Chapitre 1 EDP Classification

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

GEER1

Chapitre I
Module
Maths pour Classification des EDPs
ingénieurs I

Professeur :
Younes ABOUELHANOUNE
Chapitre 1
EDP d’ordre 2

1. Classi…cation et théorèmes
EDP Elliptique

2. Résolution des EDP à coe¢ cients constants.

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Chapitre 1
EDP d’ordre 2

1. Classi…cation et théorèmes
EDP Elliptique
EDP Parabolique

2. Résolution des EDP à coe¢ cients constants.

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Chapitre 1
EDP d’ordre 2
1. Classi…cation et théorèmes
EDP Elliptique
EDP Parabolique
EDP Hyperbolique

2. Résolution des EDP à coe¢ cients constants.

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Chapitre 2
Méthode de la transformée de Fourier pour
l’équation de la chaleur

1 1- Produit de convolution
2 2- Transformation de Fourier
3 3- Application : Equation de la chaleur
@u @2u
=a 2 ; a>0
@t @x

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Chapitre 3

Méthode de la séparation des variables


1 pour les EDP
2

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Chapitre 1

EDP d’ordre 2
Classi…cation et Méthodes de Résolution

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
De…nition 1

Equations di¤érentielles ordinaires (EDO)


Les équations di¤érentielles ordinaires se rencontrent dans tous
les domaines de la physique :
- Electricité,
- Mécanique,
- Thermique. . . .
C’est une relation entre une fonction inconnue et ses dérivées.
La fonction inconnue ne dépend que d’une seule variable, par
exemple f(x).

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Etablissement du courant électrique dans une bobine

L’évolution du courant électrique est donnée par


di(t)
L + Ri(t) = u(t) avec i(t = 0) = 0
dt

La résolution de cette équation di¤érentielle donne :


U t L
i(t) = 1 e avec =
R R

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Oscillateur en mécanique du point :

La projection de la relation sur l’axe des x nous amène à


l’équation di¤érentielle du 2ème ordre
d2 x f dx k dx
+ + x = 0 avec x(t = 0) = x0 et jt=0 = 0
dt2 m dt m dt

La résolution de l’équation di¤érentielle du second


qordre (régime
f k 2 :
pseudopériodique) donne avec 2 = m et w = ( m )

t
x(t) = x0 e cos(wt) + sin(wt)
w
Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Qu’est-ce qu’une EDP ?
De…nition 2
la modélisation des phénomènes physiques repose sur la
résolution d’équations aux dérivées partielles.
- Mécanique des ‡uides : équations de Navier-Stockes
- Electromagnétisme : équations de Maxwell
- Thermique : équation de la chaleur
- Mécanique quantique : équation de Schrödinger

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
De…nition 3

Le caractère particulier d’une équation aux dérivées partielles


(EDP) est de mettre en jeu des fonctions de plusieurs variables
(x; y; ::::) ! u(x; y; ::::): Une EDP est alors une relation entre les

variables et les dérives partielles de u.

De…nition 4

Une équation aux dérivées partielles relie une fonction inconnue


à ses dérivées. La fonction inconnue dépend de plusieurs
variables (variables d’espace et le temps).

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
EDP - Quelques exemples classiques et historiques
✔ Equation des ondes (1747)
∂2Φ 2
(t, x) − c ∆Φ(t, x) = f (t, x)
∂t2
- c vitesse de l’onde, f (t, x) force externe
- Corde vibrante : Φ(t, x) champ de déplacement vertical
- Acoustique linéaire : Φ(t, x) champ de pression du fluide
- Equation linéaire non homogène du second ordre
J. le Rond d’Alembert (1717-1783)
- Capteurs à corde vibrante : jauge de contrainte extensométrique (génie civil)
✔ Equation de la chaleur (1811)

∂Φ
(t, x) − c2 ∆Φ(t, x) = f (t, x)
∂t

- Φ(t, x) champ de température, c2 coefficient de diffusivité thermique


- f (t, x) source de chaleur
- Equation linéaire non homogène du second ordre
- Diffusion de la chaleur, surface isolée, pas de convection
- Problème aux limites unidimensionnel : Φ(t, 0) = 0 et Φ(t, l) = 0 J. Fourier (1768-1830)
EDP - Quelques exemples classiques et historiques
✔ Equation dynamique de la poutre d’Euler-Lagrange (1750 - vibrations de flexion)

∂2Φ ∂4Φ
µ 2 (t, x) + EI 4 (t, x) = f (t, x)
∂t ∂x

- µ = ρS masse par unité de longueur, E module de Young, I moment d’inertie


