Modlisation des changes par convection dans un canal rectangulaire muni dailettes

Modlisation des changes par convection dans un canal rectangulaire muni dailettes
Menasria Fouad, Laboratoire de Gnie Mcanique (LGM), Universit Mohamed Khider Biskra Email : menasria. [email protected]

Rsum : Lobjectif de cette tude, consiste tablir un model empirique de calcul du coefficient dchange thermique par convection, lors de lcoulement de lair dans un conduit rectangulaire, dont le plan infrieur est muni de chicanes de formes rectangulaires, disposes en quinconce. Par la mthode de lanalyse dimensionnelle, on a pu mettre en relation tout les paramtres physiques, thermophysiques et les caractristiques gomtriques des chicanes, et en fonction du rgime dcoulement on estime un coefficient dchange par convection dans lespace dcoulement du fluide. Mots cls : corrlation, coefficient dchange, convection, rugosit artificielle, chicane, coulement
turbulent, capteur solaire plan air. Nomenclature

: Pas entre deux ranges de[m] : Longueur de la conduite [m]

= = = =

. . . .

: Nombre de Reynolds : Diamtre hydraulique de la conduite : Nombre de Prandtl utile[m] = [m2] : Temprature lentre du capteur [C] : Temprature la sortie du capteur [C] :Facteur de Colburn : Flux solaire reu par la face suprieure du capteur [W/m] = [K] = [Pas]
1

: Nombre de Brinkman . : Nombre de Peclet : Nombre dEckert

: Surface de labsorbeur [m2]

: Section de passage de lair dans la veine

=
=
=
.

: Nombre de Nusselt

: Chaleur massique de lair [J/Kg.K] : Vitesse moyenne du fluide[m/s] : Longueur de chicane[m] : Largeur du capteur[m] : Pas entre deux chicanes[m]

: Diffrence de temprature

: Diffrence de pression

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

: Flux utile [W] : Temprature de parois [K] : Temprature du fluide : Coefficient dchange thermique par convection [W/m2.K] : Dbit massique de lair [Kg/S] : Surface transversale minimale du capteur [m2] : Temprature dentre du fluide [K] : Temprature de la sortie du fluide [K] : Longueur dune tranche fictive[m]

: Temprature du fluide moyenne (exprimentale)[K] : Primtre mouill[m] : Nombre des chicanes par range : Coefficient de conductibilit thermique de lair [w/m.K] :Coefficient de conductibilit thermique de lair [w/m.K] : Masse volumique de lair [Kg/m3] : Viscosit dynamique [Kg.m-1.s-1] : Hauteur absolue ou rugosit absolue dune chicane[m]

I- Introduction Les capteurs solaires plans air, qui sont lobjet de notre tude transforment lnergie solaire incidente solaire en nergie thermique extraite par lair en coulement dans le capteur. Cette nergie et ensuite utilise dans diverses application, savoir ; le chauffage la rfrigration, le schage, etc. Ltude du comportement thermique de notre capteur relativement aux configurations

utilise, est une tude parmi dautres dj labores dans divers laboratoire de recherche. Lobjectif principal est de mettre en vidence les paramtres de performances, noter loptimisation des changes convectifs dans le conduit utile du capteur, et la modlisation de ces changes en fonction de plusieurs variables, notamment les paramtres physiques, thermophysiques et les caractristiques gomtriques des chicanes. Plusieurs types de capteurs air on t construits et tests. Il se diffrent par leurs configurations, du mode dcoulement de lair, du type et du nombre de vitrage utilis, par la nature de labsorbeur ainsi que par le mode dadjonction de tles minces places en ligne ou en quinconce entre labsorbeur et lisolant figure (1).

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

Chicanes en ligne Chicane

Chicanes en quinconce

Pas entre deux ranges

Pas entre deux chicanes

Longueur d'une chicane

Fig. (1) Disposition en range et en quinconces des chicanes dans la veine dair mobile II- Description du dispositif exprimental Pour le droulement du protocole exprimental, un capteur solaire plan air simple passe entre labsorbeur et la plaque absorbante en acier galvanis a t ralis au dpartement de gnie mcanique de luniversit Mohamed Khider de Biskra fig.(2).

