background image
Code_Aster
®
Version
7.2
Titre :
WTNA102 - Diffusion d'air dissous (axi)
Date :
14/10/04
Auteur(s) :
S. GRANET, C. CHAVANT
Clé
:
V7.33.102-A
Page :
1/8
Manuel de Validation
Fascicule V7.33 : Thermo-hydro-mécanique en milieu poreux de structures axisymétriques
HT-66/04/005/A
Organisme(s) :
EDF-R&D/AMA















Manuel de Validation
Fascicule V7.33 : Thermo-hydro-mécanique en milieu poreux de structures
axisymétriques
Document : V7.33.102



WTNA102 - Diffusion d'air dissous (axi)




Résumé :

On considère ici un problème à température et saturation constantes. Par des conditions aux limites
appropriées on impose une pression d'eau et une pression de vapeur constantes. Une pression de gaz est
imposée sur un bord du domaine (flux nuls de l'autre côté). Seules les pressions d'air sec et d'air dissous
reliées par la loi de Henry évoluent. Ce problème se ramène en une équation pour la pression d'air sec de type
« équation de la chaleur ». La solution de référence sera alors un calcul thermique ASTER.
background image
Code_Aster
®
Version
7.2
Titre :
WTNA102 - Diffusion d'air dissous (axi)
Date :
14/10/04
Auteur(s) :
S. GRANET, C. CHAVANT
Clé
:
V7.33.102-A
Page :
2/8
Manuel de Validation
Fascicule V7.33 : Thermo-hydro-mécanique en milieu poreux de structures axisymétriques
HT-66/04/005/A
1
Problème de référence
1.1 Géométrie




Coordonnées des points (m) :
A 0 0
C
1 0,5
B
1
0
D
0
0,5
1.2
Propriétés du matériau
On ne donne ici que les propriétés dont la solution dépend, sachant que le fichier de commandes
contient d'autres données de matériau (modules d'élasticité, conductivité thermique ...) qui finalement
ne jouent aucun rôle dans la solution du problème traité.
Eau liquide
Masse volumique (kg.m
-3
)
Chaleur massique à pression constante (J.K
-1
)
Viscosité dynamique de l'eau liquide (Pa.s)
coefficient de dilatation thermique du liquide (K
-1
)
Perméabilité relative à l'eau
10
3
0.
0.001
0.
( )
5
.
0
=
S
kr
w
Vapeur
Chaleur massique (J.K
-1
)
Masse molaire(kg.mol
-1
)
0.
0,01
Gaz
Chaleur massique (J.K
-1
)
Masse molaire(kg.mol
-1
)
Perméabilité relative au gaz
Viscosité du gaz (kg.m
-1
.s
-1
)
0.
0,01
( )
5
.
0
=
S
kr
gz
0.001
Air dissous
Chaleur massique (J.K
-1
)
Constante de Henry (Pa.m
3
.mol
-1
)
0.
50000
Etat initial
Porosité
Température (K)
Pression de gaz (Pa)
Pression de vapeur (Pa)
Pression capillaire (Pa)
Saturation initiale en liquide
1
300
1.01E5
1000
1.
E
6
0,4
Constantes
Constante des gaz parfaits
8,32
Coefficients
homogénéisés
Masse volumique homogénéisée (kg.m
-3
)
Isotherme de sorption
Coefficient de Biot
Fick Vapeur (m
2.
s
-1
)
Fick air dissous (m
2.
s
-1
)
Perméabilité intrinsèque (m
2
)
2200
( )
4
.
0
=
c
P
S
0
FV=0
FA=6 . E-10
Kint = 1.E-19
y
x
A
B
C
D
background image
Code_Aster
®
Version
7.2
Titre :
WTNA102 - Diffusion d'air dissous (axi)
Date :
14/10/04
Auteur(s) :
S. GRANET, C. CHAVANT
Clé
:
V7.33.102-A
Page :
3/8
Manuel de Validation
Fascicule V7.33 : Thermo-hydro-mécanique en milieu poreux de structures axisymétriques
HT-66/04/005/A
1.3
Conditions aux limites et chargements
Sur l'ensemble du domaine, on veut :
1
)
(
0
0
1
0
0
0
0
0
=
=
=
=
=
=
=
=
=
=
=
=
=
T
cte
T
S
cte
p
S
F
p
cte
p
cte
K
p
cte
p
c
vp
vp
vp
w
w
w
w
w
ol
ad
ol
vp
ol
as
M
M
M
=
=
Sur tous les bords : Flux hydrauliques et thermiques nuls.
On va maintenant linéariser
vp
p
en fonction de
w
p
.

Ecriture de p
vp
fonction linéaire de p
w
:
La section 4.2.3 du document de référence Aster [R7.01.11] nous donne la relation
:
w
w
ol
vp
vp
vp
dp
RT
M
p
dp
=
. Si on linéarise cette expression on obtient
:




-
+
=
0
0
0
0
0
0
w
w
ol
vp
vp
vp
w
w
ol
vp
vp
vp
p
M
RT
p
p
p
M
RT
p
p
que l'on peut écrire sous la forme :
B
Ap
p
w
vp
+
=
éq 1.3-1
avec
0
0
w
ol
vp
vp
M
RT
p
A
=
et
0
0
0
0
w
w
ol
vp
vp
vp
p
M
RT
p
p
B
-
=
6
10
115000
E
pc
p
gz
=
=
:
AB
bord
le
Sur

background image
Code_Aster
®
Version
7.2
Titre :
WTNA102 - Diffusion d'air dissous (axi)
Date :
14/10/04
Auteur(s) :
S. GRANET, C. CHAVANT
Clé
:
V7.33.102-A
Page :
4/8
Manuel de Validation
Fascicule V7.33 : Thermo-hydro-mécanique en milieu poreux de structures axisymétriques
HT-66/04/005/A
2
Solution de référence
2.1
Méthode de calcul
2.1.1 Calcul de la conservation de la masse d'air
La conservation de la masse de gaz s'écrit :
0
)
(
=
+
+
ad
as
air
div
dt
dm
M
M
éq
2.1.1-1
On écrit que la masse totale d'eau et la masse totale d'air sont conservée (car il n'y a pas de flux
d'eau ni de gaz au bord) et on obtient :
)
)(
1
(
)
(
0
0
0
0
as
as
ad
ad
ad
as
air
S
S
m
m
m
-
-
+
-
=
+
=
donc
as
ad
ad
as
d
S
d
S
m
m
d
)
1
(
)
(
0
0
-
+
=
+
éq
2.1.1-2
as
ol
as
as
dP
RT
M
d
=
et
as
H
ol
ad
ad
dP
K
M
d
=
dt
dP
RT
M
S
K
M
S
dt
dm
as
ol
as
H
ol
as
air




-
+
=
)
1
(
0
.
0
Calcul des vitesses :
)
(
as
gz
as
P
-
=
as
M
éq
2.1.1-3
puisque
0
=
vp
F
et
0
=
vp
P
et
ad
ad
lq
lq
ad
d
C
F
P
-
-
=
)
(
a
M
avec
ad
ad
C
=
Comme
as
H
ad
ad
w
lq
P
K
RT
P
P
P
P
=
=
+
=
as
ad
H
ol
ad
as
H
lq
ad
d
P
F
K
M
P
K
RT
-
-
=
.
.
)
(
a
M
[éq 2.1.1-1] peut alors se simplifier sous la forme suivante :
)
(
as
as
P
Ldiv
dt
dP
C
=
avec




+
+
=
-
+
=
ad
H
ol
as
lq
ad
H
gz
as
ol
as
H
ol
as
F
K
M
K
RT
L
RT
M
S
K
M
S
C
.
)
1
(
0
0
0
.
0
Equation de la chaleur dont on connaît le résultat.
background image
Code_Aster
®
Version
7.2
Titre :
WTNA102 - Diffusion d'air dissous (axi)
Date :
14/10/04
Auteur(s) :
S. GRANET, C. CHAVANT
Clé
:
V7.33.102-A
Page :
5/8
Manuel de Validation
Fascicule V7.33 : Thermo-hydro-mécanique en milieu poreux de structures axisymétriques
HT-66/04/005/A
2.2
Résultats de référence
Avec les valeurs numériques précédentes, on trouve :
4992
10
0
0
5
=
=
=
as
H
ad
as
P
K
RT
P
P
4
.
0
0
0
=
=
as
ol
as
as
P
RT
M
et
02
.
0
0
0
=
=
ad
ol
ad
ad
P
RT
M
3
0
10
.
4
-
=
=
vp
vp
Les constantes de l'équation de la chaleur sont alors :
16
6
10
.
4
,
1
4810
,
2
-
-
=
=
L
C

2.3 Incertitudes
Les incertitudes sont assez grandes parce-que la solution analytique est une solution approchée du
fait de la linéarisation des équations.
background image
Code_Aster
®
Version
7.2
Titre :
WTNA102 - Diffusion d'air dissous (axi)
Date :
14/10/04
Auteur(s) :
S. GRANET, C. CHAVANT
Clé
:
V7.33.102-A
Page :
6/8
Manuel de Validation
Fascicule V7.33 : Thermo-hydro-mécanique en milieu poreux de structures axisymétriques
HT-66/04/005/A
3 Modélisation
A
3.1
Caractéristiques de la modélisation A
Modélisation en déformations planes. 20 éléments QUAD8.
3.2 Fonctionnalités
testées
Commande Option
AFFE_MODELE
AXIS_THH2D
DEFI_MATERIAU
THM_LIQU
THM_GAZ
THM_VAPE_GAZ
THM_AIR_DISS
THM_DIFFU
THM_INIT
ELAS
AFFE_CHAR_MECA DDL_IMPO
PRE1
PRE2
TEMP
STAT_NON_LINE COMP_INCR
RELATION KIT_THH
RELATION_KIT
ELAS
LIQU_AD_GAZ_VAPE
HYDR_UTIL
Discrétisation en temps : 100 pas de temps de 5E7 s chacun.
3.3 Résultats
X (m)
Temps (s)
PRE2 Aster
PRE2 calcul thermique
Erreur relative
0,2 3E9s 7.90E3
7.94E3 0.56%
0,2 5E9s 9.50E3
9.56E3 0.60%
Comparaison pression d'air sec, calcul thermique
9,80E+04
1,00E+05
1,02E+05
1,04E+05
1,06E+05
1,08E+05
1,10E+05
1,12E+05
1,14E+05
0
0,5
1
X
Pa
s
Pas; t=1E9
Pas; t=1E10s
Pas; t=2E10s
TEMP; t=2E9s
TEMP; t=1E10s
TEMP; t=2E10s
background image
Code_Aster
®
Version
7.2
Titre :
WTNA102 - Diffusion d'air dissous (axi)
Date :
14/10/04
Auteur(s) :
S. GRANET, C. CHAVANT
Clé
:
V7.33.102-A
Page :
7/8
Manuel de Validation
Fascicule V7.33 : Thermo-hydro-mécanique en milieu poreux de structures axisymétriques
HT-66/04/005/A
4
Synthèse des résultats
Les résultats Aster sont en très bon accord avec la solution analytique.
background image
Code_Aster
®
Version
7.2
Titre :
WTNA102 - Diffusion d'air dissous (axi)
Date :
14/10/04
Auteur(s) :
S. GRANET, C. CHAVANT
Clé
:
V7.33.102-A
Page :
8/8
Manuel de Validation
Fascicule V7.33 : Thermo-hydro-mécanique en milieu poreux de structures axisymétriques
HT-66/04/005/A


























Page laissée intentionnellement blanche.