background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
1/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Organization (S):
EDF-R & D/SINETICS















Manual of Reference
R4.03 booklet: Analyze sensitivity
Document: R4.03.02



Calculation of sensitivities in thermics


Summary:

During digital simulations obtaining a rough result is not sufficient any more. The user is moreover in
petitioning of calculation of sensitivity compared to the data input of the problem. That allows him
to estimate the uncertainty which the field result according to the law of variation of the data answers. This
derived is also the basic substrate of problems opposite (retiming of parameters…) and of problems
of optimization.
This sensitivity can be obtained “manually”, but the experiment shows that these parametric studies
are often expensive, little mutualisables and less reliable than an analytical calculation established in the software of
calculation.
In this note, one places oneself in the perimeter of use of the standard thermal operators of
Code_Aster and one are interested in this analytical sensitivity of the field of temperature and its flow by
report/ratio with the characteristics material and the loadings. One described the process allowing there to exhume
the linear system which this derivative checks. In order to minimize the overcost calculation, a particular effort was
brought to bind its resolution to that of the initial problem.
One details theoretical, numerical work and data processing which governed the establishment of these
calculations of sensitivity in the code. One specifies their properties and their limitations while connecting these
considerations with a precise parameter setting of the accused operators and with the choices of modeling of the code. One has
tried constantly to bind different the items approached while detailing, has minimum, the demonstrations a little
techniques.
Required environment, the parameter setting and the perimeter of use of this new functionality
are described. An example extracted from an official case-test is clarified.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
2/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Contents
1
Problems ....................................................................................................................................... 4
2
Linear thermics ................................................................................................................................. 8
2.1
Derived compared to voluminal heat .................................................................................. 10
2.1.1
Theoretical elements ............................................................................................................ 10
2.1.2
Establishment in Code_Aster ......................................................................................... 13
2.2
Derived compared to thermal conductivity ........................................................................... 15
2.2.1
Theoretical elements ............................................................................................................ 15
2.2.2
Establishment in Code_Aster ......................................................................................... 16
2.3
Derived compared to the source ..................................................................................................... 18
2.3.1
Theoretical elements ............................................................................................................ 18
2.3.2
Establishment in Code_Aster ......................................................................................... 18
2.4
Derived compared to the imposed temperature ............................................................................. 19
2.4.1
Theoretical elements ............................................................................................................ 19
2.4.2
Establishment in Code_Aster ......................................................................................... 20
2.5
Derived compared to normal flow imposed ................................................................................... 22
2.5.1
Theoretical elements ............................................................................................................ 22
2.5.2
Establishment in Code_Aster ......................................................................................... 22
2.6
Derived compared to the coefficient from convectif exchange ................................................................ 24
2.6.1
Theoretical elements ............................................................................................................ 24
2.6.2
Establishment in Code_Aster ......................................................................................... 25
2.7
Derived compared to the outside air temperature ........................................................................... 26
2.7.1
Theoretical elements ............................................................................................................ 26
2.7.2
Establishment in Code_Aster ......................................................................................... 27
3
Nonlinear thermics ........................................................................................................................ 28
3.1
Derived compared to voluminal heat .................................................................................. 29
3.1.1
Theoretical elements ............................................................................................................ 29
3.1.2
Establishment in Code_Aster ......................................................................................... 34
3.2
Derived compared to thermal conductivity ............................................................................. 37
3.2.1
Theoretical elements ............................................................................................................ 37
3.2.2
Establishment in Code_Aster ......................................................................................... 38
3.3
Derived compared to the source ..................................................................................................... 38
3.3.1
Theoretical elements ............................................................................................................ 38
3.3.2
Establishment in Code_Aster ......................................................................................... 39
3.4
Derived compared to the imposed temperature ............................................................................. 39
3.4.1
Theoretical elements ............................................................................................................ 39
3.4.2
Establishment in Code_Aster ......................................................................................... 40
3.5
Derived compared to linear normal flow imposed ...................................................................... 41
3.5.1
Theoretical elements ............................................................................................................ 41
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
3/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
3.5.2
Establishment in Code_Aster ......................................................................................... 42
3.6
Derived compared to non-linear normal flow imposed ............................................................... 42
3.6.1
Theoretical elements ............................................................................................................ 42
3.6.2
Establishment in Code_Aster ......................................................................................... 43
3.7
Derived compared to the coefficient from convectif exchange ................................................................ 43
3.7.1
Theoretical elements ............................................................................................................ 43
3.7.2
Establishment in Code_Aster ......................................................................................... 44
3.8
Derived compared to the outside air temperature ........................................................................... 45
3.8.1
Theoretical elements ............................................................................................................ 45
3.8.2
Establishment in Code_Aster ......................................................................................... 46
3.9
Derived compared to emissivity/constant of Stefan-Boltzmann .............................................. 46
3.9.1
Theoretical elements ............................................................................................................ 46
3.9.2
Establishment in Code_Aster ......................................................................................... 47
3.10
Derived compared to the temperature ad infinitum ...................................................................... 48
3.10.1
Theoretical elements ................................................................................................ 48
3.10.2
Establishment in Code_Aster ............................................................................ 49
4
Summary of the sensitivities of the temperature .................................................................................. 50
5
Sensitivity of the heat transfer rate ............................................................................................................... 54
6
Implementation in Code_Aster .................................................................................................... 56
6.1
Particular difficulties .................................................................................................................. 56
6.2
Environments necessary/parameter settings ................................................................................. 56
6.3
Perimeter of use .................................................................................................................... 60
6.4
Example of use ...................................................................................................................... 61
7
Perspective conclusion/....................................................................................................................... 63
8
Bibliography ........................................................................................................................................ 64
Appendix 1
Concept of derived “within the meaning of the distributions” ...................................................... 65
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
4/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
1 Problems
During digital simulation, obtaining a result starting from a data file is not sufficient any more.
Taking into account uncertainties which weigh on the evaluation of the loadings, the geometries and of
characteristics material, taking into account also the numerical approximations due to modelings
employees, with their space-time discretizations and the algorithms of resolution, the user is
increasingly petitioning of calculations of sensitivity [bib5]. One then seeks to evaluate
sensitivity of a variable compared to an input datum of the problem. It makes it possible to estimate
the uncertainty (the probabilistic function) to which answers the field result according to the law of
variation of certain data.
This derivative can be estimated “manually”, but these of parametric studies are often
expensive, little mutualisables and less reliable than an analytical calculation established in the code.
Note:
·
The sensitivities by finished differences are of course dependant on the parameters on
shift and of the mesh, but into non-linear, another worsening factor is superimposed:
degree of convergence of the solution. In any rigor, the aforementioned also intervenes on quality
analytical sensitivities, because one uses the field of temperature solution for
to assemble the linear system “derived”.
·
The third channel gathers the techniques of automatic differentiation (ODYSSEY [bib7],
[bib9]…) but they are not plantable in Code_Aster because of sound
software architecture (transmission of arguments between the routines by pointer…). Of
any way these products are still “relatively embryonic” and their use
seem fixed quotas for with model problems or parts of software well
specific. The ideal would be of course to incorporate these problems as of the compiler…
Recently, the introduction of calculations of sensitivity of mechanical thermo fields [R4.03.01] and of
rate of refund of energy [R7.02.01] compared to a variation of field, showed
relevance and the feasibility of this type of approach in Code_Aster. By coupling this last with
software PROBAN, one can thus know the probability of priming of the rupture for a distribution
of variation of field given. This type of studies mechanic-reliability engineers, for example, was carried out in
the framework of project PROMETE [bib4] to determine the probability of rupture of a tank ITEM in
considering the variability the thickness of its lining.
These sensitivities can also intervene in a crucial way in the resolution of problems opposite
(retiming of parameters…) and in many problems of optimization.
In this document, one restricts oneself with the linear and non-linear thermal problems
Code_Aster and thus with the analytical sensitivities of the field of temperature
T
(and of its flow) by
report/ratio with the characteristics material and the loadings. One places oneself in the perimeter
of use of the standard thermal operators (the loadings are supposed to be fixed, one
do not thus be interested in the phenomena of convection-dissemination in pointer of
THER_NON_LINE_MO
[R5.02.04]) for isoparametric finite elements (one does not treat it
thermal problem for the thin hulls [R3.11.01] (modeling
COQUE_ *
) and for
elements of Fourier (resp.
AXIS-FOURIER
)) (
THER_LINEAIRE
[R5.02.01] and
THER_NON_LINE
[R5.02.02]) and also in that of the operators of preprocessing of the data (
DEFI_MATERIAU
[U4.44.01],
AFFE_MATERIAU
[U4.44.03] and
AFFE_CHAR_THER
[U4.44.02]).
One is interested only in derivation of
T
and of its flow, fields depending on the variables on space
X
and of time
T
and of parameters materials and loadings, compared to one of these parameters
(which must be a constant scalar by geometrical areas (these under-parts are supposed to be distinct
and motionless, one neglects thus in particular the phenomena of dilation)). Thus let us consider, by
example, a Bi-material whose isotropic thermal conductivity is a constant reality by areas:
1
on
1
and
2
on
2
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
5/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A








Appear 1-a: Designation of thermal conductivities for a Bi-material
One thus modelizes total thermal conductivity in the form:
()
()
() (
)
2
2
1
2
2
1
1
,
:
+
=
X
X
X
I
I
éq
1-1
with
I
I
the indicating function of the ième part (
()


=
if not
0
if
1
:
I
I
X
I X
). One is interested in the sensitivity
field of temperature compared to one of the two parameters
()
T
T
I
I
I
,
0
X
=




. The problem
formulate same manner for a loading or a limiting condition.
By parameterizing advisedly the “derived loadings and materials” (cf [§6.2], [§6.4]) calculations
developed thereafter and the data-processing developments that they subtend can too
to take into account the modelings more sophisticated with several space dependences and
temporal. For example, considering a thermal source,
()
() ()
() ()
T
I
S
T
I
S
T
S
,
,
:
,
2
2
2
1
1
1
X
X
X
X
X
+
=
éq
1-2
one can calculate
()
T
S
T
I
I
S
S
I
,
0
X
=




while parameterizing


=
=
2
1
1
1
1
1
on
0
on
I
S
S
and


=
=
2
2
1
2
2
2
on
on
0
I
S
S
in the definition of the derived sources. In the remainder of the document,
we restrict ourselves with the first modeling [éq 1-1], in order not to overload them
later theoretical developments, but also because it appears closer to the real needs
users. On a case-by-case basis, we specify however more sophisticated derivations which are
accessible taking into account the new introduced functionalities and the perimeter from use of the code.
In order to more easily be able to switch over space or temporal derivation with derivation by
report/ratio with one of the parameters, one works with a derivation “within the meaning of the distributions” on
initial parabolic problem (derivation clarified in appendix 1). But the same exercise could
to be carried out starting from its version semi-discretized in time, of the variational formulation or of
linear system (one then takes the “discrete derivative” i.e. compared to the components of
discretized parameter) resulting from its discretization. In linear thermics, one shows that these
derivations, to each stage of the numerical process, lead to the same result. The problem
derived discretized being identical to the discrete derived problem, theoretical results exhumed on
continuous problem can apply to the problem actually implemented.
1
1
2
2
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
6/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
The approach adopted for the calculation of the sensitivity is thus purely analytical. Him could
to prefer an semi-analytical method (it is the step partly retained by the codes
MSC/Nastran, FEMtools and ABAQUS (cf [§9] [bib8])) : broadly analytical for the determination
derived problem and using differences finished at the local level to determine the derivative of
matrices and of the second elementary derived members. But the latter, although easier with
to establish and to maintain, is also more expensive and it introduces a dependence with respect to
parameter of shift.
For linear thermics, the “derived” problem is very similar to the initial problem.
first member of the resulting linear system is preserved. He does not have thus to be reassembled,
only the second member is to be packed by suitable a source term. Theoretical results of existence,
of unicity and convergence of the solution are not appreciably modified. Obtaining the derivative
in temperature require the same processes numerical (dualisation and inversion of the system
linear resulting).
On the other hand, in non-linear thermics, the derived problem is metamorphosed: the operator
parabolic is modified. He became linear, just like the limiting conditions. These last
are more than of two types: Dirichlet or Robin, exit conditions of Neumann and radiation.
Because of its linear character, the usual theoretical results are thus much easier with
to exhume. In addition the resolution of the derived problem is faster and more robust than that of
initial problem. One does not need to have recourse to an algorithm of Newton-Raphson for
to determine the increment of temperature between two contiguous moments. A linear solvor is enough: no one is not
need to assemble a tangent matrix with each under-iteration.
This time the two members of this equation are fundamentally different from those of the problem in
temperature. However, after each pitch of time, once determined
T
+
from
T
-
, one does not have
to reassemble all the matrix of the linear system and its second associate member. It is enough to
to supplement the first tangent matrix of the pitch of time following by the term due to the non-linearity of
thermal conductivity. One also leaves the second member of the problem in temperature for
to constitute that which interests us: one packs it by the terms of implicitation of non-linearities of
thermal conductivity and of the limiting conditions.
One details theoretical, numerical work and data processing which governed the establishment of these
calculations of sensitivities in the code. One specifies their properties and their limitations while connecting these
considerations with a precise parameter setting of the accused operators and with the choices of modeling of
code. One tried constantly to bind different the items approached while detailing, has minimum, them
a little technical demonstrations.
In short, the perimeter of use of this functionality gathers it or not thermal, linear,
isotropic or anisotropic, stationary or transitory, being pressed on isoparametric finite elements
lumpés or not. Within this framework there, it covers the same perimeter as that with the operators
accused thermics.
The request for one or more sensitivities does nothing but enrich the structure of data
thermics (
EVOL_THER
) and provides also the thermal field of which they are the derivative. In term
of performance, the calculation of an analytical sensitivity is much less expensive than a calculation
standard since the same factorized matrix is re-used.
In addition to the calculation of sensitivities in thermics, Code_Aster proposes their hanging in mechanics
statics or quasi-static [R4.03.03] and in dynamics [R4.03.04]. All these functionalities and theirs
postprocessings associated (impressions, tests…) are included in the user's documentation
[U4.50.02] and belong to deliverable project “Uncertainties of numerical calculations” [bib5].
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
7/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
This document is articulated around the following parts:
·
initially, one is interested in the various problems resulting from
derivation of the linear problem of thermics (
THER_LINEAIRE
) compared to the parameters
characteristics material and loadings,
·
then one reiterates this process on the nonlinear problem of thermics (
THER_NON_LINE
),
·
various linear systems “derived”, - directly plantable in Code_Aster
to determine such or such sensitivity -, are recapitulated in the third part,
·
in the following paragraph, one describes necessary postprocessings to obtain the sensitivities
heat transfer rate (
CALC_ELEM/CALC_NO
),
·
one concludes by approaching the practical difficulties from implementation, the environment, it
parameter setting and the perimeter of use. An example of use extracted from an official case-test
(
SENST04A
) is also detailed.