Fig.(2) : capteur solaire plan air expos au soleil (photo dgm.Univ de Biskra 2007) En vue damliorer les performances thermiques du capteur solaire, on a introduit plusieurs ranges dailettes mtalliques trs minces dpaisseur 0,4 mm confectionnes dacier galvanis, de forme rectangulaire, dont la partie suprieure est incline. Les chicanes sont
3

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

implantes sur le plan dcoulement infrieur, travers par un coulement dair stationnaire. Le plan dcoulement infrieur est aussi dacier galvanis de dimension gale de . . , lpaisseur de la veine dair mobile est de . entre le plan infrieur et labsorbeur. Lentre dair dans le capteur se fait travers des orifices circulaires localiss et centres larrire du capteur, de telle sorte assurer au moyen dune pompe muni dun variateur de vitesse, laspiration en permanence dun flux dair. En raison de linsuffisance de lchange thermique ralis dans le capteur solaire plan air entre le fluide et labsorbeur lutilisateur est amen apporter des amliorations pour de meilleures performances ou une meilleure efficacit thermique. La figure (3) suivante reprsente les lments du collecteur solaire plan air tudi :

, qui correspond lespace compris

Fig. (3) Schma dun capteur solaire plan air III- Bilan thermique Le bilan thermique dun insolateur au niveau de la veine dair mobile rapport lunit de surface, se rsume ainsi par la quantit dnergie utile dans notre cas il sagit de lair. On considre la tranche" " voir figure (8-a), La temprature est suppose uniforme au emporte par le fluide caloporteur,

niveau de la section de la veine dair mobile de linsolateur dmuni des chicanes.

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

a)capteur solaire muni de chicanes capteur chicanes. b) Tranche fictive " contrle. Fig. (4) : schma simplifier dune section longitudinale

" du volume de

La conservation de lnergie dans cette tranche nous donne avec ltant de largeur d capteur du . . Aprs simplification et intgration on obtient :
. .

(1) (2)

Supposons aussi que la temprature du fluide varie linairement le long du capteur, la valeur moyenne de la temprature est par consquent tant gale la moyenne arithmtique entre lentre et la sortie, par consquent : , (3) Ou dune manire gnrale :
/

(4)

Donc on peut calculer la temprature moyenne du fluide daprs la relation suivante : (5) A partir des deux dernires quations, on dduit : (6) IV. Effet des obstacles sur les changes convectifs Lintroduction dobstacles dans la veine dair dynamique ncessite un choix judicieux du mode de fixation, il est possible quils soient fixs sur lisolant ou sous labsorbeur, ou les , combinaisons ensemble. Dans tous les cas les performances thermiques sont nettement amliors, du fait que la surface dchange totale mise en jeu se trouve augmente et la s, section de passage offerte devant le passage du fluide qui est par consquent rduite.
5

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

Dans la veine dair mobile, le calcul du nombre de Reynolds (Re) est obtenu par le calcul de la vitesse maximale de lair ( ) correspondant la section de passage minimale ( ), qui

est exprim par la relation suivante :

.

(7)

Ou,

est le diamtre hydraulique, dfinit comme tant le rapport de 4 fois la section de

passage sur le primtre mouill, autrement : =

.

(8)

Dans le cas dune veine rectangulaire muni de chicanes figure (4), le diamtre hydraulique est calcul partir du schma quivalent suivant, ou on tient compte de la section obture transversale par une range constitue de n obstacles de chicanes disposes devant le passage du fluide.

a(disposition en ligne)

b(disposition en quinconce 2ieme ranges)

Fig. (4) Coupe transversale dun capteur muni de chicanes En se rfrant au schma de la figure (4a) on aurait : =
[ . . . . ]

(9)

Et en se rfrant au schma de la figure (4b), on aurait : =

[ . ( . . ). ]

(10) = ,

Sachant que la section minimale offerte est Par consquent la vitesse maximale = = est :

(11)

Le nombre de Reynolds serait alors : Soit dans la veine o les chicanes sont disposes en ligne est :
6

Modlisation des changes par convection dans un canal rectangulaire muni dailettes
.

. ]

(12)

Soit dans la veine o les chicanes sont disposes en quinconce est : convectif dfinit par : =
[ . ( ). ]

(13)

Le coefficient dchange thermique h entre labsorbeur et lair dpend notamment du nombre Reynolds. On dduit alors que laugmentation de la rugosit absolue entrane celle de Re

