Deformable bodies in a 3-dimensional viscous flow: Vorticity-Stream vector formulation

技術要約：渦度・ストリームベクトル定式化による3次元粘性流中の変形体のシミュレーション

問題提起
三次元（3D）の非圧縮性粘性流と、変形体または弾性体（小胞、液滴、細胞など）との相互作用をシミュレーションすることは、大きな課題を伴う。格子ボルツマン法（LBM）、境界積分法、および粒子ベースの手法である平滑化散逸粒子動力学（SDPD）を含む既存の手法は、しばしば高い計算複雑性、複雑な数値実装、あるいは複雑な界面力学の取り扱いの困難さに直面する。渦度・流線関数定式化とフェーズフィールドモデルを用いた2次元シミュレーションは堅牢であることが証明されているが、低レイノルズ数における3次元変形界面のための、これと同等に軽量で古典的な理論に基づいたフレームワークは不足していた。具体的には、カンハム–ヘルフリッヒ（膜の曲げ）またはカーン–ヒリアード（表面張力）エネルギーによって支配される3次元界面のための、渦度に基づく定式化に関する文献には空白が存在する。

手法
著者らは、特に低レイノルズ流用に調整された、渦度・ストリームベクトルアプローチに基づく新しい数値定式化を導入する。この手法は、2次元の流線関数を3次元のベクトルポテンシャル（$\vec{\psi}$）へと一般化し、速度場を $\vec{v} = \nabla \times \vec{\psi}$ と定義する。

手法の核となるのは、以下の3つの主要な観測量に関する結合系の解法である：

1. 界面化学ポテンシャル ($\mu$): フェーズフィールド秩序パラメータ ($\phi$) に対する自由エネルギー ($F$) の汎関数微分として導出される。著者らは、カンハム–ヘルフリッヒ・エネルギー（膜の曲げ用）とカーン–ヒリアード・エネルギー（表面張力用）の両方を用いて、このフレームワークを実証している。
2. 流体渦度 ($\vec{\omega}$): $\vec{\omega} = \nabla \times \vec{v}$ と定義される。
3. ストリームベクトル ($\vec{\psi}$): ゲージ条件 $\nabla \cdot \vec{\psi} = 0$ を満たすベクトルポテンシャル。

支配方程式は、一連の結合ポアソン問題へと簡略化される：

• フェーズフィールド $\phi$ の進化は、化学ポテンシャルと速度場を含む移流拡散方程式によって支配される。
• 渦度方程式は、ストークス方程式（または慣性項を含むナビエ・ストークス方程式）のカールから導かれ、ソース項がフェーズ勾配と化学ポテンシャル勾配の外積 $\nabla^2 \vec{\omega} = \frac{1}{\eta} (\nabla \phi) \times (\nabla \mu)$ に依存するポアソン方程式となる。
• ストリームベクトルは、$\nabla^2 \vec{\psi} = -\vec{\omega}$ を解くことで得られる。

数値実装
システムは、固定された、適合していない（unfitted）立方体格子上の有限差分法を用いて解かれる。界面はフェーズフィールド法を介して拡散領域として扱われるため、複雑なメッシュ変形や埋め込み境界法を必要としない。時間の進化には、フェーズフィールドに対して前向きオイラー積分を用い、ポアソン方程式に対しては、スタガード（逐次分離）スキームを利用した反復ソルバーを用いる。コードはPythonで実装されており、オープンソースのリポジトリとして公開されている。

主な貢献

• 3Dソルバーの簡素化: この定式化は、流体ソルバーを一連の結合ポアソン問題へと変換し、圧力の直接計算を回避し、構成によって厳密な非圧縮性を保証する。
• 界面エネルギーの統一フレームワーク: この手法は、異なる自由エネルギーモデル（膜の曲げのためのカンハム–ヘルフリッヒ、および表面張力のカーン–ヒリアード）を同一の数値構造内にシームレスに統合しており、膜力学の容易なカスタマイズを可能にする。
• 診断的観測量: 著者らは、未摂動の流れに対する渦度およびストリームベクトルの二乗揺らぎに基づく新しい観測量を提案している。これらの指標は、形状変形や形態変化に対して非常に高い感度を示すことが示されている。

結果
本手法は、ニュートン流体中の単一小胞および液滴のシミュレーションを通じて検証された：

• ポアズイユ流れ: シミュレーションは、赤血球（RBC）の典型的な軸対称形状（「パラシュート」型や「スリッパー」型を含む）を正確に再現し、既存の文献と一致した。診断的観測量は、緩和ダイナミクスを効果的に追跡し、形状変化（例：ディンプルの消失）に対応する明確なピークを示し、物体が定常状態に達するとプラトーに達することを示した。
• クエット流れ: 本フレームワークは、横方向移動（lateral migration）、すなわち小胞が流れの方向に対して垂直方向に移動する非慣性揚力を捉えた。シミュレーションは、平衡位置と移動ダイナミクスが、内部流体と外部流体の間の粘度コントラスト ($\zeta$) に系統的に依存することを実証した。粘度コントラストが高いほど、移動壁に向かうドリフトがより顕著になった。
• 応力解析: 本手法は面内のせん断応力分布の計算に成功し、より高いせん断速度が応力の不均一性を誘発し、形態的不安定性を引き起こし得ることを明らかにした。

意義と主張
本論文は、この定式化が、既存のメゾスコピックな手法や高複雑度のソルバーに対する、軽量で柔軟かつ計算効率の高い代替案を提供すると主張している。運動論的理論（LBMのような）ではなく、古典的な流体力学の原理に依拠することで、この手法は連続体メカニクスとの整合性の証明を必要としない。

著者らは、現在の研究がニュートン流体中の単一体のダイナミクスに焦点を当てているものの、このフレームワークが拡張性に優れていることを強調している。この手法は、体積力、慣性効果、粘弾性媒体、および多体相互作用を自然に取り込むことができる。2次元近似では失われがちな3次元面内せん断ダイナミクスを捉える能力により、本手法は細胞生物学やソフトマター物理学の研究における貴重なツールとしての地位を確立している。コードのオープンソース化は、この新興分野におけるベンチマークとさらなる発展を促進することを目的としている。