物理学上の現象
音響流は、高振幅の 音響 振動の吸収によって駆動される流体内の定常流です。この現象は、音響放射体の近く、または クント管内 の定在波で観察できます。音響流は 、1884年に レイリー卿 によって初めて説明されました。 [1]
これは、流れによる音の発生の反対ですが、あまり知られていません。
音が伝播媒体で吸収される状況は 2 つあります。
バルク流(「エッカート流」)中の伝播中。 [2] 減衰係数は であり 、 ストークスの法則(音の減衰) に従います。この効果は周波数が高いほど強くなり、水中(1 MHz で約 100 m)よりも空気中(1 MHz で約 10 cm の 特性距離で減衰が発生)の方がはるかに大きくなります。空気中では、 石英風 として知られています 。
α
=
2
η
ω
2
/
(
3
ρ
c
3
)
{\displaystyle \alpha =2\eta \omega ^{2}/(3\rho c^{3})}
α
−
1
{\displaystyle \alpha^{-1}}
α
−
1
{\displaystyle \alpha^{-1}}
境界付近では、音響流が発生します (「レイリー流」)。音が境界に到達したとき、または境界が静止した媒体内で振動しているときのいずれかです。 [3]壁が平行に振動すると、 ストークス振動境界層 内で減衰した振幅のせん断波が発生します。この効果は、1 MHz で空気と水の両方で数マイクロメートルのオーダーの減衰長に局在します 。音波とマイクロバブル、弾性ポリマー、 [4] 、さらには生物細胞 [5] との相互作用によって生成されるストリーミングフローは、境界駆動型音響流の例です。
δ
=
[
η
/
(
ρ
ω
)
]
1
/
2
{\displaystyle \delta =[\eta /(\rho \omega )]^{1/2}}
レイリーストリーミング
の 速度場に対応する平面定在音波を考えます 。問題の特性(横方向)次元を とします 。ここで説明した流れ場は非粘性流れに対応します。しかし、固体壁の近くでは粘性効果が重要になります。その場合、厚さ または侵入深さ の境界層が存在します 。レイリー流は、近似 で最もよく視覚化されます。 の場合と同様に 、速度成分は よりもはるかに小さくなります。さらに、境界層内の特性時間スケールは 、音響時間スケール と比較して 非常に大きくなります( が小さいため) 。これらの観察結果は、境界層内の流れは非圧縮性と見なせることを示唆しています。
あなた
(
x
、
t
)
=
ヴ
0
コス
け
x
コス
ω
t
=
ε
コス
け
x
ℜ
(
e
−
私
ω
t
)
{\displaystyle U(x,t)=v_{0}\cos kx\cos \omega t=\varepsilon \cos kx\Re (e^{-i\omega t})}
け
=
2
π
/
λ
=
ω
/
c
{\displaystyle k=2\pi /\lambda =\omega /c}
l
{\displaystyle l}
δ
=
(
2
ν
/
ω
)
1
/
2
{\displaystyle \delta =(2\nu /\omega )^{1/2}}
λ
≫
l
≫
δ
。
{\displaystyle \lambda \gg l\gg \delta .}
あなた
(
x
、
t
)
{\displaystyle U(x,t)}
(
あなた
、
ヴ
)
{\displaystyle (u,v)}
c
{\displaystyle c}
δ
{\displaystyle \delta}
l
/
c
{\displaystyle l/c}
非定常非圧縮 境界層 方程式は
∂
あなた
∂
t
+
あなた
∂
あなた
∂
x
+
ヴ
∂
あなた
∂
ええ
−
ν
∂
2
あなた
∂
ええ
2
=
あなた
∂
あなた
∂
x
+
∂
あなた
∂
t
{\displaystyle {\frac {\partial u}{\partial t}}+u{\frac {\partial u}{\partial x}}+v{\frac {\partial u}{\partial y}}-\nu {\frac {\partial ^{2}u}{\partial y^{2}}}=U{\frac {\partial U}{\partial x}}+{\frac {\partial U}{\partial t}}}
ここで、右辺の項は境界層に課せられる圧力勾配に対応します。この問題は、 および を満たす流れ関数を使用して解くことができます。 定義 により 、 音波の速度場は非常に小さいため 、 などの
漸近 級数を導入することで、境界層方程式の解を正式に得ることができます。
ψ
{\displaystyle \psi }
u
=
∂
ψ
/
∂
y
{\displaystyle u=\partial \psi /\partial y}
v
=
−
∂
ψ
/
∂
x
.
{\displaystyle v=-\partial \psi /\partial x.}
U
{\displaystyle U}
ε
→
0
{\displaystyle \varepsilon \rightarrow 0}
u
=
ε
u
1
+
ε
2
u
2
+
⋯
{\displaystyle u=\varepsilon u_{1}+\varepsilon ^{2}u_{2}+\cdots }
ψ
=
ε
ψ
1
+
ε
2
ψ
2
⋯
{\displaystyle \psi =\varepsilon \psi _{1}+\varepsilon ^{2}\psi _{2}\cdots }
最初の近似では、
∂
u
1
∂
t
−
ν
∂
2
u
1
∂
y
2
=
−
ω
cos
k
x
ℜ
(
i
e
−
i
ω
t
)
.
{\displaystyle {\frac {\partial u_{1}}{\partial t}}-\nu {\frac {\partial ^{2}u_{1}}{\partial y^{2}}}=-\omega \cos kx\Re (ie^{-i\omega t}).}
壁で滑りなし条件を満たし 、次の ように近づく解は 次のように与えられる。
y
/
δ
=
0
{\displaystyle y/\delta =0}
U
{\displaystyle U}
y
/
δ
→
∞
{\displaystyle y/\delta \rightarrow \infty }
u
1
=
ℜ
[
cos
k
x
(
1
−
e
−
κ
y
)
e
−
i
ω
t
]
,
ψ
1
=
ℜ
[
cos
k
x
ζ
1
(
y
)
e
−
i
ω
t
]
{\displaystyle u_{1}=\Re \left[\cos kx\,(1-e^{-\kappa y})\,e^{-i\omega t}\right],\quad \psi _{1}=\Re \left[\cos kx\,\zeta _{1}(y)\,e^{-i\omega t}\right]}
どこで そして
κ
=
(
1
−
i
)
/
δ
{\displaystyle \kappa =(1-i)/\delta }
ζ
1
=
y
+
(
e
−
κ
y
−
1
)
/
κ
.
{\displaystyle \zeta _{1}=y+(e^{-\kappa y}-1)/\kappa .}
次の式は
∂
u
2
∂
t
−
ν
∂
2
u
2
∂
y
2
=
U
∂
U
∂
x
−
u
1
∂
u
1
∂
x
−
v
1
∂
u
1
∂
y
.
{\displaystyle {\frac {\partial u_{2}}{\partial t}}-\nu {\frac {\partial ^{2}u_{2}}{\partial y^{2}}}=U{\frac {\partial U}{\partial x}}-u_{1}{\frac {\partial u_{1}}{\partial x}}-v_{1}{\frac {\partial u_{1}}{\partial y}}.}
右辺の各項は2次式なので、周波数と の項になります。これらの 項は 、 の時間に依存しない強制に対応します 。この時間に依存しない部分のみに対応する解を見つけましょう。これは、 が 方程式 [6] を満たす場合に得られます。
ω
+
ω
=
2
ω
{\displaystyle \omega +\omega =2\omega }
ω
−
ω
=
0.
{\displaystyle \omega -\omega =0.}
ω
=
0
{\displaystyle \omega =0}
u
2
{\displaystyle u_{2}}
ψ
2
=
sin
2
k
x
ζ
2
(
y
)
/
c
{\displaystyle \psi _{2}=\sin 2kx\,\zeta _{2}(y)/c}
ζ
2
{\displaystyle \zeta _{2}}
2
δ
ζ
2
‴
=
1
−
|
ζ
1
′
|
2
+
ℜ
(
ζ
1
ζ
1
″
)
{\displaystyle 2\delta \zeta _{2}'''=1-|\zeta _{1}'|^{2}+\Re (\zeta _{1}\zeta _{1}'')}
レイリーストリーミング
ここでプライムは に関する微分を表す。 壁面の境界条件は 、 が 有限でなければならないことを意味する。上記の式を2回積分すると、
y
.
{\displaystyle y.}
ζ
(
0
)
=
ζ
′
(
0
)
=
0.
{\displaystyle \zeta (0)=\zeta '(0)=0.}
y
/
δ
→
∞
{\displaystyle y/\delta \rightarrow \infty }
ζ
2
{\displaystyle \zeta _{2}}
ζ
2
′
=
3
8
−
1
8
e
−
2
y
/
δ
−
e
−
y
/
δ
[
sin
y
δ
+
1
4
cos
y
δ
+
y
4
δ
(
sin
y
δ
−
cos
y
δ
)
]
.
{\displaystyle \zeta _{2}'={\frac {3}{8}}-{\frac {1}{8}}e^{-2y/\delta }-e^{-y/\delta }\left[\sin {\frac {y}{\delta }}+{\frac {1}{4}}\cos {\frac {y}{\delta }}+{\frac {y}{4\delta }}\left(\sin {\frac {y}{\delta }}-\cos {\frac {y}{\delta }}\right)\right].}
なので 、 結果として となる。 したがって、境界の端では、振動運動に重ね合わせた定常流体運動が発生する。この速度強制により、境界層の外側に定常流動運動が発生する。興味深い結果は、 が に依存しないため 、境界層の外側で発生する定常流動運動は粘性にも依存しないということである。ただし、定常流動運動の発生源は粘性境界層である。
y
/
δ
→
∞
{\displaystyle y/\delta \rightarrow \infty }
ζ
′
(
∞
)
=
3
/
8
{\displaystyle \zeta '(\infty )=3/8}
v
2
(
x
,
∞
,
t
)
=
(
3
/
8
c
)
sin
2
k
x
.
{\displaystyle v_{2}(x,\infty ,t)=(3/8c)\sin 2kx.}
v
2
(
∞
)
{\displaystyle v_{2}(\infty )}
ν
{\displaystyle \nu }
外側の定常流動非圧縮運動は、問題の幾何学的形状に依存します。2つの壁があり、1つはとの場合 、 解は
y
=
0
{\displaystyle y=0}
y
=
2
h
{\displaystyle y=2h}
ψ
2
=
3
16
c
sin
2
k
x
[
−
(
y
−
h
)
+
(
y
−
h
)
3
/
h
2
]
{\displaystyle \psi _{2}={\frac {3}{16c}}\sin 2kx\,[-(y-h)+(y-h)^{3}/h^{2}]}
これは、図に示すように、反対方向に回転する渦の周期的な配列に対応します。
原因: 流体中の音響吸収による体積力
音響流は非線形効果である。
[7]
速度場を振動部分と定常部分に分解することができる 。振動部分は 音によるもので、定常部分は音響流速度(平均速度)である。 ナビエ・ストークス方程式は 音響流速度について次式を意味する。
u
=
v
+
u
¯
{\displaystyle {u}=v+{\overline {u}}}
v
{\displaystyle v}
ρ
¯
∂
t
u
¯
i
+
ρ
¯
u
¯
j
∂
j
u
¯
i
=
−
∂
p
¯
i
+
η
∂
j
2
u
¯
i
−
∂
j
(
ρ
v
i
v
j
¯
/
∂
x
j
)
.
{\displaystyle {\overline {\rho }}{\partial _{t}{\overline {u}}_{i}}+{\overline {\rho }}{\overline {u}}_{j}{\partial _{j}{\overline {u}}_{i}}=-{\partial {\overline {p}}_{i}}+\eta {\partial _{j}^{2}{\overline {u}}_{i}}-{\partial _{j}}({\overline {\rho v_{i}v_{j}}}/{\partial x_{j}}).}
定常流は、右側に現れる 定常体積力から生じます。この力は、 乱流における レイノルズ応力 と呼ばれるものの関数です。レイノルズ応力は音の振動の振幅に依存し、体積力はこの音の振幅の減少を反映します。
f
i
=
−
∂
(
ρ
v
i
v
j
¯
)
/
∂
x
j
{\displaystyle f_{i}=-{\partial }({\overline {\rho v_{i}v_{j}}})/{\partial x_{j}}}
−
ρ
v
i
v
j
¯
{\displaystyle -{\overline {\rho v_{i}v_{j}}}}
この応力は速度振幅において非線形( 2 次 )であることがわかります。速度振幅が変化する場所でのみ、応力はゼロではありません。流体の速度が音によって のように振動する場合 、2 次非線形性により に比例する一定の力が発生します
。
ϵ
cos
(
ω
t
)
{\displaystyle \epsilon \cos(\omega t)}
ϵ
2
cos
2
(
ω
t
)
¯
=
ϵ
2
/
2
{\displaystyle \scriptstyle {\overline {\epsilon ^{2}\cos ^{2}(\omega t)}}=\epsilon ^{2}/2}
音響流速度の大きさのオーダー
たとえ粘性が音響流の原因であったとしても、境界付近の音響流の場合、結果として生じる流速から粘性の値は消えます。
流速の大きさの順序は以下のとおりです。 [8]
U
∼
−
3
/
(
4
ω
)
×
v
0
d
v
0
/
d
x
,
{\displaystyle U\sim -{3}/{(4\omega )}\times v_{0}dv_{0}/dx,}
音の振動速度 に伴って 壁の境界に沿って流れます。流れは音の振動が減少する方向(振動の節)に向けられます。
v
0
{\displaystyle v_{0}}
x
{\displaystyle x}
静止半径aの振動泡 [9] の近くでは、その半径は相対振幅 (または )で脈動し、その質量中心も相対振幅 (または )で周期的に移動する。位相シフトで
ϵ
=
δ
r
/
a
{\displaystyle \epsilon =\delta r/a}
r
=
ϵ
a
sin
(
ω
t
)
{\displaystyle r=\epsilon a\sin(\omega t)}
ϵ
′
=
δ
x
/
a
{\displaystyle \epsilon '=\delta x/a}
x
=
ϵ
′
a
sin
(
ω
t
/
ϕ
)
{\displaystyle x=\epsilon 'a\sin(\omega t/\phi )}
ϕ
{\displaystyle \phi }
U
∼
ϵ
ϵ
′
a
ω
sin
ϕ
{\displaystyle \displaystyle U\sim \epsilon \epsilon 'a\omega \sin \phi }
壁から遠い [10] 流れの起源から遠い( 音響パワー、 動粘性、 音速)。流れの起源に近づくにつれて、速度はの平方根に比例します 。
U
∼
α
P
/
(
π
μ
c
)
{\displaystyle U\sim \alpha P/(\pi \mu c)}
P
{\displaystyle P}
μ
{\displaystyle \mu }
c
{\displaystyle c}
P
{\displaystyle P}
生物種、例えば接着細胞でさえも、音波にさらされると音響流を呈することが示されています。表面に接着した細胞は、表面から離れることなく、mm/sオーダーの音響流を生成できます。 [11]
参照
参考文献
^ レイリー、L. (1884)。クントの管内で観察される空気の循環と、それに関連する音響上の諸問題について。ロンドン王立協会哲学論文集、175、1-21。
^ http://lmfa.ec-lyon.fr/spip.php?article565&lang=en のビデオをご覧ください
^ Wan, Qun; Wu, Tao; Chastain, John; Roberts, William L.; Kuznetsov, Andrey V.; Ro, Paul I. (2005). 「振動する圧電バイモルフによって確立された狭いチャネル内の音響ストリーミングによる強制対流冷却」。 流れ 、 乱 流、燃焼 。74 (2): 195–206。CiteSeerX 10.1.1.471.6679。doi : 10.1007 /s10494-005-4132-4。S2CID 54043789 。
^ Nama, N., Huang, PH, Huang, TJ, および Costanzo, F., 振動する鋭いエッジの周りの音響流パターンの調査, Lab on a Chip, Vol. 14, pp. 2824-2836, 2014
^ Salari, A.; Appak-Baskoy, S.; Ezzo, M.; Hinz, B.; Kolios, MC; Tsai, SSH (2019) 細胞と踊る: 振動する細胞によって生成される音響マイクロフロー。https://doi.org/10.1002/smll.201903788
^ Landau, LD, & Lifshitz, EM (2000). 流体力学(理論物理学講座、第6巻)。
^ サー・ ジェームズ・ライトヒル (1978)「音響ストリーミング」、61、391、音響振動ジャーナル
^ Squires, TM & Quake, SR (2005) マイクロ流体工学: ナノリットルスケールの流体物理学、Review of Modern Physics、第 77 巻、977 ページ
^ Longuet-Higgins, MS (1998). 「振動する球状泡からの粘性流」 Proc. R. Soc. Lond. A . 454 (1970): 725–742. Bibcode :1998RSPSA.454..725L. doi :10.1098/rspa.1998.0183. S2CID 123104032.
^ Moudjed, B.; V. Botton; D. Henry; Hamda Ben Hadid; J.-P. Garandet (2014-09-01). 「音響ストリーミングジェットのスケーリングと次元解析」 (PDF) . 流体物理学 . 26 (9): 093602. Bibcode :2014PhFl...26i3602M. doi :10.1063/1.4895518. ISSN 1070-6631.
^ Salari, A.; Appak-Baskoy, S.; Ezzo, M.; Hinz, B.; Kolios, MC; Tsai, SSH (2019) 細胞と踊る: 振動する細胞によって生成される音響マイクロフロー。https://doi.org/10.1002/smll.201903788