6. Modélisation D

6.1. Caractéristiques de la modélisation D

Les résultats présentés sont issus de la modélisation en 3D_THH2MS. Il s’agit de la même modélisation que la modélisation A mais avec un terme d’échange non linéaire en HR et non en pression.

Pour compléter, les calculs sont également réalisés en 3D_HH2MS, en 3D_THH2S et en en 3D_HH2S.

Les déplacements et températures étant bloqués, les résultats doivent être les mêmes dans tous les cas. Le but est ici de vérifier la bonne prise en compte de ce terme d’échange.

Le maillage est composé de 640 éléments HEXA20 .

6.2. Valeurs testées et résultats

On vérifie que le flux d’eau (“FH11X”+”FH12X” récupéré par SIEF_NOEU ) obtenu sur le bord d’imposition correspond bien à son expression analytique, à savoir \(q_{w}^{ext}{.}n= \alpha.\rho_{vp}.\left( exp\left[ \frac{Pc.M_{H2O}}{\rho_{l}.R.T}\right] -{HR}^{ext} \right)\). Une tolérance de 10% est permise du fait de la taille du maillage et de l’erreur générée par l’interpolation au noeud.

Points \((x,y)\)	Temps ( \(s\) )	“FH11X”+”FH12X”	Valeur analytique
\((0,2;0)\)	1800	7.3754E-09	8.09e-09

Pour des raisons de gain de temps les modélisations autres que 3D_THH2MS ne seront testées que pour 60 s. On vérifie que les pressions à 60 s sont toutes les mêmes quelles que soient les modélisations.

Points \((x,y)\)	Temps ( \(s\) )	PRE1
\((0,2;0)\)	1800	10017