Cours de physique IV (PHYS-206) – Prof.

Paolo Ricci – SPC

21 mars 2025

Corrigé 5 : Ondes stationnaires

1 Ondes stationnaires dans une colonne d’eau
Un long cylindre vertical de rayon R et de longueur L, ouvert à son extrémité supérieure, est rempli
d’une colonne d’eau de hauteur h. La hauteur de la colonne d’air est notée H = L − h. Une pompe
de débit Q permet d’ajuster h. Un diapason, dont la fréquence propre est ν, est placé au sommet du
cylindre. Dans un premier temps, on suppose que l’eau est un milieu opaque (qui ne permet pas la
transmission d’ondes sonores).
(a) Déterminer les hauteurs de la colonne d’air Hn pour lesquelles on observe une résonance associée
au ne mode propre de l’onde sonore du diapason. On note uair la vitesse de propagation des
ondes sonores dans l’air.
On commence par s’intéresser aux modes propres dans la colonne d’air pour une hauteur H
fixée. On exprime la fonction d’onde stationnaire à partir de l’indication,

ξ(x, t) = ξ0 sin(kx + α) cos(ωt + φt ), (1)

où φt correspond à un déphasage temporel. On choisit l’origine x = 0 à l’interface entre

l’eau et l’air, et x = H > 0 à l’extrémité supérieure, ouverte, du cylindre. Les conditions
aux limites à imposer sont donc
— un nœud en x = 0 (l’eau étant considérée comme un milieu opaque), ce qui implique
que α = 0, et
— un ventre en x = H, puisque l’extrémité supérieure du cylindre est ouverte.
La dérivée spatiale de la fonction d’onde est nulle au ventre, ce qui donne
∂ξ
(H, t) = ξ0 k cos(kH) cos(ωt) = 0 =⇒ cos(kH) = 0. (2)
∂x
On obtient ainsi les nombres d’onde kn et les fréquences propres νn des modes :

π kn u (2n + 1)
kn H = (2n + 1) , νn = = uair , (3)
2 2π 4H
où n = 0, 1, 2, ... et uair est la vitesse de propagation.
On remarque cependant que, selon les instructions de l’exercice, la fréquence ν du diapason
est fixée, alors que la hauteur de la colonne d’air H est variable. On transforme donc la
condition de résonance (3) afin de fixer ν pour obtenir un ensemble de hauteurs Hn pour
lesquelles une onde stationnaire existe,

(2n + 1)
Hn = uair . (4)
4ν

(b) À l’aide de la pompe, on fait monter le niveau d’eau. Déterminer le temps ∆t séparant deux
instants où le cylindre entre en résonance avec le diapason.
On a une résonance lorsque la fréquence du diapason est égale à une des fréquences propres.
La distance entre deux hauteurs de liquide pour deux résonances consécutives est donc :
uair
∆H = Hn+1 − Hn = . (5)
2ν

1
 La variation de hauteur ∆H dans un temps ∆t est : ∆H = (Q/πR2 )∆t. Finalement :

πR2 uair
∆t = = 237 s. (6)
2Qν

Dans un second temps, on ne considère plus l’eau comme étant un milieu opaque. On considère
d’abord une onde progressive incidente sinusoïdale d’amplitude ξi0 et de fréquence ω. À l’interface
entre l’eau et l’air, l’onde incidente se décompose en une onde réfléchie, se propageant dans l’air, et
une onde transmise dans l’eau.
(c) En posant les conditions de continuité de l’amplitude du déplacement et de la pression à
l’interface, dériver les amplitudes ξr0 et ξt0 des ondes transmises et réfléchies en fonction de
l’amplitude incidente ξi0 . On note κair et κeau les coefficients de compressibilités des milieux
respectifs.
À l’interface entre les deux milieux, l’onde transmise ne change pas de fréquence ω, mais
uniquement de longueur d’onde, caractérisée par le nombre d’onde k. Il en découle que

