E.D.
P linéaires du second ordre : Généralisations
Définition et classification : Equation générale
z est une fonction inconnue des n variables (x1,x2,…,xn) = x,
les coefficients Aij et Bj étant, comme C et F, des fonctions
données de ces n variables sur un domaine Ω dans .
Notons que
On pose M matrice carré d'ordre n
Et on définit la forme quadratique Q à n variables (f1,f2,…,fn) = f par
Trois cas se présentent en un point x:
1) La matrice M est régulière et Q garde un signe constant quelque soit f
vecteur de IRn non nul. Dans ce cas Q est décomposable en une somme de n
carrés. L'E.D.P est dite elliptique dans Ω.
2) La matrice M est régulière, mais le signe de Q n'est pas constant lorsque f ϵ
IRn. Q est alors décomposable en une somme algébrique de n carrés. L'E.D.P
est dite hyperbolique dans Ω.
3) La matrice M est singulière, donc Q est la somme algébrique de (n-1)
carrés. L'E.D.P est dite parabolique dans Ω.
Exemples :
a) Equation de Laplace dans IR3: x=(x1,x2,x3)
La matrice M est donnée par
Notons que les coefficients Bi, C et F sont tous nuls.
Comme le déterminant de cette matrice est non nul alors
elle est régulière. Calculons donc la forme quadratique.
f=(f1,f2,f3) étant un vecteur de IR3 non nul
Q est positive quelque soit f non nul, donc l'équation de
Laplace est elliptique.
2) Equation des ondes acoustiques
L'équation des ondes acoustiques est celle qui régit la
propagation des vibrations de faible amplitude dans un
fluide compressible dont la loi d'état est donné par P =
f(ρ) où P est la pression, ρ étant la densité volumique.
Cette équation est de la forme, pour x=(x1,x2,x3,t)
c représente la vitesse du son à travers le fluide tandis
que z représente le potentiel de vitesse c'est à dire le
vecteur vitesse est égale au gradient de z.
La matrice M est donnée par
Les coefficients Bi, C et F Sont tous nuls. C'est une matrice
régulière son déterminant est (-1/c2). Calculons à présent la
forme quadratique, f= (f1,f2,f3,f4) étant un vecteur de IR4, on a
Q change de signe, donc I'E.D.P est hyperbolique.
3) Equation de la diffusion de la chaleur
soit z (x1,x2,x3,t) la température d'un solide ou
d'un fluide dans un domaine Ω de IR3.
L'équation est donné par
k est appelé la diffusivité thermique.
La matrice M est donnée par
son déterminant est nul, la matrice est donc singulière.
Quelque soit f= (f1,f2,f3,f4) de IR4 non nul, la forme
quadratique est égale
Q est la somme de (n-1) carrés, n= 4. L'équation est
donc parabolique.
Remarque importante :
Revenons maintenant à l'équation à deux
variables
La matrice M est donnée par
a) Si S2-RT est non nul, l'équation est hyperbolique ou
elliptique. On calcule alors la forme quadratique, à
savoir Q = tf M f avec f=(f1,f2)
Il apparaît que si
S2-RT < 0 Q a un signe constant (celui de R), alors
l'équation est elliptique.
S2-RT > 0 Q change de signe, alors l'équation est
hyperbolique.
b) Enfin si S2-RT est nul, l'équation est parabolique.
Autre méthode de classification
La classification de l'équation envisagée se fait aussi en travaillant
sur les valeurs propres de la matrice M. M étant symétrique et à
éléments réels, ses valeurs propres λi sont toutes réelles.
a) Si les valeurs propres λi sont strictement positives ou bien
strictement négatives, quelque soit i, alors l'équation est de type
elliptique.
Voir exemple de l'équation de Laplace dans IR3, les valeurs
propres sont 1, 1, 1, toutes strictement positives.
b) Si au moins une valeur propres est de signe
opposé au signe d'au moins une autre valeur propre
alors l'équation est de type hyperbolique.
Voir exemple de l'équation des ondes acoustiques,
les valeurs propres sont 1, 1, 1, (-1/c2) .
c) Si une valeur propre au moins est nulle alors
l'équation est de type parabolique.
Voir exemple de l’équation de la diffusion de la
chaleur, les valeurs propres sont 1, 1, 1, 0.
Problèmes aux conditions initiales et aux conditions aux limites
En général, on travaille sur des équations aux dérivées partielles linéaires du
second ordre à quatre variables indépendantes, trois d'espace et une de
temps, situation qui se présente le plus souvent en physique mathématique.
Comme on vient de le constater sur des exemples, une E.D.P admet en
général une infinité de solutions dans lesquelles interviennent des fonctions
arbitraires, donc l'unicité de la solution n'est pas satisfaite.
La question donc ne consiste pas à trouver une solution quelconque u(x,y,z,t)
de l'équation, mais à choisir parmi toutes les solutions possibles, une d'entre
elles qui satisfait à des conditions accessoires convenablement données.
Ainsi u(x,y,z,t) doit vérifier I'E.D.P dans un certain domaine fermé Ω et
satisfaire sur la frontière de Ω c'est à dire sur ∂Ω à certaines conditions.
On rencontre le plus généralement des conditions
accessoires de deux espèces.
Les conditions initiales
Ces dernières marquent qu'à l’instant initial, il existe dans
un certain domaine Ω certaines relations entre les
variables de l’espace, la fonction u(x,y,z,t) ou une ou
plusieurs de ses dérivées partielles premières.
Les conditions aux frontières ou aux limites
Elles expriment que, quel que soit l'instant t, certaines
conditions sont satisfaites sur la frontière du domaine Ω:
∂Ω, par exemple qu'il existe une relation donnée entre x,
y, z; u(x,y,z,t), u’t, …
Considérons l'équation des ondes à deux variables
a étant une constante. On peut la considérer comme définissant
l'élongation u(x,t) au point d'abscisse x, au temps t, d'une corde AB, de
longueur L, fixée à ses extrémités A et B.
On peut se donner les valeurs,
pour t=0, de u et de u’t avec 0 ≤ x ≤ L.
De plus ajouter à ces conditions les relations qui expriment que la
corde est fixée à ses extrémités,
u(0,t) = 0 et u(L,t) = 0; t > 0.
On se donne ainsi les conditions initiales et les conditions aux limites.
Exemples de problèmes aux conditions initiales et aux
conditions aux limites
De façon générale, on rencontre l’un des problèmes suivants:
a) Problème de Dirichlet
Il consiste à résoudre une E.D.P de type elliptique définie
sur le domaine Ω. La valeur de la solution étant imposée
sur la frontière du domaine. C'est à dire
L étant un opérateur différentiel linéaire du second ordre
dans Ω.
b) Problème de Neuman
Il s'agit du même type de problème, mais cette
fois, c'est la dérivée normale de la solution qui
est imposé sur la frontière
N.B:
c) Problème Mixte
Il existe des problèmes de type mixte, où par exemple
les conditions aux limites sont de
Dirichlet sur une partie de la frontière ∂Ω1 et de type
Neumann sur la partie
complémentaire ∂Ω2 : ∂Ω = ∂Ω1 ∂Ω2. u vérifie donc le
système
d) Problème de Cauchy
Ce problème consiste à déterminer une solution
u(x1,x2,x3,x4) d'une E.D.P dans un
domaine Ω satisfaisant aux conditions suivantes,
dans lesquelles a désigne une valeur
d'une des variables xi (qui peut être le temps)
e) Problème extérieur
Il se pose lorsque l’on cherche à définir une
solution d'une E.D.P dans un domaine
extérieur à un obstacle δ (écoulement autour
d'un fuselage par exemple) il faut alors
outres les conditions aux limites sur ∂δ, des
conditions à l’infini.
