Modélisation B
==============

Comportement du fluide : THMC = LIQU_GAZ_ATMavec une saturation constante :math:`S=1`


Caractéristiques de la modélisation B
-------------------------------------

Modélisation plane : D_PLAN_HM


.. image:: images/100010E20000159A000013A3E04FB40CCC36EFED.svg
    :width: 288
    :height: 288


.. _RefImage_100010E20000159A000013A3E04FB40CCC36EFED.svg:

1 maille DPQ8 de la modélisation D_PLAN_HM : HM_ DPQ8


Grandeurs testées et résultats
------------------------------

Discrétisation en temps : Plusieurs pas de temps (16) pour étudier l’évolution de la pression pendant la phase transitoire jusqu’à se stabiliser. Le schéma en temps est implicite :math:`(\vartheta =1)` .

Liste des instants de calcul en secondes :

:math:`1,5,10,50,100,500,{10}^{3},{5.10}^{3},{10}^{4},5.{10}^{4},{10}^{5},5.{10}^{5},{10}^{6},5.{10}^{6},{10}^{7},{10}^{10}`

Les inconnues nodales de pression de fluide évaluées dans Code_Aster sont des variations par rapport aux pressions initiales de référence définies sous le mot clé THM_INIT, c’est pourquoi ce tableau présente des variations de pression dans notre comparaison entre le calcul Code_Aster et la solution de référence.


.. csv-table::

    "**Nœud / point**", "**Numéro d’ordre**", "**Pression**", "**Référence**
    :math:`(\mathrm{Pa})`", "**Tolérance**
    :math:`(\text{\%})`"
    ":math:`\mathrm{N1}/A`", "1 (t=1 s)", ":math:`\mathrm{PRE1}`", "-3,98.10-2", "1.0"
    "", "2 (t=5 s)", ":math:`\mathrm{PRE1}`", "-1,99.10-1", "1.0"
    "", "3 (t=10 s)", ":math:`\mathrm{PRE1}`", "-3,98.10-1", "1.0"
    "", "4 (t=50 s)", ":math:`\mathrm{PRE1}`", "-1,99", "1.0"
    "", "8 (t=5.103s)", ":math:`\mathrm{PRE1}`", "-1,95.10+2", "1.0"
    "", "16 (t=1010s)", ":math:`\mathrm{PRE1}`", "-5.10+3", "1.0"
    ":math:`\mathrm{N3}/C`", "1 (t=1 s)", ":math:`\mathrm{PRE1}`", "3,98.10-2", "1.0"
    "", "2 (t=5 s)", ":math:`\mathrm{PRE1}`", "1,99.10-1", "1.0"
    "", "3 (t=10 s)", ":math:`\mathrm{PRE1}`", "3,98.10-1", "2.0"
    "", "4 (t=50 s)", ":math:`\mathrm{PRE1}`", "1,99", "2.0"
    "", "8 (t=5.103s)", ":math:`\mathrm{PRE1}`", "1,95.10+2", "1.0"
    "", "16 (t=1010s)", ":math:`\mathrm{PRE1}`", "5.10+3", "1.0"



Remarques
---------

On remarque que les pressions calculées pour les deux comportements précédents (THMC=LIQU_SATU (modèle :math:`A` ) et THMC=LIQU_GAZ_ATM (modèle :math:`B` )) sont égales en valeurs absolues. La différence de signes est due au fait que :


*   la pression :math:`\mathrm{PRE1}` évaluée dans le code est la pression de l’eau pour le comportement THMC=LIQU_SATU,


*   :math:`\mathrm{PRE1}` est égale à la pression capillaire pour le comportement THMC=LIQU_GAZ. La pression capillaire est égale à la différence entre la pression de gaz et la pression liquide. Dans le cas particulier où la pression d’air sec est la pression atmosphérique (THMC=LIQU_GAZ_ATM), la pression capillaire a pour valeur l’opposé de la pression liquide.