[論文レビュー] Numerical study of 3D vesicles under flow: discovery of new peculiar behaviors

Thierry Biben, Alexander Farutin+1 more|ArXiv.org|Dec 23, 2009

Blood properties and coagulation被引用数 4

ひとこと要約

本研究は、せん断流れ下における赤血球（RBC）のモデルとしての3次元数値解析を実施し、2つの新規な挙動を明らかにした。1つは、せん断速度が増加するにつれて、ふらつき呼吸（VB）モード領域が著しく広がることであり、もう1つは、たるみ（TB）およびVB振動周期がせん断速度に強く依存しない、強固な独立性を示すことである。境界積分シミュレーションと解析的モデルを用いて、著者らは、大きな変形下では高次調波が力学的挙動において重要になることを示し、強い流れ条件下では古典的ケラーサクラック理論の有効性が疑問視されることを示した。

ABSTRACT

The study of vesicles under flow, a model system for red blood cells (RBCs), is an essential step in understanding various intricate dynamics exhibited by RBCs in vivo and in vitro. Quantitative 3D analyses of vesicles under flow are presented. The regions of parameters to produce tumbling (TB), tank-treating (TT), vacillating-breathing (VB) {and even Kayaking (K)} modes are determined. New qualitative features are found: (i) a significant widening of the VB mode region in parameter space upon increasing shear rate $\dotγ$ and (ii) a striking robustness of period of TB and VB with $\dotγ$. Analytical support is also provided. These findings shed new light on the dynamics of blood flow.

研究の動機と目的

• せん断流れ下における3次元ヘモリティック細胞の動的挙動を体系的・包括的に調査すること。
• たるみ（TB）、タンクトレッド（TT）、ふらつき呼吸（VB）、カイキング（K）モードのパラメータ領域を特定・特徴づけること。
• 特に、VBモードのせん断速度依存的拡大と、TBおよびVB振動周期のせん断速度不変性といった、これまで観察されていなかったヘモリティック細胞の挙動を解明すること。
• 強い流れや大きな変形下におけるケラーサクラック（KS）理論の限界を評価すること。
• 観測された数値結果を説明する解析的理論を構築・検証すること、特に形状動的挙動における高次調波の役割を解明すること。

提案手法

• 単純せん断流れ下におけるヘモリティック細胞の Stokes 方程式を境界積分（BI）法で数値的に解くシミュレーション。
• 膜の弾性を表現するため、ヘルドリッヒの曲げ力と、曲げ剛性定数 $ \kappa $、体積不変制約を組み込む。
• ヘモリティック細胞の形状変化を特徴付ける次元なし変形パラメータ $ C_{\bar{k}} $ を用いる。
• ケラーサクラック（KS）理論との比較を通じて、ヘモリティック細胞の主軸の角度変化 $ \dot{\theta} $ を分析する。
• 大変形下で高次調波（例：$ \cos(4\theta) $）が出現するメカニズムを説明するための摂動的解析理論を適用する。
• 長時間シミュレーションにより、過渡的挙動および安定状態モードへの緩和挙動を評価し、特にTT-VBやTB-Kの遷移線付近の挙動を重点的に分析する。

実験結果

リサーチクエスチョン

• RQ1 パラメータ空間におけるふらつき呼吸（VB）モード領域の広がりは、せん断速度 $ \dot{\gamma} $ の増加に伴いどのように変化するか？
• RQ2 なぜ、たるみ（TB）およびVBモードの振動周期が、過去のシミュレーションとは異なり、せん断速度にほとんど依存しないのか？
• RQ3 強い流れまたは大きなヘモリティック細胞の変形下では、ケラーサクラック理論がどの程度崩壊するのか？
• RQ4 高次調波（例：$ \cos(4\theta) $）は、非常に変形したヘモリティック細胞の挙動において、どのような役割を果たすのか？
• RQ5 モード遷移線（例：TT-VB）付近では、ヘモリティック細胞が安定状態への緩和に要する時間はどの程度で、その時間スケールは何か？

主な発見

• パラメータ空間におけるふらつき呼吸（VB）モード領域は、せん断速度 $ \dot{\gamma} $ の増加に伴い著しく広がり、これは過去の理論やシミュレーションでは捉えられていなかった挙動である。
• たるみ（TB）およびVBモードの振動周期は、せん断速度の10倍の範囲にわたりほぼ一定であり、$ \dot{\gamma} $ に強く依存しない強固な動的挙動を示している。
• ケラーサクラック理論は、変形が小さい場合または粘度比 $ \lambda $ が大きい場合にはTBおよびTT動的挙動を正確に記述できるが、高 $ \dot{\gamma} $ または中程度の $ \lambda $ では顕著なずれが生じる。
• 高次調波、特に $ \cos(4\theta) $ は、変形度 $ C_{\kappa} \geq 8 $ の大変形域において、形状動的挙動に顕著に現れ、TB-VB や TB-K の分岐的遷移の始まりを示す兆候となる。
• TT-VB遷移線付近での安定状態への緩和時間は、$ 30 \eta_{in}R^3 / \kappa $ のオーダーであり、約125回の振動に相当し、緩和が遅い過渡的状態が続くことが示された。
• 本研究で開発した解析的理論は、数値的観察をうまく説明できており、特に高次項の出現と振動周期の安定性の両方を的確に解釈している。