[論文レビュー] A simple artificial damping method for total Lagrangian smoothed particle hydrodynamics

Chi Zhang, Yujie Zhu+3 more|arXiv (Cornell University)|Feb 5, 2021

Fluid Dynamics Simulations and Interactions参考文献 28被引用数 7

ひとこと要約

本論文は、ボン・ニューマン＝リヒトマイヤーにインspiredしたスケーリング係数を有するケルビン＝ヴォイクト型ダンパーに基づく単純な人工ダンピング手法を提案し、全ラグランジュ型滑らか粒子法（TL-SPH）シミュレーションの安定化を図る。この手法は鋭い勾配領域における不自然な振動を効果的に抑制し、大変形を伴う固体動力学、特にバイオメカニクス分野における複雑なステント構造の安定的かつ高精度なシミュレーションを可能にする。

ABSTRACT

In this paper, we present a simple artificial damping method to enhance the robustness of total Lagrangian smoothed particle hydrodynamics (TL-SPH). Specifically, an artificial damping stress based on the Kelvin-Voigt type damper with a scaling factor imitating a von Neumann-Richtmyer type artificial viscosity is introduced in the constitutive equation to alleviate the spurious oscillation in the vicinity of the sharp spatial gradients. After validating the robustness and accuracy of the present method with a set of benchmark tests with very challenging cases, we demonstrate its potentials in the field of bio-mechanics by simulating the deformation of complex stent structures.

研究の動機と目的

• 鋭い空間的勾配近くで生じる不自然な振動によって引き起こされるTL-SPHにおける数値的不安定性を解消すること。
• 大変形を伴う困難な固体動力学問題に対するTL-SPHのロバスト性と精度を向上させること。
• SPHを用いて、心血管ステントを含む複雑なバイオメカニカル構造の信頼性のあるシミュレーションを可能にすること。
• 既存のTL-SPHフレームワークにスムーズに統合可能な、簡単で実装が容易な安定化技術を提供すること。
• ベンチマークテストおよび実際のバイオメカニカル応用を通じて、本手法の有効性を実証すること。

提案手法

• ケルビン＝ヴォイクトダンパーのモデルに基づき、構成方程式に人工ダンピング応力項を導入する。
• ボン・ニューマン＝リヒトマイヤーの人工粘性にインspiredしたスケーリング係数を用いて、ダンピングの大きさを制御する。
• 最初のピオラ＝クリーフ応力テンソルにケルビン＝ヴォイクトのダンピング寄与を加えることで、応力の時間発展を修正する。
• 元のTL-SPH定式化を保持しつつ、この単純な追加により安定性を向上させる。
• 変形勾配の計算における一次の整合性を維持するために、カーネル補正行列を用いる。
• 再現可能性およびアクセス性を高めるために、本手法をオープンソースのSPHinXsysライブラリに実装する。

実験結果

リサーチクエスチョン

• RQ1 鋭い空間的勾配近くで生じる不自然な振動は、どのようにしてTL-SPHシミュレーションで効果的に抑制できるか？
• RQ2 単純なケルビン＝ヴォイクトに基づく人工ダンピング手法は、精度を損なわせることなくTL-SPHを安定化できるか？
• RQ3 極端な変形および高いひずみ勾配を伴うベンチマークケースにおいて、本手法はどのように性能を発揮するか？
• RQ4 本手法は、ステントのような複雑な現実世界のバイオメカニカル構造を信頼性高くシミュレートできるか？
• RQ5 スケーリング係数は、さまざまな材料モデルにおいて数値的安定性および精度にどのように影響を与えるか？

主な発見

• 提案されたTL-SPH-KV手法は、ソードチューブおよびショックチューブテストを通じて、鋭い空間的勾配領域における物理的でない振動を完全に除去したことが実証された。
• 極端な変形下でも、弾性リングおよびプレート衝突テストを含むベンチマーク問題において、高い精度を維持した。
• PS形状ステントのケースでは、ひし形変形に起因するセルのコーナー部にボン・ミーゼス応力集中が正確に捉えられ、現実的な応力分布が再現された。
• C形状ステントのケースでは、曲がったストラット部に応力が局在化し、自由端の拡張が正しく再現された。これは、反対称的制約下でも本手法のロバスト性を示している。
• 本手法により、現実的なネオ＝フック型心血管ステントのSPHベースのシミュレーションが初めて成功裏に実現された。これにより、バイオメカニカルモデリングの新たな道が開かれた。
• すべてのシミュレーション結果は、期待される物理的挙動と整合しており、本手法は多様な材料および幾何的配置において一貫した性能を示した。