파라메트릭 표면파가 유도하는 혼합 메커니즘

초록

수직 진동으로 발생한 파라메트릭(파라볼릭) 파동, 즉 파라데이 파가 밀도 차이가 작은 두 가용성 유체 사이의 계면을 교란시켜 혼합을 촉진한다. 실험과 직접수치시뮬레이션(DNS)을 통해 파동이 계면을 진동·전단시키고, 큰 파동이 붕괴하면서 기포와 가벼운 유체가 하부로 침투하는 과정을 확인하였다. 혼합 속도는 초기에는 계면 장벽이 약해지면서 증가하지만, 표면으로부터 전달되는 에너지가 감소하면서 다시 감소하는 비단조적 거동을 보인다. 저자들은 깊이에 따라 변하고 내부 리차드슨 수에 의존하는 난류 확산계수를 포함한 1차원 모델을 제시해 관측된 현상을 정성적으로 재현하였다.

상세 분석

본 연구는 파라데이(Faraday) 파가 가용성 두 층 사이의 미세한 밀도 구배를 어떻게 교란하고, 장기적인 혼합을 유도하는지를 체계적으로 탐구한다. 실험에서는 20 rad s⁻¹의 고정 진동수와 가변 진폭(a ω²/g)으로 수직 진동을 가해, 파라데이 불안정에 의해 서브하모닉(ω/2) 파동을 유도하였다. 파동의 파장은 탱크 폭에 맞춰 선택되었으며, 얕은 경우(초기 계면 깊이 h_init≈2 cm)와 깊은 경우(h_init≈27 cm)를 포함한 다양한 초기 구배를 시험하였다. 영상 처리와 임계값 기반 세분화 기법을 통해 자유표면(ζ_surf)과 내부 계면(ζ_mix)의 위치를 실시간으로 추적했으며, 기포 발생·파고 붕괴 현상을 시각화하였다.

수치적으로는 다상 흐름을 부피분율(f) 변수와 연계한 비압축성 Navier‑Stokes 방정식(∇·u=0, ρ(∂u/∂t+u·∇u)=−∇p+∇·τ+g(t))을 3차원 DNS로 풀었다. 표면 장력과 점성, 그리고 분자 확산(D_mix≈1.3×10⁻⁹ m² s⁻¹)을 포함해 실제 실험 조건을 그대로 재현했으며, 특히 가속도 항 a ω² cos(ωt) 를 중력에 가산함으로써 파라데이 파의 강도를 조절하였다. DNS 결과는 실험과 일치하게, 작은 진폭에서는 계면이 단순히 파동에 따라 진동하고, 임계 진폭을 초과하면 파고가 붕괴하면서 기포와 가벼운 물이 하부로 급격히 투입되는 현상을 보여준다.

혼합 속도는 두 상호작용 메커니즘에 의해 비단조적이다. 첫 단계에서는 계면의 밀도 차이(Atwood 수)가 감소함에 따라 에너지 장벽이 낮아져 난류가 쉽게 발생하고, 혼합층 두께가 급격히 증가한다. 그러나 파동 에너지는 자유표면에서 깊이로 전달될수록 급감하는데, 이는 입사 제트·플룸과 유사한 에너지 감쇠 메커니즘이다. 결과적으로 계면이 충분히 약해지면 표면 파동이 더 이상 충분한 전단을 제공하지 못해 혼합 속도가 감소하고, 최종적으로는 계면이 표면 강제와 완전히 분리(decoupling)되는 정착 상태에 도달한다.

이러한 현상을 설명하기 위해 저자들은 깊이(z)와 내부 리차드슨 수(Ri(z)=N²/(∂u/∂z)²) 에 의존하는 난류 확산계수 κ_t(z)=κ₀ exp