ω = kair uair = keau ueau . (7)

On pose donc la forme générale des ondes progressives pour l’onde incidente, réfléchie et
transmise :

Onde incidente : ψi (x, t) = ξi0 sin(kair x − ωt) (8)

Onde réflechie : ψr (x, t) = ξr0 sin(−kair x − ωt) (9)
Onde transmise : ψt (x, t) = ξt0 sin(keau x − ωt) (10)

Sans perte de généralité, on pose x = 0 la position de l’interface. Il est nécessaire d’imposer

des conditions de continuité à l’interface. La première condition est celle de continuité des
amplitudes, imposant que l’amplitude des ondes résultantes dans l’air ψ1 = ψi + ψr et dans
l’eau ψ2 = ψt soient égales à l’interface :

ψ1 (0, t) = ξi0 sin(−ωt) + ξr0 sin(−ωt) = ξt0 sin(−ωt) = ψ2 (0, t) =⇒ ξi0 + ξr0 = ξt0 . (11)

La deuxième condition est la condition de continuité sur la pression. La pression est donnée
par l’expression suivante (vue en cours) :
∂ψ
p(x, t) = p0 − κ (x, t), (12)
∂x
où p0 est la pression d’équilibre qui est uniforme dans l’eau et l’air. Ceci donne donc la
condition suivante :
∂ψ1 ∂ψ2
κair = κeau . (13)
∂x ∂x
Explicitement, ceci donne une deuxième condition sur les amplitudes ξi0 , ξr0 et ξt0 , valable
en tout temps t :
κair kair (ξi0 − ξr0 ) = κeau keau ξt0 . (14)
En résolvant les équations (11) et (14) pour ξr0 et ξt0 , on obtient les expressions suivantes :
κeau keau − κair kair 0
ξr0 = ξ = Rξi0 , (15)
κeau keau + κair kair i
2κair kair
ξt0 = ξ 0 = T ξi0 , (16)
κeau keau + κair kair i
où R et T sont appelés « coefficients de Fresnel ».

2
Finalement, on s’intéresse à la possibilité d’observer des ondes stationnaires qui se développent dans
l’eau et l’air.
(d) Déterminer les conditions nécessaires pour qu’une onde stationnaire soit présente dans les deux
milieux en fonction des nombres d’onde kair et keau , des coefficients de compressibilité κair et
κeau , de H ainsi que L.
L’onde transmise à l’interface possède une longueur d’onde différente, mais une fréquence
est égale à l’onde incidente. Par conséquent, par la relation d’onde dans l’air et l’eau, on
trouve que les nombres d’onde kair et keau doivent satisfaire la relation donnée en (7). On
écrit la solution d’onde
dans l’air : ψ1 (x, t) = A sin(kair x + α1 ) cos(ωt + φt ), (17)
et dans l’eau : ψ2 (x, t) = B sin(keau x + α2 ) cos(ωt + φt ). (18)
Pour simplifier les développements, on choisit de poser x = 0 le fond du cylindre, et donc
x = L l’extrémité ouverte du cylindre ainsi que x = L − H la position de l’interface. Les
conditions aux limites sont :
— un nœud en x = 0 (fond du cylindre), ce qui donne : α2 = mπ, m ∈ N,
— un ventre en x = L, car l’extrémité supérieure du cylindre est ouverte, ce qui donne
α1 + kair L = π/2 + nπ, n ∈ N.
Selon l’indication, on se concentre sur les modes fondamentaux par simplicité, donc n =
m = 0 pour obtenir
dans l’air : ψ1 (x, t) = A cos(kair (L − x)) cos(ωt + φt ), (19)
dans l’eau ψ2 (x, t) = B sin(keau x) cos(ωt + φt ). (20)

À l’interface, on impose la continuité des solutions d’onde. En effet, ψ représente le dépla-

