background image
Code_Aster
®
Version
6.4
Titre :
TPNA01 - Problème axisymétrique stationnaire avec rayonnement
Date
:
02/06/03
Auteur(s) :
J.P. LEFEBVRE
Clé
:
V4.41.001-C
Page :
1/8
Manuel de Validation
Fascicule V4.41 : Thermique stationnaire avec rayonnement
HT-66/03/008/A
Organisme(s) :
EDF-R&D/AMA













Manuel de Validation
Fascicule V4.41 : Thermique stationnaire avec rayonnement
Document : V4.41.001



TPNA01 - Problème axisymétrique stationnaire
avec rayonnement




Résumé :

Ce test élémentaire permet de traiter un problème axisymétrique en thermique stationnaire avec une condition
aux limites de type rayonnement. La solution est analytique. Le problème est traité en axisymétrique et en
volumique.

Pour les modélisations présentées ici, les écarts des résultats obtenus par le Code_Aster se situent entre 1 et
2% de la référence calculée analytiquement.
background image
Code_Aster
®
Version
6.4
Titre :
TPNA01 - Problème axisymétrique stationnaire avec rayonnement
Date
:
02/06/03
Auteur(s) :
J.P. LEFEBVRE
Clé
:
V4.41.001-C
Page :
2/8
Manuel de Validation
Fascicule V4.41 : Thermique stationnaire avec rayonnement
HT-66/03/008/A
1
Problème de référence
1.1 Géométrie
z
Re = 0.300 m
Re = 0.391 m
r
Le cylindre creux est supposé infiniment long.
1.2
Propriétés de matériaux
Seul le coefficient de conductivité intervient. Le Code_Aster rend obligatoire la fourniture d'une fonction
représentant l'enthalpie volumique déterminée à partir du coefficient de chaleur volumique.
chaleur volumique
1.00 J / m
C
conductivité thermique
40 W / m C
C
p
=
=
°
°
3
k
1.3
Conditions aux limites et chargements
Condition de type rayonnement sur la surface interne du cylindre, condition de type convection
(échange avec le milieu extérieur) sur la surface externe.
Pas de condition aux limites sur les extrémités du cylindre (ce qui revient à imposer un flux nul).
(
)
(
)
[
]
C
20.0
et
C
142.0
avec
externe
surface
Celsius
en
exprimée
C,
500.0
et
/
10
5.73
,
0.6
avec
273.15
273.15
interne
surface
2
4
2
8
4
4
°
°
°
°
=
=
-
=
=
=
=


