Projet M2 Mathematiques Appliquées UAD

UNIVERSITÉ ALIOUNE DIOP ANNÉE 2022-2023
MASTER II MATHÉMATIQUES ET APPLICATIONS : UFR SATIC
EXAMEN SIMULATION NUMÉRIQUE DES EDP

Problème 1 (Approximation du modèle EDP de type elliptique)
On considère le problème :
−∆u = f ∀x ∈ ]0, 1[
(1)
u(0) = u(1) = 0
Partie I : Formulation variationnelle : existence et unicité de solution
Soit Th : {jh}M +1
j=0 avec (M + 1)h = 1 une partition uniforme de l’intervalle [0, 1] en sous-intervalles
Ij = ((j − 1) h, jh) , j = 1, 2, ..., M + 1.
soit Vh0 l’espace de dimension finie telle que
Vh0 = {v ∈ C(0, 1) : v est une f onction linéaire par morceaux sur Th , v(0) = v(1) = 0},
avec des fonctions de base {ϕj }M j=1 .
Z 1
1
On définit V0 = H0 (0, 1) = {w : (w(x)2 + w0 (x)2 )dx < ∞, w(0) = w(1) = 0}.
0
1. Montrer que la formulation variationnelle du problème (1) est de trouver u ∈ V0 telle que :
Z 1 Z 1
u0 (x)v 0 (x)dx = f (x)v(x)dx; ∀v ∈ V0
0 0
2. Montrer que la stabilité de cette formulation variationnelle est contrôlée si u ∈ H 1 (0, 1) solution
du problème (1) vérifie k u kH 1 (0,1) ≤ C k f kL2 (0,1) où C est une constante.
Z 1 Z 1
0 0
3. On pose a(u, v) = u (x)v (x)dx et l(v) = f (x)v(x)dx.
0 0
a. Montrer inégalités suivantes :
-
a(u, v) ≤k u kH 1 (0,1) k v kH 1 (0,1) ;
-
1
a(v, v) ≥ k v k2H 1 (0,1) ;
1 + c2p
-
l(v) ≤k f kL2 (0,1) k v kH 1 (0,1) .
b. En déduire que ∀u, v ∈ H01 (0, 1, l’existence et l’unicité de la solution du problème (1) est assurée.
Partie II : Interpolation-méthode des éléments finis
On donne la méthode des éléments finis de Galerkin du problème (1) comme suit :
Z 1 Z 1
Trouver U (x) ∈ Vh0 telle que U 0 (x)v 0 (x)dx = f (x)v(x)dx; ∀v(x) ∈ Vh0 qui permet de déter-
0 0
miner
ξj = U (xj )l’approximation de u(x) aux nœuds xj = jh, 1 ≤ j ≤ M.
1
4. Montrer que l’approximation de la formulation variationnelle s’exprime
M Z 1 Z 1
X 0
0
ξj ϕj (x)v (x)dx = f (x)v(x)dx, ∀v(x) ∈ Vh0
j=1 0 0
On rappelle que v(x) est une combinaison linéaire de fonctions de base ϕi (x), pour i = 1, ..., M .
5. Montrer que le problème est de résoudre le système M × M équations
A= {aij }M T
i,j=1 b= (b1 , b2 , ...., bn ) à déterminer.
 x−(i−1)h
 h (i − 1)h ≤ x ≤ ih
On définit ϕi (x) = (i+1)h−x
 h ih ≤ x ≤ (i + 1)h
0 sinon
6.
a. Donner les éléments de ϕ0i (x) pour (i − 1)h ≤ x ≤ ih et ih ≤ x ≤ (i + 1)h
b. Pour |i − j| > 1, ϕi et ϕj ont des supports disjoints, vérifier alors que aij =0.
c. Montrer que aii= h2 , ∀i = j.
d. Montrer que ai,i+1 = ai+1,i = − h1 , ∀j = i ± 1.
e. En tenant de la symétrie de A, donner la forme générale de la matrice de rigidité.
Partie III : Implémentation sur Freefem+ et Matlab
On considère toujours le même modèle :

−∆u = f dans Ω
u = 0 sur ∂Ω
1. Implémenter le problème Elliptique pour Ω = ]0, 1[, f (x) = 1
2. Ω = [0, 1] × [0, 1], f (x, y) = sin(πx) sin(πy), h = 1/100.
Problème 2 (Approximation de l’équation de Stokes )
soit Ω ⊂ R (N = 2ou3) un domaine borné et régulier.

On considère le problème suivant

 ν∆u + ∇p = f dans Ω
(P1 ) ∇.u = 0 dans Ω
u = 0 sur ∂Ω

On introduit les formes bilinéaires suivantes :

a : H01 (Ω)N × H01 (Ω)N −→ R
b : H01 (Ω)N × L2 (Ω) −→ R
on introduit aussi l’espace des fonctions de L2 (Ω) à moyenne nulle par : L20 (Ω) = {q ∈ L2 (Ω), Ω qdx = 0}
R
1. Montrer que la formulation du problème (P1 ) est de trouver u ∈ H01 (Ω)N , p ∈ L20 (Ω) tels que
a(u, v) + b(v, p) = (f, v) ∀v ∈ H01 (Ω)N
b(u, q) = 0 ∀q ∈ L20 (Ω)
2
où les formes bilinéaires sont à expliciter.
On pose X = [H01 (Ω)]N , Y = L20 (Ω),
2. Vérifier que la forme a est continue et coercive sur [H01 (Ω)]N × [H01 (Ω)]N .
3. Vérifier que la forme b est continue sur [H01 (Ω)]N × L2 (Ω).
Rappels
Lemme : Il existe une constante C telle que pour tout p ∈ L2 (Ω), il existe v ∈ [H01 (Ω)]N tel que
∇.v = p et kvkH 1 (Ω)N ≤ Ckpk2L (Ω).
pdx = 0 alors on peut prendre v ∈ [H01 (Ω)]N .

R
De plus, si P vérifie Ω
Condition ’inf-sup’ : D’après le Lemme, il existe w ∈ [H01 (Ω)]N , q = −∇.w, b(w, q) =

− Ω (∇.w)dx = kqk2L2 (Ω) vérifiant les mêmes conditions que v tels que
R
b(w, q) 1 b(w, q)
≥ = βkqkL2 (Ω) .
kwkH 1 (Ω)N C kqkL2 (Ω)
4. Montrer qu’il existe une constante C > 0 indépendante de h (pas de discrétisation) telle que

ku − uh kX ≤ inf ku − vh kX + inf kp − qh kY
vh ∈Vh qh ∈Yh
où Xh ∈ X, Yh ∈ Y et Vh = {vh ∈ Xh |b(vh , qh ) = 0, ∀qh ∈ Yh }.

5. Montrer que
inf ku − vh kX ≤ C inf ku − wh kX
vh ∈Vh wh ∈Xh
6. Montrer que (choix pour bonus +2)

kp − ph kY ≤ C ku − uh kX + inf kp − qh kY .
qh ∈Yh
7. En déduire que (choix pour bonus +2)

ku − uh k + kp − ph kY ≤ C inf ku − vh kX + inf kp − qh kY .
vh ∈Vh qh ∈Yh
8. Implémenter le Stokes sur Freefem+ en utilisant les données de vos choix.
Indication disponibles dans les exemples du dossier Freefem+ et dans plusieurs ressources.
BON N E CHAN CE

Projet M2 Mathematiques Appliquées UAD

Transféré par

Droits d'auteur :

Formats disponibles

Projet M2 Mathematiques Appliquées UAD

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Projet M2 Mathematiques Appliquées UAD

Transféré par

Droits d'auteur :

Formats disponibles

UNIVERSITÉ ALIOUNE DIOP ANNÉE 2022-2023

MASTER II MATHÉMATIQUES ET APPLICATIONS : UFR SATIC

EXAMEN SIMULATION NUMÉRIQUE DES EDP

Partie II : Interpolation-méthode des éléments finis

On considère toujours le même modèle :

2. Ω = [0, 1] × [0, 1], f (x, y) = sin(πx) sin(πy), h = 1/100.

Problème 2 (Approximation de l’équation de Stokes )

soit Ω ⊂ R (N = 2ou3) un domaine borné et régulier.

On introduit les formes bilinéaires suivantes :

a(u, v) + b(v, p) = (f, v) ∀v ∈ H01 (Ω)N

b(u, q) = 0 ∀q ∈ L20 (Ω)

On pose X = [H01 (Ω)]N , Y = L20 (Ω),

∇.v = p et kvkH 1 (Ω)N ≤ Ckpk2L (Ω).

pdx = 0 alors on peut prendre v ∈ [H01 (Ω)]N .

Condition ’inf-sup’ : D’après le Lemme, il existe w ∈ [H01 (Ω)]N , q = −∇.w, b(w, q) =

où Xh ∈ X, Yh ∈ Y et Vh = {vh ∈ Xh |b(vh , qh ) = 0, ∀qh ∈ Yh }.

6. Montrer que (choix pour bonus +2)

7. En déduire que (choix pour bonus +2)

8. Implémenter le Stokes sur Freefem+ en utilisant les données de vos choix.

Vous aimerez peut-être aussi