II-2.2.2a Températures imposées
Les températures étant données, on a comme
dans le cas des problèmes de diffusion:
? q ( L , t
) ? T ( L , t
)
q ( L , t
) ?? ( T ) ? ? k (
T )[
? t ? y
a1
1
a2
?
?
0
a ?
1
0
b
?
?
T1
t
??t
(II-31a)
a M
1
a M ?
?
1
?
0
a M
?
?
0
1
t
(II-31b)
38
II-2.2.2b Flux imposés
Nous rappelons que le flux thermique est donné par la loi
de Cattanéo-Vernotte (I-15) dans le
cas d'une onde thermique. Ainsi aux frontières on a:
]en (II-32a)
]en (II-32b)
Nous nous proposons également dans le cas des
équations de propagation de l'onde thermique d'examiner en exploitant
l'équation de conservation et la relation constitutive du flux thermique
aux extrémités du mur.
CAS DES EQUATIONS DE CONSERVATION
Dans le cas où la propagation est modélisée
par (II-23a) et à partir des développements faits dans les cas
des problèmes de diffusion, on a :
HEUGANG NDJANDA Audrey Steven
Thèse de Master of science,
Option physique, Spécialité
Mécanique-Energétique/2012 Laboratoire de Mécanique et de
Modélisation des Systèmes Physiques, Université de
Dschang
2
k1
a2
?y1 2
0
a?1
?
?t
? ?y
a1 C
2?t
??? ( T ) ? 1
? 1 ?
? ? t ?
?
?
t
??t
k1
2
?y1 2
?
? T ?
( )
? ?t ?
?y
(II-33a)
??t
12
?
? ?? t ?? t ?
q t
q M M
?
)
2
aM?
aM
1
b? ? ?
0
1
?2 1 ?
kM
?
?1
?
kM
?1 2
?1 2
?yM
?yM
?
12
? ?
?
q1
t
??t
?
?
(T1)
?t
?y
a ? C
t
? ? t ?
?? t
?
?
? ?
M
?
q
?
(TM
?
??1 t
C T t??t
t??t
M
?t
?
)?y
(TM
? ??
q t??t ? q t
1 1
?y
2?t C TMt ?t
2?t
(II-33b)
Tandis que dans le cas où la propagation
modélisée par (II-23b), on a plutôt :
2
k1
?
q
?
a2
?
a?1
t
?y1 2
0
t ?? t t 1
? ? t ?
b ? C ?
C
2?t ?t 2?t
? ( T ) ? y
1
T t t t
? C T
1 1
?t
t
k1
2
b ? ? 2 ?
1 ?
?
? ? ( T ) ?
M
?(T1)C
?y
??t
C
? ?
?
2
a1
?
?
?t
?t
??
?y12
?? t
?
( T 1 )
?
t
1
? t
? ? ( T ) ?
? ( ) ? C
?
?
? ?? q t ?? t ?
q t
1 1
?
?y
(II-
? ?33a)
39
HEUGANG NDJANDA Audrey Steven
Thèse de Master of science,
Option physique, Spécialité
Mécanique-Energétique/2012 Laboratoire de Mécanique et de
Modélisation des Systèmes Physiques, Université de
Dschang
?y ?
?(T ) ?
1
a ? ? Ct??t
? C ?
M 2?t ? ?t ?
aM?
aM
1
0
kM
?1 2
?1 2
?1
?yM
t ?t
2?t
2?
1
?
kM
? 1 2
?yM
?1 2
t
? ?t
Ct TM
?
? ?? t??t ? q
t qM M
?
??t
?
(T1)
?
q
?
? ?
?
t
M
?t
(II-34b)
II-2.2.2c Convection aux frontières
Elle se traduit aux frontières du mur des problèmes
de propagations par les équations
suivantes :
Avec
Dans le cas où la propagation est modélisée
par (II-23a), on a :
2
k1
a2
?y12
0
a?1
?
??? ( T ) ? 1
? 1 ?
? ? t ?
t
??t
C
?y
2
a1
?t
? ? h1 ?
2
? y ? t
1 2
)h t
k1
?
(T1
??t
1
t
?
?
? ?
2?t ?t 2?t
h1
1
h1
?
?t
Tt
f
,
? T t??t
? ?
f,1
?
? T
? ? ?
( ) ?
1
b ?2 1 ?
? ?t ?
Ct Tt ? ?(T1)
?y Ct??tTt??t
?y
(II-35a)
40
HEUGANG NDJANDA Audrey Steven
Thèse de Master of science,
Option physique, Spécialité
Mécanique-Energétique/2012 Laboratoire de Mécanique et de
Modélisation des Systèmes Physiques, Université de
Dschang
?
? ?
|
|
h t??t
?(TM) ? h
?t
|
? ?? t t
? ? h T
?
|
|
??? ( T ) ?
M ? 1
?
k M?1 2
? ? h ? h
M M
aM
? ?(T ) ?
M
b? ?2 ?
1 ?
? ?t ?
?y ? (T
)?y
t t M C T
t??t t ??t
C T ?
M M
2?t ?t 2?t
?
t
??t
C
?y
t
)
??t
(T M
?t
(II-35b)
? t ?
2 ? t ?
?y M ? 12
Tandis que dans le cas où elle modélisée
par (II-23b), on a plutôt :
?
?
? ?
M
? ?? t t
? ? h T
?
? ??t
h t??t
?(TM) ? h
? t
k1
a2
?y1 2
0
a?1
?
t
?
?
2
)h1
(T1
??t
h1
?
?
? ?
?t
t
h1
1
Tt
f
,
? t ?
2?t ?
a1
?
M M
T t t
? ?(T ) ?
1 C t ??t
? ? C ?
? ?t ?
2
) h t
? ?t
? Tt??t
? ?
f,1
?
? ? h1 ??y1 2
k1
?y
f , M M f ,
M
?y
2 ? t ? t 2 ?
t
?
??t
(T1
?t
1
0
? 1 2
kM
1
2
?yM
?
(C
t?
C
t
)
?
M?
(T
t
Tt
2?t ?t 2?t
b?
t
t
kM
?
??t
12
?
)
??t
?
(T M
?y
?
?(T1) C
? hM ?
a M
hM
C
? ?
?
2
?t
?
b??(Ct??t?Ct)?(T1)?C
?
?t
?t
? ?
?
?
?yM
? 1
aM
1
t
? ?t
Ct TM
?1 2
?y
(II-36a)
(II-36b)
41
M M
T t t
f,M M f , M
CAS DES RELATIONS CONSTITUTIVES DES
FLUX
Lorsqu'il s'agit des problèmes de propagation avec flux
imposés aux frontières, exprimée par la relation
constitutive (loi de Cattanéo-Vernotte) donnée par:
HEUGANG NDJANDA Audrey Steven
Thèse de Master of science,
Option physique, Spécialité
Mécanique-Energétique/2012 Laboratoire de Mécanique et de
Modélisation des Systèmes Physiques, Université de
Dschang
Les équations discrétisées se
réécrivent respectivement en (y=0) et (y=1) quelle que soit la
k k ? ? ( T )
?
? M 1 2
T t t M
?? 1 2 ? ?? ??
? ? ? ? ??
t t
? ?? ? q t
T t t q t t M q
M M ? 1 M M M
? y ? y ? ? t
M ? 1 2 M ? 1
2 ?
forme considérée des problèmes de
propagation (d'enthalpie et de température):
soient
(II-47b)
Les coefficients des systèmes algébriques sont
alors donnés respectivement en (y=0) et
(y=L) par :
(II-48b)
(II-48a)
2
k1
1 2 (II-46)
a1
6y12
2
k1
a2
a-1
6y1 2
0
?
b
q
t+At + z (T )
1 At
?
(q1+At q1 )
?
?1 2
kM
?1 2
aM
?yM
0
aM?
1
?1 2
kM
1
?1 2
aM
?
?yM
?
?
)
??t
?
(TM
?
q
t
M
?
?t
b
42
? ?
? ?
?
? ??
t??t ? q t qM M
II-2.3. DISCRETISATION DES CONDITIONS INITIALES
La condition initiale des problèmes de diffusion est
donnée par la relation (I-29a) du chapitre précédent. Sa
forme discrétisée est:
TP ? Tinitiale(II-45)
Pour les problèmes de propagation, il faut en plus de
(II-36) ajouter la forme discrétisée de l'équation (I-30b)
du précédent chapitre.
TP ? TP
HEUGANG NDJANDA Audrey Steven
Thèse de Master of science,
Option physique, Spécialité
Mécanique-Energétique/2012 Laboratoire de Mécanique et de
Modélisation des Systèmes Physiques, Université de
Dschang
|