Warning:
The reader in a hurry and/or not very interested by the theoretical springs genesis of these
sensitivities and the details of modeling of the code can, from the start, to jump to [§4] and [§6] which
the main theoretical and practical contributions recapitulate preceding chapters.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
8/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
2 Thermics
linear
One considers a limited open body occupying related
R
Q
(
Q
=2 or 3) of border
lipschitzienne characterized by its voluminal heat with constant pressure
C
p
(X)
(the variable
vectorial
X
symbolize the couple here
(X, y)
(resp.
(X, y, Z)
) for
Q
=2 (resp.
Q
=3)) and its coefficient of
isotropic thermal conductivity
(X)
. These data materials are supposed to be independent of
time (modeling
THER
of Code_Aster) and constants by element (discretization
P
0
).
Note:
With modeling
THER_FO
these characteristics can depend on time. As of
first versions of the code and before the installation of
THER_NON_LINE
, it allowed
to simulate “pseudo” non-linearities. Taking into account its rather marginal use, us
we will not interest, initially, with its derivation.
One is interested in the changes of the temperature in any point
X
opened and at any moment
[[
(
)
0
,
0
>
T
, when the body is subjected to limiting conditions and loadings
independent of the temperature but being able to depend on time. It is about voluminal source
S (X, T),
boundary conditions of imposed the temperature type
F (X, T)
(on the external portion of surface
1
),
imposed normal flow
G (X, T)
(on
2
) and exchanges convectif
H (X, T)
and
T
ext.
(X, T)
(on
3
).
One places oneself thus within the framework of application of the operator
THER_LINEAIRE
[R5.02.01] of
Code_Aster by retaining only the conductive aspects of this linear thermal problem.
This problem in extreme cases interfered (type Cauchy-Dirichlet-Neumann-Robin (also called condition of
Fourier) inhomogenous, linear and with variable coefficients) is formulated
(
)
] [
] [
] [
] [






=
×
=
+
×
=
×
=
×
=
-
)
(
)
0
,
(
,
0
,
0
,
0
,
0
div
0
3
2
1
X
X
T
T
HT
HT
N
T
G
N
T
F
T
S
T
T
T
C
ext.
p
éq
2-1
Note:
·
The condition of Robin modelizing the convectif exchange (key word
EXCHANGE
) on a portion of
edges of the field, can be duplicated to take account of exchanges between two pennies
parts of the border in opposite (key word
ECHANGE_PAROI
). This limiting condition
modelize a thermal resistance of interface
] [
] [


×
=
+
×
=
+
=
=
,
0
,
0
,
21
1
2
2
12
2
1
1
21
12
3
HT
HT
N
T
HT
HT
N
T
T
T
ij
I
has
one
With
éq
2-2
·
The condition of Dirichlet can spread in the form of linear relations between the ddls
(key word
LIAISON_ *
) to simulate, in particular, of geometrical symmetries of
structure.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
9/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
()
() ()
] [
()
()
] [
DLL)
(LIAISON_D
ROUP)
(LIAISON_G
,
0
,
,
,
0
,
,
,
,
1
1
2
2
1
1
21
12
1
×
=
×
=
+
=
=
on
T
T
T
on
T
T
T
T
T
T
T
I
I
I
J
J
J
I
I
I
I
ij
X
X
X
X
X
simply
more
or
has
one
With
éq
2-3
·
One will not speak about the functionalities
LIAISON_UNIF
and
LAISON_CHAMNO
who allow
to impose the same temperature (unknown) on a whole of nodes, because they are not
that a surcouche of the preceding conditions imposing of the couples (
,
)
private individuals.
·
When the material is anisotropic (modeling
THER_ORTH
), conductivity is
modelized by a diagonal matrix expressed in the reference mark of orthotropism of material.
That basically does not change calculations according to which only hold account
isotropic case. Guard should just be taken not to switch over more, under the conditions limit
Neumann and of Robin, the scalar product with the normal and multiplication by
conductivity. In practice, in elementary calculations, one is not interested in the derivative
normals. The problem thus arises only in the preliminary theoretical part.
The sensitivity compared to one of the components of anisotropic conductivity is not
not yet available. These calculations were set up in the accused subroutines
(elementary calculations YOU.), they await nothing any more but the software evolution consisting with
to extend the taking into account of the anisotropy to the functions (a modeling
THER_ORTH_FO
). Indeed, from a point of view (cf structures.
DEFI_PARA_SENSI
[bib6]),
variable ASTER representing the significant parameter must be a data-processing object of
function type.
·
In all following calculations of sensitivities, one calculates only it derived by report/ratio
with a constant parameter by area. If not, it would be necessary to introduce a concept of derivative
directional!
This does not exclude a temporal or space dependence from characteristics material or
loadings. By parameterizing advisedly the “derived” loadings and materials
in the command file, one can also have access to some derivatives made up
(cf [§6.2]/[§6.4]).
·
For a transitory calculation, the initial temperature can be selected in three manners
different: by carrying out a stationary calculation over the first moment, by fixing it at one
uniform or unspecified value created by one
AFFE_CHAM_NO
and by carrying out a recovery with
to start from a preceding transitory calculation. This choice will affect the initialization of
derived problem.
·
We will not treat the case where (almost) all the loadings are multiplied by one
even function dependant on time (option
FONC_MULT
(this well adapted functionality
for certain mechanical problems is disadvised in thermics, because it can re-enter in
conflict with the temporal dependence of the loadings and, in addition, it applies
selectively with each one of them. It was not included besides in
THER_NON_LINE
)).
In order to be able to consider different the derivative from the temperature in configurations
multimatériau and multichargement one introduces the following notations:
(
)
3
or
2
1,
I
=
=
=
U
U
I
ij
I
I
I
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
10/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
The characteristics material and the loadings will be subscripted to notify their membership of such
or such opened or portion of border. Thus, if one is interested in a Bi-material,
1
modelize heat
voluminal of material occupying the open first
1
and
2
that of material occupying the second
opened
2
.







Appear 2-a: Designation of voluminal heats for a Bi-material

2.1
Derived compared to voluminal heat
2.1.1 Elements
theoretical
One modelizes total voluminal heat penny the form
()
() (
)
=
I
I
I
I
p
I
C
X
X
. All them
opened being solidified, one has
I
I
p
I
C
=
the indicating function of the ième part
I
. The derivation of
[éq 2-1] then leads us “trivialement” (cf Appendix 1) to the new problem in extreme cases of which is
solution sought sensitivity, noted
I
T
U
=
,
(
)
] [
] [
] [
] [






=
×
=
+
×
=
×
=
×
=
-
0
3
2
1
~
)
0
,
(
,
0
~
,
0
~
,
0
~
,
0
~
div
U
U
H
hu
N
U
G
N
U
F
U
S
U
T
U
C
p
X
éq
2.1.1-1
with the new source voluminal and the new limiting conditions and initial
0
~
~
,
0
~
~
~
0
=
-
=
=
=
=
U
T
T
I
S
H
G
F
I
and
éq
2.1.1-2
One is thus brought to solve a homogeneous problem in
U
similar to that which answers
T
. One fixes
a pitch of time
T
such as
T
that is to say an entirety
NR
. Semi-discretization in times of [éq 2.1.1-1],
[éq 2.1.1-2] by
- method leads to the following problem: to find a continuation.
1
1
2
2
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
11/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Note:
By using an adaptation of the theorem of Lax-Milgram to the parabolic problems
([R4.10.03 §1] or [bib1] theorems 1 & 2 chap XVIII pp615-620 or [bib3] pp220-241) one
show, under certain conditions of regularity on opened, materials, the loadings and
the initial condition, that this problem admits a single solution.
()
()
{
}
0
/
1
1
0
0
=
=
U
H
U
V
U
NR
N
N
éq
2.1.1-3
This space comprises also the conditions of “generalized” Dirichlet of linear relations type
between ddls when they exist.
such as:
(
)
(
)
(
)






=
-
=
+
-
=
-
=
-
-
=
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
0
0
3
1
1
1
1
2
1
1
1
1
1
1
1
1
~
)
(
1
0
~
1
0
~
1
0
~
1
0
~
~
div
1
div
U
U
NR
N
H
U
H
N
U
NR
N
G
N
U
NR
N
F
U
NR
N
T
S
S
U
U
T
U
U
C
N
N
N
N
N
N
N
N
N
N
N
N
N
N
p
X
éq 2.1.1-4
while posing:
()


=
=
=
=


-
=


=
NR
N
H
H
H
G
F
NR
N
T
C
S
NR
N
U
U
N
N
N
N
I
p
N
N
,
0
~
~
~
,
,
~
,
,
X
X
X
X
and
éq 2.1.1-5
By applying the theorem of Green to [éq 2.1.1-3], [éq 2.1.1-4], [éq 2.1.1-5] and by introducing them
following notations
(
)
{
}
1
0
,
,
,
1
,
-


=


+
=
-
+
NR
N
H
T
U
NR
N
NR
N
and
with
and
X
X
one is brought to solve the following variational problem:
()
()




=
±
+
+
+
-
±
±
v
L
v
U
has
V
v
V
U
U
T
H
,
0
0
,
that
such
To calculate
and
given
Being
éq
2.1.1-6
with the bilinear form depending on the current moment (via
H
+
)
()
+
+
+
+
+
+
+
+
=
3
1
,
D
v
U
H
dx
v
U
dx
v
U
C
T
v
U
has
p
éq
2.1.1-7
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
12/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
and the linear form parameterized by the moments running and precedent (via
H
,
U
and
T
)
()
(
)
(
)
(
)
+
-
-
-
-
-
±
-
+
-
+
-
+
=
dx
v
T
T
I
T
D
v
U
H
dx
v
U
dx
v
U
C
T
v
L
I
p
1
1
1
1
3
éq 2.1.1-8
Note:
·
Contrary to the initial problem, the unknown field and the function test belong to same
functional space, which is more comfortable from a theoretical and numerical point of view.
·
The first member of this equation is formally identical to that of the equation in
temperature. After each pitch of time, once determined
T
+
from
T
-
, one rests
also on
H
-
, H
+
and
U
-
to determine
U
+
.
The matrix of the linear system corresponding does not have
not with being reassembled. Only the second member is to be packed by suitable the source term.
·
By deriving the variational formulation (cf [§5.1.3]) from the problem in temperature [R5.01.02]
one finds well [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8].
·
For a condition of exchange between walls, the term of usual exchange is of course replaced
by (by taking again the notations of [éq 2-2])
(
)
±
±
±
-
ij
D
v
U
U
H
J
I
.
·
These problems of suitable initialization of the problem derived are recurring. It was found
in the calculation of the derivative Lagrangian of the field of temperature compared to one
variation of field (cf [R4.03.01] operand
DEUL_INIT
).
·
The reader interested by a theoretical study of the thermal problem actually put in
place in the code, which underlines these holding and bordering and their links with the choices of
modeling, will be able to refer to [§1] Doc. R: “Indicating of error in residue for
transitory thermics “[R4.10.03]. It relates to a related field of improvement and of
calibration of the studies, that of the space errors due to the mesh finite elements.
To solve this problem numerically one spatially discretizes it by considering one
subspace
H
V
0
of
0
V
of dimension finished
()
{
}
K
U
K
V
U
V
NR
U
U
K
K
H
H
H
H
I
N
I
I
H
=
=
±
±
=
±
±
/
:
0
0
1
while noting (
T
H
)
H
a regular family of triangulations of the polygonal or polyhedric field discretized
H
,
P
K
(
K)
the space of the polynomials of degree <
(k+1)
on
K
and
NR
I
the function of form associated with
node n°i. From where the discretized variational problem
Note:
This property of total continuity of the elements and maximization of their characteristics
geometrical (which ensures the convergence of the finite element method) is checked for
all isoparametric elements of the code: segment, triangle, quadrangle, tetrahedron,
pentahedron and hexahedron.
(
)
()




=
±
+
+
+
-
±
±
H
H
H
H
H
H
H
H
H
H
H
H
v
L
v
U
has
V
v
V
U
U
T
H
,
0
0
,
that
such
To calculate
and
given
Being
éq
2.1.1-9
leading to the linear system of command
N.
L
U
With
=
+
éq
2.1.1-10
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
13/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
p
boundary conditions of the Dirichlet type are taken into account in Code_Aster by one
technique of double dualisation [R3.03.01] via “ddls (degrees of freedom)” of Lagrange
()
p
I
I
,
1
=
=
. That is to say
J
the whole of the nodes belonging to the portion of border on which
apply the condition of Dirichlet (card (
J
) =
p
), the new unknown vector is considered
[
]
T
+
=
U
U
~
and the operator
B
(of command
p
X
N
) checking
(
)
J
I
U
I
I
=
+
+
with
U
B
The homogeneous condition of Dirichlet of [éq 2.1.1-1], [éq 2.1.1-2] is realized while imposing
0
C
U
B
=
=
+
The dualized problem then consists in reversing the system of command
n+2p




=








-
-
=
+
+
0
0
L
U
Id
Id
B
Id
Id
B
B
B
With
L
U
With
T
T
~
~
~
éq
2.1.1-11
Note:
·
One carries out the possible taking into account of limiting conditions of Dirichlet generalized like
in the problem in temperature (here
p
=1) but with a second null member
0
=
=
=
+
+
I
J
J
J
J
J
U
C
U
B
and
with
We now will see how these calculations are declined in the code.

2.1.2 Establishment in Code_Aster
The matrix of this system results from the assembly of the following elementary terms, due to
contribution of the nodes
(
I, J)
at the point of gauss
G
(of weight
G
(this weight gathers in fact “truth
weight “of the formula of quadrature multiplied by the jacobien of the element considered and possibly
by the radius of the point of gauss
R
G
(in modeling
AXIS
or
AXIS_DIAG
))) current element
K
(in Code_Aster the conditions limit are affected on particular elements of skin of
dimension
q-1
. No confusion not being really possible, one will make formal distinction here
between those and the elements of volume which support them).
(
)
(
)
=
=
3
1
,
,
I
G
I
ij
G
ij
K
With
K
With
éq
2.1.2-1
with
·
the thermal term of mass (calculated by the option
MASS_THER)
(
)
()
() ()
G
I
G
J
p
G
G
ij
NR
NR
K
C
T
K
With
=
,
1
·
the thermal term of rigidity (
RIGI_THER)
(
)
()
()
()
G
I
G
J
G
G
ij
NR
NR
K
K
With
=
.
,
2
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
14/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
·
the term of rigidity due to the conditions limit of exchange (
RIGI_THER_COEF_F/R)
(
)
(
)
() ()
G
I
G
J
G
G
ij
NR
NR
K
H
K
With
3
3
3
,
=
+
In the event of exchange between walls this term is replaced by (
RIGI_THER_PARO_F/R
)
(
)
(
)
()
()
()
(
) ()
G
I
G
J
F
G
J
G
G
ij
NR
NR
NR
K
H
K
With
-
=
+
3
3
3
,
while noting
F
the bijection putting in opposite the two walls.
The second member is written, with the same notations,
(
)
(
)
=
=
3
1
,
,
I
G
I
J
G
J
K
L
K
L
éq
2.1.2-2
where
·
the term resulting from the implicitation of the matrix of rigidity and mass (new option
CHAR_SENS_EVOL
, copy of
CHAR_THER_EVOL
with
U
-
instead of
T
-
, the field material
derived and standard and the new source term)
(
)
()
() ()
(
)
()
()
()
G
J
G
G
G
J
G
p
G
G
J
NR
U
K
NR
U
K
C
T
K
L
-
+
=
-
-
.
1
,
1
·
the term resulting from the implicitation of the conditions limit of exchange
(
CHAR_THER_TEXT_F/R
with
T
ext.
=0
and
U
-
instead of
T
-
)
(
)
() (
)
() ()
G
J
G
G
G
J
NR
U
K
H
K
L
-
-
-
=
3
3
2
1
,
In the event of exchange between walls this term is replaced by (
CHAR_THER_PARO_F/R
with
U
-
instead of
T
-
)
(
)
() (
)
()
()
(
)
(
) ()
G
J
G
G
G
G
J
NR
F
U
U
K
H
K
L
-
-
-
-
-
=
3
3
2
1
,
·
the term due to the “new source” comprising the field derived material (cf.
CHAR_SENS_EVOL
above)
(
)
()
() ()
(
) ()
G
J
G
G
I
G
G
J
NR
T
T
K
I
T
K
L
-
+
-
-
=
,
3
As one already specified all the elementary terms of the matrix are the subject of an option
of calculation and will already have been evaluated for the calculation of
T
+
.
It thus remains to estimate the second member
by re-using (with a different parameter setting) the existing options of calculation or by introducing one
news (
CHAR_SENS_EVOL
). This new option is common with the other derived material
(thermal conductivity) and it redirects towards the same routine of elementary calculation (
TE
.). The channel
of character power station (
SENS
instead of
THER
) joined to a detection of the nullity of the field material
derived, allows to parameterize this routine towards one of its three possible orientations: calculation of
term of implicit mass and standard rigidity, idem of sensitivity compared to one of both
characteristics material which thus adds a new source term.
In accordance with the principles of architecture set up in the code to treat calculations
of sensitivity [bib6], the assembly and the resolution of [éq 2.1.1-11] are started by the analysis of
table of correspondence associated with the significant variable. It was seen that this calculation is very close to one
standard linear thermal calculation, only the initial condition and the loadings are modified
0
,
0
,
0
,
0
0
=
-
=
=
=
=
U
T
T
I
S
T
G
F
I
ext.
and
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
15/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
This feedback is carried out via the succession of controls
v = DEFI_PARA_SENSI (VALE = < value of
I
>)
my = DEFI_MATERIAU (THER = _F
(RHO_CP = v))
affe = AFFE_MATERIAU (AFFE = _F
(GROUP_MA = < definition of
I
>,
MATER =
my))
one = DEFI_CONSTANTE (VALE = 1. )
MEMO_NON_SENSI (NAME =_F
(NOM_SD = “my”, PARA_SENSI = “v”,
NOM_COMPOSE = `
ma_v'))
ma_v = DEFI_MATERIAU (THER = _F
(RHO_CP = one))
MEMO_NON_SENSI (NAME = _F
(NOM_SD = “affe”, PARA_SENSI = “v”,
NOM_COMPOSE
= `
affe_v'
))
affe_v = AFFE_MATERIAU (AFFE = _F
(GROUP_MA = <
I
>, MATER =
ma_v))
resu = THER_LINEAIRE (CHAM_MATER = affe,
SENSITIVITY
=
(
v)
…)
Note:
·
In his command file, the user will not have soon any more but to specify the first and it
third blocks of instruction. The block of the medium will be generated automatically by
supervisor thanks to the shaft of dependence which it builds between the various controls.
·
The essential data of this calculation, derived voluminal heat
I
p
I
C
I
=
, is provided by
ma_v
.
·
To take into account a more sophisticated modeling of voluminal heat
()
() () (
)
=
I
I
I
I
I
p
I
C
X
X
X
in these calculations of sensitivity, it is enough to substitute for the indicating function the new term
source
I
I
I
p
I
C
=
in the definition of the field derived material
ma_v
.
We will unroll the same process for the various sensitivities, to start with that
concerning the other characteristic material: thermal conductivity.

2.2
Derived compared to thermal conductivity
2.2.1 Elements
theoretical
One poses
()
() (
)
=
I
I
I
I
I
X
X
and
I
T
U
=
required sensitivity. All the open ones
being solidified, one has
I
I
I
=
the indicating function of the ième part
I
. The derivation of [éq 2-1] us
conduit with a problem in extreme cases identical to [éq 2.1.1-1] but with a nonnull voluminal source
and of the new conditions of Neumann and Robin
(
)
T
I
S
N
T
I
H
G
I
I
=
-
=
=
div
~
~
~
and
éq
2.2.1-1
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
16/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
It is a homogeneous problem of Cauchy-Dirichlet and inhomogenous Neumann-Robin similar to that
which answers
T
. Its semi-discretization in time led to seek a continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 2.1.1-4] whose first relation is rewritten
(
)
(
)
(
)
(
)
1
0
~
1
~
div
1
div
1
1
1
-
-
+
=
-
-
-
-
+
+
+
NR
N
S
S
U
U
T
U
U
C
N
N
N
N
N
N
p
éq 2.2.1-2
with the new source term
(
)
{
}
1
,
div
~
+
=
N
N
m
T
I
S
m
I
m
éq
2.2.1-3
and new limiting conditions
N
T
I
H
G
N
I
N
N
-
=
=
+
+
1
1
~
~
éq
2.2.1-4
From where a variational problem identical to [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] comprising
even bilinear form [éq 2.1.1-8] joined to the linear form [éq 2.1.1-7] of which only the fourth
integral is modified to adapt to the new source
()
(
)
(
)
-
+
±
-
+
-
=
dx
v
T
T
I
v
L
I
1
L
éq
2.2.1-5
Note:
·
The first member of this equation is formally identical to that of the equation in
temperature. After each pitch of time, once determined
T
+
from
T
-
, one rests on
H
-
, H
+
and
U
-
to determine
U
+
.
The matrix of the linear system corresponding does not have to be
reassembled. Only the second member is to be packed by suitable the source term.
·
By deriving the variational formulation (cf [§5.1.3]) from the problem in temperature [R5.01.02]
one finds well [éq 2.2.1-5].
·
In stationary regime this complementary source term is tiny room to
()
-
=
dx
v
T
I
v
L
I
.
L
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].

2.2.2 Establishment in Code_Aster
Compared to [§2.1.2], only the term due to the new source is to be modified
(
)
()
()
(
)
()
(
)
()
G
J
G
G
I
G
G
J
NR
T
T
K
I
K
L
-
+
-
=
-
+
1
,
3
éq
2.2.2-1
What is made in the new option of calculation
CHAR_SENS_EVOL
with the field derived material and
standard.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. It was seen that this calculation is very close to a calculation
standard linear thermics, only the initial condition and the loadings are modified
(
)
0
div
,
,
0
0
=
=
-
=
=
=
U
T
I
S
N
T
I
HT
G
F
I
I
ext.
and
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
17/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
This feedback is carried out via the same succession of controls as with [§2.1.2] in
substituent
LAMBDA
(or
LAMBDA_L/T/NR
into orthotropic) with
RHO_CP
in
DEFI_MATERIAU
.
Note:
·
The essential data of this calculation, derived thermal conductivity
I
I
I
=
, is provided
by
ma_v
.
·
When the material is anisotropic the thermal matrix of conductivity is expressed in
identify orthotropism of material: it is thus diagonal (in our example
Q
=2)


