Code_Aster ®
Version
7.4
Titrate:
WTNP110 - Orthotropic Ecoulement saturated 2D
Date:
20/10/04
Author (S):
Key S. GRANET
:
V7.32.110 - A Page:
1/4
Organization (S): EDF-R & D/AMA
Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
Document: V7.32.110
WTNP110 - Orthotropic Ecoulement saturated 2D
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 orthotropic saturated medium. This test corresponds to test 1.2 of the plan of qualification
project Al [bib1].
The reference solution is an analytical solution.
Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/04/005/A
Code_Aster ®
Version
7.4
Titrate:
WTNP110 - Orthotropic Ecoulement saturated 2D
Date:
20/10/04
Author (S):
Key S. GRANET
:
V7.32.110 - A Page:
2/4
1
Problem of reference
1.1 Geometry
D C
(0;0)
WITH B
Co-ordinates of the points (m):
With - 0,1
- 0,1
C 0,1 0,1
B 0,1 - 0,1
D - 0,1 0,1
1.2
Properties of material
One gives here only the properties whose solution depends, knowing that the command file
contains other data of material (moduli of elasticity,…) who finally do not play any part
in the solution of the dealt with problem.
Liquid water
Density (kg.m-3)
1
Viscosity
1
Parameters
Intrinsic permeability (component in X)
1
homogenized
Intrinsic permeability (component in Y)
3/4
Intrinsic permeability (component in Z)
1
Initial State
Porosity
1
Pressure of liquid
0
1.3
Boundary conditions and loadings
1.3.1 Modeling
With
On AB P = P (X) = - 45X+30.5
On BC P = P (X) = - 80Y+18
On CD P = P (X) = - 45X+14.5
On DA P = P (X) = - 80Y+27
Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/04/005/A
Code_Aster ®
Version
7.4
Titrate:
WTNP110 - Orthotropic Ecoulement saturated 2D
Date:
20/10/04
Author (S):
Key S. GRANET
:
V7.32.110 - A Page:
3/4
2
Reference solution
The analytical solution in pressure is a polynomial of degree in X and y, speed is constant and
horizontal:
P (X, y) = - 45x - 80y +
5
,
22
V (X, y) = (45K 80
; K) = (
)
60
;
45
X
y
2.1
Results of reference
One gives the value of the pressure and speed in 3 points:
X - 0,05 0,00 +0,05
y - 0,05 0.00 0.05
P 28,8 22.5 16.3
Vx 45 45 45
Vy 60 60 60
3 Modeling
With
3.1
Characteristics of modeling A
Modeling in plane deformations. 20 * 20 Q8 elements
3.2 Functionalities
tested
Order Option
AFFE_MODELE
D_PLAN_HMD
DEFI_MATERIAU
THM_LIQU
THM_DIFFU
PERMIN_X
PERMIN_Y
PERMIN_Z
THM_INIT
ELAS
AFFE_CHAR_MECA DDL_IMPO
PRE1
DX
DY
STAT_NON_LINE COMP_INCR RELATION
KIT_HM
RELATION_KIT
ELAS
LIQU_SATU
HYDR_UTIL
Discretization in time: only one step of time: 1 S.
Handbook of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/04/005/A
Code_Aster ®
Version
7.4
Titrate:
WTNP110 - Orthotropic Ecoulement saturated 2D
Date:
20/10/04
Author (S):
Key S. GRANET
:
V7.32.110 - A Page:
4/4
3.3 Values
tested
Standard node
of
Moment (S)
Reference
Aster Difference
value
(analytical)
(%)
N441
PRE1 1 16,3
16,25 0,3
%
(0.05;0.05)
N689
PRE1 1 28,8
28,5 0.96
%
(- 0.05;- 0.05)
N341
PRE1 1 22,5
22,5 0.
%
(0;0)
Meshs Standard
of
Moment (S)
Reference
Aster Difference
value
(analytical)
(%)
DROIT FH11X
1 45
45 0
%
DROIT FH11Y
1 60
60 0
%
4
Summary of the results
Very good agreement with the reference solution.
5 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/04/005/A
Outline document