디스크 충돌에 의한 물 표면 공동 역학

초록

본 연구는 선형 모터로 속도를 정밀 제어한 원판이 물 표면에 충돌할 때 형성되는 일시적인 공기 공동의 형성·붕괴 과정을 고속 영상, 입자 이미지 유속(PIV) 및 경계 적분 시뮬레이션으로 조사한다. 프루드 수(Fr)가 0.6~200 범위에서 공동의 반경·깊이·최소 반경·기포 부피가 Fr¹⁄²에 비례하는 스케일링을 보이며, 무한 원통 붕괴 모델이 실험·수치와 매우 좋은 일치를 보여준다. 또한 공기 밀봉(surface seal) 현상이 프루드 수가 높을 때 발생해 공동 붕괴를 방해함을 확인한다.

상세 분석

이 논문은 디스크 반경 h₀와 속도 V를 선형 모터로 정확히 제어함으로써 Froude 수 Fr = V²/(g h₀)를 0.6 ~ 200 사이에서 조절한다. 실험에서는 고속 카메라(1000 fps)로 공동의 형성·수축 과정을 촬영하고, PIV를 이용해 공동 주변의 축대칭 유동장을 측정하였다. 수치적으로는 무점성·무회전 가정하에 경계 적분법을 적용했으며, 축대칭성을 이용해 1차원 라인 적분으로 변환하였다. 시간 스텝은 현지 속도와 노드 간 거리에 기반한 적응형 크루크-니콜슨 스킴으로 결정하고, 고곡률 영역에서는 노드 밀도를 자동으로 증가시켜 핀치‑오프 근처의 미세 구조를 정확히 포착한다.

실험 영상과 시뮬레이션 결과는 Fr ≤ 13.6 정도에서는 거의 차이가 없으며, 특히 공동의 외곽 형태와 수축 속도가 정량적으로 일치한다. 그러나 Fr ≈ 13.6 ~ 200 구간에서 공기 흐름에 의한 surface seal 현상이 나타나, 실험에서는 공동 상부가 공기 흐름에 의해 조기에 봉인되어 공동 깊이가 얕아지는 반면, 시뮬레이션은 공기 모델을 배제했기 때문에 이러한 현상이 재현되지 않는다. 이는 프루드 수가 클수록 공기 역학이 공동 붕괴에 중요한 역할을 함을 시사한다.

공동 반경의 시간 진화는 무한 원통 붕괴 모델을 기반으로 한 간단한 ODE : R̈ = −g · (R₀/R)² 로 기술된다. 이 모델은 확장 단계와 수축 단계 사이의 비대칭성을 자연스럽게 설명한다. 모델 해석을 통해 폐쇄 깊이 z_c와 전체 깊이 z_max이 모두 ∝ Fr¹⁄²·h₀·f(Fr) 형태로 스케일링됨을 보였으며, 여기서 f(Fr) 는 프루드 수에 약하게 의존하는 전처리 상수이다. 최소 반경 R_min의 시간적 감소율은 t → t_c 근처에서 R_min ∝ (t_c − t)^{2/3} 와 같은 특이 지수를 갖는다. 또한, 공기 기포 부피 V_bubble는 V_bubble / h₀³ ≈ (1 + 0.26 Fr^{1/2}) Fr^{1/2} 로 정량화되었다. 이는 실험적으로 측정된 기포 부피와 시뮬레이션 결과가 모두 이 식에 잘 부합함을 의미한다.

이러한 결과는 기존 연구에서 제시된 “무한 원통 붕괴” 가정이 실제 물-공기 시스템에서도 유효함을 입증한다. 특히, 프루드 수가 중간 범위(Fr ~ 1 ~ 100)에서는 점성·표면장력 효과가 무시될 정도로 Reynolds 수와 Weber 수가 충분히 커서, 관성·중력 주도 흐름이 지배한다는 점을 확인하였다. 반면, 프루드 수가 매우 높아지면 공기 동역학(표면 밀봉, 공기 유입)과 미세한 캡illary 파동이 비정상적인 공동 형태를 유발한다는 점에서, 향후 연구에서는 공기-액체 상호작용을 포함한 전산 유체역학(CFD) 모델이 필요함을 제시한다.