단일 SPH 입자 낙하의 물리적 타당성 분석

초록

본 논문은 다상 유동 및 자유표면 흐름을 SPH로 시뮬레이션할 때 발생하는 단일 입자(액체 방울·고체)의 움직임이 물리적으로 올바른지, 그 유효 직경은 어떻게 정의되는지, 그리고 해석 파라미터(h, h/Δr, 커널 형태)가 결과에 미치는 영향을 실험적으로 검증한다. 낙하 실험을 통해 드래그 계수와 속도를 비교한 결과, 입자 직경을 스무딩 길이 h로 정의할 때 가장 일치함을 확인하였다.

상세 분석

이 연구는 SPH(Smoothed Particle Hydrodynamics) 방식이 다상 흐름을 다룰 때 종종 발생하는 ‘단일 입자 현상’—즉, 한 상의 입자가 다른 상 속으로 고립되어 이동하는 상황—의 물리적 의미를 체계적으로 탐구한다. 저자들은 먼저 기존 다상 SPH 공식 중 Hu‑Adams(2006) 방식을 변형하여 연속 방정식을 입자 간 상대 변위 기반으로 재정의하고, 이를 WCSPH(Weakly Compressible SPH)와 결합했다. 이 새로운 연속식은 입자 간 거리 변화에만 밀도가 반응하도록 설계돼, 자유표면 근처에서 입자 수가 감소해도 밀도 계산이 크게 왜곡되지 않는다.

검증 사례로는 Rayleigh‑Taylor 불안정과 다상 댐‑브레이크 실험을 수행했으며, 특히 댐‑브레이크 시뮬레이션에서 두 종류의 단일 입자(밀도 차이에 의해 생성된 고체 입자와 자유표면 파편)가 나타나는 현상을 관찰했다. 여기서 핵심 질문은 “이 입자들의 이동이 실제 물리 현상을 재현하는가?”와 “입자의 유효 직경은 무엇인가?”이다.

이를 답하기 위해 저자들은 고정된 부피의 단일 입자를 중력에 의해 낙하시키는 간단한 실험을 설계했다. 입자 직경을 정의하는 세 가지 가설을 검증하였다: (1) 입자 간 초기 간격 Δr, (2) 스무딩 길이 h, (3) 질량과 밀도로부터 계산된 부피 기반 직경 DΩ. 각각에 대해 다양한 밀도비(%S/%L)와 해상도(h/L=1/16, 1/32)에서 시뮬레이션을 수행하고, 얻어진 속도와 드래그 계수를 이론적 Stokes식(저레인즈수) 및 Schiller‑Naumann 실험식과 비교했다.

결과는 다음과 같다. Δr을 직경으로 사용할 경우, 시뮬레이션 속도가 이론값보다 크게 편차를 보이며, 특히 저속 영역에서 표준편차가 크게 나타났다. DΩ를 사용하면 약간 개선되지만 여전히 오차가 존재한다. 반면 h를 직경으로 가정하면, 시뮬레이션 결과가 이론 곡선에 가장 근접하고, 드래그 계수 역시 실험 데이터와 일치하는 경향을 보였다. 이는 SPH 입자의 물리적 ‘크기’가 실제 입자 간 거리보다 스무딩 길이와 더 연관이 깊다는 것을 시사한다.

또한, 입자 속도가 빨라질수록(즉, 높은 레이놀즈수) 드래그 계수의 표준편차가 감소하고, 결과가 보다 안정적으로 수렴한다는 점을 확인했다. 이는 입자와 주변 유체 입자 간 상호작용이 충분히 많은 경우(해상도·h 비율이 충분히 작을 때) SPH가 연속 매체의 저항을 적절히 포착한다는 의미이다. 반대로, 저속·저해상도 상황에서는 입자와 주변 입자 간 간격이 커서 ‘입자‑입자 충돌’에 의한 인공적인 진동이 발생, 물리적 의미가 약해진다.

이러한 분석을 통해 저자들은 단일 입자 현상이 반드시 수치적 오류라기보다, 해상도와 커널 파라미터에 따라 물리적 의미를 부여할 수 있음을 강조한다. 특히, 입자 직경을 h로 정의하고, h/Δr 비율을 2.5 이상 유지하며, 부드러운 Wendland 커널을 사용하면 단일 입자 낙하 실험에서 신뢰할 만한 drag coefficient와 terminal velocity를 얻을 수 있다.

요약하면, SPH에서 단일 입자(방울·고체)의 물리적 움직임은 해석 파라미터에 크게 좌우되며, 적절한 스케일링(특히 h 기반 직경 정의)과 충분한 입자 밀도 확보가 필수적이다. 이 결과는 다상 흐름에서 미세 입자(예: 에멀전, 기포) 모델링 시 SPH 파라미터 선택에 실질적인 가이드라인을 제공한다.