Modélisation A
==============

Caractéristiques de la modélisation et du maillage
--------------------------------------------------

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

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

On réalise un test de non-régression sur les flux hydrauliques et thermiques sur l'arête :math:`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 TEMP du mot-clé RESULTANTE de l'opérateur POST_RELEVE_T. On notera que ces composantes sont négatives, car les flux sont pris positifs dans le sens de la normale sortante au barreau.



.. csv-table::

    "**Composante**", "**Type de référence**", "**Valeur**"

    "PRE1", "NON_REGRESSION", ":math:`-0.09639395` kg/s"
    "TEMP", "NON_REGRESSION", ":math:`-6893720649` J/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 d'entropie dans l'ensemble du barreau. Ceux-ci convergent vers deux valeurs asymptotiques qui s'approchent par :


.. math::
    \left\lbrace
    \begin{array}{l}
    \Delta M_{lq,\infty} = \cfrac{\phi^0\rho_{lq}^0}{K_{lq}}\Delta p_{lq}|\Omega_0| \\
    \Delta S_\infty = \left[r_0C_\sigma^S+\phi^0\rho_{lq}^0(C_{lq}^p-C_\sigma^S)\right]\Delta T|\Omega_0|
    \end{array}
    \right.
    :label: eq-1

avec :math:`\Delta p_{lq} = 100` kPa, :math:`\Delta T = 100` K, et :math:`|\Omega_0|=L\times W\times 1\text{ m} = 5` m3.

En effet, les relations :eq:`eq-1` sont exactes dans l'hypothèse des petites perturbations pour le fluide et le squelette. Dans le cas présent, la figure suivante montre donc bien que les conditions du cas-test ne s'en éloignent que très peu.


.. figure:: images/resu_modelisation_a.svg
 :align: center
 :width: 960
 :height: 380