기술 요약: 액적-풀 상호작용의 탄성 유체 역학 (Elasto-hydrodynamics)

문제 정의

뉴턴 유체의 경우 액적이 액체 풀에 충돌할 때 발생하는 유체 역학(공동 형성, 크라운 확장, 모세관 파동 및 워싱턴 제트 발생 등)은 잘 규명되어 있으나, 탄성(또는 보거, Boger) 유체의 경우 그 거동에 대한 연구는 거의 이루어지지 않았다. 기존 문헌은 주로 뉴턴 시스템이나 통제된 기하학적 구조 내에서의 점탄성 필라멘트 가늘어짐(viscoelastic filament thinning)에 집중되어 있다. 충돌하는 액적과 대상이 되는 풀이 모두 탄성 고분자 유체일 때, 공동 형성, 붕괴 및 후속 제트 분출의 비정상적이고 공간적으로 진화하는 역학을 이해하는 데에는 상당한 공백이 존재한다. 특히, 관성, 모세관 및 탄성력 사이의 상호작용이 이러한 사건들의 시공간적 진화를 지배하는 방식과 그로 인해 발생하는 독특한 형태학적 체계(morphological regimes)를 규명하기 위한 체계적인 조사가 필요하다.

방법론

본 연구는 실험적, 이론적, 계산적 접근 방식을 결합한 삼중 구조를 채택한다:

1. 실험 설정:

   • 유체: 고분자량 폴리에틸렌 옥사이드(PEO, ~5–6 백만 Da)의 수용액을 농도별로(0, 250, 500, 1000, 2000 ppm) 사용하여 탄성(보거) 유체를 제조하였다. 대조군으로는 순수 물과 물-글리세롤 혼합물(탄성 유체와 점도를 맞춘 것)을 사용하였다.
   • 절차: 액적을 준정상 상태(quasi-steadily)로 생성하여 깊은 풀 위로 자유 낙하시켰다. 충돌 속도는 1.5에서 4.3 m/s 범위였다.
   • 진단: 고속 촬영(4000 fps)을 통해 공동 확장, 붕괴, 크라운 형성 및 제트 분출의 시간적 진화를 포착하였다.
   • 무차원 매개변수: 역학적 분석은 Weber 수($We$)와 Deborah 수($De$)의 범위 내에서 수행되었으며, 구형의 충돌 전 액적을 보장하기 위해 Bond 수($Bo$)는 0.27로 고정하였다. Ohnesorge 수($Oh$)는 0.0046 미만으로 유지하여 점성 효과가 무시될 수 있도록 하였다.

2. 이론적 프레임워크:

   • 공동 진화: 포텐셜 흐름(potential-flow) 프레임워크 내에서 에너지 보존 모델을 도출하였다. 초기 액적의 운동 에너지는 중력 위치 에너지, 표면 에너지, 풀의 운동 에너지, 그리고 저장된 탄성 에너지 항의 균형을 맞추도록 설계되었다.
   • 제트 가늘어짐: FENE-P(Finitely Extensible Nonlinear Elastic - Peterstein) 구성 방정식 모델을 사용하여 워싱턴 제트 반경의 진화를 설명하는 이론적 모델을 개발하였다. 이 모델은 고분자 사슬의 유한한 신장성과 탄성 및 모세관 응력 사이의 균형을 고려한다.

3. 수치 시뮬레이션:

   • 솔버: 시뮬레이션은 COMSOL Multiphysics를 사용하여 수행되었다.
   • 지배 방정식: 계면을 추적하고 탄성 응력을 모델링하기 위해 Cahn-Hilliard 페이즈 필드(phase-field) 모델을 FENE-P 구성 방정식과 결합하였다.
   • 검증: 시뮬레이션 결과는 공동 형태, 제트 높이, 속도장을 중심으로 뉴턴 및 탄성 사례 모두에 대해 실험 데이터와 비교 검증되었다.

주요 기여 및 결과

1. 점성에 대한 탄성의 지배력

탄성 유체(2000 ppm PEO)와 점도가 일치하는 뉴턴 유체(물-글리세롤) 간의 비교 실험을 통해, 관찰된 형태학적 변화가 점성이 아닌 순수하게 유체의 탄성에 의해 발생함을 입증하였다.

• 공동 형태: 뉴턴 유체는 반구형 공동을 생성하는 반면, 탄성 유체는 수축 단계에서 뚜렷한 사다리꼴 공동 형태를 보인다. 이는 비균일한 고분자 신장과 그에 따른 탄성 응력 분포 때문이며, 이로 인해 모세관 파동의 전파가 지연되고 공동 바닥이 평평해진다.
• 점도 독립성: 본 연구는 테스트된 범위 내에서 점성이 파쇄 과정에 미치는 영향은 미미하며, 탄성 응력 분포가 공동 형태와 제트 역학을 결정한다는 것을 확인하였다.

2. 에너지 전환 및 저장

에너지 보존 모델을 통해 운동 에너지가 탄성 에너지로 전환되는 과정을 정량화하였다.