=
2
1
0
0
One can then dissociate derivation compared to a value of one of his diagonal terms
derivation compared to a value of this diagonal. The formulas detailed here are
identical whatever the configuration selected. Only the evaluation of
kl
ij
K
ij
or
must
to take account of these characteristics.
·
In practice, one has access to the sensitivity compared to an isotropic conductivity
constant by area. Sensitivity compared to a component of conductivity
anisotropic is not yet available. These calculations were set up in
subroutines accused (elementary calculations YOU.), they await nothing any more but the evolution
software consisting in extending the taking into account of the anisotropy to the functions (one
modeling
THER_ORTH_FO
). Indeed, from a point of view structures
(cf.
DEFI_PARA_SENSI
[bib6]), the variable ASTER representing the significant parameter
must be a data-processing object of function type.
·
It was reasoned here as if the condition of initial Cauchy of the problem had been fixed
uniform or unspecified. If it is given by carrying out a stationary calculation on
first moment it is necessary to reiterate this process with the derived problem. On the other hand, if it results
of a recovery starting from a preceding transitory calculation, the derived problem must be initialized with
to leave the value of the same derivative at the same moment of recovery.
In short, two initializations (that of the problem in temperature and that of the problem
derived) must be homogeneous. On the other hand, contrary to a calculation of thermics
standard, one cannot thus modify the conditions limit and it is necessary to carry out a recovery with
to start from a calculation of comparable nature.
·
To take into account a more sophisticated modeling of thermal conductivity
()
() () (
)
=
I
I
I
I
I
I
X
X
X
in these calculations of sensitivity, it is enough to substitute for the indicating function the new one
source term
I
I
I
I
=
in the definition of the field derived material
ma_v
.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
18/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
2.3
Derived compared to the source
2.3.1 Elements
theoretical
One poses
()
() (
)
=
I
I
I
I
S
I
S
S
X
X
and
I
S
T
U
=
required sensitivity. All the open ones being
solidified, one has
I
I
I
S
S
=
the indicating function of the ième part
I
. The derivation of [éq 2-1] leads us
with a problem in extreme cases identical to [éq 2.1.1-1] but with another voluminal source
I
I
S
=
~
éq
2.3.1-1
It is a homogeneous problem similar to that which answers
T
. Its semi-discretization in time
conduit to seek a continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 2.1.1-4] of which
first relation is rewritten
(
)
(
)
(
)
1
0
~
div
1
div
1
1
-
=
-
-
-
-
+
+
NR
N
S
U
U
T
U
U
C
N
N
N
N
p
éq
2.3.1-2
From where a variational problem identical to [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] comprising
even bilinear form [éq 2.1.1-7] joined to the linear form [éq 2.1.1-8] of which only the fourth
integral is modified to adapt to the new source
()
±
+
=
dx
v
I
v
L
I
L
éq
2.3.1-3
Note:
·
The first member of this equation is formally identical to that of the equation in
temperature. After each pitch of time, once determined
T
+
from
T
-
, one rests
also on
H
-
, H
+
and
U
-
to determine
U
+
.
The matrix of the linear system corresponding does not have
not with being reassembled. Only the second member is to be packed by suitable the source term.
·
By deriving the variational formulation (cf [§5.1.3]) from the problem in temperature [R5.01.02]
one finds well [éq 2.3.1-3].
·
In stationary regime this complementary source term is not modified.
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].

