공동 흐름 장치에서의 흐름 안정성 재검토

초록

공동 흐름(co‑flow) 마이크로채널에서 발생하는 액체 제트는 전통적인 전역 선형 안정성 분석으로는 드리핑 전이 조건을 정확히 예측하지 못한다. 실험과 전이(Transient) 시뮬레이션 결과, 제트는 콘형 매니스크스가 안정된 상태에서도 먼저 파괴되며, 초기 섭동의 형태에 따라 파괴 길이가 크게 달라진다. 이는 감쇠 고유모드들의 비정상적(Non‑modal) 중첩이 단기적으로 에너지를 증폭시켜 제트를 붕괴시키기 때문이다.

상세 분석

본 논문은 공동 흐름 장치에서 발생하는 액체 제트의 안정성을 전역 선형 안정성(Global Linear Stability)과 전이(Transient) 시뮬레이션을 결합해 재검토하였다. 먼저, 축대칭 Navier‑Stokes 방정식을 비차원화하고, 내부·외부 유체의 밀도·점도비(ρ=0.962, μ=36.1)와 Ohnesorge 수(Oh_i=0.0264)를 고정한 뒤, 반경비 R=R_o/R_i=6.67이라는 실험적 파라미터 집합을 사용하였다. 전역 선형 안정성 분석에서는 정상 흐름(base flow)을 가정하고, 선형화된 연산자에 대한 고유값 문제를 Chebyshev 스펙트럼 방법으로 풀어 고유주파수 ω=ω_r+iω_i를 얻었다. 이때 가장 큰 성장률을 가진 지배 고유모드의 실수부가 양(ω_i>0)이면 불안정, 음이면 안정으로 판단한다.

분석 결과, 매니스크스(콘형 팁) 근처에서는 어떠한 파라미터 조합에서도 지배 고유모드의 진폭 δF(z)가 눈에 띄게 증가하지 않았다. 즉, 전역 선형 분석상 매니스크스는 항상 안정적인 것으로 나타났다. 반면 실험에서는 Q (내부/외부 유량비) 가 일정 임계값 Q_min≈1.7×10⁻⁴ 이하로 감소하면 매니스크스가 진동하면서 다분산 드리핑이 발생한다는 점을 확인하였다.

전이 시뮬레이션에서는 초기 섭동을 자유면(F(z,0)−F₀(z)=β e^{−α|z−z₀|}) 형태로 가정하고, 선형화된 방정식을 직접 적분하였다. β=10⁻³, α=10이라는 작은 변형을 주었을 때, 감쇠 고유모드들의 중첩이 단기적으로 에너지 증폭을 일으켜 제트가 매니스크스보다 먼저 파괴되는 현상이 관찰되었다. 특히, 초기 섭동의 위치(z₀)와 형태에 따라 파괴 길이가 크게 달라졌으며, 이는 전통적인 고유모드 분석만으로는 설명되지 않는다.

이러한 비정상 성장(Non‑modal growth)은 감쇠 고유모드들의 비직교성에 기인한다. 감쇠 모드들이 동시에 활성화되면, 순간적인 에너지 증폭이 발생해 비선형 효과가 나타나기 전에 제트가 파열된다. 따라서 전역 선형 안정성 분석은 “지배 고유모드가 안정”이라는 결론을 내리지만, 실제 흐름은 초기 섭동에 민감한 비정상적인 동역학에 의해 제어된다.

결론적으로, 공동 흐름 장치에서의 제트‑드리핑 전이는 매니스크스의 절대 불안정성보다는 제트 자체의 전도성(convective) 불안정성과 초기 섭동에 의한 비정상 성장에 의해 결정된다. 따라서 설계 단계에서 단순히 전역 선형 안정성 결과만을 의존하는 것은 위험하며, 전이 분석이나 비정상 모드 해석을 병행해야 한다.