- Φ(t, x) déformée verticale, f (t, x) charge distribuée par unité de longueur
- Equation linéaire non homogène d’ordre 4
- Théorie de l’élasticité isotrope pour un matériau homogène
L. Euler (1707-1783)
- Pas de cisaillement et pas de gauchissement
✔ Equations de l’élastodynamique (1821)

2 3  
∂ Φ ∂ X ∂Φ ∂Φj
ρ(x) − λ(x)divΦ − µ(x) + = fi (t, x), i = 1, 2, 3
∂t2 ∂xi j=1
∂xj ∂xi

- ρ densité de masse, λ, µ coefficients de Lamé


- Φ(t, x) champ de déplacement, f (t, x) force externe par unité de longueur
- Equation linéaire non homogène d’ordre 2
- Théorie de l’élasticité isotrope pour un matériau homogène
C.L.M.H. Navier (1785-1836)
- Ondes sismiques, vibrations des structures
De…nition 5
Un problème est bien posé (ou stable) si une « petite »
variation des conditions de bord entraîne une « petite »
variation de la solution.

De…nition 6
Soit E un espace vectoriel de fonctions. Un opérateur linéaire
sur E est une application A dé…nie sur E et telle que :
1- A( u) = A(u) pour tout u 2 E
2- A(u + v) = A(u) + A(v) pour tout u; v dans E.

Exemples 7

Soit E l’espace vectoriel des fonctions réelles de classe 1


X 2 C sur R.
@u @ u
* L’application A dé…nie par Au = @xi et u = @x2i
sont des
opérateurs linéaires sur E.
Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Exemples d’équations linéaires

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Classi…cation des EDP
Considérons la forme générale d’une (EDP) de second Ordre
suivant les deux variables indépendanres (x; y)
@2 @2 @2 @ @
A 2
+ B + C 2
+D +E +F =f
@x @x@y @y @x @y

Une classi…cation assez simple de cette Equation peut etre faite


sur la base des coe¢ cients A, B et C.
= B2 4AC

l’Equation et dite de type :


1 Elliptique si <0
2 Parabolique si =0
3 Hyperbolique si >0
Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
EDP elliptique : Laplace
Cette catégorie d’EDP est associée aux problèmes de nature
stationnaire ou d’équilibre tels que l’écoulement stationnaire
d’un ‡uide visqueux, la répartition stationnaire du champ de
température ou la distribution d’un potentiel.

@2u @2u @2u


u= + 2 + 2 = f (u)
@x2 @y @z

La principale caractéristique de ce type d’équation elliptique est


qu’une perturbation introduite en un point quelconque à
l’intérieur du domaine de calcul in‡ue sur la totalité du
domaine.
Ceci implique que pour résoudre un problème de type elliptique
il est impératif de poser les conditions aux limites sur toutes les
frontières du domaine.
Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
EDP parabolique : Transport
C’est le cas d’un problème de propagation associé à un
mécanisme de dissipation tel que la conduction thermique
(Problème de Di¤usion).

@u @2u
= f (u)
@t @x2

C’est une équation linéaire d’ordre 2, parabolique par rapport à


la variable du temps t.
La propagation en avant dans le temps et la di¤usion dans
l’espace, font que la solution en un point P peut in‡uencer
n’importe qu’elle point pour t ti

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
EDP hyperbolique : Propagation des
Ondes
Cette catégorie d’EDP peut être considérée comme extension
des équations elliptiques pour lesquels certaines valeurs
critiques des paramètres doivent être déterminées en même
temps que la distribution d’équilibre correspondante.

@2u @2u
= f (u)
@t2 @x2

Exemples : L’équation de propagation d’une onde , la résonance


de circuit électrique ou d’enceintes acoustiques ainsi que la
détermination des fréquences propres des structures élastiques
constituent des exemples de ce type d’équations.

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Conditions aux limites :
Soit un problème dans un domaine , limité par la frontière @ .
Les conditions aux li,ites peuvent etre de trois natures :
1 Type Dirichlet : les conditions de la valeur de lq vqriqble

dépendant est imposée sur la frontière du domaine de calcul


=f sur @

2 Type Newman : La varibale dépendant n’est pas connue


sur la frontière mais sa dérivée est bien dé…nit
@
= f sur @
@n

3 Type mixte : une combinaison des deux conditions est


imposée sur la frontière
@
+ k = f ; k > 0 sur @
@n
Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Examples 8

1- Déformation d’un …l élastique


