1. Reference problem#
1.1. Geometry#
The field studied consists of two porous masses (\(\mathrm{AEFD}\) and \(\mathrm{EBCF}\)) separated by a hydraulic joint. The play between the structure and the gas injection wall (\(\mathrm{AGIB}\)) is also modelled.

Point coordinates (in meters):
1.2. Material properties#
Properties of intersticial fluid (dihydrogen):
Molar mass |
\(\mathrm{0,002}{\mathrm{kg.m}}^{-3}\) |
Viscosity |
\({9.10}^{-6}\mathrm{Pa.s}\) |
Properties of the upper rock matrix:
The matrix is elastic and has the following properties:
Young’s module |
\(\mathrm{3,0}\mathrm{GPa}\) |
Poisson’s Ratio |
\(\mathrm{0,12}\) |
Porosity |
\(\mathrm{0,18}\) |
Intrinsic permeability |
\(\mathrm{2,75}{.10}^{-20}{m}^{2}\) |
Properties of the lower rock matrix:
The lower rock matrix has the same mechanical characteristics as the upper matrix, but has an intrinsic permeability that is ten times lower.
Young’s module |
\(\mathrm{3,0}\mathrm{GPa}\) |
Poisson’s Ratio |
\(\mathrm{0,12}\) |
Porosity |
\(\mathrm{0,18}\) |
Intrinsic permeability |
\(\mathrm{2,75}{.10}^{-21}{m}^{2}\) |
Properties of the discontinuity:
The mechanical behavior of discontinuity is described by Bandis’s law. Hydraulic flow is given by cubic law.
Initial normal stiffness |
\({1.10}^{9}{\mathrm{Pa.m}}^{-1}\) |
Initial asymptotic opening |
\(\mathrm{0,4}\mathrm{mm}\) |
Coefficient \(\gamma\) |
|
Game Properties:
The play between the emission wall of the dihydrogen flow and the rock is an elastic medium that is very inrigid and has very high permeability.
Young’s module |
\(\mathrm{3,0}\mathrm{MPa}\) |
Poisson’s Ratio |
\(\mathrm{0,12}\) |
Porosity |
\(1\) |
Intrinsic permeability |
\({1.10}^{-8}{m}^{2}\) |
1.3. Boundary conditions and loading#
The hydraulic boundary conditions are as follows:
On \([\mathrm{GD}]\) zero hydraulic flow
On \([\mathrm{DC}]\) \({p}_{g}={p}_{0}=\mathrm{0,1}\mathrm{MPa}\)
On \([\mathrm{CI}]\) zero hydraulic flow
On \([\mathrm{IG}]\) gas flow \({F}_{g}={1.10}^{-10}\mathrm{kg.}{s}^{-1}\mathrm{.}{m}^{-2}\)
The mechanical boundary conditions are given by the figure below.

1.4. Initial conditions#
The initial conditions are as follows:
initial opening: \(\mathrm{3,48}{.10}^{-6}m\)
initial pressure in the massif: \({p}_{0}=\mathrm{0,1}\mathrm{MPa}\)
compressive stress in direction \(y\): \(\mathrm{12,3}\mathrm{MPa}\)
temperature: \(303°K\)