Chapitre I

Equations de la physique mathmatique

I.1 Classification des quations
I.1.1 Notions gnrales
La relation reliant des variables indpendantes dune fonction inconnue et ses drives
partielles sappelle quation aux drives partielles :
(

la fonction u(x,y) est continue dans un domaine D, ainsi que ses drives. En
transformant cette quation en identit, elle sappelle la solution rgulire ou classique de
cette quation. On va considrer lquation de deuxime ordre suivante :
(

Cette quation est quasi-linaire si A, B, C dpendent de x, y, u, ux, uy et linaire si A, B,

C dpendent seulement de x et y avec F- est une fonction linaire par rapport u, ux, uy.
Soit lquation linaire de deuxime ordre sous sa forme gnrale :

o A, B, C,, F, - des fonctions de x et y. Lquation est linaire et non homogne si

et linaire homogne si
I.1.2 Equation linaire homogne de premier ordre
Soient n variables
indpendantes x=(x1,xn). Considrons lquation linaire
homogne aux drives partielles de premier ordre :

Admettons que
alors :

est continue, diffrentiable dans le domaine D de lespace E n,

Le systme dquations diffrentielles ordinaires est :

On appelle

la solution (courbe intgrale) du systme (I.2). Si la fonction

(

est constante

sur chaque courbe intgrale du systme (I.2), alors lquation

premire intgrale de ce systme.

est la

Thorme. La fonction
sera la solution de lquation linaire
homogne aux drives partielles si et seulement si
est lintgrale des
quations diffrentielles ordinaires du systme correspondant.
Dmonstration.
La fonction

est la solution de lquation (I.1) signifie que :

Par ailleurs, la fonction

dpend seulement de t :
(

par consquent

- solution du systme (I.2), do

mais

Si

- la premire intgrale du systme (I.2) ; sur la courbe intgrale

)
, on a (

par consquent

mais

- la solution du systme (I.2), ce qui signifie

c'est--dire que

est la solution de (I.1)

I.1.3 Types dquations de second ordre

Les courbes de deuxime ordre

sont de trois types (hyperbolique, parabolique et elliptique) selon le signe du discriminant

de lquation
Ainsi, pour une quation de type :

qui au point (x, y) sera hyperbolique si

et elliptique si

, parabolique si

Si les coefficients A, B, C de (I.4) dans le domaine D sont constants, (x, y) = const et

le type dquation sera le mme tous les points du domaine.
Si dpend de (x, y), on pourra avoir (x, y) > 0 et lquation sera hyperbolique, alors
que dans lautre partie
elle sera elliptique et sur la courbe (x, y) = 0, se produit
une dgnration du type dquation. Alors, la courbe (x, y) = 0 sappelle ligne de
dgnrescence parabolique.
Ainsi, lquation

a le descriminant
, c'est--dire que lquation sera de type
hyperbolique si 13y > 0, elliptique si 13y < 0 et sur la ligne (1-3y) se produira une
dgnrescence du type parabolique.
Dans le cas de n variables indpendantes (n>2) :
lquation quasi-linaire est :

, la forme gnrale de

ici
- sont fonctions de x1, , xn.
Au point fixe (x1, , xn) soit la forme quadratique :

On peut passer, laide de la transformation linaire de (I.6), la forme canonique :

2
j

. Lquation (I.5) sera alors de type elliptique au point (x1, , xn) si

, j=1, , n , hyperbolique si
pour j = 1, , m < n et
pour j = m+1,
,n et parabolique si au moins un des coefficients
.
, o

I.2 Transformations des quations de deuxime ordre

I.2.1 Invariance du type dquation
Soit dans le domaine D, lquation
(

qui est linaire par rapport aux drives dordre suprieur. Admettons, que les coefficients A,
B, C sont continus dans le domaine D et montrons que le type dquation est invariant par
rapport la transformation des variables indpendantes.
Introduisons de nouvelles variables :

Admettons que (I.8) est rsolvable aux alentours du point

|
Si

, ce qui signifie que :

est compose, alors :

En remplaant les drives dans (I.7), on peut crire :

(

En calculant le discriminant de (I.11) :

(
(

(
)(

)
)

(
(
(

)
(

((

)
)

)
(

Alors :
(

Donc, si A, B, C sont continus aux alentours du point (x, y), alors que la transformation (I.8)
nest pas dgnre, donc les signes des discriminants et 1 aux alentours correspondant
aux points (x, y) et
sont les mmes et par consquent le type dquation est conserv.
I.2.2 Forme canonique
de faon simplifier lquation aux
pour que dans lquation (I.11) les

On peut choisir les fonctions

drives partielles. Alors, choisissons
coefficients soient nuls pour
:

c'est--dire que les fonctions

soient des solutions de :

(

Lquation diffrentielle ordinaire

correspondant (I.14) est la caractristique et ses intgrales sont les quations caractristiques
aux drives partielles (I.7). Supposons que A0 et rsolvons (I.15) par rapport y' :

Considrons trois cas :

Type hyperbolique. Le discriminant
(I.16) sont diffrentes, do les intgrales :

. Les parties droites de

de ces quations qui sont indpendantes. Posons les nouvelles variables indpendantes
suivantes :

Alors, on peut crire les caractristiques (I.17) de lquation (I.7) sous la forme :
(

En divisant lquation par 2B1, on obtient :

(

qui est la premire forme canonique de lquation hyperbolique.

Si on pose :

alors, lquation (I.18) sous la deuxime forme canonique de lquation hyperbolique scrit :
(

Type elliptique. Le discriminant

lquation (I.16) sont complexes et conjugues, cest pourquoi :

, les parties droite de

qui sera complexe et conjugue.

Si lon passe de nouvelles variables :

alors, on a :
(
o

et

sont des variables complexes. Si lon passe de nouvelles variables

alors,
(
et lquation (I.20) scrit sous la forme :

qui est la forme canonique de lquation elliptique.

Type parabolique. Le discriminant
et on a un ensemble de
caractristiques
. En introduisant de nouvelles variables, alors :

o
- les caractristiques de lquation (I.7) et
quelconque ne dependant pas de
.

une fonction diffrentielle

est la solution de (I.14), do :

La fonction

Le type dquations ne change pas pendant la transformation et par consquent,

. Cependant,
, signifie que
Montrons, que C 0 aux alentours du point (x, y). Admettons linverse, soit

A condition que les coefficients A, B, C ne sannulent pas en mme temps, cest pourquoi si
lon pose A = 0, C = 0, alors de
, A = B = C = 0, donc on a soit A 0
soit
.
Si A > 0 aux alentours du point (x, y), alors de
. Alors,

et on peut crire

Analogiquement, de

, on trouve

Par consquent,
sont dpendantes, ce qui contredit la condition de choix de la
fonction
et signifie, que
aux alentours du point (x, y). Donc, on obtient enfin
, et la transformation de lquation (I.7), aprs division par C1, scrit :
(
qui sappelle la forme canonique de lquation parabolique.

Donc, lquation hyperbolique

indpendantes

a deux familles de caractristiques

.

Lquation elliptique
de nouvelles variables
sous la forme canonique.

ne possde pas de caractristiques. Si lon passe

, alors lquation se transformera

Lquation parabolique ( (x, y) = 0) a une seule famille de caractristiques

. Si lon passe de nouvelles variables
, o ne dpend pas de
, alors, lquation se transformera sous la forme canonique.
Si lon a plus de deux variables indpendantes, on peut transformer chaque point de
lquation sous la forme canonique, et donc on a :
le type hyperbolique :

le type parabolique :

le type elliptique :

Afin de simplifier davantage lquation, on peut introduire une nouvelle fonction

inconnue. Soit lquation hyperbolique linaire sous la forme canonique :

On admet dans le domaine D, que les coefficients

diffrentiables. Introduisons une nouvelle fonction (x, y) :

restent pour le moment quelconques. De (I.24), on trouve

En substituant les drives dans (I.123), on a :

sont continus et

tant arbitraires, on peut annuler un des coefficients dans (I.25), et par consquent
simplifier davantage lquation.
Thorme 1. La forme canonique de lquation hyperbolique peut scrire, en posant
, sous la forme :

si et seulement si
Dmonstration.
Soit lquation (I.23) transforme sous forme de (I.26), alors :

En diffrentiant

par rapport x et b1 par rapport y, on a :

En remplaant dans ces galits

on obtient enfin
Soit

partir de (I.27), on a :

. La fonction

est choisie de sorte que dans (I.25) le coefficient :

- une fonction quelconque

En diffrentiant (I.27) par rapport x et en utilisant

Choisissons

afin que

Donc, de (I.29), on a :

, on obtient :

En remplaant

dans (I.28), on a

c'est--dire que lquation (I.23) se ramne la

Alors, les coefficients seront

forme (I.26).
Consquence.
Si les coefficients de lquation

satisfont les conditions

alors, on peut crire

Dmonstration.
En prenant

selon la formule (I.30), on aura :

do

((
(

Ainsi, on suit la mme procdure pour lquation elliptique

Thorme 2. Lquation canonique de type elliptique, en posant

scrire, sous la forme :

si et seulement si
sous la forme

. Si lon pose
.

, peut

, alors on peut crire aussi lquation

Si les coefficients sont constants dans lquation de deuxime ordre donne sous forme
canonique, alors la condition
est accomplie. En introduisant une nouvelle fonction
, o et sont constantes, alors :

En introduisant de nouvelles variables indpendantes :

on aura les coefficients

- hyperbolique
- elliptique
- parabolique
gaux un.

I.2.3 Solution gnrale

La solution gnrale dune quation diffrentielle ordinaire de deuxime ordre dpend
de deux constantes arbitraires. La notion de solution gnrale pour une quation aux drives
partielles peut se ramener une solution dpendant de deux fonctions arbitraires.
Exemple. Trouver la solution gnrale de :

Solution.
Transformons cette quation sous la forme canonique
(

Si
, lquation sera hyperbolique et si
dgnre.
Lquation caractristique est :

En intgrant pour

)
, elle sera paraboliquement

, on a

ce qui reprsente une famille de droites et une famille dhyperboles.

En introduisant les variables caractristiques

En substituant les drives dans lquation de dpart, on a :

, on trouve

, on a lquation aux variables caractristiques :

En simplifiant par

En considrant le facteur entre parenthse comme une nouvelle fonction inconnue, en

intgrant suivant , on trouve :

o
est une fonction arbitraire de . Lquation trouve peut tre considre comme une
fonction diffrentielle linaire ordinaire de premier ordre ( - une variable indpendante, fixe). En lintgrant, on a la solution gnrale:

est arbitraire. Dsignons par

et revenons aux anciennes variables,

on aura la solution gnrale de lquation initiale :

qui dpend de

et

I.3 Equations principales de la physique mathmatique

Soit le point x = (x1, , xn) dans un espace multidimensionnel et soit t le temps. u(x, t)
est la fonction cherche. Les principales quations de la physique mathmatique sont les :
- Equations doscillations
x

2 u x, t n
u x, t
p x

q x u x, t F x,t ,
2
xk
t
k 1 x k

avec x , px , qx - les coefficients du milieu, F x, t - caractrise lintensit de laction

(force) extrieure.
Si

x , px const , qx 0, alors

F x,t
p
a2 ,
f x,t ,

et on obtient lquation donde dans lespace En:

2
2u
2u
2 u

f x,t ;

...

x 2
t 2
xn2

Introduisons loprateur de Laplace :

2
2

...

,
x12
x n2

lquation donde devient alors :

2u
a 2 u f x,t .
2
t
Les processus doscillations sont dcrits par lquation hyperbolique.
- Equation de la chaleur
La conductibilit thermique et la diffusion sont exprimes par lquation parabolique :
x

u x,t n

t
k 1 x k

u x,t
px
qx u x,t F x,t .
xk

x , px const , qx 0,

si, en particulier
alors, on aura :

u
a 2 u f x, t .
t

- Equations stationnaires
Si les processus dcrits par lquation doscillations ou de la chaleur ne dpendent pas
du temps, c'est--dire sils sont stationnaires, alors on passe lquation de type elliptique
suivante :
n

u x
px
q(x)u(x) F(x)
xk
k

x
k 1

et en particulier lquation de Poisson, si p(x) = const, q (x) = 0,

2u
f(x);
2
k 1 x k
n

u x

Si laction extrieure est absente f x 0, alors on obtient lquation de Laplace :

2u
0,
2
k 1 x k
n

u x

- Equation de Helmholtz
Pour chaque quation de type elliptique, on cherche la solution de lquation donde :
utt a 2 u

sous forme dun produit ux, t ( x)T (t ) . Alors, en la remplaant dans lquation, on
obtient :
( x) T " (t ) T (t ) a 2 ( x);

En divisant par : a 2 (x)T (t) , on obtient lgalit

(x) T" (t)

2
,
(x)
a T(t)

qui nexistera que si le rapport des termes constants situs la partie gauche ne dpend pas de
t, et celui de droite ne dpend pas de x. Si lon galise ces rapports pour T (t) et (x),
alors :
T " (t ) a 2 T (t ) 0,

(x) (x) 0.
Lquation elliptique pour la dtermination de (x) sappelle quation de Helmoholtz.
Considrons maintenant la dduction dquations et des conditions aux limites pour
certains cas concrets.
Vibrations dune corde
On considre les quations de vibrations transversales dune corde fine, flexible et
lastique.
Le principe de base est celui de dAlembert : la force rsultante, y compris les forces
dinerties agissant sur le systme, est nulle.
Considrons que la corde, la position dquilibre, soit tendue selon laxe ox et soit u
(x, t) llongation dun point de cette corde au temps t. u = u (x, t) est alors lquation
dcrivant la forme de la corde linstant t.
Soient x la densit de la corde, F (x, t) - la densit des forces exterieures agissant
u
suivant la direction
, F x - la densit des forces de rsistance- proportionnelle la
t
vitesse ut. Soient langle dinclinaison de la corde par rapport
au point x et ds la
diffrentielle de la longueur de larc de la corde. Alors, approximativement, on a :

ds 1 u x2 dx 1 u x2 ...dx dx ds dx ;
2

1
1
tg x u x , cos x
1 u x2 ... 1 ;
2
2
1 tg x

sin x

tgx
1 tg 2 x

u x 1 u x2 ... u x ;
2

1 u x2
ux

ici, on dcompose 1 u x2

1
2

en srie de Taylor.
u

X
Fig I.1

Imaginons quune partie de la corde de projections (x1, x2) et de tensions

correspondantes T(x1, t), T(x2, t) est dirige tangentiellement par rapport la corde. La
longueur de cette partie considre est :

S x1 , x2

x2

x2

1 u dx dx x2 x1 ,
2
x

x1

x1

La projection sur ox de la force agissant sur le morceau (x1, x2), selon le principe de
dAlembert, est :

T x1 cos x1 Tx2 cos x2 0 Tx1 cos x1 Tx2 cos x2 .

mais cos x1 cos x2 1 , alors, Tx1 Tx2 , c'est--dire que la tension T x, t T ne

dpend pas de x. Ainsi, on a T x, t T , o T est la longueur du vecteur tension de la corde

ltat de repos.
La projection des forces de tension sur ox sera :

T cos x1 T cos x 2 T sin x1 sin x 2 T u x x1 , t u x x 2 , t

2

x2

Tu xx x, t dx.
x1

Les forces exterieures, les rsistances et linertie agissent travers ox et la projection de

leur rsultante sur ox sera alors :
x2

Tu

x1

xx

x2

x2

x2

x1

x1

x1

dx F x,t dx k x ut dx x utt dx 0,

En admettant, que la fonction sous intgrale est continue et si lon applique le thorme
de la moyenne et quon divise par (x2 x1) et en passant par la suite la limite
x1 x, x 2 x, on obtiendra lquation de vibrations de la corde sous forme de :
Tu xx F x,t k x ut x utt 0.

(I.31)

Posons :

F x,t
T
k x
a2 ;
k1 ;
f x,t ,
x
x
x
rcrivons lquation de vibrations de la corde dans un milieu rsistant :
utt k1ut a 2 u xx f x,t .

(I.32)

Si, en particulier, k (x) = 0, on obtiendra lquation donde dans un espace homogne E 1 pour
un milieu non rsistant :
utt a 2 u xx f x,t ,

(I.33)

avec f(x, t) - la densit des forces externes par unit de masse de la corde, si x const ,
alors, a 2 const .
Lquation donde homogne f x, t 0 dans E1 :
utt a 2 u xx .

(I.34)

dcrit les vibrations libres de la corde.

Conditions aux limites
Le caractre de fixation des extrmits de la corde, c'est--dire les conditions aux
limites, influe sur sa vibration, autrement dit, sur lquation (1.32).
On considre, quun point dune corde est fix dans une direction si ce point ne bouge
pas dans la direction ou dans la direction inverse. Il est mou- si son dplacement est libre, et
lastique- si le dplacement dune variable dans cette direction produit dans le nud de
fixation une raction :

R k1, k1 0.
Si
on a une fixation rigide, et si
dans ce qui suit la fixation des extrmits

alors on a une fixation moue. On considre

=l le long de laxe ox - rigide.

Soit une partie de la corde 0 x x1 x 2 x l , contenant lune de ses extrmits

(fig.I.2). Cherchons les conditions dajustement des nuds de fixation sur la vibration de la
corde. La raction dans un nud de fixation est
(la raction
. La tension dans lautre extrmit est T. La projection des forces agissant
sur la partie 0, x1 (sur la partie x2 ,l ) sera aux alentours du point x = 0 :

x1
x1
x1

k1 u 0,t T x1 cos T x1 , Ou F dx k ut dx utt 0 ;

0
0
0

Aux alentours du point x = l, elle sera :

l
l
l

k1u l, t T x2 cos T x2 , Ou F dx k ut dx utt dx 0.

x2
x2
x2

Mais

(T x1 , Ou ) x1 cos T x1 , Ou sin x1 u x x1 ,t .
2

(T x2 , Ou) x2 cos T x2 , Ou sin x2 u x x2 ,t .

2

En remplaant dans cette galit les cosinus et en passant la limite x1 0 pour x 2 l ,

on obtient la condition aux limites au point x = 0 (au point x = l).
k1 u 0, t T u x 0, t 0 ;
k1 u l , t T u x l , t 0.

0
R(0)

R(l )
Fig.I.2

Conditions aux limites de premier type.

Posons k1 , on obtient les conditions de fixation rigide des extrmits :

u0,t 0 ; ul, t 0.

(I.35)

Si les extrmits se dplacent suivant les lois de 1 (t ), 2 t , alors, les conditions aux
limites seront :

u0, t 1 t ;
Conditions aux frontires de deuxime type

ul , t 2 t .

Posons k1 = 0, nous trouvons les conditions dune fixation moue (les extrmits sont
libres suivant ox) :
ux (0, t) = 0 ;

ux (l, t) = 0,

(I.36)

cela signifie gomtriquement, que langle dinclinaison de la corde par rapport ox aux
points de fixation est nul.
Condition aux frontires de troisime type
Les extrmits sont fixes lastiquement ( k1 0, k1 ). Dsignons par :
pour lextrmit gauche de la corde et

k1
h1 , T

k1
h2 -pour son extrmit droite. Alors, les
T

conditions aux limites pour une fixation lastique sont :

u x 0, t h1 u(0, t ) 0 ; u x (l , t ) h2 u(l , t ) 0.

(I.37)

Ces trois conditions aux limites peuvent tre crites sous la forme :

1 u x o, t 1u0, t 1 t ; 2 u x l , t 2 ul , t 2 (t ),
o

1 ,..., 2 const ; 12 12 0 ; 22 22 0.

Vibrations dune tige

Lquation donde dans E1 :

u tt a 2 u xx f x, t ;

a2

dcrit aussi des petites vibrations longitudinales dune tige lastique. Ici nous avons est la
densit volumique de la tige, E le module dlasticit du matriau, f(x, t) la densit des
forces exterieures, u(x, t) le dplacement suivant ox. Alors :

u ( x x, t ) u ( x, t )
u
lim
x x 0
x

0
Fig.I.3

est le dplacement relatif au point x. Selon la loi de Hooke, la contrainte au point x est
u
T ES , o E est le module dYoung, S - la section. Les conditions aux limites pour la
x
tige de longueur l sont :
Pour les extrmits x = 0,

x = l fixes rigidement :

u (0, t ) 0;

u (l , t ) 0,

Pour les extrmits libres :

u x (0, t ) 0;

u x (l , t ) 0.

Si les extrmits x =0, x = l sont libres, alors les contraintes aux sections x = 0, x = l
sont nulles suivant la loi de Hooke. Pour les extrmits fixes lastiquement, on a :
u x (0, t ) h1u(0, t ) 0; u x (l , t ) h2 (l , t ) 0

Vibrations dune membrane

Soit une membrane dans un tat de contraintes nulles situe dans un plan ox1x2, et soit
u( x1 , x2 , t ) linclinaison du point ( x1 , x2 ) de la membrane de sa position dquilibre. Les
vibrations traversales faibles de la membrane sont dcrites par lquation donde non
homogne dans E 2 ( deux variables gomtriques) :

u tt a 2 (u x1 x1 u x2 x2 ) f ( x1 , x 2 , t ); a 2

o la densit superficielle de la membrane, T- la tension calcule par unit de longueur,

f ( x1 , x2 , t ) - la densit des forces externes agissant perpendiculairement par rapport la

position dquilibre et calcule par unit de masse.
Vibrations dun milieu tridimensionnel
Lquation donde non homogne dans E 3 est :
utt a(u x1x1 u x2 x2 u x3 x3 ) f ( x1 , x2 , x3 , t )

Elle dcrit les petites vibrations lastiques des corps solides, les vibrations de gaz, du son, et
les oscillations lectromagntiques.
Equation de la chaleur
Soit un corps de densit ( x), x ( x1 , x2 , x3 ), de permabilit thermique c(x) et de
conductibilit thermique interieure k (x). On considre que le corps est isotrope. Soit u ( x, t ) la temprature au point x linstant t. Selon le principe de conservation de lnergie et la loi
de Fourier, la quantit de chaleur dQ passant pendant dt travers la surface d suivant le
vecteur de la normale v est :

dQ k(x)

u
d dt,
v

k(x) 0.

La quantit de chaleur passant suivant la direction v sera ngative (ou positive) si la

temprature u ( x, t ) diminue (ou augmente) suivant cette direction, c'est--dire que
u
u
0( 0) . En prsence dune source de chaleur de densit F(x, t) calcule par unit de
v
v
volume dans un volume D limite par la surface , nous supposons une quantit de chaleur
Q1 dgage par cette source pendant lintervalle de temps t 2-t1, Q2 - la chaleur servant
augmenter la temprature du corps et Q3- la chaleur passant par la surface suivant la
normale extrieure v pendant t2- t1, alors :
Q1 Q2 Q3

La source dgage pendant dt dans llment du corps de volume dx dx1 dx2 dx3 la
chaleur dQ F (x , t ) dx dt . Alors, lintgrale suivant le volume D de t 1 t2 donne :
t2

Q1 dt F (x , t ) dx .
t1

Selon la loi de Fourier, la quantit de chaleur Q2 passant travers la surface suivant

la normale extrieure v sera :
t2

Q2 dt k (x )
t1

u (x , t )
d.
v

La quantit de chaleur d Q3 servant augmenter la temprature de llment dx de

u ( x, t1 ) u ( x, t2 ) est :

u (x , t )
dt ;
t
t1

t2

d Q3 ( (x )dx ).c (x ).(u (x , t 2 ) u (x , t 1 )) c (x ) (x )dx

Lintgrale suivant tout le volume est :
t2

Q3 dt c (x ) (x )
t1

u (x , t )
dx .
t

Remplaons Q1 , Q2 , Q3 dans lgalit Q1 Q2 Q3 ,

alors, on trouve ;
t2

t2

dt F ( x, t )dx dt k ( x)
t1

t1

Ainsi, lintgrale par rapport donne :

2
u ( x, t )
u ( x, t )
d dt c( x) ( x)
dx.
v

t
t1
D

u
u
d k
cos(v , x j )d
j 1 x
v
j

Si lon pose dans la formule de Gauss-Ostrogradsky :

3

k Aj cos(v, x j )d
j 1

Aj

j 1

x j

dx

la relation :

u
x j

Aj k

j 1, 2,3,

alors,

3
u

u
d A j
(k
)dx
j

1
v

x
j
j
D

(I.38)

En remplaant lintgrale suivant par son expression dans (I.38), on obtient lquation
intgrale u (x , t ) :

u
u
(k
) c( x) ( x) dx 0
j 1 x
x j
t
j

t2

dt F ( x, t )
t1

En considrant lexpression sous intgrale continue en fonction de (x, t), et en utilisant le

thorme de la moyenne des intgrales, alors :

u
u
(k
) c( x) ( x) ) x xcp .(t2 t1 ) D 0
j 1 x
x j
t t tcp
j
3

F ( x, t )

En simplifiant par (t2 t1 ) D , o D est le volume de D, en passant la limite pour

t1 , t 2 t , et en concentrant le corps D au point x = (x1, x2, x3), alors on obtient lquation

diffrentielle de la chaleur pour la fonction u( x, t ) :
3

F ( x, t )

j 1

u
u
(k
) c( x) ( x) 0
x j x j
t

Si la densit ( x), la permabilit c(x) et la conductibilit thermiques k (x) sont

constantes et si lon pose :

k
a2
c
lquation de la chaleur scrira :

F ( x, t )
f ( x, t ),
c

u
2u 2u 2u
a 2 ( 2 2 2 ) f ( x1 , x2 , x3 , t ),
t
x1 x2 x3

(I.39)

Elle est parabolique.

Si la temprature u( x1 , x2 , x3 , t ) ne dpend pas de x3, on aura :
u
2u 2u
a 2 ( 2 2 ) f ( x1 , x2 , t ),
t
x1 x2

(I.40)

Si la temprature u( x1 , x2 , x3 , t ) ne dpend pas de x2 et x3, alors ;

u
2u
a 2 2 f ( x1 , t ),
t
x1

(I.41)

Lquation (I.39) peut exprimer le processus de diffusion en cas o :

a2

k
;
c

f (x , t )

F (x , t )
c

o k est le coefficient de diffusion, c- le coefficient de porosit, F(x, t) la densit des sources

de la substance de diffusion.
Conditions aux limites de premier type
La temprature sur la surface dun corps D est :
u( x, t ) ( x, t );

Conditions aux limites de deuxime ordre

La quantit de chaleur passant par la surface unitaire et par unit de temps (flux
thermique) est :

u ( x, t )
v

( x, t ),

avec v la normale extrieure la surface . Si le corps est isol thermiquement, alors :

u
0
v
Conditions aux limites de troisime type
Le milieu dans lequel se trouve le corps D a une temprature ( x, t ) ; lchange de
chaleur entre le corps et le milieu environnant se ralise par la loi de Newton : le flux de
chaleur Q travers la surface suivant la normale extrieure v est proportionnel la
diffrence de temprature :

Q (u( x, t ) ( x, t )),

x , ( 0),

o - le coefficient dchange thermique.

Ou encore selon la loi de Fourier Q = - k
k

u
, on a la condition aux limites :
v

u ( x, t )
(u ( x, t ) Q( x, t ) )
v

u ( x, t )
(hu ( x, t )) ( x, t ), x , (h 0)
v

(I.42)

Processus stationnaires
Il sagit ici dquations de type elliptique.
1. Si le processus de vibrations dcrit par lquation donde :
2
2u
2u 2u
2 u
a ( 2

) f1 ( x1 , x2 , x3 )
t 2
x1
x2 2 x32

u 2u
sarrte et ne dpend plus du temps, alors 2 0 et lquation donde se transformera en
t t
quation de Poisson :
2u
2u 2u
u ( x1 , x2 , x3 ) 2

f1 ( x1 , x2 , x3 ),
x1
x2 2 x32

Si f1 ( x1 , x2 , x3 ) 0, c'est--dire en absence de perturbations externes, alors lquation de

Laplace sera :
2u
2u 2u
u ( x1 , x2 , x3 ) 2

0
x1
x2 2 x32

2. Si la distribution de temprature est invariable par rapport au temps c'est--dire que

u
0, et lquation (I.39) se transforme en quation de Poisson u( x1 , x2 , x3 ) f1 ( x1 , x2 , x3 )
t
et par consquent, on na pas de source de chaleur et f1 ( x1 , x2 , x3 ) 0 , on obtiendra lquation
de Laplace :

u( x1 , x2 , x3 ) 0

3. La tendance stationnaire est parfaite, c'est--dire un liquide incompressible et non

visqueux se caractrise par la vitesses :

(x 1 , x 2 , x 3 ) grad u (x 1 , x 2 , x 3 ).

o u est le potentiel en cas dabsence de sources, alors : div (x 1 , x 2 , x 3 ) 0 et lquation de

Laplace est :
div grad u 0

u ( x1 , x2 , x3 ) 0

4. Dans un milieu homogne et conducteur en absence de sources, on a :

u 0.

I.4 Position des problmes de la physique mathmatique

Considrons un domaine D , dlimit par une surface dans un espace n- dimensions.
Soit une fonction f ( x), x ( x1,..., xn ), dfinie dans ce domaine. On comprend par la limite de
:
f ( x) x f ( x) ( x)

comme tant la limite :

1

lim f ( x) ( x)
1

x .

x x
1

xD

On dit que f (x) est continue dans le domaine D jusqu' sa frontire, si

1) f (x) - est continue dans D
2) il existe la limite continue :
1

lim f ( x) ( x)
1

x .

x x
1

xD

f ( x) x D
et la fonction F ( x)
( x) x

est continue dans D D

Alors, on crit dans ce cas, que f ( x) C ( D) .

Si f (x) est continue et diffrentiable S fois dans le domaine D, et elle a sur la frontire
du domaine des drives jusqu' lordre k S , alors on pourra crire f ( x) C S ( D) C R ( D) .
En particulier, si f ( x) C 2 ( D) C ( D) , donc, on aura une fonction continue dans le
domaine D jusqu' sa frontire et elle sera doublement diffrentiable.
La surface est dite lisse ou appartenant la classe C 01 , si aux alentours de chaque
0

point x , elle est reprsente par une quation ( x) 0 , et la fonction (x) - est
continue, diffrentiable et grad ( x) 0 . Alors, lorsque la fonction (x) est continue, ainsi

que ses drives jusqu' lordre S 0 , on dit que la surface appartient la classe C S et on
crit C S .
Si la surface est compose dun nombre limit de surfaces lisses, on lappelle surface
lisse par morceau.
On va considrer ci-aprs, que la surface est ferme et lespace E n est partag en
deux domaines : D - fini et contenu dans la sphre x R , o R- est fini, et D - infini,
contenant un point infiniment loign.
Si la surface est considre comme la frontire du domaine D, alors le vecteur de la
normale extrieure y la surface au point y , est un vecteur issu de ce domaine D.
Lensemble des points ( x, t ) de sorte que x D , 0 t T , est un cylindre C D 0,T .
La frontire du cylindre est forme par la base infrieure D 0 , la base suprieure D T

et la surface latrale 0, T . Ainsi, si n 2 , le domaine D (a, b) est un intervalle (a, b) ,

le cylindre C est un rectangle dans le plan O t . Si n 1, le domaine D est la surface dans
le plan Ox1 x2 , le cylindre C- est le domaine dans lespace trimensionnel Ox1 x2t limit par les
domaine D dans le plan t 0, t T , les axes parallles Ot et la direction (fig.I.4)
t

t
T

T
D
C

0 a

Fig.I.4

Les quations diffrentielles ont un nombre indtermin de solutions. Afin de dgager,

de lensemble, la solution la plus correcte qui dcrit un processus donn, on pose des
conditions complmentaires qui dpendent du type dquation.
Les conditions complmentaires seront tablies dans ce paragraphe pour :
Les quations de vibrations (hyperboliques) :

( x)

2u
div ( P( x) grad u ) q( x)u F ( x, t );
t 2

Lquation de la chaleur (parabolique) :

(I.43)

 ( x)

u
div ( P( x) grad u ) q( x)u F ( x, t );
t

(I.44)

Lquation stationnaire (elliptique) :

div ( ( x) grad u) q( x)u F ( x) 0;

(I.45)

ici
n

div ( ( x) grad u )

k 1

u
( ( x)
),
xk
xk

n 1,2,3.

Les coefficients de lquation suivant leur sens physique seront :

( x) 0, P( x) 0, q( x) 0, ( x),
q( x) C ( D), P( x) C1 ( D)

o D- est le domaine de variation des variables gomtriques dans les quations (I.44-I.45),
c'est--dire le domaine o se droule le processus.

I.5 Problmes pour les quations hyperboliques

On va formuler le problme dabord pour un espace E1 ( une variable gomtrique), et
par la suite, ltaler sur lespace n.

2u
u
( x) 2 ( P( x) ) q( x)u F ( x, t ),
t
x
x

(I.46)

Problme mixte. Dans un rectangle C (0, l ) (0, T ) , il faut trouver la solution

rgulire de lquation continue, diffrentiable jusqu' la frontire, satisfaisant la condition
aux limites :

u
u
1u ) x 0 1 (t ); ( 2
2u ) x l 2 (t )
x
x

(I.47)

et la condition initiale :

u ( x, t ) t 0 ( x);

u ( x, t )
( x).
t t 0

(I.48)

Pour que la position du problme ne soit pas contradictoire, on doit satisfaire les conditions de
lissage :
F ( x, t ) C (u ); ( x) C1 (a, B);

( x) C (a, b); 1 (t ), 2 (t ) C (0, T )

et les conditions complmentaires :

 1 (0) 1 (0) u1 (0), 2 (l ) 2 (l ) u2 (l ).

Les conditions aux limites dterminent le comportement de la fonction sur la frontire
du domaine quelque soit t. Les conditions initiales donnent la fonction cherche et la vitesse
de sa variation x (0, l ) au moment initial.
Le problme (I.46), (I.47), (I.48) sera homogne, si F(x, t ), u1 (t ), u 2 (t ) 0, et il sera
non homogne lorsque au moins lune de ces fonctions est diffrente de zro.
En relation avec les conditions de lissage, on dit que le problme mixte est de :
premier type, si :

k 0, k 1 : u(0, t ) 1 (t ); u(l , t ) 2 (t );
deuxime espce, si :

k 1,
u (0, t )
1 (t );
x

k 0 :
u (l , t )
2 (t );
x

troisime espce, si :

k 1, 1 h1, 2 h2 , hk 0 :

u (0, t )
u (l , t )
h1u (0, t ) 1 (t );
h2u (l , t ) 2 (t ).
x
x

On formule le problme mixte, pour lquation (I.43), comme suit :

Dans un cylindre C D 0, T, il faut trouver la solution de lquation (I.43), continue,
diffrentiable jusqu' la frontire, et satisfaisant la condition aux limites :
(

u ( x, t )
u ( x, t )) x ( x, t )
v

et la condition initiale :

u ( x, t ) t 0 ( x);

u ( x, t )
( x).
t t 0

Alors, on doit satisfaire les conditions dune surface lisse et de concordance :

F ( x, t ) C ( D);

( x) C ( D).

( x) C1 ( D);
u ( x, t ) C ( 0, T );

( x)
( x) u ( x, t ) t 0
v

Problme de Cauchy. A la diffrence du problme mixte, on propose ici, que

lquation (I.43) dcrit le processus dans tout lespace des variables gomtriques
En ( D En ) t 0, cest pourquoi nous navons pas de conditions aux limites suivant ces
variables.
Le problme de Cauchy est formul comme suit.
Dans

un

demi

espace

t 0

de

variables ( x1 ,..., xn , t ), on

cherche

la

fonction,

u( x, t ) C (t 0) C1 (t 0) satisfaisant lquation (I.43) et les conditions initiales :

2

u ( x, t )
( x);
t t 0

u ( x, t ) t 0 ( x);
(I.49)

supposons pour cela les conditions de lissage suivantes :

F ( x, t ) C (t 0);

( x) C1 ( En ); ( x) C ( En ).

Le problme de Cauchy admet une solution gnralise dans ce cas, lorsque les
conditions initiales ne sont pas donnes pour t 0, alors pour la surface lisse, on a :

: t ( x);

u ( x, t )
;
v

u ( x, t ) ;

on tabli ici la solution dans une partie de lespace appartenant la surface .

Problme de Goursat. Il existe dautres conditions complmentaires pour lquation
hyperbolique. Pour le cas du problme de Goursat, on donne des conditions complmentaires
sur les caractristiques. Formulons le problme de Goursat pour lquation :

2u
u
u
a b cu F ( x, t ),
xt
x
t
(I.50)
Il faut trouver dans un rectangle C (0, l ) (0, T ) (fig.I.5) la fonction u ( x, t ) continue
jusqu' la frontire et qui est une solution rgulire de lquation (I.50) et satisfaisant les
conditions complmentaires :

u ( x,0) ( x), 0 x l;
u (0, t ) (t ),0 0 t T
sur les caractristique x 0, t 0
t
T
C
0

l X

Fig.I.5

I.6 Problmes pour les quations paraboliques

Le problme mixte et le problme de Cauchy pour lquation parabolique
u n
u
( x)

( p( x)
) q( x)u F ( x, t )
k

1
t
xk
xk
diffrent du problme pour lquation hyperbolique par le fait davoir une seule condition
initiale (y-compris la drive premire par rapport t) :
u( x, t ) t 0 ( x) .

De plus, il faut le limiter linfini, pour lunicit de la solution :

y( x, t ) A

x En

I.7 Problmes pour les quations elliptiques

Dans lquation elliptique :
n

k 1

u ( x)
( p ( x)
) q( x)u ( x) F ( x)
xk
xk

(I.51)

le temps t nest pas inclu et le domaine D de variation des variables gomtriques

x ( x1 ,..., xn ) est un domaine de lquation donne. Supposons, quune surface ferme et
lisse par morceau divise son domaine E n en domaine fini D et infini D . On considre
suivant D ou D et les condition aux limites les paramtres :

Problme intrieur aux limites. On cherche dans D une solution rgulire de

lquation (I.51), continue, diffrentiable jusqu' la frontire et satisfaisant la condition aux
limites :
(

u ( x)
u ( x)) x ( x),
v

2 2 0;

(I.52)

on suppose pour cela la condition F ( x) C ( D ) ; ( x) C ( ).

Problme extrieur aux limites. Dans le domaine D, il faut chercher la fonction
A
u( x) C 2 ( D ) C1 ( D ), rgulire linfini ( u ( x ) n 2 pour ( x : n 2,3),
x
satisfaisant lquation (I.51) et la condition aux limites (I.52).

La fonction u (x) doit tre rgulire linfini pour lunicit de la solution.

On considre suivant les valeurs des cfficients : , :
Les conditions aux limites de premire espce ( 0, 1) :

( x) x ( x)
et le problme correspondant sappelle le premier problme aux limites.
En particulier, lorsque lquation (I.51) est une quation de Poisson :
2u
u(x) 2 F(x).
k 1 x
k
n

le problme aux limites dans D (problme dans D ) avec les conditions aux limites de
premire espce sappelle le problme intrieur de Dirichlet et scrit : le problme D (le
problme exterieur de Dirichlet scrit: problme D )
Conditions aux limites de deuxime espce ( 1, 0) :

u ( x)
( x)
v x
et le problme correspondant sappelle le deuxime problme aux limites
Pour lquation de Poisson, on appelle problme aux limites dans D (problme D )
avec conditions aux limites de deuxime espce le problme intrieur (extrieur) de
Newmann et scrit : problme N (problme N ).
Conditions aux limites de troisime espce ( 0, 0) :
(

u ( x)
u ( x)) x ( x)
v

et le problme correspondant sappelle le troisime problme aux limites.

Notons, que pour rsoudre le problme de Newmann dans D , il faut une condition
complmentaire ( y )d 0.

I.8 Convenance du problme

On dit quun problme est pos correctement, si la solution : 1- existe, 2- est singulire,
3- est stable.
Exemple : Montrons, que le problme de Cauchy nest pas correct pour lquation
elliptique. Soit la fonction u ( x, y) dans un demi espce y 0 qui est une solution de
lquation de Laplace :

u ( x, y )

2u 2u

0
x 2 y 2

et satisfaisant les conditions de Cauchy

u ( x, y ) y 0 ( x);
La fonction ( x, y) u ( x, y)

u ( x, y )
y

( x) .
y 0

sin nxshny
sera ainsi une solution de lquation de
n2

Laplace,
donc :

xx yy (uxx sin nxshny) (u yy sin n shny) 0.

Par ailleurs, la fonction ( x, y) satisfait les conditions :
sin nxshny
) y 0 ( x);
n2
u ( x, y ) sin nxchny
sin nx
(

) y 0 ( x)
y
n
n

( x, y ) y 0 (u ( x, y )
( x, y )
y

y 0

Evaluons la diffrence des conditions initiales u( x, y), ( x, y) :

( x, y) u ( x, y) y 0 0,

( x, y) u ( x, y)

y
y

y 0

sin nx 1

n
n

Ainsi, pour n suffisamment grand, les conditions initiales pour rsoudre u ( x, y) , ( x, y) sont
proches et la diffrence des solutions sera:

( x, y) u ( x, y)

sin nxshny
n2

pour x 0, y 0 et par consquent :

shny e Ay e ny

n2
2n 2

pour

n.

Donc, la solution nest pas stable et le problme de Cauchy nest pas correct pour
lquation de Laplace.

I.9 Problmes non homognes

En physique mathmatique, on ramne les problmes complexes des problmes plus simples
pour pouvoir les rsoudre. Alors, il est utile dexaminer dans cette partie le principe de
superposition pour le cas de problmes mixtes de lquation hyperbolique.

Principe 1. Si uk (x, y) est la solution dun problme mixte :

Lu x u tt px u x x qx u Fk x, t ;
'

(I.43)

1u 1 u x 0,t 1u(0,t) 0 ; 2 u 2 u x (lt) 2 u(l,t) 0 ;

u(x,0) k (k); ut (x,0) k (x),

(I.44)
(I.45)

alors, la fonction :
m

u(x,t) Ak (x,t),

Ak const

k 1

(I.46)

sera la solution du problme :

m

Lu Ak Fk (x,t),

1u 0,

k 1

2u 0 ;
m

u(x,0) Ak k (x),

ut (x,0) Ak k (x).

k 1

k 1

Dmonstration.
En substituant la solution propose (I.46) dans (I.43)-(I.45) et en utilisant la linarit des
oprateurs diffrentiels L, 1 et 2 , on obtient :
m

Lu L Ak u k Ak L u k Ak Fk (x, t) ;
j u j Ak u k Ak j u k 0,
m

j 1, 2 ;

u (x,0) Ak u k (x,0) Ak k (x) ;

u (x,0) m
u(x,0) m
Ak k
Ak k (x),
t
t
1
1
Principe 2. Si u k ( x, t ) (k = 1, 2, 3, m)- est la solution de lquation homogne L u =
0, satisfaisant les conditions homognes aux limites 1u 0, 2 u 0 , la combinaison linaire
m

des

solutions

u ( x, t ) A k u k ( x, t ) sera
1

Lu 0, 1u 0, 2 u 0.

Dmonstration.

de

nouveau

la

solution

du

problme

En substituant la solution propose dans lquation, ainsi que dans les conditions aux
limites, on obtient :
m

k 1

k 1

Lu L k u k Ak L u k 0 ;
m

A u
k

k 1

Ak j u k 0,

j 1, 2.

k 1

Principe 3 (gnralis). Si u k ( x, t ), (k 1, 2, 3,...)est la solution de lquation L u =

satisfaisant les conditions aux limites 1 u 0, 2 u 0 , alors la srie
m

u(x,t) Ak u k (x,t) converge, admet une application des oprateurs L, 1 , 2 et elle sera de
1

nouveau la solution du problme Lu 0, 1u 0, 2 u 0.

Dmonstration.
En substituant la solution propose dans lquation, les conditions aux limites, en utilisant
les oprateurs L, 1 , 2 et la somme

, on obtient :
1

Lu L Ak u k Ak u k 0 ;

k 1

j u j Ak u k Ak j

u k 0,

j 1, 2.

Principe 4 (gnralis). Si la fonction u (x, t, ) dpendant du paramtre

1/ est la solution du problme Lu 0, 1u 0,
2/ lintgrale

K ()u( x, t, ) d

2u 0 ;

converge uniformment ;

3/ les oprateurs L, 1 2 peuvent tre mis sous intgrale ;

L K(()u(t, )d K(( u(x, t )d ;

j K u(x, t, )d K(( j u(x, t, )d,

alors, la fonction :
(x,t) K(()u(t, )d

sera encore une solution.

j 1, 2,

Dmonstration.
La convergence uniforme de lintgrale conduit la continuit de la fonction ( x, t ) et
la possibilit de mettre sous intgrale les oprateurs L, 1 , 2 permettent ( x, t ) de devenir
diffrentiable.
En appliquant les oprateurs L, 1 , 2 lgalit
( x, t ) K()u ( x, t , ) d .

et en tenant compte que u (x, t) est la solution de Lu 0, 1u 0, 2 u 0 , on arrive ainsi la

dmonstration demande.

I.10 Fonction delta de Dirac

Ltude des processus physiques de rpartition des masses, des forces et des sources
de chaleur, se ralise par la considration de leur caractre ponctuel.
Soit, par exemple, une masse unitaire distribue sur lintervalle [x0, x0 +
linaire ( x, x0 ) , alors sa rpartition unitaire est exprime par :
x0

(x, x

] de densit

)dx 1.

x0

La fonction ( x, x0 ) tant nulle lextrieur de lintervalle [x0, x0 +

(x, x0 ) 0

x [x 0 , x0 ];

(x, x

], alors :

)dx

x0

(x, x

)dx 1.

x0

Si la masse unitaire se concentre au point x0, alors le processus de transition limite est
obtenu pour la densit de la masse ponctuelle unitaire comme suit :
0 x x0
lim (x, x0 )
x x0 .
0
x x0

La masse reste inchange durant le processus de cette concentration :

lim (x, x0 )dx lim

0

x0

x0

(x, x0 )dx 1 x x0 dx.

Rappelons quune fonction u (x) sur lintervalle a, b est une limite faible des
fonctions intgrables u (x) pour 0 , si quelque soit f ( x) c a,b , il existe la limite :

lim u (x) f(x)dx u(x)f(x)dx .

0

(I.47)

Alors, si f(x) est continue sur laxe ox, donc pour ( x, x0 ) 0 ,on a :
x0

lim
0

x0

(x, x0 )f(x)dx lim f(x cp ) (x, x0 )dx f(x0 ).

x0

x0

Donc, la limite faible x x0 existe et peut scrire sous la forme :

x x f(x)dx

x0

x x f(x)dx f(x
0

0,

x0

Si lon pose : f(x) 1, alors,

x x0 dx

x0

x x dx 1
0

0.

x0

La fonction gnralise x x0 satisfaisant les conditions :

x x0,

1/

0
x x0

2/

f(x
x x f(x)dx 0 ,

x x0 ;

),

x0 a,b,

x0 a,b

quelque soit f (x) continue aux alentours du point x0, sappelle la fonction
fonction dimpulsion.

de Dirac ou la

La fonction ( x, x0 ), en respectant les conditions (1) et (2), peut tre tablie

comme :
1
, x x0 , x0 ,
x, x0

0, x x0 , x0

Srie de Fourier pour la fonction

. Dmontrons, que la fonction :

n (x, x0 )

2 n
k
k
sin x sin x0 .

l k 1
l
l

forme une fonction sur lintervalle [0, l].

Soit f (x)- une fonction arbitraire se dveloppant en srie de Fourier, alors :

2 l
k
k
lim n (x, x0 )f(x)dx lim f(x) sin xdx sin x0 .
n
l
l
n
k 1 l 0
0

Si f (x) est continue au point x0, alors :

l

lim n x, x0 f(x)dx f(x0 ),

Ainsi, on obtient pour la fonction de Dirac la reprsentation sous forme dune srie de
Fourier :

x x0

2
k
k
sin x sin x0 .

l k 1
l
l

Sur lintervalle [-l, l] , la srie de Fourier pour la fonction

x x0

sera :

1 1
k
k
k
k
cos
x cos
x0 sin x sin x0 .
2l l k 1
l
l
l
l

(I.48)