v7.32.112 WTNP112 - Resaturation d’une colonne#
Résumé:
Le test présenté ici permet de vérifier le bon fonctionnement des opérateurs utilisés pour la résolution des équations d’un écoulement en milieu insaturé. Ce test correspond au test 3.1 du plan de qualification du projet ALLLIANCES [bib1]. Il représente la resaturation d’une colonne.
Solution de référence#
La solution de référence est la solution semi-analytique de Phillips [bib1].
Références bibliographiques#
Projet Alliances plan de qualification, note ANDRA CNT-ASCS 02-075B
Modélisation A#
Caractéristiques de la modélisation A#
Modélisation THHD en déformations planes. 40 élément \(\mathit{Q8}\)
Résultats#
Le tableau ci-dessous présente les profils de pressions capillaires et saturation le long du barreau et à chaque heure: On intercale entre les 2 profils, celui correspondant à la solution semi-analytique de Philips.
Valeur testée#
\(X\) ( \(m\) ) |
Temps ( \(s\) ) |
\(\mathit{PRE1}\) Aster |
Pression capillaire (Phillips) |
Erreur relative |
36000s |
0.461 |
0.460 |
\(\text{0.2\%}\) |
Modélisation B#
Caractéristiques de la modélisation B#
Modélisation THHD en déformations planes. 80 éléments \(\mathit{Q8}\)
Résultats#
Le tableau ci-dessous présente les profils de pressions capillaires et saturation le long du barreau et à chaque heure, comme pour la modélisation A on intercale la solution semi-analytique de Phillips:
Valeur testée#
\(X\) ( \(m\) ) |
Temps ( \(s\) ) |
\(\mathit{PRE1}\) Aster |
Pression capillaire (Phillips) |
Erreur relative |
36000s |
0.459 |
0.460 |
0.1% |
Modélisation C#
Caractéristiques de la modélisation C#
Modélisation THH2D en déformations planes. 80 éléments Q8. C’est exactement le même cas que la modélisation B mais avec une structure THH2D et un coefficient de Henry infini. Cette modélisation a uniquement pour but de se ramener à une structure de donnée THH2D qui est celle connue par Alliances.
Résultats#
Les résultats sont bien évidemment les même que pour la modélisation B.
Modélisation D#
Caractéristiques de la modélisation D#
Modélisation HMD en déformations planes. 80 éléments Q8. C’est exactement le même cas que la modélisation B mais avec la loi de couplage LIQU_GAZ_ATM spécifique à la modélisation de Richards (équivalente à modélisation non saturé avec pression de gaz imposée). On est bloqué en mécanique.
Résultats#
Les résultats sont bien évidemment les même que pour la modélisation B.