대칭·비대칭 부유 방울의 에너지 비교와 수치 해법

초록

본 논문은 부유 방울의 평형 상태를 기술하는 자유경계 문제를 뉴턴 방법과 체비셰프 스펙트럴 콜로케이션을 결합한 수치 알고리즘으로 해결한다. 9개의 물리 파라미터를 탐색하여 2차원·3차원에서 대칭·비대칭 구성의 잠재 에너지를 비교하고, 해의 비유일성 및 에너지 최소화 해의 다중성 현상을 제시한다.

상세 분석

이 연구는 세 개의 불혼합 유체가 제한된 용기 안에서 평형을 이루는 ‘부유 방울’ 문제를 수치적으로 분석한다. 먼저, 각 유체 사이의 계면에 대한 표면 장력 σij와 접촉각 γjk를 정의하고, 중력에 의한 체적 에너지와 습윤 에너지를 포함한 전체 에너지 함수 E를 구성한다. 변분 원리를 적용하면 Young‑Laplace 방정식 2H = κz + λ/σ 형태의 비선형 경계값 문제가 도출되며, 여기서 κ는 각 계면의 캡릴러 상수, λ는 부피 제약을 위한 라그랑주 승수이다.

저자들은 λ를 자유 경계 Γ에서 계면을 매칭함으로써 제거하고, 실제 물리 계면 u, v, w를 구한다. 이때 중심 대칭 방울은 반지름 r=0을 통과하는 형태로, 벽에 접촉하는 방울은 annular 형태로 모델링한다. 두 경우 모두 방정식은 아크 길이 s를 매개변수로 하는 1차 ODE 시스템( dr/ds = cosψ, dU/ds = sinψ, dψ/ds = κU−sinψ·r 등)으로 변환된다.

수치 해법은 Newton 반복을 통해 비선형성을 처리하고, 각 반복 단계에서 Chebyshev 스펙트럴 콜로케이션을 이용해 선형 시스템을 고정 정확도로 푼다. 14개의 Chebyshev 점을 사용해 10⁻¹⁰ 수준의 상대 오차를 달성했으며, 알고리즘은 기존의 shooting 방법보다 안정적이다.

파라미터 공간은 9차원(부피 V, 용기 반경 R, 두 유체 밀도 ρ1, ρ2, 세 표면 장력 σ01, σ02, σ12, 두 접촉각 γ10^p, γ20^p)으로 정의된다. 저자들은 이 공간을 체계적으로 샘플링하여, (1) 중심에 위치한 방울과 벽에 붙은 방울의 에너지 차이를 계산하고, (2) 동일 파라미터에서 여러 해가 존재함을 확인하며, (3) 경우에 따라 비대칭 구성이 대칭 구성보다 낮은 에너지를 갖는 대칭 파괴 현상을 발견한다. 특히, 일부 파라미터 조합에서는 두 구성이 정확히 같은 에너지를 가져 에너지 최소화 해가 비유일함을 강력히 시사한다.

이러한 결과는 기존 이론이 보장하는 대칭성(예: Treinen의 존재 정리)과는 달리, 실제 물리적 조건 하에서는 다중 최소화 해와 비대칭 안정 상태가 공존할 수 있음을 보여준다. 또한, 2차원 모델에서도 동일한 현상이 재현되어, 차원 축소가 본질적인 비유일성 메커니즘을 제거하지 않음을 확인한다.