Etude des phénomènes du transfert de chaleur[1]

Etude des phénomènes du transfert de chaleur par
simulation numérique
Introduction :
Le transfert de chaleur, qui regroupe la conduction, la convection et le rayonnement, joue un rôle clé dans de
nombreux systèmes thermiques. La simulation numérique de ce Tp permet de modéliser ces phénomènes avec
précision, offrant une alternative économique et flexible aux expériences physiques.
But de la Simulation :
Cette simulation vise à comprendre et optimiser les mécanismes de transfert de chaleur. Elle permet d’analyser
les performances thermiques, de réduire les coûts d’expérimentation, et de concevoir des systèmes
énergétiquement efficaces et durables.
Les étapes a suivre pour la simulation :

 Modélisation géométrique : Création de la poutre carrée dans l'environnement CAO de
SolidWorks.
 Définition des matériaux : Attribution des propriétés thermiques (conductivité thermique) au

matériau constituant la poutre.
 Application des conditions aux limites : Imposition de températures différentes aux deux
extrémités de la poutre pour créer un gradient thermique.
 Maillage : Discrétisation de la géométrie en éléments finis pour la résolution numérique.
 Simulation : Lancement du calcul numérique pour obtenir la distribution de température dans

la poutre.
 Résultats : Visualisation et analyse des résultats, comparaison avec la solution analytique.
Caractéristique de la pièce :
Matériau : matériau spécifique
 Masse volumique : 7870

 Coefficient de chaleur : 5
 Chaleur spécifique : 480
Dimensions : 1000*1000*5000 (mm)
Apres avoir la forme géométrique de la pièce on a mis les températures :

o Température 100°C
o Température 25°C
 On a créé le maillage en choisissant le maillage standard et la taille du maillage

Résultat obtenu après la simulation :
o Distribution de la chaleur :
o Evolution de la température tout au long de la poutre (en °C) :
Interprétation : D'après le résultat qu’on a obtenu, le problème étudié concerne le transfert de

chaleur par conduction dans une poutre de section carrée. La simulation numérique a permis de
visualiser la distribution de température au sein de la poutre et de valider la solution analytique.
D’après La loi de fourrière :
Φ = -λ * S * (dT/dx) (W) avec :
o Φ est le flux thermique (en W/m²)

o λ est la conductivité thermique du matériau (en W/m.K)
o S est la surface de la section (en m²)
o dT/dx est le gradient de température (en K/m)
Calculons T2 , T3 et T4 ?
Données numériques :
 Φ = 75 W/m² (constant)
 S = 1 m²
*Calcul de la température T2:
En appliquant la loi de Fourier, on obtient :
dT = - Φ * dx / (λ * S)
En intégrant cette équation, on trouve :
T2 = T1 + (Φ * Δx) / (-λ * S)
AN :
T2 = 100 + (75 * 1.25) / (-5) = 81.25 °C
* T3 et T4:
Par la meme méthode on va chercher T3 et T4
T3 = T2 + (Φ * Δx) / (-λ * S)
A.N : T3= 81.25 + (75*1.25/-5) = 81.25 - 18.75 = 62.5 °C
T4 = T3 + (ΦΔx) / (-λS)
A.N : T4 = 62.5 + (75 * 1.25) / (-5) = 62.5 - 18.75 = 43.75 °C
PARTIE 2 : simulation numérique de la convection naturelle autour d’un cube

But :
Modélisation et analyse du phénomène de convection naturelle de l’air autour d’un objet solide (le
cube) en utilisant le module Flow Simulation du logiciel SolidWorks.
A. Méthodologie et démarche de l’étude :

1. Modèle géométrique :
2. Méthodologie et démarche de l’étude :
 Ouvrir un nouveau projet appelé <<convection naturelle autour un cube >>dans le module
flow simulation
 Ouvrir assistant en choisissant le système d’unité S.I.(m-kg-s)
 Choisir type d’analyse ‘externe’ avec un régime permanant, cocher la case gravité sans oublier
de vérifier la valeur, direction et sens de force gravitationnelle
 Choisir le type de fluide ‘air’, et le type d’écoulement laminaire uniquement
 Définir les conditions thermiques de paroi « température de paroi = 313.2K »
 Création du maillage volumique avec un rapport d’expansion égale 1,2.
 Après le raffinage des cellules , on lance le calcule

 Création du maillage volumique sue le plan (x,y) avec un décalage de 0.05m
3. Visualisation des résultats de la simulation :
 Distribution de la température de l’air dans le plan de visualisation (x,y)
 Distribution de la masse volumique de l’air dans le plan de visualisation (x,y)
 Champs de vitesse de l’air autour du cube dans le plan de visualisation ( x , y)

 Ligne de courant autour du cube
 Interprétation des résultats :

 D'après le plan de visualisation 1, on remarque que la température est plus élevée
