On the development of OpenFOAM solvers for simulating MHD micropolar fluid flows with or without the effect of micromagnetorotation

本論文は、外部磁場下での血液やフェロ流体などのマイクロポーラー流体の挙動を解析するため、マイクロマグネトロテーション（MMR）の影響を考慮した新しい OpenFOAM ソルバーを開発し、その高精度な検証と生体医学応用における可能性を示したものである。

要約

この論文は、**「血液の流れをより正確にシミュレーションするための、新しいコンピュータープログラム（ソルバー）」**の開発について書かれたものです。

専門用語を避け、日常のイメージを使って説明しましょう。

1. 血液は「ただの液体」ではない

まず、血液は水のように均一な液体ではありません。赤血球という「小さな硬い粒子」が水（血漿）の中に浮いている**「微粒子入り液体（マイクロポーラー流体）」**です。

• イメージ: 水に小さなビー玉を混ぜたような状態。
• 特徴: 普通の水と違い、このビー玉（赤血球）は、流れの中で**「自分自身で回転」**することができます。

2. 磁石の力で「回転」が止まる現象（MMR）

この研究の核心は、**「磁石」**の存在です。
MRI（磁気共鳴画像法）のように強い磁場をかけると、血液の中の赤血球（ヘモグロビンが鉄を含んでいるため）は磁石のように反応します。

• 普通の考え方（これまでの研究）: 磁場をかけると、電流が流れて「ローレンツ力」という力が働き、流れが少し変わるだろうと考えられていました。
• この論文の発見（MMR）: しかし、実はそれよりも重要な現象が起きていました。それは**「マイクロマグネトロテーション（MMR）」**と呼ばれるものです。
  • アナロジー: 磁場をかけると、赤血球たちは**「磁石の方向に揃って、回転することを拒否する」**ようになります。
  • 結果: 赤血球が「回転」しなくなると、血液全体の流れが**「かたまり」になり、速度が劇的に落ちます（最大 40% 減）。また、赤血球の回転自体も99.9% 以上停止**してしまいます。
  • 重要な点: 血液は電気をあまり通さないため、従来の「ローレンツ力」による影響はほとんど無視できるほど小さいですが、「回転が止まる（MMR）」という効果は非常に巨大です。

3. 開発された「魔法の道具」：2 つの新しいプログラム

研究者たちは、この複雑な現象を計算するために、オープンソースの CFD（数値流体力学）ソフト「OpenFOAM」を使って、2 つの新しいプログラムを作りました。

1. epotMicropolarFoam（エポット・マイクロポーラー・フォーム）:
   • 赤血球の「回転」は考慮するが、磁場による「回転の停止（MMR）」は考慮しない、基本的なプログラム。
2. epotMMRFoam（エポット・MMR・フォーム）:
   • ここが新しさです。磁場によって赤血球の回転がどう止まるか（MMR）まで計算できる、より高度なプログラム。

これらは、**「磁気熱療法（がん治療など）」や「標的薬物送達（薬を患部に届ける技術）」**などの医療技術に役立つと期待されています。

4. 実験結果：動脈と「風船」のような瘤（こぶ）

このプログラムを使って、実際の血管のようなシミュレーションを行いました。

• 単純な血管（3D 動脈）:
  • 磁場をかけると、赤血球が整列して回転しなくなり、血液の流れがスムーズに止まる（遅くなる）ことが確認されました。
• 動脈瘤（2D アニエリスム）:
  • イメージ: 血管の一部が風船のように膨らんでいる状態。ここには通常、**「渦（うず）」**ができて、血液がグルグルと循環しています。
  • 発見: 磁場をかけると、この**「渦」が驚くほど小さくなり、安定しました。**
  • 理由: 赤血球が回転しなくなったため、渦を作るエネルギーが失われ、血液が「静まり返った」状態になったからです。

5. まとめ：なぜこれがすごいのか？