u00 (x) + c(x)u(x) = f (x) 8x 2 ]0; 1[
u(0) = u(1) = 0

Condition Direchlet

2- Vibration d’une corde


8 9
> @ 2 u(x;t) @ 2 u(x;t) >
>
> @t2 @x2
= f (t; x) 8t 2 R + ; 8x 2 ]0; 1[ >
>
< =
u(0; t) = u(1; t) = 0 8t 2 R +
>
> u(x; 0) = u0 (x) 8x 2 ]0; 1[ >
>
>
: @u(x;0) >
;
@t = u1 (x) 8x 2 ]0; 1[

Conditions Mixte
Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Examples 9

Donner le type de chacune des équations suivantes :

@2u @2u
=F
@x2 @y 2
@2u @2u
c =F
@t2 @x2
@2u @2u @ 2 u @u
+ 4 + 2 + =F
@x2 @x@y @y 2 @x

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Résolution des EDP à coe¢ cients
constants
Une équation di¤érentielle linéaire du second ordre à coe¢ cients
constants est du type :
ay 00 (x) + by(x) + c(x) = 0 (Ec )

où les réels a; b et c sont donnés dans R :R:R


On appelle équation caractéristique associée à (Ec ) :
ar2 + br + c = 0

Notons = b2 ac son discriminant.


1- Si > 0 alors l’équation caractéristique admet deux
solutions r1 6= r2 réelles. La solution générale de (Ec ) s’écrit :

y(x) = C1 er1 x + C2 er2 x ; C1 ; C2 2 R

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
2- Si = 0 alors l’équation caractéristique admet une unique
b
solution r = 2a (racine double) et les solutions de (Ec )
s’écrivent :
y(x) = (A + Bx)erx

où A et B sont deux constantes arbitraires réelles.


3- Si < 0 alors l’équation caractéristique admet deux
_
solutions complexes conjuguées r = + i et r = i .
Comme précédemment on a donc les solutions complexes :
_
y(x) = C1 erx + C2 e r x ; C1 ; C2 2 C:

ou bien si on veut les solutions réelles de (Ec ) :

y(x) = e x (Acos( x) + Bsin( x)) ; A; B 2 R:

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Recherche d’une solution particulière
Soit l’ED complète :

ay 00 (x) + by(x) + c(x) = f (x) (EG )

On cherche une solution particulière par appliquation l’une des


méthodes suivante :
1 Si le second membre est de la forme
f (x) = cos(x) + sin(x) alors on cherche une solution de
la forme
y0 (x) = c1 cos(x) + c2 sin(x)

Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Recherche d’une solution particulière
1 Si le second membre est de la forme f (x) = e x Pn (x) avec
Pn (x) Polynome de degré n: Alors on cherche une solution
de la forme :
i 1ére cas : Si a 2 + b + c 6= 0 tel que 6= r1
et 6= r2 ; on cherche une solution sous la forme
y0 (x) = e x Qn (x)

ii 2éme cas : Si a 2 + b + c = 0 tel que = r1 ou


= r2 avec r1 6= r2 ; on cherche une solution sous
la forme
y0 (x) = e x xQn (x)

iii 3éme cas : Si = r1 = r2 ; on cherche une


solution sous la forme
y0 (x) = e x x2 Qn (x)
Y . A b o u e lh a n o u n e F iliè re G é n ie M é c a n iq u e
Exercices

Résoudre sur ℝ les équations différentielles suivantes:

𝑦 ! + 2𝑦 = 𝑥 " 𝐸#

𝑦 !! − 2𝑦 ! + 𝑦 = 0 (E2)

𝑦 !! + 𝑦 = 2 cos " 𝑥 (E3)


Solution

On a 𝑦 ! + 2𝑦 = 𝑥 " 𝐸#

On commence par résoudre l'équation homogène associée 𝑦 ! + 2𝑦 = 0

les solutions sont les 𝑦(𝑥) = 𝜆𝑒 "#$ , 𝜆 ∈ ℝ


Il suffit ensuite de trouver une solu1on par1culière de (E1). Le second membre étant
polynomial de degré 2, on cherche une solu1on par1culière de la même forme:
𝑦% (𝑥) = 𝑎𝑥 # + 𝑏𝑥 + 𝑐 est solution de 𝐸&
⟺ ∀𝑥 ∈ ℝ, 𝑦%! (𝑥) + 2𝑦% (𝑥) = 𝑥 #
⟺ ∀𝑥 ∈ ℝ, 2𝑎𝑥 # + (2𝑎 + 2𝑏)𝑥 + 𝑏 + 2𝑐 = 𝑥 #
Ainsi, en identifiant les coefficients, on voit que

& & &


𝑦% (𝑥) = # 𝑥 # − # 𝑥 + ' convient.

Les solu1ons de (E1) sont


Exercices

𝑦 !! − 2𝑦 ! + 𝑦 = 0 (E2)

𝑦 !! + 𝑦 = 2 cos " 𝑥 (E3)

Vous aimerez peut-être aussi