ATLAS Ouvrir dans l'atlas →

Proposition 1Équations de Laplace et de la chaleur

Le noyau de la chaleur

L'équation de la chaleur sur ℝ{\mathbb R} est ∂tu=∂x2u\partial_t u = \partial_x^2 u, pour une fonction u(t,x)u(t, x) définie pour t>0t > 0 et xx réel. Le noyau de la chaleur est la fonction G(t,x)=14πtexp⁡(−x24t),G(t, x) = \frac1{\sqrt{4\pi t}}\exp\Big(-\frac{x^2}{4t}\Big), pour t>0t > 0 et xx réel, avec l'exponentielle réelle du nœud Exponentielle.

Énoncé

(i) GG est de classe C∞C^\infty et strictement positive sur ]0,+∞[×ℝ]0, +\infty[ \times \mathbb R, et ∂tG=∂x2G\partial_t G = \partial_x^2 G.

(ii) ∫ℝG(t,x) dx=1\int_{\mathbb R} G(t, x)\,dx = 1 pour tout t>0t > 0.

(iii) Pour tout δ>0\delta > 0, ∫∣x∣⩾δG(t,x) dx→0\int_{|x| \geqslant \delta} G(t, x)\,dx \to 0 quand t→0+t \to 0^+.

(iv) Pour 0<t0<t10 < t_0 < t_1 et j,k∈ℕj, k \in {\mathbb N}, il existe CC tel que ∣∂tj∂xkG(t,x)∣⩽C (1+∣x∣)2j+kexp⁡(−x24t1)\begin{aligned} &|\partial_t^j\partial_x^k G(t, x)| \\ &\quad \leqslant C\,(1 + |x|)^{2j + k}\exp\Big(-\frac{x^2}{4t_1}\Big) \end{aligned} pour t0⩽t⩽t1t_0 \leqslant t \leqslant t_1 et xx réel.

Démonstration

(i) L'exponentielle réelle est de classe C∞C^\infty, strictement positive, et exp⁡′=exp⁡\exp' = \exp (Exponentielle, (i) et (ii)) ; t↦t−1/2t \mapsto t^{-1/2} est de classe C∞C^\infty sur ]0,+∞[]0, +\infty[ (Puissances réelles, (iii)). Par composition et produit, GG est de classe C∞C^\infty (Classe Cᵏ). On calcule ∂xG=−x2t G,∂x2G=(x24t2−12t)G=∂tG.\begin{aligned} \partial_x G &= -\frac x{2t}\,G, \\ \partial_x^2 G &= \Big(\frac{x^2}{4t^2} - \frac1{2t}\Big)G = \partial_t G. \end{aligned}

(ii) Le changement de variable x=2t vx = 2\sqrt t\,v (Changement de variables, (ii)), de jacobien 2t2\sqrt t, et l'intégrale de Gauss (Intégrale de Gauss, (ii)) donnent ∫exp⁡(−x24t) dx=2t∫e−v2 dv=2πt\int \exp(-\frac{x^2}{4t})\,dx = 2\sqrt t\int e^{-v^2}\,dv = 2\sqrt{\pi t}.

(iii) Le même changement de variable donne ∫∣x∣⩾δG(t,x) dx=1π J(δ2t),J(R)=∫∣v∣⩾Re−v2 dv.\begin{aligned} \int_{|x| \geqslant \delta} G(t, x)\,dx &= \frac1{\sqrt\pi}\,J\Big(\frac\delta{2\sqrt t}\Big), \\ J(R) &= \int_{|v| \geqslant R} e^{-v^2}\,dv. \end{aligned} La fonction JJ est décroissante, et J(k)→0J(k) \to 0 quand l'entier kk tend vers l'infini, par convergence dominée, avec le majorant intégrable e−v2e^{-v^2}. Donc J(R) \to 0 quand R→+∞R \to +\infty.

(iv) Par récurrence sur j+kj + k, ∂tj∂xkG=Qj,k(1t,x) G\partial_t^j\partial_x^k G = Q_{j,k}(\frac1t, x)\,G, où Qj,kQ_{j,k} est un polynôme en deux variables de degré au plus 2j+k2j + k en xx : c'est vrai pour Q0,0=1Q_{0,0} = 1, et ∂x(QG)=(∂xQ−x2tQ)G,∂t(QG)=(∂tQ−12tQ+x24t2Q)G,\begin{aligned} \partial_x(QG) &= \Big(\partial_x Q - \frac x{2t}Q\Big)G, \\ \partial_t(QG) &= \Big(\partial_t Q - \frac1{2t}Q + \frac{x^2}{4t^2}Q\Big)G, \end{aligned} où ∂tQ\partial_t Q, calculé par la règle de la chaîne, ne change pas le degré en xx. Pour t0⩽t⩽t1t_0 \leqslant t \leqslant t_1, 1t⩽1t0\frac1t \leqslant \frac1{t_0}, d'où ∣Qj,k(1t,x)∣⩽C′(1+∣x∣)2j+k|Q_{j,k}(\frac1t, x)| \leqslant C'(1 + |x|)^{2j+k} ; 14πt⩽14πt0\frac1{\sqrt{4\pi t}} \leqslant \frac1{\sqrt{4\pi t_0}}, et exp⁡(−x24t)⩽exp⁡(−x24t1)\exp(-\frac{x^2}{4t}) \leqslant \exp(-\frac{x^2}{4t_1}), l'exponentielle étant croissante.

□
Remarque

G(t,⋅)G(t, \cdot) est la température à l'instant tt d'une barre infinie qui a reçu à l'instant 00 une quantité de chaleur 11 concentrée en 00. Le point (ii) dit que la chaleur totale se conserve ; un calcul formel donne ∫x2G(t,x) dx=2t\int x^2 G(t, x)\,dx = 2t et G(t,0)=12πtG(t, 0) = \frac1{2\sqrt{\pi t}} : la chaleur s'étale sur une largeur de l'ordre de t\sqrt t, et le pic s'abaisse d'autant.

x t = 0,05 t = 0,25 t = 1 −2 2
Le noyau G(t,⋅)G(t, \cdot) aux instants t=0,05t = 0{,}05 (rouge), t=0,25t = 0{,}25 (bleu) et t=1t = 1 (vert). Les trois courbes enferment la même aire 11 ; leurs maxima en 00, 12πt≈1,26\frac1{2\sqrt{\pi t}} \approx 1{,}26, 0,560{,}56 et 0,280{,}28, décroissent comme 1t\frac1{\sqrt t}, et les points d'inflexion, en x=±2tx = \pm\sqrt{2t}, marqués sur chaque courbe, s'écartent comme t\sqrt t.
Remarque

Le noyau est un semi-groupe : ∫G(t,x−y) G(s,y) dy=G(t+s,x)\int G(t, x - y)\,G(s, y)\,dy = G(t + s, x), égalité vérifiée par un calcul formel, qui se démontre en complétant le carré dans l'exposant. Laisser diffuser la chaleur pendant ss, puis pendant tt, revient à la laisser diffuser pendant t+st + s. Par la transformée de Fourier (Transformée de Fourier), l'égalité devient e−tξ2e−sξ2=e−(t+s)ξ2e^{-t\xi^2}e^{-s\xi^2} = e^{-(t+s)\xi^2}.

Remarque

Fourier a fondé l'étude mathématique de la propagation de la chaleur dans sa Théorie analytique de la chaleur, publiée en 1822 (W. W. Rouse Ball, A Short Account of the History of Mathematics, 1908) ; c'est pour la résoudre sur un intervalle qu'il a développé les fonctions en séries trigonométriques.

Repose sur

C'est un point de départ.

Sert à

Ouvrir dans l'atlas interactif →