Code_Aster
®
Version
4.0
Titrate:
Thermal TTLV100 Shock 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
Manual of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A
Organization (S)
: EDF/IMA/MN
Manual of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
Document: V4.25.100
TTLV100 - Thermal shock 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 modelized by a linear slope
T = 269°C in 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 '').
Code_Aster
®
Version
4.0
Titrate:
Thermal TTLV100 Shock 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
Manual of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A
1
Problem of reference
1.1 Geometry
Z
Re
IH
IH = 417 mm
Re = 496 mm
R
With
B
M2
M1
Z
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
C p
= 4.89488 106 J/m
3
°C
1.3
Boundary conditions and loadings
()
()
Exchange
with
T
N
H T
T R T
H
W m
C
T
With
C
T
With
C
R R
ext.
ext.
T
ext.
T
S
I
=
=
=
=
-
=
°
=
°
=
°
(
(,))
.
/
.
.
40000
289
20
2
0
12
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
Code_Aster
®
Version
4.0
Titrate:
Thermal TTLV100 Shock 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
Manual of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A
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.).
Code_Aster
®
Version
4.0
Titrate:
Thermal TTLV100 Shock 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
Manual of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A
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.
With
B
D
C
H
M
1
Me
2
Me
1
M
2
Limiting conditions
on [AB], [BC] and [CD]:
null flow
on [AD]: exchange H, T
ext.
Initial conditions
T = 289 °C
points
M1
M2
nodes
N10
N5
M'2
M'1
N7
N11
Z
0.0
0.0
40.
40.
R
443.33
469.67
443.33
469.67
T
0
12 S
289 °C
T
ext.
20°C
3.2
Characteristics of the mesh
A number of nodes: 12
A number of meshs and types: 3
QUAD4
, 6
TRIA6
3.3 Functionalities
tested
Controls
Keys
AFFE_MODELE
MODELING
AXIS_DIAG
[U4.22.01]
AFFE_CHAR_THER_F
EXCHANGE
[U4.25.02]
THER_LINEAIRE
[U4.62.01]
Code_Aster
®
Version
4.0
Titrate:
Thermal TTLV100 Shock 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
Manual of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A
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
Overall dimension memory:
8 megawords
Time CPU To use: 7.6 seconds
Code_Aster
®
Version
4.0
Titrate:
Thermal TTLV100 Shock 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
Manual of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A
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.
With
B
D
C
H
M
1
Me
2
Me
1
M
2
Limiting conditions
on [AB], [BC] and [CD]:
null flow
on [AD]: exchange H, T
ext.
Initial conditions
T = 289 °C
points
M1
M2
nodes
N25
N9
M'2
M'1
N1
N28
Z
0.0
0.0
40.
40.
R
443.33
469.67
443.33
469.67
T
0
12 S
289 °C
T
ext.
20°C
5.2
Characteristics of the mesh
A number of nodes: 35
A number of meshs and types: 3
QUAD9
, 6
TRIA6
5.3 Functionalities
tested
Controls
Keys
CREA_MAILLAGE
MODI_MAILLE
OPTION:
“QUAD8_9”
[U4.12.06]
AFFE_MODELE
MODELING
AXIS_DIAG
[U4.22.01]
AFFE_CHAR_THER_F
EXCHANGE
[U4.25.02]
THER_LINEAIRE
[U4.62.01]
Code_Aster
®
Version
4.0
Titrate:
Thermal TTLV100 Shock 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
Manual of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A
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
Overall dimension memory:
8 megawords
Time CPU To use: 8.0 seconds
Code_Aster
®
Version
4.0
Titrate:
Thermal TTLV100 Shock 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
Manual of Validation
V4.25 booklet: Transitory thermics of the voluminal structures
HI-75/98/040 - Ind A
7
Summary of the results
Modeling
“AXIS_DIAG”
give rather satisfactory results. Although mesh
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.