3. Modélisation A#

3.1. Caractéristiques de la modélisation#

Cette modélisation permet de valider l’évolution de la température lors d’un échange thermique non-linéaire (rayonnement).

On n’impose pas de conditions de Dirichlet mais une condition d’échange sur les nœuds aux bords.

CHARGE = AFFE_CHAR_THER_F(
    MODELE=MODELE,
    RAYONNEMENT=_F(
        GROUP_MA="LINE_EXT", SIGMA=C_sigma, EPSILON=C_epsilon, TEMP_EXT=TEMP_EXT
    ),
)

La température sur les 4 nœuds de la maille devrait être la même à un instant donné.

La température en régime stationnaire devrait être la même que la température externe.

Pour imposer l’utilisation de Fortran avec THER_NON_LINE, on rajout le mot-clé :

OBSERVATION=_F(NOM_CHAM="TEMP", NOM_CMP="TEMP", TOUT="OUI"),

3.2. Caractéristiques du maillage#

Nombre de nœuds : 4

Nombre de mailles et types : 1 QUAD4

3.3. Valeurs testées#

Le nœud observé a pour coordonnées \((0.0, 0.0)\), mais d’autres points sont possible (puisque la température devrait être la même).

Avec fortran
temps	température	référence	error
100	54.533188379364574	54.56435728271139	5.712319E-02%
600	99.70542232072329	99.70768743887123	2.271759E-03%
STAT	100.00000002900663	100.0	2.900663E-08%

Avec Python
temps	température	référence	error
100	54.53443932011566	54.56435728271139	5.483060E-02%
600	99.70542525641507	99.70768743887123	2.268814E-03%
STAT	100.00000002900661	100.0	2.900661E-08%

où STAT fait référence au régime stationnaire.

Fig 3.1. Comparaison des résultats (avec Python) avec référence théorique