Model Compte Rendu PDF
Model Compte Rendu PDF
Model Compte Rendu PDF
Université de Ghardaïa
Thème
Simulation de l'effet d'un Obstacle dans un
Ecoulement Laminaire dans une Conduite
Par :
Mr. Allami abdelwahab
Mr.Elkhir younes
Mr. Okbaoui youcef
Nom de prof :
V.M.Chergui
تاثير وضعية العنصر المخضوع و المسافة بن عنصرين. عدة اشكال للعناصر المخوضة قد خضعت للدراسة
النتائج الرقمية المتحصل عليها تثبت.كما تائثير تغير المخضوع على الجريان المائع قد درسوا غي هذا العمل
كما انه لو حظ. بوضح ان كفاءات الحركية و الحرارية متعلقة بعدد و بعد و وضعية و شكل العناصر المخوضة
تطور بنيات دوار نية بجوار العناصر المخوضة هذة البينات تبدوا السبب االساسي لزيادة التبادل الحراري
. كما ان ضياع الطاقة ناتج اساسا عن تعطيل الجريان بسبب هذه العناصر المثبتة على الجدار الداخلي االنبوب
Résumé :
Une simulation numérique du transfert de chaleur et l’analyse de pertes de charge sont présentation
pour des tubulateur de différance forme insérés un conduite de section circulaire. Un flux constant
de chaleur es appliqué du coté externe de la conduite. Les calculs numériques sont réalisées à l’aide
du code COMSOL Multiphysics Déférentes formes de turbulateur sont utilisées. L’effet de la
position de turbulateur et la distance entre deux tubulateur, ainsi que l’effet d’augmentation de la
hauteur sont étudiés .Les résultats numérique obtenus montrent clairement que les performances
dynamique et thermique sont dépendantes du nombre, des dimension, de la position et de la
configuration des tubulateur . On remarque aussi le développement de structures tourbillonnaires,
en aval et en amont des turbulateur, Ces structure semblent être à l’origine de l’augmentation du
transfert de chaleur, ainsi que des pertes de charge dues essentiellement au blocage de l’écoulement
créé par l’emplacement des ces élément
Abstract:
A numerical simulation heat-transfer study and pressure drops analysis are presented for turbulators
of different shapes inserted in a circular pipe. A constant heat flux is applied to the outer side of the
pipe .Numerical are conducted with the COMSOL Multiphysics .The turbulators shapes are
different. The effete of the turbulators position. the distance between two turbulators and the
effect of turbulatos height are studied .The numeriel obtained resulta show clearly that the dynamic
and thermal performances are dependent on the number, the dimensions, the position and the shope
of the tabulators .We can observe a development of vortex structures, behind turbulators , which
seem to be responsible on the heat transfer increase , as pressure drops awed essentially blocking of
the flow by these element
7. Modélisation Numérique par COMSOL MULTIPHYSICS version 4.3, Transfert de
chaleur fluide (eau) et solide dans un conduit présence un obstacle :
Peut être séparée en quatre parties. à gauche, on retrouve le bibliothèque de model où il est possible
de définir des variables et autres paramètres du problème (Parameters), le modèle (Model) où se
retrouvent la géométrie (Geometry), les propriétés des matériaux formant la géométrie, le ou les
modèles de physique s’appliquant au problème étudié et les paramètres de maillage (Mesh). Le
bibliothèque de model comporte aussi le type de problèmes et les paramètres de solveur (Study) et
les options d’affichage et de post traitement des données (Results). La colonne directement à
droite comporte les options relatives à ce qui est sélectionné dans Le bibliothèque de model, par
exemple, les dimensions d’un objet qui vient d’être créé dans (Geometry). C’est aussi à cet endroit
que les paramètres initiaux de la simulation et les modèles physiques nécessaires sont choisis. En
haut à droite, la fenêtre d’affichage graphique (Graphics) permet de visualiser la géométrie, le
maillage ou les résultats. Au haut de cette fenêtre se retrouvent diverses options permettant de
changer le grossissement de l’affichage, l’orientation d’un objet tridimensionnel, etc. Les options
permettant de sélectionner des objets, des domaines, des frontières ou des points se retrouvent aussi
au haut de cette fenêtre.
Finalement, directement au-dessous de la fenêtre d’affichage graphique, il y a une fenêtre
permettant de visualiser les messages d’erreurs, le progrès des simulations, la liste des opérations
effectuées lors du calcul de la solution ainsi que des résultats numériques calculés une fois la
simulation terminée
Les différents éléments de l’interface utilisateur de COMSOL 4.3 a sont présentés à la figure (III-1)
20
Figure (III-1)
7.2. Création du modèle de simulation:
Figure (III-2)
Après choisir le Model Transfert de chaleur et choisir Transfert de chaleur dans les fluides (ht)
voir la figure (III-3)
21
Figure (II
Après cliquer sur flache de terminé on se trouve cette fenêtre suivant voir la figure (III-5) choisir le
22
Figure (III-5)
Et (clic droit de la souris) de Global Définitions et cliquez sur Paramètres voir Figure (III-6).
Lorsque vous définissez une nouvelle variable, il est possible d’associer des unités à cette dernière.
Les unités doivent être mises entre crochets
Diamètre de (D) = 0.25cm
Largeur (L) = 1cm
La masse volumique
La viscosité dynamique
La vitesse d’entrée
Le nombre Reynolds
23
Figure (III-6)
7.3. Géométrie :
Figure (III-7)
24
Après cliquiez sur adroit sur géométrie choisissez sur opération booléennes et choisissez
Différance pour créé un obstacle
Figure (III-8)
25
7.4. Ajout des propriétés du matériel composant le conduit:
Faites apparaître le menu contextuel de Matériels dans le Model transfert de chaleur puis
sélectionnez ouvert Matériel Choisissez liquide choisissez water (eau).
Figure (III-9)
26
7.5. Ajout des conditions aux limites:
Dans transfert de chaleur (ht), ajoutez les conditions aux limites.
Les paramètres des autres conditions sont :
• température d étirée
Figure (III-10)
• le flux de chaleur
Figure (III-11)
27
• la vitesse d’entrée
Figure (III-12)
28
7.6. Maillage:
Cliquez sur Mesh et choisissez comme taille de maillage laissez le maillage.
Cliquez sur construire tout, il est maintenant possible de voir le maillage triangulaire.
Il est possible de définir manuellement les paramètres de maillage.
Figure (III-12)
7.7. Résultats :
29
Figure (III-13)
Figure (III-14)
Vous devez maintenant tracer le profil de température sous l’onglet Résulta et créez une ligne de
coupe en deux dimensions (Line 2D). Définissez la position des points à chaque bout de la ligne de
manière à placer cette dernière à .
Pouvez cliquer sur construire pour vérifier que la ligne de coupe se situe à l’endroit désiré.
Maintenant, créez un groupe de graphiques en une dimension en effectuant un clic droit sur Résulta
et en choisissant 1D Group graphique. Une fois le groupe de graphiques créé, faites un clic droit sur
ce dernier et choisissez l’option Graphique sur line . Une fois le Line Graph sélectionné, sous
l’onglet Data, choisissez la ligne de coupe que vous venez de créer comme Data Set. Sous -Axis
30
Data, entrez dans expression. En cliquant sur Replace Expression ( ), vous avez accès à la listes
de différentes variables calculées par COMSOL si vous désirez tracer des courbes représentant
d’autres variables. Sous -Axis Data, choisissez expression plutôt que Arc Longueur et entrez x
afin d’afficher la courbe en fonction de la position en plutôt qu’en fonction de la longueur de la
ligne de coupe. Il est aussi possible si vous le désirez de changer les titres du graphique et des axes.
La figure ….. Présente la ligne de coupe et la figure 6 présente le graphique du profil de
température dans centre de conduit.
Figure (III-15)
Les figures (III-16) et (III-17) au-dessous présentent les graphes évaluaion de la vitesses et la
température
31
Figure (III-16)
Figure (III-17)
32
Figure (IV-6) : Evolution de la température [K] le long de l’axe de la conduite sans présence de
turbulateurs
Figure (IV-7) : Contours de vitesse [m/s] en présence d’un turbulateur de forme carrée
37
Figure (IV-8) : Contours de température [K] en présence d’un turbulateur de forme carrée
Dans la figure (IV-8), on observe des températures plus basses en aval de l’obstacle en comparant
avec celles en amont, et ça à cause des tourbillons créées par le turbulateur où il semble que des
forts gradients augmentent le transfet de chaleur au tube.
La figure (IV-9) illustre l’augmentation brusque de la vitesse à l’axe de la conduite juste au niveau
parallèle à la position du turbulateur à cause du blocage de l’écoulement à cette position en présence
de l’obstacle carrée.
38
Figure (IV -10) : Evolution de la température [K] le long de l’axe de la conduite en présence
de turbulateur de forme carrée
Par contre dans la figure (IV -10) la température continue son augmentation le long de l’axe de la
conduite au fur à mesure avec l’échauffement de l’écoulement en déplacement suite au flux de
chaleur imposé aux parois, et l’obstacle ne perturbe pas l’allure de cette évolution car il fait partie
de la paroi où la chaleur est dégagée vers l’écoulement.
Ici on fait varier la hauteur H du turbulateur par rapport à la largeur (Diamètre) D de la conduite. Le
turbulateur est placé à la position x=0.475 m, y=0.
La figure (IV-11) montre les différentes hauteurs choisies du turbulateur.
39
H=0.025 m
La figure (IV-14) : Contours de vitesse [m/s] au sein de la conduite pour différentes hauteurs du
turbulateur.
Dans la figure (IV-14), on observe une chute de la vitesse au niveau du tabulateur qui augmente
avec l’augmentation de la hauteur du tabulateur, du essentiellement à l’augmentions de la surface
de blocage de l’écoulement caractérisé pour la hauteur du tabulateur.
La figure (IV -15) présente les contours de témpérature au sein de la conduit pour différentes
hauteurs du turbulateur.
H=0.1 m
H=0.05 m
41
H=0.025 m
La figure (IV-15) : Contours de température [K] au sein de la conduite pour différentes hauteurs
du turbulateur.
Dans ce paragraphe on étudie l’effet de la forme des turbulateurs sur le transfert thermique et les
profils de vitesse. On propose d’étudier trois formes distinctes, à savoir : triangulaire et demi-cercle.
42
Conclusion générale :
Notre étude nous a permis de développer une simulation simple et efficace destinée à améliorer la
conception des échangeurs de chaleur. Elle nous a permis, en outre, d’avoir une connaissance plus
précise des échangeurs thermiques fluide-paroi, des facteurs conditionnant le temps de séjour du
fluide et du profil de températures le long de l’écoulement.
Les objectifs visés par cette modeste contribution ont été atteints dans leur majorité. L’analyse des
comportements dynamiques et thermiques de l’échangeur nous a permis d’aboutir à des résultats
qui pourraient être exploitables en matière de fonctionnement, de dimensionnement et de sensibilité
aux paramètres géométriques constructifs, ainsi qu’aux caractéristiques dynamiques de
l’écoulement du fluide.
Les résultats numériques montrent clairement que les performances dynamiques et thermiques sont
dépendantes du nombre, des dimensions et de la configuration des turbulateurs.
L’étude a montré que le champ dynamique de l’écoulement est fortement perturbé en présence du
turbulateur, le transfert de chaleur et les pertes de charge sont d’autant plus importants que le
nombre et la hauteur des turbulateurs augmentent.
Il apparaît aussi que le champs thermique, lui aussi, d’une grande sensibilité à la présence du
turbulateur. En effet le gradient de la température du fluide à la paroi est plus élevé que celui
correspondant au tube vide sans turbulateur.
Compte tenu des résultats obtenus, on peut déjà souligner l’utilité du modèle numérique élaboré. En
effet, il pourrait être simplement utilisé comme référence et cas test car il ne permet qu’une
approche du coefficient d’échange global métal-fluide. Cependant il nous permettra, d’une part, de
mettre au point d’autres modèles à partir de résultats connus et, d’autre par, de définir le domaine
ou il y a concordance entre ces résultats et ceux des autres modèles.
Il est serait plus judicieux de trouver un compromis entre les optimums dynamiques et thermiques
afin de ne pas altérer la résistance mécanique des éléments constitutifs par un temps de séjour
excessif du fluide, et afin d’augmenter les puissances thermiques tout en réduisant le nombre et les
dimensions des composants des échangeurs à tube.
48
Référence bibliographique
49
[14] : M. LEBID , « Etude expérimental des paramétre contribuant à l’amélioration des
performances dynamique et thermique des échangeur tubulaires » Thése de magistére, Université de
Blida, 1995
[15] : PATRICK CHASSAING, « Tubulateur en mécanique des fluides », ENSSEEIHT,
Collection POLYTECH de L.N.P de Toulouse.
[16] : PH .MARTY «Transfert thermique convectifs », cour .IPU .Génie des system industrielle de
Grenoble 2001.
[17] : ROLAND SCHIESTEL ; « Les écoulement turbulent modélisation et simulation », Edition
HERMES, Paris 1993.
[18]: R.A. WIRTZ and WAYNE MCAULIFFE : «Expermentant mdelling of convecton
upstream form a row electronic packages » Intersociety on thermal phenomena IEEE, CH2798,
Reno, 1990
[19]: THORSTEN STOESSER, FABRICE MATHEY: JOCHEN FROHLICH AND
WOLFCGANG RODI « Les of flow over multiple cubes », ERCOFTAC Bulletin Germany March
2003.
[20]: S.T. CHAN AND D.E. STEVENS: , «An evaluation of two advanced turbulence models for
simulation the flow and dispersion around buildings », Millennium NATO/CCMS intermationnal
technical meeting on air pollution modeling and its application ,UCRL-135557,2000
[21] : YVES JANNOT « Cours de transferts thermiques » 2005
[22] : www.comsol.com
[23] : www.wikipedia.org
50