압력 강인 재구성 연산자를 이용한 Stokes 최적 제어 문제의 사후오차 추정 개선

초록

본 논문은 발산이 없는 재구성 연산자를 도입한 압력-강인 유한요소법을 Stokes 방정식에 기반한 분산형 최적 제어 문제에 적용하고, 압력에 독립적인 새로운 잔차 기반 사후오차 추정기를 제안한다. 제안된 추정기는 전역적인 신뢰성(upper bound)과 효율성(lower bound)을 보이며, 속도와 압력 오차를 분리해 속도 오차가 압력에 영향을 받지 않도록 설계되었다. 수치 실험을 통해 이론적 결과와 비균질 경계조건에서도 적용 가능함을 확인하였다.

상세 분석

이 연구는 Stokes 최적 제어 문제에서 기존 사후오차 추정기가 압력 항에 의존해 저점성(ν→0) 상황에서 속도 오차가 과도하게 커지는 문제점을 정확히 짚어낸다. 이를 해결하기 위해 저자들은 세 번째 접근법인 ‘발산이 없는 재구성 연산자(R_h)’를 채택한다. R_h는 V_h(속도 유한요소 공간)에서 H(div)‑정합성을 만족하는 W_h(RT1) 공간으로 사상하면서 (i) 발산이 0, (ii) 저차 다항식에 대해 L^2‑직교, (iii) 근사 오차가 h‑스케일에 비례하도록 설계된다. 이러한 연산자를 우변에 적용함으로써 원래 연립식 a(y_h,v_h)−b(v_h,p_h)=(f,R_h v_h) 로 변형하고, 이는 연속 문제의 ‘압력-불변성(invariance)’을 이산화 수준에서도 보존한다.

논문은 먼저 기존 비압력‑강인 추정기(정리 3.1)와 그 한계를 서술한다. 기존 추정식은 볼륨 항에 ∇p_h가 포함돼 ν⁻¹·‖p‖가 속도 오차에 직접 기여한다. 이어서 압력‑강인 재구성 연산자를 도입한 새로운 이산화 방식을 제시하고, 이를 기반으로 잔차 기반 추정기 ˆε를 재정의한다. 핵심은 (f−∇p_h+νΔy_h) 대신 (f+ν div_h(∇y_h)) 형태로 바꾸어 압력 항을 완전히 소거하는 것이다. 또한, 상태와 대입 변수 각각에 대해 압력‑제거된 잔차를 별도로 구성해, 전체 추정기가

‖u−u_h‖ + ‖∇(y−y_h)‖ + ‖∇(z−z_h)‖ ≤ C·ˆε + OSC

와 같은 형태의 상한을 만족함을 증명한다. 여기서 OSC는 데이터 진동(oscillation) 항으로, 고차 근사에 의한 오차를 정량화한다. 하한(효율성) 증명에서는 거품 함수와 역추정 기법을 활용해 ˆε가 실제 오차를 과소평가하지 않음을 보인다.

수치 실험에서는 2차 및 3차 Lagrange 요소, RT1 기반 재구성 연산자를 사용해 다양한 ν값(10⁻¹~10⁻⁴)과 비균질 경계조건을 시험한다. 결과는 압력‑강인 추정기가 ν에 무관하게 수렴률을 유지하고, 기존 추정기에 비해 오류 상수가 현저히 작아짐을 보여준다. 또한, 적응적 메쉬 정제(adaptive refinement) 과정에서 제안된 추정기가 효과적인 마커 역할을 수행함을 확인한다.

이 논문의 주요 기여는 다음과 같다. (1) 압력‑강인 재구성 연산자를 활용한 Stokes 최적 제어 문제의 이산화 체계 구축, (2) 압력에 독립적인 잔차 기반 사후오차 추정기 설계 및 전역 신뢰·효율성 증명, (3) 다양한 실험을 통한 이론 검증 및 비균질 경계조건에 대한 확장성 입증. 이러한 결과는 저점성 유동 시뮬레이션이나 압력‑강인 요구가 높은 공학·과학 응용 분야에서 보다 안정적이고 정확한 적응형 유한요소 해석을 가능하게 할 것으로 기대된다.