v7.32.110 WTNP110 - Ecoulement orthotrope saturé 2D#
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 saturé orthotrope. Ce test correspond au test 1.2 du plan de qualification du projet ALLIANCES [bib1].
On valide également les estimateurs d’erreur a posteriori en résidu développés pour la THM.
La solution de référence est une solution analytique.
Solution de référence#
Méthode de calcul#
La solution analytique en pression est un polynôme de degré en \(x\) et \(y\) , la vitesse est constante et horizontale:
\(\begin{array}{}P(x,y)=-45x-80y+22.3\\ V(x,y)=(45{K}_{x};80{K}_{y})=(45;60)\end{array}\)
Grandeurs et résultats de référence#
On donne la valeur de la pression et de la vitesse en 3 points:
\(x\) |
-0,05 |
0,00 |
+0,05 |
\(y\) |
-0,05 |
0.00 |
0.05 |
\(P\) |
28,8 |
22.5 |
16.3 |
\(\mathit{Vx}\) |
45 |
45 |
45 |
\(\mathit{Vy}\) |
60 |
60 |
60 |
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 en déformations planes. \(20\times 20\) éléments \(\mathit{Q8}\)
Grandeurs testées et résultats#
Nœud |
Type de valeur |
Instant ( \(s\) ) |
Référence (analytique) |
Aster |
Différence ( \(\text{\%}\) ) |
\(\mathit{N441}\) \((0.05;0.05)\) |
\(\mathrm{PRE1}\) |
1 |
16,3 |
16,25 |
0,3 % |
\(\mathrm{N241}\) \((-0.05;-0.05)\) |
\(\mathrm{PRE1}\) |
1 |
28,8 |
28,5 |
0.96 % |
\(\mathit{N341}\) \((\mathrm{0 };0)\) |
\(\mathit{PRE1}\) |
1 |
22,5 |
22,5 |
|
On valide également les estimateurs d’erreur en temps et en espace par des tests de non-régression.
Type de valeur |
Instant ( \(s\) ) |
Aster |
Différence ( \(\text{\%}\) ) |
ERRE_TPS_GLOB |
1 |
1,5E-8 |
0 % |
Nœud |
Type de valeur |
Composante |
Instant ( \(s\) ) |
Aster |
Différence ( \(\text{\%}\) ) |
\(\mathrm{N313}\) \((-0.09;-0.01)\) |
ERME_NOEU |
ERHMHY_G |
1 |
1,9E-12 |
|
Synthèse des résultats#
Très bon accord avec la solution de référence.