Internal flow and concentration in neighbouring evaporating binary droplets and rivulets

이 논문은 잉크젯 프린팅 및 분무 냉각 등 다양한 응용 분야에서 관찰되는 인접한 이성분 액적과 리블렛의 증발 과정에서 접촉각, 간격, 마랑고니 수와 같은 주요 인자들이 유동 및 농도 분포의 대칭성과 정지점에 미치는 영향을 수치 및 이론적으로 규명하였다.

핵심

1. 핵심 상황: "두 친구가 나란히 서서 땀을 흘리는 상황"

상상해 보세요. 두 명의 친구 (방울) 가 나란히 서서 뜨거운 햇빛 아래에서 땀을 흘리고 있습니다.

• 혼자 있을 때: 한 친구가 땀을 흘리면, 그 친구 주변 공기는 습해지지만, 다른 친구가 없으니 땀이 잘 날아갑니다.
• 함께 있을 때: 두 친구가 너무 가까이 서 있으면, 서로가 땀을 날려보낸 습한 공기가 서로의 땀을 막아줍니다. 이를 **'가림 효과 (Shielding Effect)'**라고 합니다.
  • 결과: 두 친구 모두 땀이 더 천천히 마르게 되고, 서로의 땀 흐름이 서로에게 영향을 미쳐 흐름이 비대칭적으로 변합니다.

이 연구는 바로 이 '가림 효과'가 액체 (방울) 내부의 흐름을 어떻게 뒤틀어 놓는지를 수학적으로 분석한 것입니다.

2. 액체 내부의 두 가지 흐름: "커피 스테인" vs "마랑고니 흐름"

방울 안에는 두 가지 힘이 서로 싸우고 있습니다.

1. 커피 스테인 효과 (Coffee-stain effect):
   • 비유: 커피가 마르면 가장자리 (테두리) 로 커피 입자가 몰려가서 고리 모양으로 남는 현상입니다.
   • 원리: 방울 가장자리에서 수분이 더 빨리 증발하기 때문에, 액체가 가장자리로 끌려가면서 입자들을 밀어냅니다.
2. 마랑고니 흐름 (Marangoni flow):
   • 비유: 비눗방울에 알코올을 떨어뜨리면 비눗방울이 급격히 변형되거나 회전하는 것처럼, 액체 표면의 '표면 장력' 차이 때문에 생기는 흐름입니다.
   • 원리: 이 연구에서는 물과 1,2-헥산디올 (비휘발성 성분) 이 섞인 액체를 다룹니다. 물이 먼저 증발하면, 비휘발성 성분이 가장자리에 남게 되어 표면 장력이 변합니다. 이 차이가 마치 액체 내부에 강력한 바람을 불어넣어 액체를 강하게 회전시킵니다.

연구의 핵심: 보통 '커피 스테인' 효과가 지배적이지만, 이 연구에서는 마랑고니 흐름이 훨씬 강력해서 액체 내부의 흐름을 완전히 뒤집어 놓는 상황을 다룹니다.

3. 두 가지 실험: "긴 물줄기 (리뷰렛)" vs "둥근 방울 (드롭렛)"

연구진은 복잡한 3 차원 문제를 단순화하기 위해 두 가지 모양을 비교했습니다.

A. 긴 물줄기 (리뷰렛, Rivulets) - "2 차원 길쭉한 물방울"

• 상황: 두 개의 긴 물줄기가 나란히 있는 상태입니다.
• 발견:
  • 물줄기 내부의 흐름은 대칭이 깨져서 한쪽으로 쏠리게 됩니다.
  • 중요한 점: 이 흐름의 중심 (정체점) 이 어디로 이동하는지는 물줄기 사이의 거리와 **방울의 각도 (접촉각)**에만 의존합니다.
  • 놀라운 사실: 마랑고니 힘 (표면 장력 차이) 이 아무리 세게 불어도, 물줄기 내부의 흐름 중심은 거의 변하지 않습니다. 마치 거대한 기차 (마랑고니 흐름) 가 달리고 있어도, 기차의 방향은 궤도 (거리와 각도) 에만 결정되는 것과 같습니다.

B. 둥근 방울 (드롭렛, Droplets) - "3 차원 구형 물방울"

• 상황: 일반적인 둥근 물방울 두 개가 나란히 있는 상태입니다.
• 발견:
  • 물줄기와는 다르게, 둥근 방울은 3 차원 공간에서 움직일 수 있습니다 (위아래, 좌우, 그리고 회전).
  • 중요한 차이: 둥근 방울은 마랑고니 힘의 세기에 따라 흐름의 중심이 움직입니다.
  • 이유: 둥근 방울은 물줄기와 달리 **회전 방향 (방사형이 아닌 원주 방향)**으로도 흐를 수 있는 '여분의 자유도'가 있기 때문입니다. 마랑고니 힘이 강해질수록 이 회전 흐름이 더 복잡해지고, 그 결과 흐름의 중심이 마랑고니 힘의 세기에 따라 달라집니다.

