Scribd Logo
Importer

Bienvenue sur Scribd !

  • 22 vues

    Examdec 08

    Transféré par

    zizou Saada

    Droits d'auteur :

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

    Examdec 08

    Transféré par

    zizou Saada
    22 vues3 pages

    Titre original

    Examdec08

    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é ?

    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
    22 vues3 pages

    Examdec 08

    Transféré par

    zizou Saada

    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
    sur 3

    Université Pierre et Marie Curie Année universitaire 2008-2009

    M1 Master mention mathématiques et applications - Cours MM006


    Bases des méthodes numériques - E. Godlewski, M. Vohralik

    Examen du 15 décembre 2008


    durée 3 heures - aucun document n’est autorisé, pas de calculatrice
    Il est demandé de justifier précisément les résultats et de rédiger soigneusement.

    Le barême est indiqué à titre indicatif

    Il est demandé de rédiger la question Scilab sur une feuille séparée

    Question 1 (9 points) On considère le problème aux limites : trouver u ∈ C 1 ([0, 1]) ∩


    C 2 (]0, 1[) tel que

    d2 u du

     − 2 + b(x)
     = f (x) dans ]0, 1[,
    (P) dx dx


    u(0) = 0, u(1) = β,

    avec b ∈ C 1 ([0, 1]), b0 ≤ 0 f ∈ C 0 ([0, 1]).


    1) Définir une formulation variationnelle (FV) de (P) dans un espace de Hilbert V
    bien choisi contenu dans H 1 (0, 1). Donner l’expression de la forme bilinéaire a et de
    la forme linéaire ` (on remarquera que a n’est pas symétrique). Montrer que (FV) a
    une solution, et que celle ci est unique (énoncer le théorème permettant de l’affirmer et
    vérifier les hypothèses).
    2) Dans quel espace (contenant C 0 ([0, 1])) peut-on prendre f en gardant un sens à la
    formulation variationnelle ? Discuter la régularité de u.
    3) Utiliser la méthode des éléments finis P1 pour définir une approximation varia-
    tionnelle du problème dans un espace de dimension finie Vh qu’on explicitera. Préciser
    une base wi de Vh . Montrer que ce problème approché a une solution unique notée uh .
    4) Rappeler le résultat fondamental dans le cadre abstrait permettant de majorer
    l’erreur ||u − uh || pour la norme de V .
    5) Ecrire le problème discret sous la forme d’un système linéaire à résoudre. Préciser
    la matrice, le vecteur inconnu et le second membre de ce système : on donnera des
    expressions théoriques en fonction de a, ` et des wi .
    6) On suppose f = 0 et b constante. Calculer explicitement la solution de (P). Calculer
    explicitement la solution du système linéaire précédent en supposant le maillage uniforme
    de pas h.

    Question 2 (6 points) On se donne un pas de temps ∆t > 0, un pas d’espace ∆x > 0,


    ∆t
    on pose λ = ∆x , on suppose λ constant, c est une constante donnée. On considère le

    1
    schéma d’approximation
    (λc)2
     n+1
     vj = vjn − λc 2
    n
    (vj+1 n
    − vj−1 )+ 2
    n
    (vj+1 n
    − 2vjn + vj−1 )), n ≥ 0, j ∈ Z,
    (1)
    vj0 donné , ∀j ∈ Z.

    1) Si dans (1), vjn désigne une approximation d’une solution exacte u au point xj =
    j∆x, tn = n∆t, quelle est l’équation aux dérivées partielles vérifiée par u pour que (1) soit
    un schéma consistant d’ordre au moins un ?(justifier le résultat par des développements
    précis).
    2) Est ce que (sous une éventuelle condition sur ∆t, ∆x) le schéma (1) est stable dans
    L2 ? (on pourra ou redémontrer ou s’appuyer sur un résultat du cours).
    3) Ecrire vjn+1 donné dans (1) comme une combinaison linéaire des trois valeurs
    n
    vj+1 n
    , vjn , vj−1 avec des coefficients notés respectivement a1 , a0 , a−1 qu’on précisera en
    fonction de c, λ. Que valent a−1 + a0 + a1 et a−1 − a1 ? Le schéma est-il stable dans
    L∞ (le principe du maximum est-il vérifié) ?
    4) On se donne vj0 = 1, ∀j ≤ −1 et vj0 = 0, ∀j ≥ 0. Calculer vj1 , ∀j ∈ Z. Représenter
    vj0 , vj1
    (en fonction des valeurs de j placées sur un axe), conclure.

    Question Scilab (5 points) à rédiger sur une copie séparée


    On considère ici un problème de réaction–diffusion sur l’intervalle ]0, 1[

     −u00 (x) + c(x)u(x) = f (x),
    (RD) u(0) = 0,
    u(1) = 0,

    où les fonctions f et c satisfont f ∈ L2 (0, 1) et c ∈ C 0 ([0, 1]). Rappelons que sa formu-
    lation variationnelle est
    u ∈ H01 (0, 1) telle que, pour tout v ∈ H01 (0, 1), on ait

     Trouver
    Z 1 Z 1 Z 1
    (F VRD ) 0 0
     u (x)v (x)dx + c(x)u(x)v(x)dx = f (x)v(x)dx
    0 0 0

    et que la méthode des éléments finis pour ce problème consiste à trouver



     Trouver uh ∈ V1h telle que, pour tout vh ∈ V1h , on ait
    Z 1 Z 1 Z 1
    (EFRD ) 0 0
     uh (x)vh (x)dx + c(x)uh (x)vh (x)dx = f (x)vh (x)dx.
    0 0 0

    On prend N ∈ N∗ , h = 1/(N +1), xk = kh, k = 0, . . . , N +1, Ik =]xk , xk+1 [, k = 0, . . . , N


    et V1h ⊂ H01 (0, 1) l’espace des fonctions linéaires par morceaux et continues. La méthode
    des éléments finis amène au système linéaire
    Ah ũh = bh ,
    où ũh = (uh (x1 ), . . . , uh (xN ))T .

    2
    1. Ecrire la matrice Ah comme Ah = Dh + Rh . Calculer explicitement (à la main)
    les coefficients de la matrice Rh , due au terme de réaction, jusqu’à la forme des
    intégrales.
    2. Ecrire une fonction Scilab calculant la matrice Rh . L’entête de cette fonction sera
    function R=CalculMatR(N,c).
    Arguments d’entrée : le nombre de points N et le nom de la fonction c. Pour la
    fonction c, on pourrait considérer

    si 0 ≤ x ≤ 21 ,
    
    c(x) = x2
    c(x) = x3 si 12 ≤ x ≤ 1.

    Argument de sortie : un tableau carré R contenant Rh . Pour évaluer les intégrales,


    utilisez la formule d’intégration numérique suivante :
    b M
    b−aX
    Z
    g(x)dx ' Hi g(ξi ),
    a 2 i=1

    où ξi sont les points d’intégration, M ∈ N définit leur nombre et Hi donnent des
    poids pour chaque point d’intégration, avec :

    M = 3,
    H1 = 0.555555555555556, H2 = 0.888888888888889, H3 = 0.555555555555556,
    ξ1 = a + (1 − q) ∗ (b − a)/2, ξ2 = (a + b)/2, ξ3 = a + (1 + q) ∗ (b − a)/2,
    q = 0.774596669241483.

    Ecrivez cet algorithme dans le langage Scilab comme si vous le programmiez dans
    l’ordinateur.
    Ceci est une question algorithmique – cela n’est pas la syntaxe exacte du
    langage Scilab, mais l’algorithmique qui compte !

    Vous aimerez peut-être aussi