[論文レビュー] Computational study of shock waves propagating through air-plastic-water interfaces

Mauricio J. del Razo, Randall J. LeVeque|arXiv (Cornell University)|Mar 31, 2015

Combustion and Detonation Processes参考文献 11被引用数 4

ひとこと要約

本研究では、爆風による脳震盪（TBI）を想定し、1次元有限体積法とTammann状態方程式、およびハイブリッドHLLC-正確なリーマン解法を用いて、空気-プラスチック-水界面における衝撃波の伝播をモデル化している。主な発見は、薄いプラスチック界面はシミュレーションにおいて無視可能であり、その影響は衝撃波の伝播にほとんど及たないため、TBI研究における複雑な生物学的界面のモデル化を簡略化できるということである。

ABSTRACT

The following study is motivated by experimental studies in traumatic brain injury (TBI). Recent research has demonstrated that low intensity non-impact blast wave exposure frequently leads to mild traumatic brain injury (mTBI); however, the mechanisms connecting the blast waves and the mTBI remain unclear. Collaborators at the Seattle VA Hospital are doing experiments to understand how blast waves can produce mTBI. In order to gain insight that is hard to obtain by experimental means, we have developed conservative finite volume methods for interface-shock wave interaction to simulate these experiments. A 1D model of their experimental setup has been implemented using Euler equations for compressible fluids. These equations are coupled with a Tammann equation of state (EOS) that allows us to model compressible gas along with almost incompressible fluids or elastic solids. A hybrid HLLC-exact Eulerian-Lagrangian Riemann solver for Tammann EOS with a jump in the parameters has been developed. The model has shown that if the plastic interface is very thin, it can be neglected. This result might be very helpful to model more complicated setups in higher dimensions.

研究の動機と目的

• 低強度の爆風衝撃波が軽度脳震盌傷（mTBI）を引き起こすメカニズムを理解すること。特に、従来の画像診断では損傷が検出されない場合のメカニズムを解明すること。
• 実験データだけではこれらの界面における複雑なダイナミクスを解明できないため、1次元モデルを用いて空気-プラスチック-水界面を通過する衝撃波の伝播をシミュレートすること。
• 実験装置（例：トランスウェル容器）におけるプラスチック界面が衝撃波伝播に顕著な影響を及ぼすかどうかを特定し、計算モデルにおいてその界面を単純化または省略可能かどうかを評価すること。
• 特に空気から水への界面を通過する際の圧力増幅と波形歪みを定量化すること。これは、脳内皮細胞への損傷リスクを評価する上で重要である。
• 爆風誘発性mTBIの文脈において、将来の3次元シミュレーションや二相流モデルの基盤を提供すること。

提案手法

• 圧縮流体のEuler方程式を解くために、保存則に基づく有限体積法を用い、圧縮性ガスおよびほぼ不圧縮性流体／固体をモデル化するためのTammann状態方程式（EOS）を併用している。
• Tammann EOSに特化した、材料パラメータの不連続およびジャンプを扱えるハイブリッドHLLC-正確なEuler-Lagrangianリーマン解法を開発した。
• 1次元モデルは、空気、薄いプラスチック層、水を含む実験的ショックチューブ設定を再現し、衝撃波が界面に作用する領域に焦点を当てている。
• 計算実験では、プラスチック界面の厚さを変化させ、その影響が伝播する衝撃波の振幅および形状に与える影響を評価した。
• 線形音響理論との比較により、プラスチックの厚さが小さくなるにつれて、結果が期待される解析解に収束することを確認し、モデルの妥当性を検証した。
• 全プレート（空気-プラスチック-水）ケースと比較して、簡略化された空気-水界面モデルの精度を評価した。

実験結果

リサーチクエスチョン

• RQ1 プラスチック界面の厚さが、空気-プラスチック-水界面における衝撃波伝播に与える影響は何か？
• RQ2 直接的な空気-水界面と比較して、プラスチック層が伝播する衝撃波の振幅および形状にどの程度の影響を及ぼすか？
• RQ3 特に水中における衝撃波増幅の文脈において、計算モデルからプラスチック界面を無視しても、精度に顕著な損失がないか？
• RQ4 空気から水に衝撃波が伝播する際の圧力増幅の大きさは何か？また、脳組織における損傷メカニズムとどのように関連するか？
• RQ5 薄いプラスチック層内での複数の反射および干渉パターンが、下流の時間分解圧力プロファイルにどのように影響するか？

主な発見

• ゲージ4（水中）での最大圧力振幅は、空気-水界面ケースで184.06 KPaの初期空気衝撃波から54％増加した284.26 KPaに達した。
• プラスチック界面の厚さを0.1 mに減少させた場合、水中の伝播する衝撃波波形は、空気-水界面ケースとほぼ区別できず、最大振幅の誤差は0.38％にとどまった。
• プラスチックの厚さが小さくなるにつれて、プラスチック層内での多重反射の干渉が極めて速くなり、圧力プロファイルは単一の衝撃波形状に収束する。
• 最も薄いプラスチックケース（0.1 m）とプラスチックなしケース（空気-水）の圧力ゲージプロットは視覚的に区別できず、プラスチック層の影響が無視できることが示された。
• シミュレーション結果は、線形音響理論の予測と一致しており、薄いプラスチック層を通る伝播圧力は、直接的な空気-水界面の期待値に収束することが確認された。
• モデリングの観点から、空気-プラスチック-水界面は、特にプラスチックが非常に薄い場合、空気-水界面に単純化しても顕著な精度の損失がないことが示された。