v7.32.118 WTNP118 - Rééquilibrage gravitaire de la saturation d’une colonne#
Résumé:
Le test présenté ici permet de vérifier la bonne prise en compte de la gravité pour la modélisation des écoulements non saturés. Ce test représente une colonne initialement uniformément désaturé. Il n’y a aucun chargement: seule la gravité est moteur de l’évolution.
Ce test fera l’objet d’un cas test Alliances.
Modélisation A#
Caractéristiques de la modélisation A#
La modélisation étudiée est HHDen déformations planes. Le maillage est composé de 80 éléments Q8.
Fonctionnalités testées#
Cette modélisation permet en particulier de vérifier l’utilisation de la fonction temporelle PESA_MULT introduite dans la carte THM_DIFFU. Cette fonction est ajoutée en facteur des termes PESA_X, PESA_Y et PESA_Z, et permet de rendre la gravité dépendante du temps.
Résultats#
Les figures ci-dessous présentent les profils de pressions capillaires, pressions de gaz et saturation le long du barreau pour différents instants. On observe bien le rééquilibrage pour atteindre un état stationnaire (plus rien ne bouge après \(0,5s\) ).
A cet état stationnaireet en \(\mathrm{Pa}\) :
\(\mathrm{Pgz}(x=0)=0,991\mathrm{Pa}\) et \(\mathrm{Pgz}(x=3)=1,014\mathrm{Pa}\)
\(\mathrm{Pc}(x=0)=22,8\mathrm{Pa}\) et \(\mathrm{Pc}(x=3)=-6,62\mathrm{Pa}\)
Ce qui donne \({\mathrm{\Delta P}}_{\text{lq}}=29,44\mathrm{Pa}\) qui correspond bien à l’équilibre hydrostatique.
Figure 2.3-a : profils de pression capillaire
Figure 2.3-b : profils de pression de gaz
Figure 2.3-c : profils de saturation
Valeur testée#
\(X(m)\) |
Temps \((s)\) |
PRE1Aster \((\mathrm{Pa})\) |
Erreur relative autorisée |
0,02 |
13,94 |
0.1% |
|
0,1 |
21,92 |
0.1% |
|
1 |
22,79 |
0.1% |
Modélisation B#
Caractéristiques de la modélisation B#
La modélisation étudiée est HH2Den déformations planes. Le maillage est composé de 80 éléments Q8. Cette déclinaison n’a pour but que de passer par la loi de comportement LIQU_AD_GAZ_VAPEavec gravité. C’est exactement le même test que précédemment mais avec une coefficient de Henry \({K}_{H}=\infty\) et tous les coefficients de Fick pris égaux à 0.
Résultats#
Les résultats sont bien sûr les mêmes que précédemment:
\(X(m)\) |
Temps \((s)\) |
PRE1Aster |
Erreur relative autorisée |
0,02 |
13,94 |
0.1% |
|
0,1 |
21,92 |
0.1% |
|
1 |
22,79 |
0.1% |
Modélisation C#
Caractéristiques de la modélisation C#
La modélisation étudiée est THH2Den déformations planes. Le maillage est composé de 80 éléments Q8. C’est exactement le même cas que la modélisation B mais avec une structure THH2Det donc de la thermique (bloquée). Cette modélisation a uniquement pour but de se ramener à une structure de donnée THH2Dafin de pouvoir en faire un cas test Alliances (seule la modélisation THH2Dest connue) .
Résultats#
Les résultats sont bien sûr les mêmes que précédemment:
\(X(m)\) |
Temps \((s)\) |
PRE1Aster \((\mathrm{Pa})\) |
Erreur relative autorisée |
0,02 |
13,94 |
0.1% |
|
0,1 |
21,92 |
0.1% |
|
1 |
22,79 |
0.1% |
Synthèse#
Ce cas test permet de valider la prise en compte de la gravité pour les cas saturés. Les résultats obtenus sont cohérents.