[논문 리뷰] A modified SPH approach for fluids with large density differences

F. Ott, Erik Schnetter|arXiv (Cornell University)|2003. 03. 27.

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

한 줄 요약

이 논문은 다상유동에서 큰 밀도 차이가 있는 경우 시뮬레이션을 안정화하기 위해 질량 밀도 대신 입자 밀도를 이산화하는 수정된 스무딩 입자 유체역학(SPH) 방법을 제안한다. 입자 밀도를 사용하는 밀도 재구성 방식으로 전환함으로써, 비정상적인 압력 진동을 제거하고 충격파관, 음파 전파, 이송 문제에서 정확도를 향상시킨다. 특히 입자의 질량이나 공간 해상도가 크게 다를 경우에도 효과적이다.

ABSTRACT

We introduce a modified SPH approach that is based on discretising the particle density instead of the mass density. This approach makes it possible to use SPH particles with very different masses to simulate multi-phase flows with large differences in mass density between the phases. We test our formulation with a simple advection problem, with sound waves encountering a density discontinuity, and with shock tubes containing a contact discontinuity between air and Diesel oil. For all examined problems where particles have different masses, the new formulation yields better results than standard SPH. This is also the case for problems in which different spatial resolutions are needed while the mass density does not change.

연구 동기 및 목표

• 유체 경계에서 큰 밀도 차이가 있는 다상유동을 시뮬레이션할 때 표준 SPH가 보이는 불안정성과 부정확성을 해결한다.
• 표준 SPH에서 입자 간 질량 차이가 커지면 불연속성 근처에서 밀도 및 압력 계산이 부정확해지는 근본적인 문제를 해결한다.
• 입자의 질량과 공간 해상도가 다를 경우에도 안정적인 시뮬레이션을 가능하게 하며, 단일 상에서 질량 밀도가 일정한 경우에도 적용 가능하다.
• 경직된 상태방정식(예: 액체-기체 또는 액체-액체 경계)을 가진 압축성 유동에서 안정성과 정확도를 유지하는 형식을 개발한다.
• 접촉 불연속성, 충격파, 음파를 포함한 높은 밀도 대비 문제에 대해 표준 SPH의 수치적 안정성 있는 대안을 제공한다.

제안 방법

• 표준 SPH의 질량 밀도 기반 밀도 재구성 방식을 입자 밀도 기반으로 대체한다. 입자 밀도는 입자 부피의 역수로 정의된다.
• SPH 커널 합산을 수정하여 질량 밀도 대신 입자 밀도를 계산하도록 재해석한다. 수식은 $ \tilde{\rho}_{\text{part}}(\mathbf{x}_i) = \sum_j V_j W_{ij} $ 로, 여기서 $ V_j $ 는 입자 부피이다.
• 입자 질량 $ m_i $ 와 입자 밀도 $ \rho_{\text{part},i} $ 를 사용하여 실제 유체 질량 밀도를 $ \rho_i = m_i / V_i $ 로 계산함으로써 입자 운동과의 일致성을 확보한다.
• 운동 방정식(모멘텀 및 에너지 진화 포함)에 수정된 형식을 적용하면서도 메서드의 라그랑주 성격을 유지한다.
• 충격파를 안정화하기 위해 인공 점성 계수 $ \alpha = 0.5 $ 를 사용하며, 약 8개의 스무딩 길이에 걸쳐 분포시키되, 초기 불연속성을 악화시키지 않는다.
• 필요한 경우에만 초기 자료에 스무딩을 적용하지만, 불연속적인 초기 조건이 새로운 형식을 불안정하게 하지 않음을 보여준다.

실험 결과

연구 질문

• RQ1 질량 밀도 대신 입자 밀도 기반의 수정된 SPH 형식이 큰 밀도 대비를 보이는 다상유동에서 안정성과 정확도를 향상시키는가?
• RQ2 표준 SPH와 비교해 볼 때, 새로운 형식은 밀도 불연속성의 이송을 어떻게 시뮬레이션하는가?
• RQ3 수정된 SPH 방법은 공기와 디젤 기름(밀도 비율 약 1:80)을 포함한 고밀도 대비 충격파관에서 충격파 및 팽창파를 정확하게 포착하는가?
• RQ4 다른 질량과 공간 해상도를 가진 입자를 사용할 때도 비정상적인 압력 진동 없이 처리할 수 있는가?
• RQ5 초기 자료에 불연속성이나 고주파 모드가 포함되어 있어도 새로운 형식은 견고한가?

주요 결과

• 수정된 SPH 형식은 밀도 불연속성이 있는 1차원 이송 문제를 성공적으로 시뮬레이션하여 비정상적인 압력 진동 없이 정확한 경계 유지가 가능했다.
• 음파 시험에서는 새로운 방법이 전파 및 반사 계수를 정확히 포착하는 반면, 표준 SPH는 잘못된 파동 행동과 압력 피크를 발생시켰다.
• 공기와 디젤 기름(밀도 비율 약 1:80)을 포함한 충격파관에서 접촉 불연속성이 안정적으로 유지되고 잘 해상화되었으며, 시간이 지남에 따라 압력 진동이 증가하지 않았다.
• 충격파 속도($ v_s = 1056 \pm 20 \, \text{m/s} $), 충격 후 속도($ v_D = 5.94 \pm 0.02 \, \text{m/s} $), 압력 상승($ \Delta p_D = 4.94 \pm 0.02 \, \text{MPa} $)의 시뮬레이션 결과가 해석적 해와 매우 유사하게 일치했다.
• 모든 주요 충격 파rameter에 대해 해석적 충격 관계와 뛰어난 일치를 보였으며, 상대 오차는 2% 이하로 매우 높은 정확도를 입증했다.
• 표준 SPH는 경직된 유체에서 실제 압력보다 수개의 주기수를 초과하는 비물리적 압력 오차로 인해 동일 조건에서 안정적인 시뮬레이션을 수행하지 못했다.