Effect of velocity, fluid properties and drop shape on coalescence and neck oscillation

본 연구는 축대칭 수치 시뮬레이션을 통해 액적의 속도, 유체 특성 및 모양이 액적 합체 및 목 (neck) 진동에 미치는 영향을 분석하여 부분 합체와 완전 합체의 경계를 규명하고, 타원형 액적에서 2 차 액적 형성이 가장 두드러지며 레이leigh-Plateau 불안정성이 주요 동인이 아님을 밝혔습니다.

핵심

🌧️ 1. 연구의 배경: 물방울의 두 가지 운명

비 오는 날, 물웅덩이에 물방울이 떨어지는 장면을 상상해 보세요.

• 운명 A (완전 합체): 물방울이 물웅덩이와 완전히 섞여 사라집니다. (완전 합체)
• 운명 B (부분 합체): 물방울이 물웅덩이와 닿자마자, 마치 요술처럼 작은 물방울 하나를 만들어 위로 솟구친 뒤 다시 떨어집니다. (부분 합체)

과학자들은 오랫동안 "왜 어떤 물방울은 그냥 합쳐지고, 어떤 것은 작은 자식을 낳고 사라질까?"를 궁금해했습니다. 이 논문은 그 비밀을 풀기 위해 컴퓨터로 수천 번의 실험을 진행했습니다.

🎭 2. 주요 등장인물: 물방울의 '성격'과 '환경'

연구자들은 물방울이 어떻게 행동하는지 결정하는 네 가지 주요 요소를 분석했습니다.

1. 물방울의 모양 (타원형 vs 구형):

   • 물방울이 떨어지기 전 모양이 세로로 길쭉한지 (계란형), 가로로 납작한지 (접시형), 아니면 완전한 구형인지에 따라 결과가 달라집니다.
   • 비유: 세로로 길쭉한 물방울은 마치 긴 스프링처럼 행동해서, 합쳐질 때 더 쉽게 작은 물방울을 만들어냅니다. 반면 납작한 물방울은 바닥에 넓게 퍼지면서 바로 합쳐지려 합니다.

2. 떨어지는 속도 (충격력):

   • 물방울이 빠르게 떨어질수록 (높은 속도) 물웅덩이로 더 많이 흘러들어갑니다.
   • 비유: 천천히 떨어지면 물방울이 "잠시 멈추고" 작은 자식을 낳을 시간이 있지만, 너무 빠르게 떨어지면 물웅덩이로 쏙 빠져버려 합쳐지는 게 더 빨라집니다.

3. 점성 (끈적임):

   • 물이 얼마나 **끈적한지 (꿀처럼 끈적한지 vs 물처럼 묽은지)**도 중요합니다. 끈적하면 물방울이 합쳐지는 속도가 느려져서 작은 물방울을 만들기 쉽습니다.

4. 중력과 표면 장력:

   • 중력은 물방울을 아래로 당기고, 표면 장력은 물방울을 둥글게 유지하려는 힘입니다. 이 두 힘의 싸움 결과가 합체 방식을 결정합니다.

🌊 3. 발견된 비밀: '목 (Neck)'의 춤과 진동

가장 흥미로운 발견은 물방울이 합쳐지는 순간, 물방울과 물웅덩이 사이에서 생기는 '목 (Neck)' 부분의 움직임입니다.

• 목의 진동 (Neck Oscillation):
  물방울이 닿으면 물이 위로 솟구치는데, 이때 목 부분이 마치 스프링처럼 수축하고 팽창하며 진동합니다.
  • 한 번 진동 후 합쳐짐: 물방울이 바로 합쳐집니다.
  • 두 번 진동 후 합쳐짐: 물방울이 합쳐지는 과정이 조금 더 길어지거나, 작은 물방울이 만들어집니다.
  • 연구자들은 이 진동 횟수를 세어 "1 단계 부분 합체", "2 단계 부분 합체" 등으로 새로운 분류를 만들었습니다. 마치 춤을 추는 단계에 따라 결과가 달라진다고 볼 수 있죠.

🚫 4. 오해를 깨다: "레이리 - 플레이트 불안정성"은 죄가 없다?

과거 과학자들은 작은 물방울이 만들어지는 원인을 '레이리 - 플레이트 불안정성' (기름기 있는 물줄기가 여러 조각으로 부서지는 현상) 때문이라고 생각했습니다. 마치 호스 물줄기가 여러 방울로 부서지는 것처럼요.

하지만 이 연구는 **"아니요, 그건 주된 원인이 아닙니다"**라고 말합니다.

• 비유: 물줄기가 부서지는 건 부수적인 현상일 뿐, 진짜 원인은 물방울이 위로 솟구치는 힘과 아래로 떨어지는 힘의 싸움에서 비롯된다고 설명합니다. 특히 물방울이 세로로 길쭉할 때만 이 불안정성이 작은 물방울을 여러 개로 쪼개는 데 도움을 줍니다.

