[논문 리뷰] A simple artificial damping method for total Lagrangian smoothed particle hydrodynamics

Chi Zhang, Yujie Zhu외 3명|arXiv (Cornell University)|2021. 02. 05.

Fluid Dynamics Simulations and Interactions참고 문헌 28인용 수 7

한 줄 요약

이 논문은 총 라그랑주 스플린트 화학역학(SPH, Total Lagrangian Smoothed Particle Hydrodynamics) 시뮬레이션의 안정성을 향상시키기 위해 콜드버그-보이트 타입 댐퍼에 바탕을 두고 비온-뉴먼-리히트마이어가 영감을 줬던 척도 인자를 도입한 단순한 인공 감쇠 방법을 제안한다. 이 방법은 날카로운 기울기 근처에서 발생하는 비물리적 진동을 효과적으로 억제하여, 복잡한 stent 기하구조를 포함한 대변형 고체 역학의 견고하고 정확한 시뮬레이션을 가능하게 한다.

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를 사용하여 심장혈관 stent과 같은 복잡한 생체역학적 구조를 신뢰할 수 있게 시뮬레이션할 수 있도록 한다.
• 기존의 TL-SPH 프레임워크에 부드럽게 통합될 수 있는 단순하고 쉽게 구현할 수 있는 안정화 기법을 제공한다.
• 기준 시험과 실제 생체역학적 적용 사례를 통해 이 방법의 효과성을 입증한다.

제안 방법

• 콜드버그-보이트 댐퍼 모델에 기반한 구성 방정식에 인공 감쇠 응력 항을 도입한다.
• 비온-뉴먼-리히트마이어의 인공 점성에서 영감을 얻은 척도 인자를 적용하여 감쇠 크기를 제어한다.
• 제1 피아올라-키르히호프 응력 텐서에 콜드버그-보이트 감쇠 기여를 추가하여 응력 진화를 수정한다.
• 기존의 TL-SPH 수식을 유지하면서도 이 간단한 추가를 통해 안정성을 향상시킨다.
• 변형률 계산에서 일阶 일致성을 유지하기 위해 커널 보정 행렬을 사용한다.
• 재현 가능성과 접근성을 확보하기 위해 이 방법을 오픈소스 라이브러리 SPHinXsys에 구현한다.

실험 결과

연구 질문

• RQ1 TL-SPH 시뮬레이션에서 날카로운 공간 기울기 근처의 비물리적 진동을 효과적으로 억제할 수 있는가?
• RQ2 간단한 콜드버그-보이트 기반 인공 감쇠 방법이 정확도를 훼손하지 않고 TL-SPH를 안정화시킬 수 있는가?
• RQ3 극한의 변형과 높은 응력 기울기 조건을 가진 기준 시험 케이스에서 이 방법의 성능은 어떠한가?
• RQ4 이 방법은 심장혈관 stent과 같은 복잡한 실제 생체역학적 구조를 신뢰성 있게 시뮬레이션할 수 있는가?
• RQ5 척도 인자가 다양한 재료 모델에서 수치적 안정성과 정확도에 미치는 영향은 어떠한가?

주요 결과

• 제안된 TL-SPH-KV 방법은 Sod 튜브 및 충격파 튜브 시험을 통해 날카로운 공간 기울기 영역에서 비물리적 진동을 성공적으로 제거하였다.
• 극한의 변형 조건에서도 탄성 링 및 판 충격 시험과 같은 기준 문제에서 높은 정확도를 유지하였다.
• PS형 stent의 경우, 마름모꼴 변형으로 인해 세포 모서리에서 베르누이 응력 집중이 정확하게 포착되었으며, 현실적인 응력 분포를 확인하였다.
• C형 stent의 경우, 곡선 끝부분과 자유단 팽창 영역에서 높은 응력 집중이 정확하게 시뮬레이션되었으며, 이는 반대칭 경계 조건 하에서도 견고함을 보여주었다.
• 이 방법을 통해 실제 neo-Hookean 심장혈관 stent의 첫 번째 성공적인 SPH 기반 시뮬레이션이 가능해졌으며, 생체역학 모델링의 새로운 길을 열었다.
• 모든 시뮬레이션 결과는 예상되는 물리적 거동과 일치했으며, 다양한 재료 및 기하구조 조건에서 일관된 성능을 보였다.