Modélisation D
==============

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


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


*   Modélisation volumique: 3D_HM


*   1 maille HEXA20 de la modélisation 3D_HM : HM_HEXA20


.. image:: images/10010F2C00001C0000001C3388288620E469299A.svg
    :width: 277
    :height: 290


.. _RefImage_10010F2C00001C0000001C3388288620E469299A.svg:


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**", "**Numéro d’ordre**", "**Pression**", "**Référence**
    :math:`(\mathrm{Pa})`", "**Tolérance**
    :math:`(\text{\%})`"
    ":math:`\mathrm{NO20}`", "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{NO1}`", "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
---------

De même que pour la modélisation bidimensionnelle, on remarque que les pressions calculées pour les deux comportements précédents (THMC=LIQU_SATU (modèle :math:`C` ) et THMC=LIQU_GAZ_ATM (modèle :math:`D` )) sont égales en valeurs absolues. La différence de signes est due au fait que :


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


*   et pour le comportement THMC=LIQU_GAZ, :math:`\mathrm{PRE1}` est égale à la pression capillaire. 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.