Code_Aster ®
Version
5.0
Titrate:
TTLP301 Transfert of heat in a perforated plate

Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.301-A Page:
1/8

Organization (S): EDF/AMA, EDF/UTO/LOCATED, Delta CAD
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
V4.23.301 document

TTLP301 - Transfert of heat in a plate
perforated

Summary:

This test, industrialist, result from the validation independent of version 3 in linear transitory thermics.

It is about a problem plane 2D represented by two modelings, one planes, the other voluminal one.

The functionalities tested are as follows:
· plane thermal element,
· voluminal thermal element,
· transitory algorithm of thermics,
· boundary conditions of exchange and flow.


Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A

Code_Aster ®
Version
5.0
Titrate:
TTLP301 Transfert of heat in a perforated plate

Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.301-A Page:
2/8

1
Problem of reference

1.1 Geometry
Ø 12.7
Coolant
D
E
Ø 15.87
C
Ø 18.79
Graphite
With
B
Fuel


1.2
Properties of material


=
0.1 W/cm °C
Thermal conductivity
CP = 1.0 J/Cm3 °C
Voluminal heat

1.3
Boundary conditions and loadings

· [ED] Coefficient of convection H = 1 W/cm2 °C Text = 0°C,
· [AC] Density of flow Q = 1 W/cm2,
· [AB], [BE], [CD.] = 0.

1.4 Conditions
initial

T (T = 0) = 0

Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A

Code_Aster ®
Version
5.0
Titrate:
TTLP301 Transfert of heat in a perforated plate

Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.301-A Page:
3/8

2
Reference solution

2.1
Method of calculation used for the reference solution

The reference solution is a numerical solution obtained by the finite element method. This
solution is based on a linear triangular network presented below. Calculations were
carried out by considering an increment of time T = 0.01 S.

2.2
Results of reference

Temperature at the point C for t=0.1, 0.2,…, 0.9, 1.0, 1.1, 1.2 S

2.3 References
bibliographical

[1]
J. Donea, “One the accuracy off finite element solutions to the transient heat-conduction
equation ", Int. J. num. Meth. Engng, flight 8, n°1, p 103-110, 1974

Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A

Code_Aster ®
Version
5.0
Titrate:
TTLP301 Transfert of heat in a perforated plate

Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.301-A Page:
4/8

3 Modeling
With

3.1
Characteristics of modeling

PLAN (TRIA6)

Limiting conditions:
- dimensioned AC
Q
= 1 W/cm ²
- dimensioned ED
H
= 1 W/cm ² °C
Text = 0°C
- dimensioned AB, BE, cd.
= 0
y
D
E
N1
C
With
B
X


3.2
Characteristics of the grid

A number of nodes:
718
A number of meshs and types: TRIA6: 335 (SEG3: 22)

3.3
Functionalities tested

Commands




AFFE_MODELE
THERMIQUE
PLAN
TOUT
AFFE_CHAR_THER
FLUX_REP

ECHANGE

THER_LINEAIRE
TEMP_INIT
VALE

INCREMENT
LIST_INST

RECU_CHAMP
INST


3.4 Remarks

Limiting condition, Flux = 0., not modelled (implicit).

The discretization in step of time is as follows:

10 steps for [0. , 0.2D+0] is T = 2.D-2
10 steps for [0.2D+0, 1.2D+1] are T = 1.D-1
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A

Code_Aster ®
Version
5.0
Titrate:
TTLP301 Transfert of heat in a perforated plate

Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.301-A Page:
5/8

4
Results of modeling A

4.1 Values
tested

Identification Reference
Aster
% difference
Tolerance
Node Temps (S)
T (°C)
T (°C)