cement infinitésimal du milieu. Puisque le déplacement de l’interface doit être le même
dans les deux milieux, cela implique donc que ψ1 (H, t) = ψ2 (H, t). Ceci donne l’équation
suivante :
A cos(kair (L − H)) = B sin(keau H). (21)
Comme au point (c), la pression doit également être continue à l’interface. On applique de
la même manière qu’en (14)
κair kair A sin(kair (L − H)) = κeau keau B cos (keau H). (22)
Ces deux équations (21) et (22) peuvent s’écrire sous forme matricielle :
! ! !
cos(kair (L − H)) − sin(keau H) A 0
= (23)
κair kair A sin(kair (L − H)) −κeau keau cos (keau H) B 0

Pour que (A, B) donne une solution non triviale, il faut imposer que le déterminant de la
matrice soit nul,
κair kair sin(kair (L − H)) sin(keau H) − κeau keau cos(kair (L − H)) cos (keau H) = 0. (24)
Pour montrer que cette équation est cohérente, il est intéressant de considérer le cas où les
deux milieux sont les mêmes, c.-à.-d. k = kair = keau ainsi que κ = κair = κeau . L’équation
donne
κk sin(k(L − H)) sin(kH) − κk cos(k(L − H)) cos (kH) = 0 (25)
cos(kL − 2kH) − cos(kL) − cos(kL) − cos(kL − 2kH) = 0 (26)
cos(kL) = 0, (27)
ce qui est la condition pour une onde stationnaire dans une cavité semi-ouverte.

3
2 Timbre d’un instrument à cordes
Dans cet exercice, on propose d’étudier les vibrations d’une corde de longueur L dans les cas où elle
est initialement pincée ou frappée. On considère la corde sujette à une tension T et de masse linéique
µ. La corde est fixe à chacune de ses extrémités.
Pour commencer, on considère une corde pincée en son milieu, c’est-à-dire que la corde est initialement
de forme triangulaire, avec le sommet situé à égale distance des extrémités fixes. Le sommet est
déplacé de A par rapport à la position au repos de la corde.
(a) Déterminer la forme générale de la solution de l’équation d’onde. Identifier les modes propres
du système.
On remarque que la condition initiale exhibe un profil qui ne correspond pas à la forme
sinusoïdale habituelle de la solution à l’équation d’onde. En revanche, ce profil peut être
décomposé en une superposition de modes propres, c’est-à-dire l’ensemble de fonctions
solutions de l’équation d’onde satisfaisant aux conditions aux bords. Afin de déterminer
l’évolution de la fonction d’onde totale ψ(x, t), on procède par l’approche suivante,
1. On cherche une solution générale à l’équation d’onde par séparation des variables afin
d’identifier les modes propres,
2. On projette la condition initiale sur la base fonctionnelle formée par les modes propres
trouvés au point précédent,
3. Les coefficients obtenus par projection de la condition initiale sont insérés dans
l’expression de la solution totale en tenant compte de la dépendance temporelle de
chaque mode propre déterminé à l’étape 1.
On résout l’équation d’onde unidimensionnelle par séparation de variables. On utilise
l’indication en écrivant ψn (x, t) = Xn (x)Tn (t), qu’on introduit dans l’équation d’onde afin
d’obtenir des équations harmoniques, associées au ne mode propre, séparées pour Xn (x) et
Tn (t),
Xn00 (x) + kn2 Xn (x) = 0, (28)
Tn00 (t) + u2 kn2 Tn (t) = 0. (29)
En résolvant les équations (28) et (29), on obtient les solutions de forme générale suivantes :
Xn (x) = A(n) (n)
x cos(kn x) + Bx sin(kn x), (30)
(31)
(n) (n)
Tn (t) = At cos(kn ut) + Bt sin(kn ut).
Les conditions aux bords fixes s’expriment
ψn (0, t) = ψn (L, t) = 0. (32)
En imposant (32), on obtient des conditions sur les paramètres Ax , et kn :
(n) (n)
Bx
Xn (0) = 0 =⇒ A(n)
x = 0, (33)
Xn (L) = 0 =⇒ Bx(n) sin(kn L) = 0. (34)
Cette dernière égalité impose que sin(kn L) = 0. Par conséquent, on obtient que kn = nπ/L.
Cette condition sur le nombre d’onde permet d’identifier les modes propres du système,
(35)
 
