스무스 입자 유체역학의 혼합 문제 해결

초록

본 논문은 기존 스무스 입자 유체역학(SPH)에서 발생하는 경계면 혼합 실패 원인을 두 가지로 규명한다. 첫째는 이웃 입자 수가 늘어나도 사라지지 않는 운동량 방정식의 1차 오차이며, 이는 커클링과 밴딩 불안정성에 의해 커널이 비정규적으로 샘플링되기 때문이다. 둘째는 엔트로피 보존으로 인한 국부 혼합 불안정성(LMI)으로, 압력 불연속을 초래한다. 저자는 적절한 커널 선택과 가중 밀도 추정법을 도입한 최적화 SPH(OSPH) 방식을 제시하고, 켈빈-헬름홀츠 불안정성 테스트와 고밀도 구체 ‘블롭 테스트’에서 기존 Eulerian 코드와 동등한 성능을 보임을 입증한다.

상세 분석

스무스 입자 유체역학(SPH)은 입자 기반 라그랑지안 형식으로 연속 방정식과 운동량 방정식을 커널 함수를 통해 근사한다. 그러나 기존 SPH는 서로 다른 엔트로피를 가진 유체가 접촉하는 경계면에서 혼합이 억제되는 현상을 보이며, 이는 수치적 불안정성과 물리적 모델링 한계가 복합적으로 작용한다는 점이 본 연구의 핵심 분석이다. 첫 번째 문제는 운동량 방정식에 존재하는 𝑂(𝑁⁻¹) 수준의 편차이다. 이 편차는 이웃 입자 수(N) 가 증가하면 이론적으로 사라져야 하지만, 실제 시뮬레이션에서는 커널이 불규칙하게 샘플링되면서 유지된다. 저자는 이를 ‘클러밍 불안정성(clumping instability)’과 ‘밴딩 불안정성(banding instability)’이라는 두 가지 메커니즘으로 구분한다. 클러밍은 입자들이 서로 뭉쳐 커널 중심에 과도한 밀도가 집중되는 현상이며, 밴딩은 입자 배열이 주기적인 파동 형태로 변형돼 커널의 대칭성이 깨지는 현상이다. 두 불안정 모두 커널 함수의 형태와 스무딩 길이 선택에 민감하게 반응한다. 특히, 전통적인 골든 스플라인 커널은 고밀도 대비 저밀도 영역에서 샘플링이 불균형해지면서 위 두 현상이 증폭된다. 저자는 고차원 다중 피크 커널(예: 웨이브렛 기반 커널)과 적절한 스무딩 길이 조정을 통해 입자 간 거리를 균등하게 유지하고, 불안정성을 억제함으로써 이론적인 𝑂(𝑁⁻¹) 수렴을 회복한다. 두 번째 문제인 ‘국부 혼합 불안정성(LMI)’은 입자들이 커널 스케일에서 혼합하려 할 때 엔트로피 보존 조건이 압력의 단일값성을 깨뜨리는 현상이다. 엔트로피가 서로 다른 두 입자가 같은 부피를 공유하면, 동일한 압력 대신 엔트로피에 비례한 압력 차이가 발생해 인위적인 압력 장벽이 형성된다. 이 장벽은 물리적으로는 존재하지 않으며, 경계면에서 유체가 서로 밀어내는 효과를 만든다. 저자는 이 문제를 해결하기 위해 ‘가중 밀도 추정(weighted density estimate)’을 도입한다. 기존 SPH는 각 입자의 질량을 커널 가중 평균으로 나누어 밀도를 계산하지만, 가중 밀도는 각 입자의 엔트로피 가중치를 포함해 압력 연속성을 보장한다. 구체적으로, ρ_i = Σ_j m_j W(r_ij, h) (P_i + P_j)/(2P_ref) 형태의 가중식을 사용해 압력이 단일값이 되도록 설계한다. 이 방식은 또한 입자 부피 추정의 정확성을 향상시켜 연속 방정식과 운동량 방정식의 오차를 동시에 감소시킨다. 최종적으로 저자는 이러한 커널 선택과 가중 밀도 추정을 결합한 ‘Optimised SPH(OSPH)’ 알고리즘을 구현하고, 켈빈-헬름홀츠(KHI)와 블롭 테스트에서 기존 SPH가 보였던 인공적인 표면 장력과 파동 억제를 극복함을 입증한다. 결과는 OSPH가 Eulerian 고정격자 코드와 거의 동일한 성장률과 혼합 효율을 보이며, 특히 1:10 밀도 대비 상황에서도 정확한 충격파와 전단 흐름을 재현한다는 점에서 큰 의미를 가진다.