9. Modélisation G#
Comportement du fluide : THMC = LIQU_SATU
9.1. Caractéristiques de la modélisation G#
Modélisation volumique
1 maille HEXA20 de la modélisation 3D_THM : THM_HEXA20
9.2. Grandeurs testées et résultats#
Discrétisation en temps : Plusieurs pas de temps (16) pour étudier l’évolution de la pression pendant la phase transitoire jusqu’à se stabiliser. Le schéma en temps est implicite \((\vartheta =1)\) .
Liste des instants de calcul en secondes :
\(1,5,10,50,100,500,{10}^{3},{5.10}^{3},{10}^{4},5.{10}^{4},{10}^{5},5.{10}^{5},{10}^{6},5.{10}^{6},{10}^{7},{10}^{10}\)
Les inconnues nodales de pression de fluides évaluées dans Code_Aster sont des variations par rapport aux pressions initiales de référence définies sous le mot clé THM_INIT, c’est pourquoi ce tableau présente des variations de pression dans notre comparaison entre le calcul Code_Aster et la solution de référence.
Nœud |
Numéro d’ordre |
Pression |
Référence \((\mathrm{Pa})\) |
Tolérance \((\text{\%})\) |
\(\mathrm{NO20}\) |
1 (t=1 s) |
\(\mathrm{PRE1}\) |
3,98.10-2 |
1.0 |
2 (t=5 s) |
\(\mathrm{PRE1}\) |
1,99.10-1 |
1.0 |
|
3 (t=10 s) |
\(\mathrm{PRE1}\) |
3,98.10-1 |
1.0 |
|
4 (t=50 s) |
\(\mathrm{PRE1}\) |
1,99 |
1.0 |
|
8 (t=5.103s) |
\(\mathrm{PRE1}\) |
1,95.10+2 |
1.0 |
|
16 (t=1010s) |
\(\mathrm{PRE1}\) |
5.10+3 |
1.0 |
|
\(\mathrm{NO1}\) |
1 (t=1 s) |
\(\mathrm{PRE1}\) |
-3,98.10-2 |
1.0 |
2 (t=5 s) |
\(\mathrm{PRE1}\) |
-1,99.10-1 |
1.0 |
|
3 (t=10 s) |
\(\mathrm{PRE1}\) |
-3,98.10-1 |
2.0 |
|
4 (t=50 s) |
\(\mathrm{PRE1}\) |
-1,99 |
2.0 |
|
8 (t=5.103s) |
\(\mathrm{PRE1}\) |
-1,95.10+2 |
1.0 |
|
16 (t=1010s) |
\(\mathrm{PRE1}\) |
-5.10+3 |
1.0 |