près du cube (Tmax = 313,28 K), et diminue à mesure que l'on s'éloigne du cube
(Tmin = 293,20 K). C'est l'effet de la convection naturelle.
 D'après le plan de visualisation 2, on observe que la densité du fluide est plus faible
près de la surface du cube (d = 1,13 kg/m³), puis augmente progressivement à
mesure que l'on s'éloigne du cube (d = 1,20 kg/m³).
Cette variation de la densité est liée à l'effet de la convection naturelle : près de la
surface du cube, qui est chauffée, l'air se dilate en raison de la température plus
élevée, ce qui réduit sa densité. En revanche, à mesure que l'on s'éloigne du cube,
l'air se refroidit, se contracte et voit sa densité augmenter.
 Les lignes de courant montrent comment l’air se déplace autour et au-dessus du
cube.
Devant le cube : L’air ralentit
Au-dessus : L’air chaud monte rapidement, visible par les lignes montant vers le haut.
La loi de NEWTON :
Q=h× s×(T1−T2)
Nombre de Rayleigh Ra :
Ra=Gr × Pr
g×B×∆T ×h³ μ × Cp
Avec Gr = v²
et Pr = λ
On a :
ρ=1.2kg/m³
λ=0.0318W/m/°C
μ=1.851×10ˉ⁵ Pa/s
Cp=1006.92J/kg/K
Calculs :
Calcul de la vitesse cinématique V :

U 1.851× 10ˉ ⁵
V=ρ = =1.541× 10ˉ ⁵ m²/s
1.2
Calcul de Gr (nombre de Grashof) :
(313−293.2)× 0.1³
Gr = 981× = 2821587.97
(1.541× 10ˉ 5 )²
Calcul de Pr (nombre de Prandtl) :
U × Cp 1.541× 10ˉ 5 ×1006.92

Pr = λ = =0.4878
0.0318
 Ra = 2821587.97× 0. 4 8 78=1377272.6
 Il s’agit d’un écoulement turbulent
Le nombre de Nusselt :
1 1
Nu=0.59× Ra 4 =0.59×1377272.6 4 =13.33
Calcul du coefficient de transfert thermique (h) :

h ×l Nu ×l 13.33× 0.0318
Nu= λ => h= = =4.24 W /m ² ° C
λ 0.1

Etude des phénomènes du transfert de chaleur[1]

Transféré par

Droits d'auteur :

Formats disponibles

Etude des phénomènes du transfert de chaleur[1]

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Etude des phénomènes du transfert de chaleur[1]

Transféré par

Droits d'auteur :

Formats disponibles

Etude des phénomènes du transfert de chaleur par

Les étapes a suivre pour la simulation :

 Définition des matériaux : Attribution des propriétés thermiques (conductivité thermique) au

 Maillage : Discrétisation de la géométrie en éléments finis pour la résolution numérique.

 Simulation : Lancement du calcul numérique pour obtenir la distribution de température dans

 Résultats : Visualisation et analyse des résultats, comparaison avec la solution analytique.

 Masse volumique : 7870

Apres avoir la forme géométrique de la pièce on a mis les températures :

 On a créé le maillage en choisissant le maillage standard et la taille du maillage

o Evolution de la température tout au long de la poutre (en °C) :

Interprétation : D'après le résultat qu’on a obtenu, le problème étudié concerne le transfert de

Φ = -λ * S * (dT/dx) (W) avec :

o Φ est le flux thermique (en W/m²)

*Calcul de la température T2:

En appliquant la loi de Fourier, on obtient :

En intégrant cette équation, on trouve :

T2 = 100 + (75 * 1.25) / (-5) = 81.25 °C

Par la meme méthode on va chercher T3 et T4

A.N : T3= 81.25 + (75*1.25/-5) = 81.25 - 18.75 = 62.5 °C

A.N : T4 = 62.5 + (75 * 1.25) / (-5) = 62.5 - 18.75 = 43.75 °C

PARTIE 2 : simulation numérique de la convection naturelle autour d’un cube

A. Méthodologie et démarche de l’étude :

 Création du maillage volumique avec un rapport d’expansion égale 1,2.

 Après le raffinage des cellules , on lance le calcule

 Distribution de la masse volumique de l’air dans le plan de visualisation (x,y)

 Champs de vitesse de l’air autour du cube dans le plan de visualisation ( x , y)

 Interprétation des résultats :

Calcul de la vitesse cinématique V :

Calcul de Gr (nombre de Grashof) :

Calcul de Pr (nombre de Prandtl) :

U × Cp 1.541× 10ˉ 5 ×1006.92

Calcul du coefficient de transfert thermique (h) :

Vous aimerez peut-être aussi