Simulations

Présentation du logiciel Jadim

Jadim est un logiciel de recherche développé par le groupe Interface de l'Institut de Mécanique des Fluides de Toulouse. Ce code permet l'étude des paramètres intervenant au sein d'écoulements diphasiques. Il résout les équations de Navier-Stokes en 3D dans le cadre de l'hypothèse de fluides incompressibles et pour des cas instationnaires.

Jadim utilise un maillage structuré et permet d'effectuer les simulations en DNS ou en LES

Hypothèses simplificatrices

Nous n'étudierons pas directement les phénomènes de rupture et de coalescence. En effet, l'étude de ces phénomène demande le suivi lagrangien des gouttes. Nous nous restreindrons dans le cadre de cette étude à l'étude eulerienne de l'écoulement.  Les gouttes ne sont pas ici modélisées et l'écoulement simulé est donc monophasique. Il s'agit donc ici de prendre les propriétés physique d'un fluide équivalent. Nous ne prendrons en compte aucun effet thermique. 

​Propriétés physiques

Les propriétés du fluide sont prises égales à celles de l'eau, c'est à dire : 

µf = 10-3

$\rho_{f}$= 103

​Contraintes

Il existe un certain nombre de contraintes inhérentes à l'outil numérique. Notamment, le respect des critères de stabilité entraîne une contrainte sur les dimensions des mailles et sur le pas de temps, ce qui peut entraîner des temps de calcul très important. 

Nous garderons également à l'esprit que le maillage doit permettre de prendre en compte toutes les échelles impliquées dans les phénomènes que l'on souhaite observer. Un compromis entre précision et temps de calcul doit donc être trouvé.

 

Première simulation - jet centré

Présentation

Présentation de la simulation

Cette simulation est une première approche du problème nous permettant de savoir s'il est possible de simuler directement l'écoulement dans les conditions industrielles et dans le cadre des hypothèses d'un écoulement monophasique incompressible. Ce premier cas nous a permis de poser par la suite d'autres hypothèses simplificatrices essentielles. Le calcul est effectué en 2D axisymétrique et tourne en DNS.

Une première simulation aux conditions industrielles demande cependant plusieurs étapes intermédiaires. Un premier cas pour V = 1m.s-1 nous permet de vérifier que le maillage généré permet au code de tourner. Il s'agit ensuite d'augmenter la vitesse jusqu'à la vitesse réelle d'écoulement. Cependant, l'essentiel restera d'obtenir un cas qui ne diverge pas, ce qui peut demander certains ajustements.

Géométrie

On gardera pour cette étude les dimensions du système industriel. La géométrie du système est la suivante : 

Géométrie industrielle

On ne modélisera cependant que la cavité dans laquelle s'écoule le jet au niveau de l'élargissement brusque

Domaine d'écoulement du jet

On imposera un profil de type Poiseuille en entrée du domaine plutôt que de simuler l'écoulement en amont.

Maillage

Le maillage utilisé pour cette simulation est raffiné en arctangente au niveau de la paroi d'entrée et autour du jet. La zone de jet, de même épaisseur que l'entrée, est maillée de manière uniforme de façon à obtenir 50 mailles dans les 10 micromètres de l'entrée.  Ce raffinage doit permettre de capter les forts gradients de vitesse en entrée causées par les contraintes de cisaillement et au sein du jet où la vitesse décroît rapidement et  dans lequel on observe une forte intensité de la turbulence.

 

Maillage du domaine

Maillage de l'entrée

Résultats et commentaires

Validation du maillage

Afin de valider le maillage que nous avons généré à l'aide du logiciel matlab, nous avons procédé à une première simulation très brèves pour une vitesse d'entrée de 1m.s-1. Nous avons ainsi pu vérifier le bon comportement du fluide aux limites du domaine, notamment en observant la présence d'un Poiseuille en entrée.

 

Champ de vitesse - test des conditions limites

 

   

  Vitesse en entrée du système - t=0.1 ms     Vitesse en sortie de système - t=0.1 ms

Le profil observé en entrée est une extrapolation entre le profil Poiseuille en première maille (non-observable) et la maille calculée suivantes. Ainsi on n'observe pas directement un profil de Poiseuille.

A la paroi limite nord du domaine, on observe bien une vitesse nulle à la paroi et une vitesse sortante à la sortie.

Test aux conditions industrielles

Ainsi assurés de la validité du maillage, nous avons pu lancer une seconde simulation pour une vitesse de 200m.s-1 correspondant à la vitesse d'entrée du fluide dans le système industriel. Le pas de temps est extrêmement faible, de l'ordre de 10-10s. Le calcul tourne alors 1 060 000 itérations (0.104 ms) avant de diverger.

Champ de vitesse dans le domaine - le calcul diverge

Dans le but de faire converger le calcul, nous avons donc décidé de diminuer la vitesse d'entrée du fluide. Cependant, pour conserver la physique de l'écoulement qui nous intéresse, il nous a fallu conserver le même nombre de Reynolds pour toutes les simulations.

Conservation du Reynolds

Afin de diminué la vitesse d'entrée du fluide tout en conservant le Reynolds de l'écoulement, nous avons choisi de diviser la vitesse et la viscosité du fluide par un même facteur. Ainsi nous avons choisi arbitrairement de diviser ces deux quantités par 10. Toutes les simulations réalisées par la suite l'ont été pour une vitesse d'entrée de 20m.s-1

Par ailleurs, on observe que le calcul à V=200 m.s-1 diverge lorsque les vitesses deviennent importantes en sortie. Dans cette zone, les vitesses passent d'une valeur nulle en paroi à une valeur très importante à "l'angle" de la conduite de sortie.  Dans la suite de notre travail, nous avons donc cherché à réduire l'importance de cette contrainte, lourde pour un code CFD.

