Code_Aster
®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell
Date:
03/01/05
Author (S):
S. GRANET
Key
:
V7.32.113-A
Page:
1/8
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A
Organization (S):
EDF-R & D/AMA
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
Document: V7.32.113
WTNP113 - Resaturation of an cell
Summary:
The test presented here makes it possible to check the correct operation of the operators used for the resolution of
equations of a flow in unsaturated medium. This test corresponds to test 3.2 of the plan of qualification of the project
Al [bib1].
It represents the evolution of water saturation of the medium around an cell of storage. Two phases are
takings into account, a phase of desaturation at the time of the exploitation of the underground work and a phase of
resaturation after the fill of gallery of the cells.
Code_Aster
®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell
Date:
03/01/05
Author (S):
S. GRANET
Key
:
V7.32.113-A
Page:
2/8
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A
1
Problem of reference
1.1 Geometry
The studied field represents a cut of ground around an cell of storage.
Co-ordinates of the points (m):
To 0
- 500
C 10 - 400
B 10 - 500
D 0 - 400
O 10 - 450
Radius of the cell: 5.6 m
1.2
Properties of material
One gives here only the properties whose solution depends, knowing that the command file
contains other data of material (temperatures,…) who finally do not play any part in
solution of the dealt with problem.
Liquid water
Density (kg.m
- 3
)
Viscosity
1000
1
Parameters
homogenized
Permeability K
Isotherm of sorption
Relative permeability
Porosity
Storage
10
- 18
m
- 2
()
[
]
33
.
0
49
.
1
8
)
10
.
5
.
6
(
1
85
.
0
15
.
0
C
C
P
P
S
+
+
=
()
3
S
P
Kr
C
W
=
0,14
4. 10
- 10
m
- 1
WITH B
C
D
O
Code_Aster
®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell
Date:
03/01/05
Author (S):
S. GRANET
Key
:
V7.32.113-A
Page:
3/8
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A
1.3
Initial conditions
The problem comprises two phases:
·
First a 15 years phase of desaturation corresponding to the exploitation of the work
underground
·
One second phase of resaturation after the fill of the cell corresponding to
exploitation (one initializes the saturation of the cell with 0,7).
The initial conditions are as follows:
For phase 1
·
Cell
C
P
= 9,4.10
7
AP (S = 0,49)
·
Geological barrier
C
P
= 1.10
5
AP (S = 0,999)
For phase 2 (T > 15 years)
·
Cell
C
P
= 3,015.10
7
AP (S = 0,7)
1.4
Boundary conditions
They are expressed on the capillary pressure.
Phase 1:
On [AB]
C
P
= 1.10
5
AP
On [CB] hydraulic Flow no one
On [CD]
C
P
= 1.10
5
AP
On [A01] U [02D] hydraulic Flow no one
On the whole of the cell
C
P
=9,4.10
7
AP (S=0,49).
Phase 2:
On [AB]
C
P
= 1.10
5
AP
On [CB] hydraulic Flow no one
On [CD]
C
P
= 1.10
5
AP
On [AD] hydraulic Flow no one
Code_Aster
®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell
Date:
03/01/05
Author (S):
S. GRANET
Key
:
V7.32.113-A
Page:
4/8
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A
2 Modeling
With
2.1
Characteristics of modeling A
The results presented here result from modeling in plane deformations carried out with
2988 elements TRI3.
2.2 Functionalities
tested
Order Option
AFFE_MODELE
D_PLAN_THHD
DEFI_MATERIAU
THM_LIQU
THM_GAZ
THM_DIFFU
THM_INIT
AFFE_CHAR_MECA DDL_IMPO
PRE1
PRE2
TEMP
STAT_NON_LINE COMP_INCR RELATION KIT_THH
RELATION_KIT
LIQU_GAZ
HYDR_UTIL
CREA_CHAMP OPERATION
EXTR
CREA_CHAMP OPERATION
AFFE
CREA_CHAMP OPERATION
ADZE
STAT_NON_LINE COMP_INCR
RECH_LINEAIRE
RELATION KIT_THH
RELATION_KIT
LIQU_GAZ
HYDR_UTIL
Discretization in time: 34 pitches of time for a 10000 years simulation.
2.3 Results
One presents the profiles of capillary pressure and saturation on a horizontal cut (y = 450m) and
verticals (X = 7m).
Code_Aster
®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell
Date:
03/01/05
Author (S):
S. GRANET
Key
:
V7.32.113-A
Page:
5/8
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A
Appear 2.3-a: Profiles of capillary pressure Y = - 450m
Appear 2.3-b: Profiles of saturation Y = - 450m
Capillary pressures Y = - 450 m
0,00E+00
2,00E+07
4,00E+07
6,00E+07
8,00E+07
1,00E+08
0
1
2
3
4
5
6
7
8
9 10
X (m)
PC
(
P
has
)
7,5 years
15 years
50 years
100 years
200 years
500 years
1000 years
10000 years
Saturation Y = - 450m
0,4
0,5
0,6
0,7
0,8
0,9
1
0
1
2
3
4
5
6
7
8
9
10
X (m)
S
7,5 years
15 years
50 years
100 years
200 years
500 years
1000 years
10000 years
Code_Aster
®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell
Date:
03/01/05
Author (S):
S. GRANET
Key
:
V7.32.113-A
Page:
6/8
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A
Appear 2.3-c: Profiles of capillary pressure X = 7 m
Appear 2.3-d: Profiles of saturation X = 7 m
2.4 Values
tested
X (m)
Y (m)
Time (years)
PRE1 (AP) Aster
5 - 450
27
2.9E7
5 - 450
10000 4.5E6
7 - 450
27
4.2E6
7 - 450
10000 4.5E6
These results are qualitatively in conformity with those which one finds in the literature of Alliances
(isovaleurs).
Capillary pressures X = 7m
0,00E+00
2,00E+06
4,00E+06
6,00E+06
8,00E+06
1,00E+07
- 500
- 480
- 460
- 440
- 420
- 400
Y (m)
PC
7,5 years
15 years
50 years
100 years
200 years
500 years
1000 years
10000 years
Saturations X = 7m
0,8
0,82
0,84
0,86
0,88
0,9
0,92
0,94
0,96
0,98
1
- 500
- 480
- 460
- 440
- 420
- 400
Y (m)
S
7,5 years
15 years
50 years
100 years
200 years
500 years
1000 years
10000 years
Code_Aster
®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell
Date:
03/01/05
Author (S):
S. GRANET
Key
:
V7.32.113-A
Page:
7/8
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A
3 Modeling
B
3.1
Characteristics of modeling B
It acts of the same modeling as above but with modeling THH 2D (with a coefficient
of infinite Henry). The awaited results must thus be exactly the same ones.
3.2 Functionalities
tested
Order Option
AFFE_MODELE
D_PLAN_THH 2D
DEFI_MATERIAU
THM_LIQU
THM_GAZ
THM_VAPE_GAZ
THM_AIR_DISS
THM_DIFFU
THM_INIT
AFFE_CHAR_MECA DDL_IMPO
PRE1
PRE2
TEMP
STAT_NON_LINE COMP_INCR RELATION KIT_THH
RELATION_KIT
LIQU_AD_GAZ_
VAPE
HYDR_UTIL
CREA_CHAMP OPERATION
EXTR
CREA_CHAMP OPERATION
AFFE
CREA_CHAMP OPERATION
ADZE
STAT_NON_LINE COMP_INCR
RECH_LINEAIRE
RELATION KIT_THH
RELATION_KIT
LIQU_AD_GAZ_
VAPE
HYDR_UTIL
Discretization in time: 34 pitches of time for a 10000 years simulation.
3.3 Values
tested
X (m)
Y (m)
Time (years)
PRE1 (AP) Aster
5 - 450
27.2.9E7
5 - 450
10000
4.5E6
7 - 450
27.4.2E6
7 - 450
10000
4.5E6
The values are of course the same ones as previously. The purpose of this modeling, heavier, is not
to correspond to model the THM most complete as required by Alliances.
4 Bibliography
[1]
Project Alliances plan of qualification, notes ANDRA CNT-ASCS 02-075B
Code_Aster
®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell
Date:
03/01/05
Author (S):
S. GRANET
Key
:
V7.32.113-A
Page:
8/8
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A
Intentionally white left page.