3. Modélisation A#
3.1. Caractéristiques de la modélisation#
La modélisation est réalisée en déformation plane avec 100 éléments QU4 pour le jeu, 3202 éléments TRI3 pour le massif et 100 éléments QU4 pour la discontinuité.
Discrétisation en temps : 24 pas de temps pour une simulation de 1000 ans.
Remarque: Le test est de type D_PLAN_THMS mais la partie thermique est paramétrée de manière à ce que le résultat soit identique au cas sans thermique, afin de vérifier informatiquement que les éléments D_PLAN_THMS sont bien corrects.
3.2. Grandeurs testées et résultats#
La figure montre les profils de pression sur la coupe verticale \(x=4m\) à différents instants. On y voit bien l’influence hydraulique de la fissure. En effet, au voisinage de celle-ci on observe des gradients de pression très importants qui correspondent à des flux dirigés vers chacun des deux massifs. Par ailleurs, on voit également l’effet de la différence de perméabilité entre les deux massifs. Conformément à ce qu’on attend, la pression s’uniformise plus rapidement dans le massif supérieur, de perméabilité plus élevée.
Figure 3.2.1: Profils de pression de gaz le long d’une coupe verticale (x=4m)
On effectue les test de non-régression suivants :
\(X(m)\) |
\(Y(m)\) |
Temps (années) |
\(\mathrm{PRE1}(\mathrm{Pa})\) Aster |
4,23 |
2,35 |
200 |
1,1159E+05 |
4,23 |
2,35 |
1000 |
1,1447E+05 |
3,92 |
-2,58 |
200 |
1,0135E+05 |
3,92 |
-2,58 |
1000 |
1,0929E+05 |
4,0 |
0,0 |
200 |
1,1470E+05 |
4,0 |
0,0 |
1000 |
1,1477E+05 |