Code_Aster ®
Version
4.0
Titrate:
Thermal TTLV100 Choc in a pipe with condition of exchange
Date: 01/12/98
Author (S)
:
J.M. PROIX, Mr. A. REDON
Key:
V4.25.100-A Page:
1/8
Organization (S): EDF/IMA/MN
Handbook of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
Document: V4.25.100
TTLV100 - Thermal Choc in a pipe with
condition of exchange
Summary:
This test of transitory linear thermics consists in imposing a cold thermal shock on a hollow roll
presumedly infinite using a limiting condition of exchange.
The shock is modelled by a linear slope T = ­ 269°C into 12 S.
The problem is dealt with into axisymmetric.
The reference solution is obtained on a fine network.
The test is carried out on 2 modelings: (TRIA3, QUAD4) and (TRIA6, QUAD9).
One tests the algorithm of linear thermics transitory when the matrix of mass is diagonalisée
(modeling AXIS_DIAG with `'farmhouse lumping '').
Handbook of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A

Code_Aster ®
Version
4.0
Titrate:
Thermal TTLV100 Choc in a pipe with condition of exchange
Date: 01/12/98
Author (S)
:
J.M. PROIX, Mr. A. REDON
Key:
V4.25.100-A Page:
2/8
1
Problem of reference
1.1 Geometry
Z
IH
With
B
R
Re
M1 m2
Z
IH = 417 mm
Re = 496 mm
R (A) = 417 mm
R (B) = 496 mm
R (M1) = 443.43 mm
R (m2) = 469.67 mm
1.2
Material properties
= 19.97 W/m °C
CP = 4.89488 106 J/m3 °C
1.3
Boundary conditions and loadings
T
Exchange

= (
-

H T
T (R, T))
N
ext.
R =ri
2
with
H =
.
40000 W/m ° C
T

()
With
0 =
.
289 ° C
ext.
T =
T () A 12 =.
20 ° C
ext.
T =
S
1.4 Conditions
initial
T (R, 0) = 289°C for any R
Discretization in time (T):
12
not for
[0.
,
12.]
that is to say
T = 1. S
2
not for
[12.
,
20.]
that is to say
T = 4. S
4
not for
[20.
,
100.]
that is to say
T = 20. S
2
not for
[100.
,
200.]
that is to say
T = 50. S
2
not for
[200.
,
400.]
that is to say
T = 100. S
8
not for
[400.
,
2000.]
that is to say
T = 200. S
Handbook of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A

Code_Aster ®
Version
4.0
Titrate:
Thermal TTLV100 Choc in a pipe with condition of exchange
Date: 01/12/98
Author (S)
:
J.M. PROIX, Mr. A. REDON
Key:
V4.25.100-A Page:
3/8
2
Reference solution
2.1
Method of calculation used for the reference solution
The reference solution is obtained on a fine network comprising 99 quadratic elements QUAD8
in the thickness without option of matrix of diagonal thermal mass.
2.2
Results of reference
Temperatures at the points M1 (R = 443.33) and m2 (R = 469.67),
and at various moments (T = 12., 100., 600. and 2000.).
Handbook of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A

Code_Aster ®
Version
4.0
Titrate:
Thermal TTLV100 Choc in a pipe with condition of exchange
Date: 01/12/98
Author (S)
:
J.M. PROIX, Mr. A. REDON
Key:
V4.25.100-A Page:
4/8
3 Modeling
With
3.1
Characteristics of modeling
TRIA3, QUAD4
The infinite cylinder being supposed, one nets only one section height H = 40 mm with 2 layers
elements.
M'1
M'2
D
C
Limiting conditions

on [AB], [BC] and [CD]: null flow

H
on [AD]: exchange H, Text
Text
With
B
M1
M2
289 °C
20°C
T
0
12 S
Initial conditions
points
nodes
R
Z
T = 289 °C

M1 N10
443.33 0.0
M2 N5
469.67
0.0
M'1
N11
443.33 40.
M'2
N7
469.67
40.
3.2
Characteristics of the grid
A number of nodes: 12
A number of meshs and types: 3 QUAD4, 6 TRIA6
3.3 Functionalities
tested
Commands
Keys
AFFE_MODELE
MODELISATION
AXIS_DIAG
[U4.22.01]
AFFE_CHAR_THER_F
ECHANGE
[U4.25.02]
THER_LINEAIRE
[U4.62.01]
Handbook of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A

Code_Aster ®
Version
4.0
Titrate:
Thermal TTLV100 Choc in a pipe with condition of exchange
Date: 01/12/98
Author (S)
:
J.M. PROIX, Mr. A. REDON
Key:
V4.25.100-A Page:
5/8
4
Results of modeling A
4.1 Values
tested
Identification
Reference
Aster
% difference
M1 (R = 443.33)
T = 12.
288.64
282.63
­ 2.08
T = 20.
202.76
199.84
­ 1.44
T = 600.
93.027
92.821
­ 0.22
T = 2000.
29.419
29.865
+1.51
M2 (R = 469.67)
T = 12.
289.00
288.84
­ 0.06
T = 20.
275.04
268.63
­ 2.32
T = 600.
143.00
142.74
­ 0.18
T = 2000.
35.858
36.629
+2.15
M'1 (R = 443.33)
T = 12.
288.64
283.09
­ 1.92
T = 20.
202.76
206.66
+1.92
T = 600.
93.027
93.731
+0.76
T = 2000.
29.419
29.988
+1.93
M'2 (R = 469.67)
T = 12.
289.00
288.82
­ 0.06
T = 20.
275.04
267.66
­ 2.68
T = 600.
143.00
141.57
­ 1.00
T = 2000.
35.858
36.470
+1.71
4.2 Parameters
of execution
Version: 4.01.12
Machine: CRAY C98
System: UNICOS
Obstruction memory:
8 megawords
Time CPU To use: 7.6 seconds
Handbook of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A

Code_Aster ®
Version
4.0
Titrate:
Thermal TTLV100 Choc in a pipe with condition of exchange
Date: 01/12/98
Author (S)
:
J.M. PROIX, Mr. A. REDON
Key:
V4.25.100-A Page:
6/8
5 Modeling
B
5.1
Characteristics of modeling
TRIA6, QUAD9
The infinite cylinder being supposed, one nets only one section height H = 40 mm with 2 layers
elements.
M'1
M'2
D
C
Limiting conditions

on [AB], [BC] and [CD]: null flow

H
on [AD]: exchange H, Text
Text
With
B
M1
M
289 °C
2
20°C
T
0
12 S
Initial conditions
points
nodes
R
Z
T = 289 °C

M1 N25
443.33 0.0
M2 N9
469.67
0.0
M'1
N28
443.33 40.
M'2
N1
469.67
40.
5.2
Characteristics of the grid
A number of nodes: 35
A number of meshs and types: 3 QUAD9, 6 TRIA6
5.3 Functionalities
tested
Commands
Keys
CREA_MAILLAGE
MODI_MAILLE
OPTION:
“QUAD8_9”
[U4.12.06]
AFFE_MODELE
MODELISATION
AXIS_DIAG
[U4.22.01]
AFFE_CHAR_THER_F
ECHANGE
[U4.25.02]
THER_LINEAIRE
[U4.62.01]
Handbook of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A

Code_Aster ®
Version
4.0
Titrate:
Thermal TTLV100 Choc in a pipe with condition of exchange
Date: 01/12/98
Author (S)
:
J.M. PROIX, Mr. A. REDON
Key:
V4.25.100-A Page:
7/8
6
Results of modeling B
6.1 Values
tested
Identification
Reference
Aster
% difference
M1 (R = 443.33)
T = 12.
288.64
286.80
­ 0.63
T = 20.
202.76
202.25
­ 0.25
T = 600.
93.027
92.955
­ 0.08
T = 2000.
29.419
29.524
+0.36
M2 (R = 469.67)
T = 12.
289.00
288.99
+0.00
T = 20.
275.04
273.35
­ 0.61
T = 600.
143.00
142.99
­ 0.00
T = 2000.
35.858
36.050
+0.54
M'1 (R = 443.33)
T = 12.
288.64
287.13
­ 0.52
T = 20.
202.76
205.35
+1.28
T = 600.
93.027
93.378
+0.38
T = 2000.
29.419
29.580
+0.55
M'2 (R = 469.67)
T = 12.
289.00
288.99
+0.00
T = 20.
275.04
272.65
­ 0.87
T = 600.
143.00
142.39
­ 0.43
T = 2000.
35.858
35.972
+0.32
6.2 Parameters
of execution
Version: 4.01.12
Machine: CRAY C98
System: UNICOS
Obstruction memory:
8 megawords
Time CPU To use: 8.0 seconds
Handbook of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A

Code_Aster ®
Version
4.0
Titrate:
Thermal TTLV100 Choc in a pipe with condition of exchange
Date: 01/12/98
Author (S)
:
J.M. PROIX, Mr. A. REDON
Key:
V4.25.100-A Page:
8/8
7
Summary of the results
Modeling “AXIS_DIAG” gives rather satisfactory results. Although grid
comprise that 3 elements in the thickness, the variation on the temperatures remains lower than 2.7%.
In spite of the violence of the thermal shock, the diagonalisation of the matrix of mass makes it possible to obtain one
solution in temperature which does not oscillate during the transient.
Handbook of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A