4. 연구의 결론과 의미

이 논문은 복잡한 수학적 모델과 컴퓨터 시뮬레이션을 통해 다음과 같은 결론을 내렸습니다.

1. 가림 효과가 핵심: 인접한 방울들은 서로의 증발을 방해하여 흐름을 비대칭으로 만듭니다.
2. 모양이 중요:
   • 길쭉한 물줄기라면, 흐름의 중심은 거리의 영향을 많이 받지만 마랑고니 힘의 세기에는 무관합니다.
   • 둥근 방울이라면, 흐름의 중심은 마랑고니 힘의 세기에 따라 민감하게 반응합니다.
3. 실용적 가치: 이 연구는 잉크젯 프린터가 잉크를 얼마나 정밀하게 뿌릴지, 농약이 식물에 얼마나 골고루 퍼질지 예측하는 데 도움을 줍니다. 특히 여러 방울이 뭉쳐 있을 때 어떤 현상이 일어나는지 이해함으로써, 더 나은 기술 개발이 가능해집니다.

요약: 한 문장으로 정리하면?

"두 개의 액적 (방울) 이 나란히 있을 때, 서로가 서로의 증발을 방해하며 내부 흐름을 뒤틀어 놓는데, 이 흐름의 중심이 어디로 갈지는 액적의 모양 (길쭉한지 둥근지) 에 따라 마랑고니 힘의 영향을 받거나 받지 않는다는 것을 밝혀냈다."

이 연구는 마치 복잡한 춤 (유체 흐름) 을 추는 두 명의 무용수가 서로의 동작을 어떻게 방해하고 조율하는지, 그리고 그들이 **긴 줄을 타고 춤추는지 (리뷰렛) 공중에서 춤추는지 (드롭렛)**에 따라 춤의 패턴이 어떻게 달라지는지를 수학적으로 해독한 것이라고 볼 수 있습니다.

기술 요약

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

• 배경: 잉크젯 프린팅, 스프레이 냉각, 농약 전달 등 다양한 공학 응용 분야에서 인접한 다성분 액적 (droplets) 이나 리벳 (rivulets, 좁은 띠 모양의 액체) 의 증발 현상이 빈번하게 발생합니다.
• 핵심 문제:
  • 차폐 효과 (Shielding Effect): 인접한 액체들 사이의 거리가 가까우면 증기 농도가 높아져 개별 증발 속도가 감소합니다.
  • 비대칭성 (Asymmetry): 이 차폐 효과는 증발 속도의 불균형을 초래하여, 액체 내부의 농도 분포와 유동장의 대칭성을 파괴합니다.
  • 복잡성: 단일 액적의 증발은 잘 연구되었으나, 인접한 다성분 (예: 물과 1,2-헥사디올) 액적의 증발 시 발생하는 마랑고니 (Marangoni) 유동과 선택적 증발의 상호작용은 3 차원적이며 계산 비용이 매우 높아 분석이 어렵습니다.
• 목표: 인접한 이원계 액적/리벳의 증발 과정에서 내부 유동과 농도 분포의 비대칭성을 결정하는 주요 물리적 인자 (접촉각, 액적 간 거리, 마랑고니 수) 를 규명하고, 이를 효율적으로 모델링하는 것입니다.

2. 방법론 (Methodology)

연구팀은 수치 시뮬레이션과 이론적 모델을 결합하여 다음과 같은 접근법을 사용했습니다.

• 물리적 가정:
  • 비휘발성 성분 (A) 과 휘발성 성분 (B) 으로 구성된 이원계 혼합물 (물 + 1,2-헥사디올) 을 가정.
  • 증발이 느리고, 마랑고니 유동이 충분히 강하여 농도 균질화가 일어난다고 가정 (준정상, Quasi-stationary).
  • 온도 변화는 무시하고 조성 변화에 따른 표면 장력 변화만 고려.
• 수치 모델:
  1. 완전 과도 수치 시뮬레이션 (Transient DNS): 2 차원 리벳에 대해 Navier-Stokes 방정식, 확산 - 대류 방정식, 증기 확산 방정식을 풀어서 기준 데이터를 생성.
  2. 준정상 모델 (Quasi-stationary Model): 시간에 따른 변화가 느리다고 가정하여 방정식을 단순화. 마랑고니 수 (Ma), 접촉각 ($\theta$), 리벳 간 거리 ($b$) 를 무차원 파라미터로 사용.
  3. 윤활 근사 모델 (Lubrication Model): 접촉각이 작다는 가정 하에 3 차원 문제를 2 차원 (수직 평균) 문제로 축소. Taylor-Aris 분산 효과를 포함하여 농도 혼합을 모델링.
