Conduction thermique - Définition

Source: Wikipédia sous licence CC-BY-SA 3.0.
La liste des auteurs de cet article est disponible ici.

Conduction en régime dynamique

La résolution de l'équation de la chaleur en régime dynamique est beaucoup plus délicate. Elle fait appel aux notions de transformées de Fourier, de produit de convolution et de distributions. Nous donnons quelques exemples de résolution.

Cas d'un domaine illimité

Principe général

Ecrivons l'équation de la chaleur sous la forme :

\frac{\partial T}{\partial t} - D \Delta T = P

où D est le coefficient de diffusivité thermique et P représente des sources de chaleur. P peut être une fonction du temps et de la position de la source de chaleur, mais aussi une distribution. Par exemple, l'injection instantanée et ponctuelle d'une quantité de chaleur peut se représenter par le produit δ(t)δ(x) d'une distribution de Dirac à l'instant t = 0 par une distribution de Dirac en x = 0, x étant l'abscisse dans le cas d'un problème unidimensionnel ou le vecteur position dans le cas général.

On se donne également l'état initial du domaine T0 = T(0,x), qui peut être également une fonction de x ou une distribution. On considère que T est nulle pour t < 0.

La méthode de résolution consiste à :

  • Appliquer une transformée de Fourier relative à la variable x, à tous les termes de l'équation différentielle. Cela transforme la dérivation par rapport à x par un produit. Si on prend F(T)(p, t) = \int T(x, t) \exp(-2i\pi px) dx, alors l'équation devient :
\frac{\partial F(T)}{\partial t} + D 4\pi^2p^2F(T) = F(P)

où plutôt, au sens des distributions :

\frac{\partial F(T)}{\partial t} + D 4\pi^2p^2F(T) = F(P) + F(T_0)\delta(t)

pour tenir compte de la condition initiale.

  • Reconnaître dans cette équation un produit de convolution :
(\frac{\partial \delta(t)}{\partial t} + 4\pi^2Dp^2 \delta(t)) * F(T) = F(P) + F(T_0)\delta(t)

L'opérateur qu'on applique à F est un produit de convolution relatif à la variable t.

F(T) = F(P) * H(t)exp( − 4π2Dp2t) + F(T0)H(t)exp( − 4π2Dp2t)

Si F(P) est une fonction et non une distribution, cette relation devient, pour t > 0 :

F(T) = \int_0^t F(P)(\tau)\exp(- 4\pi^2Dp^2(t-\tau)) d\tau + F(T_0)\exp(- 4\pi^2Dp^2t)
  • Prendre la transformée de Fourier inverse pour en déduire T.

Cas particulier

Propagation par conduction dans le plan à partir d'un point chaud. La hauteur en un point donné indique la valeur de la température en ce point

Si on prend T0 = 0 et P = δ(t)δ(x) (injection instantanée de chaleur en un point donné), la méthode décrite ci-dessus conduit à :

F(P) = δ(t)

donc, pour t > 0 :

F(T) = exp( − 4π2Dp2t)

dont la transformée de Fourier inverse est, pour t > 0 :

 T = \frac{\exp(- \frac{x^2}{4Dt})}{2\sqrt{\pi tD}} dans le cas unidimensionnel.
 T = \frac{\exp(- \frac{r^2}{4Dt})}{8\sqrt{\pi tD}^3} dans le cas tridimensionnel.

Domaine illimité sans source de chaleur

Si on se donne seulement la température initiale T0 du milieu sans source de chaleur (P = 0), alors on trouve que :

T = \frac{1}{2\sqrt{\pi tD}} \int_{-\infty}^{\infty} \exp(- \frac{(x - u)^2}{4tD}) T_0(u)\, du dans le cas unidimensionnel.
T = \frac{1}{8\sqrt{\pi tD}^3} \int_{\mathbb R^3} \exp(- \frac{(r - s)^2}{4tD}) T_0(r)\, dx_s dy_s dz_s dans le cas tridimensionnel.

