H(curl) 기반 스토크스 방정식의 슬립 경계조건 처리와 최적 유한요소 설계

초록

본 논문은 스토크스 방정식을 H(curl) 공간에서 정의하고, Navier 슬립 경계조건을 로빈 형태로 변환하여 연속 문제의 잘 정의됨을 증명한다. 이후 H(curl)‑준수 유한요소를 이용한 비정합 Galerkin 방법을 제시하고, 안정성·오차 추정식을 얻는다. 수치 실험을 통해 이론적 수렴률이 확인된다.

상세 분석

이 연구는 기존의 H¹ 기반 스토크스 해석에서 벗어나, 속도를 1‑형식으로 해석함으로써 H(curl)‑준수 유한요소를 적용할 수 있는 새로운 프레임워크를 제시한다. 핵심 아이디어는 Navier 슬립 조건을 “u·n=0, ‑n×(∇×u)+2W(uₜ)=0” 형태의 로빈 경계조건으로 재표현하는 것이다. 여기서 W는 Weingarten map이며, 곡률에 따라 양·음·불정형이 될 수 있어 α=2W가 부호가 변할 가능성을 고려한다. 연속 변분식은 a(u,v)= (f,v) 로 정의되며, a는 eₐ(∇×·,∇×·)+⟨γ∥u,γ∥v⟩와 (α−I)γ∥u·γ∥v의 합으로 구성된다. eₐ는 H(curl)에서 코에르시브함을 보이고, α가 부호가 바뀌어도 k(·,·)는 컴팩트 연산으로 처리한다. 따라서 Lax‑Milgram과 Fredholm 대안을 이용해 해의 존재·유일성을 확보한다.

또한, X=H(curl)∩H₀(div 0) 공간을 정의하고, Hodge‑Helmholtz 분해 X=H₁⊕⊥Z를 이용해 Poincaré‑Steklov 부등식과 정규성 결과를 도출한다. 이때, 도메인이 Lipschitz 다면체이면 X⊂Hˢ(Ω) (s>½) 로 연속 삽입되며, 이는 경계 트레이스 γ∥의 L²‑제어에 필수적이다.

이산화 단계에서는 Nédélec 1차 요소와 Lagrange 압축 요소를 결합한 비정합 공간 V_h⊂H(curl) 를 사용한다. V_h는 연속 공간 X와 “가깝다”는 의미에서 연속성 추정과 보강 연산자를 도입해 A_h의 coercivity와 K_h의 컴팩트성을 보인다. 결과적으로, 이산 문제도 연속 문제와 동일한 Fredholm 구조를 가지며, a_h(u_h,v_h)=(f,v_h) 가 유일해를 갖는다. 오류 분석은 Céa‑type 부등식을 기반으로, ‖u−u_h‖{H(curl)} ≤ C h^k (‖u‖{H^{k+1}}+‖p‖_{H^{k}}) 형태의 최적 차수 수렴을 얻는다. 압력은 스터프된 라그랑주 승수를 통해 별도 추정이 가능하다.

수치 실험에서는 2D 원형 및 3D 구형 도메인, 그리고 곡률이 변하는 다면체에서 정확한 Weingarten map을 사용한 경우와 근사식(Neun­teufel et al.)을 사용한 경우를 비교한다. 모든 테스트에서 이론적 수렴률(1차 요소는 O(h), 2차 요소는 O(h²))이 관측되었으며, α가 음의 값을 가질 때도 안정적으로 수렴한다. 특히, 비볼록 도메인에서 기존 H¹ 기반 방법이 해를 복원하지 못하는 반면, 제안된 H(curl) 접근법은 물리적 해를 정확히 재현한다는 점이 강조된다.

이러한 결과는 MHD 시뮬레이션에서 헬리시티 보존, 압력‑강건성 확보, 그리고 복잡한 경계조건을 가진 다중 물리 문제에 대한 새로운 수치 해법의 가능성을 열어준다.