et h. La section minimale A est fonction de la forme des obstacles, de leurs dimensions et de leur disposition. Dans la veine dair mobile, la convection force dans le cas des conduites rectangulaires, lchange entre labsorbeur et le fluide caloporteur est caractris par un coefficient dchange convectif dfinit par :
_

(14)

Dautre part, le coefficient dchange entre le fluide caloporteur et la plaque infrieure isolante est pratiquement peu diffrent de celui entre labsorbeur et lair convectif du ct isolant
_ _

, gnralement

considr de mme valeur pour simplifier le calcul. Par consquent, le coefficient dchange et labsorbeur : =
_ _ _

est donn par les mmes corrlations que celui entre le fluide (15) est considr comme

Pour le cas dun capteur muni de rugosit artificielle (chicanes),

le coefficient dchange moyen, calcul par la relation prcdente. Le nombre de Nusselt est dtermin par la corrlation de Kays (pour un absorbeur avec chicane) suivante [1] : = 0,0158.
.

Pour (

) >10

(16)

Pour les capteurs absorbeurs munis de chicanes, on utilise le facteur de Colburn J . Pour 120 < capteurs solaires V- Etablissement dun modle de calcul du coefficient dchange On essayera dtablir des relations qui permettent de mettre en vidence des modles de calcul du coefficient dchange par convection , dans le conduit rectangulaire dun < 10000, avec = 0,66 0,7 pour lair comme fluide caloporteur dans les

capteur solaire, dont le plan suprieur est parfaitement rayonn par un flux solaire incident. Le conduit utile est muni de plusieurs ranges de chicanes, qui sont disposes en quinconce. Ces chicanes sont insres directement sur le plan infrieur dcoulement de la veine fluide.
7

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

Avant de procds llaboration du modle de calcul du coefficient convectif, qui est dans notre cas bas sur la mthode de lanalyse dimensionnelle, Il fallait pralablement dfinir linventaire de tous les paramtres dont dpend la densit de flux de chaleur qui est donne par lexpression ( forme : = ( , , . T), ce qui nous a permis dtablir une relation gnrale de la , , , , , , , , , , , ) (17)

Daprs le thorme de

(de Vashy-Bukingham), on ne peut voir que 09 groupements

indpendants, adoptons alors lcriture suivante, avec ( , ) qui sont respectivement la longueur et la largeur du canal rectangulaire constituant le conduit du capteur, par consquent on obtient :
.

= . .

. .

(18)

Lcriture de lquation prcdente compte tenu les dimensions fondamentales, aprs dveloppement et identification on obtient un systme de 04 quations dont la rsolution aboutit une expression gnrale de la forme :
. . . . . . . .

.(

.(

(19)

On introduisant les nombres adimensionnels, lexpression (20) devient : =

. . .

(20)

Ou bien rcrite sous la forme : =

. . .

(21)

Pour dterminer les corrlations exprimant le coefficient dchange , il faut identifier les valeurs exacts de ( , , , , , , , ). Une srie de mesures exprimentales ont t ralises au dpartement de Gnie Mcanique de luniversit de Biskra, effectu sur un capteur solaire plan air, muni de plusieurs ranges de chicanes de forme rectangulaire de longueur 7 cm, disposes en quinconce. Les

paramtriques gomtriques considres sont montrs daprs le tableau (1) suivant :

Configurations Premier cas Paramtres gomtriques Longueur dune chicane Pas entre deux chicanes Pas entre deux ranges Deuxime cas

=7

= 10

=4

=7

= 20

=4

Tableau.(1) : Paramtres gomtriques des chicanes tudies

8

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

Dans cette tude exprimentale, on sest intress deux configurations ou le pas entre deux ranges successives est respectivement 10 cm dans le premier cas et 20 cm dans le second, par contre les autres paramtres sont fixs constant. En considrant des rgimes dcoulements variables, on a pu obtenir les rsultats suivants : A) Rgime laminaire : Pour le premier cas, ( Pe r = 10 cm), lexpression empirique du coefficient dchange est :
. . .
.

Pour le premier cas, ( Pe r = 20 cm), le coefficient dchange convectif est :

= . . .
.

B) En rgime turbulent : Pour le premier cas, ( Pe r = 10 cm) :

= . . .
.