• 탄성 저장: 특정 고농도 케이스에서 고분자 사슬의 신장을 통해 초기 운동 에너지의 **약 35%**가 탄성 에너지로 저장될 수 있음을 나타낸다. 구체적으로, 500D-500P 케이스에서 에너지 분할 계수($\beta$)가 뉴턴 유체의 0.80에서 탄성 유체의 0.45로 감소하는 것은 초기 운동 에너지의 약 **35%**가 탄성 에너지로 저장되었음을 의미한다.
• 역학에 미치는 영향: 이러한 에너지 저장은 공동 확장에 사용할 수 있는 에너지를 감소시켜, 뉴턴 유체에 비해 탄성 유체에서 최대 공동 크기가 작아지고 성장 속도가 느려지게 만든다. 모델은 에너지 분할 계수($\beta$)를 조정함으로써 평균 공동 반경의 시간적 진화를 성공적으로 예측하였다.

3. 레짐 맵(Regime Maps) 및 형태학적 전이

$[De, We]$ 매개변수 공간에서 네 가지 뚜렷한 레짐을 확립하였다:

1. 핀치오프(Pinch-off): 이차 액적의 완전한 분리 (주로 뉴턴 시스템에서 관찰됨).
2. 위성 액적 + 필라멘트(Satellite Droplet + Filament): 위성 액적이 필라멘트에 의해 연결된 형태. 여기에는 다음이 포함된다:
   • 비즈 온 어 스트링(Beads-on-a-String, BOAS): Rayleigh 불안정성 성장률이 역 이완 시간(inverse relaxation time)보다 느리거나 비슷할 때 관찰되며, 고분자 신장이 지배하기 전에 비대칭성이 발달할 수 있게 한다.
   • 가는/굵은 필라멘트(Thin/Thick Filaments): 탄성 응력이 모세관 압력과 균형을 이루어 완전한 파쇄를 방지할 때 관찰된다.
3. 제트 없음(No Jet): 낮은 $We$또는 높은$De$에서 모세관 파동의 수렴 지연으로 인해 제트 형성이 억제됨.
4. 핀치오프 없음(No Pinch-off): 높은 $De$에서 탄성 응력이 모세관 수축(necking)을 지속적으로 상쇄하여, 파쇄되지 않는 안정적인 워싱턴 제트가 형성됨.

4. 탄성-모세관 제트 가늘어짐

워싱턴 제트의 가늘어짐에 대한 이론적 및 실험적 분석은 다음과 같다:

• 지수적 감쇠: 탄성-모세관 레짐에서 제트 반경은 시간 $t$에 따라 $r(t) \sim \exp(-t/2\lambda)$의 스케일링 법칙을 따르며 지수적으로 감소한다.
• 스케일링 차이: 이는 표준 점탄성 필라멘트 가늘어짐에서 관찰되는 $r(t) \sim \exp(-t/3\lambda)$ 감쇠 법칙과 다르며, 이는 워싱턴 제트 특유의 경계 조건과 축 방향 인장력에 기인한다.
• 안정성: 높은 탄성은 Rayleigh-Plateau 불안정성을 억제하여 더 굵고 안정적인 필라멘트를 형성하며, 핀치오프를 지연시키거나 방지한다.

5. 응력 및 속도장에 대한 수치적 통찰

시뮬레이션은 실험적으로 접근할 수 없는 내부 흐름장을 시각화하였다:

• 응력 분포: 탄성 응력($\tau_{zz}$)은 제트의 축 대칭 중심에서 최대값을 가지며 인터페이스를 향해 감소한다. 이는 높은 축 방향 속도와 중심부의 큰 고분자 신장에 대응한다.
• 속도 억제: 탄성 유체에서는 공동 확장 속도가 순수 물의 약 절반 수준이다. 이는 반대 방향의 탄성 응력 때문이다. 결과적으로 공기 유입(air entrainment)이 크게 감소하며, 뉴턴 케이스와 달리 수축 중 공동 내부에서 와류(vortex)가 관찰되지 않는다.

의의 및 주장

저자들은 본 연구가 관성-모세관-탄성 상호작용이라는 미개척 패러다임의 액적-풀 상호작용 이해를 크게 진전시켰다고 주장한다.

• 기초 물리학: 본 연구는 점성 효과로부터 탄성 효과를 분리하여, 유체 탄성이 어떻게 공동 형태(사다리꼴 vs 반구형)를 근본적으로 변화시키고 제트 파쇄를 억제하는지 입증하였다.
• 정량적 통찰: 공동 반경 진화와 제트 가늘어짐에 대한 이론적 모델을 도출함으로써, 액적-풀 상호작용 연구를 통해 탄성 유체의 연신 이완 시간(extensional relaxation time)을 직접 추정할 수 있는 방법을 제공한다.
• 레짐 정의: 네 가지 형태학적 레짐의 식별은 페인트, 잉크, 생물학적 유체와 같은 비뉴턴 유체가 사용되는 응용 분야에 중요한 예측 프레임워크를 제공한다.

결론적으로, 본 논문은 운동 에너지가 탄성 에너지로 전환되는 것이 탄성 유체에서 크라운 파쇄를 억제하고, 공동 크기를 줄이며, 워싱턴 제트를 안정화하는 핵심 메커니즘임을 밝히고 있다.