Seconde simulation - jet décalé

Présentation

Démarche

Afin d'éviter d'obtenir de trop fort gradients de vitesse en sortie, nous avons décidé de décaler la position du jet de manière à éloigner la zone impactante de la sortie.

Maillage

Nous avons raffiné d'avantage le maillage au niveau des parois et notamment au niveau de la paroi opposée à l'entrée où le jet vient impacter. Le nombre de mailles selon y a été changé pour passer de 64 à 128 mailles. Ainsi nous obtenons le maillage suivant :

Maillage raffiné aux parois - jet décalé

 

 

Résultats et commentaires

Simulation du jet décalé

Le fait d'avoir décalé le jet plus loin de la sortie et d'avoir modifié le maillage nous a permis d'obtenir un écoulement beaucoup plus établi que précédement. Le calcul diverge au bout de 228 000 itérations seulement, mais le pas de temps calculé par Jadim oscille autour de 10-8 secondes au lieu de 10-10 précédement.

Champ de vitesse dans le domaine au moment où le calcul diverge

 

On peut ici observer les grands tourbillons de recirculation après impact à la paroi. A nouveau, le calcul diverge lorsque les gradients de vitesse deviennent importants au niveau de la sortie.

Cette simulation nous conforte dans l'idée que la divergence du calcul provient de variations brutales de la vitesse au niveau de la sortie. Ainsi, il s'agira dans le suite du projet d'atténuer ces brusques variations au niveau de la condition de sortie.

Troisième simulation - sortie modifiée

Présentation

Démarche

Afin d'éviter d'obtenir des gradients de vitesse en sortie trop importants, nous avons décidé de modéliser une partie de la conduite de sortie afin que le profil de vitesse ait le temps de s'homogénéiser en partie avant d'atteindre la condition de sortie.

 

Maillage

Toute la partie du domaine dont la géométrie reste inchangée demeure maillée de la même manière. La partie ajoutée en sortie apporte donc des mailles supplémentaires, raffinée en arctangente pour ne pas entraîner de "saut" dans les dimensions des mailles en entrée de conduite de sortie. On obtient ainsi le maillage suivant :

Maillage du domaine - sortie modifiée

Toute la partie du maillage qui a été ajoutée et qui n'est pas comprise dans la conduite de sortie n'est pas prises en compte par Jadim. Les équations ne sont pas résolues dans ce domaine. Ainsi le maillage "effectif" à considérer est le suivant :

Maillage du domaine de résolution des équations

 

Résultats et commentaires

La modification de la sortie a permis au calcul de converger. Nous avons cependant du l'arreter au bout d'un temps relativement court à cause des limites de temps du projet. Nous avons donc laisser le calcul tourner 764 000 itérations, ce qui correspond à 7.5ms de temps physique. Le résultat de la simulation et visible sur la vidéo suivante :

On observe bien les boucles de recirculation qui s'étendent jusqu'aux extrémités du domaine. D'autre part, nous pouvons nous intéresser aux longueurs de décroissance de la vitesse selon x et y afin d'estimer les différentes échelles de la turbulence comme expliqué au cours de l'étude théorique. 

Considérons tout d'abord la décroissance avec y de la vitesse. On se place en x égale au x du centre du jet. On trace ensuite la vitesse moyenne du fluide en fonction de y ; on obtient alors la courbe suivante : 

U moyen en fonction de Y en X = X(centre du jet)

On observe alors que la vitesse décroît de moitié sur une longueur Ydiss​ = 3.59 e-3 m. 

Pour observer la décroissance de la vitesse avec x, la question de la position en y se pose. En effet, la longueur selon x que met le fluide à perdre la moitié de sa vitesse dépend de la position en y choisie. Afin de nous placer au centre de notre domaine de dissipation, nous avons choisit de nous placer en Y=Ydiss/2. En traçant la vitesse moyenne du fluide en fonction de x, on obtient alors : 

U moyen en fonction de X en Y<Ydiss/2 (bleu), Y=Ydiss​/2 (rouge) et Y>Ydiss/2 (vert)

On trouve alors une valeur de Xdiss​=1e-5 m, c'est à dire égale à la largeur du jet. 

En reprenant les formules énoncées dans l'étude théorique, on trouve alors : 

  • Adiss=3.59 e-8 m2
  • $\epsilon$ = 5.58 e4 m2.s-3
  • $\eta$ = 2.06 e-6 m

Le maillage, en entrée du domaine, comprend des mailles de 2e-7 m selon x et de 3e-6 m selon y. Il devrait donc permettre de capter de manière satisfaisante les plus petites échelles de la turbulence. On se propose a présent d'observer le champ de vitesse du fluide en entrée du domaine : 

Champ de vitesse en entrée du domaine - visualisation des petits tourbillons (en noir)

On devine alors certaines structures tourbillonnaires, cependant la finesse ne permet pas une bonne visualisation des plus petites structures. 

D'autre part, en nous basant sur le Reynolds d'entrée et de la valeur de $\eta$ que nous venons de trouver, nous avons pu évaluer la dimension des plus grandes structures de la turbulence. On trouve alors : l0=1.1e-4 m. Visualisons à présent le champ de vitesse à l'échelle du domaine : 

Champ de vitesse dans le domaine - visualisation des grands tourbillons (en noir) - axes en mm

Les dimensions des plus grands tourbillons sont ici de quelques millimètres. Cet écart peut être du à la différence entre le Reynolds d'entrée que nous avons pris pour évaluer l0 et le Reynolds de l'écoulement, difficilement évaluable.