Cas de domaines limités, sans source de chaleur

Cas d'un domaine limité par un plan. Le problème de Kelvin

Problème de Kelvin. L'axe des x est orienté vers la droite. Le demi-espace x>0, dont la température initiale est T0, possède une frontière x=0 dont la température est constamment nulle. La hauteur en un point donné indique la valeur de la température en ce point

Supposons le domaine limité par le plan x=0. Si on se donne pour condition aux limites supplémentaire T(0,t) = 0 pour tout t, alors, il suffit de prolonger la répartition initiale de température T0 par une fonction impaire en x et d'appliquer le résultat précédent.

Le cas le plus célèbre est celui du problème de Kelvin. Ce dernier a considéré dans les années 1860 que la Terre était initialement à une température constante T0 de l'ordre de 3000° et qu'elle s'est refroidie par simple conduction. Utilisant la valeur actuelle du gradient de température en fonction de la profondeur, il en a déduit une estimation de l'âge de la Terre. On peut appliquer la méthode de résolution précédente en considérant la Terre comme plate et infiniment profonde, limitée par le plan de sa surface. Le calcul conduit à :

T= \frac{T_0}{\sqrt{\pi t D}} \int_0^x \exp(- \frac{u^2}{4tD})\, du = T_0 \,{\rm erf}(\frac{x}{2\sqrt{Dt}})

où erf est dite fonction d'erreur de Gauss.

Le gradient de température à la surface est :

\frac{\partial T}{\partial x}= \frac{T_0}{\sqrt{\pi t D}}

Connaissant \frac{\partial T}{\partial x} de l'ordre de 3 °C pour 100 mètres de profondeur et D estimé à 10^{-6} \,m^2s^{-1}, on trouve que t vaut 100 millions d'années. Ce résultat est largement sous-estimé car Kelvin ignorait les phénomènes de convection au sein du manteau terrestre.

Cas d'un domaine limité par deux plans parallèles

Propagation de la chaleur par conduction dans un domaine limité 0 < x < L. Les deux frontières du domaine sont maintenues à température constante. La hauteur en un point donné indique la valeur de la température en ce point

Considérons un domaine limité par les deux plans x=0 et x=L. Supposons qu'on se donne comme conditions aux limites T(0,t) = T(L,t) = 0. On utilise une méthode de résolution basée sur les séries de Fourier, en cherchant T sous la forme :

T = \sum_{n=1}^\infty b_n \sin(n\pi \frac{x}{L}) \exp(- \frac{Dn^2\pi^2t}{L^2})

Cette expression vérifie à la fois l'équation de la chaleur et les conditions aux limites. Si on se donne la répartition de température initiale T0, il suffit de développer celle-ci en série de Fourier pour déterminer les bn.

Par exemple, si on prend T0 constant, on obtient :

T = \frac{4T_0}{\pi} \sum_{n=0}^\infty \frac{1}{2n+1} \sin(\frac{(2n+1)\pi x}{L}) \exp(- \frac{D(2n+1)^2\pi^2t}{L^2})

En faisant tendre L vers l'infini, on retrouve la solution de Kelvin du , la somme précédente étant considérée comme une somme de Riemann convergeant vers l'intégrale.

Cas d'un domaine à géométrie sphérique

Propagation de la chaleur par conduction dans un domaine circulaire. Problème de Kelvin : la température initiale est uniforme, la température sur le cercle frontière est maintenue nulle. La hauteur en un point donné indique la valeur de la température en ce point

Dans le cas où la propagation se fait dans un domaine sphérique, et où la température ne dépend que de la distance r au centre, l'équation de la chaleur devient, compte tenu de l'expression du laplacien en sphérique :

\frac{\partial T}{\partial t} = D(\frac{2}{r} \frac{\partial T}{\partial r} + \frac{\partial^2 T}{\partial r^2})

