Version

Code_Aster default
Titre : TPLV06 - Dégagement de puissance dans une sphère c[...] Date : 02/02/2011 Page : 1/5
Responsable : HAELEWYN Jessica Clé : V4.04.006 Révision :
5a5d2302fdd3

TPLV06 - Release of power in a hollow sphere

Summary:

This test is resulting from the validation independent of version 3 in linear stationary thermics.

It is about a three-dimensional problem which aims to validate the voluminal thermal element subjected to an
imposed temperature and a source of heat.

This case test understands a modeling 3D. The results are compared with an analytical solution (VPCS).

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2021 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : TPLV06 - Dégagement de puissance dans une sphère c[...] Date : 02/02/2011 Page : 2/5
Responsable : HAELEWYN Jessica Clé : V4.04.006 Révision :
5a5d2302fdd3

1 Problem of reference
1.1 Geometry
 
Re 
Ri  Sphère creuse: 
­ Rayon intérieur   Ri =  1 m 
­ Rayon extérieur   Re = 2 m 

O  r 

Q 

1.2 Properties of material

=1. W /m ° C Thermal conductivity

1.3 Boundary conditions and loadings

• T i=T  r=R i=20 ° C ,
• T e =T r =Re =20 ° C ,
• Q=100 W / m 3 .

1.4 Initial conditions

Without object.

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2021 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : TPLV06 - Dégagement de puissance dans une sphère c[...] Date : 02/02/2011 Page : 3/5
Responsable : HAELEWYN Jessica Clé : V4.04.006 Révision :
5a5d2302fdd3

2 Reference solution
2.1 Method of calculating used for the reference solution
The reference solution is that given in card TPLV06/89 of guide VPCS.

● Temperature according to r :

[[ [ ]
]
2 2 1 1
 Re − Ri  −
Q Ri r 2 2
T =T i −r −Ri 
6 1 1
−
R i Re ]
● Density flux according to r :

=−4  r 
2 dT
dr
=−
2Q
3
2 2
 Re −Ri 
[
1 1
−
R i Re
−2 r
3
[ ] ]
2.2 Results of reference
Temperature in r =1.25 ; 1.5 and 1.75 m

2.3 Uncertainty on the solution

Analytical solution.

2.4 Bibliographical references

[1] Guide of validation of the software packages of structural analysis. French company of the
Mechanics, AFNOR 1990 ISBN 2-12-486611-7

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2021 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : TPLV06 - Dégagement de puissance dans une sphère c[...] Date : 02/02/2011 Page : 4/5
Responsable : HAELEWYN Jessica Clé : V4.04.006 Révision :
5a5d2302fdd3

3 Modeling A
3.1 Characteristics of modeling
3D (HEXA8)

  G 
N101 
 
 Conditions limites:  C 
  N106  N 105 
­ faces ABCD, E FGH,  
           ABFE, DCGH   =     0 
N111  N110 
­ faces ADH E, BCGF    T =   20°C 
F  N1 
N116  N115 

N6  x 
N120 
N121  N11  B 
H  N 5 
z  N16  N125 
D  N10 
N15 
E  N21 
N20 
C  N25   
A 
30°  30° 
y 

3.2 Characteristics of the grid

Many nodes: 125
Many meshs and types: 64 HEXA8

3.3 Sizes tested and results

Identification Reference Aster % difference tolerance
Temperature ° C

R = 1.25 (N16) 30,625 30,471  0,504 1%

R = 1.25 (N116) 30,625 30,471  0,504 1%
R = 1.25 (N20) 30,625 30,462  0,532 1%
R = 1.25 (N120) 30,625 30,462  0,532 1%
R = 1.50 (N11) 32,500 32,337  0,500 1%
R = 1.50 (N111) 32,500 32,337  0,500 1%
R = 1.50 (N15) 32,500 32,335  0,507 1%
R = 1.50 (N115) 32,500 32,335  0,507 1%
R = 1.75 (N6) 28,482 28,379  0,362 1%
R = 1.75 (N106) 28,482 28,379  0,362 1%
R = 1.75 (N10) 28,482 28,382  0,351 1%
R = 1.75 (N110) 28,482 28,382  0,351 1%

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2021 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)
 Version
Code_Aster default
Titre : TPLV06 - Dégagement de puissance dans une sphère c[...] Date : 02/02/2011 Page : 5/5
Responsable : HAELEWYN Jessica Clé : V4.04.006 Révision :
5a5d2302fdd3

4 Summary of the results

The got results are satisfactory, the maximum change obtained is of 0.53%.

Modeling 3D used to model this sphere is correct.

The quality of the results could be still improved in:

• carrying out a finer grid of the portion of sphere,

• choosing quadratic elements for better approximating the reference solution.

Warning : The translation process used on this website is a "Machine Translation". It may be imprecise and inaccurate in whole or in part and is
provided as a convenience.
Copyright 2021 EDF R&D - Licensed under the terms of the GNU FDL (http://www.gnu.org/copyleft/fdl.html)