Transfert de chaleur par conduction

Physique de la conduction thermique

Au niveau atomique, la conduction thermique résulte de deux mécanismes. Dans les solides non métalliques, la chaleur est transportée principalement par les vibrations du réseau cristallin — des paquets quantifiés d'énergie vibratoire appelés phonons. Lorsqu'une région du réseau cristallin est plus chaude, ses atomes vibrent plus vigoureusement, et ces vibrations se propagent aux atomes voisins, transportant l'énergie. Dans les métaux, un second mécanisme — généralement dominant — entre en jeu : les électrons libres. La mer d'électrons de conduction délocalisés qui confère aux métaux leur conductivité électrique sert également de vecteur extrêmement efficace d'énergie thermique. C'est pourquoi les métaux sont à la fois de bons conducteurs électriques et thermiques, une connexion formalisée par la loi de Wiedemann-Franz.

La conductivité thermique $k$ d'un matériau quantifie la facilité avec laquelle il conduit la chaleur. Le cuivre, avec $k \approx 400 \; \mathrm{W/(m \cdot K)}$, figure parmi les meilleurs conducteurs. Le diamant, grâce à son réseau cristallin extrêmement rigide, atteint $k \approx 2000 \; \mathrm{W/(m \cdot K)}$. À l'autre extrême, l'aérogel — un matériau à base de silice constitué principalement de poches d'air piégées — atteint $k \approx 0{,}015 \; \mathrm{W/(m \cdot K)}$, ce qui en fait l'un des meilleurs isolants thermiques connus.

Loi de Fourier de la conduction thermique

La pierre angulaire du transfert de chaleur conductif est la loi de Fourier, proposée par Jean-Baptiste Joseph Fourier en 1822. En une dimension, elle stipule que le flux de chaleur — le taux de transfert de chaleur par unité de surface — est proportionnel au gradient de température négatif :

$${\displaystyle q = -k \, \frac{dT}{dx}}$$

Où :

• $q$ est le flux de chaleur $\left[\mathrm{W/m^2}\right]$.
• $k$ est la conductivité thermique du matériau $\left[\mathrm{W/(m \cdot K)}\right]$.
• $dT/dx$ est le gradient de température dans la direction du flux de chaleur $\left[\mathrm{K/m}\right]$.

Le signe négatif reflète le second principe de la thermodynamique : la chaleur s'écoule spontanément du chaud vers le froid. Le taux total de transfert de chaleur à travers une surface d'aire $A$ est simplement :

$${\displaystyle \dot{Q} = q \, A = -k \, A \, \frac{dT}{dx}}$$

En trois dimensions, la loi de Fourier se généralise en une équation vectorielle impliquant le gradient de température $\nabla T$ :

$${\displaystyle \vec{q} = -k \, \nabla T}$$

Cette expression compacte est le point de départ de pratiquement tous les traitements analytiques et numériques de la conduction thermique.

L'équation de la chaleur

En combinant la loi de Fourier avec la conservation de l'énergie, on obtient l'équation de la chaleur (également appelée équation de diffusion), l'équation aux dérivées partielles fondamentale régissant la conduction thermique transitoire dans un solide :

$${\displaystyle \rho \, c_p \, \frac{\partial T}{\partial t} = \nabla \cdot \left( k \, \nabla T \right) + \dot{q}_{\text{gen}}}$$

Où :

• $\rho$ est la masse volumique du matériau $\left[\mathrm{kg/m^3}\right]$.
• $c_p$ est la capacité thermique massique $\left[\mathrm{J/(kg \cdot K)}\right]$.
• $\dot{q}_{\text{gen}}$ est le taux de génération interne de chaleur volumique $\left[\mathrm{W/m^3}\right]$, par exemple dû à un chauffage par résistance électrique ou à des réactions chimiques.

Pour un matériau homogène avec une conductivité thermique constante, cela se simplifie en :

$${\displaystyle \frac{\partial T}{\partial t} = \alpha \, \nabla^2 T + \frac{\dot{q}_{\text{gen}}}{\rho \, c_p}}$$

