Séries Infinies : Trouver Les Constantes Pour L'équation De Diffusion

Salut les amis ! Aujourd'hui, on plonge dans un sujet super intéressant qui peut faire transpirer certains : comment résoudre pour les constantes dans une équation impliquant une série de sommes infinies. C'est un défi qu'on rencontre souvent quand on utilise la méthode de séparation des variables pour résoudre des équations aux dérivées partielles comme l'équation de diffusion. Si vous avez déjà jeté un œil à l'équation de diffusion et que vous vous êtes retrouvés face à une expression du type $c(x,t) = som_{n=0}^{\infty}A_n\exp

left(-\frac{Dn^2\pi2t}{(N+1)^2\delta2}\right)\sin

left(\frac{n\pi x}{N+1}\right)$, alors vous savez de quoi je parle. Ces fameuses constantes $A_n$ peuvent sembler un peu mystérieuses au début, mais ne vous inquiétez pas, on va démystifier tout ça ensemble. L'objectif est de trouver ces valeurs de $A_n$ pour que notre solution corresponde aux conditions initiales de notre problème. Préparez votre café, car ça va être une petite exploration dans le monde fascinant des séries infinies et de la physique mathématique !

Démystifier les constantes $A_n$ : Le cœur du problème

Alors les gars, le vrai casse-tête ici, c'est de déterminer ces fichues constantes $A_n$ dans notre série infinie. Souvenez-vous, on est partis de l'équation de diffusion et on a utilisé la séparation des variables. Ça nous a menés à une solution qui est une somme de termes, chacun étant le produit d'une constante $A_n$, d'une fonction exponentielle qui dépend du temps $t$, et d'une fonction sinus qui dépend de la position $x$. Notre équation ressemble à ça : $c(x,t) = som_{n=0}^{\infty}A_n\exp

left(-\frac{Dn^2\pi2t}{(N+1)^2\delta2}\right)\sin

left(\fracn\pi x}{N+1}\right)$. Le truc, c'est que cette formule générale est super utile, mais elle ne vaut vraiment que lorsqu'on connaît les valeurs précises des $A_n$. Ces constantes sont directement liées à la condition initiale de notre problème de diffusion. Autrement dit, c'est ce que notre système (par exemple, la température ou la concentration) ressemble au tout début, quand $t=0$. Si on pose $t=0$ dans notre équation, le terme exponentiel devient $e^0 = 1$. Ça simplifie drastiquement notre formule pour devenir $c(x,0) = som_{n=0^{\infty}A_n\sin

left(\frac{n\pi x}{N+1}\right)$. Et là, mes amis, on a affaire à une série de Fourier ! La fonction $c(x,0)$ représente notre distribution initiale, et on cherche à l'exprimer comme une somme de fonctions sinus. Les constantes $A_n$ sont en fait les coefficients de cette série de Fourier. Pour trouver ces coefficients, on utilise une propriété super cool des séries de Fourier et des fonctions orthogonales. En gros, on multiplie notre équation par $\sin

left(\fracm\pi x}{N+1}\right)$ (où $m$ est un autre indice entier) et on intègre sur l'intervalle de $x$ qui nous intéresse, qui est généralement de $0$ à $N+1$ (ou $L$ dans une notation plus générale). Grâce à l'orthogonalité des fonctions sinus, la plupart des termes de la somme s'annulent, sauf celui où $n=m$. Ça nous donne une formule explicite pour calculer chaque $A_m$ $A_m = \frac{2{N+1}\int_{0}^{N+1}c(x,0)\sin

left(\frac{m\pi x}{N+1}\right)dx$. Voilà, le mystère est levé ! Les $A_n$ ne sont plus des inconnus, mais des valeurs que l'on peut calculer directement à partir de notre état initial. C'est comme si on décomposait notre état initial complexe en une somme de composantes sinus plus simples, et les $A_n$ nous disent quelle est l'amplitude de chaque composante. C'est une technique super puissante en physique mathématique.

