4. Modélisation B#
Comportement du fluide : THMC = LIQU_GAZ_ATMavec une saturation constante \(S=1\)
4.1. Caractéristiques de la modélisation B#
Modélisation plane : D_PLAN_HM
1 maille DPQ8 de la modélisation D_PLAN_HM : HM_ DPQ8
4.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 / point |
Numéro d’ordre |
Pression |
Référence \((\mathrm{Pa})\) |
Tolérance \((\text{\%})\) |
\(\mathrm{N1}/A\) |
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{N3}/C\) |
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 |
4.3. Remarques#
On remarque que les pressions calculées pour les deux comportements précédents (THMC=LIQU_SATU (modèle \(A\) ) et THMC=LIQU_GAZ_ATM (modèle \(B\) )) sont égales en valeurs absolues. La différence de signes est due au fait que :
la pression \(\mathrm{PRE1}\) évaluée dans le code est la pression de l’eau pour le comportement THMC=LIQU_SATU,
\(\mathrm{PRE1}\) est égale à la pression capillaire pour le comportement THMC=LIQU_GAZ. 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.