Modélisation A ============== Caractéristiques de la modélisation A ------------------------------------- Modélisation HHS en déformations planes. Couplage LIQU_GAZ. 100 éléments QUAD8. Discrétisation en temps: * :math:`\mathrm{1000s}` en 20 pas de temps * :math:`5000s` en 20 pas de temps * :math:`10000s` en 10 pas de temps * :math:`50000s` en 20 pas de temps * :math:`100000s` en 20 pas de temps Résultats --------- Les figures ci-dessous présentent les profils de pressions de gaz, pressions capillaires et saturations le long du bi-matériaux à différents temps: .. image:: images/100002000000022600000144B9004D8AD4F20AA2.png :width: 5.7189in :height: 3.3665in .. _RefImage_100002000000022600000144B9004D8AD4F20AA2.png: .. image:: images/10000200000001EA0000012FBA65E18634443F87.png :width: 5.0984in :height: 3.15in .. _RefImage_10000200000001EA0000012FBA65E18634443F87.png: On observe bien la désaturation de la barrière géologique par la barrière ouvragée. Les profils de pression de gaz sont eux aussi caractéristiques de ce type de problème : on observe un pic de pression de gaz à l'interface des deux matériaux qui est du au fait que le gaz est comprimé par l'eau qui resature la barrière ouvragée. Si l'existence de ce pic a une réalité physique, son ampleur est en revanche due à un problème de schéma numérique (problème connu) comme nous le verrons pour la modélisation :math:`C` . Les oscillations observées sur les saturations sont également bien connues (on se réfèrera à [ :ref:`1 ` ]). Les matériaux sont en effet définis aux éléments et non aux nœuds. Ce qui signifie que sur un nœud situé à l’interface, 2 courbes de pressions capillaires sont définies et différentes. Ce point explique que si la continuité de pression capillaire est assurée il n’en est pas de même pour la saturation. Valeurs testées --------------- Ce cas test n’a pas de valeur de référence, on en fait donc un cas de non régression. On effectue des tests sur deux valeurs : .. csv-table:: ":math:`X` **(** :math:`m` **)**", "**Temps (** :math:`s` **)**", ":math:`\mathrm{PRE1}` **Aster**", "**Erreur relative autorisée**" ":math:`0.03`", ":math:`5000`", ":math:`7.5{10}^{5}`", ":math:`1\text{\%}`" ":math:`0.03`", ":math:`100000`", ":math:`4.48{10}^{6}`", ":math:`1\text{\%}`"