2.3.2 Establishment in Code_Aster
Compared to [§2.1.2], only the term due to the new source is to be modified
(
)
()
()
G
J
I
G
G
J
NR
K
I
K
L
=
,
3
éq
2.3.2-1
It is enough to re-use the standard option
CHAR_THER_SOUR_F/R
with the derived field source.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. It was seen that this calculation is very close to a calculation
standard linear thermics, only the initial condition and the loadings are modified
0
,
0
0
=
=
=
=
=
U
I
S
T
G
F
I
ext.
and
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
19/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
This feedback is carried out via the succession of controls
v = DEFI_PARA_SENSI (VALE = < value of
S
I
>)
chth = AFFE_CHAR_THER_F (SOURCE = _F
(GROUP_MA = < definition of
I
>,
SOUR =
v))
one = DEFI_CONSTANTE (VALE = 1. )
MEMO_NON_SENSI (NAME =_F
(NOM_SD = “chth”, PARA_SENSI = “v”,
NOM_COMPOSE = `
chth_v'))
chth_v = AFFE_CHAR_THER_F (SOURCE = _F (GROUP_MA = <
I
>, SOUR =
one))
resu = THER_LINEAIRE (EXCIT = chth,
SENSITIVITY
=
(
v)
…)
Note:
·
The essential data of this calculation, the derived field source
I
I
S
S
I
=
, is provided by
chth_v
.
·
This calculation is independent of the three types of modeling of the source: constant by mesh
(
AFFE_CHAR_THER + SOUR
), constant by point of Gauss (
AFFE_CHAR_THER +
SOUR_CALCULEE
) and constant by mesh and dependant on time (
AFFE_CHAR_THER_F +
SOUR
). These considerations do not even re-enter in line of account during the effective calculation of
I
S
S
in
chth_v
, because this size represents the derivation of a function parameterized by
one of its constant parameters. One is not interested here in derivative of the type
()
()
()
()
T
T
S
S
U
T
S
S
J
I
J
G
I
,
,
,
X
X
X
O
.
·
To take into account a more sophisticated modeling of the source
()
() () (
)
=
I
I
I
I
I
S
T
I
S
T
S
,
,
X
X
X
in these calculations of sensitivity, it is enough to substitute for the indicating function the new term
source
I
I
I
I
S
S
=
in the definition of the derived loading
chth_v
.
2.4
Derived compared to the imposed temperature
2.4.1 Elements
theoretical
One poses
()
() (
)
=
I
I
I
I
F
I
F
F
X
X
and
I
F
T
U
=
required sensitivity. Portions of
external border
J
1
being solidified, one has
I
I
I
F
F
=
the indicating function of the ième portion
I
1
.
derivation of [éq 2-1] leads us to a problem in extreme cases identical to [éq 2.1.1-1] but with one
another voluminal source and a new condition of Dirichlet
0
~
~
=
=
S
I
F
I
and
éq
2.4.1-1
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
20/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
It is an inhomogenous problem of Dirichlet and Cauchy- Neumann-Robin homogeneous similar to that
which answers
T
. Its semi-discretization in time led to seek a continuation (this space can
to comprise, if necessary, also conditions of “generalized” Dirichlet of linear relations type
between ddls)
()
()
{
}
I
NR
N
N
I
U
H
U
V
U
=
=
1
/
1
1
0
éq
2.4.1-2
checking a system similar to [éq 2.1.1-3] whose first relation is rewritten
(
)
(
)
(
)
1
0
0
div
1
div
1
1
-
=
-
-
-
-
+
+
NR
N
U
U
T
U
U
C
N
N
N
N
p
éq
2.4.1-3
with the new limiting condition
I
N
I
F
=
+
1
~
éq
2.4.1-4
From where a variational problem identical to [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] with the same form
bilinear [éq 2.1.1-7] joined to the linear form [éq 2.1.1-8] whose fourth integral is null.
Note:
·
The first member of this equation is formally identical to that of the equation in
temperature. After each pitch of time, once determined
T
+
from
T
-
, one rests
also on
H
-
, H
+
and
U
-
to determine
U
+
.
The matrix of the linear system corresponding does not have
not with being reassembled. On the other hand this time it is necessary to constitute the Lagrangian part of
second member dualized in order to approximate new functional space
V
1
.
·
By deriving the variational formulation (cf [§5.1.3]) from the problem in temperature [R5.01.02]
one finds well [éq 2.4.1-3], [éq 2.4.1-4].
Space discretization in a subspace
H
V
1
and the taking into account of the condition of Dirichlet
inhomogenous lead to a linear system dualized similar to [éq 2.1.1-11]




=








-
-
=
+
+
C
C
L
U
Id
Id
B
Id
Id
B
B
B
With
L
U
With
T
T
~
~
~
éq
2.4.1-5
with
ij
I
J
K
F
F
C
=
=
éq
2.4.1-6
while noting
F
J
the value of the condition of Dirichlet to the node n°
K
(classification room) of
1
.
2.4.2 Establishment in Code_Aster
Compared to [§2.1.2], only the second member is modified since
(
)
0
,
3
=
G
J
K
L
éq
2.4.2-1
There is not thus an option of calculation particular to envisage, it is just necessary to assemble the linear system
dualized associated the derived imposed temperatures.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. It was seen that this calculation is very close to a calculation
standard linear thermics, only the initial condition and the loadings are modified
0
0
,
0
=
=
=
=
=
U
S
T
G
I
F
ext.
I
and
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
21/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
This feedback is carried out via the succession of controls
v = DEFI_PARA_SENSI (VALE = < value of
F
I
>)
chth = AFFE_CHAR_THER_F (TEMP_IMPO = _F (GROUP_MA = < definition of
I
1
>,
TEMP =
v))
one = DEFI_CONSTANTE (VALE = 1. )
MEMO_NON_SENSI (NAME =_F
(NOM_SD = “chth”, PARA_SENSI = “v”,
NOM_COMPOSE = `
chth_v'))
chth_v = AFFE_CHAR_THER_F (TEMP_IMPO = _F (GROUP_MA = <
I
1
>, TEMP =
one))
resu = THER_LINEAIRE (EXCIT = chth,
SENSITIVITY
=
(
v)
…)
Note:
·
The essential data of this calculation, the derived field imposed temperature
I
I
F
F
I
=
, is
provided by
chth_v
.
·
This calculation is independent of the three types of modeling of the imposed temperature:
constant by mesh (
AFFE_CHAR_THER + TEMP
), more dependence in time provided by
a function (
AFFE_CHAR_THER_F + TEMP
) or provided by a structure of data
“EVOL_THER”
(
AFFE_CHAR_THER_F + EVOL_THER + “TEMP”
). These considerations
do not even re-enter in line of account during the effective calculation of
I
F
F
in
chth_v
, because this
size represents the derivation of a function parameterized by one of its parameters
constant. One is not interested here in the derivative
()
()
T
T
F
F
J
I
,
,
X
X
.
·
The calculation of derived from
T
compared to one of the parameters of the relations of Dirichlet
generalized [éq 2-3] same manner would be carried out. A distinction only appears
on the level of the Lagrangian components of the dualized system:
+
+
+
-
=
=
=
J
J
I
J
J
J
J
I
T
J
J
U
T
U
C
U
B
and
with
With
éq
2.4.2-2
I
J
J
J
I
J
J
U
T
U
=
=
=
+
+
C
U
B
and
with
With
éq
2.4.2-3
The taking into data-processing account of these calculations would be carried out, like above, via the word
keys
COEF_MULT_1/2
and
COEF_IMPO
key words factors
LAISON_GROUP
and
LIAISON_CHAMNO
.
·
Sensitivity compared to a multiplying coefficient
of this condition of Dirichlet
generalized [éq 2.4.2-2 is not available] because it has little direction with coefficients
often discrete. One does not have access, by parameterizing advisedly the condition of Dirichlet
generalized derived, that with derivation compared to the total coefficient
[éq 2.4.2-3].
·
To take into account a more sophisticated modeling of a condition of Dirichlet
()
() () (
)
=
I
I
I
I
I
F
T
I
F
T
F
,
,
X
X
X
in these calculations of sensitivity, it is enough to substitute for the indicating function the new one
source term
I
I
I
I
F
F
=
in the definition of the derived loading
chth_v
.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
22/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
2.5
Derived compared to imposed normal flow
2.5.1 Elements
theoretical
One poses
()
() (
)
=
I
I
I
I
G
I
G
G
X
X
and
I
G
T
U
=
required sensitivity. Portions
J
2
being solidified, one has
I
I
I
G
G
=
the indicating function of the ième portion
I
2
. The derivation of [éq 2-1]
us leads to a problem in extreme cases identical to [éq 2.1.1-1] but with another voluminal source
and a new condition of Neumann
0
~
~
=
=
S
I
G
I
and
éq
2.5.1-1
It is a problem of inhomogenous Neumann and homogeneous Cauchy-Dirichlet-Robin similar to that
which answers
T
. Its semi-discretization in time led to seek a continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 2.1.1-3] whose first relation is rewritten
(
)
(
)
(
)
1
0
0
div
1
div
1
1
-
=
-
-
-
-
+
+
NR
N
U
U
T
U
U
C
N
N
N
N
p
éq
2.5.1-2
with the new limiting condition
I
N
I
G
=
+
1
~
éq
2.5.1-3
From where a variational problem identical to [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] comprising
even bilinear form [éq 2.1.1-7] joined to the linear form [éq 2.1.1-8] of which only the fourth
integral is modified to adapt to the new “surface” source
()
±
+
=
2
dx
v
I
v
L
I
L
éq
2.5.1-4
Note:
·
The first member of this equation is formally identical to that of the equation in
temperature. After each pitch of time, once determined
T
+
from
T
-
, one rests
also on
H
-
, H
+
and
U
-
to determine
U
+
.
The matrix of the linear system corresponding does not have
not with being reassembled. Only the second member is to be packed by suitable the source term.
·
By deriving the variational formulation (cf [§5.1.3]) from the problem in temperature [R5.01.02]
one finds well [éq 2.5.1-4]
·
In stationary regime this source term is not modified.
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].
2.5.2 Establishment in Code_Aster
Compared to [§2.1.2], only the term due to the new source is to be modified
(
)
(
)
()
G
J
I
G
G
J
NR
K
I
K
L
2
2
3
,
=
éq
2.5.2-1
It is enough to re-use the standard option
CHAR_THER_FLUN_F/R
with the field derived flow.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. It was seen that this calculation is very close to a calculation
standard linear thermics, only the initial condition and the loadings are modified
0
,
0
0
=
=
=
=
=
U
I
G
T
S
F
I
ext.
and
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
23/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
This feedback is carried out via the succession of controls
v = DEFI_PARA_SENSI (VALE = < value of
G
I
>)
chth = AFFE_CHAR_THER_F (FLUX_REP = _F (GROUP_MA = < definition of
I
2
>,
FLUN =
v))
one = DEFI_CONSTANTE (VALE = 1. )
MEMO_NON_SENSI (NAME =_F
(NOM_SD = “chth”, PARA_SENSI = “v”,
NOM_COMPOSE = `
chth_v'))
chth_v = AFFE_CHAR_THER_F (FLUX_REP = _F (GROUP_MA = <
I
2
>, FLUN =
one))
resu = THER_LINEAIRE (EXCIT = chth,
SENSITIVITY
=
(
v)
…)
Note:
·
The essential data of this calculation, the field derived normal flow
I
I
G
G
I
=
, is provided by
chth_v
.
·
This calculation is independent of the three types of modeling of the condition of Neumann:
constant by mesh (
AFFE_CHAR_THER + FLUN
), more dependence in time provided by
a function (
AFFE_CHAR_THER_F + FLUN
) or components of vectorial flow
dependant on time and constants by meshs (
AFFE_CHAR_THER_F + FLUN_X/Y/Z
).
One can then dissociate derivation compared to the q-uplet components of derivation
compared to one of its components. The formulas detailed here are identical whatever
that is to say configuration selected. Only possibly the evaluation of
{
}
(
)
Z
y
X
L
G
G
G
G
J
L
I
K
J
L
I
,
,
or
must take account of these characteristics.
·
The dependence in time of the condition of Neumann is not taken into account at the time of
calculation of
I
G
G
because this size represents the derivation of a function parameterized by one
of its constant parameters. One is not interested here in the derivative
()
()
T
T
G
G
J
I
,
,
X
X
.
·
In practice, there is thus access to the sensitivity compared to a scalar flow or vector,
constant by area. By parameterizing advisedly the vector derived flow, one can too
to obtain the sensitivity compared to one of its components
()
T
G
T
J
I
,
X
.
·
To take into account a more sophisticated modeling of a normal flow
()
() () (
)
=
I
I
I
I
I
G
T
I
G
T
G
,
,
X
X
X
in these calculations of sensitivity, it is enough to substitute for the indicating function the new one
source term
I
I
I
I
G
G
=
in the definition of the derived loading
chth_v
.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
24/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
2.6
Derived compared to the coefficient from convectif exchange
2.6.1 Elements
theoretical
One poses
()
() (
)
=
I
I
I
I
H
I
H
H
X
X
and
I
H
T
U
=
required sensitivity. Portions
J
3
being
solidified, one has
I
I
I
H
H
=
the indicating function of the ième portion
I
3
. The derivation of [éq 2-1] us
conduit with a problem in extreme cases identical to [éq 2.1.1-1] but with another voluminal source and
a new condition of Robin
(
)
0
~
~
=
-
=
S
T
T
I
H
ext.
I
and
éq
2.6.1-1
It is an inhomogenous problem of Robin and homogeneous Cauchy-Dirichlet-Neumann similar to that
which answers
T
. Its semi-discretization in time led to seek a continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 2.1.1-3] whose first relation is rewritten
(
)
(
)
(
)
1
0
0
div
1
div
1
1
-
=
-
-
-
-
+
+
NR
N
U
U
T
U
U
C
N
N
N
N
p
éq
2.6.1-2
with the new limiting condition
(
)
1
1
1
~
+
+
+
-
=
N
N
ext.
I
N
T
T
I
H
éq
2.6.1-3
From where a variational problem identical to [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] comprising
even bilinear form [éq 2.1.1-7] joined to the linear form [éq 2.1.1-8] of which only the fourth
integral is modified to adapt to the new “surface” source
()
(
)
(
)
(
)
{
}
-
-
+
+
±
-
-
+
-
+
=
3
1
dx
v
T
T
T
T
I
v
L
ext.
ext.
I
L
éq
2.6.1-4
Note:
·
The first member of this equation is formally identical to that of the equation in
temperature. After each pitch of time, once determined
+
T
from
-
T
, one rests
also on
-
H
,
+
H
and
-
U
to determine
+
U
. The matrix of the linear system corresponding
does not have to be reassembled. Only the second member is to be packed by suitable the source term.
·
By deriving the variational formulation (cf [§5.1.3]) from the problem in temperature [R5.01.02]
one finds well [éq 2.6.1-4].
·
These calculations spread without sorrow in the condition of exchange between walls of [éq 2-2]. It
is enough to replace the terms
(
)
±
±
±
-
I
D
v
T
T
ext.
3
by
(
)
(
)
{
}
I
I
J
I
I
I
H
D
v
U
U
ij
-
-
±
±
±
1
,
,
with
.
·
In hover this complementary source term is tiny room to
()
(
)
-
+
=
3
dx
v
T
T
I
v
L
ext.
I
L
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
25/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
2.6.2 Establishment in Code_Aster
Compared to [§2.1.2], only the term due to the new source is to be modified
(
)
(
)
(
)
()
(
)
(
) (
)
()
(
) ()
G
J
G
ext.
G
ext.
I
G
G
J
NR
T
K
T
T
K
T
K
I
K
L






-
-
+
-
=
-
-
+
+
3
3
3
3
3
1
,
éq2.6.2-1
What is made in the new option of calculation
CHAR_SENS_TEXT_F
with the fields coefficient
of standard replacement and derived.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. It was seen that this calculation is very close to a calculation
standard linear thermics, only the initial condition and the loadings are modified
(
)
0
,
0
0
=
-
=
=
=
=
U
T
T
I
HT
G
S
F
ext.
I
ext.
and
This feedback is carried out via the succession of controls
v = DEFI_PARA_SENSI (VALE = < value of
H
I
>)
chth = AFFE_CHAR_THER_F (EXCHANGE = _F (GROUP_MA = < definition of
I
3
>,
COEF_H =
v, TEMP_EXT = W))
one = DEFI_CONSTANTE (VALE = 1. )
zero = DEFI_CONSTANTE (VALE = 0. )
MEMO_NON_SENSI (NAME =_F
(NOM_SD = “chth”, PARA_SENSI = “v”,
NOM_COMPOSE = `
chth_v'))
chth_v = AFFE_CHAR_THER_F (EXCHANGE = _F (GROUP_MA = <
I
3
>,
COEF_H =
one, TEMP_EXT = zero))
resu = THER_LINEAIRE (EXCIT = chth,
SENSITIVITY
=
(
v)
…)
Note:
·
The essential data of this calculation, the field coefficient of derived exchange
I
I
H
H
I
=
, is
provided by
chth_v
.
·
To take into account a more sophisticated modeling of a coefficient of exchange
()
() () (
)
=
I
I
I
I
I
H
T
I
H
T
H
,
,
X
X
X
in these calculations of sensitivity, it is enough to substitute, with the indicating function, the new one
source term
I
I
I
I
H
H
=
in the definition of the derived loading
chth_v
.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
26/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
2.7
Derived compared to the outside air temperature
2.7.1 Elements
theoretical
In the case of convectif exchange with the external medium, one poses
()
()
(
)
=
I
ext.
I
I
I
ext.
ext.
T
I
T
T
X
X
and
I
ext.
T
T
U
=
required sensitivity. As previously
I
I
ext.
ext.
I
T
T
=
the indicating function of
the ième portion
I
3
. The derivation of [éq 2-1] leads us to a problem in extreme cases identical to
[éq 2.1.1-1] but with another voluminal source and a new condition of Robin
0
~
~
=
=
S
I
H
H
I
and
éq
2.7.1-1
It is an inhomogenous problem of Robin and homogeneous Cauchy-Dirichlet-Neumann similar to that
which answers
T
. Its semi-discretization in time led to seek a continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 2.1.1-3] whose first relation is rewritten
(
)
(
)
(
)
1
0
0
div
1
div
1
1
-
=
-
-
-
-
+
+
NR
N
U
U
T
U
U
C
N
N
N
N
p
éq
2.7.1-2
with the new limiting condition
I
N
N
I
H
H
1
1
~
+
+
=
éq
2.7.1-3
From where a variational problem identical to [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] comprising
even bilinear form [éq 2.1.1-7] joined to the linear form [éq 2.1.1-8] of which only the fourth
integral is modified to adapt to the new “surface” source
()
(
)
{
}
-
+
±
-
+
+
=
3
1
dx
v
H
H
I
v
L
I
L
éq
2.7.1-4
Note:
·
The first member of this equation is formally identical to that of the equation in
temperature. After each pitch of time, once determined
+
T
from
-
T
, one
also rest on
-
H,
+
H
and
-
U
to determine
+
U
. The matrix of the linear system
corresponding does not have to be reassembled. Only the second member is to be packed by the term
suitable source.
·
By deriving the variational formulation (cf [§5.1.3]) from the problem in temperature [R5.01.02]
one finds well [éq 2.7.1-4].
·
In stationary regime this source term is tiny room to
()
+
=
3
dx
v
H
I
v
L
I
L
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
27/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
2.7.2 Establishment in Code_Aster
Compared to [§2.1.2], only the term due to the new source is to be modified
(
)
(
)
(
) (
) (
)
{
} ()
G
J
I
G
G
J
NR
K
H
K
H
K
I
K
L
3
3
3
3
3
1
,
-
+
=
-
+
éq
2.7.2-1
It is enough to re-use the standard option
CHAR_THER_TEXT_F
with the fields coefficient of exchange
standard and derived, and while it “bluffant” with a field T
-
=0.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. It was seen that this calculation is very close to a calculation
standard linear thermics, only the initial condition and the loadings are modified
0
,
0
0
=
=
=
=
=
U
I
T
G
S
F
I
ext.
and
This feedback is carried out via the same succession of controls as with [§2.6.2] in
substituent
TEMP_EXT
with
COEF_H
in
AFFE_CHAR_THER_F
.
Note:
·
The essential data of this calculation, the field derived outside air temperature
I
ext.
ext.
I
T
T
I
=
, is
provided by
chth_v
.
·
To take into account a more sophisticated modeling of the outside air temperature
()
() ()
(
)
=
I
ext.
I
I
I
I
ext.
ext.
T
T
I
T
T
T
,
,
X
X
X
in these calculations of sensitivity, it is enough to substitute for the indicating function the new one
source term
I
I
I
ext.
ext.
I
T
T
=
in the definition of the derived loading
chth_v
.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
28/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
3
Nonlinear thermics
In nonlinear thermics the characteristics materials
C
p
(X, T)
and
(X, T)
can depend on
temperature (in non-linear thermics, one cannot define anisotropic materials. There does not exist
of modeling
THER_NL_ORTH
. Modeling
THER_NL)
. The body is subjected to the same types of
limiting conditions and of loadings which the linear problem to which are added two conditions
non-linear: imposed normal flow
I (X, T)
(on
4
) and radiation ad infinitum of a gray body (on
5
).
This last condition is modelized (
P
0
) by its emissivity
(X, T)
, the constant of Stefan-
Boltzmann
(X, T)
and the temperature ad infinitum
()
T
T
,
X
.
The operator of Code_Aster dedicated to this type of problem is
THER_NON_LINE
[R5.02.02]. It
allows to solve the problem in extreme cases interfered according to (Cauchy-Dirichlet-Neumann-Robin type
Inhomogenous, non-linear radiation and with variable coefficients)
()
()
(
)
] [
] [
()
] [
()
] [
()
()
] [
()
(
) (
)
[
]
] [












=
×
+
-
+
=
×
=
×
=
+
×
=
×
=
×
=
-
)
(
)
0
,
(
,
0
15
.
273
15
.
273
,
0
,
0
,
0
,
0
,
0
div
0
5
4
4
4
3
2
1
X
X
T
T
T
T
N
T
T
T
I
N
T
T
HT
HT
N
T
T
G
N
T
T
F
T
S
T
T
T
T
T
C
ext.
p
éq
3-1
Non-linearities pose theoretical problems to show the existence and the unicity of
the solution [bib2]. They can be also prejudicial with the numerical resolution itself.
Thus, with regard to the modeling of voluminal heat
C
p
(X, T)
, during an iteration,
either because the thermal transient is violent, or because the range of phase shift is
very small (for example, for a pure substance), two the reiterated successive ones of the temperature can
to locate on both sides one of its discontinuities. One then missed a large part of information
relating to the phase shift.
To free itself from this type of problem one rewrites the first equation of [éq 3-1] while introducing
a voluminal function enthalpy which will smooth these non-linearities (dependant on
T (X, T)
and while noting
T *
=
T (X, T *)
the value of the temperature at one moment
T * < T
arbitrary)
()
()
(
)
] [
()
()
=
×
=
-
T
T
p
D
C
T
S
T
T
T
T
*
,
0
div
with
éq
3-2
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
29/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Note:
·
Of share its definition (difference in a primitive (in temperature) of voluminal heat
between the temperature considered and a temperature
*
T
at one arbitrary moment), the function
enthalpy is known except for a constant of integration. It will thus have to be taken care that this
constant is eliminated in all the handled expressions.
·
The generalized conditions of type exchanges linear walls [éq 2-2] or relations between ddls
[éq 2-3] are also usable. As in [§2] one will not be interested in derived by
report/ratio with the parameters of the functionalities
LIAISON_ *
.
·
For a transitory calculation, as the problem of thermics linear, three strategies
can govern the choice of the field of initial temperature and they have an angle of attack on
the initialization of the derived problem.
·
Implicitly,
THER_NON_LINE
must tolerate conductivities rather badly strongly
non-linear. Because the tangent matrices and the initial predictive phase does not comprise it
term representing their derivative compared to the temperature. Moreover this term
complementary would have the bad taste to make nonsymmetrical the tangent matrix of
standard system and the matrix of the derived system! What is problematic to take in
count by means of computer in the handling of the structures of data.
In an any state of cause, with respect to the non-linear characteristics material
generally used, the perimeter of use of calculations of sensitivity is the same one
that that of the standard problem. It does not take into account rigorously
nonlinear thermal conductivity.
·
In addition, as in linear thermics, one has calculates only the sensitivity by report/ratio
with a constant parameter by geometrical area. What does not exclude a dependence
temporal, space or not-linear of characteristics material or loadings
nonconcerned by derivation. By parameterizing advisedly the loadings and
“derived” materials in the command file, one can also have access to
some derivatives made up (cf [§6.4]).
·
Into non-linear, obtaining a sensitivity by finished differences is even less reliable
that into linear, because it can be very sensitive to the degree of convergence of the solution. In
any rigor, that also influences the quality of the analytical sensitivity, the field of
temperature solution intervening in the assembly of the linear system “derived”.
·
Within the framework of non-linear thermics, the derivative compared to the enthalpy will not thus have
no the direction (a enthalpy is not constant!).
·
On the other hand, knowing the sensitivity
U
field of temperature compared to a parameter,
one reaches that of the enthalpy easily
v
compared to this same parameter, via the formula
()
U
C
T
C
T
T
T
v
p
p
=
=
=
=
:

3.1
Derived compared to voluminal heat
3.1.1 Elements
theoretical
That is to say
()
() (
)
=
I
I
I
I
p
I
C
X
X
, one will derive compared to the parameter
I
the formulation
[éq 3-2]. The required sensitivity is noted
I
T
U
=
. By noticing that the enthalpy can
to modelize like a function of time
T
and of
J
(considered makes some like functions
indicatrixes of the type of material)
(
)
()
()
(
)
(
)
*
,
,
,
*
,
,
,
:
T
T
C
D
C
T
p
T
T
T
J
p
J
J
J
-
=
X
X
X
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
30/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
its derivative is written
(
)




-
+
-
=




-
+
=
I
I
p
I
I
I
p
T
T
I
p
I
T
T
C
T
T
I
T
T
C
D
C
*
*
*
*
Indeed quantity
I
p
I
C
I
=
is independent of the temperature (and thus time) because one
suppose that the portions
I
are fixed (the under-parts of the body are supposed to be motionless, one
neglect in particular the phenomena of dilation). The derivative in time of the first term in enthalpy
of [éq 3-2] is thus worth
()
T
T
T
U
T
T
I
T
I
I


+
=




éq
3.1.1-1
In addition one has
(
)
(
)
{
{
{
{}
(
)
(
)
{
{
{
(
)
U
T
T
T
T
T
T
T
T
I
U
T
T
T
T
T
T
T
T
T
J
ij
I
J
J
I
I
I
J
ij
I
J
J
I
I
I
J
3
3
4
15
.
273
4
0
0
15
.
273
4
15
.
273
,
0
0
,
,
+
=








+
+
+
=
+
=








+
+
=
X
X
X
with
éq 3.1.1-2
The derivation of [éq 3-2] leads us to the new problem in extreme cases in
U
] [
] [
] [
] [
] [
(
)
] [














=
×
=
+




+
+
×
=
+




-
×
=
+


+
×
=
+
×
=
×
=


+
-


0
5
3
4
3
2
1
~
)
0
,
(
,
0
~
15
.
273
4
,
0
~
,
0
~
,
0
~
,
0
~
,
0
~
div
U
U
N
U
U
T
N
I
N
U
U
T
I
N
H
N
U
U
H
N
G
N
U
U
N
F
U
S
U
T
U
T
T
U
T
X
éq
3.1.1-3
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
31/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
with the new source voluminal and the new limiting conditions and initial
0
~
0
~
~
~
~
~
,
~
0
=
=
=
=
=
=
-
=
U
I
H
G
F
T
T
I
S
I
and
éq
3.1.1-4
Note:
·
The derivation of the formulation in voluminal heat [éq 3-1] led of course to same
result because voluminal heat is a function of
J
(considered makes some like
indicating functions of the type of material) and the temperature (depending itself on
time, of the variable of space and of
J
!)
(
)
(
)
J
J
p
T
T
C
,
,
,
X
Its derivative compared to
I
is thus written
{
{
{
{
I
T
T
T
T
T
T
C
C
C
I
J
ij
I
J
J
I
I
p
J
ij
I
J
J
p
I
p
+
=








+
+
+
=
0
0
X
from where
T
T
I
T
U
T
T
U
C
U
T
C
T
T
I
T
U
C
T
T
C
I
p
p
I
p
I
p
+


=
+
+
=
+
one finds well the formulation [éq 3.1.1-1] of the only term which distinguishes these two
modelings.

Contrary to linear thermics, the derived problem is completely metamorphosed.
The parabolic operator is modified and it became linear in
U
. All the problem is besides
become linear because the conditions limit underwent the same processing. The condition of Dirichlet
is from now on homogeneous and those of Neumann and radiation left room to conditions
of Robin. The condition of Cauchy became homogeneous. Contrary to the problem in
temperature, theoretical results of existence and unicity of the solution
U
are thus easier with
to exhume.
From a practical point of view, one does not need to have recourse, as in
THER_NON_LINE
, with one
algorithm of crossbred Newton-Raphson of a predictive phase to determine the increment of
temperature between two contiguous moments. A linear solvor is enough. No one is not need to assemble one
stamp tangent with each under-iteration. The resolution of the “derived” problem is thus more
rapid and more robust than that of the initial problem.
Semi-discretization in times of [éq 3.1.1-3] by
- method leads to the following problem:
to find a continuation
()
0
0
V
U
NR
N
N
such as (with the notations of [§2.1.1])
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
32/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
(
)
(
)


















=
-
=
+




+
+
-
=
+




-
-
=
+




+
-
=
+
-
=
-
-
=




+
-
-




+
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
0
0
5
1
1
1
1
3
1
1
1
1
4
1
1
1
1
1
1
3
1
1
1
1
1
1
2
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
1
~
)
0
(
1
0
~
15
.
273
4
1
0
~
1
0
~
1
0
~
1
0
~
1
0
~
~
div
1
div
U
U
NR
N
N
U
U
T
N
NR
N
I
N
U
U
T
I
N
NR
N
H
N
U
U
H
N
NR
N
G
N
U
U
N
NR
N
F
U
NR
N
T
S
S
U
T
U
T
U
T
U
T
T
U
T
U
T
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
éq 3.1.1-5
while posing
()
()
()
()


=


=


=


=
=
=
=
=
=
=
=
=
-
=


=


=
NR
N
T
T
NR
N
NR
N
NR
N
H
H
T
T
T
I
I
I
H
G
F
T
I
S
NR
N
T
T
NR
N
U
U
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
I
N
N
N
,
,
,
,
,
,
,
,
,
,
,
,
,
0
~
~
~
~
~
,
~
,
,
,
,
X
X
X
X
X
X
X
X
X
X
éq 3.1.1-6
By applying the theorem of Green to [éq 3.1.1-5], [éq 3.1.1-6] and by introducing the notations
following
(
)
{
}
(
)
{
}
I
NR
N
T
NR
N
T
NR
N
T
T
H
T
U
NR
N
NR
N
ext.
,
,
,
,
1
,
,
1
0
,
,
,
,
,
,
,
1
,






=






+
=
-


=


+
=
-
+
-
+
with
and
with
and
X
X
X
X
X
X
one is brought to solve the following variational problem
()
()




=
±
+
+
+
-
±
±
±
±
±
±
±
v
L
v
U
has
V
v
V
U
U
T
H
I
,
0
0
,
,
,
,
,
,
that
such
To calculate
and
given
Being
éq
3.1.1-7
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
33/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
with the bilinear form depending on the current moment
()
(
)
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
+




+
+
=
4
4
3
3
15
.
273
4
.
1
,
D
v
U
T
D
v
U
T
I
D
v
U
H
dx
v
U
U
T
T
dx
v
U
T
T
v
U
has
éq 3.1.1-8
and the linear form parameterized by the moments running and precedent
()
(
)
(
)
(
)
(
)
(
)
(
)
-
-
-
-
-
-
-
-
-
+
-
-
-
-
-
-
-
±
+
-
+
-
+
-
+
-
-




+
-
+
=
4
5
3
3
15
.
273
1
4
1
1
1
.
1
1
D
v
U
T
D
v
U
T
I
D
v
U
H
dx
v
T
T
I
T
dx
v
U
U
T
T
dx
v
U
T
T
v
L
I
éq 3.1.1-9
Note:
·
While posing in [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9]
=
=
=
=
=
±
±
±
5
4
,
0
,
and
T
C
T
p
one finds well the formulation [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] of the linear problem.
In addition, by deriving the variational formulation [éq 4.2-1] from the problem in temperature
[R5.02.02] one finds well [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9].
·
Contrary to the initial problem, the unknown field and the function test belong to
even functional space, which is more comfortable from a theoretical point of view and
numerical.
·
This time the two members of this equation are fundamentally different from those from
problem in temperature. However, after each pitch of time, once determined
T
+
with
to leave
T
-
, one does not have to reassemble all the matrix of the linear system and his second
associate member. It is enough to supplement the first tangent matrix (by taking them again
notations of [R5.02.02]) of the pitch of next time (allowing to pass from
+
+
+
=
1
T
T
with
+
+
2
T
)
by the term due to the non-linearity of thermal conductivity. One also leaves the second
member of the problem in temperature
(
)
±
-
T
T,
L
to constitute that which interests us: one
fabric by the terms of implicitation of non-linearities of the thermal conductivity and of
limiting conditions.
·
Concerning the initialization of the problem derived the remarks from the linear case apply in
extenso (cf [§2.2.2].
·
Contrary to calculation in temperature (at the time of the predictive phase of
THER_NON_LINE
,
the elimination of this constant imposes a suitable reformulation of the elementary term
CHAR_THER_EVOLNI
), the constant of integration of the enthalpy does not appear here because one
handle that its derivative in temperature.
·
For a condition of exchange between walls, the term of usual exchange is of course replaced
by (by taking again the notations of [éq 2-2])
(
)
±
±
±
-
ij
D
v
U
U
H
J
I
.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
34/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Space discretization of [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] provides the variational problem
discretized
(
)
()




=
±
+
+
+
±
±
-
±
±
±
±
±
H
H
H
H
H
H
H
H
H
H
H
H
H
H
H
H
H
v
L
v
U
has
V
v
V
U
I
U
T
H
,
0
0
,
,
,
,
,
,
that
such
To calculate
and
given
Being
éq
3.1.1-10
The taking into account of the conditions of Dirichlet leads then to the dualized linear system [éq 2.1.1-11].
We now will see how these calculations are declined in the code.
3.1.2 Establishment in Code_Aster
With the same notations that with [§2.1.1] one can break up the elementary matrix into
(
)
(
)
=
=
5
1
,
,
I
G
I
ij
G
ij
K
With
K
With
éq
3.1.2-1
with
·
the term of mass and thermal rigidity (option of calculation
MTAN_RIGI_MASS
by estimating them
characteristics material in T
-
)
(
)
() () ()
() ()
()
G
I
G
J
G
G
G
I
G
J
G
G
G
ij
NR
NR
NR
NR
T
T
K
With
+
=
+
+
.
,
1
·
a term of thermal rigidity due to the non-linearity of thermal conductivity (not introduced
because it is supposed, in the perimeter of use of the sensitivities, that
is independent of
T
)
(
)
()
() ()
()
G
I
G
J
G
G
G
G
ij
NR
NR
T
T
K
With
=
+
+
.
,
2
·
the term of rigidity due to the conditions limit of exchange (
MTAN_THER_COEF_F/R)
(
)
(
)
() ()
G
I
G
J
G
G
ij
NR
NR
K
H
K
With
3
3
3
,
=
+
In the event of exchange between walls this term is replaced by (
RIGI_THER_PARO_F/R
)
(
)
(
)
()
()
()
(
) ()
G
I
G
J
F
G
J
G
G
ij
NR
NR
NR
K
H
K
With
-
=
+
3
3
3
,
while noting
F
the bijection putting in opposite the two walls.
·
the term of rigidity due to the condition of non-linear Neumann (
MTAN_THER_FLUXNL
in
considering flow non-linear in T
-
)
(
)
() () ()
G
I
G
J
G
G
G
ij
NR
NR
T
I
K
With
-
=
+
,
4
4
·
the term of rigidity due to the condition of radiation (
MTAN_THER_RAYO_F/R)
(
)
() (
)
()
(
) () ()
G
I
G
J
G
G
G
ij
NR
NR
T
K
K
With
3
5
5
5
15
.
273
4
,
+
=
+
+
The second member is written, with the same notations,
(
)
(
)
=
=
6
1
,
,
I
G
I
J
G
J
K
L
K
L
éq
3.1.2-2
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
35/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
where
·
the term resulting from the implicitation of the matrix of rigidity and mass (new option
CHAR_SENS_EVOLNI
copy
CHAR_THER_EVOLNI
with
U
-
instead of
T
-
, fields
derived and standard material and the new source term)
(
)
() () ()
(
)
() ()
()
G
J
G
G
G
G
J
G
G
G
G
J
NR
U
NR
U
T
T
K
L
-
+
=
-
-
-
-
.
1
,
1
·
a term resulting from the implicitation of the non-linearity of thermal conductivity (not
introduced because it is supposed, in the perimeter of use of the sensitivities, that
is independent
of
T
)
(
)
(
)
() ()
()
()
G
I
G
G
G
G
G
J
NR
T
U
T
K
L
-
=
-
-
-
.
1
,
2
·
the term due to the “new source” (cf.
CHAR_SENS_EVOLNI)
(
)
()
() ()
(
) ()
G
J
G
G
I
G
G
J
NR
T
T
K
I
T
K
L
-
+
-
-
=
,
3
·
the term resulting from the implicitation of the conditions limit of exchange (
CHAR_THER_TEXT_F/R
with
T
ext.
=0
and
U
-
instead of
T
-
)
(
)
() (
)
() ()
G
J
G
G
G
J
NR
U
K
H
K
L
-
-
-
=
3
3
4
1
,
In the event of exchange between walls this term is replaced by (
CHAR_THER_PARO_F/R
with
U
-
with
place of
T
-
)
(
)
() (
)
()
()
(
)
(
) ()
G
J
G
G
G
G
J
NR
F
U
U
K
H
K
L
-
-
-
-
-
=
3
3
4
1
,
·
the term resulting from the implicitation of the condition of non-linear Neumann (new option
CHAR_SENS_FLUNL)
(
)
(
)
() () ()
G
J
G
G
G
G
J
NR
U
T
I
K
L
-
-
-
=
1
,
4
5
·
the term resulting from the implicitation of the condition of radiation (new option
CHAR_SENS_RAYO_F)
(
)
(
) () (
)
()
(
) () ()
G
J
G
G
G
G
J
NR
U
T
K
K
L
-
-
-
+
-
=
3
5
5
6
15
.
273
1
4
,
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
36/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Note:
·
The elementary terms of the matrix of the “derived” problem are identical to those of
stamp tangent initial problem. Indeed the problem in temperature can be written in
revealing a vector residue
R
[R5.02.02]
and a vectorial function test
V
(
)
(
)
±
-
+
+
=
T
T
T
T
,
.
,
.
T
L
V
T
R
V
While deriving compared to a parameter from calculation from the residue (for example a characteristic
material) noted
this vectorial relation and by translating it in indicielle form, one has
(
)
(
)
{
()
{{
(
)
ki
K
I
I
J
J
ij
ij
K
K
J
kj
J
K
K
K
K
K
K
V
L
L
U
T
With
K
L
V
T
K
T
R
V
L
V
R
V
=
=
=
=
+
+
+
+
+
±
+
+
±
+
with
4
3
42
1 T
One thus finds well the indicielle formulation of the derived linear system [éq 2.1.1-10].
·
Compared to the “derived” problem linear, the elementary matrix and the second member are
supplemented by the terms incorporating non-linearities of the thermal conductivity and of
limiting conditions. On the other hand, the new source is identical in both cases (one will be able
thus mutualiser the option of calculation).
·
To take into account a more sophisticated modeling of voluminal heat
()
() () (
)
=
I
I
I
I
I
p
T
I
T
C
,
,
X
X
X
in these calculations of sensitivity, it is enough to substitute for the indicating function the new term
source
I
I
I
p
I
C
=
in the definition of the field derived material
ma_v
.
As one already specified, all the elementary terms (except the second) of the matrix make
the object of an option of calculation and will already have been evaluated for the calculation of
+
+
2
T.
It thus remains with
to estimate the second member by re-using (with a different parameter setting) the options of calculation
existing or by introducing news (
CHAR_SENS_EVOLNI/FLUNL/RAYO_F
). These last
redirect towards the same routine of elementary calculation (
TE
.) that their associated options standards
(
CHAR_THER_EVOLNI/FLUNL/RAYO_F
).
The central character string (
SENS
instead of
THER
) joined to a detection of the nullity of the field
derived material, makes it possible to parameterize these routines towards one or the other their orientations
possible: calculation of an elementary term derived (new source term) or of a term “spectator”
had with the presence of a condition of exchange, radiation or a non-linear flow.
In accordance with the principles of architecture set up in the code to treat calculations
of sensitivity [bib6], the assembly and the resolution of [éq 2.1.1-11] are started by the analysis of
table of correspondence associated with the significant variable. This feedback is carried out via
even succession of controls that with [§2.1.2] in substituent
THER_NL
with
THER
in
DEFI_MATERIAU
and by replacing of course the operator
THER_LINEAIRE
by
THER_NON_LINE
.
We will unroll the same process for the various sensitivities, to start with that
concerning the other characteristic material.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
37/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
3.2
Derived compared to thermal conductivity
3.2.1 Elements
theoretical
By taking again the notations of [§2.2.1], one carries out derivation compared to the parameter
I
, with
I
T
U
=
sought sensitivity,
(
)
()
()
T
U
T
T
T
T
T
U
C
T
U
C
T
T
U
T
C
T
T
C
I
p
p
p
I
p


=


=
+
=


The derivation of [éq 3-2] delivers to us a problem in extreme cases identical to [éq 3.1.1-1] but with one
another voluminal source and of new conditions of Robin
(
)
N
T
I
I
H
G
T
I
S
I
I
-
=
=
=
=
=
~
~
~
~
,
div
~
éq
3.2.1-1
It is a homogeneous problem of Dirichlet-Cauchy and inhomogenous Robin. One can thus take back
same remarks concerning the linear character of the derived problem and simplifications
theoretical and numerical that that implies. Its semi-discretization in time results in seeking one
continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 3.1.1-5], [éq 3.1.1-6] of which the first relation
rewrites itself
(
)
(
)
1
0
~
1
~
div
1
div
1
1
1
1
1
1
1
1
-
-
+
=




+
-
-




+
-
-
+
+
+
+
+
+
+
+
NR
N
S
S
U
T
U
T
U
T
U
T
T
U
T
U
T
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
N
éq
3.2.1-2
with the new source term
(
)
{
}
1
,
div
~
+
=
N
N
m
T
I
S
m
I
m
éq
3.2.1-3
and new limiting conditions
N
T
I
I
H
G
N
I
N
N
N
N
-
=
=
=
=
+
+
+
+
+
1
1
1
1
1
~
~
~
~
éq
3.2.1-4
From where a variational problem identical to [éq 3.1.1-7] comprising the same bilinear form
[éq 3.1.1-8] joined to the linear form [éq 3.1.1-9] whose only third integral is modified for
to adapt to the new source
()
(
)
(
)
L
L
-
+
±
-
+
-
=
dx
v
T
T
I
v
L
I
1
éq
3.2.1-5
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
38/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Note:
·
While posing in the new variational problem [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.2.1-5]
=
=
=
=
=
±
±
±
5
4
,
0
,
and
T
C
T
p
one finds well the formulation [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] + [éq 2.2.1-5] of
linear problem. In addition, by deriving the variational formulation [éq 4.2-1] from
problem in temperature [R5.02.02] one finds well [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.2.1-5].
The new source term is identical into linear and non-linear.
·
For the moment, the operator
DEFI_MATERIAU
does not allow to modelize a conductivity
orthotropic non-linear thermics. There are not thus the particular cases of the linear problem.
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].

3.2.2 Establishment in Code_Aster
Compared to [§3.1.2], only the term due to the new source is to be modified. The same one is taken
(
)
G
J
K
L
,
3
that for the linear case (cf [éq 2.2.2-1]). What is made in the new option of calculation
CHAR_SENS_EVOLNI
with the fields derived and standard material.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. This feedback is carried out via the same one
succession of controls that with [§2.1.2] in substituent
THER_NL/LAMBDA
with
THER/RHO_CP
in
DEFI_MATERIAU
and by replacing of course the operator
THER_LINEAIRE
by
THER_NON_LINE
.
Note:
·
To take into account a more sophisticated modeling of thermal conductivity
()
() () (
)
=
I
I
I
I
I
T
I
T
,
,
X
X
X
in these calculations of sensitivity, it is enough to substitute for the indicating function the new term
source
I
I
I
I
=
in the definition of the field derived material
ma_v
.

3.3
Derived compared to the source
3.3.1 Elements
theoretical
By applying the formulas of the preceding paragraphs, the derivation of [éq 3-2] compared to
parameter
I
S
(cf [§2.3.1]) a problem in extreme cases identical to [éq 3.1.1-1] but with one delivers to us
another voluminal source
I
I
S
=
~
éq
3.3.1-1
One can thus renew the same remarks concerning the linear character of the derived problem and
theoretical and numerical simplifications that that implies. Its semi-discretization in time
conduit to seek a continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 3.1.1-5], [éq 3.1.1-6]
whose first relation is rewritten, while noting
I
S
T
U
=
sought sensitivity,
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
39/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
(
)
1
0
div
1
div
1
1
1
1
1
1
1
-
=




+
-
-




+
-
-
+
+
+
+
+
+
+
NR
N
I
U
T
U
T
U
T
U
T
T
U
T
U
T
I
N
N
N
N
N
N
N
N
N
N
N
N
N
N
éq 3.3.1-2
The application of the theorem of Green to [éq 3.3.1-2] led to solve a variational problem
identical to [éq 3.1.1-7] comprising the same bilinear form [éq 3.1.1-8] joined to the linear form
[éq 3.1.1-9] whose only third integral is modified to adapt to the new source
()
L
L
±
+
=
dx
v
I
v
L
I
éq
3.3.1-3
Note:
·
While posing in the new variational problem [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.3.1-3]
=
=
=
=
=
±
±
±
5
4
,
0
,
and
T
C
T
p
one finds well the formulation [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] + [éq 2.3.1-3] of
linear problem. In addition, by deriving the variational formulation [éq 4.2-1] from
problem in temperature [R5.02.02] one finds well [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.3.1-3].
·
The new source term is identical into linear and non-linear.
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].

3.3.2 Establishment in Code_Aster
Compared to [§3.1.2], only the term due to the new source is to be modified. The same one is taken
(
)
G
J
K
L
,
3
that for the linear case (cf [éq 2.3.2-1]) and the same option of calculation.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. This feedback is carried out via the same one
succession of controls that with [§2.3.2] by replacing the operator
THER_LINEAIRE
by
THER_NON_LINE
.
Note:
The taking into account of more sophisticated modeling of the source term is carried out as in
linear.

3.4
Derived compared to the imposed temperature
3.4.1 Elements
theoretical
By applying the formulas of the preceding paragraphs, the derivation of [éq 3-2] compared to
parameter
I
F
(cf [§2.4.1]) a problem in extreme cases identical to [éq 3.1.1-1] but with one delivers to us
null voluminal source and another condition of Dirichlet
I
I
F
S
=
=
~
0
~
and
éq
3.4.1-1
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
40/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
It is an inhomogenous problem of Dirichlet-Robin and homogeneous Cauchy. One can thus take back
same remarks concerning the linear character of the derived problem and simplifications
theoretical and numerical that that implies. Its semi-discretization in time results in seeking one
continuation
()
1
0
V
U
NR
N
N
checking a system similar to [éq 3.1.1-5], [éq 3.1.1-6] of which the first relation
rewrites itself, while noting
I
F
T
U
=
sought sensitivity,
(
)
1
0
0
div
1
div
1
1
1
1
1
1
1
-
=




+
-
-




+
-
-
+
+
+
+
+
+
+
NR
N
U
T
U
T
U
T
U
T
T
U
T
U
T
N
N
N
N
N
N
N
N
N
N
N
N
N
N
éq
3.4.1-2
with the new limiting condition
I
N
I
F
=
+
1
~
éq
3.4.1-3
The application of the theorem of Green to [éq 3.4.1-2], [éq 3.4.1-3] resulted in solving a problem
variational identical to [éq 3.1.1-7] comprising the same bilinear form [éq 3.1.1-8] joined to the form
linear [éq 3.1.1-9] whose only third integral is modified to adapt to the news
source: here this integral is null.
Note:
·
While posing in the new variational problem [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
“null source” + [éq 3.4.1-3]
=
=
=
=
=
±
±
±
5
4
,
0
,
and
T
C
T
p
one finds well the formulation [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] + [éq 2.4.1-2],
[éq 2.4.1-3], [éq 2.4.1-4] of the linear problem. In addition, by deriving the formulation
variational [éq 4.2-1] of the problem in temperature [R5.02.02] one finds [éq 3.1.1- well
7], [éq 3.1.1-8], [éq 3.1.1-9] + “null source” + [éq 3.4.1-3].
·
The new source term is identical into linear and non-linear.
·
One finds the same linear system not dualized as for the derivative in enthalpy. Only
taken into account of the condition of Dirichlet via of Lagranges will make the difference.
The space discretization and the taking into account of the inhomogenous condition of Dirichlet lead to
dualized linear system [éq 2.4.1-5]. Components of its second associate member with
Lagranges are nonnull.

3.4.2 Establishment in Code_Aster
Compared to [§3.1.2], only the term due to the new source is to be modified. The same one is taken
(
)
0
,
3
=
G
J
K
L
that for the linear case (cf [§2.4.2]) and the same option of calculation to assemble it
dualized linear system associated the “derived” conditions of Dirichlet.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. This feedback is carried out via the same one
succession of controls that with [§2.4.2] by replacing the operator
THER_LINEAIRE
by
THER_NON_LINE
.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
41/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Note:
·
Concerning various modelings of the imposed temperature and the conditions of
Dirichlet generalized, there is no difference between the non-linear one and the linear one. For more
information one can thus refer to [§2.4.2].
·
The taking into account of more sophisticated modeling of the imposed temperature is carried out
as into linear.
3.5
Derived compared to linear imposed normal flow
3.5.1 Elements
theoretical
By applying the formulas of the preceding paragraphs, the derivation of [éq 3-2] compared to
parameter
I
G
(cf [§2.5.1]) a problem in extreme cases identical to [éq 3.1.1-1] but with one delivers to us
null voluminal source and another condition of Robin on
2
I
I
S
G
=
=
~
0
~
and
éq
3.5.1-1
It is a homogeneous problem of Dirichlet-Cauchy and inhomogenous Robin. One can thus take back
same remarks concerning the linear character of the derived problem and simplifications
theoretical and numerical that that implies. Its semi-discretization in time results in seeking one
continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 3.1.1-5], [éq 3.1.1-6] of which the first relation
rewrites itself, while noting
I
G
T
U
=
sought sensitivity,
(
)
1
0
0
div
1
div
1
1
1
1
1
1
1
-
=




+
-
-




+
-
-
+
+
+
+
+
+
+
NR
N
U
T
U
T
U
T
U
T
T
U
T
U
T
N
N
N
N
N
N
N
N
N
N
N
N
N
N
éq
3.5.1-2
with the new limiting condition
I
N
I
G
=
+
1
~
éq
3.5.1-3
From where a variational problem identical to [éq 3.1.1-7] comprising the same bilinear form
[éq 3.1.1-8] joined to the linear form [éq 3.1.1-9] whose only third integral is modified for
to adapt to the new “surface” source
()
L
L
±
+
=
2
dx
v
I
v
L
I
éq
3.5.1-4
Note:
·
While posing in the new variational problem [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.5.1-4]
=
=
=
=
=
±
±
±
5
4
,
0
,
and
T
C
T
p
one finds well the formulation [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] + [éq 2.5.1-4] of
linear problem. In addition, by deriving the variational formulation [éq 4.2-1] from
problem in temperature [R5.02.02] one finds well [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.5.1-4].
·
The new source term is identical into linear and non-linear.
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
42/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
3.5.2 Establishment in Code_Aster
Compared to [§3.1.2], only the term due to the new source is to be modified. The same one is taken
(
)
G
J
K
L
,
2
3
that for the linear case (cf [§2.5.2]) and the same option of calculation.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. This feedback is carried out via the same one
succession of controls that with [§2.5.2] by replacing the operator
THER_LINEAIRE
by
THER_NON_LINE
.
Note:
·
Concerning various modelings of this condition of linear Neumann there is not
no difference between the linear and non-linear derived problem. For more information one
can thus refer to [§2.5.2].
·
The taking into account of more sophisticated modeling of imposed normal flow is carried out
as into linear.

3.6
Derived compared to non-linear imposed normal flow
3.6.1 Elements
theoretical
One poses
()
() (
)
=
I
I
I
I
I
I
I
I
X
X
and
I
I
T
U
=
required sensitivity. Portions of border
external
J
4
being solidified, one has
I
I
I
I
I
=
the indicating function of the ième portion
I
4
. By applying them
formulas of the preceding paragraphs, the derivation of [éq 3-2] compared to the parameter
I
I
us delivers
a problem in extreme cases identical to [éq 3.1.1-1] but with a null voluminal source and another
condition of Robin on
4
I
I
S
I
=
=
~
0
~
and
éq
3.6.1-1
It is a homogeneous problem of Dirichlet-Cauchy and inhomogenous Robin. One can thus take back
same remarks concerning the linear character of the derived problem and simplifications
theoretical and numerical that that implies. Its semi-discretization in time results in seeking one
continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 3.1.1-5], [éq 3.1.1-6] of which the first relation
rewrites itself
(
)
1
0
0
div
1
div
1
1
1
1
1
1
1
-
=




+
-
-




+
-
-
+
+
+
+
+
+
+
NR
N
U
T
U
T
U
T
U
T
T
U
T
U
T
N
N
N
N
N
N
N
N
N
N
N
N
N
N
éq
3.6.1-2
with the new limiting condition
I
N
I
I
=
+
1
~
éq
3.6.1-3
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
43/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
From where a variational problem identical to [éq 3.1.1-7] comprising the same bilinear form
[éq 3.1.1-8] joined to the linear form [éq 3.1.1-9] whose only third integral is modified for
to adapt to the new “surface” source
()
L
L
±
+
=
4
dx
v
I
v
L
I
éq
3.6.1-4
Note:
·
While posing in the new variational problem [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.6.1-4]
=
=
=
=
=
=
=
±
±
±
±
±
5
2
4
,
2
,
0
,
,
0
,
and
G
I
T
I
T
C
T
p
one finds well the formulation [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] + [éq 2.5.1-4] of
linear problem. In addition, by deriving the variational formulation [éq 4.2-1] from
problem in temperature [R5.02.02] one finds well [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.6.1-4].
·
This new source term is identical, except for the borders, at the end source of the derivative
compared to linear normal flow.
·
In stationary regime, it is of course unchanged.
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].
3.6.2 Establishment in Code_Aster
Compared to [§3.1.2], only the term due to the new source is to be modified. The same one is taken
(
)
G
J
K
L
,
4
3
that for the linear case by replacing the border
2
by
4
(cf [éq 2.5.2-1]). What
is made in the new option of calculation
CHAR_SENS_FLUNL
with the field derived non-linear flow.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. This feedback is carried out via the same one
succession of controls that with [§2.5.2] in substituent
FLUX_NL
with
FLUX_REP
in
AFFE_CHAR_THER_F
and by replacing of course the operator
THER_LINEAIRE
by
THER_NON_LINE
.
Note:
To take into account a truly non-linear modeling
()
() () (
)
=
I
I
I
I
I
I
T
I
I
T
I
,
,
X
X
X
in these calculations of sensitivity, it is enough to substitute for the indicating function the new term
source
I
I
I
I
I
I
=
in the definition of the derived loading
chth_v
.
3.7
Derived compared to the coefficient from convectif exchange
3.7.1 Elements
theoretical
By applying the formulas of the preceding paragraphs, the derivation of [éq 3-2] compared to
parameter
I
H
(cf [§2.6.1]) a problem in extreme cases identical to [éq 3.1.1-1] but with one delivers to us
null voluminal source and another condition of Robin on
3
(
)
T
T
I
S
ext.
I
H
-
=
=
~
0
~
and
éq
3.7.1-1
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
44/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
It is a homogeneous problem of Dirichlet-Cauchy and inhomogenous Robin. One can thus take back
same remarks concerning the linear character of the derived problem and simplifications
theoretical and numerical that that implies. Its semi-discretization in time results in seeking one
continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 3.1.1-5], [éq 3.1.1-6] of which the first relation
rewrites itself, while noting
I
H
T
U
=
sought sensitivity,
(
)
1
0
0
div
1
div
1
1
1
1
1
1
1
-
=




+
-
-




+
-
-
+
+
+
+
+
+
+
NR
N
U
T
U
T
U
T
U
T
T
U
T
U
T
N
N
N
N
N
N
N
N
N
N
N
N
N
N
éq
3.7.1-2
with the new limiting condition
(
)
1
1
1
~
+
+
+
-
=
N
N
ext.
I
N
T
T
I
H
éq
3.7.1-3
From where a variational problem identical to [éq 3.1.1-7] comprising the same bilinear form
[éq 3.1.1-8] joined to the linear form [éq 3.1.1-9] whose only third integral is modified for
to adapt to the new source
()
(
)
(
)
(
)
{
}
L
L
-
-
+
+
±
-
-
+
-
+
=
3
1
dx
v
T
T
T
T
I
v
L
ext.
ext.
I
éq
3.7.1-4
Note:
·
While posing in the new variational problem [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.7.1-4]
=
=
=
=
=
±
±
±
5
4
,
0
,
and
T
C
T
p
one finds well the formulation [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] + [éq 2.6.1-4] of
linear problem. In addition, by deriving the variational formulation [éq 4.2-1] from
problem in temperature [R5.02.02] one finds well [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.7.1-4].
·
The new source term is identical into linear and non-linear.
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].

3.7.2 Establishment in Code_Aster
Compared to [§3.1.2], only the term due to the new source is to be modified. The same one is taken
(
)
G
J
K
L
,
3
3
that for the linear case (cf [éq 2.6.2-1]) and the same option of calculation.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. This feedback is carried out via the same one
succession of controls that with [§2.6.2] by replacing the operator
THER_LINEAIRE
by
THER_NON_LINE
.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
45/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Note:
·
The taking into account of more sophisticated modeling of the coefficient of exchange is carried out
as into linear.
·
The condition of imposed normal flow non-linear (cf [§3.6]) makes it possible to modelize an exchange
convectif with a non-linear coefficient of exchange via one
DEFI_FONCTION
adequate
()
() ()
(
)
T
T
T
T
T
J
T
I
ext.
,
,
)
,
(
,
X
X
X
X
-
=
Same manner, knowing
I
I
T
U
=
one could have access “easily” to
I
J
T
W
=
. Indeed,
{
{








-
+
-
=




J
T
J
T
I
T
T
J
I
I
T
J
T
ext.
from where
(
)
(
)
ju
T
T
U
W
ext.
+
-
=
1

3.8
Derived compared to the outside air temperature
3.8.1 Elements
theoretical
By applying the formulas of the preceding paragraphs, the derivation of [éq 3-2] compared to
parameter
I
ext.
T
(cf [§2.7.1]) a problem in extreme cases identical to [éq 3.1.1-1] but with one delivers to us
null voluminal source and another condition of Robin on
3
I
hI
S
H
=
=
~
0
~
and
éq
3.8.1-1
It is a homogeneous problem of Dirichlet-Cauchy and inhomogenous Robin. One can thus take back
same remarks concerning the linear character of the derived problem and simplifications
theoretical and numerical that that implies. Its semi-discretization in time results in seeking one
continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 3.1.1-5], [éq 3.1.1-6] of which the first relation
rewrites itself, while noting
I
ext.
T
T
U
=
sought sensitivity,
(
)
1
0
0
div
1
div
1
1
1
1
1
1
1
-
=




+
-
-




+
-
-
+
+
+
+
+
+
+
NR
N
U
T
U
T
U
T
U
T
T
U
T
U
T
N
N
N
N
N
N
N
N
N
N
N
N
N
N
éq
3.8.1-2
with the new limiting condition
I
N
N
I
H
H
1
1
~
+
+
=
éq
3.8.1-3
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
46/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
From where a variational problem identical to [éq 3.1.1-7] comprising the same bilinear form
[éq 3.1.1-8] joined to the linear form [éq 3.1.1-9] whose only third integral is modified for
to adapt to the new “surface” source
()
(
)
{
}
L
L
-
+
±
-
+
+
=
3
1
dx
v
H
H
I
v
L
I
éq
3.8.1-4
Note:
·
While posing in the new variational problem [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.8.1-4]
=
=
=
=
=
±
±
±
5
4
,
0
,
and
T
C
T
p
one finds well the formulation [éq 2.1.1-6], [éq 2.1.1-7], [éq 2.1.1-8] + [éq 2.6.1-4] of
linear problem. In addition, by deriving the variational formulation [éq 4.2-1] from
problem in temperature [R5.02.02] one finds well [éq 3.1.1-7], [éq 3.1.1-8], [éq 3.1.1-9] +
[éq 3.8.1-4].
·
The new source term is identical into linear and non-linear.
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].

3.8.2 Establishment in Code_Aster
Compared to [§3.1.2], only the term due to the new source is to be modified. The same one is taken
(
)
G
J
K
L
,
3
3
that for the linear case (cf [éq 2.7.2-1]) and the same option of calculation.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. This feedback is carried out via the same one
succession of controls that with [§2.6.2] in substituent
TEMP_EXT
with
COEF_H
in
AFFE_CHAR_THER_F
and while replacing, of course, the operator
THER_LINEAIRE
by
THER_NON_LINE
.
Note:
The taking into account of a more sophisticated modeling of the outside air temperature
be carried out as into linear.

3.9
Derived compared to emissivity/constant of Stefan-Boltzmann
3.9.1 Elements
theoretical
One poses
()
() (
)
=
I
I
I
I
I
X
X
(resp.
()
() (
)
=
I
I
I
I
I
X
X
) the parameter and
I
T
U
=
(resp.
I
T
U
=
) required sensitivity. Portions
J
5
being solidified, one has
I
I
I
=
(resp.
I
I
I
=
) the indicating function of the ième portion
I
5
. By applying the formulas of
preceding paragraphs, the derivation of [éq 3-2] delivers a problem in extreme cases identical to us to
[éq 3.1.1-3] but with another voluminal source and another condition of Robin on
5
(resp. in
inverting the role of
and of
)
(
) (
)
{
}
4
4
15
.
273
15
.
273
~
0
~
+
-
+
=
=
T
T
I
I
S
I
and
éq
3.9.1-1
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
47/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
One can thus renew the same remarks concerning the linear character of the derived problem and
theoretical and numerical simplifications that that implies. Its semi-discretization in led time
to seek a continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 3.1.1-5], [éq 3.1.1-6] of which
first relation is rewritten
(
)
1
0
0
div
1
div
1
1
1
1
1
1
1
-
=




+
-
-




+
-
-
+
+
+
+
+
+
+
NR
N
U
T
U
T
U
T
U
T
T
U
T
U
T
N
N
N
N
N
N
N
N
N
N
N
N
N
N
éq
3.9.1-2
with the new condition (resp limits. by inverting the role of
and of
)
(
) (
)
{
}
4
1
4
1
1
1
15
.
273
15
.
273
~
+
-
+
=
+
+
+
+
N
N
N
I
N
T
T
I
I
éq
3.9.1-3
From where a variational problem identical to [éq 3.1.1-7] comprising the same bilinear form
[éq 3.1.1-8] joined to the linear form [éq 3.1.1-9] whose only third integral is modified for
to adapt to the new “surface” source (resp. by inverting the role of
and of
)
()
(
) (
)
{
}
(
)
(
) (
)
{
}
dx
T
T
v
I
T
T
v
I
v
L
I
I
4
4
4
4
15
.
273
15
.
273
1
15
.
273
15
.
273
5
+
-
+
-
-
+
-
+
+
=
-
-
-
+
+
+
±
L
éq 3.9.1-4
Note:
·
Derivation compared to the constant of Stefan-Boltzmann has certainly only one interest
practical miner. But its overcost of numerical establishment was moderate and it allows
cross validations.
·
In stationary regime, this source term is reduced to (resp. by inverting the role of
and of
)
()
(
) (
)
{
}
L
L
+
-
+
+
=
5
4
4
15
.
273
15
.
273
dx
v
T
T
I
v
L
I
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].

3.9.2 Establishment in Code_Aster
Compared to [§3.1.2], only the term due to the new source is to be modified (resp. by inverting it
role of
and of
)
(
)
() (
)
(
)
(
)
()
(
)
{
}
()
(
) () (
)
(
)
(
)
()
(
)
{
}
()
G
J
G
I
G
G
J
G
I
G
G
J
NR
T
K
T
K
I
NR
T
K
T
K
I
K
L
4
4
5
5
4
4
5
5
5
3
15
.
273
15
.
273
1
15
.
273
15
.
273
,
+
-
+
-
+
+
-
+
=
-
-
-
+
+
+
éq 3.9.2-1
What is made in the new option of calculation
CHAR_SENS_RAYO_F
with the fields of
standard and derived radiation.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. This feedback is carried out via
succession of controls.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
48/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
v = DEFI_VALEUR_SENSI (VALE = < value of
I
>)
(resp.
I
)
chth = AFFE_CHAR_THER_F (RAYONNEMENT= _F (GROUP_MA = < definition of
I
5
>,
SIGMA =
v, EPSILON = W, TEMP_EXT = Z))
one = DEFI_CONSTANTE (VALE = 1. )
zero = DEFI_CONSTANTE (VALE = 0. )
MEMO_NON_SENSI (NAME =_F
(NOM_SD = “chth”, PARA_SENSI = “v”,
NOM_COMPOSE = `
chth_v'))
chth_v = AFFE_CHAR_THER_F (RADIATION = _F (GROUP_MA = <
I
5
>,
SIGMA =
one, EPSILON = zero, TEMP_EXT = zero))
resu = THER_NON_LINE (EXCIT = chth,
SENSITIVITY
=
(
v)
…)
Note:
·
The essential data of this calculation, the field derived radiation
0
,
0
,
=
=
=
I
I
I
I
T
I
(resp.
I
I
I
=
), is provided by
chth_v
.
·
To take into account a more sophisticated modeling of emissivity (resp. constant
of Stefan-Boltzmann)
()
() () (
)
=
I
I
I
I
I
T
I
T
,
,
X
X
X
in these calculations of sensitivity, it is enough to substitute for the indicating function the new term
source
I
I
I
I
=
in the definition of the derived loading
chth_v
.

3.10 Derived compared to the temperature ad infinitum
3.10.1 Theoretical elements
One poses
()
()
(
)
=
I
I
I
I
T
I
T
T
X
X
the parameter and
I
T
T
U
=
required sensitivity.
portions
J
5
being solidified, one has
I
I
I
T
T
=
the indicating function of the ième portion
I
5
. By applying them
formulas of the preceding paragraphs, the derivation of [éq 3-2] delivers a problem in extreme cases to us
identical to [éq 3.1.1-3] but with another voluminal source and another condition of Robin on
5
(
)
3
15
.
273
4
~
0
~
+
=
=
T
I
I
S
I
and
éq
3.10.1-1
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
49/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
One can thus renew the same remarks concerning the linear character of the derived problem and
theoretical and numerical simplifications that that implies. Its semi-discretization in time
conduit to seek a continuation
()
0
0
V
U
NR
N
N
checking a system similar to [éq 3.1.1-5], [éq 3.1.1-6]
whose first relation is rewritten
(
)
1
0
0
div
1
div
1
1
1
1
1
1
1
-
=




+
-
-




+
-
-
+
+
+
+
+
+
+
NR
N
U
T
U
T
U
T
U
T
T
U
T
U
T
N
N
N
N
N
N
N
N
N
N
N
N
N
N
éq
3.10.1-2
with the new limiting condition
(
)
3
1
15
.
273
4
~
+
=
+
T
I
I
I
N
éq
3.10.1-3
From where a variational problem identical to [éq 3.1.1-7] comprising the same bilinear form
[éq 3.1.1-8] joined to the linear form [éq 3.1.1-9] whose only third integral is modified for
to adapt to the new “surface” source
()
()
(
)
(
) ()
(
)
{
}
L
L
-
-
+
+
±
+
-
+
+
+
=
5
3
3
15
.
273
1
15
.
273
4
dx
v
T
T
I
v
L
I
éq 3.10.1-4
Note:
In stationary regime, this source term is unchanged.
The space discretization and the taking into account of the homogeneous condition of Dirichlet lead to
dualized linear system [éq 2.1.1-11].
3.10.2 Establishment in Code_Aster
Compared to [§3.1.2], only the term due to the new source is to be modified
(
)
(
)
(
) (
)
(
) ()
(
)
(
)
(
) (
)
(
) ()
G
J
I
G
G
J
I
G
G
J
NR
K
T
K
I
NR
K
T
K
I
K
L
3
5
5
3
5
5
5
3
15
.
273
1
4
15
.
273
4
,
+
-
+
+
=
-
-
-
+
+
+
éq
3.10.2-1
What is made in the new option of calculation
CHAR_SENS_RAYO_F
with the fields of
standard and derived radiation.
The assembly and the resolution of the linear system are started by the analysis of the table of
correspondence associated with the significant variable. This feedback is carried out via the same one
succession of controls that with [§3.9.2] in substituent
TEMP_EXT
with
SIGMA
.
Note:
Ad infinitum to take into account a more sophisticated modeling of the temperature
()
() ()
(
)
=
I
I
I
I
I
T
T
I
T
T
T
,
,
X
X
X
in these calculations of sensitivity, it is enough to substitute for the indicating function the new term
source
I
I
I
I
T
T
=
in the definition of the derived loading
chth_v
.
We now will recapitulate all the linear systems “derived” to assemble according to
desired sensitivities.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
50/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
4
Summary of the sensitivities of the temperature
That one is in or not linear linear thermics, the sensitivity required to the current moment
U
+
is
solution of a dualized system of the type (with the notations of [§2.1])




=








-
-
=
+
+
C
C
L
U
Id
Id
B
Id
Id
B
B
B
With
L
U
With
T
T
~
~
~
where
(
)
(
)
=
I
G
I
ij
G
ij
K
With
K
With
,
,
(
)
(
)
=
I
G
I
J
G
J
K
L
K
L
,
,
One will draw up the nomenclature of the potential elementary terms of this system and of their
option of calculation in Code_Aster. In order to be more synthetic, one agrees here on a classification
different from those practiced in the preceding paragraphs.
For the matrix, there are the seven following possible terms:
·
MASS_THER
(cf [§2.1.2])
(
)
()
() ()
G
I
G
J
p
G
G
ij
NR
NR
K
C
T
K
With
=
,
1
·
RIGI_THER
(cf [§2.1.2])
(
)
()
()
()
G
I
G
J
G
G
ij
NR
NR
K
K
With
=
.
,
2
·
RIGI_THER_COEF_F/R
or
MTAN_THER_COEF_F/R
(cf [§2.1.2], [§3.1.2])
(
)
(
)
() ()
G
I
G
J
G
G
ij
NR
NR
K
H
K
With
3
3
3
,
=
+
In the event of exchange between walls this term is replaced by
RIGI_THER_PARO_F/R
or
MTAN_THER_PARO_F/R
(cf [§2.1.2], [§3.1.2])
(
)
(
)
()
()
()
(
) ()
G
I
G
J
F
G
J
G
G
ij
NR
NR
NR
K
H
K
With
-
=
+
3
3
3
,
·
MTAN_RIGI_MASS
(cf [§3.1.2])
(
)
()
(
) () ()
()
(
) ()
()
G
I
G
J
G
G
G
I
G
J
G
G
G
ij
NR
NR
T
NR
NR
T
T
T
K
With
+
=
+
+
.
,
4
·
Not coded. (cf [§3.1.2])
(
)
()
(
)
() ()
()
G
I
G
J
G
G
G
G
ij
NR
NR
T
T
T
K
With
=
+
+
.
,
5
·
MTAN_THER_FLUXNL
(cf [§3.1.2])
(
)
()
(
) () ()
G
I
G
J
G
G
G
ij
NR
NR
T
T
I
K
With
+
-
=
,
4
6
·
MTAN_THER_RAYO_F/R
(cf [§3.1.2])
(
)
() (
)
()
(
) () ()
G
I
G
J
G
G
G
ij
NR
NR
T
K
K
With
3
5
5
7
15
.
273
4
,
+
=
+
+
For the second member, there are sixteen possibilities:
·
CHAR_SENS_EVOL
(cf [§2.1.2])
(
)
()
() ()
(
)
()
()
()
G
J
G
G
G
J
G
p
G
G
J
NR
U
K
NR
U
K
C
T
K
L
-
+
=
-
-
.
1
,
1
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
51/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
·
CHAR_THER_TEXT_F/R
(cf [§2.1.2])
(
)
() (
)
() ()
G
J
G
G
G
J
NR
U
K
H
K
L
-
-
-
=
3
3
2
1
,
In the event of exchange between walls this term is replaced by
CHAR_THER_PARO_F/R
(cf.
[§2.1.2])
(
)
() (
)
()
()
(
)
(
) ()
G
J
G
G
G
G
J
NR
F
U
U
K
H
K
L
-
-
-
-
-
=
3
3
2
1
,
·
CHAR_SENS_EVOL
(cf [§2.1.2])
(
)
()
() ()
(
) ()
G
J
G
G
I
G
G
J
NR
T
T
K
I
T
K
L
-
+
-
-
=
,
3
·
CHAR_SENS_EVOL
(cf [§2.1.2])
(
)
()
()
(
)
()
(
)
()
G
J
G
G
I
G
G
J
NR
T
T
K
I
K
L
-
+
-
=
-
+
1
,
4
·
CHAR_THER_SOUR_F/R
(cf [§2.3.2], [§3.3.2])
(
)
()
()
G
J
I
G
G
J
NR
K
I
K
L
=
,
5
·
CHAR_THER_FLUN_F/R
(cf [§2.5.2], [§3.5.2])
(
)
(
)
()
G
J
I
G
G
J
NR
K
I
K
L
2
2
6
,
=
·
CHAR_SENS_TEXT_F
(cf [§2.6.2], [§3.7.2])
(
)
(
)
(
)
()
(
)
(
) (
)
()
(
) ()
G
J
G
ext.
G
ext.
I
G
G
J
NR
T
K
T
T
K
T
K
I
K
L






-
-
+
-
=
-
-
+
+
3
3
3
3
7
1
,
·
CHAR_THER_TEXT_F
(cf [§2.7.2], [§3.8.2])
(
)
(
)
(
) (
) (
)
{
} ()
G
J
I
G
G
J
NR
K
H
K
H
K
I
K
L
3
3
3
3
8
1
,
-
+
=
-
+
·
CHAR_THER_TEXT_F
(cf [§2.7.2], [§3.8.2])
(
)
()
(
) () ()
(
)
()
(
) ()
()
G
J
G
G
G
G
J
G
G
G
G
J
NR
U
T
NR
U
T
T
T
K
L
-
+
=
-
-
-
-
.
1
,
9
·
CHAR_SENS_EVOLNI
(cf [§3.1.2])
(
)
(
)
()
(
) () ()
()
G
I
G
G
G
G
G
J
NR
T
U
T
T
K
L
-
=
-
-
-
.
1
,
10
·
CHAR_SENS_FLUNL
(cf [§3.1.2])
(
)
(
)
()
(
) () ()
G
J
G
G
G
G
J
NR
U
T
T
I
K
L
-
-
-
=
1
,
4
11
·
CHAR_SENS_RAYO_F
(cf [§3.1.2])
(
)
(
) () (
)
()
(
) () ()
G
J
G
G
G
G
J
NR
U
T
K
K
L
-
-
-
+
-
=
3
5
5
12
15
.
273
1
4
,
·
CHAR_SENS_FLUNL
(cf [§3.6.2])
(
)
(
)
()
G
J
I
G
G
J
NR
K
I
K
L
4
4
13
,
=
·
CHAR_SENS_RAYO
(cf [§3.9.2])
(
)
() (
)
(
)
(
)
()
(
)
{
}
()
(
) () (
)
(
)
(
)
()
(
)
{
}
()
G
J
G
I
G
G
J
G
I
G
G
J
NR
T
K
T
K
I
NR
T
K
T
K
I
K
L
4
4
5
5
4
4
5
5
5
14
15
.
273
15
.
273
1
15
.
273
15
.
273
,
+
-
+
-
+
+
-
+
=
-
-
-
+
+
+
(
)
() (
)
(
)
(
)
()
(
)
{
}
()
(
) () (
)
(
)
(
)
()
(
)
{
}
()
G
J
G
I
G
G
J
G
I
G
G
J
NR
T
K
T
K
I
NR
T
K
T
K
I
K
L
4
4
5
5
4
4
5
5
5
15
15
.
273
15
.
273
1
15
.
273
15
.
273
,
+
-
+
-
+
+
-
+
=
-
-
-
+
+
+
·
CHAR_SENS_RAYO
(cf [§3.10.2])
(
)
(
)
(
) (
)
(
) ()
(
)
(
)
(
) (
)
(
) ()
G
J
I
G
G
J
I
G
G
J
NR
K
T
K
I
NR
K
T
K
I
K
L
3
5
5
3
5
5
5
16
15
.
273
1
4
15
.
273
4
,
+
-
+
+
=
-
-
-
+
+
+
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
52/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
The operator
B
allowing to check the condition of Dirichlet is written
(
)
J
I
U
I
I
=
+
+
with
U
B
with its second associate member
0
C
=
Note:
·
The possible taking into account of limiting conditions of Dirichlet generalized requires
to rewrite this operator in the form
J
J
U
J
J
J
=
+
+
with
U
B
·
As one rigorously does not take into account nonlinear thermal conductivity,
terms
With
5
and L
10
are not programmed yet.
The eighteen dualized systems exhumed in the preceding paragraphs can then gather
in the following table:
Type of
sensitivity
Variable
Stamp
Second
member
Key word
Characteristics
Thermics
linear
Stamp
identical to
that of
problem
direct
THER_LINEAIRE
Sensitivity compared to
a constant parameter
by area
(possibly derived
composed cf [§6.4])
Heat
voluminal
I
T
U
=
With
1
+
With
2
+
With
3
L
1
+
L
2
+
L
3
THER/RHO_CP
Conductivity
thermics
I
T
U
=
Idem
L
1
+
L
2
+
L
4
THER/LAMBDA
Sensitivity compared to
a characteristic
orthotropic not
accessible
Source
I
S
T
U
=
Idem
L
1
+
L
2
+
L
5
SOURCE
Three types of
modeling
Temperature
imposed
I
F
T
U
=
Idem
L
1
+
L
2
TEMP_IMPO
ij
I
J
K
F
F
C
=
=
Three types of
modeling.
Sensitivity compared to
a coefficient
multiplier of one
condition of Dirichlet
generalized not
accessible
Normal flow
imposed
I
G
T
U
=
Idem
L
1
+
L
2
+
L
6
FLUX_REP
Three types of
modeling.
Sensitivity compared to
a component of flow
vectorial accessible.
Exchange
convectif
I
H
T
U
=
Idem
L
1
+
L
2
+
L
7
EXCHANGE/COEF_H
or
ECHANGE_PAROI
Not
ECHANGE_PAROI
in lumpé.
Temperature
external
I
ext.
T
T
U
=
Idem
L
1
+
L
2
+
L
8
EXCHANGE/
TEMP_EXT

background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
53/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Type of
sensitivity
Variable
Stamp
Second
member
Key word
Characteristics
Thermics
nonlinear
Stamp
near to
stamp
tangent of
problem
direct.
The new ones
sources terms
are similar
into linear and in
nonlinear.
THER_NON_LINE
Sensitivity compared to
a constant parameter
by area
(possibly derived
composed cf [§6.4]).
No sensitivity by
report/ratio with the enthalpy.
Account is not held
of a conductivity not
linear
With
5
and L
10
not taken in
count
Heat
voluminal
I
T
U
=
With
3
+
With
4
+
With
5
+
With
6
+
With
7
L
2
+
L
3
+
L
9
+L
10
+
L
11
+
L
12
THER_NL/RHO_CP
Conductivity
thermics
I
T
U
=
Idem
L
2
+
L
4
+
L
9
+L
10
+
L
11
+
L
12
THER_NL/LAMBDA
Idem linear case
Source
I
S
T
U
=
Idem
L
2
+
L
5
+
L
9
+L
10
+
L
11
+
L
12
SOURCE
Idem linear case
Temperature
imposed
I
F
T
U
=
Idem
L
2
+
L
9
+L
10
+
L
11
+
L
12
TEMP_IMPO
Idem linear case
Normal flow
imposed
linear
I
G
T
U
=
Idem
L
2
+
L
6
+
L
9
+L
10
+
L
11
+
L
12
FLUX_REP
Idem linear case.
Normal flow
imposed not
linear
I
I
T
U
=
Idem
L
2
+
L
13
+
L
9
+L
10
+
L
11
+
L
12
FLUX_NL
Allows to modelize
sensitivity compared to
a convectif exchange
nonlinear
Exchange
convectif
linear
I
H
T
U
=
Idem
L
2
+
L
7
+
L
9
+L
10
+
L
11
+
L
12
EXCHANGE/COEF_H
or
ECHANGE_PAROI
Not
ECHANGE_PAROI
in lumpé.
Temperature
external
I
ext.
T
T
U
=
Idem
L
2
+
L
8
+
L
9
+L
10
+
L
11
+
L
12
EXCHANGE/
TEMP_EXT
Emissivity
I
T
U
=
Idem
L
2
+
L
15
+
L
9
+L
10
+
L
11
+
L
12
RADIATION/
EPSILON
Constant of
Stefan-
Boltzmann
I
T
U
=
Idem
L
2
+
L
14
+
L
9
+L
10
+
L
11
+
L
12
RADIATION/
SIGMA
Temperature with
the infinite one
I
T
T
U
=
Idem
L
2
+
L
16
+
L
9
+L
10
+
L
11
+
L
12
RADIATION/
TEMP_EXT
Table 4-1: Summary of the linear systems “derived”
We will now close the theoretical part this document by calculating the sensitivities of one
related size with the field of temperature: heat transfer rate.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
54/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
5
Sensitivity of the heat transfer rate
After having carried out a thermal calculation, one has access to the heat transfer rate (calculation is carried out either with
nodes (field with the nodes by elements via
OPTION=' FLUX_ELNO_TEMP'
), that is to say at the points of Gauss
(field at the points of Gauss by elements via
OPTION=' FLUX_ELGA_TEMP'
)) via the operators
CALC_ELEM/CALC_NO
. It is given starting from the field of temperature by using the Fourier analysis
who is written in linear thermics with isotropic materials
()
() ()
T
T
T
,
,
X
X
X
Q
-
=
With anisotropic materials, thermal conductivity is modelized by a diagonal matrix
(expressed in its reference mark of orthotropism). Into nonlinear (inevitably isotropic in Code_Aster) one
has on the other hand
()
() ()
T
T
T
T
,
,
,
X
X
X
Q
-
=
By taking again the notations and the step developed in the preceding paragraphs, one
synthesize the calculation of derived from the heat transfer rate in the following table. The column entitled
“formula” expresses the relation to be set up in the operators
CALC_ELEM/CALC_NO
for
to determine the sensitivity of the heat transfer rate required. It depends on the field of temperature and
of its sensitivity compared to the same parameter, all the two resulting one from a thermal calculation (via
THER_LINEAIRE
or
THER_NON_LINE
).
Note:
In practice, calculations of sensitivity of the heat transfer rate do not take account of the possible one
non-linearity of thermal conductivity. In non-linear thermics, the first term is not
thus not programmed.
Type of
sensitivity
Variables
of exit of
the operator
thermics
Sensitivity
sought
Formulate/
Options of calculation in Code_Aster
Thermics
linear
Heat
voluminal
T
T
U
I
,
=
I
=
Q
v
()
() ()
T
U
T
,
,
X
X
X
v
-
=
FLUX_ELGA/NO_TEMP
Conductivity
thermics
T
T
U
I
,
=
I
=
Q
v
()
() () () ()
T
U
T
T
I
T
I
,
,
,
X
X
X
X
X
v
-
-
=
FLUX_ELGA/NO_SENS
Source
T
S
T
U
I
,
=
I
S
=
Q
v
()
() ()
T
U
T
,
,
X
X
X
v
-
=
FLUX_ELGA/NO_TEMP
Temperature
imposed
T
F
T
U
I
,
=
I
F
=
Q
v
Idem
Normal flow
imposed
T
G
T
U
I
,
=
I
G
=
Q
v
Idem
Exchange
convectif
T
H
T
U
I
,
=
I
H
=
Q
v
Idem
Temperature
external
T
T
T
U
I
ext.
,
=
I
ext.
T
=
Q
v
Idem



background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
55/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Type of
sensitivity
Variables
of exit of
the operator
thermics
Sensitivity
sought
Formulate/
Options of calculation in Code_Aster
Thermics
nonlinear
Heat
voluminal
T
T
U
I
,
=
I
=
Q
v
()
() (
) ()
() ()
T
U
T
T
T
U
T
T
T
,
,
,
,
,
X
X
X
X
X
v
-
-
=
FLUX_ELGA/NO_TEMP
Conductivity
thermics
T
T
U
I
,
=
I
=
Q
v
()
() () ()
()
() ()
T
U
T
T
T
T
I
T
U
T
T
T
I
,
,
,
,
,
,
,
X
X
X
X
X
X
X
v
-


+
-
=
FLUX_ELGA/NO_SENS
Source
T
S
T
U
I
,
=
I
S
=
Q
v
()
() (
) ()
() ()
T
U
T
T
T
U
T
T
T
,
,
,
,
,
X
X
X
X
X
v
-
-
=
FLUX_ELGA/NO_TEMP
Temperature
imposed
T
F
T
U
I
,
=
I
F
=
Q
v
Idem
Normal flow
imposed
linear
T
G
T
U
I
,
=
I
G
=
Q
v
Idem
Normal flow
imposed not
linear
T
I
T
U
I
,
=
I
I
=
Q
v
Idem
Exchange
convectif
linear
T
H
T
U
I
,
=
I
H
=
Q
v
Idem
Temperature
external
T
T
T
U
I
ext.
,
=
I
ext.
T
=
Q
v
Idem
Emissivity
T
T
U
I
,
=
I
=
Q
v
Idem
Constant
of Stefan-
Boltzmann
T
T
U
I
,
=
I
=
Q
v
Idem
Temperature
ad infinitum
T
T
T
U
I
,
=
I
T
=
Q
v
Idem
Table 5-1 Sensitivities of the heat transfer rate
Note:
While posing in the formulas of nonlinear thermics
()
() ()
X
X
X
=
=
T
T
T
,
0
,
and
one finds the formulations of linear thermics well.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
56/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
6
Implementation in Code_Aster
6.1 Difficulties
particular
The main difficulty of these calculations of sensitivity is to detect the presence of parameters
sensitive in the standard loadings and materials, and, to associate the fields to them
suitable derivatives. To be done, a whole architecture (cf [bib6] and [§6.2], [§6.4]) has being put in
place in order to notify to the supervisor the affiliation of a derived field with such or such significant variable. Via
controls
DEFI_PARA_SENSI
and
MEMO_NOM_SENSI
and their associated utilities FORTRAN, one
can thus make the joint within a thermal operator, between a significant parameter, the field where
it intervenes and the associated derived field.
With a particular mention for the materials which are known only in one coded form. One
detect the significant characteristic of the aforesaid material by testing the nullity of derived material (a problem
related is posed for the characteristics of the radiation).
Inter alia installations, it A was also necessary:
·
To insert, within the thermal operators, a loop on the sensitivities requested,
wedged between the temporal loops and those on the loadings.
·
To set up the resolution of a linear system in the non-linear process of
THER_NON_LINE
.
·
To organize and manage the mutualisation of the matrices between the standard problem and it (or them)
problem (S) derived (S).
·
To take into account, in robust and fast manner, the insertion of possible parameters
insensitive (creation of one
CHAM_NO
null components for the thermal operators and
of one
CHAM_ELEM
null components for postprocessings of calculation of flow).
·
To avoid an inopportune proliferation of option of calculations (and theirs
TE000
FORTRAN
associated) while mutualisant and “bluffant” existing it (one can visualize the exact routing of
calculation (at the macroscopic level of the options) and easy ways deployed to re-use
existing the cf [§6.2]).
Beyond these fastidious developments, a large effort of validation “numérico-data processing”
was deployed on all the meshs supports, all modelings, all the loadings, all them
types of initialization of the thermal solveurs and for all the sensitivities. These hard tests on
small cases model tests (
TPLL01A/H
for
2D
PLAN
and
3D
and
TPNA01A
for
2D
AXIS
) are
revealed profitable (including for standard thermal calculation in lumpé!) and essential.
Because one seldom has theoretical values allowing to validate a thermal calculation
complicated: “nothing resembles more one sensitivity… than another value of sensitivity! ”. In
crossing with converged finished differences, one thus tried to release a maximum of
confidence in all building blocs constituting the derived system.
To be exhaustive on the aspect validation, let us note that several case-tests were delivered
(
SENST
<00> [V1.01.15 *]) including one analytical (
SENST04
cf [§6.4]).

6.2 Environments
necessary/parameter settings





Appear 6.2-a Désignation of thermal conductivities for a Bi-material
C
p
=
v=5
1
C
p
=10
2
v
T
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
57/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
To calculate the sensitivity of the thermal field compared to a constant parameter by area (here one
particular value of conductivity in one of the areas of a Bi-material during a stationary calculation),
it is first of all necessary to notify this significant parameter in the command file via one
DEFI_PARA_SENSI
v = DEFI_PARA_SENSI (VALE = 5. )
as well as the constant functions equal to zeros and the unit which will be used to define the essential one
indicating function
zero = DEFI_CONSTANTE (VALE = 0. )
one = DEFI_CONSTANTE (VALE = 1. )
Then, as for a standard problem, one definite the fields material (or the conditions limit)
associated. One does not have to modify the other standard loadings.
MA1 = DEFI_MATERIAU (THER = _F
(RHO_CP = v))
MA2 = DEFI_MATERIAU (THER = _F
(RHO_CP = 10. ))
affe = AFFE_MATERIAU (AFFE = (_F
(GROUP_MA = <
1
>, MATER =
ma1),
_F
(GROUP_MA = <
2
>, MATER =
ma2))
….
other insensitive data material and limiting conditions with respect to
v
It is then necessary to define the “derived” fields material (or the conditions limit “derived”) and
only those concerned with this derivation. It is them which will provide information
essential with the assembly of the derived problem: the indicating function of the derived field


=
=
2
1
1
on
0
on
1
v
p
C
I
.
One notifies to the supervisor their affiliation with the initial field via the control
MEMO_NOM_SENSI
(it is it
who filled a structure of data suitable sensitivity allowing to make the joint, within
thermal operators, between the significant parameter, the field where it intervenes and the derived field
associated)
MEMO_NON_SENSI (NAME =_F
(NOM_SD = “my”,
PARA_SENSI = `
v',
NOM_COMPOSE = `
ma_v'))
ma_v = DEFI_MATERIAU (THER = _F
(RHO_CP = one))
mazero = DEFI_MATERIAU (THER = _F
(RHO_CP = zero))
MEMO_NON_SENSI (NAME = _F
(NOM_SD = “affe”,
PARA_SENSI = `
v',
NOM_COMPOSE
= `
affe_v'
))
affe_v = AFFE_MATERIAU (AFFE = (_F
(GROUP_MA = <
1
>, MATER =
ma_v),
_F
(GROUP_MA = <
2
>, MATER =
mazero))
One concludes by requesting from the thermal operator the aforementioned calculation via the key word sensitivity and notifying to
supervisor the presence of a field sensitivity in
EVOL_THER
resu
(accessible in the field
symbolic system
resu_v
SD result)
MEMO_NON_SENSI (NAME =_F
(NOM_SD = “resu”,
PARA_SENSI = `
v',
NOM_COMPOSE = `
resu_v'))
resu = THER_LINEAIRE (CHAM_MATER = affe,
….
SENSITIVITY = (
v))
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
58/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
One can of course make same postprocessings on the derived field as on the standard field, it
is enough to add with the suitable controls (
CALC_ELEM
for the calculation of the flow of it
CHAM_ELEM
with
nodes,
CALC_NO
for its transformation into one
CHAM_NO
,
POST_RELEVE
for the statement of value,
IMPR_RESU
for the impression and
TEST_RESU
for the comparison of values) the key word
SENSITIVITY
and to inform it good value of parameter of derivation. The operator
the good field sensitivity sharing will select then the same SD as the standard field. By
example, one filled out here the table
rlresuv
certain components of derived from
T
compared to
v
rlresuv = POST_RELEVE_T (ACTION=_F (RESULTAT= resu,
….
SENSITIVITY = (
v))
By printing the table one can then visualize the components of the sensitivity specifically
wished.
--------------------------------------------------------------------------
ASTER 6.02.24
CONCEPT RLRESUV CALCULATES THE 27/03/2002 A 14:56:43 OF TYPE
TABL_POST_RELE
ENTITLE NODE RESU NOM_CHAM
PAR_SENS NUME_ORDRE INST COOR_X/Y/Z TEMP
TEMPERAT NO1
RLRESUV TEMP v 0 0. 0. 1.1 - 1.5 0. 1.69334E+01
TEMPERAT NO2
RLRESUV TEMP v 0 0. 0. 1.1 - 1.0 0. 1.88532E+01
TEMPERAT NO3
RLRESUV TEMP v 0 0. 0. 1.1 - 0.5 0. 1.19433E+01
TEMPERAT NO4
RLRESUV TEMP v 0 0. 0. 1.1 - 0.0 0. 2.15978E+01
Example 6.2-1: Layout, via
IMPR_TABLE
, in the file result
For further information, the reader will be able to consult [bib6], [U4.50.02] or the case-tests of
type
SENST
<00> [V1.01.15 *].
Note:
·
By parameterizing advisedly the loadings and materials “derived” in the file from
order, one can also have access to some derivatives made up. Cf the remark of
[§6.3].
·
The syntax of the specific controls of the sensitivity (
DEFI_PARA_SENSI
,
MEMO_NOM_SENSI
) can be brought to change according to the architectural choices which will be taken
(automation of the process of detection and parameter setting of the fields derived by
supervisor or not). Methodology will not be on the other hand modified and one will always be able
to pass in “manual”. The user then remaining only judge of the relevance of these calculations
“except limits”.
·
During the construction of the derived fields, in
AFFE_MATERIAU
or
AFFE_CHAR_THER
, one
need does not have to specify the “nullity” of the insensitive fields. They are initialized by defect
to zero.
To be completely exhaustive on the parameter setting, let us conclude by a functionality which can
to interest of future mainteneurs/developers or the “pointed” users. Unfolding
exact of calculation (options calculated with their field
IN
and
OUT
, loadings and materials taken in
count…) as well as the easy ways deployed to re-use existing it are traced in the file
message if one notifies
INFO=2
in the thermal solvor accused and in
CALC_ELEM
(for
options
FLUX_ELGA/NO_TEMP)
.
*******************************************
THERMAL CALCULATION OF SECOND MEMBER: NXACMV

TYPESE/INST: 3 5.0000000000000
LINEAR CALCULATION: F
......
-->
COMPLEMENTARY CALCULATION IN SENSITIVITY
--> BLUFF OF THE OPTION: T EAST REPLACES BY (DT/DS) -
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
59/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
--> AND ADDITION Of a NEW SOURCE TERM
--> OPTION:CHAR_SENS_EVOLNI
LPAIN/LCHIN:PGEOMER MALL .COORDO
LPAIN/LCHIN:PTEMPER &&OP0186VAR_____001
......
NO. LOADS: 6
LOOP ON THE LOADS OF THE TYPE NEUMANN FLAX
CHARGE:CH3
EXICHA/NOMCHS: 1????????
CHARGE:CH4
EXICHA/NOMCHS: 1????????
-->
COMPLEMENTARY CALCULATION IN SENSITIVITY
--> EXCHANGE/NUMCHM: 2
--> BLUFF OF THE OPTION: CARD TEXT+/- NULL
--> K: 1
--> OPTION:CHAR_THER_TEXT_F
LPAIN/LCHIN:PT_EXTF &&VECHTH.T_EXTNUL
LPAIN/LCHIN:PGEOMER MALL .COORDO
Example 6.2-2: Layout of
THER_LINAIRE
or THER_NON_LINE
,
via
INFO=2
, in the file message.
*******************************************
CALCULATION OF HEAT FLUXES
OPTION OF CALCULATION FLUX_ELGA_TEMP
MODEL MOTH
SD EVOL_THER GIVEN TH
RESULT TH
MATERIAL TAKEN INTO ACCOUNT CMAT
SEQUENCE NUMBER 1 NUMBERS
A NUMBER OF SIGNIFICANT PARAMETERS 7
*******************************************

OP0058 **********
INST/IAUX/IORDR 0. 1 0
NRPASS/TYPESE/NOPASE 1 3
PS1
CHTEMP/CHTESETH_PS1 .001.000000/TH .001.000000
--> OPTION:FLUX_ELGA_TEMP
LPAIN/LCHIN:PGEOMER MALL .COORDO
LPAIN/LCHIN:PMATERC CMAT .MATE_CODE
LPAIN/LCHIN:PTEMPER
TH_PS1 .001.000000
LPAIN/LCHIN:PTEMPSR &&MECHTI.CH_INST_R
LPAIN/LCHIN:PTEREF
LPAIN/LCHIN:PMATSEN
CMAT_PS1.MATE_CODE
LPAIN/LCHIN:PTEMSEN TH .001.000000
Example 6.2-3: Layout of
CALC_ELEM
,
via
INFO=2
, in the file message.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
60/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
6.3 Perimeter
of use
The perimeter of use of thermal calculations of sensitivities can be formulated in some points:
1. Thermal calculations of sensitivity is pressed on the operators
THER_LINEAIRE
and
THER_NON_LINE
. They deal with thermal, linear problems or not, isotropic or
anisotropic, stationary or transitory.
2. They were not still extended to the problems of drying and hydration which are too
treaty by
THER_NON_LINE
(option
THER_HYDR
or
SECH_ *
of
COMP_THER_NL/RELATION
).
If this type of problems were retained for a calculation of sensitivity, calculation stops in
ERREUR_FATALE
after having meant the reason of it.
3. In an exhaustive way, these calculations of sensitivities relate to all the meshs supports
(
TRIA3/6, QUAD4/8/9, TETRA4/10, PENTA6/13/15, PYRAM5/10
and
HEXA8/20/27
)
and all isoparametric modelings (
PLAN, PLAN_DIAG, AXIS, AXIS_DIAG, 3D
and
3d_DIAG
). If other modelings (only
THER_LINEAIRE
accept other modelings
that the usual isoparametric elements) are present in the mesh (
COQUE_
or
AXIS_FOURIER
) calculation stops (option ­ 1 front the line of the suitable options of the catalog
elements) in
ERREUR-FATALE,
after having meant the reason of it.
4. While being pressed on the perimeter of use of the code, one did not envisage the calculation of sensitivity
(of the temperature and its flow) in the presence of non-linear conductivity.
5. One has access to the sensitivity of the field of temperature (and its flow cf [§5]) compared to
all parameters of loadings and all the characteristics material. Those
having to depend that variables on space.
6. The only exceptions are the entropy (a constant entropy according to the temperature,
that would not have any direction!), multiplying coefficients of the conditions of Dirichlet
generalized (calculation of sensitivity not established because it has little direction with coefficients
often discrete) and conductivities anisotropic (data-processing limitations, absence of
modeling
THER_ORTH_FO
).
7. One can, on the other hand, calculate the sensitivity compared to one of the components of a normal flow
vectorial.
8. Calculation provides, at the same time, the field of standard temperature and the sensitivities of this same
field compared to the parameters provided to the key word
SENSITIVITY
of the thermal operator.
On the other hand, for the calculation of flow (via
CALC_ELEM
and/or
CALC_NO
), only the part is calculated
sensitivity. To also obtain that of the field of standard temperature, it is necessary to reiterate
the operation without the key word
SENSITIVITY
.
9. If the user, involuntarily or not, request a calculation of sensitivity compared to one
parameter known as “insensitive”, i.e. nonconcerned by thermal calculation in progress, one
message of alarm prevents it. No calculation is produced, the structure of data planned for
the greeting of it
CHAMNO
is initialized to zero, and calculation continues with the other sensitivities
asked.
10. The request for one or more sensitivities does nothing but enrich the structure of data
thermics (
EVOL_THER
) and provides also the field of temperature of which they are the derivative
(cf n°8). In term of performance, the calculation of an analytical sensitivity is much less
expensive that a standard calculation since the same factorized matrix is re-used.
11. During the calculation of the sensitivity of the heat flux via
CALC_ELEM
, it is necessary to specify not
only the field material used (as for the standard problem) but also them
loadings accused by the parameters of derivation.
12. To avoid any confusion (at the supervisory level, but also on the level user!), it
is not better to re-use several times a significant parameter in loadings or
different materials.
13. The calculation of the sensitivity of the heat flux was not developed, as for the problem
standard, in
COQUE_PLAN
and in
COQUE_AXIS
.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
61/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Note:
By parameterizing advisedly the loadings and materials “derived” in the file from
order, one can also have access to some derivatives made up. Thus if one
loading or a characteristic depends, explicitly or implicitly, of the time or of
space and that one is able to exhume this dependence, one can then calculate the derivative
of
T
(and of its flow) compared to this coefficient of dependence.
The main thing being which it also does not depend on the calculated solution (the field of temperature).
In which case, it will be necessary to develop a true made up calculation of sensitivity.

In short, the perimeter of use of this functionality gathers it or not thermal, linear,
isotropic or anisotropic, stationary or transitory, being pressed on finite elements
isoparametric lumpés or not. Within this framework there, it covers the same perimeter as that
accused thermal operators.

6.4 Example
of use
To familiarize oneself with the use of this new functionality, one can take as a starting point this
expurgée version of the case test
SENST04A
[V1.01.154]. In this case analytical test, it is about
to make sure of the validity of derived compared to the coefficient from exchange-wall and conductivity, in
a calculation of transitory thermal response linear of two plates separated by a play in which
a transfer of heat is carried out.
The problem is two-dimensional, but the limiting conditions make that the field of temperature reached
quickly the stationary state and depends analytically only on the X-coordinate and the data. One in
then deduced easily the analytical expressions from the sensitivities of the field of temperature and sound
flow compared to the thermal parameters which interest us.
# 1. Definition/memorizing of the constant functions
PS_UN=DEFI_CONSTANTE (VALE=1.0,);
MEMO_NOM_SENSI (NOM_UN =
PS_UN);
PS_ZERO=DEFI_CONSTANTE (VALE=0.0,);
MEMO_NOM_SENSI (NOM_ZERO =
PS_ZERO);
.........
# 2. Definition of the significant parameters and other parameters
PS1=DEFI_PARA_SENSI (VALE=80.0,);
PS2=DEFI_PARA_SENSI (VALE=40.0,);
.........
# 3.1. Installation of the materials (std and sensitive)
MAT=DEFI_MATERIAU (THER_FO=_F (LAMBDA=PS2, RHO_CP=A5))
MEMO_NOM_SENSI (NOM=_F (NOM_SD='
MAT',
PARA_SENSI=
PS2,
NOM_COMPOSE=
“MAT_PS2”));
MAT_PS2=DEFI_MATERIAU (THER_FO=_F (LAMBDA = PS_UN,
RHO_CP =
PS_ZERO))
CMAT=AFFE_MATERIAU (MAILLAGE=MAIL, AFFE=_F (TOUT=' OUI', MATER=MAT))
MEMO_NOM_SENSI (NOM=_F (NOM_SD='
CMAT',
PARA_SENSI=
PS2,
NOM_COMPOSE=
“CMAT_PS2”));
CMAT_PS2=AFFE_MATERIAU (MAILLAGE=MAIL, AFFE=_F (TOUT=' OUI', MATER=MAT_PS2))
.........
# 3.2.2 Loadings significant (echange_paroi)
CH_1=AFFE_CHAR_THER_F (MODELE=MODTHER,
ECHANGE_PAROI=_F (GROUP_MA_1=' INTERG',
GROUP_MA_2=' INTERD', COEF_H=
PS1));
MEMO_NOM_SENSI (NOM=_F (NOM_SD=
“CH_1”,
PARA_SENSI=
PS1,
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
62/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
NOM_COMPOSE=
“CH_1_PS1”));
CH_1_PS1=AFFE_CHAR_THER_F (MODELE=MODTHER,
ECHANGE_PAROI=_F (GROUP_MA_1=' INTERG',
GROUP_MA_2=' INTERD', COEF_H=
PS_UN));
.........
# 4.1 standard Calculation + 2 calculations of sensitivity
MEMO_NOM_SENSI (NOM= (_F (NOM_SD='
RESU',
PARA_SENSI=
PS1,
NOM_COMPOSE=
“RESU_PS1”),
_F (NOM_SD='
RESU',
PARA_SENSI=
PS2,
NOM_COMPOSE=
“RESU_PS2”)));
RESU=THER_LINEAIRE (MODELE=MODTHER, CHAM_MATER=CMAT,
EXCIT= (_F (CHARGE=CH_0), _F (CHARGE=
CH_1)),
TEMP_INIT=_F (CHAM_NO=TEMPINIT),
SENSIBILITE= (
PS1, PS2),
INCREMENT=_F (LIST_INST=LINST))

# 4.2 Calculation of flow to the nodes
RESU=CALC_ELEM (reuse=RESU, MODELE=MODTHER, CHAM_MATER=CMAT,
RESULTAT=RESU, OPTION=' FLUX_ELNO_TEMP',
EXCIT= (_F (CHARGE=CH_0), _F (CHARGE=
CH_1)),
SENSIBILITE= (
PS1, PS2))
RESU=CALC_NO (reuse=RESU,
RESULTAT=RESU, OPTION=' FLUX_NOEU_TEMP',
SENSIBILITE= (
PS1, PS2))
Example 6.4-1: Installation of a calculation of sensitivity in linear thermics
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
63/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
7 Outline Conclusion/
During digital simulations obtaining a rough result is not sufficient any more. The user is of
more in petitioning of calculation of sensitivity compared to the data input of the problem. That
he makes it possible to estimate the uncertainty which the field result according to the law of variation answers
data. This derivative is also the basic substrate of problems opposite (retiming of
parameters…) and of problems of optimization.
This sensitivity can be obtained “manually”, but the experiment shows that these studies
parametric are often expensive, little mutualisables and less reliable than an analytical calculation
established in the computation software.
In this note, one places oneself in the perimeter of use of the thermal operators
standards of Code_Aster and one is interested in this analytical sensitivity of the field of
temperature and of its flow compared to the characteristics material and the loadings. One y
described the process allowing to exhume the linear system which this derivative checks. In order to
to minimize the overcost calculation, a particular effort was brought to bind its resolution to that of
initial problem.
One details theoretical, numerical work and data processing which governed the establishment
of these calculations of sensitivity in the code. One specifies their properties and their limitations all in
connecting these considerations to a precise parameter setting of the accused operators and to the choices of
modeling of the code. One tried constantly to bind different the items approached while detailing,
has minimum, the a little technical demonstrations.
Required environment, the parameter setting and the perimeter of use of this news
functionality are described. An example extracted from an official case-test is clarified.
Thereafter, the prospects for this work are several commands:
·
From a functional point of view, one could extend truly the perimeter of use
thermal calculations of sensitivity (and also that of the standard problem on which they
lean) with non-linear and/or anisotropic conductivities and the conditions of exchange
wall in lumpé. With few expenses, one could also deal with the problems of hydration, of
drying and of convection-dissemination.
·
In addition, of the developments still remain to establish to be able to treat
suitably a broad spectrum of derivative made up and, in particular, to carry out
chainings thermomechanical (which are the genuine target of the current developments).
One will be able to then obtain the mechanical sensitivities of variables (displacements, deformations
and stresses) compared to loadings or characteristics material of the problem
thermics.
·
From a theoretical point of view, it remains to carry out a study “numérico-functional calculus”
similar to that of this document, to exhume the same thermal sensitivities in
modeling
HULL
and
FOURIER
. That will open another field of investigation then:
sensitivities compared to the geometrical characteristics of the structural elements.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
64/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
8 Bibliography
[1]
R. DAUTRAY & J.L. LIONS and Al Analyzes mathematical and numerical calculation for
sciences and techniques. ED. Masson, 1985.
[2]
J.L. LIONS. Some methods of resolution of the problems in extreme cases non-linear.
ED. Dunod, 1969.
[3]
H. BREZIS. Analyze functional, theory and applications. ED. Masson, 1983.
[4]
V. VENTURINI and probabilistic Al Study of the tank by a coupling mechanic-reliability engineer. Assessment of
P1-97-04 project: PROMETE. Note HP-26/99/012, Nov. 1999.
[5]
G. NICOLAS. Feasibility study of the I7-01-02 project. Note HI-72/00/020, 2000.
[6]
G. NICOLAS & J. PELLET. Structure for the sensitivity. Note HI-72/01/009.
[7]
G. CORLISS and Al AD off algorithms: from simulation to optimization. ED. Springer Verlag,
2002.
[8]
S. CAMBIER. Vibratory sensitivities of sizes. Theory and algorithms for one
establishment in Code_Aster. Note HT-62/01/011.
[9]
C. DUVAL and Al Applicability of the automatic differentiation to a system of governing EDP
thermohydraulic phenomena in a heating tube. Note 96NJ00017, 01/02/96.
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
65/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A
Appendix 1 Concept of derived “within the meaning of the distributions”

Let us place within the framework of a Bi-material and a derivation compared to the voluminal heat of one of these
materials, in linear thermics. The same reasoning can be led into non-linear and compared to
do not matter that it another characteristic material, loading or limiting condition.

One considers the modeling of following total voluminal heat (
I
I
the indicating function of ième
part
I
)
()
()
() (
)
2
2
1
2
2
1
1
,
:
+
=
X
X
X
I
I
C
p
éq
A1.1






Appear A1-a: Designation of voluminal heats for a Bi-material

That is to say the family of distributions parameterized by
()
] [
(
)
T
D
T
,
0
'
*
×
+
éq A1.2
such as
T
is the solution of the problem in extreme cases
(
P
)
(
)
(
)
(
)
] [
] [
] [
] [






=
×
=
+
×
=
×
=
×
=
-
+
+
)
(
)
0
,
(
,
0
,
0
,
0
,
0
div
1
0
3
2
1
2
2
1
1
X
X
T
T
HT
HT
N
T
G
N
T
F
T
S
T
T
T
I
I
ext.
éq
A1.3

One of course recognizes the solutions of as many slightly shifted problems of linear thermics (cf [éq 2
1] in order to be able to approximate gradually, when the small parameter
tends towards zero,
1
T
.
1
1
2
2
background image
Code_Aster
®
Version
6.0
Titrate:
Calculation of sensitivities in thermics
Date:
01/07/03
Author (S):
O. BOITEAU
Key
:
R4.03.02-A
Page
:
66/66
Manual of Reference
R4.03 booklet: Analyze sensitivity
HI-23/03/001/A

By writing that
T
0
and
T
check, respectively, (
P
0
) and (
P
), and by withdrawing members from members both
EDP one obtains
(
P
-
P
0
)
(
)
(
)
(
)
] [
] [
(
)
] [
(
) (
)
] [
(
)






=
-
×
=
-
+
-
×
=
-
×
=
-
×
-
=
-
-
-
0
)
0
,
(
,
0
0
,
0
0
,
0
0
,
0
div
0
3
0
0
2
0
1
0
1
1
0
0
X
T
T
T
T
H
N
T
T N
T
T
T
T
T
T
I
T
T
T
T
T
éq
A1.4
It any more but does not remain to divide by
1
this new EDP and to make tend towards zero the parameter
. The problem
in extreme cases then becomes a simple problem interfered similar homogeneous Cauchy-Dirichlet-Neumann-Robin type
with the initial problem and comprising however a different source term
(
)
0
1
0
1
lim
P
P
-
(
)
] [
] [
] [
] [






=
×
=
+
×
=
×
=
×
-
=
-
0
)
0
,
(
,
0
0
,
0
0
,
0
0
,
0
div
3
2
1
0
1
X
U
hu
N
U
N
U
U
T
T
I
U
T
U
éq
A1.5
which one notes
(
)
0
1
0
1
lim
:
T
T
U
-
=
the solution. The aforementioned exists and is single because it inherits all them
good properties of the initial problem. By definition, one then indicates it under the term of “distribution
derived compared to the characteristic
1
voluminal heat “
1
:
=
T
U
.

Taking into account the process of obtaining of this type of derived (translation, passage in extreme cases then definition
of a distribution derived as soon as one is ensured of his existence and his unicity), it appears that that returns
“formally” to directly derive the terms from the initial problem in extreme cases [éq A1-2], [éq A1-3].
It is in this way that one proceeded for all the derivative exhumed in this document. Moreover, for
each one of them, one checked that this “formal” derivation produced the same result as it
process describes previously.