고점성 액체 시트의 후퇴 역학: 얇은 막 모델과 자가유사 팁 흐름의 통합 해석

초록

본 연구는 Oh ≫ 1·l₀/h₀ ≫ 1인 한정된 평면 액체 시트의 후퇴를 1차원 얇은 막 모델과 팁 영역의 자가유사 Stokes 흐름으로 분리하고, 두 영역을 매칭해 유효 경계조건을 도출한다. 결과적으로 두 변수(두께·속도)를 하나의 열방정식으로 축소시켜 단일 무차원 파라미터 ℒ = l₀/(4h₀ Oh) 로 기술한다. 초기에는 후퇴 속도가 T^{1/2} 로 증가하고, ℒ≫1인 경우 Taylor‑Culick 속도에 접근한다. T ≈ ℒ에서 급격히 감속해 T^{-2} 로 늦은 시간 해석에 이른다.

상세 분석

본 논문은 고점성(Oh ≫ 1)·긴(l₀/h₀ ≫ 1) 액체 시트의 후퇴 현상을 두 개의 비동질 영역으로 나누어 분석한다. 첫 번째는 전역적인 얇은 막 영역으로, 전형적인 질량·운동 방정식(연속식과 1차원 모멘텀식)을 적용한다. 여기서는 표면 장력이 직접적인 항으로 등장하지 않으며, 대신 팁 영역에서 전달되는 경계조건을 통해 구동된다. 두 번째는 자유 팁 근처의 O(h₀) 규모 영역으로, 전형적인 2차원 Stokes 흐름을 가정하고 자가유사 해를 구한다. 매칭 과정에서 얻은 유효 경계조건 4 μ h ∂v/∂x = −γ는 팁 영역의 점성·표면 장력 균형을 의미한다. 이 조건은 Oh ≫ 1, 즉 관성 항이 무시될 정도로 점성이 지배적인 경우에만 유효함을 논증한다.

스케일링 분석을 통해 얇은 막 영역에서 점성 응력이 표면 장력 응력보다 우세함을 확인하고(γ h² l⁻³ ≪ γ l*⁻¹), 따라서 모멘텀식에서 표면 장력 항을 제거한다. 그러나 관성 항은 무시할 수 없으며, 점성·관성 항의 균형으로부터 길이·시간 스케일 l_I = 4 h₀ Oh, t_I = 4 μ h₀/γ 를 도출한다. 이 스케일은 Taylor‑Culick 속도 u_TC = γ/(ρ h₀)와 직접 연결된다.

핵심적인 수학적 발견은 얇은 막 방정식이 보존량 Q = Hv − ∂H/∂X 를 갖는다는 점이다. 이를 이용해 속도와 두께를 선형적으로 연결하면, 원래의 연속·운동 방정식은 단일 열방정식 ∂H/∂T = ∂²H/∂X² 로 축소된다(시간 의존 경계조건 포함). 동일한 보존량을 이용하면 Burgers 방정식 형태로도 기술 가능하다.

축소된 문제는 무차원 파라미터 ℒ = l₀/(4 h₀ Oh) 만을 남기며, ℒ가 크면 초기 후퇴 속도가 Taylor‑Culick 속도에 근접하고, ℒ가 작으면 전체 시트 길이에 의해 제한된다. 수치 해석은 열방정식 풀이와 기존 전산 유체역학(N‑S) 시뮬레이션, 그리고 이전 얇은 막 모델(예: Brenner‑Gueyffier, Savva‑Bush, Deka‑Pierson) 결과와 정량적으로 일치함을 보여준다.

시간에 따른 속도 스케일은 세 단계로 구분된다. (1) 초기 확산‑주도 단계: u ∝ T^{1/2}. (2) ℒ≫1인 경우의 중간 단계: u ≈ u_TC (Taylor‑Culick). (3) 후기 단계: T ≈ ℒ에서 급격히 감속해 u ∝ T^{-2} 로 감소한다. 특히 ℒ≫1일 때는 속도가 거의 일정하게 유지되다가, 시트 끝에 도달하면서 급격히 감소하는 현상이 관찰된다. 이러한 전이 현상은 경계조건(고정된 끝에서 v=0, ∂H/∂X=0)과 보존량에 의해 자연스럽게 발생한다.

본 연구는 기존 얇은 막 모델이 팁 근처의 급격한 기울기와 곡률을 제대로 다루지 못한다는 한계를 극복하고, 점성·관성·표면 장력의 역할을 명확히 구분함으로써 고점성 시트 후퇴 현상의 전반적인 물리적 메커니즘을 통합적으로 제시한다.