Scribd Logo
Importer

Bienvenue sur Scribd !

  • 14 vues

    Sujet 16

    Transféré par

    Sarra Atrous
    Ce document décrit un problème d'optimisation de chauffage d'un four rectangulaire comportant des résistances. L'objectif est de trouver les valeurs des résistances pour que la température dans une partie du four soit proche d'une valeur cible, en résolvant l'équation de la chaleur. Le document présente les étapes pour résoudre le problème direct puis inverse à l'aide du logiciel FreeFem++.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd

    Sujet 16

    Transféré par

    Sarra Atrous
    14 vues2 pages
    Ce document décrit un problème d'optimisation de chauffage d'un four rectangulaire comportant des résistances. L'objectif est de trouver les valeurs des résistances pour que la température dans une partie du four soit proche d'une valeur cible, en résolvant l'équation de la chaleur. Le document présente les étapes pour résoudre le problème direct puis inverse à l'aide du logiciel FreeFem++.

    Titre original

    sujet16

    Copyright

    © © All Rights Reserved

    Formats disponibles

    PDF, TXT ou lisez en ligne sur Scribd

    Avez-vous trouvé ce document utile ?

    Ce contenu est-il inapproprié ?

    Ce document décrit un problème d'optimisation de chauffage d'un four rectangulaire comportant des résistances. L'objectif est de trouver les valeurs des résistances pour que la température dans une partie du four soit proche d'une valeur cible, en résolvant l'équation de la chaleur. Le document présente les étapes pour résoudre le problème direct puis inverse à l'aide du logiciel FreeFem++.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    14 vues2 pages

    Sujet 16

    Transféré par

    Sarra Atrous
    Ce document décrit un problème d'optimisation de chauffage d'un four rectangulaire comportant des résistances. L'objectif est de trouver les valeurs des résistances pour que la température dans une partie du four soit proche d'une valeur cible, en résolvant l'équation de la chaleur. Le document présente les étapes pour résoudre le problème direct puis inverse à l'aide du logiciel FreeFem++.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    Vous êtes sur la page 1sur 2

    Mini-Projet de MAP 431, Ecole Polytechnique

    sujet proposé par F. Nataf


    nataf@ann.jussieu.fr

    Optimisation du chauffage d’un four


    On considère un four de forme rectangulaire, comportant des résistances. On se propose de chercher
    les valeurs à donner aux résistances de sorte que la température dans une partie du four soit proche d’une
    valeur fixée à l’avance.

    On notera Ω l’ouvert représentant le four et S la pièce à chauffer. Le four Ω est un carré de 2m de côté.
    L’origine est placée au centre du carré. La pièce S à chauffer est un rectangle de 1m sur 0.4m placée au
    centre du four. On notera Ci (i = 1, . . . , 4) les résistances modélisées comme des cercles de rayon 0.05m et
    placées aux points de coordonnées (±0.75m, ±0.75m) (cf. Figure 1).

    Le bord supérieur du four est supposé maintenu à une température fixe de 50o C, le bord inférieur à une
    température fixe de 10o C. Les deux bords latéraux sont isolés, le flux de chaleur y est nul. La température
    dans le four est régie par l’équation de la chaleur

    −∇ · (k∇T ) = f. (1)

    Dans cette équation, k est le coefficient de diffusion thermique qui est variable et vaut 1 dans la pièce à
    X 4
    chauffer et 10 dans le reste du four, f est la source de chaleur et est de la forme f = αi 1Ci , où les
    i=1
    coefficients αi sont à déterminer et 1C est la fonction caractéristique du domaine C.

    Travail à réaliser

    1 Problème direct
    On suppose connus les coefficients (αi )i=1,...,4 et l’on cherche à calculer la température qui en résulte.
    Question 1. Ecrire la formulation variationnelle correspondante. On notera T (α) la solution du problème
    correspondant, où α est le vecteur (αi )1≤i≤4 .
    Question 2. Ecrire le programme FreeFem++ correspondant, en prenant pour valeurs numériques
    ∀i = 1, . . . , 4, αi = 10. Calculer la température moyenne Tmoy dans la pièce S.

    2 Problème indirect
    On se donne une température souhaitée TC dans la pièce S. On cherche les coefficients αi de sorte que
    T (α) soit proche de TC dans S. Plus précisément, on cherche à minimiser
    1
    Z
    J(α) := |T (α) − TC |2 . (2)
    2 S
    Question 3. Montrer que par linéarité, il existe cinq champs de température Ti , i = 0, . . . , 4, tels que
    4
    X
    T (α) = T0 + αi Ti .
    i=1

    Calculer les champs de température Ti (i = 0, . . . , 4) à l’aide de FreeFem++.

    1
    y

    T = 50o C
    C2 C1

    Ω
    S x
    ∂T k=1 ∂T
    ∂n =0 ∂n =0
    k = 10
    C3 C4
    T = 10o C

    Figure 1 – Four Ω et pièce à chauffer S

    Question 4. Montrer que le calcul des coefficients optimaux αi peut se mettre sous la forme d’un système
    linéaire à résoudre : Aα = b, où A est une matrice 4 × 4 et b un vecteur colonne. Expliciter A et b.

    Question 5. Calculer effectivement les coefficients αi pour TC = 40o C.

    Question 6. Contrôler à l’aide d’une simulation directe sous FreeFem++ que votre solution est acceptable.
    En particulier, calculer l’erreur relative entre la température moyenne obtenue et la température moyenne
    attendue. Dans tous les cas, on cherchera un maillage convenable en expérimentant avec l’adaptation de
    maillage de FreeFem++.

    Vous aimerez peut-être aussi