[논문 리뷰] Further investigation of the spurious interface fragmentation in multiphase Smoothed Particle Hydrodynamics

Kamil Szewc, Michał T. Lewandowski|arXiv (Cornell University)|2016. 02. 25.

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

한 줄 요약

이 논문은 다상 유동을 위한 스플라인 입자 방법(SPH)에서 널리 사용되는 인공적 표면계면 보정 항목의 사용을 도전하며, 표면장력이 무시될 경우 발생하는 '오류 있는 계면 분열' 현상이 실제로는 켈빈-헬름홀츠 불안정성의 물리적으로 타당한 표현임을 주장한다. 안정성 분석과 SPH/VOF 시뮬레이션을 통해 이러한 보정 항목이 비물리적 해를 유도할 수 있음을 입증하고, 이전 연구에서 관측된 보정 계수 ε와 스무딩 길이 h 사이의 역관계가 유체역학적 안정성 이론의 직접적 결과임을 설명한다.

ABSTRACT

This article presents results of further investigation of the problem of spurious interface fragmentation in the multiphase SPH. In order to remove arising instabilities, many authors introduced the artificial interface correction procedure. In the present paper we show that the interface instabilities are physical and the introduction of the interface correction procedure may leads to non-physical solutions. We also explain the puzzling relation between the parameters $\varepsilon$ and $h$. The analysis is performed on the basis of the stability analysis and numerical calculations using SPH and Volume Of Fluid (VOF) approach.

연구 동기 및 목표

• 표면장력이 무시될 경우 다상 SPH에서 발생하는 계면 불안정성의 물리적 기원을 조사하기 위해.
• 일반적으로 시뮬레이션 안정화를 위해 사용되는 인공적 계면 보정 항목의 타당성을 도전하기 위해.
• 이전 연구에서 관측된 경험적 관계 ε_min ∝ 1/h를 명확히 하기 위해.
• 서스만의 버블 시험 케이스에서 표면장력이 비록 작지만, 핵심적인 안정화 역할을 한다는 것을 입증하기 위해.
• SPH 시뮬레이션에서 계면 보정 절차의 사용에 있어 신중함을 권장하기 위해.

제안 방법

• 균일한 속도와 밀도 대비를 갖는 이중 유체 계면에 대해 선형화된 켈빈-헬름홀츠 안정성 분석을 수행하였다.
• 계면 보정을 모의하기 위해 인공적 법선 응력 항목을 포함한 분산 관계를 유도하였다 (ε-항).
• 이전 기준과의 일관성을 확보하기 위해 비차원 단위(반지름 R = 1 cm, 중력 가속도 g = 1 cm/s², 밀도 ρ = 1 g/cm³)를 사용하였다.
• 계면 보정 항목이 포함된 경우와 포함되지 않은 경우의 SPH 시뮬레이션을 수행하고, VOF 방법 결과와 비교하였다.
• 수치 실험을 통해 최소 안정 ε(ε_min)가 스무딩 길이 h에 따라 어떻게 변화하는지 분석하였다.
• 안정 조건에서 유도된 이론적 관계 ε_min ∝ 1/h를 통해 수치적 관측 결과를 확인하였다.

실험 결과

연구 질문

• RQ1 표면장력이 무시될 경우 SPH 시뮬레이션에서 관측되는 계면 분열 현상은 수치적 오류인지, 물리적 불안정성인지?
• RQ2 인공적 계면 보정 항목(ε-항)이 SPH 계면 안정화에 수행하는 진정한 물리적 역할은 무엇인가?
• RQ3 왜 최소 안정 ε 값이 경험적으로 관측된 것처럼 스무딩 길이 h와 반비례하는가?
• RQ4 서스만 시험 케이스에서 비록 작지만 표면장력은 상당한 안정화 역할을 하는가?
• RQ5 관측된 ε_min ∝ 1/h 관계는 수압 안정성 이론에서 분석적으로 유도될 수 있는가?

주요 결과

• 표면장력이 무시된 SPH 시뮬레이션에서 관측된 계면 불안정성은 물리적으로 타당하며, 켈빈-헬름홀츠 불안정성의 시작을 나타낸다.
• 인공적 계면 보정 항목(ε-항)은 물리적 결함을 수정하는 것이 아니라 실제 불안정성을 가리기 때문에, 비물리적 해를 유도할 수 있다.
• 이전 연구에서 경험적으로 관측된 scaling ε_min ∝ 1/h는 켈빈-헬름홀츠 분석에서 도출된 안정 조건에 의해 분석적으로 설명된다.
• 서스만 시험 케이스에서 비록 작지만 표면장력은 계면을 안정화시키고 분열을 방지하는 데 충분하다.
• 계면 보정 항목은 주로 수치적 '미적 보정' 역할를 하며, 물리적 안정화 요소는 아니며, 사용 시 신중함이 필요하다.