これまでの研究では、「磁場をかけると血液の流れが変わる」というのは、電気の力（ローレンツ力）によるものだと考えられていましたが、実は**「赤血球の回転が磁場で止まること（MMR）」**が、流れを劇的に変える主役だったのです。

この新しいプログラムを使えば、医師や研究者は**「磁石を使って、どうすれば血液の流れを制御できるか」**を、より現実的にシミュレーションできるようになります。

一言で言うと：
「血液の中の小さな赤血球が、磁石の力で『回転するのをやめて』整列することで、流れが劇的に変わる現象を、初めて正確に計算できるプログラムを作ったよ！」という研究です。

技術概要

論文の技術的サマリー：OpenFOAM による MMR 効果を考慮した MHD 微極流体のシミュレーション

本論文は、磁気粒子（赤血球やフェロ流体など）を含む微極流体（Micropolar fluid）の磁気流体力学（MHD）シミュレーションを対象として、OpenFOAM 環境で開発された 2 つの新しいソルバ（epotMicropolarFoam および epotMMRFoam）の提案、実装、検証、および生物医学的応用に関する結果を報告したものである。特に、外部磁場下での「マイクロマグネトロテーション（MMR: Micromagnetorotation）」効果を考慮することの重要性と、それが血流の速度や微極回転に与える劇的な影響に焦点を当てている。

1. 背景と課題 (Problem)

• 微極流体理論: 微極流体は、母液中を独立して運動する微小な剛体粒子（赤血球など）を含む流体を記述する理論であり、角運動量保存則を考慮する必要がある。従来のナビエ - ストークス方程式にはない「カップル応力（couple stress）」と「スピン粘性（spin viscosity）」が導入される。
• MHD と血液: 血液は赤血球内のヘモグロビン（酸化鉄）が磁性粒子として振る舞うため、外部磁場を印加するとフェロ流体として扱える可能性がある。従来の MHD 研究では、ローレンツ力（電流と磁場の相互作用）が主要な力として扱われてきた。
• 未解決の課題: しかし、血液のような導電率が低い流体では、ローレンツ力の影響は微小である。一方で、磁性粒子の磁化ベクトルと外部磁場の不一致によって生じる「マイクロマグネトロテーション（MMR）」による磁気トルクの影響は、従来の MHD 研究や血液流体力学において十分に検討されてこなかった。
• 数値コードの欠如: 外部磁場下での微極流体、特に MMR 項を含む構成式（磁化方程式）を完全に結合したオープンソースの数値コードは存在しなかった。

2. 手法とソルバ開発 (Methodology)

著者らは、OpenFOAM の既存ソルバ（icoFoam, epotFoam, micropolarFoam）を基盤として、以下の 2 つの非定常ソルバを開発した。

2.1 開発されたソルバ

1. epotMicropolarFoam:
   • 従来の MHD 微極流体シミュレーション用。
   • MMR 項を含まない。
   • 運動量方程式に微極回転と渦度の差（$2\mu_r \nabla \times (\omega - \Omega)$）を源項として追加。
   • 内部角運動量方程式を求解。
2. epotMMRFoam:
   • 上記ソルバを拡張し、MMR 効果を含む。
   • Shizawa-Tanahashi モデルに基づく磁化方程式（構成式）を求解。
   • 内部角運動量方程式に磁気トルク項（$\mathbf{M} \times \mathbf{H}$）を追加。
   • 磁気流束密度 $\mathbf{B}$ を磁場強度 $\mathbf{H}$ と磁化 $\mathbf{M}$ の関数として計算。

2.2 数値アルゴリズム

• 圧力 - 速度結合: PISO (Pressure-Implicit with Splitting of Operators) アルゴリズムを使用。
• MHD 近似: 低磁気レイノルズ数近似（$Re_m \ll 1$）を採用。誘導方程式を無視し、電気ポテンシャル（Electric Potential）定式化を用いる。
• 離散化: 有限体積法（FVM）。
• 実装: C++ で実装され、運動量方程式、角運動量方程式、電気ポテンシャルのポアソン方程式を反復的に求解する構造となっている。

