[PDF] Analyse numérique des EDP TD 1





Previous PDF Next PDF



Feuille dexercices : Formulations Faibles

Donnez sa formulation variationnelle. 2. Peut-on appliquer le Théorème de Lax-Milgram ? Exercice 6. Soit ?: H1(?)? L2(??) l 



Corrigé de la Séance 2 : Formulations variationnelles

Dans la suite ? est un ouvert borné de R3



Sorbonne Université Année 2019-2020 Master MPE mention

Méthodes pour les EDP P. Frey. Corrigé de l'examen du 8 janvier 2020. Exercice 1. partie



Analyse numérique des EDP TD 1

29 janv. 2016 Avec certains corrigés ... Exercice 1 (Défaut de coercivité dans C1) ... Formulation variationnelle et existence de la solution.



Chapitre 5 ´ETUDE MATH´EMATIQUE DES PROBL`EMES

Exercice 5.2.1 A l'aide de l'approche variationnelle démontrer l'existence intégration par partie on obtient la formulation variationnelle suivante :.



Méthodes variationnelles

Définition 3.5 (Formulation variationnelle) Soit f ? L2(?); Exercice 39 (Conditions aux limites de Fourier et Neumann) Corrigé en page 130.



Exercices Corrigés - Analyse numérique et optimisation Une

29 août 2012 FORMULATION. VARIATIONNELLE DES. PROBL`EMES ELLIPTIQUES. Exercice 2.1.1 Si f est une fonction continue sur [01]



METHODES DAPPROXIMATION DES EQUATIONS AUX

1.1 Définition d'une équation aux dérivées partielles (e.d.p) . . . . . . . . 7 2.4 Exercices . ... 7 Etablissement de formulations variationnelles.



Université Lyon 1 Année 2013-2014 Master Mathématiques