Pour le premier cas, ( Pe r = 20 cm) :

= . . .
.

a)Cas (

= 10

b) Cas (

= 10

Fig. (5) : Evolution de coefficient dchange en fonction de nombre de Reynolds (Re<2100)

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

a)Cas (

10

b) Cas (

10

Fig. (6) : Evolution de coefficient dchange en fonction de nombre de Reynolds (Re>2100)

a)Cas ( 10 ) b) Cas ( 10 ) Fig. (7) : Evolution de coefficient dchange en fonction de nombre de Prandtl (Re<2100)

a)Cas ( 10 ) b) Cas ( 10 ) Fig. (8) : Evolution de coefficient dchange en fonction de nombre de Prandtl (Re>2100)

10

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

a)Cas ( 10 ) b) Cas ( 10 ) Fig. (9) : Evolution de nombre de Nusselt en fonction de nombre de Reynolds (Re<2100)

a)Cas (

10

b) Cas (

10

Fig. (10) Evolution de nombre de Nusselt en fonction de nombre de Reynolds (Re>2100)

a)Cas (

10

b) Cas (

10

Fig. (11) Evolution de facteur de Colburn en fonction de nombre de Reynolds (Re<2100)

11

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

a)Cas (

10

b) Cas (

10

Fig. (12) Evolution de facteur de Colburn en fonction de nombre de Reynolds (Re>2100) IV- Discussion et Interprtation des Rsultats : En rgime dcoulement laminaire ( < 2100), concernant linfluence du nombre de

Reynolds sur le coefficient dchange thermique convectif, on constate que le coefficient dchange thermique ( ) varie dans lintervalle de 4 45 [W/m2.K], pour le premier cas ou le pas entre range Pe r = 10 cm, comme le montre la Fig. (5-a), dautre part il varie de 4 35 [W/m2], dans le deuxime cas, ou le pas entre range est le double Pe r = 20 cm , voir Fig. (5b), ce qui montre leffet significatif du pas entre range, ou on pourrait amliorer lchange convectif lorsque le pas entre range diminue. En Rgime turbulent ( > 2100), lcart est plus significatif, le coefficient dchange thermique varie de 40 70 [W/m2.K] dans le premier cas Fig. (6-a), par contre il ne varie que de 35 55 [W/m2.K] dans le deuxime cas Fig. (6-b). Lanalyse des rsultats numriques, montre que les amliorations apports au coefficient dchange convectif sont plus importantes et devient presque le double en rgime turbulent dans les deux cas tudies, notamment en prsence de la configuration des chicanes, ou le pas entre range est de 10 cm. Quand au nombre de Nusselt, les figures (9 et 10)(a, b), en rgime laminaire montre que sa variation en fonction du nombre de Reynolds, Ou le Nusselt varie de 05 28 pour le premier cas, et de 2 24 dans le deuxime cas. En rgime turbulent, le nombre de Nusselt varie de 22 42 dans le cas ou le pas entre range est de 10 cm, alors quil ne varie que de 20 35, lorsque le pas entre range est de 20 cm.
12

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

Concernant lvolution du coefficient dchange en fonction du nombre du Prandtl figures (7, 8)-(a-b), ou on constate que lchange thermique est pratiquement dpendant du nombre de Prandtl, leffet de laugmentation du coefficient dchange commence tre ressenti pour des nombre de Prandtl comprise entre 0,69 et 0,75 pour lair, on remarque qu partir de cet intervalle, la qualit de lchange thermique est nettement meilleure pour les grandes valeurs de Prandtl. VII- Simulation dynamique et thermique de lcoulement de lair : Dans la prsente tude, on considre la simulation numrique du transfert thermique par convection force tridimensionnelle, dans une veine rectangulaire horizontale muni des chicanes (particulirement dans les capteurs solaire plan), traverse par lair, dont la vitesse varie suivant le rgime dcoulement soit laminaire ou turbulent. La simulation numrique des paramtres considrs, a t ralise par Fluent. Ce programme est bas sur le travail dcrit par Patankar. Le recours au logiciel Fluent, est pour de multiples raisons, particulirement la possibilit de visualiser le contour des lignes de courants, de temprature et de la dpression entre lamont et laval du conduit occup par des obstacles La dfinition de la gomtrie et la gnration du maillage ont t ralises laide du code gnrateur de mailles Gambit. VII-1. Modlisation Mathmatique : Les principes physiques qui rgissent ce problme, sont modliss par les quations de conservation de la masse, des quantits de mouvements et de lnergie. Equation de continuit : 0 Equations de transport de quantit de mouvement : Suivant laxe y : Suivant laxe z- : (25) (24) (22)

Suivant laxe x- le mouvement des particules fluides est rgit par (23)

13

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

Avec :
(

=,

= 0,

= 0,

Equation de conservation de lnergie :

)

)=

