EMD AnalyseNumérique GCM1 9JUIN 2021

Télécharger au format docx, pdf ou txt
Télécharger au format docx, pdf ou txt
Vous êtes sur la page 1sur 5

Université Saad Dahlab de Blida Juin 2021

Département de génie des procédés


Analyse Numérique GCM1 2020-2021 Examen de moyenne durée N°1
1 heure

A. Donner la formule générale d’une équation différentielle aux dérivées partielles du second
ordre de la fonction u=u(x,y) 1pt

B.Donner pour chacune des équations suivantes : l’ordre, linéaire (homogène ou non) ou non
linéaire. 7 pts

1. uxy -6uxx+7uy-ux+10x =0
3
2. x 2uxyz -8uxx+6y 4uy-ux+9xe x =0

3. uxy+uy+u2-xy = 0

4. uxuxyy-2xuuyy-3xyuy-u4=0

5. uxyy-2xyuyy = -xy4ux

6 .log(uy u xy)- log(uxuyx)= c ; c = constante

∂ 2 ∂ 2
7.
∂y
(u x +u y ) + 4 u xy u xy −4 u x u y −2 x= ( u −u ) +4 u yx u yx−u x −2 x =0
∂y x y

C-Classer les équations aux dérivées partielles suivantes par type : parabolique, elliptique,
hyperbolique. 5pts

1. uy - 3uxx = 0

2. uy=x3 uxx

3. uxx+2uxy+5uyy-32u=0

4. (x+1)uxx+2uxy+(x+1)uyy+4ux+5uy=0

5. uxx+2(x+1) uxy-uyy+2ux+2uy=0

D) Soit uxx+ uyy = f(x,y) donner le résultat final à 10-2 7 pts

L’équation aux dérivées partielles est valide sur un domaine rectangulaire (R1), de l= 10 cm
de largeur, L=20 cm de longueur

1. Discrétiser l’équation pour ∆x=h, ∆y=k, après prendre h=k=5 cm.

2. Classer par leur type ces conditions u(x,0) =0 , u(x,10)=0 , u(0,y)=40 , u(20,y)=100 .

3. Représenter les précédentes conditions sur les frontières puis résoudre E.D.P pour les cas
suivants : f(x,y)=0 , f(x,y )= -1
Université Saad Dahlab de Blida Juin 2021
Département de génie des procédés
Analyse Numérique GCM1 2020-2021 Examen de moyenne durée N°1
1 heure

Solution

A.

N ordre Linéaire Homogène


1 2 oui non
2 3 oui non
3 2 non
4 3 non
5 3 oui oui
6 1 oui oui
7 2 non

B. On prendra ∆=B2− A . C ou bien ∆=(2∗B)2−4 A .C

∆ Type d’équation
1 0 Parabolique
2 0 Parabolique
3 <0 Elliptique

4. ∆=-x^2-2x

x
-∞ -2 0 +∞
∆ - + -
elliptique hyperbolique elliptique
∆=0 parabolique parabolique

5. ∆=x^2+2x+2 cette équation n’admet pas de racines dans R, ∆ > 0 équation


hyperbolique.

Forme générale d’une équation différentielle partielle du second ordre d’une fonction
u=u(x,y)

F (x , y , u ,u ¿ ¿ x , u y ,u xx , u xy ,u yy )=0 ¿

Auxx+2Buxy+Cuyy+Dux+Euy+Fu=G(x,y)
Université Saad Dahlab de Blida Juin 2021
Département de génie des procédés
Analyse Numérique GCM1 2020-2021 Examen de moyenne durée N°1
1 heure

( )
1
1
D. uxx+ uyy = f(x,y) équation de Poisson ; Après discrétisation 2
1 −4 1 ui , j=f i , j
h 1

Sous forme matricielle fi,j= 0 équation de Laplace avec condition du premier type ou
bien Dirichlet.

( )( ) ( )
−4 1 0 u1 −40
1 −4 1 . u 2 = 0
0 1 −4 u 3 −100

( )( ) ( )
−4 1 0 u1 −40+ h2 . fij
2
Si fi,j≠ 0 , h2= 25 1 −4 1 . u 2 = 0+h . fij
0 1 −4 u 3 2
−100+h . fij

f ij=0 u 1=12 , 50u 2=10 ,00 u 3=27 , 50

f ij=−1 u1=21 , 43 u 2=20 , 71u 3=36 , 43


Université Saad Dahlab de Blida Juin 2021
Département de génie des procédés
Analyse Numérique GCM1 2020-2021 Examen de moyenne durée N°1
1 heure

En plus

∂z ∂z 3 y x
9. Soit a +x + y z + x e = y définir pour l’E.D.P l’ordre, linéaire (homogène ou
∂x ∂y
non) ou bien non linéaire
2 2
3 ∂ z ∂ z δz ∂ ❑z ∂ z
a=constante , 2, xy , z y , , , , , cos ( z ) , sin ⁡( ❑ )
∂ x ∂ y δx ∂ x ∂ y ∂x ∂ y

Résolution Matlab f=0

uxo =0 ;uxl =0 ; uoy= 40 ; uLy=100 ;h=5 ; f= 0 ; h2=h^2*f; a=[-4 1 0 ;1 -4 1 ;0 1 -4]; b=[-


uoy+h2; 0+h2 ; - uLy+h2] ; x=a\b

x=

12.5000

10.0000

27.5000

>> uxo =0 ;uxl =0 ; uoy= 40 ; uLy=100 ;h=5 ; f= -1 ; h2=h^2*f; a=[-4 1 0 ;1 -4 1 ;0 1 -4];


b=[-uoy+h2; 0+h2 ; - uLy+h2] ; x=a\b

x=

21.4286

20.7143

36.4286
Université Saad Dahlab de Blida Juin 2021
Département de génie des procédés
Analyse Numérique GCM1 2020-2021 Examen de moyenne durée N°1
1 heure

Vous aimerez peut-être aussi