Modeling A
==============


Characteristics of modeling
-----------------------------------

Here we validate the behavior of the MACR_ECREVISSE command in the presence of a horizontal through crack. A modeling in plane deformations is chosen for the mechanics D_ PLAN, and a thermal modeling PLAN_DIAG, with a diagonalized mass matrix for the thermal one.

The flow characteristics have been described with the geometry in §1.1.


Characteristics of the mesh
----------------------------

The mesh used is shown:


.. image:: images/10000000000002040000014A44769FA4A4278A52.png
    :width: 4.822in
    :height: 2.9917in


.. _RefImage_10000000000002040000014A44769FA4A4278A52.png:

**Figure** 3.2-a **: Full mesh**

The dimensions are the same as for the geometry.

Number of knots: 1260

Number of meshes: 1156 QUAD4

The growth parameters used to perform the mesh are as follows: from :math:`0.002m` (intrados) to :math:`0.05m` (extrados), and idem from BFISH and BFISB to ENCASTRE.


Sizes and results
----------------------

The significant quantities of this calculation are tested, namely: the flow rate, as well as the displacement and the temperature at the inlet and outlet of cracks. Two moments are tested.

At :math:`t=500s`:


.. csv-table::

    "**Identification**", "**Tolerance**"
    "Move :math:`\mathit{DY}` to node :math:`{E}_{\mathit{HD}}` :math:`\mathrm{[}m\mathrm{]}` ", "1.0E-06"
    "Material temperature at node :math:`{E}_{\mathit{HD}}` :math:`\mathrm{[}°C\mathrm{]}` ", "1.0E-06"
    "Move :math:`\mathit{DY}` to node :math:`{I}_{\mathit{HD}}` :math:`\mathrm{[}m\mathrm{]}` ", "1.0E-06"
    "Temperature at node :math:`{I}_{\mathit{HD}}` :math:`\mathrm{[}°C\mathrm{]}` ", "1.0E-06"
    "Flow", "1.0E-06"


At :math:`t=10000s`:


.. csv-table::

    "**Identification**", "**Tolerance**"
    "Move :math:`\mathit{DY}` to node :math:`{E}_{\mathit{HD}}` :math:`\mathrm{[}m\mathrm{]}` ", "1.0E-06"
    "Material temperature at node :math:`{E}_{\mathit{HD}}` :math:`\mathrm{[}°C\mathrm{]}` ", "1.0E-06"
    "Move :math:`\mathit{DY}` to node :math:`{I}_{\mathit{HD}}` :math:`\mathrm{[}m\mathrm{]}` ", "1.0E-06"
    "Temperature at node :math:`{I}_{\mathit{HD}}` :math:`\mathrm{[}°C\mathrm{]}` ", "1.0E-06"
    "Flow", "1.0E-06"



notes
---------

At :math:`t=10000s`, the I_HD node moved about :math:`\mathrm{-}70\mu m`, (in particular due to thermal expansion), which corresponds to the half-opening of the crack. This means that the lips are in contact (keyword DEFI_CONTACT activated) and the crack is closed. On the other hand, for the hydraulic part (managed by *Crayvisse*), a persistent opening is fixed, so that the flow is never zero.