진화 표면 위 Navier–Stokes 방정식의 완전 이산 유한요소법

초록

본 논문은 정상 속도가 미리 주어진 3차원 진화 표면 위에서의 비압축성 Navier–Stokes 방정식을 대상으로, 일반화된 Taylor‑Hood 유한요소쌍(Pₖᵤ–Pₖₚᵣ–Pₖ_λ)과 후진 오일러 시간스키마를 이용한 두 가지 완전 이산 스키마를 제안한다. 라그랑주 승수 λ를 통해 정상 속도 제약을 약하게 강제하고, 기하학적 정보 요구량에 따라 ‘정보‑집약형’과 ‘정보‑경량형’ 두 스키마를 구분한다. 각 스키마에 대해 안정성 및 최적 오차 추정(속도 L²ₐₕ‖·‖, 압력 L²‖·‖ 및 이중 공간 L²‖·‖×L²H⁻¹ₕ‖·‖)을 증명하고, 수치 실험으로 이론을 검증한다.

상세 분석

이 연구는 진화 표면 유체역학 분야에서 아직 충분히 다루어지지 않은, 정상 속도가 외부에서 강제되는 상황(eNSEL)을 대상으로 한다. 기존 연구들은 주로 표면 자체의 변형을 연동하거나, Eulerian‑CutFEM 접근을 사용해 시간‑연속성 문제를 회피했지만, 본 논문은 Lagrangian‑parametric ESFEM을 채택해 표면의 움직임을 정확히 추적한다. 핵심 아이디어는 두 개의 라그랑주 승수(압력 p와 정상 속도 제약 λ)를 동시에 도입하고, 이를 각각 Pₖₚᵣ와 Pₖ_λ 다항식 공간으로 근사함으로써 inf‑sup 안정성을 확보하는 것이다.

특히, k_λ의 선택에 따라 두 가지 스키마가 정의된다. k_λ = kᵤ인 경우, 라그랑주 승수 λ를 속도와 동일 차수로 근사해 ‘정보‑집약형(eNSW_dh)’과 ‘정보‑경량형(eNSW_ch)’ 모두에서 최적 속도 L²ₐₕ 오차를 얻을 수 있다. 여기서는 iso‑parametric 격자를 사용해 기하학적 근사오차를 최소화한다. 반면, k_λ = kᵤ−1인 경우에는 λ의 차수가 낮아지므로 기하학적 정보가 더 많이 필요하고, super‑parametric 격자를 도입해야만 최적 압력 L²‖·‖ 오차를 확보할 수 있다.

새롭게 도입된 ‘표면 Ritz‑Stokes 사영 R_h’는 연속 문제의 해를 이산 공간에 투사하면서 물질미분(∂∘)과 사영이 교환되지 않는 문제를 해결한다. 이 사영의 근사 특성을 정량화하기 위해, 물질미분을 포함한 H^{kᵤ+1} 정규성 가정 하에 서브‑옵티멀한 L²‖∂∘(u−R_h u)‖ 추정식을 유도한다. 이러한 결과는 이후 속도와 압력에 대한 최적 오차 증명에 핵심적인 역할을 한다.

압력 안정성 분석에서는 ‘Leray 시간 사영’ ˆu^{n−1}_h 를 도입한다. 이는 이전 시간 단계의 이산 속도 u^{n−1}_h 를 현재 표면에 투사해 약한 발산 자유(divergence‑free) 조건을 만족하도록 만든다. 이를 통해 두 시간 단계 사이에 발생할 수 있는 발산 불일치를 억제하고, L²‖·‖×L²H⁻¹ₕ‖·‖ 형태의 압력 오차 추정에 필요한 이산 역 Stokes 연산자를 정의한다. 또한, λ에 대한 이산 inf‑sup 조건을 L²H⁻¹ₕ‖·‖(또는 L²‖·‖) 공간에서 확보함으로써, λ 자체의 오차는 약한 이중 공간에서만 제어되지만 압력 p는 강한 L²‖·‖ 오차를 얻는다.

수치 실험에서는 두 스키마 모두에 대해 수렴률을 확인하고, 기존의 penalty 방식과 비교해 라그랑주 승수 기반 접근이 더 높은 정확도와 안정성을 제공함을 보여준다. 특히, 복잡한 표면 변형과 비정상적인 외부 힘이 가해지는 테스트 케이스에서도 이론적으로 예측한 수렴 차수를 유지한다는 점이 강조된다.

전체적으로, 본 논문은 진화 표면 위 Navier–Stokes 방정식에 대한 완전 이산 ESFEM 해석을 최초로 제공하며, 라그랑주 승수 기반 정상 속도 제약, 새로운 Ritz‑Stokes 사영, Leray 시간 사영 등 여러 혁신적인 수학적 도구를 결합해 최적 오차와 안정성을 동시에 달성한다는 점에서 학술적·실용적 의의가 크다.