Si on pose F = rT, l'équation devient :

\frac{\partial F}{\partial t} = D \frac{\partial^2 F}{\partial r^2}

On peut alors appliquer les méthodes précédentes pour déterminer F, puis en déduire T en divisant par r.

Ainsi, la résolution du problème de Kelvin dans le cas d'une boule de rayon R (température initiale uniformément égale à T0, la surface étant maintenue à une température nulle) conduit à l'expression suivante de T :

T(r,t) = 2T_0 \sum_{n=1}^\infty (-1)^{n+1} \, {\rm sinc}(n\pi \frac{r}{R}) \exp(- \frac{Dn^2\pi^2t}{R^2})

où sinc est la fonction sinus cardinal.

Cas de domaines limités, avec source de chaleur

On considère l'équation :

\frac{\partial T}{\partial t} - D \Delta T = P

avec P non nul. On cherche en général une solution particulière à cette équation, de façon à ce que, une fois retranchée à T, on puisse se ramener à une équation sans second membre. Voici quelques exemples, dans le cas où P représente une densité de source de chaleur constante, indépendante de la position et du temps.

Domaine limité par deux plans parallèles

Propagation de la chaleur par conduction dans un domaine limité 0 < x < L. Le domaine est chauffé uniformément, alors que la température initiale est nulle, et que les deux bords restent à cette température nulle. La hauteur en un point donné indique la valeur de la température en ce point.

Considérons un domaine limité par les deux plans x=0 et x=L. On suppose qu'à l'instant initial, la température du domaine est égale à une température de référence nulle, et que les bords du domaine resteront en permanence à cette température nulle. T vérifie donc :

\frac{\partial T}{\partial t} - D \frac{\partial^2 T}{\partial x^2} = P
T(0,t) = T(L,t) = 0 pour tout t positif.
T(x,0) = 0 pour tout x entre 0 et L.

La fonction \frac{Px(L-x)}{2D} indépendante de t vérifie les deux premières relations, de sorte que, si on pose G = T - \frac{Px(L-x)}{2D}, alors G vérifie :

\frac{\partial G}{\partial t} - D \frac{\partial^2 G}{\partial x^2} = 0
G(0,t) = G(L,t) = 0
G(x, 0) = - \frac{Px(L-x)}{2D}

On peut appliquer la méthode en cherchant G sous la forme d'une série :

G = \sum_{n=1}^\infty b_n \sin(\frac{n\pi x}{L}) \exp(- \frac{n^2\pi^2Dt}{L^2})

qui vérifie les deux premières relations. Comme, pour des raisons de symétrie, on s'attend à ce que G(x) = G(Lx), on peut supposer que les coefficients bn sont nuls lorsque n est pair, de sorte que :

G = \sum_{n=0}^\infty b_{2n+1} \sin(\frac{(2n+1)\pi x}{L}) \exp(- \frac{(2n+1)^2\pi^2Dt}{L^2})

Pour t = 0, on a :

- \frac{Px(L-x)}{2D} = \sum_{n=0}^\infty b_{2n+1} \sin(\frac{(2n+1)\pi x}{L})

On trouve les b2n + 1 en développant - \frac{Px(L-x)}{2D} en série de Fourier. On trouve :

b_{2n+1} = - \frac{4PL^2}{(2n+1)^3\pi^3D}

D'où G, puis finalement :

T = \frac{Px(L-x)}{2D} - \frac{4PL^2}{D\pi^3} \sum_{n=0}^\infty  \frac{1}{(2n+1)^3} \sin(\frac{(2n+1)\pi x}{L}) \exp(- \frac{(2n+1)^2\pi^2Dt}{L^2})

Lorsque t tend vers l'infini, la température du domaine tend vers \frac{Px(L-x)}{2D}, l'échauffement thermique dans le milieu étant alors en équilibre avec l'évacuation de la chaleur par les deux bords.