La magie de l'orthogonalité des fonctions sinus

Parlons un peu plus en détail de ce truc génial qu'est l'orthogonalité des fonctions sinus. Sans elle, trouver nos constantes $A_n$ serait un vrai cauchemar, genre devoir résoudre un système infini d'équations ! Heureusement, les fonctions sinus, dans un intervalle donné, possèdent cette propriété d'orthogonalité. Qu'est-ce que ça veut dire concrètement ? Eh bien, si vous prenez deux fonctions sinus différentes, disons $\sin

left(\frac{n\pi x}{L}\right)$ et $\sin

left(\frac{m\pi x}{L}\right)$ avec $n \neq m$, et que vous calculez leur produit scalaire sur l'intervalle $[0, L]$ (ici, $L$ est remplacé par $N+1$ dans notre cas), le résultat est zéro. Mathématiquement, ça s'écrit : $

som_{n \neq m} \int_{0}^{L}\sin

left(\frac{n\pi x}{L}\right)\sin

left(\fracm\pi x}{L}\right)dx = 0$ . Ce résultat est fondamental. Reprenons notre équation pour la condition initiale à $t=0$ $c(x,0) = som_{n=1^{\infty}A_n\sin

left(\frac{n\pi x}{N+1}\right)$ (Notez que j'ai commencé la somme à $n=1$ car pour $n=0$, le sinus est nul, donc ce terme ne contribue pas). Maintenant, imaginez qu'on veuille trouver $A_m$. On va multiplier toute l'équation par $\sin

left(\frac{m\pi x}{N+1}\right)$ : $ c(x,0)\sin

left(\frac{m\pi x}{N+1}\right) = som_{n=1}^{\infty}A_n\sin

left(\frac{n\pi x}{N+1}\right)\sin

left(\fracm\pi x}{N+1}\right)$. Et ensuite, on intègre des deux côtés, de $x=0$ à $x=N+1$ $\int_{0^{N+1}c(x,0)\sin

left(\frac{m\pi x}{N+1}\right)dx = \int_{0}^{N+1}\left( som_{n=1}^{\infty}A_n\sin

left(\frac{n\pi x}{N+1}\right)\sin

left(\fracm\pi x}{N+1}\right)\right)dx$ . On peut intervertir l'intégrale et la somme (sous certaines conditions de convergence, mais dans nos problèmes physiques, c'est généralement le cas) $\int_{0^{N+1}c(x,0)\sin

left(\frac{m\pi x}{N+1}\right)dx = som_{n=1}^{{\infty}A_n\int_{0}}{N+1}\sin

left(\frac{n\pi x}{N+1}\right)\sin

left(\frac{m\pi x}{N+1}\right)dx$ . Et c'est là que l'orthogonalité fait des miracles ! Pour tous les termes où $n \neq m$, l'intégrale $\int_{0}^{N+1}\sin

left(\frac{n\pi x}{N+1}\right)\sin

left(\fracm\pi x}{N+1}\right)dx$ vaut zéro. Il ne reste donc que le terme où $n=m$ $\int_{0^{N+1}c(x,0)\sin

left(\frac{m\pi x}{N+1}\right)dx = A_m\int_{0}^{N+1}\sin2

left(\frac{m\pi x}{N+1}\right)dx$ . L'intégrale restante $\int_{0}^{N+1}\sin2

left(\fracm\pi x}{N+1}\right)dx$ a une valeur connue, qui est $\frac{N+1}{2}$ (pour $m \neq 0$). Donc, on obtient $\int_{0^{N+1}c(x,0)\sin

left(\fracm\pi x}{N+1}\right)dx = A_m \frac{N+1}{2}$ . En réarrangeant pour isoler $A_m$, on tombe sur la formule qu'on a vue précédemment $A_m = \frac{2{N+1}\int_{0}^{N+1}c(x,0)\sin

left(\frac{m\pi x}{N+1}\right)dx$ . C'est cette propriété d'orthogonalité qui nous permet de