où $\alpha = k / (\rho \, c_p)$ est la diffusivité thermique $\left[\mathrm{m^2/s}\right]$ — une propriété unique qui indique la rapidité de propagation des variations de température dans un matériau. Le cuivre, avec $\alpha \approx 1{,}17 \times 10^{-4} \; \mathrm{m^2/s}$, répond thermiquement bien plus rapidement que la brique à $\alpha \approx 5{,}2 \times 10^{-7} \; \mathrm{m^2/s}$.

Conduction en régime permanent : résistance thermique

Lorsque les conditions ne varient pas dans le temps ($\partial T / \partial t = 0$) et qu'il n'y a pas de génération interne de chaleur, l'équation de la chaleur se réduit à l'équation de Laplace : $\nabla^2 T = 0$. Pour le cas simple mais extrêmement utile de la conduction unidimensionnelle en régime permanent à travers une paroi plane d'épaisseur $L$, le profil de température est linéaire et le taux de transfert de chaleur est :

$${\displaystyle \dot{Q} = \frac{k \, A}{L} \left(T_1 - T_2\right) = \frac{T_1 - T_2}{R_{\text{cond}}}}$$

où la résistance thermique pour la conduction à travers la paroi est définie comme :

$${\displaystyle R_{\text{cond}} = \frac{L}{k \, A}}$$

Cette analogie avec la loi d'Ohm ($V = IR$) est extraordinairement puissante. De même que les résistances électriques peuvent être combinées en série et en parallèle, les résistances thermiques peuvent être enchaînées pour modéliser les parois composites, les couches d'isolation, les résistances de contact et les conditions aux limites convectives dans un cadre unifié unique.

Parois composites

Considérons une paroi composée de trois couches (par exemple plâtre, brique, isolation) avec des conditions aux limites convectives des deux côtés. La résistance thermique totale est :

$${\displaystyle R_{\text{total}} = \frac{1}{h_1 A} + \frac{L_1}{k_1 A} + \frac{L_2}{k_2 A} + \frac{L_3}{k_3 A} + \frac{1}{h_2 A}}$$

et le taux de transfert de chaleur global s'obtient directement par $\dot{Q} = \Delta T_{\text{global}} / R_{\text{total}}$. Cette approche est le fondement des calculs énergétiques des bâtiments et est utilisée quotidiennement par les architectes et les ingénieurs en mécanique pour évaluer les performances d'isolation et les valeurs U.

Conduction radiale : cylindres et sphères

De nombreux composants d'ingénierie — tuyaux, câbles, cuves de réacteurs — ont une géométrie cylindrique ou sphérique. Pour la conduction en régime permanent à travers un cylindre creux (rayon interne $r_1$, rayon externe $r_2$, longueur $L_{\text{cyl}}$), le profil de température est logarithmique et le taux de transfert de chaleur est :

$${\displaystyle \dot{Q} = \frac{2\pi \, k \, L_{\text{cyl}} \left(T_1 - T_2\right)}{\ln\left(r_2 / r_1\right)}}$$

La résistance thermique correspondante pour la conduction radiale à travers une coque cylindrique est :

$${\displaystyle R_{\text{cyl}} = \frac{\ln\left(r_2 / r_1\right)}{2\pi \, k \, L_{\text{cyl}}}}$$

Une conséquence intéressante survient lorsqu'on ajoute de l'isolation à un tuyau ou un fil mince : il existe un rayon critique d'isolation $r_{\text{cr}} = k_{\text{ins}} / h$ en dessous duquel l'ajout d'isolation augmente en fait la perte de chaleur, car la plus grande surface extérieure améliore la convection davantage que l'isolation ne réduit la conduction. Ce résultat contre-intuitif a des implications pratiques pour l'isolation des câbles électriques et des tubulures de petit diamètre.

Conduction transitoire et nombre de Biot

Lorsque les températures varient dans le temps — une pièce trempée, un four en préchauffage, un disque de frein absorbant de la chaleur par frottement — nous entrons dans le domaine de la conduction transitoire. La première question est de savoir si la température intérieure du corps peut être traitée comme spatialement uniforme ou si des gradients internes significatifs se développent. Le nombre de Biot fournit la réponse :

$${\displaystyle \mathrm{Bi} = \frac{h \, L_c}{k_s}}$$

