5. Modélisation C

5.1. Caractéristiques de la modélisation

Modélisation PLAN_DIAG :

5.2. Caractéristiques du maillage

20 QUAD8

5.3. Remarque

La chaleur latente de fusion est fournie par l’intermédiaire de l’enthalpie sur un intervalle de \(0.01°C\) .

La matrice de masse est diagonalisée pour limiter les oscillations

5.4. Valeurs testées

Les nœuds observés ont pour coordonnée \(y=0.0\)

Identification température	Référence
t = 0.5 s N6 (x = 0.005)	682.43
t = 1.0 s N6 (x = 0.005)	661.33
t = 3.0 s N6 (x = 0.005)	628.20
t = 6.0 s N6 (x = 0.005)	614.25
t = 0.5 s N11 (x = 0.010)	726.05
t = 1.0 s N11 (x = 0.010)	705.75
t = 3.0 s N11 (x = 0.010)	669.63
t = 6.0 s N11 (x = 0.010)	647.49
t = 0.5 s N16 (x = 0.015)	738.11
t = 1.0 s N16 (x = 0.015)	728.70
t = 3.0 s N16 (x = 0.015)	696.06
t = 6.0 s N16 (x = 0.015)	673.22
t = 0.5 s N21 (x = 0.020)	739.86
t = 1.0 s N21 (x = 0.020)	737.22
t = 3.0 s N21 (x = 0.020)	714.94
t = 6.0 s N21 (x = 0.020)	692.06

Le calcul par éléments finis nécessite une discrétisation en temps de \(\Delta t=5.{10}^{-4}s\) au moins pour les premiers pas. La condition au limite imposée à l’origine faisant passer brusquement la température de \(740.°C\) à \(580.°C\) . On observe au niveau des premiers pas de temps quelques oscillations qui se stabilisent ensuite assez rapidement, malgré tout la température maximum est dépassée, il n’y a pas respect du maximum discret. Ce phénomène est observé lors des chocs thermiques, seul un traitement numérique particulier au niveau de la matrice de masse peut remédier à ce dernier.