Code_Aster ®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell


Date:
03/01/05
Author (S):
Key S. GRANET
:
V7.32.113-A Page:
1/8

Organization (S): EDF-R & D/AMA

Handbook 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.

Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A

Code_Aster ®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell


Date:
03/01/05
Author (S):
Key S. GRANET
:
V7.32.113-A Page:
2/8

1
Problem of reference

1.1 Geometry

The studied field represents a cut of ground around an cell of storage.

D
C

O

WITH B



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)
1000
Viscosity
1
Parameters
Permeability K
10-18 m2
homogenized
Isotherm of sorption
S (P =
+
c)
0 15
.

0 85
.


[
8
1 + (6 5
. 10
.
P)
C
] 33.0
49
.
1
Relative permeability
Kr
=

W (PC)
3
S

Porosity
0,14
Storage
4. 10-10m-1
Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A

Code_Aster ®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell


Date:
03/01/05
Author (S):
Key S. GRANET
:
V7.32.113-A Page:
3/8

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 P = 9,4.107 Pa (S = 0,49)
C

· Geological barrier P = 1.105 Pa (S = 0,999)
C

For phase 2 (T > 15 years)

· Cell P = 3,015.107 Pa (S = 0,7)
C

1.4
Boundary conditions

They are expressed on the capillary pressure.

Phase 1:

On [AB] P = 1.105 Pa
C
On [CB] hydraulic Flux no one
On [CD] P = 1.105 Pa
C
On [A01] U [02D] hydraulic Flux no one

On the whole of the cell P =9,4.107 Pa (S=0,49).
C

Phase 2:

On [AB] P = 1.105 Pa
C
On [CB] hydraulic Flux no one
On [CD] P = 1.105 Pa
C
On [AD] hydraulic Flux no one

Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A

Code_Aster ®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell


Date:
03/01/05
Author (S):
Key S. GRANET
:
V7.32.113-A Page:
4/8

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
ASSE



STAT_NON_LINE COMP_INCR RELATION KIT_THH
RECH_LINEAIRE

RELATION_KIT
LIQU_GAZ

HYDR_UTIL

Discretization in time: 34 steps 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).

Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A

Code_Aster ®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell


Date:
03/01/05
Author (S):
Key S. GRANET
:
V7.32.113-A Page:
5/8

Capillary pressures Y = - 450 m

1,00E+08

8,00E+07
7,5 years

15 years

has
) 6,00E+07

P
50 years

(
100 years

4,00E+07
PC

200 years

2,00E+07
500 years

0,00E+00
1000 years

0
1
2
3
4
5
6
7
8
9 10
10000 years

X (m)

Appear 2.3-a: capillary Profils of pressure Y = - 450m

Saturation Y = - 450m


1
7,5 years

15 years

0,9

50 years

0,8
100 years

200 years

S 0,7

500 years

0,6
1000 years

0,5
10000 years

0,4

0
1
2
3
4
5
6
7
8
9
10

X (m)

Appear 2.3-b: Profils of saturation Y = - 450m

Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A

Code_Aster ®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell


Date:
03/01/05
Author (S):
Key S. GRANET
:
V7.32.113-A Page:
6/8

Capillary pressures X = 7m

1,00E+07

7,5 years

8,00E+06

15 years

6,00E+06

50 years

PC
100 years

4,00E+06

200 years

2,00E+06

500 years

0,00E+00
1000 years

- 500
- 480
- 460
- 440
- 420
- 400
10000 years

Y (m)

Appear 2.3-c: Profiles of capillary pressure X = 7 m

Saturations X = 7m

1

0,98
7,5 years

0,96
15 years

0,94
50 years

0,92
100 years

S
0,9

0,88
200 years

500 years

0,86

0,84
1000 years

0,82

10000 years

0,8

- 500
- 480
- 460
- 440
- 420
- 400

Y (m)

Appear 2.3-d: Profiles of saturation X = 7 m

2.4 Values
tested

X (m)
Y (m)
Time (years)
PRE1 (Pa) 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).
Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A

Code_Aster ®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell


Date:
03/01/05
Author (S):
Key S. GRANET
:
V7.32.113-A Page:
7/8

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
ASSE



STAT_NON_LINE COMP_INCR RELATION KIT_THH
RECH_LINEAIRE

RELATION_KIT
LIQU_AD_GAZ_

VAPE
HYDR_UTIL

Discretization in time: 34 steps of time for a 10000 years simulation.

3.3 Values
tested

X (m)
Y (m)
Time (years)
PRE1 (Pa) 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]
Plane Alliances project of qualification, notes ANDRA CNT-ASCS 02-075B
Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A

Code_Aster ®
Version
7.3
Titrate:
WTNP113 - Resaturation of an cell


Date:
03/01/05
Author (S):
Key S. GRANET
:
V7.32.113-A Page:
8/8

Intentionally white left page.
Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/05/005/A

Outline document