📊 5. 결론: 물방울의 지도 만들기

연구자들은 이 모든 요소 (속도, 점성, 중력, 모양) 를 조합하여 **3 차원 지도 (Phase Diagram)**를 만들었습니다.

• 이 지도를 보면, "이런 조건에서는 작은 물방울이 생기고, 저런 조건에서는 그냥 합쳐진다"를 정확히 예측할 수 있습니다.
• 특히 세로로 길쭉한 물방울이 가장 작은 물방울을 잘 만들어내고, 빠르게 떨어지거나 끈적한 액체일수록 합쳐지는 경향이 강해진다는 사실을 밝혀냈습니다.

💡 요약

이 논문은 **"물방울이 물에 떨어질 때, 왜 어떤 것은 작은 자식을 낳고 사라지는가?"**에 대한 답을 찾았습니다.
그 답은 물방울의 모양, 떨어지는 속도, 그리고 목 부분의 진동에 달려 있었습니다. 마치 물방울이 물웅덩이와 춤을 추는데, 그 춤의 리듬과 스타일에 따라 결과가 달라진다는 것이죠. 이 연구는 잉크젯 프린팅, 비가 땅에 떨어지는 현상, 기름과 물 분리 등 다양한 분야에서 유용하게 쓰일 수 있습니다.

기술 요약

1. 연구 배경 및 문제 정의 (Problem)

• 배경: 액적 (droplet) 이 액체 풀 (liquid pool) 에 충돌하여 병합 (coalescence) 되는 현상은 분무, 연소, 잉크젯 프린팅, 유체 분리 등 다양한 공학적 응용 및 자연 현상 (빗방울 충돌 등) 에서 중요합니다.
• 핵심 현상: 충돌 후 얇은 공기막이 파열되면 두 액체가 접촉하며 병합이 시작됩니다. 이때 두 가지 결과가 발생할 수 있습니다.
  1. 부분 병합 (Partial Coalescence): 액적이 완전히 흡수되지 않고, 수직 기둥이 형성된 후 목 (neck) 이 조여져 2 차 액적 (secondary droplet) 이 분리됩니다.
  2. 완전 병합 (Complete Coalescence): 액적이 풀로 완전히 흡수되어 2 차 액체가 형성되지 않습니다.
• 연구 동기: 기존 연구들은 점성 (Ohnesorge 수, Oh) 과 중력 (Bond 수, Bo) 의 영향을 주로 다루었으나, 관성 (Weber 수, We) 의 영향과 낙하 방울의 모양 (구형, 편평형, 길쭉형) 이 병합 과정, 특히 목 진동 (neck oscillation) 과 2 차 액체 형성 메커니즘에 미치는 체계적인 영향은 명확히 규명되지 않았습니다. 또한, Rayleigh-Plateau 불안정성이 부분 병합의 주된 원인인지에 대한 논쟁이 존재합니다.

2. 방법론 (Methodology)

• 수치 해석: 축대칭 (axisymmetric) 수치 시뮬레이션을 수행했습니다.
• 솔버: 오픈소스 유체 역학 솔버인 Gerris를 사용했습니다. 이는 유한 체적법 (Finite-Volume Method) 기반이며, VOF (Volume-of-Fluid) 방법을 사용하여 공기 - 액체 계면을 추적합니다.
• 격자 및 검증:
  • 동적 적응 격자 세분화 (Dynamic Adaptive Mesh Refinement) 를 사용하여 계면 근처와 와도 (vorticity) 가 높은 영역의 해상도를 높였습니다.
  • 기존 실험 데이터 (Blanchette & Bigioni, 2006; Tran et. al., 2013) 와 비교하여 격자 수렴성 및 결과의 정확성을 검증했습니다.
• 제어 변수:
  • 무차원 수: Weber 수 (We, 관성/표면장력), Ohnesorge 수 (Oh, 점성/관성 및 표면장력), Bond 수 (Bo, 중력/표면장력).
  • 방울 모양: 종횡비 (Aspect Ratio, AR) 를 변화시켜 편평형 (Oblate, AR<1), 구형 (Spherical, AR=1), 길쭉형 (Prolate, AR>1) 을 구현했습니다.
  • 영역: 깊은 액체 풀을 가정하고, 초기 방울과 풀 표면 사이에 0.1Deq 의 간격을 두고 충돌을 시뮬레이션했습니다.

3. 주요 기여 및 결과 (Key Contributions & Results)

가. 3 차원 위상도 (Phase Diagram) 구축

