모서리 모세관 상승의 보편성: t¹⁄³ 스케일링과 옴스거 원리 기반 해석

초록

옴스거 원리를 이용해 곡선형 모서리( y = c xⁿ, n>1)에서의 모세관 상승 동역학을 기술하는 편미분 방정식을 도출하고, 수치·자기유사 해를 통해 전단 높이가 t¹⁄³ 로 증가함을 확인하였다. 전단 전진 계수 C는 n에 대해 약 10 % 정도만 변해, 모서리 형상에 관계없이 보편적인 거동을 보인다.

상세 분석

본 논문은 두 곡면이 y = ±E(x)/2 로 정의되는 파워‑law 코너에서, 습윤액이 중력 하에 상승하는 과정을 옴스거 원리(Rayleighian 최소화)로 정형화한다. 자유에너지 변화율 ˙F와 점성 손실 Φ를 각각 (2.6)·(2.14)식으로 구한 뒤, 변분 조건 δR/δQ=0을 적용해 유량 Q(z,t)를 얻는다. 보존식(2.3)과 결합해 meniscus 프로파일 G(z,t)의 진화 방정식(2.17)을 도출하고, 이를 파워‑law 형태 E(x)=c xⁿ (n>1)으로 특수화한다. 핵심 파라미터 B(G), A′(G), ∂(L′/A′)/∂z 를 (3.2)~(3.4)식으로 전개해, 최종 비선형 편미분식(3.6)을 얻는다.

스케일링 변수를 H_c와 t_c 로 정의해 차원 없는 형태(3.8)·(3.9)로 정리하면, 자가유사 해 G̃(ẑ, t̃)=F(χ) t̃^{α}, χ=ẑ t̃^{β} 를 가정했을 때 α=−1/(3n), β=−1/3 이 도출된다. 이는 전단 위치 Z_m ∝ t^{1/3} 임을 의미한다. ODE (3.15)와 경계조건(3.16, 3.23)을 수치적으로 풀어 χ₀ 값을 구하면, 전단 전진 계수 C = χ₀ (n² cos²θ/3)^{1/3} 가 된다. χ₀는 n에 대해 약 1~1.1 정도로 변동해, n=1(선형 코너)부터 n=5까지 전단 속도는 10 % 이내 차이만 보인다.

선형 코너(n=1)와 2차 코너(n=2)에 대해 각각 (4.2), (4.10)식의 차원 없는 형태를 수치해석했으며, 시간에 따른 meniscus 형태와 Z_m(t) 곡선이 실험 데이터( Higuera 2008, Ponomarenko 2011)와 일치함을 확인했다. 특히, 장기(large‑t) 영역에서 t^{1/3} 스케일링과 앞서 제시된 C값이 정확히 재현된다.

이러한 결과는 기존 Lucas‑Washburn t^{1/2} 스케일링(중력 무시)과는 달리, 중력과 점성 손실이 동시에 작용하는 경우에도 코너 형상에 관계없이 동일한 지배 방정식과 스케일링을 갖는다는 강력한 보편성을 제시한다. 옴스거 원리를 통한 변분 접근법은 경계 움직임이 복잡한 유동 문제에 적용 가능성을 넓히며, 실험적 관측과 이론적 예측을 일관되게 연결한다는 점에서 학문적·공학적 의의가 크다.