A Stabilized Mixed Hybrid Finite Element Framework for Large-Deformation Poromechanics with Torsion: Application to Transmural Perfusion and Vascular Hemodynamics

人間の心臓は、単に血液を押し出すだけのポンプではなく、鼓動のたびにねじれ、伸び、自らの形を変えるポンプである。左心室の厚い筋壁の内部では、複雑な流体力学の問題が展開されている。心筋は単なる固体の組織の塊ではなく、およそ15パーセントが微細な管を通じて移動する流体で満たされた、スポンジのような素材である。心臓が血液で満たされ、それを射出するためにねじれる際、この内部の流体は絞られ、押し出され、心臓自身の微細な血管を通る血流を助けることもあれば、妨げることもある圧力を作り出す。もし筋組織内の圧力が高くなりすぎると、これらの繊細な血管を押しつぶし、心臓組織から酸素を奪ってしまう可能性がある。心臓のねじれる動きが、この内部の流体の流れとどのように相互作用するかを理解することは極めて重要であるが、従来のコンピュータモデルはこれを正確にシミュレートすることに苦慮することが多く、組織がほとんど圧縮不可能な状態になると、混乱を招くエラーを生じさせたり、「ロック（固着）」したりすることが頻繁にある。

アイントホーデンド工科大学の研究者は、この挙動をシミュレートするための新しい手法を開発し、心壁の大規模なねじれや伸びを、一般的な罠に陥ることなく捉える数学的枠組みを構築した。心筋を流体を含んだスポンジとして扱い、固体組織の動きと流体の流れを個別に追跡する特殊な手法を用いることで、この研究は、心臓の自然なねじれ運動がいかに保護メカニズムとして機能しているかを明らかにしている。シミュレーションによれば、心臓が充満する際に自由にねじれることが許されると、内壁の繊維へのストレスが大幅に軽減され、血管の崩壊を脅かす内部圧力が低下することが示された。逆に、高血圧によって心壁が異常に厚くなった場合、流体が閉じ込められ、圧力が危険なレベルまで上昇するが、標準的なコンピュータモデルではこの現象を正しく予測できないことが多い。

この研究の核心は、心筋のような柔らかく湿った素材のモデリングにおける長年の課題を解決することにある。標準的なコンピュータシミュレーションでは、材料がほぼ非圧縮的（つまり、体積を小さく押しつぶすことができない状態）である場合、数学的な計算が破綻し、現実とは異なる不自然な圧力の急上昇や、ギザギザとした振動的な結果が生じることがある。これは「体積ロッキング（volumetric locking）」として知られている。これを解決するために、研究者は「混合ハイブリッド（mixed hybrid）」アプローチを考案した。すべてを一括で計算しようとするのではなく、この手法は問題をより小さな断片に分解し、各断片の境界を越えて流体が連続的に流れるという厳格なルールを適用することで、現実の体内で行われているのと同様の挙動を強制する。この技術により、コンピュータは組織の動きとシミュレーションブロックの端における圧力を別々に解くことができ、流体の保存が正確に行われ、極端な変形下においても結果が滑らかで物理的に現実的なものとなる。

この新しいフレームワークをテストするために、研究者はまず、膨張し、かつねじられた厚肉円筒に関する精密な数学的解を導出し、これを完璧なベンチマークとした。次に、この手法を左心室のモデル、すなわち人間の解剖学的構造に一致する寸法を持つ心筋の円筒部分に適用した。シミュレーションでは、3つの特定のシナリオを検証した。第一に、心臓自身の冠血管内の血圧レベルの違いが、壁の剛性にどのように影響するかを調査した。その結果、血管が血液で十分に満たされているとき、内部の流体圧力が上昇し、それが弛緩した心壁を硬化させることを発見した。「ホースの園芸（garden-hose effect）」として知られるこの現象は、心臓が充填されるためにより高い圧力を必要とし、取り込める血液量を制限することを意味する。新しいモデルは、この変化を正確に予測し、高い内部灌流圧が筋壁の固体応力を10キロパスカルから3キロパスカルへと減少させ、流体と組織の間で負荷を効果的に分担していることを示した。

第二のシナリオでは、心臓のねじれ運動の重要な役割を検証した。健康な心臓では、筋線維は内表面から外表面に向かって螺旋状に配置されており、進むにつれて角度が変化する。心臓が充満する際、この螺旋構造によって壁はねじれる。シミュレーションでは、自由にねじれることが許された心臓と、回転が固定された心臓を比較した。ねじれが阻止されると、内壁の筋線維にかかるストレスは57.20キロパスカルへと急上昇し、危険な力の集中を引き起こした。しかし、モデルが自然にねじれることが許されると、そのピークストレスは45.5パーセント減少して31.20キロパスカルとなった。この力の再分配は、壁を通じた流体の排出も容易にし、内部の孔隙圧を約128ミリメートル水銀柱から122ミリメートル水銀柱へと低下させた。この研究は、心臓のねじれが単なる収縮の副産物ではなく、ストレスを緩和し、筋層への血流を維持するための不可欠なメカニズムであることを示している。

第三の、そして最も重要なケーススタディでは、高いポンプ圧に耐えるために壁が著しく肥厚する、重度の高血圧状態にある心臓を取り上げた。このシナリオでは、心臓の内半径を2.00ミリメートル、外半径を5.00ミリメートルとし、腔内圧を180ミリメートル水銀柱としてモデル化した。このような条件下では、厚くなった壁によって、流体が脱出するための経路が長くなるため、内表面付近に流体が閉じ込められる。新しい混合ハイブリッドモデルは、この閉じ込められた流体が、初期の負荷の後でも持続する高い圧力帯を作り出すことを示した。対照的に、研究者が標準的なコンピュータモデリング技術を用いて同じシナリオを実行したところ、結果は混沌としたものとなった。標準的な手法は、内壁付近で非物理的な圧力の振動を生じさせ、真の応力分布を捉えることができず、12.75パーセントという膨大な誤差を生じさせた。しかし、新手法は、わずか0.01パーセントの誤差で滑らかな単調増加の圧力曲線を描き、厚い壁とねじれ運動が、排水を遅らせ、高い内部圧力を持続させる機械的な障害物となることを正しく特定した。

これらの知見は、心臓がねじれ、また元に戻る能力が、その水力学的な健康状態の根本的な部分であることを示唆している。心臓が自由にねじれるとき、ストレスは壁全体に均一化され、流体の排出が促進され、微細血管の崩壊から保護される。壁が厚くなったり、ねじれが制限されたりすると、流体が閉じ込められ、血流を損なう高圧状態を招く。本研究は、新しい計算フレームワークが、これらの微妙かつ生命に関わる相互作用を捉えることができることを裏付けており、機械的な力と流体の流れがどのように結合しているかを理解するための信頼できるツールを提供している。現在の研究は概念を証明するために簡略化された円筒形状に焦点を当てているが、大きな変形や複雑な流体・構造相互作用を扱うこの手法の成功は、こうしたモデルがより複雑な三次元の心臓形状に応用され、心不全や高血圧の影響をより深く理解できる未来を示唆している。