• We, Oh, Bo 의 3 차원 공간에서 부분 병합과 완전 병합의 경계를 정의하는 위상도를 제시했습니다.
• 관찰:
  • Oh 증가: 점성력이 증가하면 표면장력에 의한 교란이 억제되어 수직 붕괴가 가속화되므로 완전 병합이 촉진됩니다.
  • Bo 증가: 중력이 지배적이 되어 액적의 배수가 빨라지며 완전 병합으로 전환됩니다.
  • We 증가: 관성 효과가 커지면 액적의 풀로의 배수가 가속화되어 수직 붕괴율이 수평 수축율보다 빨라지며, 이는 완전 병합을 유도합니다.
  • 임계값: 낮은 We 조건에서 완전 병합을 위한 임계 Oh 값 ($Oh_c$) 은 약 0.021 로 추정되었으며, 이는 기존 연구 결과와 일치합니다.

나. 목 진동 (Neck Oscillation) 과 병합 단계 분류

• 부분 병합과 완전 병합 사이의 전이 영역에서 목 진동이 발생하는 것을 발견하고, 이를 기준으로 4 가지 병합 regimes 를 체계적으로 분류했습니다.
  1. 1 단계 부분 병합: 첫 번째 진동 주기 동안 2 차 액체가 분리됨.
  2. 2 단계 부분 병합: 두 번째 진동 주기 동안 2 차 액체가 분리됨.
  3. 1 단계 완전 병합: 첫 번째 진동 주기 동안 완전히 흡수됨.
  4. 2 단계 완전 병합: 두 번째 진동 주기 동안 완전히 흡수됨.
• 메커니즘: 병합의 결과는 수평 방향의 수축 운동량 (목이 조여지는 힘) 과 수직 방향의 붕괴 운동량 (액주가 아래로 내려가는 힘) 간의 경쟁에 의해 결정됩니다.
  • 수평 수축이 수직 붕괴보다 빠르면 2 차 액체가 분리됩니다.
  • We, Oh, Bo 가 증가할수록 수직 붕괴가 가속화되어 완전 병합으로 전환됩니다.

다. 방울 모양 (Shape) 의 영향

• 길쭉형 (Prolate, AR>1): 목이 좁게 형성되어 배수가 지연되고 수평 수축이 우세하여 2 차 액체 형성 확률이 가장 높습니다.
• 편평형 (Oblate, AR<1): 넓은 목이 형성되어 배수가 빨라지고 수직 붕괴가 촉진되므로 완전 병합이 쉽게 일어납니다.
• 구형 (Spherical): 두 극단 사이의 중간 성향을 보입니다.
• 다중 2 차 액체 형성: 길쭉형 방울의 경우, 형성된 액주가 길게 늘어져 Rayleigh-Plateau 불안정성을 유발하여 2 차 액체가 다시 여러 개의 작은 액적으로 분열되는 현상이 관찰되었습니다.

라. Rayleigh-Plateau 불안정성의 역할 재검토

• 기존 일부 연구는 Rayleigh-Plateau 불안정성이 부분 병합의 주원인이라고 주장했으나, 본 연구는 액주 높이와 목 직경의 비율 ($\psi = h/\pi w_n$) 분석을 통해 이를 반박했습니다.
• 결론: 본 연구에서 조사된 파라미터 공간 내에서는 Rayleigh-Plateau 불안정성이 부분 병합을 주도하지 않으며, 오히려 수평/수직 운동량의 경쟁과 목 진동이 핵심 메커니즘임을 규명했습니다. 불안정성은 주로 2 차 액체가 형성된 후 추가 분열을 일으킬 때만 관여합니다.

4. 의의 및 결론 (Significance)

• 이론적 기여: 관성, 점성, 중력, 그리고 방울 모양이 복합적으로 작용하는 병합 역학에 대한 포괄적인 이해를 제공했습니다. 특히, 단순한 2 차원 위상도를 넘어 3 차원 (We-Oh-Bo) 공간에서의 전이 경계를 명확히 했습니다.
• 메커니즘 규명: 목 진동 (Neck oscillation) 이 병합 단계 (1 단계/2 단계, 부분/완전) 를 구분하는 강력한 지표임을 제시했습니다.
• 응용 가능성: 미세유체공학, 연소 효율 개선, 스프레이 공정 최적화 등 액적 - 액체 상호작용이 중요한 분야에서 예측 모델을 개선하는 데 기여할 수 있습니다.
• Rayleigh-Plateau 설의 수정: 부분 병합의 주된 동인이 Rayleigh-Plateau 불안정성이 아님을 증명하여 해당 분야의 기존 이론적 논쟁에 중요한 통찰을 제공했습니다.

이 연구는 다양한 유체 조건과 방울 형태 하에서 액적 병합의 복잡한 역학을 체계적으로 규명하여, 액적 충돌 현상에 대한 예측 정확도를 높이고 새로운 제어 전략을 모색하는 데 기초 자료를 제공합니다.