Code_Aster ®
Version
5.0
Titrate:
TPNL300 unidimensional Transfert of heat with radiation
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.42.300-A Page:
1/8
Organization (S): EDF/AMA, EDF/UTO/LOCATED, Delta CAD
Handbook of Validation
V4.42 booklet: Nonlinear stationary thermics of the linear structures
V4.42.300 document
TPNL300 - Unidimensional Transfert of heat
with radiation
Summary:
This test results from the validation independent of version 3 in nonlinear stationary thermics.
It is about a linear problem 1D represented by two modelings, one planes, the other voluminal one.
The functionalities tested are as follows:
· plane thermal element,
· voluminal thermal element,
· limiting conditions: (imposed temperatures, radiation).
The interest of the test lies in the taking into account of the radiation.
The results are compared with those provided by NAFEMS.
Handbook of Validation
V4.42 booklet: Nonlinear stationary thermics of the linear structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPNL300 unidimensional Transfert of heat with radiation
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.42.300-A Page:
2/8
1
Problem of reference
1.1 Geometry
With
B
X
0.1 m
1.2
Properties of material
= 55.6 W/m.°C
Thermal conductivity
C
= 460 J/kg.°C
Specific heat
= 7850 kg/m3 Masse
voluminal
1.3
Boundary conditions and loadings
· temperature imposed on point a: MT = 726.85°C,
· exchange by radiation at point b:
- outside temperature = 26.85°C,
- = 0.98 emissivity,
- = 5.67x10-8 W/m ² K4 (constant of Stefan-Boltzmann).
1.4 Conditions
initial
Without object.
Handbook of Validation
V4.42 booklet: Nonlinear stationary thermics of the linear structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPNL300 unidimensional Transfert of heat with radiation
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.42.300-A Page:
3/8
2
Reference solution
2.1
Method of calculation used for the reference solution
The reference solution is that given in card “TEST n°2” of the tests of reference published by
NAFEMS.
2.2
Results of reference
Temperature at point b: T = 653.85°C
2.3
Uncertainty on the solution
Nonavailable on card NAFEMS.
2.4 References
bibliographical
[1]
NAFEMS (the National Agency for Finite Element Methods and Standard (the U.K.)): “The
standard NAFEMS Benchmarcks ", TNSB rév 3, October 1990.
Handbook of Validation
V4.42 booklet: Nonlinear stationary thermics of the linear structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPNL300 unidimensional Transfert of heat with radiation
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.42.300-A Page:
4/8
3 Modeling
With
3.1
Characteristics of modeling
PLAN (TRIA3, QUAD4)
y
C
D
N55
0.02m
N53
N51
With
B
X
Boundary conditions:
- Dimensioned AD:
T = 726.85°C
- Dimensioned AB, CD: = 0
- Dimensioned BC:
Text=26.85°C
= 0.98
3.2
Characteristics of the grid
A number of nodes:
55
A number of meshs and types: 60: (20 QUAD4, 40 TRIA3)
3.3 Functionalities
tested
Commands
AFFE_MODELE
THERMIQUE
PLAN
TOUT
DEFI_MATERIAU
THER_NL
AFFE_CHAR_THER_F
FLUX_NL
TEMP_IMPO
THER_NON_LINE
TEMP_INIT
STATIONNAIRE
“OUI”
EXCIT
CHARGE
RECU_CHAMP
INST
0.
Handbook of Validation
V4.42 booklet: Nonlinear stationary thermics of the linear structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPNL300 unidimensional Transfert of heat with radiation
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.42.300-A Page:
5/8
4
Results of modeling A
4.1 Values
tested
Identification Reference
Aster %
difference
tolerance
Temperature at the point B
in °C
N51 653.85
653.87
0.003
2%
N53 653.85
653.87
0.003
2%
N55 653.85
653.87
0.003
2%
4.2 Parameters
of execution
Version: 5.03
Machine: SGI - ORIGIN 2000 - R12000
Obstruction memory:
8 megawords
Time CPU To use: 2.64 seconds
Handbook of Validation
V4.42 booklet: Nonlinear stationary thermics of the linear structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPNL300 unidimensional Transfert of heat with radiation
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.42.300-A Page:
6/8
5 Modeling
B
5.1
Characteristics of modeling
3D (HEXA8)
Z
0.02m
y
N91
0.02m
N97
Boundary conditions:
X
N92
- face X = 0
T=726.85°C
N98
- face X = 0.1
Text=26.85°C
= 0.98
- others faces
= 0
5.2
Characteristics of the grid
A number of nodes:
99
A number of meshs and types: 40 HEXA8
5.3 Functionalities
tested
Commands
AFFE_MODELE
THERMIQUE
PLAN
TOUT
DEFI_MATERIAU
THER_NL
AFFE_CHAR_THER_F
FLUX_NL
TEMP_IMPO
THER_NON_LINE
TEMP_INIT
STATIONNAIRE
“OUI”
EXCIT
CHARGE
RECU_CHAMP
INST
0.
Handbook of Validation
V4.42 booklet: Nonlinear stationary thermics of the linear structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPNL300 unidimensional Transfert of heat with radiation
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.42.300-A Page:
7/8
6
Results of modeling B
6.1 Values
tested
Identification Reference
Aster %
difference
tolerance
Temperature at the point B
in °C
N91 653.85
653.87
0.003
2%
N92 653.85
653.87
0.003
2%
N97 653.85
653.87
0.003
2%
N98 653.85
653.87
0.003
2%
6.2 Parameters
of execution
Version: 5.03
Machine: SGI - ORIGIN 2000 - R12000
Obstruction memory:
8 megawords
Time CPU To use: 2.66 seconds
Handbook of Validation
V4.42 booklet: Nonlinear stationary thermics of the linear structures
HT-66/02/001/A
Code_Aster ®
Version
5.0
Titrate:
TPNL300 unidimensional Transfert of heat with radiation
Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.42.300-A Page:
8/8
7
Summary of the results
This test is recommended by NAFEMS (but with another type of grid).
Two modelings give very satisfactory results, the maximum change obtained is of
0.003%.
For modeling PLAN, in spite of the not-symmetry of the grid, one notes that the temperature with
nodes (points of observation) pertaining to the TRIA3 and the QUAD4 is identical.
The limiting condition of radiation was imposed via a nonlinear loading of flow (flow function
temperature). In this test the taking into account of the radiation is completely correct.
This test with license to test command AFFE_CHAR_THER_F (associate with the operand FLUX_NL which
allows to affect a flow non_linéaire in the case of) modelings PLAN and 3D.
Handbook of Validation
V4.42 booklet: Nonlinear stationary thermics of the linear structures
HT-66/02/001/A
Outline document