Introduction aux transferts thermiques/Équation de la chaleur

On a vu au [[../Concepts généraux#Équation de la chaleur, cas simple|chapitre n°1]] une mise en équation locale du phénomène de transfert thermique dans un corps. Cette approche ne traitait qu'une partie des questions liées à cette mise en équation. On traite ici un cas plus général.

Le système considéré, de volume $V$ et de surface externe $\mathrm{\Sigma }$, est indéformable. Il en résulte que la capacité thermique massique à volume constant est égale à la capacité thermique massique à pression constante (J kg^-1 K^-1) : $c_v=c_p=c$.

Le taux de variation d'énergie interne $U$ du système s'écrit : $\underset{V}{\iiint }\rho c\frac{\partial T}{\partial t}\mathrm{d}V$, où $\rho$ est la masse volumique (kg m^-3)
L'énergie reçue par le système à ses frontières est : $\int \subset \supset \underset{\mathrm{\Sigma }}{\int }-\overrightarrow{\phi }\cdot \overrightarrow{n}\mathrm{d}\mathrm{\Sigma }=\underset{V}{\iiint }-\mathrm{d}\mathrm{i}\mathrm{v}\overrightarrow{\phi }\mathrm{d}V$. On convient du signe positif quand le système reçoit effectivement cette énergie.

Le système peut éventuellement avoir des sources internes d'énergie (effet Joule, réaction chimique, etc.) qui s'écrivent : $\underset{V}{\iiint }q\mathrm{d}V$.Selon le premier principe de la thermodynamique :

$\underset{V}{\iiint }\rho c\frac{\partial T}{\partial t}\mathrm{d}V=\int \subset \supset \underset{\mathrm{\Sigma }}{\int }-\overrightarrow{\phi }\cdot \overrightarrow{n}\mathrm{d}\mathrm{\Sigma }+\underset{V}{\iiint }q\mathrm{d}V$.

Le théorème de la divergence permet de ramener cette équation de bilan à une équation locale en chaque point, ci-dessous.

Hypothèses supplémentaires :

• il n'y a aucun transfert d'énergie ne se produisant par déplacement de matière, il n'y a pas de transfert thermique par convection ;
• le milieu est opaque, il n'y a pas de transfert thermique par rayonnement ;
• le milieu est isotrope et homogène, les propriétés thermophysiques sont les même dans toutes les directions et dans tout le système.

Il en résulte que le flux thermique n'est dû qu'à la conduction thermique :

$\overrightarrow{\phi }=-\lambda \overrightarrow{\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}}T$ , où $\lambda$ est la conductivité thermique (W K^-1 m^-1).

Dans le cas général, les propriétés thermophysiques $\lambda$, $\rho$, $c$ sont dépendantes de la température. Elle ne sont cependant pas dérivées spatialement. On peut alors écrire :

$-\mathrm{d}\mathrm{i}\mathrm{v}\overrightarrow{\phi }=\mathrm{d}\mathrm{i}\mathrm{v}(\lambda (T)\overrightarrow{\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}}T)$;

$-\mathrm{d}\mathrm{i}\mathrm{v}\overrightarrow{\phi }=\lambda (T)\mathrm{d}\mathrm{i}\mathrm{v}\left(\overrightarrow{\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}}T\right)+\overrightarrow{\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}}T\overrightarrow{\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}}\lambda$ ;

$-\mathrm{d}\mathrm{i}\mathrm{v}\overrightarrow{\phi }=\lambda (T)\mathrm{\Delta }T+\frac{\partial \lambda (T)}{\partial T}{\left(\overrightarrow{\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}}T\right)}^2$.

L'équation de la chaleur devient :

$\rho C\frac{\partial T}{\partial t}=\lambda (T)\mathrm{\Delta }T+\frac{\partial \lambda (T)}{\partial T}{\left(\overrightarrow{\mathrm{g}\mathrm{r}\mathrm{a}\mathrm{d}}T\right)}^2+q$

Si, pour simplifier les calculs, on peut considérer la conduction thermique constante, on a Sans la thermodépendance on a :

$\mathrm{\Delta }T-\frac{\rho c}{\lambda }\frac{\partial T}{\partial t}=-\frac{q}{\lambda }$.

On peut introduire la diffusivité thermique $a=\frac{\lambda }{\rho C}$ en m² s^-1. Modèle:Encart