10. Modélisation H#
Comportement du fluide : THMC = LIQU_GAZ_ATMavec une saturation constante \(S=1\)
10.1. Caractéristiques de la modélisation H#
Modélisation volumique : 3D_THM
1 maille HEXA20 de la modélisation 3D_THM : THM_HEXA20
10.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 fluide é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 |
10.3. Remarques#
De même que pour la modélisation bidimensionnelle, on remarque que les pressions calculées pour les deux comportements précédents (THMC=LIQU_SATU (modèle \(G\) ) et THMC=LIQU_GAZ_ATM (modèle \(H\) )) sont égales en valeurs absolues. La différence de signes est due au fait que :
pour le comportement THMC=LIQU_SATU, la pression \(\mathrm{PRE1}\) évaluée dans le code est la pression de l’eau ,
et pour le comportement THMC=LIQU_GAZ, \(\mathrm{PRE1}\) est égale à la pression capillaire. La pression capillaire est égale à la différence entre la pression de gaz et la pression liquide. Dans le cas particulier où la pression d’air sec est la pression atmosphérique (THMC=LIQU_GAZ_ATM), la pression capillaire a pour valeur l’opposé de la pression liquide.