ψn (x, t) = Xn (x)Tn (t) = sin(kn x) bn cos(kn ut) + cn sin(kn ut) ,

où on a rassemblé les constantes bn = Bx At et cn = Bx Bt . Une solution générale,

(n) (n) (n) (n)

satisfaisant les conditions aux bords, s’écrit donc

∞
(36)
X
ψ(x, t) = ψn (x, t).
n=1
Ceci conclut la première partie du raisonnement (étape 1. de la liste en page 4).

4
(b) Imposer que la solution de l’équation d’onde trouvée au point (a) satisfasse les conditions
initiales, en déduire que
∞
8A(−1)n (2n + 1)πx
X  
ψ(x, t = 0) = 2 2
sin .
n=0
(2n + 1) π L

On écrit formellement les conditions initiales et les conditions aux bords de la corde pincée.
Aux conditions initiales, la corde est pincée en son centre et immobile. Ceci ce traduit par
une fonction d’onde « chapeau » dont la dérivée temporelle est nulle :
2Ax

si 0 ≤ x ≤ L/2,


∂ψ
et L (37)

(x, 0) = 0, ∀x ∈ [0, L], ψ(x, 0) =
∂t 2A
(L − x) sinon.



L
On cherche dès à présent à déterminer les coefficients bn et cn de telle sorte que les conditions
initiales (37) soient satisfaites, notamment,
∞
(38)
X
ψ(x, 0) = bn sin(kn x),
n=1
∞
∂ψ
(39)
X
(x, 0) = cn kn u sin(kn x) = 0.
∂t n=1

L’équation (39) doit être satisfaite pour n’importe quel choix de kn = nπ/L et x ∈ [0, L], ce
qui implique en particulier que les coefficients cn s’annulent, cn = 0. Il reste donc à projeter
la condition chapeau initiale (37) sur une base de fonctions sinus (38). On remarque que
cette base représente un sous-ensemble de celle utilisée pour une série de Fourier, c’est-à-dire
que l’expression (38) est une série de Fourier dont les coefficients associés aux termes cosinus
s’annulent (voir note en fin de question). Similairement à une série de Fourier, on note que
les fonctions de base sont orthogonales,
2
Z L 
πnx πmx
  
sin sin dx = δnm . (40)
L 0 L L
Cette propriété peut être exploitée afin de déterminer les coefficients bn : par multiplication
de l’éq. (38) par (2/L) sin(πnx/L) suivi d’une intégration entre x = 0 et L, on obtient :
2
Z L 
πnx

bn = ψ(x, 0) sin dx. (41)
L 0 L
On substitue l’expression de la fonction chapeau pour poursuivre l’intégration,
2
Z L/2 2A nπx

2
 Z L 2A nπx
 
bn = x sin dx + (L − x) sin dx (42)
L 0 L L L L/2 L L
4A L/2
Z
nπx
 L/2 
nπ
Z   
= 2 x sin dx + y sin (L − y) dy (avec y = L − x) (43)
L 0 L 0 L
Z L/2
4A nπx
 
= 2 1 − (−1) n

x sin dx (44)
L 0 L
Z 1/2
= 4A 1 − (−1)n (45)


z sin(nπz) dz (avec z = x/L).
0

Entre l’équation (43) et (44), on a substitué

nπy nπy nπy
     
sin nπ − n
= (−1) sin − = −(−1) sin n
. (46)
L L L

5
 On remarque que les coefficients bn tels que n soit pair s’annulent, b2m = 0. Ce résultat