où $h$ est le coefficient convectif de surface, $L_c$ est une longueur caractéristique (typiquement volume/surface), et $k_s$ est la conductivité thermique du solide. Lorsque $\mathrm{Bi} < 0{,}1$, les gradients de température internes sont négligeables et la méthode de la capacité globale s'applique — le corps entier est traité comme étant à une seule température dépendant du temps, décrite par la loi de refroidissement de Newton :

$${\displaystyle T(t) = T_{\infty} + \left(T_0 - T_{\infty}\right) e^{-t/\tau}}$$

avec une constante de temps $\tau = \rho \, V \, c_p / (h \, A_s)$. Lorsque $\mathrm{Bi} > 0{,}1$, les variations spatiales de température à l'intérieur du corps deviennent importantes et l'équation complète de la chaleur doit être résolue — analytiquement pour les géométries simples (par séparation de variables et diagrammes de Heisler) ou numériquement par la méthode des éléments finis pour les formes complexes.

Applications en ingénierie

1. Isolation des bâtiments et efficacité énergétique

La conduction à travers les murs, toitures et planchers est la voie principale de perte de chaleur dans les bâtiments. Les ingénieurs utilisent l'approche du réseau de résistances thermiques pour évaluer la valeur U globale des assemblages de construction et déterminer où l'isolation supplémentaire offre le meilleur retour. Des matériaux comme le polystyrène expansé ($k \approx 0{,}035 \; \mathrm{W/(m \cdot K)}$), la laine minérale et les panneaux d'isolation sous vide sont sélectionnés spécifiquement pour leur faible conductivité thermique.

2. Gestion thermique de l'électronique

En microélectronique, la chaleur générée par les transistors doit être conduite à travers la puce de silicium, à travers les matériaux d'interface thermique (TIM), puis dans le dissipateur et le radiateur avant de pouvoir être évacuée par convection. Chaque couche présente une résistance thermique, et la résistance totale jonction-ambiance détermine la température de fonctionnement de la puce. Réduire n'importe quel maillon de cette chaîne — par exemple en utilisant des dissipateurs en cuivre ou en diamant à haute conductivité — abaisse directement la température de jonction et améliore la fiabilité.

3. Traitement thermique métallurgique

Les processus tels que le recuit, la trempe et le revenu dépendent de manière critique du contrôle de la vitesse et de l'uniformité du changement de température dans une pièce métallique. La conduction interne détermine la rapidité avec laquelle le cœur d'une forge épaisse en acier atteint la température cible par rapport à sa surface. Une erreur entraîne des contraintes résiduelles, des déformations ou des microstructures indésirables. L'analyse thermique prédictive par FEA est largement utilisée pour concevoir des cycles de traitement thermique qui atteignent les propriétés souhaitées tout en minimisant les défauts.

4. Systèmes de protection thermique

Les engins spatiaux en rentrée atmosphérique subissent un chauffage extrême — les températures de surface peuvent dépasser 1500 °C. Les boucliers thermiques ablatifs et les systèmes de tuiles céramiques exploitent des matériaux à faible conductivité thermique pour maintenir la structure du véhicule au frais, gagnant du temps pour que l'impulsion thermique passe avant de pénétrer jusqu'aux composants critiques. La conception de ces systèmes implique la résolution de problèmes de conduction transitoire avec des propriétés de matériaux dépendant de la température et une ablation de surface.

Étude de cas : perte de chaleur à travers une canalisation isolée

Un tuyau en acier ($k_{\text{acier}} = 50 \; \mathrm{W/(m \cdot K)}$) transporte de la vapeur à 200 °C. Le tuyau a un rayon interne $r_1 = 25 \; \mathrm{mm}$, un rayon externe $r_2 = 30 \; \mathrm{mm}$, et est recouvert de 40 mm de laine minérale isolante ($k_{\text{iso}} = 0{,}04 \; \mathrm{W/(m \cdot K)}$, rayon externe $r_3 = 70 \; \mathrm{mm}$). La température ambiante est de 20 °C avec un coefficient convectif extérieur $h_o = 10 \; \mathrm{W/(m^2 \cdot K)}$. Par mètre de longueur de tuyau ($L = 1 \; \mathrm{m}$) :

