v7.32.113 WTNP113 - Resaturation d’une alvéole#
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.2 du plan de qualification du projet ALLLIANCES [bib1].
Il représente l’évolution de la saturation en eau du milieu autour d’une alvéole de stockage. Deux phases sont prises en compte, une phase de désaturation lors de l’exploitation de l’ouvrage souterrain et une phase de resaturation après le remblaiement de galerie des alvéoles.
Solution de référence.#
On s’appuie sur les résultats du plan de qualification du projet ALLIANCES [bib1].
On effectue également des tests de non régression.
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#
Les résultats présentés ici résultent de la modélisation en déformations planes D_PLAN_THHD effectuée avec 2988éléments TRI3.
Résultats#
On présente les profils de pression capillaire et de saturation sur une coupe horizontale (\(y=\mathrm{450m}\) ) et verticales (\(x=\mathrm{7m}\) ).
Figure 3.3-a : Profils de pression capillaire \(Y=-\mathrm{450m}\)
Figure 3.3-b : Profils de saturation \(Y=-\mathrm{450m}\)
Figure 3.3-c : Profils de pression capillaire \(X=7m\)
Figure 3.3-d : Profils de saturation \(X=7m\)
Grandeurs testées et résultats#
\(X\) ( \(m\) ) |
\(Y\) ( \(m\) ) |
Temps (ans) |
\(\mathrm{PRE1}\) ( \(\mathrm{Pa}\) ) Source Externe |
Tolérance ( \(\text{\%}\) ) |
5 |
-450 |
27 |
2.9E7 |
20 |
5 |
-450 |
10000 |
4.35E6 |
20 |
7 |
-450 |
27 |
3.6E6 |
20 |
7 |
-450 |
10000 |
4.34E6 |
20 |
Les résultats sont qualitativement conformes à ceux que l’on trouve dans la littérature d’Alliances:on ne possède que des isovaleurs en référence, aussi la tolérance d’erreur est assez grande (\(\text{20\%}\) ).
Modélisation B#
Caractéristiques de la modélisation B#
Il s’agit de la même modélisation que ci-dessus mais avec la modélisation miscible THH2D (avec un coefficient de Henry infini). Les résultats attendus doivent donc être exactement les mêmes.
Grandeurs testées et résultats#
\(X\) ( \(m\) ) |
\(Y\) ( \(m\) ) |
Temps (ans) |
\(\mathit{PRE1}\) ( \(\mathit{Pa}\) ) Source Externe |
Tolérance ( \(\text{\%}\) ) |
5 |
-450 |
27 |
2.9E7 |
20 |
5 |
-450 |
10000 |
4.35E6 |
20 |
7 |
-450 |
27 |
3.6E6 |
20 |
7 |
-450 |
10000 |
4.34E6 |
20 |
Les résultats sont qualitativement conformes à ceux que l’on trouve dans la littérature d’Alliances:on ne possède que des isovaleurs en référence, aussi la tolérance d’erreur est assez grande (\(\text{20\%}\) ).
Modélisation C#
Caractéristiques de la modélisation C#
Même modélisation que la A mais en sélective.
Grandeurs testées et résultats#
\(X\) ( \(m\) ) |
\(Y\) ( \(m\) ) |
Temps (ans) |
\(\mathrm{PRE1}\) ( \(\mathrm{Pa}\) ) Source Externe |
Tolérance ( \(\text{\%}\) ) |
5 |
-450 |
27 |
2.9E7 |
20 |
5 |
-450 |
10000 |
4.35E6 |
20 |
7 |
-450 |
27 |
3.6E6 |
20 |
7 |
-450 |
10000 |
4.34E6 |
20 |
Les résultats sont qualitativement conformes à ceux que l’on trouve dans la littérature d’Alliances:on ne possède que des isovaleurs en référence, aussi la tolérance d’erreur est assez grande (\(\text{20\%}\) ).
Modélisation D#
Caractéristiques de la modélisation D#
Même modélisation que la B mais en sélective.
Grandeurs testées et résultats#
\(X\) ( \(m\) ) |
\(Y\) ( \(m\) ) |
Temps (ans) |
\(\mathit{PRE1}\) ( \(\mathit{Pa}\) ) Source Externe |
Tolérance ( \(\text{\%}\) ) |
5 |
-450 |
27 |
2.9E7 |
20 |
5 |
-450 |
10000 |
4.35E6 |
20 |
7 |
-450 |
27 |
3.6E6 |
20 |
7 |
-450 |
10000 |
4.34E6 |
20 |
Les résultats sont qualitativement conformes à ceux que l’on trouve dans la littérature d’Alliances:on ne possède que des isovaleurs en référence, aussi la tolérance d’erreur est assez grande (\(\text{20\%}\) ).
Synthèse des Résultats#
D’une façon générale, les résultats obtenus avec Code_Aster sont qualitativement conformes à ceux que l’on trouve dans la littérature d’ALLIANCES.