[論文レビュー] Flow field from transient bubble oscillation in a narrow gap: numerical simulations and effect on biological cells

Milad Mohammadzadeh, Fenfang Li+1 more|arXiv (Cornell University)|May 13, 2016

Ultrasound and Cavitation Phenomena被引用数 7

ひとこと要約

本研究では、細い円形ギャップ内でのレーザー誘発キャビテーションバブルが生成する一時的流れ場を3次元CFDシミュレーションでモデル化し、壁近くで流れの逆転および強い境界層せん断応力が生じることを明らかにした。結果は、バブルの膨張および収縮に伴い、高いせん断速度（約2×10⁶ s⁻¹）が生じるため、赤血球が引き伸ばされる実験的観察を説明している。一方、より大きな細胞は主に流れの慣性により運ばれる。

ABSTRACT

The flow driven by a rapidly expanding and collapsing cavitation bubble in a narrow cylindrical gap is studied with the volume of fluid method. The simulations reveal a developing plug flow during the early expansion followed by flow reversal at later stages. An adverse pressure gradient leads to boundary layer separation and flow reversal, causing large shear stress near the boundaries. Analytical solution to a planar pulsating flow shows qualitative agreement with the CFD results. The shear stress close to boundaries has implications to deformable objects located near the bubble: experiments reveal that thin, flat biological cells entrained in the boundary layer become stretched, while cells with a larger cross-section are mainly transported with the flow.

研究の動機と目的

• 一時的キャビテーションバブルが細い円形ギャップ内で生成する複雑な3次元流れ構造を理解すること。
• 境界層および流れの逆転が近接する生物学的細胞にせん断応力を及ぼす役割を調査すること。
• 理想化されたポテンシャル流れモデルとマイクロ流体キャビテーション応用における現実的で粘性のある壁付き流れの間のギャップを埋めること。
• 定量的流れ場解析を通じて、実験的観察された赤血球の引き伸ばしを説明すること。

提案手法

• 細い円形ギャップ内での脈動するバブルの気液界面を追跡するため、体積フラッド（VOF）法を用いた数値シミュレーション。
• 時間依存の膨張および収縮ダイナミクスを示すレーザー誘発の一時的バブルのシミュレーション。
• 粘性で非定常な流れをモデル化するため、軸対称かつ非圧縮性流れの仮定に基づいたナビエ＝ストークス方程式の適用。
• 流れの逆転およびせん断応力の傾向を検証するため、平面脈動流れに対する解析解とCFD結果を比較。
• 変形可能な生物学的細胞（例：赤血球およびRKO細胞）の実験データを組み込み、流れの特徴と観察された細胞行動の相関を分析。
• 壁近くおよびコア流れ領域におけるせん断ひずみ率および変形度の定量的評価。

実験結果

リサーチクエスチョン

• RQ1 キャビテーションバブルの一時的膨張および収縮が、狭いギャップ内でどのように複雑な流れ構造を生成するか？
• RQ2 境界層および逆圧勾配が、壁近くでの流れの逆転および渦度を引き起こす役割は何か？
• RQ3 生じるせん断応力は、赤血球のような薄く平たい生物学的細胞を変形させるのに必要な値と比べてどの程度か？
• RQ4 なぜ一部の細胞は顕著な引き伸ばしを受けるのに対し、他の細胞は主に流れの慣性によって運ばれるのか？
• RQ5 平面脈動流れの簡略化された解析モデルは、3次元CFDシミュレーションで観察された主要な特徴をどの程度再現できるか？

主な発見

• 初期のバブル膨張段階では発達するプラグ流れが観察され、その後壁近くの逆圧勾配により流れの逆転が生じる。
• 遅い膨張および初期の収縮段階で境界層の剥離が発生し、制限された表面近くに強い渦度および循環領域が生成される。
• 壁近くでの最大せん断ひずみ率は約2×10⁶ s⁻¹に達し、実験的に観察された赤血球の引き伸ばしを説明するのに十分である。
• 境界層に位置する細胞は顕著なせん断誘発変形を経験するが、コア流れに位置する細胞は低減の変形を伴いながら主に流れに従って運ばれる。
• 平面脈動流れに対する解析解はCFD結果と定性的に一致しており、流れの逆転およびせん断応力生成の物理的解釈を支持する。
• 境界層外側のせん断応力は顕著に低く（1 µm厚さで約10％）、中心に位置する細胞では変形が最小限に抑えられる。