A Modified Suspension-Balance Model for Deformable Particle Suspensions: Application to Blood Flows with Cell-Free Layer

本論文は、微小血管内における血液の流れを効率的にシミュレートするために流体力学的揚力を組み込んだ修正懸濁平衡モデルを提案しており、細胞フリー層の形成、ヘマトクリットおよび速度プロファイル、そしてファレウス効果およびファレウス・リンドクリスト効果といった主要な現象を捉えることに成功している。

要約

忙しい高速道路を想像してみてください。そこでは、車は実は小さくて柔らかい赤血球であり、道路は微小な血管です。通常の高速道路では、交通量は均等に分散していると予想されるでしょう。しかし、私たちの体内では、これらの「車」には奇妙な癖があります。彼らは壁の近くにいることを嫌うのです。彼らは道路の中央に集まることを好み、路肩のすぐ横には、がらんとした空のレーンを残していきます。

この空のレーンのことは、**細胞フリー層（Cell-Free Layer: CFL）**と呼ばれます。これは、血液がより速く、摩擦を少なく流れるのを助ける、私たちの血液流における極めて重要な特徴です。

問題点：古い地図には曲がり角が欠けていた

科学者たちは、血液の流れをシミュレートするためのコンピュータモデルを長年構築しようとしてきました。彼らは「懸濁液バランスモデル（Suspension Balance Model: SBM）」と呼ばれるものを使用しています。このモデルを、交通シミュレーション・ソフトウェアだと考えてください。

旧バージョンのこのソフトウェアは、細胞同士がぶつかり合うことで、車が道路の中央へと移動することを予測することには長けていました。しかし、なぜ車がこれほどまでに壁を離れたがるのかという「理由」を説明することには失敗していました。それは、あの「細胞フリーのレーン」を端の部分に作り出すことができなかったのです。それはまるで、車が動いていることは分かっているものの、車が自ら進んで縁石を避けていることまでは理解できていないGPSのようなものでした。

解決策：新しい「プッシュ」ボタン

Hugo Castillo-SánchezとLeonardo Liu率いる著者たちは、このソフトウェアを修正することに決めました。彼らは、赤血球は柔らかい（変形能がある）ため、壁に近づきすぎると特別な「目に見えない力」を生み出すことに気づいたのです。

彼らはこれを**揚力（Lift Force）**と呼んでいます。

• 比喩： プールのサイド付近で泳いでいるところを想像してください。あなたが動くと、水があなたを壁から少し押しやります。赤血球の場合、彼らが壁をすり抜ける際に形を変えながら押しつぶされるため、この「押し」は非常に強力になります。
• 修正方法： チームはこの「揚力」をコンピュータモデルに追加しました。彼らは**修正懸濁液バランスモデル（Modified Suspension Balance Model: MSBM）**を作り上げました。今や、このソフトウェアはただ車を見守るだけでなく、水が泳ぎ手を押し出すように、能動的に細胞を壁から遠ざけることができるようになったのです。

シミュレーションを実行した結果はどうなったか？

彼らがコンピュータの中でこの「揚力」をオンにしたとき、結果は劇的に変化しました。

1. 空のレーンの出現： シミュレーションは、実生活で見られるような、壁の近くのクリアなゾーン（CBF）を正常に作り出しました。
2. 中央での交通渋滞： 赤血球は中央に密集し、高密度のコア（核）を作り出しました。
3. 流れの形状： 細胞が中央に集まり、端がクリアになったため、血液は（通常の川のように）滑らかな曲線を描いて流れるのではなく、平らな上面を持つ、固体のプラグ（栓）やピストンのように流れるようになりました。これは、実際の微小血管内で起きている現象そのものです。

新しいモデルのテスト

チームは単に推測したのではなく、実世界のデータや他の複雑なシミュレーションと比較して、新しいモデルをテストしました。

• タイムトラベル（時間の経過）： 彼らは「空のレーン」がどのように形成されるかを観察しました。最初は細胞があらゆる場所に存在していましたが、次第に「揚力」が細胞を壁から押し退け、レーンがクリアになりました。これは、高速カメラによる実験で見られる速度や挙動と一致していました。
• 「ファラエウス効果（Fåhræus Effect）」： これは単純な観察に対する専門的な名称です。血液は予想よりも小さな管の中を速く流れ、中央の細胞濃度は出口での濃度とは異なります。彼らの新しいモデルは、これを完璧に予測しました。
• 「ファラエウス・リンドクイスト効果（Fåhræus-Lindqvist Effect）」： これは、血液が非常に小さな管の中を流れるとき、「より薄く（粘り気が少なく）」なるという観察結果です。彼らのモデルもこれも捉えており、壁の近くの空のレーンが摩擦を軽減し、血液の流れを容易にしていることを示しました。

まとめ

この論文は、シンプルな「押し（揚力）」をコンピュータモデルに加えることで、血液が微小血管内でどのように振る舞うかを正確にシミュレートできるようになったと主張しています。

• それができること： 細胞フリー層の形成、プラグ状の流れ、そして血液の流れを効率的にしている有名な「ファラエウス」効果を捉えることができます。
• まだできていないこと（現時点では）： 著者らは、非常に大きな管（40マイクロメートル以上）において、モデルが細胞を押しすぎる傾向があると認めています。彼らは、密集している際に細胞が互いに壁からの影響を「遮蔽（シールド）」する仕組みが、まだモデルに含まれていないためではないかと推測しています。彼らは将来の研究でこれを修正する予定です。

要するに、彼らは、赤血球が単なる受動的な乗客ではなく、高速道路をクリアに保つために自ら壁から離れようとする「能動的な泳ぎ手」であることを理解した、より優れたデジタルツインを構築したのです。

技術概要

