초기 레일리 모드가 급가속 낙하액 방울 변형·파열에 미치는 영향

초록

본 연구는 축대칭 레일리 진동 모드( n,0 )와 그 위상(0, π)을 미리 지정해 부피 보존 상태로 초기화한 뒤, 다양한 밀도·점도·웨버수 조건에서 급가속을 가했을 때 방울의 변형·파열 메커니즘을 직접수치 시뮬레이션(VOF)으로 조사한다. 결과는 초기 표면에 저장된 진동 에너지와 외부 작업이 어떻게 회복 가능한 진동 에너지, 질량 중심 이동, 점성 소산으로 분배되는가에 따라 파열 여부가 결정되며, 점도·밀도비가 모드 결합과 에너지 흡수 시간창을 조절한다는 것을 보여준다.

상세 분석

이 논문은 기존에 대부분 구형 상태를 가정해 온 낙하액 파열 연구와 달리, 실제 산업·자연 현장에서 흔히 관찰되는 비구형, 진동이 내재된 초기 상태를 체계적으로 탐구한다. 저자는 레일리 모드 전개를 이용해 (2,0), (3,0), (4,0) 축대칭 ‘존알’ 모드를 정의하고, 진폭 A = 0.3 D 로 설정한 뒤 부피 보존을 위해 n = 0 모드(기저 반경 보정)를 추가하는 수학적 절차를 제시한다. 이렇게 만든 초기 형태는 표면 에너지와 내부 운동 에너지(진동 에너지)의 비율을 정확히 제어할 수 있다.

시뮬레이션은 Basilisk 기반 VOF 방법으로, 두 상의 연속성·운동량 방정식과 CSF 방식의 표면장력을 풀어 급가속(초기 속도 V₀) 상황을 재현한다. 파라미터 공간는 물‑공기(ρ≈815, Oh_d≈1.3×10⁻³, We≈12), 고점도 물(Oh_d≈1.3×10⁻¹, We≈18), 물‑다른 액체(ρ≈10, Oh_d≈1.3×10⁻³, We≈10) 세 경우로 구성돼, 각각 밀도비·점도비가 파열 임계값에 미치는 영향을 비교한다.

핵심 발견은 다음과 같다. (1) 초기 모드와 위상이 외부 공기·유체 압력에 의해 유도되는 변형 시간스케일 τ_D와 레일리 진동 주기 2π/ω_n이 비슷할 때, 두 현상이 건설적(증폭) 혹은 파괴적(완화)으로 상호작용한다. 예를 들어, (2,0) 모드의 프로레이트 위상(φ = 0)에서는 공기 압력에 의해 추가적인 신축이 일어나 파열이 촉진되지만, 역위상(φ = π)에서는 초기 수축이 압력에 맞물려 방울을 일시적으로 안정시킨다. (2) 점도는 진동 감쇠와 모드 결합을 조절한다. 고점도 경우에는 진동 에너지가 빠르게 소산돼 건설적 상호작용이 억제되고, 파열 임계 We가 크게 상승한다. (3) 밀도비가 클수록 질량 중심 이동 에너지가 증가해 외부 작업이 더 많이 질량 이동에 전환되며, 이는 파열을 지연시키는 효과를 만든다.

또한, 에너지 흐름 분석을 통해 외부 작업이 회복 가능한 진동 에너지(ΔE_osc), 질량 중심 이동 에너지(ΔE_cm), 점성 소산(ΔE_visc)으로 어떻게 분배되는지를 정량화하였다. 파열이 일어나는 경우는 ΔE_osc/ΔE_total 비율이 높을 때이며, 반대로 ΔE_cm 혹은 ΔE_visc 비중이 클수록 방울은 파열 없이 변형 후 재구성한다. 이러한 관점은 기존의 We_cr 기반 임계값 모델을 확장해, 초기 형태·에너지 상태를 포함한 새로운 파열 예측 지표를 제시한다는 점에서 의의가 크다.