Exercices de cours du chapitre IV : Formulations intégrales - approximation
Exercice IV-2: problème de thermique
Objectifs : Utiliser la formulation variationnelle et l’approximation sur un problème de physique.
On considère un problème stationnaire de conduction thermique. Les équations pour un milieu isotrope Ω de
frontière ∂Ω sont :
G
⎧divq − r = 0 dans Ω
⎪ GG G JJJJG
(P) ⎨ [Link] = Φ d sur ∂ 2 Ω ∂ avec le vecteur flux de chaleur q = −λ grad T dans Ω
⎪T =T sur ∂ 1Ω
⎩ d
et r une source interne de production de chaleur,
Td température imposée sur sa frontière ∂ 1Ω ,
Φ d un flux de chaleur imposé sur son bord ∂ 2 Ω
Montrer que le problème (P) est équivalent à la formulation variationnelle suivante :
JJJJJJJG JJJJJJJJJG
∫ gradT .gradδ T − rδ T ) dV +
(λ ∫ Φ d δ TdS = 0 ∀δ T tel que δ T = 0 sur ∂ 1Ω
Ω ∂Ω
Et T = Td sur ∂1Ω (thermiquement admissible)
G G G JJJJJJG G GG
On rappelle que div (α f ) = α divf + f . gradα et ∫ divf dV = ∫ f .next dS .
Ω ∂Ω
Formulation unidimensionnelle
Pour étudier le transfert de chaleur dans un mur de hauteur Φd
r (x)
Td
infinie, on étudie le problème unidimensionnel où la température x = −l / 2 x =l/2
ne dépend plus que d’un paramètre spatial x .
On considère le domaine représenté sur la figure ci-contre
Identifier les domaines : Ω , ∂ 1Ω et ∂ 2Ω .
Montrer que, dans ce cas, la formulation précédente se réduit à :
dT d (δ T )
l/2 l/2
∫ λ dx dx dx = ∫ r ( x) δ Tdx + Φ d δ T x= −
l
−l / 2 −l / 2 2
Préciser les conditions à satisfaire pour établir cette équation.
On considère un élément fini linéaire à deux nœuds i et j de longueur le .
L’approximation de la température sur cet élément est alors : T ( x ) = Tiϕi ( x ) + T j ϕ j ( x)
1 x 1 x
où les fonctions d’interpolation linéaires sont : ϕi ( x ) = − et ϕ j ( x ) = + .
2 le 2 le
Discrétiser la formulation précédente et calculer :
[ ]
La matrice raideur K ;
Le vecteur flux généralisé { Fr } dû à une production uniforme de chaleur r ( x ) = r0 sur l’élément ;
Le vecteur flux généralisé { Fd } dû à un flux imposé Φ d au nœud i de l’élément.
En utilisant cet élément fini, résoudre le problème présenté sur la figure suivante.
Φd 1 2 3
Td = 0
Il n’y a pas de source de chaleur, l’élément (1-2) est de longueur l et a pour conductivité thermique λ , l’élément (2-3)
est de longueur l / 3 et a pour conductivité thermique 2λ .
Déterminer la répartition de température dans les éléments, et calculer le flux de chaleur inconnu au point 3.
5
Exercices de cours du chapitre IV : Formulations intégrales - approximation
Corrigé de l’exercice IV-2: problème de thermique
Formulation variationnelle d’un problème de thermique stationnaire
G
⎧divq − r = 0 dans Ω
⎪ GG r >0
Équations : ⎨ [Link] = Φ d sur ∂ 2 Ω Source de chaleur T
⎪T =T sur ∂ 1Ω
⎩ d
G G
Où q (flux de chaleur) est donné par : q = −λ grad T dans Ω
G G
divq + r = 0 ⇔ ∀δ T ∫ (divq − r ) δ T dV = 0 Représente l’erreur pondérée sur le domaine
Ω
G G G JJJJJJG
On rappelle que : div (α f ) = α divf + f . gradα
G G JJJJJJJJJG
D’où ∀δ T ∫ (div(δ Tq ) − [Link] δ T − rδ T ) dV = 0
Ω
G GG
Et que : ∫ divf dV = ∫ f .next dS
Ω ∂Ω
JJJJJJJG JJJJJJJJJG GG
D’où ∀δ T ∫ (λ gradT .gradδ T − rδ T ) dV + ∫ q.nδ T dS = 0
Ω ∂Ω
GG
En tenant compte de la condition de flux sur la frontière∂ 2Ω : q.n = Φ d
Et en utilisant un champ virtuel thermiquement admissible δ T = 0 sur ∂ 1Ω
JJJJG JJJJG
Nous obtenons ∀δ TTh − admissible ∫ λ gradT . gradδ T dV − ∫ r δ TdV + ∫ Φ d δ TdS = 0
Ω Ω ∂ 2Ω
Il faudra satisfaire la condition T = Td sur ∂ 1Ω
Cas unidimensionnelle JJG G
φD = −φD n Avec φD > 0 flux entrant
φd r(x)
−A / 2
T
A/2 d
x
Pour ce problème la formulation se réduit à :
dT d (δ T )
l/2 l/2
∫ λ
dx dx
dx = ∫ r ( x) δ Tdx + Φ d δ T x= −
l
−l / 2 −l / 2 2
⎧δ T = 0
Avec en x = A / 2 ⎨
⎩ T = Td
Élément fini
⎧⎪ Ti ⎫⎪ 1 x 1 x ⎧⎪ Ti ⎫⎪
L’approximation de la température est : T ( x ) =< ϕ ( x ) > ⎨ ⎬ =< − + >⎨ ⎬
⎩⎪ j ⎭⎪
T 2 Ae 2 A e ⎩⎪T j ⎭⎪
Ae / 2
λ ⎡1 −1⎤
La matrice raideur [K ] = ∫ < ϕ, x >T λ < ϕ, x > dx =
−A e / 2
A e ⎢⎣ −1 1 ⎥⎦
Ae / 2
⎧A / 2 ⎫
Le vecteur flux généralisé {FR } = ∫ < ϕ >T ro dx = ro ⎨ e ⎬
−A e / 2 ⎩A e / 2 ⎭
⎧Φ d ⎫
Le vecteur flux généralisé { Fd } = Φ d < ϕ ( − A e / 2) > = ⎨
T
⎬
⎩ 0 ⎭
6
Exercices de cours du chapitre IV : Formulations intégrales - approximation
Application
Soit le problème présenté sur la figure suivante
Φd 1 A 2 A/3 3
Td = 0
λ 2λ
⎡ 1 −1 0 ⎤
λ⎡1 −1⎤ 6λ ⎡ 1 −1⎤ λ ⎢ −1 7 −6 ⎥
[ K1 ] = ⎢−1 1 ⎥ et [ K 2 ] = ⎢ −1 1 ⎥ Î [ K ] = ⎢ ⎥
A⎣ ⎦ A ⎣ ⎦ A
⎢⎣ 0 −6 6 ⎥⎦
⎡ 1 −1 0 ⎤ ⎧ T1 ⎫ ⎧Φ d ⎫ ⎧ 0 ⎫
λ⎢ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪
Le système d’équations est :
⎢ −1 7 −6 ⎥⎥ ⎨T2 ⎬ = ⎨ 0 ⎬ + ⎨ 0 ⎬
A ⎪ ⎪ ⎪ ⎪ ⎪ ⎪
⎣⎢ 0 −6 6 ⎦⎥ ⎩ 0 ⎭ ⎩ 0 ⎭ ⎩−Φ 3 ⎭
λ ⎡ 1 −1⎤ ⎧T1 ⎫ ⎧Φ d ⎫ ⎧0 ⎫ Les températures
La résolution de : ⎢ ⎥ ⎨ ⎬ = ⎨ ⎬+⎨ ⎬ sont positives car le
A ⎣ −1 7 ⎦ ⎩T2 ⎭ ⎩ 0 ⎭ ⎩0 ⎭ flux est entrant
7A A
Nous donne les températures T1 = Φ d et T2 = Φ d
6λ 6λ
Puis la dernière équation nous donne le flux de chaleur inconnu Φ 3 = Φ d