-
+
=
-
+
e
ext
e
e
ext
e
iext
iext
T
W/m
h
T
T
h
n
T
k
T
T
K
m
W
T
T
T
k
n
background image
Code_Aster
®
Version
6.4
Titre :
TPNA01 - Problème axisymétrique stationnaire avec rayonnement
Date
:
02/06/03
Auteur(s) :
J.P. LEFEBVRE
Clé
:
V4.41.001-C
Page :
3/8
Manuel de Validation
Fascicule V4.41 : Thermique stationnaire avec rayonnement
HT-66/03/008/A
2
Solution de référence
2.1
Méthode de calcul utilisée pour la solution de référence
On dispose d'une solution analytique dans le cas d'un cylindre de longueur infinie :
Log(Ri)
-
Log(Re)
Log(Ri)
T
-
Log(Re)
T
r
Log
Log
T
T
T(r)
e
i
i
e
+
)
(


-
=
Ri
Re

2.2
Résultats de référence
r (m)
T (°C)
.30000
105.55
.32275
99.21
.34550
93.30
.36825
87.76
.39100
82.56
Valeur de la température en fonction de r

r (m)
(W / m )
2
.300
11577.49
.391
8822.98
Valeur du flux en fonction de r

2.3
Incertitude sur la solution
Solution exacte.

2.4 Références
bibliographiques
[1]
Guide de validation des progiciels de calcul de structures. Société Française des Mécaniciens
AFNOR 1990 ISBN 2-12-486611-7
background image
Code_Aster
®
Version
6.4
Titre :
TPNA01 - Problème axisymétrique stationnaire avec rayonnement
Date
:
02/06/03
Auteur(s) :
J.P. LEFEBVRE
Clé
:
V4.41.001-C
Page :
4/8
Manuel de Validation
Fascicule V4.41 : Thermique stationnaire avec rayonnement
HT-66/03/008/A
3 Modélisation
A
3.1
Caractéristiques de la modélisation
Modélisation 2D :
N1
N3
N6
N11
N9

3.2
Caractéristiques du maillage
2 QUAD8

3.3
Fonctionnalités testées

Commande Mot-clé
facteur
Mot-clé
simple
Argument
DEFI_MATERIAU THER_NL
LAMBDA
BETA
THER_NON_LINE TEMP_INIT
STATIONNAIRE 'OUI'
CONVERGENCE
RESI_GLOB_RELA
1.E-2
ITER_GLOB_MAXI
9
CRIT_LAGR_RELA
1.E-3
OPTION
'FLUX_ELNO_TEMP'
background image
Code_Aster
®
Version
6.4
Titre :
TPNA01 - Problème axisymétrique stationnaire avec rayonnement
Date
:
02/06/03
Auteur(s) :
J.P. LEFEBVRE
Clé
:
V4.41.001-C
Page :
5/8
Manuel de Validation
Fascicule V4.41 : Thermique stationnaire avec rayonnement
HT-66/03/008/A
4
Résultats de la modélisation A
4.1 Valeurs
testées
Les noeuds observés ont pour coordonnée z = 0.0
Identification
température
Référence
Aster %
différence
N1 (r=.30000)
105.55
105.52
­0.03
N3 (r=.32275)
99.21
99.17
­0.04
N6 (r=.34550)
93.30
93.26
­0.04
N11 (r=.36825)
87.76
87.73
­0.04
N9 (r=.39100)
82.56
82.53
­0.04
Identification
flux
Référence
Aster %
différence
maille M1 noeud N1
11577.49
11533.37
­0.38
maille M2 noeud N9
8822.98
8855.97
­0.37

4.2 Remarques
La condition aux limites de type rayonnement est fournie sous la forme d'une fonction de la
température interpolée linéairement entre chaque point (on a discrétisé ici la courbe à l'aide de 101
points).
background image
Code_Aster
®
Version
6.4
Titre :
TPNA01 - Problème axisymétrique stationnaire avec rayonnement
Date
:
02/06/03
Auteur(s) :
J.P. LEFEBVRE
Clé
:
V4.41.001-C
Page :
6/8
Manuel de Validation
Fascicule V4.41 : Thermique stationnaire avec rayonnement
HT-66/03/008/A
5 Modélisation
B
5.1
Caractéristiques de la modélisation
Modélisation 3D :


5.2
Caractéristiques du maillage
32 HEXA20

5.3
Fonctionnalités testées
Commande Mot-clé
facteur
Mot-clé
simple
Argument
DEFI_MATERIAU THER_NL
LAMBDA
BETA
THER_NON_LINE TEMP_INIT
STATIONNAIRE
'OUI'
CONVERGENCE
RESI_GLOB_RELA
1.E-2
ITER_GLOB_MAXI
10
CRIT_LAGR_RELA
1.E-3
OPTION
'FLUX_ELNO_TEMP'
background image
Code_Aster
®
Version
6.4
Titre :
TPNA01 - Problème axisymétrique stationnaire avec rayonnement
Date
:
02/06/03
Auteur(s) :
J.P. LEFEBVRE
Clé
:
V4.41.001-C
Page :
7/8
Manuel de Validation
Fascicule V4.41 : Thermique stationnaire avec rayonnement
HT-66/03/008/A
6
Résultats de la modélisation B
6.1 Valeurs
testées
Les noeuds observés ont pour coordonnées : y = z = 0.0

Identification
température
Référence
Aster
% différence
NO106 (x=.30000)
105.55
105.51
­0.03
NO105 (x=.32275)
99.21
99.17
­0.04
NO115 (x=.34550)
93.30
93.26
­0.04
NO125 (x=.36825)
87.76
87.73
­0.04
NO123 (x=.39100)
82.56
82.52
­0.04

Identification
flux
Référence
Aster
% différence
maille MA17 noeud NO106
11577.49
11680.50
+0.89
maille MA16 noeud NO123
8822.98
8968.85
+1.65

6.2 Remarques
La fonction de flux utilisée dans la modélisation A est aussi utilisée ici.
background image
Code_Aster
®
Version
6.4
Titre :
TPNA01 - Problème axisymétrique stationnaire avec rayonnement
Date
:
02/06/03
Auteur(s) :
J.P. LEFEBVRE
Clé
:
V4.41.001-C
Page :
8/8
Manuel de Validation
Fascicule V4.41 : Thermique stationnaire avec rayonnement
HT-66/03/008/A
7
Synthèses des résultats
La prise en compte des conditions de rayonnement est tout à fait correcte dans ce cas stationnaire.
Notons que ce test fait intervenir un coefficient de conductivité thermique constant, la seule
non-linéarité porte donc sur les conditions aux limites.
Les erreurs sont plus élevées sur le calcul du flux, ce que l'on peut expliquer par l'utilisation des
lissages aux noeuds, effectués à partir des valeurs calculées aux points d'intégration (points de
GAUSS).