[論文レビュー] A one-dimensional model for aspiration therapy in blood vessels

Michaël Herty, Niklas Kolbe+1 more|arXiv (Cornell University)|Jan 1, 2024

Gastric Cancer Management and Outcomes被引用数 2

ひとこと要約

本稿では、血管とカテーテルの間の界面で緩和型リーマンスolverを用いた1次元双曲型モデルを提示し、脳血管における吸引血栓吸引術の流れをモデル化する。このモデルは、吸引による流れの逆転、圧力低下、血管の狭窄を捉え、より大きなカテーテルと強い吸引が血行動態的乱れを増幅することを示しており、数値結果では臨床的吸引条件における閉塞状況下で圧力波の伝播と血流の減少が確認されている。

ABSTRACT

Aspiration thrombectomy is a treatment option for ischemic stroke due to occlusions in large vessels. During the therapy a device is inserted into the vessel and suction is applied. A new one-dimensional model is introduced that is capable of simulating this procedure while accounting for the fluid-structure interactions in blood flow. To solve the coupling problem at the tip of the device a problem-suited Riemann solver is constructed based on relaxation of the hyperbolic model. Numerical experiments investigating the role of the catheter size and the suction forces are presented.

研究の動機と目的

• 血管内の流体-構造相互作用を捉える計算効率の良い1次元モデルの開発を目的とする。
• 血管とカテーテルの流れが衝突するカテーテル先端における複雑な結合問題に対処することを目的とする。
• カテーテルサイズと吸引強度が流れの逆転、圧力低下、血管の狭窄といった血行動態的変化に与える影響を調査することを目的とする。
• 双曲型偏微分方程式系における界面結合に、緩和型リーマンスolverを用いた数値的に安定なフレームワークを提供することを目的とする。
• 臨床応用の前に、シミュレーションによる血管内治療戦略の検証を可能にすることを目的とする。

提案手法

• 血管内での血液の1次元双曲型系を、質量保存則と運動量保存則に従い、血管弾性に基づく圧力則を用いて定式化する。カテーテルは剛体と仮定する。
• カテーテル内流れを粘性のないバーガース方程式でモデル化し、血管内流れの式とは独立に扱う。
• カテーテル-血管界面での結合を処理するため、補助変数と緩和率εを用いた緩和系を導入する。
• 界面における断面積と質量フラックスの連続性を保証するように、問題特有の緩和型リーマンスolverを構築する。
• 安定性と精度を確保するため、適応的CFL条件と特徴量に基づく波速推定を用いた有限体積法を適用する。
• カテーテル半径と吸引速度を変化させた数値実験を実施し、血行動態的反応を分析する。

実験結果

リサーチクエスチョン

• RQ1 吸引カテーテルの挿入が1次元血管モデルにおける圧力、血流速度、血管断面積に与える影響は何か？
• RQ2 カテーテルサイズが吸引時の圧力波の伝播および流れの逆転に与える影響は何か？
• RQ3 吸引力の変動がカテーテル先端における圧力低下と血管の狭窄に与える影響は何か？
• RQ4 緩和型リーマンスolverは、カテーテル付加部と自由流れ部の双曲型系を正確かつ安定に結合できるか？
• RQ5 この1次元モデルは、閉塞状況下で流れの逆転や波の反射といった臨床的関連血行動態現象をどの程度再現できるか？

主な発見

• カテーテルの挿入により、圧力では前方に伝播するショック波、血流では後方へ伝播する希薄波が発生し、血管断面積の減少と血流速度の低下を引き起こす。
• カテーテル半径を0.25 cm（0.1 cmと比較）に増大させると、圧力低下と血流減少の効果が顕著に増幅され、血管断面積の狭窄も顕著に強まる。
• 強い吸引（w = -10,000 cm/s）により流れの逆転が発生し、カテーテル先端で負のネット流量に起因するkink型の血流速度最小値が観測される。
• 遠位閉塞状況下でも波の反射と持続的高圧状態をモデルが的確に捉えており、圧力は最大1×10⁵ dyne/cm²に達する。
• w = -1,000 cm/sでの吸引は、閉塞血管における高圧状態の持続時間を短縮するため、血行動態的ストレスの軽減という治療的利点があると示唆される。
• 緩和型リーマンスolverは、強い吸引や大きな断面積の不連続性が生じる状況下でも、界面結合における安定性と収束性を確保している。