est effectivement attendu en remarquant que la condition initiale est une fonction paire
sur l’intervalle d’intégration (autour de x = L/2), alors que les fonctions de bases pour
n = 2m sont impaires. On évalue l’intégrale restante pour le cas où n = 2m + 1 est impair
en procédant par parties,
Z 1/2 z z=1/2 1/2 1 Z
z sin(nπz) dz = − cos(nπz) + cos(nπz) dz = · · ·
0 nπ z=0 0 nπ
1 z=1/2 sin(mπ + π/2) (−1)m
··· = 0 + sin(nπz) = = . (47)
(nπ)2 z=0 ((2m + 1)π)2 (2m + 1)2 π 2

On a donc, pour le cas impair, que le coefficient b2m+1 est donné par

8A(−1)m
b2m+1 = . (48)
(2m + 1)2 π 2

Ce dernier résultat conclut la deuxième partie du raisonnement proposé en page 4 (étape 2).
Note sur les séries de Fourier. On remarque que si le coefficient Ax n’était pas contraint
à être nul, éq. (33), l’équation (38) aurait pris la forme
∞
a0 X nπx nπx
    
ψ(x, 0) = f(x) = + an cos + bn sin , (49)
2 n=1
L L

qui est la forme générale de la série de Fourier d’une fonction f(x) sur un intervalle
périodique x ∈ [0, L] Le facteur 1/2 est introduit afin que la définition des coefficients an
donnés ci-dessous soit valide pour tout n. Les coefficients satisfont

2
Z L 
nπx

2
Z L 
nπx

an = f(x) cos dx, et bn = f(x) sin dx. (50)
L 0 L L 0 L

La base de fonctions trigonométriques est orthogonale, c.-à-d. que les relations (40), et (51)
sont satisfaites :
2
Z L 
mπx

nπx
 
2
Z L 
mπx nπx
  
cos cos dx = δmn , cos sin dx = 0. (51)
L 0 L L L 0 L L

(c) Dériver l’évolution temporelle du déplacement de la corde. Quelle est l’intensité de l’onde sur la
corde en fonction du temps ?
On substitue les coefficients bn et cn = 0 dérivés ci-dessus dans la solution générale de
l’équation d’onde (36) pour finalement obtenir
∞
8A(−1)m (2m + 1)πx (2m + 1)πut
   
(52)
X
ψ(x, t) = 2 2
sin cos
m=0
(2m + 1) π L L
∞
(2m + 1)πx
 
(53)
X
= βm (t) sin .
m=0
L

On calcule l’intensité de l’onde sur l’intervalle [0, L], définie comme :

1
Z L
I(t) = |ψ(x, t)|2 dx. (54)
L 0

En utilisant l’orthogonalité des modes sinusoïdaux, comme donné en (40), il est possible de

6
 montrer que, par substitution de (53),

1 L
Z
I(t) = |ψ(x, t)|2 dx (55)
L 0
∞ X ∞ Z L
1 X (2m + 1)πx (2n + 1)πx
   
= βm (t)βn (t) sin sin dx (56)
L m=0 n=0 0 L L
∞ ∞
1 X
(57)
X
2
= βm (t) = Im (t).
2 m=0 m=0

L’intensité de chaque harmonique Im (t) est donnée par :

2
1 8A (2m + 1)πut
 
Im (t) = 2 2
cos . (58)
2 (2m + 1) π L

(d) Refaire les points (a)–(c) en considérant une corde pincée au 1/3 de sa longueur.
Le raisonnement reste identique à celui du point (c). Il faut cependant adapter l’inté-
grale (42)–(45) associée au calcul des coefficients bn aux nouvelles conditions initiales
décrites par l’énoncé :
3A L

x si 0 ≤ x ≤ ,


L 3 (59)

ψ(x, 0) =
3A
(L − x) sinon,