• 분석 지표:
  • 정체점 (Stagnation Point, $x_0$): 액적 내부의 두 와류 (vortex) 가 만나는 지점 (접촉면에서의 접선 속도 0 지점). 이 점의 위치 이동을 통해 유동장의 비대칭성 정도를 정량화.

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

가. 모델 검증 (Validation)

• 준정상 모델의 정확성: 과도 (Transient) 시뮬레이션 결과와 준정상 모델 결과를 비교한 결과, 마랑고니 수가 충분히 큰 조건에서 두 모델은 매우 높은 일치도를 보였습니다. 이는 마랑고니 유동이 '커피 스테인 (coffee-stain)' 효과를 압도하여 농도 분포가 안정화됨을 의미합니다.
• 윤활 근사의 한계: 리벳의 경우 윤활 근사 모델이 초기 단계 (큰 접촉각) 에서는 증발 플럭스 (evaporation flux) 의 접촉각 의존성을 제대로 반영하지 못해 정체점 위치를 정확히 예측하지 못했으나, 접촉각이 작아지는 후기 단계에서는 잘 일치했습니다.

나. 리벳 (Rivulets) 의 유동 특성

• 마랑고니 수 (Ma) 의 비의존성: 2 차원 리벳의 경우, 정체점의 위치 ($x_0$) 는 마랑고니 수 (Ma) 에 거의 의존하지 않습니다. 오직 접촉각 ($\theta$) 과 리벳 간 거리 ($b$) 에만 의존합니다.
• 접촉각과 거리의 영향:
  • 접촉각이 증가하거나 리벳 간 거리가 감소할수록 차폐 효과가 강화되어 정체점이 중심에서 더 멀리 이동 (비대칭성 증가) 합니다.
  • 이론적 분석 (윤활 근사 + 해석적 유도) 을 통해 정체점 이동량이 농도 최대값 위치와 일치하며, 마랑고니 수와 무관함을 증명했습니다.

다. 액적 (Droplets) 의 유동 특성

• 마랑고니 수 (Ma) 의 의존성: 3 차원 액적의 경우, 리벳과 달리 정체점 위치가 마랑고니 수 (Ma) 에 강하게 의존합니다.
• 원인: 액적은 리벳과 달리 방위각 (azimuthal) 방향으로 유동이 발생할 수 있는 추가적인 자유도가 있습니다. 이로 인해 농도 구배의 비선형적 진화가 발생하고, 마랑고니 수의 변화가 유동 구조와 정체점 위치에 복잡한 영향을 미칩니다.
• 수직 구배: Ma 가 증가함에 따라 액적 표면과 기판 근처의 정체점 위치가 서로 다르게 이동하며, 이는 수직 방향의 큰 속도 구배를 의미합니다.

4. 의의 및 결론 (Significance & Conclusions)

• 계산 효율성 확보: 3 차원 과도 시뮬레이션의 고비용 문제를 해결하기 위해, 2 차원 리벳 모델을 통해 핵심 물리 현상을 규명하고 이를 기반으로 단순화된 준정상/윤활 모델을 개발했습니다.
• 물리적 통찰:
  • 리벳: 유동 비대칭성은 기하학적 인자 (접촉각, 거리) 에 의해 지배되며, 마랑고니 효과는 혼합을 촉진하지만 비대칭성 위치에는 영향을 주지 않음.
  • 액적: 3 차원성 (방위각 유동) 으로 인해 마랑고니 수가 유동 구조와 정체점 이동에 결정적인 역할을 함.
• 실용적 함의: 인접한 액적/리벳의 증발 거동을 예측할 때, 시스템의 차원 (2D 리벳 vs 3D 액적) 에 따라 마랑고니 수를 고려해야 할지 여부를 판단할 수 있는 기준을 제시했습니다. 이는 잉크젯 프린팅의 품질 제어, 스프레이 냉각 효율 최적화, 농약 코팅 균일성 향상 등에 기여할 수 있습니다.
• 한계 및 향후 과제: 현재 윤활 모델은 접촉각이 0 인 평평한 액적에 대한 증발 플럭스 식을 사용하므로, 유한한 접촉각을 가진 인접 액적의 정확한 증발 플럭스 해석식을 개발하는 것이 향후 중요한 과제로 남았습니다.

이 논문은 인접한 증발 액적 시스템의 복잡한 유동 - 확산 - 표면 장력 상호작용을 체계적으로 분해하여, 각 시스템 (리벳 vs 액적) 에 지배적인 물리 인자가 어떻게 다른지를 명확히 규명한 중요한 연구입니다.