Code_Aster
®
Version
7.4
Titrate:
WTNP110 - Orthotropic flow saturated 2D
Date:
20/10/04
Author (S):
S. GRANET
Key
:
V7.32.110 - A
Page:
1/4
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/04/005/A
Organization (S):
EDF-R & D/AMA
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
Document: V7.32.110
WTNP110 - Orthotropic flow 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.
Code_Aster
®
Version
7.4
Titrate:
WTNP110 - Orthotropic flow saturated 2D
Date:
20/10/04
Author (S):
S. GRANET
Key
:
V7.32.110 - A
Page:
2/4
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/04/005/A
1
Problem of reference
1.1 Geometry
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
)
Viscosity
1
1
Parameters
homogenized
Intrinsic permeability (component in X)
Intrinsic permeability (component in Y)
Intrinsic permeability (component in Z)
1
3/4
1
Initial State
Porosity
Pressure of fluid
1
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
(0;0)
D C
WITH B
Code_Aster
®
Version
7.4
Titrate:
WTNP110 - Orthotropic flow saturated 2D
Date:
20/10/04
Author (S):
S. GRANET
Key
:
V7.32.110 - A
Page:
3/4
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/04/005/A
2
Reference solution
The analytical solution in pressure is a polynomial of degree in
X
and
y
, speed is constant and
horizontal:
)
60
;
45
(
)
80
;
45
(
)
,
(
5
,
22
80
45
)
,
(
=
=
+
-
-
=
y
X
K
K
y
X
V
y
X
y
X
P
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
THM_INIT
ELAS
PERMIN_X
PERMIN_Y
PERMIN_Z
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 pitch of time: 1
S.
Code_Aster
®
Version
7.4
Titrate:
WTNP110 - Orthotropic flow saturated 2D
Date:
20/10/04
Author (S):
S. GRANET
Key
:
V7.32.110 - A
Page:
4/4
Manual of Validation
V7.32 booklet: Thermo hydro-mechanical in porous environment unsaturated
HT-66/04/005/A
3.3 Values
tested
Standard node
of
value
Moment (S)
Reference
(analytical)
Aster Difference
(%)
N441
(0.05;0.05)
PRE1 1 16,3
16,25 0,3
%
N689
(- 0.05;- 0.05)
PRE1 1 28,8
28,5 0.96
%
N341
(0;0)
PRE1 1 22,5
22,5 0.
%
Meshs Standard
of
value
Moment (S)
Reference
(analytical)
Aster Difference
(%)
RIGHT FH11X
1 45
45 0
%
RIGHT FH11Y
1 60
60 0
%
4
Summary of the results
Very good agreement with the reference solution.
5 Bibliography
[1]
Project Alliances plan of qualification, notes ANDRA CNT-ASCS 02-075B