Domaine limité par un plan

Propagation de la chaleur par conduction dans un domaine x > 0 limité par un bord. Le domaine est chauffé uniformément, alors que la température initiale est nulle, et que le bord reste à cette température nulle. La hauteur en un point donné indique la valeur de la température en ce point.

La résolution du même problème dans le cas où x>0 consiste à déterminer T tel que  :

\frac{\partial T}{\partial t} - D \frac{\partial^2 T}{\partial x^2} = P
T(0,t) = 0 pour tout t positif.
T(x,0) = 0 pour tout x > 0.

On peut obtenir la solution en faisant tendre L vers l'infini dans l'expression donnée dans le , en assimilant la série à une somme de Riemann. On obtient alors l'expression suivante :

T = - \frac{Px^2}{2D} + \frac{Px^2}{2D} \,{\rm erf}(\frac{x}{2\sqrt{Dt}}) + \frac{Px\sqrt{t}}{\sqrt{D\pi}} \exp(- \frac{x^2}{4Dt}) + Pt \, {\rm erf}(\frac{x}{2\sqrt{Dt}})

erf est la fonction dite fonction d'erreur de Gauss. On peut également trouver cette expression en appliquant la méthode découlant du relatif à un domaine illimité, après avoir étendu à l'espace entier les fonctions T et P en des fonctions impaires en x, de façon à ce que T s'annule en x = 0.

Quand t tend vers l'infini, T vaut environ Pt, analogue à celle d'un domaine infini. Le bord unique n'est pas suffisant pour évacuer la chaleur.

Domaine à géométrie sphérique

Propagation de la chaleur par conduction dans un domaine à bord circulaire. Le domaine est chauffé uniformément, alors que la température initiale est nulle, et que le bord reste à cette température nulle. La hauteur en un point donné indique la valeur de la température en ce point.

Dans le cas d'un domaine dont le bord est une sphère de rayon R, on utilise l'expression du laplacien en sphérique et on est amené à résoudre :

\frac{\partial T}{\partial t} = D\left(\frac{2}{r} \frac{\partial T}{\partial r} + \frac{\partial^2 T}{\partial r^2}\right) + P
Pour tout t, T(R,t) = 0
Pour tout r, T(r,0) = 0

En posant G = rT + \frac{r^3P - rR^2P}{6D}, G vérifie le système :

\frac{\partial G}{\partial t} = D\frac{\partial^2 G}{\partial r^2}
Pour tout t, G(R,t) = 0
Pour tout r, G(r,0) = \frac{r^3P - rR^2P}{6D}

La méthode des séries de Fourier suggère de chercher G sous la forme d'une série \sum_{n=1}^\infty b_n \sin(\frac{n\pi r}{R}) \exp(- \frac{n^2\pi^2 Dt}{R^2}), où les bn sont trouvés en développant \frac{r^3P - rR^2P}{6D} en série de Fourier. On obtient :

G = \frac{2P R^3}{D\pi^3} \sum_{n=1}^\infty \frac{(-1)^n}{n^3} \sin(\frac{n\pi r}{R}) \exp(- \frac{n^2\pi^2 Dt}{R^2})

et donc :

T = \frac{R^2P - r^2P}{6D} + \frac{2P R^2}{D\pi^2} \sum_{n=1}^\infty \frac{(-1)^n}{n^2} {\rm sinc}(\frac{n\pi r}{R}) \exp(- \frac{n^2\pi^2 Dt}{R^2})

où sinc est la fonction sinus cardinal.

Quand t tend vers l'infini, la température T tend vers la répartition limite \frac{R^2P - r^2P}{6D}.

Page générée en 0.123 seconde(s) - site hébergé chez Contabo
Ce site fait l'objet d'une déclaration à la CNIL sous le numéro de dossier 1037632
A propos - Informations légales | Partenaire: HD-Numérique
Version anglaise | Version allemande | Version espagnole | Version portugaise