Mise en place des calculs

Calcul en un seul bloc

Le calcul se fera comme précédemment avec le solveur buoyantBoussinesqSimpleFoam, avec cette fois-ci une modification de la géométrie ( l'entrée ne commence plus à x=1.44m mais plutôt à 0). Le solveur potentialFoam n'est plus utilisé puisqu'il n'y a plus de problème de discontinuité des conditions de vitesses, et les conditions aux limites sont résumés dans le tableau ci-dessous :

  Entrée Plaque Sortie Bord
U (m/s) nonuniform List<vector>

fixedValue uniform (0 0 0)

zeroGradient zeroGradient
p (Pa) zeroGradient zeroGradient

fixedValue uniform

fixedValue uniform
T (K) nonuniform List<scalar>

fixedValue uniform 333.15K

zeroGradient zeroGradient

 

Calculs en deux blocs

Dans le premier domaine (de x=0 à x=0.8 m), nous lançons un premier calcul sous buoyantBoussinesqSimpleFoam avec cette fois-ci une viscosité turbulente nulle pour forcer le régime laminaire, et les conditions aux limites seront exactement identiques à celles utilisées précedemment.

Une fois le calcul convergé, nous récupérons les résultats au dernier pas de temps en x= 0.8m, pour ensuite les imposer comme conditions d'entrée pour le second domaine/calcul.

La récupération de ces profils en entrée peut se faire en passant d'abord par le logiciel "pararaview", pour ensuite enregistrer les datas sous format texte (data.dat), et enfin faire une des opérations de tri via quelques commandes linux :

  • Afin d'extraire les lignes contenant uniquement les données relatifs à la position x=0.8 et z=0.001

> grep -w 0.8 data.dat | grep -w 0.001  > Y.dat

  • Pour extraire les composantes des vitesses (colonnes 3, 4 et 5)

> cut -d',' -s -f3,4,5 Y.dat > U.dat

  • Pour extraire le profil de température (colonne 2)

> cut -d',' -s -f2 Y.dat > T.dat

  • Pour extraire le profil des pressions (colonne 1)

> cut -d',' -s -f1 Y.dat > p.dat

  • Pour respecter la forme avec laquelle s'écrit un profil de vitesses de type nonuniform List<vector>

> sed 's/^/(/' U.dat > U.dat (ouvrir une parenthèse au début de chaque ligne)

> sed 's/$/)/' U.dat > U.dat (fermer les parenthèses à la fin de chaque ligne)

> sed -e 's/,/ /g' U.dat > U.dat (remplacer le virgules par des espaces)

Après avoir récupérer les profils de vitesse, pression et température en sortie du premier bloc, on passe au calcul dans le deuxième bloc où les conditions en entrée seront de type nonuniform. Ainsi, il va falloir copier les dernier profils de vitesse, température et pression, lancer un premier calcul sous potentialFoam pour résoudre les problèmes de discontinuité, et ensuite réaliser un dernier calcul via buoyantBoussinesSimpleFoam.