(26)

La mthode des volumes finis, est base sur lalgorithme SIMPLE [2], utilise pour rsoudre le systme de quations (22), (23), (24), (25), (26), on impose les conditions aux limites suivantes 1. A lentre : vitesse du fluide lentre (velocity_inlet) 2. A la sortie : Pression la sortie (pressure_outlet) (101380Pas) 3. Au niveau des autres parois : wall

Fig. (13) prsentation dune chicane en3D

Fig. (14) Prsentation du plan dcoulement muni de plusieurs ranges de chicanes Lespacement (maillage) . . Type de maillage Ttradrique Ttradrique Nombre de mailles 465252 1162694

Cas1 Cas2

Fig. (15) : visualisation des tourbillons Axes verticaux en prsence des chicanes disposes en quinconce.
14

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

Fig. (16) Tourbillons Axes horizontaux (lignes de courant)

Fig. (17) : contour de la temprature en prsence des chicanes.

Fig. (18) contour de la dpression de lamont jusqu laval du conduit en prsence des chicanes La Fig. (17) prsente le champs de temprature le long du plan dcoulement, Les caractristiques 3D sont clairement visibles sur les champs de temprature, en effet la temprature augmente sensiblement lorsque les particules fluides se trouvent juste

proximit immdiate des chicanes. Dautre part la Fig. (18) montre la dpression le long du conduit, ce qui implique quil y a une perte de charge entre lentre et la sortie du conduit utile. Cette perte de charge est induite par les obstacles interposs devant le passage du fluide (chicanes).Dautre on sest intress au comportement dynamique et thermique de la lame fluide en coulement le long dune section longitudinale du conduit, savoir que le plan de base considr avait les coordonnes spatiales suivantes :
15

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

X[m] 0 1.4 1.4 0 Y[m] 0 0 0.025 0.025 Z[m] 0.5 0.44 0.44 0.44 La Figure Fig. (19), illustre les contours des profils de tempratures dans le plan orthogonal, et longitudinalement par rapport au conduit utile. On constate que dans la rgion situe juste aprs la chicane la temprature augmente, par contre la figure Fig. (27) montre que la pression diminue.

Fig. (19) prsentation de la surface a tudie

Fig. (20) contour des vitesses au niveau de la section considre

Fig. (21) contour des profiles de tempratures au niveau de la section tudie

Fig. (22) contour de la pression au niveau de la section tudie

Fig. (23) Evolution du coefficient dchange convectif le long du conduit = , = ,

Fig. (24) Evolution du nombre de Nusselt le long du conduit = < 2100

16

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

Fig. (25) Evolution de la temprature dans Fig. (26) Evolution de la pression le long la veine fluide, en prsence des chicanes de la veine fluide, en prsence des chicanes disposes en quinconce disposes en quinconce = , = , = < 2100 Les pertes de charge gnres par FLUENT, sont reprsentes daprs figure Fig. (25) et valent environ (29.56Pa), ou on note que la valeur des perte de charge calcule numriquement, concide avec celle calcule par les approches semi-empiriques dvelopps au laboratoire. La discontinuit des graphes implique la discontinuit dcoulement cause de lexistence des obstacles. Remarque : Comme, le nombre de Reynolds est lev, cela prouve que lcoulement nest pas du tout laminaire. On a utilis par consquent un modle de turbulence. Contrairement au cas laminaire, les solutions sont plus dpendantes du maillage, et une attention particulire doit tre porte sur la vrification a posteriori de la finesse du maillage prs des parois, afin de sassurer que tous les phnomnes sont capturs. Dans cette tude on a utilis le modle de turbulence le plus simple qui existe dite, modle de VIII- Conclusion gnrale Le recours aux rugosits artificielles pour optimiser les surfaces dchanges, est en vue de remdier aux proprits thermophysiques mdiocres de lair, par rapport aux autres fluides caloporteurs. Les chicanes utilises dans ce travail, sont de formes rectangulaires avec une partie incline de longueur 07 cm, disposes en quinconce, le pas entre deux chicanes est de 04 cm, par contre on a considr deux pas entre les ranges qui sont de 10 cm et 20 cm. La configuration en quinconce offre la possibilit davoir des rtrcissements et des largissements brusques qui se manifestent par la formation dune structure tourbillonnaire axe verticale. Dautre
17

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

