血栓塞栓症の発症メカニズムに関する数値シミュレーション

2003 Fiscal Year Annual Research Report

• Project/Area Number: 14750107
• Research Institution: Tohoku University
• Principal Investigator: 伊澤 精一郎 東北大学, 大学院・工学研究科, 助手 (90333856)
• Keywords: SPH法 / 移動境界 / 粘性モデル / 密度均一化アルゴリズム / 自由表面

Research Abstract

研究の初年度は、血流の数値シミュレーションを始めるにあたって、本研究で採用した粒子法の1種であるSPH(Smoothed Particle Hydrodynamics)法について、流入・流出の境界条件や粘性の取り扱い、それに壁面の構成方法など基礎的な事柄に関して検証を進めた。SPH法はもともと圧縮性流体の運動を解析する方法として出発しており、本研究では、非圧縮性流体に適用するために、各タイムステップに粒子数密度を一様化させる独自のアルゴリズムを開発した。このアルゴリズムに加え、粒子のカーネル関数の拡散方程式に基づいた粒子間の運動量交換から粘性モデルを定義し、その評価をストークス層、平板間流れおよび円管流れに適用することで行った。この結果、解析解との比較を通して最適なモデル係数を同定することで、流れ場を正しく表現できることが明らかとなった。また、密度の評価モデルに改良を加えることで、物体後流などの非定常な渦運動を含む流れ場の様子をシミュレートすることに成功した。
上記の基礎的事項の検証に一応のめどがついた段階で、自由表面を含む流れの計算をスタートさせた。これは、血流に含まれる赤血球や白血球などの物質を模擬するために必要なモデルを開発するためのものである。本研究では、上述の密度均一化アルゴリズムに加え、界面を構成している粒子に作用する表面張力を表面粒子位置における曲率半径から算出する自由表面モデルを新たに開発した。これらのモデルを片面だけが回転するディスク間に入った空気と油の気液2相流に対する計算に適用した。その結果、液体表面の粒子スケールの小さな凹凸を自己補正するアルゴリズムを併用することで、SPH法によって移動境界と自由表面を含む非圧縮流れの計算が可能であることが明らかとなった。

Research Products

• [Publications] 伊澤 精一郎, 岡地 範明, 熊 鰲魁, 福西 祐: "SPH法を応用した渦が存在する非圧縮流れ場の計算"日本機械学会流工部門講演論文集. (CD-ROM). (2003)
• [Publications] 黒田 将史, 伊澤 精一郎, 熊 鰲魁, 福西 祐: "片面が回転するディスク間の気液2相流の計算"日本機械学会流工部門講演論文集. (CD-ROM). (2003)