J(v). Formulation variationnelle de problèmes elliptiques. Exercice 4. (Laplacien + Dirichlet). Soit ? un ouvert de Rn borné et régulier (de 



Interrogation no2 - corrigé groupe 1 Lundi 13/02/17.

13 févr. 2017 On considère l'EDP elliptique d'ordre 2 suivante : ... de telle sorte que toute solution u de (2) satisfasse la formulation variationnelle.



[PDF] Corrigé de la Séance 2 : Formulations variationnelles - ENSTA Paris

Construire la formulation variationnelle (FV1) associée `a (1) Corrigé de la question 1 : En multipliant la 1`ere équation de (1) par v ? H1(?) et en 



[PDF] Feuille dexercices : Formulations Faibles

1 Donnez sa formulation variationnelle 2 Peut-on appliquer le Théorème de Lax-Milgram ? Exercice 6 Soit ?: H1(?)? L2(??) l'application trace sur ??



[PDF] M2-0120 corrigépdf - Sorbonne Université

Méthodes pour les EDP P Frey Corrigé de l'examen du 8 janvier 2020 Exercice 1 partie on obtient la formulation variationnelle suivante :



[PDF] Exercices Corrigés

29 août 2012 · FORMULATION VARIATIONNELLE DES PROBL`EMES ELLIPTIQUES Exercice 2 1 1 Si f est une fonction continue sur [01] montrer que l'équation dif 



[PDF] Méthodes variationnelles

Définition 3 5 (Formulation variationnelle) Soit f ? L2(?); Exercice 39 (Conditions aux limites de Fourier et Neumann) Corrigé en page 130



[PDF] Analyse numérique des EDP TD 1

29 jan 2016 · Avec certains corrigés Exercice 1 (Défaut de coercivité dans C1) Formulation variationnelle et existence de la solution



[PDF] Université Lyon 1

J(v) Formulation variationnelle de problèmes elliptiques Exercice 4 (Laplacien + Dirichlet) Soit ? un ouvert de Rn borné et régulier (de 



Examen corrige Formulation variationnelle

Calcul Scientifique Séance no3 Formulations variationnelles Corrigé 29 Novembre 2005 Exercice 1 Formulation variationnelle 1 1 - Soit v Interrogation 



Formulation variationnelle - Cours et Exercices

est une forme sur V (´linéaire si les conditions aux limites de l'EDP de départ le sont) Remarque : Une formulation 



[PDF] TD 6 EDP elliptiques

EDP elliptiques Exercice 1 Identification de formulations variationnelles On considère le problème suivant : { trouver u ? V tel que

:

Analyse numérique des EDP - TD1 1

Master MAPI

31ereannée Université Paul Sabatier - Toulouse 3 Année 2015-2016

Analyse numérique des EDP

TD 1

Avec certains corrigés

Les numéros de Théorèmes, Propositions, etc ... font référence aux notes de cours.

Exercice 1 (Défaut de coercivité dansC1)On considèreX=fv2 C1([0;1]); v(0) =v(1) = 0gmuni de sa norme usuelle et l"énergieEdéfinie par

E(v) =12

Z 1 0 jv0j2dxZ 1 0 fv dx; oùfest une fonction continue qui est identiquement nulle sur un intervalle non trivial[;][0;1]. Montrer qu"il existe une suite(vn)nd"éléments deXde la forme v n(x) =1pn '(n(xx0)); avec'etx0bien choisis, telle que(E(vn))nest bornée maiskvnkX!n!11.Corrigé :

On prendx0=+2

,'une fonctionC1c(R)à support dans]1;1[non identiquement nulle. Ainsi la fonctionvnest bien nulle en0et1pour toutnassez grand. On a immédiatement, pour nassez grand,kvnkL1=k'kL1pn etkv0nkL1=pnk'0k1. On a donc bienkvnkX!

1quandn! 1.

Calculons l"éner giede vn

E(vn) =k2

Z 1 0 jv0n(x)j2dxZ 1 0 f(x)vn(x)dx:

Par hypothèse surf, et par construction devn, les supports defet devnsont disjoints et donc le second terme est

nul. Il nous reste donc à évaluer la première intégrale. On utilise la définition devnet un changement de variable

y=n(xx0)

E(vn) =k2

Z R nj'0(n(xx0))j2dx=k2 Z R j'0(y)j2dy:

Ainsi la suite(E(vn))nest bien bornée (elle est même constante sur cet exemple!!).Exercice 2 (Défaut de complétude dansC1)On considère toujours l"espaceX=fv2 C1([0;1]); v(0) =v(1) = 0gmais muni, vette fois, de la norme

kvkH1=qkvk2L2+kv0k2L2.

Vérifier que la suite(un)ndéfinie par

u n(x) =s x12 2 +1n r1 4 +1n est une suite de Cauchy non convergente dans(X;k kH1).Corrigé :

Commençons par montrer que la suite (un)nn"est pas convergente dans(X;kkH1). Pour cela, on suppose qu"elle

converge vers un certainu2X. D"après le lemmeI.8 , ceci implique que,unconverge uniformément versusur

[0;1]. En particulier, on aurait la convergence simple deunversusur[0;1]. On voit alors immédiatement que cela impliqueu(x) =x12 12 . Or, cette fonctionun"est pas de classeC1sur [0;1]à cause de la singularité enx=12 . Ceci établit donc une contradiction.

F. BOYER- VERSION DU29JANVIER2016

2 Analyse numérique des EDP - TD1

Montrons maintenant que (un)nest de Cauchy dansH1. Pour cela on commence par utiliser le lemmeI.8 pour

établir que,unetun+pétant dansX, on a

Ainsi, pour montrer le critère de Cauchy dansH1pour(un)nil suffit de vérifier que(u0n)est de Cauchy dansL2.

Calculons doncu0nu0n+p

u

0n(x)u0n+p(x) =x1=2q

(x1=2)2+1n x1=2q (x1=2)2+1n+p = (x1=2)q(x1=2)2+1n+pq(x1=2)2+1nq (x1=2)2+1n q(x1=2)2+1n+p (x1=2)

1n+p1n

q (x1=2)2+1n q(x1=2)2+1n+p q(x1=2)2+1n +q(x1=2)2+1n+p

On prend la valeur absolue de cette expression et on essaie de majorer intelligemment les différents termes. Comme

on sait qu"il n"y a pas convergence uniforme de cette suite de fonctions (car sinon la limiteu0serait continue ...), la

majoration qu"on doit obtenir doit encore dépendre dex! De façon plus précise, on utilise les inégalités p(x1=2)2+ 1=(n+p) jx1=2j et p(x1=2)2+ 1=(n+p) +p(x1=2)2+ 1=np(x1=2)2+ 1=n:

Il vient ainsi

ju0n(x)u0n+p(x)j 1n (x1=2)2+1n 1n

1(x1=2)2+1n

78

1(x1=2)2+1n

18 1n 11 n 78

1jx1=2j14

=1n

1=81jx1=2j14

Si maintenant on élève l"inégalité au carré et qu"on l"intègre entre0et1, il vient Z 1 0 ju0nu0n+pj2dx1n 1=4Z 1

01jx1=2j12

dx:

Comme l"intégrale qui apparaît dans le membre de droite est convergente (c"est une intégrale de Riemann), on a

bien montré ku0nu0n+pkL2Cn 18

ce qui prouve bien que la suite(u0n)nest de Cauchy dansL2, et donc que(un)nest de Cauchy dans(X;k kH1)

d"après la remarque initiale.Exercice 3 (FonctionsC1par morceaux)SoitIun intervalle ouvert borné deRet2I. Montrer qu"il n"existe pas de fonctiong2L2(I)telle que

'() =Z I 'g dx;8'2 C1c(I): En déduire que sif:I!Rest une fonction de classeC1par morceaux, alors f2H1(I)()fest continue:

Que vaut la dérivée faible defdans le cas où elle appartient àH1(I)?F. BOYER- VERSION DU29JANVIER2016

Analyse numérique des EDP - TD1 3

Corrigé :

Il s"agit d"un raisonnement par l"absurde. On suppose qu"il existe une fonctiongqui vérifie l"égalité de l"énoncé

'() =Z I 'g dx;8'2 C1c(I):

On va construire une suite de fonctions'n2 C1c(I)telle que'n() = 1pour toutnet telle que les intégralesR

I'ng dx

tendent vers0.

Il s"agit de les choisir de plus en plus concentrée autour de. Plus précisément, on fixe unr >0tel que[r;+r]

Iet on va fixer une fonction 2 C1c(I)vérifant () = 1etSupp [r;+r](une telle fonction existe, voir le

lemme II.12 du cours).

On définit alors

n(x) = (n(x));8x2I:

On vérifie alors que

n() = 1;8n1

Supp'n[r=n;+r=n]I;

et donc ce sont bien des fonctions à support compact dansI. De plus, comme les supports se concentrent autour de,

nous avons n(x)!n!10;8x6=:

Comme par ailleurs, nous pouvons écrire

j'n(x)g(x)j k k1jg(x)j; et queg2L1(I), on peut appliquer le théorème de convergence dominée et en déduire que Z I ng dx!n!10; ce qui constitue bien une contradiction.

Pour la seconde partie du résultat, on se donne une fonctionfde classeC1par morceaux. Pour simplifier (mais cela ne

change rien en pratique) nous supposons que la fonctionfa un seul point de discontinuité. Ainsi siI=]a;b[, on suppose

qu"il existe2Ietfg;fdde classeC1surRtout entier tels que f=fg;sur[a;[; f=fd;sur];b]:

Ceci implique (voir la figure

1 ) en particulier quefadmet des limites à gauche et à droite ennotéesf()etf(+)et qui vérifient f(+) = limx!x>f(x) =fd(); f() = limx!x

Ix=f()f(+)f

gf df

FIGURE1 - Une fonction de classeC1par morceaux

F. BOYER- VERSION DU29JANVIER2016

4 Analyse numérique des EDP - TD1

On se donne maintenant une fonction test'2 C1c(I)et on calcule à l"aide de la relation de Chasles et d"une

intégration par parties sur chacun des deux intervalles : Z I f'0dx=Z a f'0dx+Z b f'0dx Z a f g'0dx+Z b f d'0dx = [fg']aZ a f0g'dx[fd']bZ b f0d'dx = [f()f(+)]'()Z b a (f0g1[a;]+f0d1[;b])'dx: On a utilisé ici que la fonction test's"annule aux bornes de l"intervalle'(a) ='(b) = 0.

Le second terme est bien de la forme recherchée (une intégrale d"une fonction deL2contre la fonction test'). Ainsi,

la fonctionfsera dansH1si et seulement si on peut mettre sous la même forme le premier terme. Or la première partie

de l"exercice nous dit que cela n"est possible que si ce terme est nul, c"est-à-dire si le saut defenest nul, c"est-à-dire

si et seulement si f(+) =f(): Cette condition est exactement la condition de continuité def. Ainsi, sifest continue elle est bien dansH1et sa dérivée faible vaut xf(x) =( f0g(x)six < f

0d(x)six > :

On constate que celle-ci n"est pas bien définie au pointce qui n"est pas surprenant car il s"agit d"une fonction deL2

(qui n"est donc bien définie qu"à un ensemble de mesure nulle près).Exercice 4 (Dérivation d"un produit)

SoitIun intervalle ouvert borné deR. Montrer que le produit de deux fonctionsu;v2H1(I)est encore un

élément deH1(I)et on a l"identité suivante (au sens faible) x(uv) = (@xu)v+u(@xv): En particulier, on peut intégrer par parties au sens faible Z b a (@xu)v dx= (uv)(b)(uv)(a)Z b a u(@xv)dx;8[a;b]I:

Corrigé :

Il s"agit ici d"utiliser les propriétés de densité des fonctions régulières dansH1. Ainsi, siu;vsont deux éléments de

H

1(I), le ThéorèmeI.12 nous dit qu"il e xistedeux suites de fonctions (un)n;(vn)n C1(I)qui convergent dansH1

versuetvrespectivement. D"après la définition de la normeH1, cela signifie kunukL2!n!10; k@xun@xukL2!n!10; kvnvkL2!n!10; k@xvn@xvkL2!n!10: Pournfixé, nous avons bien la formule de dérivation usuelle (unvn)0=u0nvn+unv0n;

et donc, pour une fonction test'2 C1c(I)fixée, nous pouvons intégrer par parties et obtenir la formule

Z I u nvn'0dx=Z I (u0nvn+unv0n)'dx=Z I (@xunvn+un@xvn)'dx: On souhaite maintenant justifier le passage à la limite dans cette égalité.

F. BOYER- VERSION DU29JANVIER2016

Analyse numérique des EDP - TD1 5

Pour le premier terme, on peut écrire par exemple, en utilisant l"inégalité de Cauchy-Schwarz,

Z I u nvn'0dxZ I uv'0dx=Z I (unvnuv)'0dx=Z

I(unu)vn+u(vnv)'0dx

k'0kL1kunukL2kvnkL2+kvnvkL2kukL2; et cette quantité tend bien vers0d"après les hypothèses. On traite de façon analogue les deux autres termes et on obtient, à la limite que Z I uv'0dx=Z I (@xuv+u@xv)'dx: Il reste à observer que, d"après le corollaire I.14 , les fonctionsuetvsont dansL1et que donc, la fonction@xuv+u@xv est dansL2.

Au final, on a bien prouvé que le produituvest dansH1et que sa dérivée faible vaut@xuv+u@xv.

La formule d"intégration par parties en découle grâce au troisième point du Théorème

I.12 appliqué au produit uv.Exercice 5 (Dérivation faible de fonctions composées) SoitIun intervalle ouvert borné deR. Montrer que siF2 C1(R;R)etu2H1(I), alorsF(u) =Fuest aussi une fonction deH1(I)et que nous avons, au sens faible x(F(u)) =F0(u)@xu;presque partout:Corrigé :

On procède comme précédemment en introduisant une suite(un)nde fonctions de classeC1(I)qui converge,

dansH1, versu. Lesunétant régulières, on a bien la formule de dérivation des fonctions composées@x(F(un)) =

F

0(un)@xunce qui, après intégration par parties contre une fonction test'2 C1c(I)nous donne l"égalité

Z I

F(un)'0dx=Z

I

F0(un)(@xun)'dx:

On veut maintenant justifier le passage à la limite dans ces intégrales.

On va utiliser le Corollaire

I.14 qui nous dit que (un)nconverge versuuniformément. Ceci implique tout d"abord l"existence d"un intervalle compact[R;R]tel que u n(x)2[R;R];8n;8x2I:

Sur le compact[R;R], les fonctionsFetF0sont continues et donc uniformément continues (théorème de Heine). On

en déduit donc queF(un)(resp.F0(un)) converge uniformément versF(u)(resp.F0(u)). Comme par ailleurs, on a par

définition la convergence dansL2de@xunvers@xu, on voit que l"on peut bien passer à la limite dans tous les termes de

l"égalité et aboutir à Z I

F(u)'0dx=Z

I

F0(u)(@xu)';8'2 C1c(I):

Ceci justifie bien queF(u)2H1(I)et que sa dérivée faible est donnéeF0(u)@xu(qui est bien dansL2carF0(u)est

F. BOYER- VERSION DU29JANVIER2016

6 Analyse numérique des EDP - TD1

bornée). Exercice 6 (Autour de l"inégalité de Poincaré) Dans cet exercice on travaille sur un intervalleI=]0;L[avecL >0. 1. On souhaite démontr erl"inégalité de P oincarésuivante

8v2H10(]0;L[);kvkL2(I)L

k@xvkL2(I):(1) (a) Montr erque si on sait pr ouver(1)pour les fonctionsv2 C1c(]0;L[), alors on peut en déduire le résultat souhaité. (b) Démontr erl"inégalité (1)pour les fonctionsC1c(]0;L[)en utilisant les séries de Fourier. (c)

Vérifier que la constante

L estoptimale, c"est-à-dire que c"est la plus petite constante pour laquelle

l"inégalité(1)est vraie. On pourra expliciter une fonctionvpour laquelle(1)est une égalité.

2.

On intr oduitun nouvel espace fonctionnel

H

1m(]0;L[) =(

v2H1(]0;L[);Z L 0 v(x)dx= 0)

En utilisant un procédé similaire au cas précédent, montrer que l"inégalité de type Poincaré suivante est

également vraie, avec une constante optimale :

8v2H1m(]0;L[);kvkL2(I)L

k@xvkL2(I):

Montrer que l"applicationv2H1m(I)7! k@xvkL2(I)est une norme surH1m(I)équivalente àk kH1(I).Exercice 7

On travaille dans l"intervalleI=]0;1[. Soientqetfdeux fonctions continues surI. On considère le problème

suivant : (u00+q(x)u=f(x);8x2]0;1[; u(0) =u(1) = 0;(2) 1. Vérifier que si q(x) =2alors le problème est mal posé, c"est-à-dire : (a) P ourf= 0, il existe des solutionsunon nullesde(2). (b) Il e xistedes fonctions fpour lesquelles(2)n"a aucune solution.

Autrement dit, le problème(2)n"admet pas toujours des solutions, et s"il en existe elles sont nécessairement

non uniques. 2. A partir de maintenant on suppose que la fonction qvérifie inf

Iq >2:

Proposer une formulation variationnelle dansH10(I)du problème(2). Donner l"expression d"une fonc-

tionnelle "énergie"u7!E(u)dont la solution de(2), si elle existe, serait un minimiseur.3.Montr erque la fonctionnelle Eobtenue est continue pour la topologie deH10(I).

4.

Montr erque Eest minorée surH10(I). En déduire l"existence d"une suite minimisante(un)n, i.e. une suite

vérifiantE(un)!n!1inf v2H10(I)E(v). 5. Démontr erqu"il e xisteun unique u2H10(I), tel queE(u) = inf v2H10(I)E(v). Montrer que la fonctionuainsi obtenue est bien une solution faible du problème(2). 6. Démontr erque la fonction uest de classeC2et queuest bien une solution du problème(2)au sens classique. 7.

Quelle(s) étape(s) de la démonstr ationtombe(nt) en défaut quand q(x) =2(car d"après la question 1

le problème est mal posé)?F. BOYER- VERSION DU29JANVIER2016

Analyse numérique des EDP - TD1 7

Exercice 8 (Un exemple de problème aux limites non linéaire) On poseI=]0;1[. On cherche dans cet exercice à résoudre le problème suivant (u00+u3=f;dansI; u(0) =u(1) = 0;(3) oùfest une fonction continue donnée etun paramètre réelpositifdonné.

1.Unicité :Soientu1etu2deux fonctions de classeC2solutions du problème(3). Démontrer que l"on a

Z 1 0 ju01u02j2dx+Z 1 0 (u31u32)(u1u2)dx= 0:

En déduire que l"on a nécessairementu1=u2.

La suite de l"exercice est donc consacrée à la preuve de l"existence d"une solution à ce problème.

2.Formulation variationnelle et existence de la solution

(a)

Démontr erl"inégalité suivante

8a;b2R;a+b2

4 12 a4+12 b4: (b)

Montr erque toute solution éventuelle u2 C2([0;1])du problème(3)vérifie la formulation suivante

8v2H10(I);Z

1 0 u0v0dx+Z 1 0 u3v dx=Z 1 0 fv dx:(4) On pourra commencer par considérer desv2 C1c(I). (c) On s"intér essedonc maintenant au pr oblèmesuivant Trouveru2H10(I)vérifiant les équations(4):(P) i.

P ourquoine peut-on pas utiliser le théorème de Lax-Milgr ampour résoudr ele pr oblème(P)?

ii. P ourtout v2H10(I), on introduit la fonctionnelleénergiesuivante

E(v) =12

Z 1 0 jv0j2dx+4 Z 1 0 jvj4dxZ 1 0 fv dx;8v2H10(I): Expliquez rapidement pourquoi tous les termes de cette formule sont bien définis. iii. Démontr erque l"éner gieEest minorée sur l"espaceH10(I). On pourra essayer de minorerE

par une fonctionnelle dont les propriétés sont connues.iv.Démontr erqu"il e xisteune suite minimisante pour EdansH10(I). On note(un)nH10(I)une

telle suite.quotesdbs_dbs13.pdfusesText_19
[PDF] formulation variationnelle exercices corrigés pdf

[PDF] pecheur d'islande film

[PDF] madame chrysanthème

[PDF] pecheur d'islande film 1996

[PDF] ramuntcho

[PDF] aziyadé

[PDF] cours modélisation et simulation des systèmes pdf

[PDF] différence entre modélisation et simulation

[PDF] modélisation et simulation cours

[PDF] modélisation et simulation cours informatique

[PDF] modélisation et simulation pdf

[PDF] pierre et jean résumé court

[PDF] pierre et jean personnages

[PDF] fonction affine activité

[PDF] analyse pierre et jean