3. 検証と結果 (Validation and Results)

3.1 解析解との検証（ポアズイユ流れ）

• 設定: 2 次元および 3 次元のポアズイユ流れ（血管モデル）において、ヘマトクリット値（$\phi = 25\%, 45\%$）と磁場強度（1 T, 3 T, 8 T）を変化させてシミュレーション。
• 精度: 解析解との比較において、速度および微極回転の誤差は 2% 未満（多くのケースで 1% 以下）であり、ソルバの高精度が確認された。
• MMR の影響:
  • MMR 無視の場合: 血液の導電率が低いため、ローレンツ力の影響は極めて小さく、速度や微極回転への影響は無視できるレベルであった。
  • MMR 考慮の場合: 磁場強度が強いほど、赤血球の磁化によるトルクが内部回転を抑制する。
    • 速度の減少：最大で 40% まで低下。
    • 微極回転の減少：最大で 99.9% まで抑制（粒子が磁場方向に整列し、内部回転が「凍結」する状態）。
    • ヘマトクリット値が高いほど（粒子数が多いほど）、この効果は顕著になる。

3.2 3 次元動脈流シミュレーション

• 単純な円筒動脈モデルを用い、1 T と 5 T の磁場を印加。
• MMR を考慮しない場合、磁場の影響は微小（速度変化 1% 未満）。
• MMR を考慮すると、5 T で速度が約 27-39% 減少し、微極回転は 99% 以上抑制された。

3.3 2 次元動脈瘤（Aneurysm）流シミュレーション

• 166% および 200% の拡張率を持つ動脈瘤モデルでシミュレーション。
• 無磁場・MMR 無視: 瘤内部に大きな再循環渦（recirculation zones）が発生。
• MMR 考慮時:
  • 磁場強度の増加に伴い、再循環渦が強く抑制され、安定化する。
  • 渦の核心部が圧縮され、せん断が緩和される。
  • 8 T の磁場では、微極回転がほぼ完全に抑制され（99.9% 減少）、流れの構造が劇的に変化。
  • 速度プロファイルも圧縮され、最大で 42.9% の減少が見られた。

4. 主要な貢献 (Key Contributions)

1. 初のオープンソース実装: MMR 項と磁化方程式を完全に結合した MHD 微極流体用の OpenFOAM ソルバ（epotMicropolarFoam, epotMMRFoam）を初めて開発し、GitHub で公開した。
2. MMR の重要性の定量的証明: 血液のような低導電性流体において、ローレンツ力ではなく、MMR（磁気トルク）が流れ制御の主要なメカニズムであることを数値的に実証した。
3. 生物医学的応用への示唆: 強い磁場下では赤血球の回転が抑制され、流れの抵抗が増大し、再循環領域が安定化することを示した。これは、磁気ハイパーサーミアや標的薬物送達（Targeted Drug Delivery）などの医療応用において、磁場制御のメカニズム理解に不可欠な知見である。

5. 意義と将来展望 (Significance and Future Work)

• 科学的意義: 従来の MHD 研究では見落とされがちだった「微粒子の磁気的整列と回転抑制」が、巨視的な流れ場（速度分布、渦構造）に決定的な影響を与えることを明らかにした。
• 実用性: 開発されたソルバは、複雑な血管形状（動脈瘤、狭窄部など）における磁場制御血流の解析に直接適用可能であり、医療機器の設計や治療法の最適化に貢献する。
• 今後の課題:
  • 現在のシミュレーションは層流・定常流に限定されているため、拍動流（Pulsatile flow）への拡張が必要。
  • 血液の非ニュートン性や赤血球の変形性、凝集現象のモデル化。
  • 患者固有の解剖学的構造（Patient-specific geometries）への適用。

本論文は、磁気流体と微極流体理論の融合を計算流体力学の観点から飛躍させ、特に生体流体の磁気制御に関する新たな研究分野を開拓する重要な基盤を提供した。