Code_Aster
®
Version
8.2
Titre :
WTNV127 Désaturation d'un milieu poreux sans air (3D)
Date :
04/05/06
Auteur(s) :
S. GRANET
Clé
:
V7.31.127-B
Page :
1/8
Manuel de Validation
Fascicule V7.31 : Thermo-hydro-mécanique en milieu poreux de structures volumiques
HT-62/06/005/A
Organisme(s) :
EDF-R&D/AMA
Manuel de Validation
Fascicule V7.31: Thermo-hydro-mécanique en milieu poreux de structures volumiques
Document : V7.31.127
WTNV127 Désaturation d'un milieu poreux sans
air (modélisation 3D_THV)
Résumé :
On chauffe un milieu poreux dont les pores sont remplis d'un mélange d'eau et de vapeur d'eau. La saturation
initiale en liquide est de 50%, le chargement est un flux thermique uniforme sur les bords du domaine. La
modélisation faite par un seul élément cubique correspond à la modélisation d'un problème homogène en
espace.
La solution de référence est une solution analytique approchée.
Code_Aster
®
Version
8.2
Titre :
WTNV127 Désaturation d'un milieu poreux sans air (3D)
Date :
04/05/06
Auteur(s) :
S. GRANET
Clé
:
V7.31.127-B
Page :
2/8
Manuel de Validation
Fascicule V7.31 : Thermo-hydro-mécanique en milieu poreux de structures volumiques
HT-62/06/005/A
1
Problème de référence
1.1 Géométrie
Coordonnées des points (m) :
A
0
0
100
E
0
0
0
B 100
0 100 F 100
0 0
C 100 100 100 G 100 100 0
D
0 100 100 H
0 100 0
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 à pression constante (J.K
-1
)
coefficient de dilatation thermique du liquide (K
-1
)
10
3
4180
0.
Vapeur
Capacité calorifique (J.K
-1
)
Enthalpie initiale (chaleur latente de vaporisation)
Masse molaire(kg.mol
-1
)
1900
2,5E6.
0,018
Squelette
Capacité calorifique à contrainte constante(J.K
-1
) 1050
Etat initial
Porosité
Température (K)
Pression de liquide (Pa)
Pression de vapeur (Pa)
Saturation initiale en liquide
0,3
300
1E5
3700
0,5
Constantes
Constante des gaz parfaits
8,315
Coefficients
homogénéisés
Masse volumique homogénéisée (kg.m
-3
)
Isotherme de sorption
2200
( )
(
)
0
0
12
10
5
.
0
c
vp
c
c
P
P
P
P
S
-
-
-
=
-
Avec
3700
0
=
vp
P
5
0
10
-
=
c
P
D
C
A B
H
G
F
E
Code_Aster
®
Version
8.2
Titre :
WTNV127 Désaturation d'un milieu poreux sans air (3D)
Date :
04/05/06
Auteur(s) :
S. GRANET
Clé
:
V7.31.127-B
Page :
3/8
Manuel de Validation
Fascicule V7.31 : Thermo-hydro-mécanique en milieu poreux de structures volumiques
HT-62/06/005/A
1.3
Conditions aux limites et chargements
Sur toutes les faces :
Flux thermique
6
10
.
=
n
q
ext
Flux hydraulique nul
2
Solution de référence
2.1
Méthode de calcul
2.1.1 Calcul de la pression de vapeur à partir de la température
Nous supposons la courbe de saturation linéaire. Elle s'écrit donc :
c
P
S
S
S
+
=
0
éq
2.1.1-1
[R7.01.11] équation [éq 3.2.1-2] donnent alors :
(
)
(
)
c
vp
vp
vp
vp
lq
lq
c
lq
lq
P
S
S
m
S
P
S
m
-
-
-
=
-
=
0
0
0
0
0
0
0
0
1
éq
2.1.1-2
On écrit que la masse totale d'eau est conservée (car il n'y a pas de flux d'eau au bord) et on obtient :
(
)
(
)
(
)
0
1
0
0
0
=
-
-
+
-
=
+
S
P
S
m
m
vp
vp
c
vp
lq
vp
lq
éq
2.1.1-3
[R7.01.11] équation [éq 4.4-1] donne par ailleurs
(
)
(
)
-
+
-
+
-
-
+
=
1
ln
1
1
1
ln
0
0
0
0
0
0
T
T
T
T
C
C
R
M
T
T
h
h
R
M
P
RT
M
p
p
p
lq
p
vp
ol
vp
lq
vp
ol
vp
lq
ol
vp
vp
vp
lq
éq
2.1.1-4
Le couplage des équations [éq 2.1.1-3] et [éq 2.1.1-4], auquel il faut ajouter l'équation des gaz parfaits
pour la vapeur, est un système fortement non linéaire que nous résoudrons en petite perturbation, ce
qui permet de le linéariser.
Code_Aster
®
Version
8.2
Titre :
WTNV127 Désaturation d'un milieu poreux sans air (3D)
Date :
04/05/06
Auteur(s) :
S. GRANET
Clé
:
V7.31.127-B
Page :
4/8
Manuel de Validation
Fascicule V7.31 : Thermo-hydro-mécanique en milieu poreux de structures volumiques
HT-62/06/005/A
Tous calculs faits, on obtient :
(
)
(
)
(
)
(
)
(
)
-
=
-
-
=
-
-
-
+
-
2
0
0
0
0
0
2
0
0
0
0
0
0
0
1
1
T
T
h
h
R
M
P
RT
M
p
P
T
T
R
M
p
S
P
S
RT
M
S
S
P
lq
vp
ol
vp
lq
lq
ol
vp
vp
vp
ol
vp
vp
lq
vp
lq
ol
vp
vp
lq
vp
éq 2.1.1-5
2.1.2 Calcul de la température
[R7.01.11] équation [éq 3.2.4.3-1] donne :
T
C
p
T
Q
vp
m
gz
+
-
=
0
3
éq
2.1.2-1
(puisque les autres coefficients de dilatation sont nuls).
Equation [éq 3.2.4.3-2] donne :
(
)
T
S
lq
m
gz
3
1
-
=
éq
2.1.2-2
On obtient donc :
(
)
T
C
p
S
Q
vp
lq
+
-
-
=
0
1
éq
2.1.2-3
Dans ce problème,
Q
n'est rien d'autre que la chaleur apportée par unité de volume.
En appelant
Vol
le volume total de la pièce et
Surf
sa surface latérale et
t
le temps d'application
des flux :
n
q .
ext
Vol
Surf
t
Q
=
éq
2.1.2-4
2.1.3 Système à résoudre
(
)
(
)
(
)
(
)
(
)
(
)
=
-
-
-
-
-
-
-
-
-
-
+
-
n
q .
0
0
1
0
1
1
1
0
2
0
0
0
0
0
2
0
0
0
0
0
0
0
ext
lq
vp
lq
lq
vp
ol
vp
lq
ol
vp
vp
ol
vp
vp
vp
lq
ol
vp
vp
lq
Vol
Surf
t
T
P
P
C
S
T
h
h
R
M
RT
M
p
RT
M
p
S
S
RT
M
S
S
éq 2.1.2-5
Code_Aster
®
Version
8.2
Titre :
WTNV127 Désaturation d'un milieu poreux sans air (3D)
Date :
04/05/06
Auteur(s) :
S. GRANET
Clé
:
V7.31.127-B
Page :
5/8
Manuel de Validation
Fascicule V7.31 : Thermo-hydro-mécanique en milieu poreux de structures volumiques
HT-62/06/005/A
2.2
Résultats de référence
On donne la valeur de la température, de la pression de liquide et de la pression de vapeur, solution
du système (10) avec les données résumées au paragraphe [§1.2] et rappelées ci dessous. Pour le
calcul des capacités calorifiques, on utilise les relations suivantes :
(
)
(
)
0
0
0
0
0
0
0
0
1
1
vp
l
l
lq
s
S
S
r
-
-
-
=
-
(
)
(
)
p
vp
vp
l
p
lq
l
lq
s
s
C
S
C
S
C
C
-
+
+
-
=
1
1
0
0
0
= C
C
, cette dernière relation étant vraie parce que le coefficient de dilatation des grains est nul.
0
S
S
0
T
0
vp
p
0
vp
h
0
vp
(calculé)
lq
5,00E-01 -1,00E-12 3,00E+02
3,70E+03
2,50E+06
2,67E-02 1,00E+03
0
r
0
s
(calculé)
s
C
p
lq
C
l
p
vp
C
0
C
(calculé)
2,20E+03 3,00E-01 2,93E+03
1,05E+03
4,18E+03
1,90E+03 2,78E+06
n
q
.
ext
t
Surf
Vol
1,00E+06 100
6,00E+04
1,00E+06
Après résolution, on obtient les résultats suivants :
vp
P
425
lq
P
l
-1,4E+06
T
2E+00
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.
Code_Aster
®
Version
8.2
Titre :
WTNV127 Désaturation d'un milieu poreux sans air (3D)
Date :
04/05/06
Auteur(s) :
S. GRANET
Clé
:
V7.31.127-B
Page :
6/8
Manuel de Validation
Fascicule V7.31 : Thermo-hydro-mécanique en milieu poreux de structures volumiques
HT-62/06/005/A
3 Modélisation
A
3.1
Caractéristiques de la modélisation A
Modélisation en déformations planes. Un élément Q8.
3.2 Fonctionnalités
testées
Commande Option
AFFE_MODELE
3D_THVD
DEFI_MATERIAU
THM_LIQU
THM_DIFFU
THM_INIT
ELAS
AFFE_CHAR_MECA DDL_IMPO PRE1
TEMP
STAT_NON_LINE COMP_INCR
RELATION KIT_THV
RELATION_KIT
LIQU_VAPE
HYDR_UTIL
Discrétisation en temps : un seul pas de temps : 10
2
s.
3.3 Valeurs
testées
Noeud
Type de valeur
Instant (s)
Référence
(analytique)
Aster Différence
(%)
NO1
DEPL/TEMP
10
2
2 2,15 8%
NO1
DEPL/PRE1
10
2
-1,4
10
6
-1,46
10
6
4%
NO1
VARI_ELNO_ELGA/V3
10
2
425 460
8%
On trouve donc des résultats relativement proches des résultats analytiques. L'incertitude demeurant
assez large du fait de la linéarisation des équations.
Code_Aster
®
Version
8.2
Titre :
WTNV127 Désaturation d'un milieu poreux sans air (3D)
Date :
04/05/06
Auteur(s) :
S. GRANET
Clé
:
V7.31.127-B
Page :
7/8
Manuel de Validation
Fascicule V7.31 : Thermo-hydro-mécanique en milieu poreux de structures volumiques
HT-62/06/005/A
4 Modélisation
B
Même modélisation mais en sélective.
4.1
Caractéristiques de la modélisation B
Modélisation en déformations planes. Un élément Q8
4.2 Fonctionnalités
testées
Commande
Option
AFFE_MODELE
3D_THVS
DEFI_MATERIAU
THM_LIQU
THM_DIFFU
THM_INIT
ELAS
AFFE_CHAR_MECA DDL_IMPO PRE1
TEMP
STAT_NON_LINE COMP_INCR
RELATION KIT_THV
RELATION_KIT
LIQU_VAPE
HYDR_UTIL
Discrétisation en temps : un seul pas de temps : 10
2
s.
4.3 Valeurs
testées
Noeud
Type de valeur
Instant (s)
Référence
(analytique)
Aster Différence
(%)
NO1
DEPL/TEMP
10
2
2
2,15 8%
NO1
DEPL/PRE1
10
2
-1,4
10
6
-1,46
10
6
4%
NO1
VARI_ELNO_ELGA/V3
10
2
425
460
8%
On trouve donc des résultats relativement proches des résultats analytiques. L'incertitude demeurant
assez large du fait de la linéarisation des équations.
Code_Aster
®
Version
8.2
Titre :
WTNV127 Désaturation d'un milieu poreux sans air (3D)
Date :
04/05/06
Auteur(s) :
S. GRANET
Clé
:
V7.31.127-B
Page :
8/8
Manuel de Validation
Fascicule V7.31 : Thermo-hydro-mécanique en milieu poreux de structures volumiques
HT-62/06/005/A
Page laissée intentionnellement blanche.