MENUM : Thermique stationnaire
Objet : Utilisation des Méthodes Numériques dans le cadre de la MEF
Dans ce projet vous aurez l'occasion de développer un nouvel élément fini pour résoudre
un problème de thermique stationnaire. Puis vous testerez différentes techniques
d'intégration numérique pour résoudre le problème inverse qui consiste à déterminer les
consignes de chauffe pour obtenir une température donnée sur une pièce placée dans le
four.
Documents et supports
Site MEF : [Link]
Chapitre V du polycopié de cours
MEFlab : avant de débuter ce projet, vous devez avoir étudié le document de présentation
de cet outil et compris le fonctionnement des différents scripts que vous téléchargerez sur
le site.
Livres disponibles à la bibliothèque de l’école
G. Dhatt & G. Touzot : une présentation de la méthode des éléments finis.
Les presses de l'université Laval - Québec, 1981
G. Dhatt - G. Touzot & E Lefrancois : méthode des éléments finis.
Hermes Lavoisier, 2005.
Travail à réaliser : en binôme ou personnel
La durée prévue est de 12 heures en salle informatique
Scripts que vous avez développés (ils doivent être commentés)
Synthèse et analyse des résultats obtenus
À envoyer par mail dans un dossier compressé à votre nom.
Optimisation d'un four
On considère un four de forme carré comportant un certain nombre de résistances. L’objectif est
de déterminer, en régime permanent, les consignes en densité surfacique de chaleur des résistances
de sorte que la température d’une pièce placée dans le four soit proche d'une consigne Tc fixée.
Données
0,25 Γ2 Résistance On notera
Ω le domaine représentant le four
Four
P le domaine représentant la pièce
k le coefficient de diffusion thermique
k = 1 dans la pièce
Γ o 0,4
Pièce à Γo
2
chauffer k = 10 dans le four (air)
1 Pour chaque résistance i (i = 1..N r )
0,4 ri la densité surfacique de chaleur de la
Γ1
résistance (inconnues du problème)
2 Les résistances sont de diamètre φ = 0,1
Données géométriques
1
Les conditions aux limites sur la frontière du four sont :
Flux de chaleur supposés nuls sur les bords Γ o
Température imposée à T1 = 10°C sur la frontière Γ1
Température imposée à T2 = 50°C sur la frontière Γ 2
Que signifie physiquement l'hypothèse sur Γ o ?
A : Formulation EF du problème
1°: Écrire la formulation variationnelle de ce problème en se plaçant dans le cadre des hypothèses
de thermique stationnaire pour un milieu homogène isotrope (loi de Fourier).
2°: En déduire les expressions matricielles du modèle éléments finis de ce problème
Écrire les scripts "Q4_th" et "T3_th" : calcul des termes élémentaires [ K e ] et { Fe }
Écrire le script "therm2D" de résolution du problème thermique 2D
3°: Valider vos programmes sur le cas test ci-dessous
T2 = 50°C
Déterminez la solution analytique de ce problème.
Comparer avec les résultats du calcul éléments finis.
Four Maillage automatique du four
Flux Source Flux
2 surfacique nul
Étude de convergence
nul
r
A partir de la solution analytique calculer la consigne r
pour que la température au centre du four soit égale à
T1 = 10°C
100°C.
2
Décomposer ce problème en deux problèmes élémentaires :
Pb1 : Températures imposées nulles sur les frontières, et densité surfacique de chaleur
r donnée.
Pb2 : Températures imposées sur les frontières avec r = 0 .
Retrouver la consigne r pour que la température au centre du four soit égale à 100°C, à partir des
solutions de ces deux problèmes élémentaires.
B : Résolution du problème inverse
Il faut généraliser l'approche utilisée dans le cas test précédent pour passer d'un problème à une
inconnue r à un problème à N r inconnues ri , consignes des résistances de sorte que la température
dans la pièce soit proche de la température fixée Tc.
Soit TO la solution du problème correspondant à ri = 0 (∀i = 1..N r ) respectant les températures
imposées sur les frontières du four.
Par linéarité on cherche la solution T ( r ) comme une combinaison linéaire de N r problèmes
élémentaires correspondant à une température imposée nulle sur les frontières du four pour des
consignes: ri = 1 , rj = 0 (∀j ≠ i ) .
Nr
Soit T ( r ) = TO + ∑ ri Ti
i =1
1
La fonctionnelle à minimiser sur le domaine de la pièce à chauffer est: J ( r ) = ∫ (T (r ) − TC ) dS
2
2P
2
1°: Décrire les N r + 1 problèmes élémentaires à traiter sous forme directe.
1 T
2°: Montrer que J ( r ) peut s'écrire sous forme matricielle : J ( r ) = {r} [ A]{r} − {r} {b} + Cte
T
2
avec Aij = ∫ Ti T j dS et bi = ∫ Ti (TC − T0 )dS
P P
3°: Comment obtenez vous les coefficients ri optimaux ?
4°: Réalisation du script Matlab permettant de résoudre le problème inverse.
Script de type 1:
H1 : La dimension des résistances est supposée suffisamment petite, pour pouvoir les
assimiler à une source de chaleur ponctuelle de densité ri s , avec s surface de la résistance.
H2 : On calcul les coefficients Aij et bi à partir des valeurs moyennes sur les éléments
Script de type 2:
On conserve l'hypothèse H1, et on peut au choix utiliser Matlab ou une technique EF pour
l'intégration des coefficients Aij et bi .
Script de type 3:
On conserve l'hypothèse H2, mais on modélise les résistances par des éléments finis.
Pour chaque script
On comparera la prévision T ( r ) au résultat du calcul obtenu en appliquant au modèle EF les
consignes optimales ri obtenues.
Valeur des températures sur la pièce et dans le four
Température moyenne dans la pièce
On effectuera une étude de convergence pour différents maillages.
Scripts à me fournir
Synthèse et analyse des résultats obtenus, validation Samcef si vous avez le temps.