4. Modélisation C#

4.1. Caractéristiques de la modélisation et du maillage#

Il s’agit d’une modélisation D_PLAN_HHMS. Le maillage est composé de 200 mailles QUA8 et 206 mailles SEG3.

4.2. Grandeurs testées et résultats#

On réalise un test de non-régression sur les flux hydrauliques sur l’arête \(x=5\) m un jour après le début de la simulation. Ces flux sont récupérés par les composantes PRE1 et PRE2 du mot-clé RESULTANTE de l’opérateur POST_RELEVE_T.

Composante

Type de référence

Valeur

PRE1

NON_REGRESSION

\(-11.5681454\) kg/s

PRE2

NON_REGRESSION

\(0.0072103482\) kg/s

On affiche ci-dessous l’opposé de l’intégrale de ces flux. Il s’agit donc des apports de masse de liquide et de gaz dans l’ensemble du barreau. Ceux-ci convergent vers deux valeurs asymptotiques strictement égales à :

(4.1)#\[\begin{split}\left\lbrace \begin{array}{l} \Delta M_{lq,\infty} = \phi^0\rho_{lq}^0\left(1-S_{lq}^0\right)|\Omega_0|\\ \Delta M_{gz,\infty} = -\cfrac{p_{gz}^0M^{mol}_{gz}}{RT}\phi^0S_{lq}^0|\Omega_0| \end{array} \right.\end{split}\]

avec \(p_{gz}^0 = 100\) kPa, \(T = 273\) K (c’est-à-dire les valeurs de référence en pression de gaz et en température), et \(|\Omega_0|=L\times W\times 1\text{ m} = 5\) m3. On notera que l’apport en gaz est négatif, celui-ci s’échappant de l’échantillon lors de sa saturation.

_images/resu_modelisation_c.svg