n1 0.1
1.045
1.0664
2.05%
2%
“0.2
1.447
1.4515
0.31%
2%
“0.3
1.742
1.7480
0.35%
2%
“0.4
1.982
1.9847
0.14%
2%
“0.5
2.189
2.1929
0.18%
2%
“0.6
2.373
2.3757
0.11%
2%
“0.7
2.541
2.5451
0.16%
2%
“0.8
2.698
2.7010
0.11%
2%
“0.9
2.846
2.8491
0.11%
2%
“1.0
2.986
2.9889
0.10%
2%
“1.1
3.120
3.1232
0.10%
2%
“1.2
3.248
3.2517
0.11%
2%

4.2 Parameters
of execution

Version: 5.03


Machine: SGI - ORIGIN 2000 - R12000

Obstruction memory:
8 megawords
Time CPU To use: 3.89 seconds
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A

Code_Aster ®
Version
5.0
Titrate:
TTLP301 Transfert of heat in a perforated plate

Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.301-A Page:
6/8

5 Modeling
B

5.1
Characteristics of modeling

3D (PENTA15)

Limiting conditions:
- dimensioned AC
Q
= 1 W/cm ²
- dimensioned ED
H
= 1 W/cm ² °C
Text = 0°C
- dimensioned AB, BE, cd. = 0
y
D
n1
E
C
n1
n719
n911
n1629
n1821
With
B
X


5.2
Characteristics of the grid

A number of nodes:
2538
A number of meshs and types: PENTA15: 670 (QUAD8: 44)

5.3 Functionalities
tested


Commands




AFFE_MODELE
THERMIQUE
3D
TOUT
AFFE_CHAR_THER
FLUX_REP

ECHANGE

THER_LINEAIRE
TEMP_INIT
VALE

INCREMENT
LIST_INST

RECU_CHAMP
INST


5.4 Remarks

Condition limits = 0. implicit: not modelled.

The discretization in step of time is as follows:

10 steps for [0. , 0.2D+0] is T = 2.D-2
10 steps for [0.2D+0, 1.2D+1] are T = 1.D-1
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A

Code_Aster ®
Version
5.0
Titrate:
TTLP301 Transfert of heat in a perforated plate

Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.301-A Page:
7/8

6
Results of modeling B

6.1 Values
tested

Identification Reference
Aster
% difference
Tolerance
Node Temps (S)
T (°C)
T (°C)


n1 0.1
1.045
1.0665
2.05%
2%
“0.2
1.447
1.4514
0.30%
2%
“0.3
1.742
1.7480
0.35%
2%
“0.4
1.982
1.9847
0.14%
2%
“0.5
2.189
2.1929
0.18%
2%
“0.6
2.373
2.3757
0.11%
2%
“0.7
2.541
2.5451
0.16%
2%
“0.8
2.698
2.7010
0.11%
2%
“0.9
2.846
2.8491
0.11%
2%
“1.0
2.986
2.9889
0.10%
2%
“1.1
3.120
3.1232
0.10%
2%
“1.2
3.248
3.2517
0.11%
2%

6.2 Remarks

Difference between the values with the nodes n1, n719, n911, n1629, n1821 about 1.e-8.

6.3 Parameters
of execution

Version: 5.03


Machine: SGI - ORIGIN 2000 - R12000

Obstruction memory:
8 megawords
Time CPU To use: 7.74 seconds
Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A

Code_Aster ®
Version
5.0
Titrate:
TTLP301 Transfert of heat in a perforated plate

Date:
20/09/02
Author (S):
C. DURAND, E. SCREW, F. LEBOUVIER Clé
:
V4.23.301-A Page:
8/8

7
Summary of the results

Two modelings give results whose only one value exceeds the fixed tolerance
initially. The maximum change is equal to 2.05%, and is thus not very higher than the fixed tolerance (2%). it
is located on the smallest value of temperature and for the weakest moment the T (starting of
problem).

Two modelings, PLAN (TRIA6) and 3D (PENTA15) gives the same results, which is
normal since the grid and the degree of interpolation are identical.

A finer grid in the zone of the N1 node should improve quality of the results which are
regarded as acceptable taking into account modelings carried out.

Handbook of Validation
V4.23 booklet: Transitory thermics of the plane systems
HT-66/02/001/A

Outline document