技術要約：変形粒子懸濁液のための修正懸濁平衡モデル

問題提起
拘束された幾何学的形状内における赤血球（RBC）およびその他の変形粒子の流れをモデリングすることは、固有の非線形性、多体相互作用、および非ニュートン流体的なレオロジー挙動により、大きな課題となっている。微小循環における重要な現象は、細胞欠乏層（CFL）の形成であり、これは壁付近の赤血球が枯渇した領域である。この層は、見かけの粘性を低く維持し、血小板の辺縁化を促進するために生理学的に不可欠である。細胞スケールの直接的な計算モデルはCFL形成を捉えることができるが、血管ネットワークや臓器スケールのシミュレーションにおいては計算コストが極めて高い。一方、従来の懸濁平衡モデル（SBM）のような既存の連続体モデルは、粒子を壁から遠ざける変形性に起因する流体力学的リフト力を十分に考慮していないため、CFLを捉えることができないことが多い。

手法
著者らは、壁と赤血球の間の流体力学的相互作用を考慮するために、特定のリフト力項を従来のSBMフレームワークに付加した「修正懸濁平衡モデル（MSBM）」を提案している。

• 理論的枠組み: 本モデルは、懸濁液バルク相および粒子相に関する標準的な質量保存および運動量保存方程式を保持している。核心となる修正は、粒子相の運動量バランスにある。従来のSBMは、粒子相応力のダイバージェンス（$\nabla \cdot \langle \Sigma \rangle_p$）および妨げられた沈降ドラグによって駆動される粒子移動を考慮している。MSMSでは、変形ベシクルおよび赤血球の実験的観察に基づいた、体積平均化されたリフト力項 $\phi L_\perp$ を導入している。
• リフト力の定式化: リフト力は、局所的なせん断速度（$\dot{\gamma}$）、壁からの距離（$h$）、および縮小体積パラメータ（$\nu$）の関数としてモデル化されている。その形式は $L_\perp \propto \frac{\dot{\gamma}}{(h+h_0)^\beta}$ であり、ここで $\beta$ は粒子依存性（具体的には赤血球）を反映するために調整されたべき乗則係数であり、$h_0$ は壁での特異点を回避するための幾何学的因子である。
• 数値実装: MSBMは、既存の SbmFoam ソルバーを拡張した SbmBloodFoam を用いてOpenFOAMに実装された。このソルバーは有限体積法（FVM）アプローチを利用している。連続体モデルの検証として、著者らは格子ボルツマン法（LB）と赤血球の粗視化変形カプセルモデルに基づいた、実験的に検証済みのマルチスケール血液流ソルバーを用いた。

主な結果
本研究は、微小血管チャネルおよびチューブ内の血液流をシミュレートする上でのMSBMの能力を実証しており、過渡的および定常的な挙動を捉えている。

1. CFLの形成と速度プロファイル: 従来のSBMが放物線状の速度プロファイルと壁面での非ゼロ濃度を予測するのに対し、MSBMは血液に特徴的なCFLの形成と、それに伴う「プラグフロー（固形物流）」のような速度プロファイルを正常に捉えている。モデルは、リフト力パラメータ（$\beta$）が増加するとCFLの厚さが増すことを示している。
2. パラメトリック感度:
   • バルクヘマトクリット（$\phi_b$）: $\phi_b$ が増加すると、文献通りCFLの厚さは減少する。
   • 最大充填濃度（$\phi_m$）: 赤血球の凝集または充填傾向を表す $\phi_m$ の値が高くなると、より厚いCFLと高い中心部ヘマトクリットピークが生じる。
3. 時間発展: モデルはCFLの過渡的な発達を捉えている。初期状態ではヘマトクリットは均一であるが、時間の経過とともに赤血球は中心部へと移動し、壁付近に細胞欠乏層、中心部に濃縮コアを形成して、最終的に定常状態に達する。
4. 理論および実験との検証:
   • 単純せん断流: モデルは、短時間スケールにおいてCFLの厚さの3乗（$\delta_{CFL}^3$）が無次元時間（$t\dot{\gamma}$）に対して線形スケーリングすることを再現しており、これは単一粒子移動実験で観察されるスケーリング則と一致している。
   • チャネル流: MSBMによる速度およびCFL厚さの予測は、様々なチャネル高さおよびバルクヘマトクリットにおけるLosserandらやSalameの実験データとよく一致している。
   • 管状流: モデルは、ファラー効果（排出ヘマトクリットが管内ヘマトクリットを上回る現象）およびファラー・リンドキスト効果（血管径の減少に伴う見かけの粘性の低下）を正常に捉えている。予測された見かけの粘度およびヘマトクリット比は、Priesらによる経験的相関関係、および特に直径40 µm以下の血管における実験データと良好な一致を示している。

意義および主張
本論文は、拘束条件下における高濃度変形粒子懸濁液の微細構造の不均一性と非ニュートン流動挙限を効率的に捉える、新しい連続体計算フレームワークを確立している。著者らは、変形性に起因するリフト力を組み込むことにより、MSBMが計算コストの高い離散粒子シミュレーションと、精度の低い連続体モデルとの間の溝を埋めるものであると主張している。

本研究は、効率的な、大規模なバイオメディカルおよび産業への応用に向けた一歩として提示されている。著者らは、モデルが小さな血管については良好に機能するものの、より大きな管（>40 µm）においてはCFLおよびファラー・リンドキスト効果をわずかに過大評価することを控えめに述べている。これは、現在のリフト力定式化に、局所的な粒子-壁間の潤滑相互作用によるリフト力の減少を考慮する「粒子スクリーニング効果」が欠けているためであると考えている。この限界は、現在のモデルの失敗ではなく、将来の研究課題として特定されている。