Modélisation J
==============

Cette modélisation s'appuie sur l'ensemble des caractéristiques du problème de référence 2D décrit à la page  :ref:`2 <RefNumPara__3514_276907781>` ., sauf que l'on considère une modélisation axisymétrique. Il s'agit donc d'un barreau cylindrique, mais la solution analytique décrite à la page  :ref:`4 <RefNumPara__2359_165696795>` reste valable pour ce problème.


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

Pour cette modélisation, le barreau est cylindrique et ses dimensions sont les suivantes:


*   rayon :math:`\mathit{Lx}=0,5m`


*   hauteur :math:`\mathit{Ly}=\mathrm{5m}`

L'interface est introduite à mi-hauteur du barreau par une level set :math:`\mathit{LSN}` d'équation:

:math:`\mathit{LSN}` : :math:`y`

On utilise la modélisation AXIS du phénomène THERMIQUE.


Caractéristiques du maillage
----------------------------

Le maillage comporte 5 mailles de type TRIA3.


.. image:: images/10000201000000EB0000023FF4A3DC63BDBE1701.png
    :width: 1.2165in
    :height: 2.9807in


.. _RefImage_10000201000000EB0000023FF4A3DC63BDBE1701.png:

Figure 12.2-1: Maillage J


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

On teste d'abord les valeurs des degrés de liberté classique TEMP et Heaviside H1 du champ de température en sortie de l'opérateur THER_LINEAIRE, aux nœuds situés juste en dessous (2 nœuds) et au dessus de l'interface (2 nœuds).


.. csv-table::

    "**Identification**", "**Type de référence**", "**Valeur de référence**", "**Tolérance**"
    "Tous les nœuds situés juste au dessus de l'interface - :math:`\mathit{TEMP}`", "'ANALYTIQUE'", "20", "0,1%"
    "Tous les nœuds situés juste en dessous de l'interface - :math:`\mathit{TEMP}`", "'ANALYTIQUE'", "10", "0,1%"
    "Tous les nœuds situés juste en dessous / au dessus de l'interface - :math:`\mathit{H1}`", "'ANALYTIQUE'", "5", "0,1%"


On teste ensuite la valeur du degré de liberté TEMP du champ de température en sortie de POST_CHAM_XFEM, aux nœuds situés juste en dessous et au dessus de l'interface.


.. csv-table::

    "**Identification**", "**Type de référence**", "**Valeur de référence**", "**Tolérance**"
    "Tous les nœuds situés juste en dessous de l'interface - :math:`\mathit{TEMP}`", "'ANALYTIQUE'", "10", "0,1%"
    "Tous les nœuds situés juste au dessus de l'interface :math:`\mathit{TEMP}` -", "'ANALYTIQUE'", "20", "0,1%"


On teste enfin la valeur de la composante TEMP du champ TEMP_ELGA sur les points de Gauss situés en dessous et au dessus de l'interface ( *cf.*  remarque page 6).


.. csv-table::

    "**Identification**", "**Type de référence**", "**Valeur de référence**", "**Tolérance**"
    "Sur les points de Gauss situés en dessous de l'interface - :math:`\mathit{TEMP}`", "'ANALYTIQUE'", "10", "0,1%"
    "Sur les points de Gauss situés au dessus de l'interface - :math:`\mathit{TEMP}`", "'ANALYTIQUE'", "20", "0,1%"