2L
∂ψ
et (x, 0) = 0, ∀x ∈ [0, L]. (60)
∂t
On note que la condition (60) implique à nouveau cn = 0, cf. éq. (39). Pour bn ,

2 L nπx
Z  
bn = ψ(x, 0) sin dx (61)
L 0 L
Z L/3 Z L
2 3A nπx 3A nπx
    
= x sin dx + (L − x) sin dx (62)
L 0 L L L/3 2L L
1/3 (−1)n 2/3
Z Z 
= 6A z sin(nπz) dz − y sin(nπy) dy avec y = 1 − z; z = x/L (63)
0 2 0
z=1/3 y=2/3
z cos(nπz) sin(nπz) (−1)n y cos(nπy) sin(nπy)
 
= 6A − + − − + (64)
nπ n2 π 2 z=0 2 nπ n2 π 2 y=0
sin(nπ/3) − (nπ/3) cos(nπ/3) sin(2πn/3) − (2πn/3) cos(2πn/3)
 
= 6A 2 2
− (−1)n . (65)
n π 2n2 π 2

Entre (62) et (63), on a utilisé que

sin(nπ(1 − y)) = sin(nπ − nπy) = −(−1)n sin(nπy). (66)

L’amplitude sera donc donnée par

∞
nπx nπut
   
(67)
X
ψ(x, t) = bn sin cos .
n=1
L L

On remarque des différences majeures par rapport au résultat trouvé en (c). En effet, lorsque
la corde était initialement pincée au milieu, par symétrie, les coefficients b2m étaient nuls

7
 par parité : la fonction sin(2mπx/L) était impaire par rapport au point x = L/2 alors que
ψ(x, 0) était paire. Or, dans le cas de la corde pincée au 1/3 de la longueur, ces coefficients
sont non nuls et des ondes stationnaires antisymétriques contribuent donc aux vibrations
de la corde.
On aimerait également étudier le cas d’une corde frappée plutôt que pincée, comme pour un piano
par exemple. Initialement, la corde est supposée être à sa position d’équilibre, mais sur un intervalle
de longueur a centré en L/2, une vitesse initiale v0 est donnée aux éléments de la corde. Cette
condition se traduit en

si x ∈ [(L − a)/2, (L + a)/2],