part, les parties suprieures des chicanes qui sont inclines de 60, permettent davoir un espace rduit vis vis du plan dcoulement suprieur (absorbeur), simplique par la

formation des tourbillons axes horizontaux, qui se dveloppent laval des arrtes suprieures des ailettes. Lalternance conjugue de ce doublet de tourbillons axe verticaux et horizontaux assure un coulement dsordonn et par consquent un change convectif trs intense. En utilisant le thorme de lanalyse dimensionnelle, une corrlation a t dveloppe, qui permettrai de mettre en relation touts les paramtres physiques, thermophysiques et gomtriques. Un programme informatique a t dvelopp, crit sous lenvironnement Matlab 6 .5, ce qui nous a permis dobtenir plusieurs informations, en fonction du rgime dcoulement laminaire ou turbulent, et du pas relatif sparant deux ranges successives. Les illustrations graphiques, qui montrent lvolution du coefficient dchange en fonction des nombres de Reynolds et de Prandtl, ainsi que lvolution du Nusselt et du facteur de Colburn en fonction de Reynolds, respectivement en comparaison avec les rsultats

exprimentaux, ont permis de juger la validit des corrlations avec la thorie. Comme rsultats important signaler, concernant linfluence du nombre de Reynolds sur le coefficient dchange thermique convectif, nous pouvons constater quen rgime laminaire, le coefficient de transfert convectif (h) varie de 4 45 [W/m2.K], dans le cas du pas entre range , alors quil ne varie que de 4 35 [W/m2.K], dans le cas, ou , ce qui montre leffet significatif du pas entre range En rgime turbulent ( [W/m2.K] dans le second. Les rsultats obtenus par la simulation numrique assiste par Fluent 6.1, nous ont permis de comprendre et de visualiser de faon dtaille le comportement dynamique et thermique de lcoulement de lair, particulirement proximit des chicanes. Les profils de distribution de la vitesse axiale, montrent une zone de recirculation relativement intense aprs les rtrcissements et largissements brusques ; notamment en prsence de la disposition en quinconce des ailettes. IIV- Bibliographie : [1]N.Chouchane, A.Moummi, N.Moummi, K.Oues. Modelisation des pertes de charge
18

> 2100), lcart est plus significatif, le coefficient dchange

thermique varie de 40 70 [W/m2.K] dans le premier, par contre il ne varie que de 35 55

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

dans un conduit rectangulaire garni de Rugosites Artificielles avec une partie superieure Incline. Sminaire international de la physqiue nergtique SIPE'7, Bechar du 03 au 05 Octobre 2004. [2]A.Abene,V.Dubois,M. Le Ray,

[8]Dong Junqi , Chen Jiangping, Chen Zhijiu, Zhou Yimin, Zhang Wenfeng. Heat transfer and pressure drop correlations for the wavy fin and flat tube heat exchangers. Applied Thermal Engineering 27. 2066 2073, (2007) [9]B.B.Prasad et J.S.Saini, Effect of

A.Ouagued. Study of a solar air flat plate collector : use of obstacles and application for the drying of grape, Journal of Food Engineering,65 ,15-22 (2004). [3]Mohammed.Alia. Modelisation des

