Phenomenon in physics
音響流は、高振幅の 音響 振動の吸収によって駆動される流体内の定常流です。この現象は、音響放射体の近く、または クント管内 の定在波で観察できます。音響流は 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}}
レイリーストリーミング 速度場 に対応する平面定在音波を考えます 。問題の特性(横方向)次元を とします 。ここで説明した流れ場は、 非粘性流れ に対応します。しかし、固体壁の近くでは粘性効果が重要になります。その場合、厚さ(または侵入深さ) の境界層が存在します 。レイリー流は、近似 で視覚化するのが最善です。 の場合と同様に 、速度成分は よりもはるかに小さくなります。さらに、境界層内の特性時間スケールは 、音響時間スケール と比較して 非常に大きくなります( が小さいため) 。これらの観察結果は、境界層内の流れは非圧縮性と見なせることを示唆しています。 U ( x , t ) = v 0 cos k x cos ω t = ε cos k x ℜ ( e − i ω t ) {\displaystyle U(x,t)=v_{0}\cos kx\cos \omega t=\varepsilon \cos kx\Re (e^{-i\omega t})} k = 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 .} U ( x , t ) {\displaystyle U(x,t)} ( u , v ) {\displaystyle (u,v)} c {\displaystyle c} δ {\displaystyle \delta } l / c {\displaystyle l/c}
非定常非圧縮 境界層 方程式は
∂ u ∂ t + u ∂ u ∂ x + v ∂ u ∂ y − ν ∂ 2 u ∂ y 2 = U ∂ U ∂ x + ∂ U ∂ 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}''),} レイリーストリーミング ここで、プライムは に関する微分を表す。 壁における境界条件は 、 As が有限でなければならないことを意味する。上式を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つは で 、 もう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}}}}
この応力は速度振幅に対して非線形( 二次 )であることがわかります。この応力は、速度振幅が変化する部分でのみゼロではありません。流体の速度が音波によって のように振動すると 、二次非線形性によって に比例する一定の力が発生します 。 ϵ cos ( ω t ) {\displaystyle \epsilon \cos(\omega t)} ϵ 2 cos 2 ( ω t ) ¯ = ϵ 2 / 2 {\displaystyle {\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]
アプリケーション 音響流に関する研究は、特に粒子操作を中心に多くの効果的な用途を示していますが、商業化への応用はほとんどの用途において初期段階にあります。マイクロ流体工学においては、細胞操作と細胞選別に利用できます。 [12] [13] これらの用途には、細胞操作と 細胞選別 、薬物送達、反応物の均質化などが含まれます。音響流は、細胞膜透過性を高めるための ソノポレーション にも関連しています。音響流は膜プロセスにも利用されており、膜の汚れを抑制し、粒子の回収率を向上させることができます。 [14] 他の用途においても、バイオフィルムを制御できます。 [15]
参照
参考文献 ^ レイリー、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, and 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 ^ ^ ^ Barrio-Zhang, Andres; Anandan, Sudharshan; Deolia, Akshay; Wagner, Ryan; Warsinger, David M.; Ardekani, Arezoo M. (2024). 「音響的に強化された多孔質媒体は、ろ過性能を劇的に向上させる」『 分離・精製技術 』 342 126972. doi :10.1016/j.seppur.2024.126972. ^ Lin, Fangfei; Yuan, Songmei; Ji, Pengzhen; Xu, Weixian (2023). 「超音波による細菌バイオフィルム形成の制御:オートインデューサー2の役割と音響流の有限要素解析」. Ultrasound in Medicine & Biology . 49 (9): 2191– 2198. doi :10.1016/j.ultrasmedbio.2023.06.016.