(
∂ψ v0
ψ(x, 0) = 0 et =
∂t t=0 0 sinon.

(e) Dériver l’évolution temporelle de la corde frappée.

Similairement aux points précédents, on détermine les coefficients de Fourier bn et cn . On
considère pose à nouveau les équations (38) et (39),
∞
(68)
X
ψ(x, 0) = bn sin(kn x) = 0,
n=1
∞
∂ψ
(69)
X
(x, 0) = cn kn u sin(kn x).
∂t n=1

Comme ψ(x, 0) = 0, tous les coefficients bn sont cette fois nuls. Il suffit donc de déterminer
les coefficients cn en projetant la condition initiale (∂ψ/∂t)(x, 0) sur la base de fonctions
sinus. Une approche similaire à celle aboutissant à l’équation (41) donne donc

2 1 L ∂ψZ
nπx
 
cn = (x, 0) sin dx (70)
L ukn 0 ∂t L
2v0 (L+a)/2 nπx
Z  
= sin dx (71)
nπu (L−a)/2 L
x=(L+a)/2
2v0 L nπx
 
= − cos (72)
nπu nπ L x=(L−a)/2
2v0 L nπ nπa nπ nπa
    
=− 2 2
cos + − cos − (73)
n π u 2 2L 2 2L
4v0 L nπ nπa
   
= 2 2 sin sin . (74)
n π u 2 2L

On a donc à nouveau une annulation des coefficients pairs c2m = 0 en raison de la parité de
la fonction projetée sur la base sinus, alors que

4v0 L (2m + 1)πa

 
c2m+1 = (−1)m 2 2
sin . (75)
(2m + 1) π u 2L

Par conséquent, l’évolution temporelle de la corde frappée est donnée par :

∞
4v0 L (2m + 1)πa (2m + 1)πx (2m + 1)πut
X      
m
ψ(x, t) = (−1) 2π2u
sin sin sin .
m=0
(2m + 1) 2L L L
(76)

8
3 Corde dans un milieu visqueux
On considère une corde de longueur L et de masse linéique µ dans un milieu visqueux. Les bords de
la corde sont considérés comme étant fixes. Chaque élément infinitésimal de corde dx est soumis à
une force de frottement infinitésimale dF = −λv dx, où λ > 0 est le coefficient de frottement par
unité de longueur. La corde est soumise à une tension T. La vitesse de propagation de la perturbation
est notée u.
(a) Montrer que l’équation d’onde, en tenant compte des frottements, s’écrit :
∂ 2 ψ λ ∂ψ 2
2∂ ψ
+ = u .
∂t2 µ ∂t ∂x2

L’équation d’onde unidimensionnelle sans frottements est donnée par :

∂2ψ 2
2∂ ψ
= u .
∂t2 ∂x2
Cette équation peut s’obtenir à l’aide d’une approche Newtonienne, où l’interaction d’un
point avec son voisinage est due à la tension, à laquelle on ajoute la contribution de la force
de frottement. Puisque F = −λv, la force de frottement s’écrit donc F = −λ(∂ψ/∂t). Par
conséquent, la 2e loi de Newton pour la corde est donnée par :

∂2ψ X ∂ψ
µ dx 2
= Fi = −λ dx + Ty (x + dx) + Ty (x). (77)
∂t i
∂t

Comme vu en cours, la tension T est de norme constante en tout point de la corde et la

composante verticale s’écrit Ty (x) = T sin(θ(x)) ≈ T θ(x) = T(∂ψ/∂x). En divisant par µ dx,
on obtient le résultat souhaité, sachant que u2 = T /µ pour la corde :

∂ 2 ψ λ ∂ψ u2 ∂ψ ∂ψ 2
2∂ ψ
 
+ = (x + dx) − (x) ≈ u . (78)
∂t2 µ ∂t dx ∂x ∂x ∂x2

(b) En utilisant la méthode de séparation de variables, dériver les modes propres d’une corde dans
un milieu visqueux. Discuter de l’évolution temporelle des différents régimes observés pour ces
modes propres.
On utilise la méthode de séparation de variables en posant ψ(x, t) = X(x)T(t). En injectant
cette expression dans l’équation d’onde amortie et en choisissant k 2 comme solution de
l’équation, on obtient :
1 λ 0 X 00 (x)
 
T 00
(t) + T (t) = = −k 2 = const. (79)
u2 T(t) µ X(x)
Ceci donne deux équations différentielles linéaires pour X(x) et T(t) :
X 00 (x) + k 2 X(x) = 0, (80)
λ
T 00 (t) + T 0 (t) + k 2 u2 T(t) = 0. (81)
µ
L’équation (80) correspond à l’équation harmonique dont la forme générale est donnée par :
X(x) = A sin(kx) + B cos(kx). (82)
Les conditions aux bords fixes imposent que :
X(0) = 0 =⇒ B = 0, (83)
X(L) = 0 =⇒ A sin(kL) = 0. (84)

9
 L’équation (83) impose donc que kL = nπ, ce qui implique que la partie spatiale de ψ est
donnée par : X(x) = A sin(kn x) avec kn = nπ/L.
Pour la partie temporelle, on observe que l’équation (81) admet les solutions de la forme
T(t) = Aeat . En appliquant cet Ansatz, on obtient que a obéit à la relation suivante :

λ
a2 + a + u2 k 2 = 0, (85)
µ
et donc s
λ λ2 √
a=− ± 2
− k 2 u2 = −Λ ± ∆. (86)
2µ 4µ
La forme générale pour T(t) est donc donnée par :
√ √
T(t) = e−Λt (Ae ∆t + Be− ∆t ) (87)
√ √ √
On peut distinguer trois cas particuliers : ∆ ∈ R, ∆ ∈ iR et ∆ = 0. La notation iR
dénote l’ensemble des nombres purement
√ imaginaires,
√ c.-à-d. les nombres complexes dont
la partie réelle est nulle. Les cas ∆ ∈ R et ∆ ∈ iR correspondent respectivement aux
régimes sur-critiques√et sous-critiques comme vu pour l’équation de l’oscillateur harmonique
amorti. Pour le cas ∆ = 0, il faut procéder plus méthodiquement. En effet, a1 = a2 = −Λ,
ce qui implique que les deux solutions trouvées précédemment sont égales. Pour trouver
la seconde solution, on utilise la méthode de la variation de la constante. En posant
T(t) = A(t)e−Λt , on peut montrer que√le terme A(t) doit satisfaire A00 (t) = 0, ce qui
implique que A(t) = A + Bt. Le cas ∆ = 0 est appelé régime critique et la partie
temporelle est donnée par T(t) = e−Λt (A + Bt). Finalement, les modes propres peuvent
s’écrire de la manière suivante :
√ √
∆n ∈ R (sur-critique) : ψn (x, t) = e−Λt sin(kn x)(Ae + Be− (88)
p
∆n t ∆n t
),
∆n = 0 (critique) : −Λt
(89)
p
ψn (x, t) = e sin(kn x)(A + Bt),
√ √
∆n ∈ iR (sous-critique) : ψn (x, t) = e−Λt sin(kn x)(Aei| ∆n |t
+ Be−i| ∆n |t
(90)
p
),

où ∆n = Λ2 − kn2 u2 .
(c) Obtenir la relation de dispersion ω = ω(k).
On rappelle l’expression de la transformée de Fourier d’une fonction ψ(x, t)
Z Z
F(ψ(x, t)) = ψ̃(k, ω) = dx dt ψ(x, t)e−i(kx−ωt) . (91)

Les propriétés de la transformée de Fourier implique que

∂ψ ∂2ψ ∂2ψ
     
F = −iωF(ψ), F = −ω 2 F(ψ) et F = −k 2 F(ψ). (92)
∂t ∂t2 ∂x2
Ces propriétés peuvent être trouvées en utilisant l’intégration par parties et en considérant
qu’aux limites t → +∞ et x → ±∞, la fonction d’onde s’annule ψ(x, t) → 0. En considérant
la transformée de Fourier de l’équation d’onde amortie, on obtient :
∂ 2 ψ λ ∂ψ ∂2ψ λ
   
F 2
+ − u2 2 = 0 =⇒ −ω 2 − i ω + u2 k 2 ψ̃(k, ω) = 0. (93)
∂t µ ∂t ∂x µ
Puisque ψ̃(k, ω) est non nul pour une solution non triviale, on obtient donc que la relation
de dispersion peut être écrite comme suit :
λ
u2 k 2 = ω 2 + i ω. (94)
µ

10
 On remarque que cette relation de dispersion est très similaire à la relation de dispersion
pour un milieu sans dissipation ω 2 = u2 k 2 . Le terme complexe iλω/µ implique donc une
décroissance exponentielle dans le temps.
(d) En considérant ω ∈ C et k ∈ R, que peut-on dire sur l’évolution temporelle d’un paquet d’onde ?
En résolvant (94) pour la variable ω, on obtient que :
s
λ λ2
ω = −i ± u2 k 2 − . (95)
2µ 4µ2

Cette expression est analogue à (86). Cela veut donc dire que pour certaines valeurs du
nombre d’onde k, ω peut être purement imaginaire, c.-à-d. ω ∈ iR, ce qui correspondrait à
un amortissement sur-critique de ces solutions d’onde. Ce résultat coïncide donc avec les
modes de Fourier trouvés en (88).

11