artificial roughness on heat transfer and friction factor in solar air-heaters, Solar Energy 41, 6. 505-560, (1980). [10]R.M.Manglik and A.E.Bergles,Heat Transfer and Pressure Drop Correlations for the Rectangular Offset Strip Fin Compact Heat Exchanger,Experimental Thermal and Fluid Science,10.PP.171-180 (1995). [11]A.R.Wieting,Empirical Correlations for Heat Transfer and Flow Friction Characteristics of Rectangular Offset-Fin PlateFin Heat Exchangers,Trans.ASME J.Heat Transfer 97,488-490 (AUGUST1975). [12]H.M.Joshi and R.L.Webb, Heat transfer and friction in the offset strip-fin heat exchanger.Int.J.heat Mass Transfer 30,69-84(1987). [13]W.M.Kays,and A.L.London,Heat transfer and flow friction characteristics of some compact heat exchanger surfaces,J.Engng Power January,2734(1960). [14]Sen Hu and Keith E.Herold,Prandtl number effect on offset fin heat exchanger performance:experimental results.Int.J.heat Mass Transfer 38,1053-1061. [15]A.E.Bergles Heat Transfer Enhancement The Encouragement and Accommodation of High Heat Fluxes
19

capteurs solaire air distines au schage (convection force), mmoire de magistere, Universit de Tiaret, (2005), Algrie [4]Ben Slama Romdhane. The air solar collectors : Comparative study, introduction of baffles to favor the heat transfer, Renewable Energy, 81 139-149 (2007). [5]Sabri Youcef-Ali. Etude numerique et exprimentale des schoirs solaires indirects convection force : Application la pomme de terre, these de doctorat,

Universit de Valenciennes (2001 France. [6]N.Chouchane, Modlisation des pertes de charge dans un conduit rectangulaire garni de rugosits artificielles avec une partie suprieure incline, Mmoire de Magister, Universit de Biskra, Algrie 2003. [7]F.Mnasria, Modlisation des changes par convection dans un canal rectangulaire muni dailettes, Mmoire de Magister, Universit de Biskra, Algrie 2003.

Modlisation des changes par convection dans un canal rectangulaire muni dailettes

Trans.ASME J.Heat 19(FEBRUARY 1997).

Transfer

119,9-

convection force :Application la pomme de terre, these de doctorat (2001). [23]A.Ahmed-Zad,A.Moulla,M.S.Hantala et J.Y.Desmons,Amlioration des Performances des Capteurs Solaires Plan air:Application au Schage de lOignon Jaune et du Hareng.Rev.energ.Ren.4.6978(2001).(shema perte de charge). [24]Yves (2008). Jannot, Transfert Thermique

[16]FIEBIG,M ,Vortice Tools to Influence Heat Transfer, Recent Developments.Proceedings of the 2nd European Thermal-Sciences and 14th UTT National Heat Transfer Conference,Edizoioni FTS,Pisa,Italy 1.P.4156 (1996). [17]Amon,C.H ,Spectral Element Fourier Method for Transtional Flows in Complex Geometries,AlAA Journal,6.No.1,PP4248(1993). [18]N.Chouchane, A.Moummi, N.Moummi, K.Oues, Modelisation des Pertes de Charge dans un Cconduit Rectangulaire Garni de Rugosites Artificielles avec une Partie Superieure Incline. [19]A.Aben,V.Dubois,M. Le Ray,A.Ouagued,Study of a solar air flrt plate collector : use of obstacles and application for the drying of grape,Journal of Food Engineering,65,15-22(2004). [20]Mohammed.Alia,Modelisation des Capteurs Solaire air distines au schage (convection force)these de Magistre (2005). [21]Ben Slama Romdhane,The air solar collectors : Comparative study,introduction of baffles to favor the heat transfer,J Solar Energy 81,139-149 (2007). [22]Sabri Youcef-Ali,Etude numerique et exprimentale des schoirs solaires indirects

[25]Abdelmalek Zitouni,Rabie Ait Chelalet,les capteurs solaire, mmoire dingniorat,Jijel(2008) [26]Derouiche Younes,Medjouri Mohamed Ali,Etude de levolution de la difference temperature en fonction des parametres geometriques et physiques dans un capteur solaire plan a air,mmoire dingniorat,Biskra(2007). [27]Rachid Saim Said Abboudi Boumedine Benyoucef Ahmed Azzi,Analyse numriquede la convection force turbilente dans les tubes muni des chicanes transversales,13 Journes internationales de thermique,Albi,France(2007) . [28]Nazih Marzouqui,T.P.Fluent,Cours Mcanique des fluides (2006). [29]Balima Olivier, Fontfryde Thomas, George Erwin, Marinhas Sandrine, Rezgui Azdine,Etude de la convection naturelle dans une cavit carre en 2D et 3D sous Fluent et Gambit

20