Modélisation du transfert de chaleur dans un écoulement périodiquement
inversé
1-problematique
c' est un problème couplé de transfert de chaleur et d'écoulement de fluide. Un canal plan épais est chauffé
avec un flux constant depuis le fond (à y=−e ) avec un flux de chaleur constant q comme indiqué sur la
figure ci-jointe. Le canal chauffé est soumis à un écoulement alternatif de vitesse U(t) = umaxsin(2πft). la
vitesse d'écoulement a une moyenne de 0, ce qui signifie que pour la première moitié du cycle, elle atteint
umax et dans la seconde moitié atteint -umax. Dans l'expression de la vitesse, le terme f représente la
fréquence du flux oscillant. Ainsi, la période de temps de l'écoulement est τ=1/f. Cela indique que pour
0<t<τ/2, la frontière à x=0,0<y<d agit comme l'entrée, tandis que pour la seconde moitié, c'est-à-dire
τ>t>τ/2 la frontière à x=L ,d>y>0 est l'entrée. Le plan supérieur (y=d,0<x<L), la face solide gauche
(x=0,0>y>−e) et la face solide droite (x=L,−e<y<0) sont tous isolé. Les domaines solide et fluide sont
couplés par la continuité de la température et du flux à l'interface (y=0).
Dans ce scénario, après un certain temps, l'ensemble du système est supposé atteindre un état stable cyclique. Je
vais maintenant mentionner les équations gouvernantes et les conditions aux limites :
Fluide
Hydrodynamique
∂u∂x+∂v∂y=0 (1)
∂u∂t+u∂u∂x+v∂u∂y=−1ρ∂p∂x+μ(∇2u) (2)
∂v∂t+u∂v∂x+v∂v∂y=−1ρ∂p∂y+μ(∇2v) (3)
No slip and No penetration condition on y=0 and y=d
.Condition de non glissement et de non pénétration sur y=0 et y=d.
thermique
Thermal
ρcp∂T∂t+ρcpu∂T∂x+ρcpv∂T∂y=kf(∇2T)(4)
(The fluid has an inlet temperature T=Ti
. So the boundary condition will be T(x=0)=Ti,∂T∂x|x=L=0 for τ/2>t>0 and T(x=L)=Ti,∂T∂x|
x=0=0 for τ/2<0<τ )
Le fluide a une température d'entrée T=Ti. Donc la condition aux limites sera T(x=0)=Ti,∂T∂x|x=L=0 pour
τ/2>t>0 et T(x=L)=Ti,∂T∂x|x= 0=0 pour τ/2<0<τ
Couplage
A l'interface entre le solide et le fluide, voici ce qui couple le problème :
Continuité de la température
T(x,0)=Ts(x,0),0<x<L(6)
Continuité du Flux
−ks∂Ts∂y|y=0,x=−kf∂T∂y|y=0,x(7)
L'objectif est de résoudre les champs de vitesse et de température dans le solide et le fluide. De plus, le
transfert de flux du solide vers le fluide à l'interface est intéressant.
(2) Comment modéliser le couplage entre le solide et le fluide en utilisant (6) et (7) ? Dans un pas de temps,
que faut-il calculer en premier, la température du solide ou celle du fluide ? Je comprends que les deux
champs doivent être calculés simultanément.
(3) Enfin, puisque le flux s'inverse, comment modéliser le basculement des conditions aux limites à chaque
demi-cycle ?
Certaines valeurs de paramètre L=0.025, d=0.002, e=0.002, k_f=0.614, k_s=390, q=5000, rho=997,
rhos=8960, mu=8.90*10^-4, cp =4178, cps = 385. Ces paramètres représentent l'écoulement de l'eau sur le
cuivre. Un profil de vitesse d'écoulement typique peut être U = 0,3 * sin (2 * pi * 1 * t), qui est un flux de 1
Hz donnant tau = 1.
Un résultat CFD typique obtenu à partir de COMSOL
Pour les mêmes combinaisons de matériaux mais avec un e=4, f=0.5, umax=0.22875, Ti=288 et un apport
calorifique de 1W
qui se traduit par q=40000, j'attache la variation moyenne de la température d'interface (y=0) (ligne intégrée
sur la longueur de x=0 à x=L) avec le temps.
Solution1. Nous utilisons une forme non dimensionnelle d'équations d'échelle d et ts=d2/(kf/(cpρ)). Nous
définissons deux régions reg1, reg2 pour décrire l'écoulement et la température du fluide en conséquence. C'est le
code dans un cas de forme divergente de l'équation de température. Nous utilisons la forme mise à l'échelle du flux
de chaleur qn = q ts/(cp rho)/d, mais la température n'est pas mise à l'échelle
1. Introduction
1-introduction
L’ étude hydrodynamique des écoulements oscillants
2. Equations des écoulements et transfert:
On écrira les équations de Navier stockes avec les conditions aux limites
3-Resolution numeriques : on resout le système d'equations differentielles en utilisant le
logiciel mathematica
4-Analyse et discutions