$${\displaystyle R_{\text{acier}} = \frac{\ln(r_2/r_1)}{2\pi \, k_{\text{acier}} \, L} = \frac{\ln(30/25)}{2\pi \times 50 \times 1} = \frac{0{,}1823}{314{,}16} \approx 5{,}8 \times 10^{-4} \; \mathrm{K/W}}$$

$${\displaystyle R_{\text{iso}} = \frac{\ln(r_3/r_2)}{2\pi \, k_{\text{iso}} \, L} = \frac{\ln(70/30)}{2\pi \times 0{,}04 \times 1} = \frac{0{,}8473}{0{,}2513} \approx 3{,}371 \; \mathrm{K/W}}$$

$${\displaystyle R_{\text{conv}} = \frac{1}{h_o \, 2\pi \, r_3 \, L} = \frac{1}{10 \times 2\pi \times 0{,}07 \times 1} = \frac{1}{4{,}398} \approx 0{,}227 \; \mathrm{K/W}}$$

La résistance totale est $R_{\text{total}} = 5{,}8 \times 10^{-4} + 3{,}371 + 0{,}227 \approx 3{,}60 \; \mathrm{K/W}$, et la perte de chaleur par mètre est :

$${\displaystyle \dot{Q} = \frac{200 - 20}{3{,}60} \approx 50 \; \mathrm{W/m}}$$

Remarquez que la paroi en acier ne contribue presque rien à la résistance totale ($0{,}016\%$), tandis que la couche d'isolation fournit $93{,}6\%$ de la résistance thermique. Sans isolation, la perte de chaleur bondirait à environ 800 W/m — soit une augmentation d'un facteur seize. Ce type d'analyse de résistance thermique est courant dans la conception de tuyauteries industrielles, l'ingénierie des procédés et les audits énergétiques.

Méthodes numériques pour les problèmes de conduction complexes

Les solutions analytiques n'existent que pour des géométries et des conditions aux limites relativement simples. En pratique, les composants d'ingénierie ont des formes irrégulières, des propriétés de matériaux dépendant de la température, des sources de chaleur internes, des résistances de contact aux interfaces et des assemblages multi-matériaux complexes. Pour ces problèmes, la méthode des éléments finis (FEM) est l'outil numérique standard. La FEM discrétise la géométrie en un maillage d'éléments, approxime le champ de température par des fonctions polynomiales par morceaux et résout le système d'équations algébriques résultant pour produire une carte détaillée de la température dans tout le composant.

Les logiciels FEA modernes gèrent la conduction en régime permanent et transitoire, couplée à des conditions aux limites convectives et radiatives, des propriétés de matériaux non linéaires, des changements de phase, et même un couplage thermo-mécanique où les contraintes et déformations thermiques interagissent avec le problème thermique. Lorsque l'écoulement du fluide lui-même fait partie du problème — convection forcée ou naturelle, transfert thermique conjugué à travers un échangeur de chaleur — le Computational Fluid Dynamics (CFD) prend le relais, résolvant les équations d'écoulement et d'énergie simultanément. Ensemble, la FEA et le CFD constituent l'épine dorsale de l'analyse thermique moderne, permettant aux ingénieurs de prédire les températures, les flux de chaleur et les contraintes thermiques dans des composants de toute complexité.

Conclusion

La conduction thermique est le mode de transfert de chaleur le plus intuitif et le plus aisé à traiter mathématiquement, mais son importance en ingénierie est immense. La loi de Fourier et le concept de résistance thermique fournissent des outils analytiques élégants pour les estimations rapides et l'orientation de la conception, tandis que l'équation de la chaleur — résolue numériquement par la méthode des éléments finis — gère toute la complexité des problèmes réels. Que l'objectif soit de minimiser les pertes thermiques d'une enveloppe de bâtiment, de gérer les températures dans une puce, de contrôler la vitesse de refroidissement pendant un processus de trempe, ou de concevoir un système de protection thermique pour la rentrée atmosphérique, la conduction est invariablement partie intégrante du problème. Pour les défis d'ingénierie où la précision est primordiale et les enjeux élevés, des services d'analyse thermique spécialisés apportent la rigueur et les connaissances nécessaires pour transformer la compréhension thermique en conceptions fiables et optimisées.