Début de la boite de navigation du chapitre
fin de la boite de navigation du chapitre
En raison de limitations techniques, la typographie souhaitable du titre, «
Dynamique des fluides compressibles : Vitesse du son Dynamique des fluides compressibles/Vitesse du son », n'a pu être restituée correctement ci-dessus.
Équations d'onde, vitesse du son
modifier
L'écoulement est supposé irrotationnel et parfait, et on traite les champs comme des perturbations de l'état stationnaire uniforme :
P
(
r
→
,
t
)
=
P
0
+
p
(
r
→
,
t
)
{\displaystyle P({\vec {r}},t)=P_{0}+p({\vec {r}},t)}
μ
(
r
→
,
t
)
=
μ
0
+
δ
μ
(
r
→
,
t
)
{\displaystyle \mu ({\vec {r}},t)=\mu _{0}+\delta \mu ({\vec {r}},t)}
p
(
r
→
,
t
)
{\displaystyle p({\vec {r}},t)}
,
δ
μ
(
r
→
,
t
)
{\displaystyle \delta \mu ({\vec {r}},t)}
et
v
→
(
r
→
,
t
)
{\displaystyle {\vec {v}}({\vec {r}},t)}
sont des infiniment petits du premier ordre.
Les équations décrivant un tel écoulement sont :
∂
v
→
∂
t
=
−
1
μ
0
∇
→
p
;
∂
p
∂
t
=
−
1
χ
T
d
i
v
v
→
{\displaystyle {\frac {\partial {\vec {v}}}{\partial t}}=-{\frac {1}{\mu _{0}}}{\vec {\nabla }}p~~~~~~;~~~~~~{\frac {\partial p}{\partial t}}=-{\frac {1}{\chi _{T}}}\mathrm {div} {\vec {v}}}
(où
χ
T
{\displaystyle \chi _{T}~}
est le coefficient de compressibilité isotherme.)
D'où l’on tire les équations d'onde :
∂
2
p
∂
t
2
−
1
μ
0
χ
T
Δ
p
=
0
;
∂
2
v
→
∂
t
2
−
1
μ
0
χ
T
Δ
v
→
=
0
{\displaystyle {\frac {\partial ^{2}p}{\partial t^{2}}}-{\frac {1}{\mu _{0}\chi _{T}}}\Delta p=0~~~~~~;~~~~~~{\frac {\partial ^{2}{\vec {v}}}{\partial t^{2}}}-{\frac {1}{\mu _{0}\chi _{T}}}\Delta {\vec {v}}=0}
(en utilisant
r
o
t
→
v
→
=
0
→
{\displaystyle {\overrightarrow {\rm {rot}}}~{\vec {v}}={\vec {0}}}
).
On en déduit la vitesse du son :
c
=
1
μ
0
χ
T
{\displaystyle c={\frac {1}{\sqrt {\mu _{0}\chi _{T}}}}}
Pour la première équation, on part de l'équation d'Euler (on néglige la viscosité et la gravité) :
∂
v
→
∂
t
+
(
v
→
⋅
∇
→
)
v
→
=
−
1
μ
∇
→
p
{\displaystyle {\frac {\partial {\vec {v}}}{\partial t}}+({\vec {v}}\cdot {\vec {\nabla }}){\vec {v}}=-{\frac {1}{\mu }}{\vec {\nabla }}p}
Au premier ordre,
(
v
→
⋅
∇
→
)
v
→
{\displaystyle ({\vec {v}}\cdot {\vec {\nabla }}){\vec {v}}}
est négligé devant
∂
v
→
∂
t
{\displaystyle {\frac {\partial {\vec {v}}}{\partial t}}}
, et on a
1
μ
∇
→
p
=
1
μ
0
(
1
−
δ
μ
μ
0
+
.
.
.
)
∇
→
p
≈
1
μ
0
∇
→
p
{\displaystyle {\frac {1}{\mu }}{\vec {\nabla }}p={\frac {1}{\mu _{0}}}(1-{\frac {\delta \mu }{\mu _{0}}}+...){\vec {\nabla }}p\approx {\frac {1}{\mu _{0}}}{\vec {\nabla }}p}
. On en déduit l'équation utilisée, dite équation d'Euler linéarisée.
L'autre relation est issue de la définition du coefficient de compressibilité isotherme :
χ
T
=
−
1
V
(
∂
V
∂
P
)
T
=
1
μ
(
∂
μ
∂
P
)
T
≈
1
μ
δ
μ
p
⇒
δ
μ
≈
μ
χ
T
p
≈
μ
0
χ
T
p
{\displaystyle \chi _{T}=-{\frac {1}{V}}\left({\frac {\partial V}{\partial P}}\right)_{T}={\frac {1}{\mu }}\left({\frac {\partial \mu }{\partial P}}\right)_{T}\approx {\frac {1}{\mu }}{\frac {\delta \mu }{p}}~~\Rightarrow ~~\delta \mu \approx \mu \chi _{T}~p\approx \mu _{0}\chi _{T}~p}
On a ensuite, avec l'équation de conservation de la masse :
d
i
v
v
→
=
−
1
μ
∂
μ
∂
t
≈
−
1
μ
0
∂
(
δ
μ
)
∂
t
≈
−
χ
T
∂
p
∂
t
{\displaystyle \mathrm {div} {\vec {v}}=-{\frac {1}{\mu }}{\frac {\partial \mu }{\partial t}}\approx -{\frac {1}{\mu _{0}}}{\frac {\partial (\delta \mu )}{\partial t}}\approx -\chi _{T}{\frac {\partial p}{\partial t}}}