벡터 문제를 위한 사전·사후 오차 식과 볼록 이중성의 새로운 적용

초록

본 논문은 볼록 이중성을 이용해 비선형·비부드러운 스칼라 문제에서 성공적으로 활용된 사후 오차 식과 식별자를 벡터 문제, 특히 비압축성 Stokes와 Navier‑Lamé 방정식에 확장한다. 비동질 혼합 경계조건과 에너지 공간의 위상 이중에 속하는 하중까지 포함하고, Crouzeix‑Raviart와 Raviart‑Thomas 요소의 호환성을 활용해 최소 정규성 가정 하에 데이터 진동 항 없이 준최적 수렴을 증명한다.

상세 분석

이 연구는 기존의 Prager‑Synge 식을 볼록 이중성(fr convex duality) 프레임워크 안에서 일반화함으로써, 에너지 차이와 쌍대 차이가 동일한 형태의 오류 측정값을 제공한다는 점을 강조한다. 특히, 스칼라 Poisson 문제에서 사용되던 강한 이중성 관계 I(u)=D(q)를 벡터 공간으로 옮겨, 비압축성 Stokes 문제에서는 속도 u∈H₀¹(Ω)ᵈ와 응력 텐서 T∈F+H(div;Ω) 사이에 ν²‖∇v−∇u‖²+½ν‖dev τ−dev T‖²=ν²‖∇v−(1/ν)dev τ‖²와 같은 식을 도출한다. 여기서 dev는 편향 텐서의 편차 부분을 의미한다. 이 식은 임의의 시험 함수 v와 τ가 각각 연속성(div v=0, div(τ−F)=−f) 조건을 만족하면 정확한 등식이 되며, 이를 통해 사후 오차를 직접 계산할 수 있다.

Navier‑Lamé 문제에 대해서는 대칭성 제약을 완화한 확장된 쌍대 에너지 D를 정의하고, C¹/₂와 C⁻¹/₂ 연산자를 이용해 ‖C¹/₂ε(v−u)‖²+‖C⁻¹/₂(τ−σ)‖²≈‖C¹/₂(ε(v)−C⁻¹τ)‖² 형태의 등가 관계를 얻는다. 이때 σ는 실제 응력 텐서이며, 등가 관계의 상수는 지역 Korn 상수에만 의존한다는 점이 중요한데, 이는 비대칭 근사에 대한 허용 오차를 명시적으로 제어할 수 있음을 의미한다.

이론적 결과를 디스크리트 수준으로 옮기기 위해 저자는 Crouzeix‑Raviart 비정규 요소와 Raviart‑Thomas 저차원 요소를 각각 원시 변수와 쌍대 변수의 근사에 사용한다. 중요한 기술적 요점은 이산화된 원시 에너지 I_h와 쌍대 에너지 D_h 사이에 강한 이중성 I_h(u_h)=D_h(q_h)를 유지하도록, 원시 변수에 대한 요소별 그라디언트 ∇_h와 정규화 연산자 Π_h를 도입한 것이다. 이 구조 덕분에 이산 Prager‑Synge 식 ½‖∇_h u_h−∇_h v_h‖²+½‖Π_h q_h−Π_h r_h‖²=½‖Π_h r_h−∇_h v_h‖²를 얻을 수 있다.

이 식을 활용해 저자는 다음과 같은 준최적성 결과를 증명한다.
‖∇h u_h−∇u‖²+‖Π_h q_h−q‖² ≈ inf{v_h∈V_h^D}‖∇h v_h−∇u‖² + inf{r_h∈Σ_h, div r_h=−Π_h f}‖Π_h r_h−q‖².
여기서 상수는 메쉬의 형태 규칙성에 독립적이며, 데이터 진동 항이 전혀 등장하지 않는다. 이는 기존 결과(예: Verfürth 1996)의 주요 제한을 완전히 제거한 것으로, 특히 최소 정규성(u∈H₀¹)만을 가정해도 동일한 수렴 속도를 확보한다는 점에서 혁신적이다.

Stokes 문제에 대해서는 추가적으로 dev‑div 부등식(‖dev τ‖≤C‖div τ‖)을 이용해 비발산 시험 함수 v가 아닌 일반 v에 대해서도 등가 관계를 얻으며, 상수 C는 도메인에만 의존한다. Navier‑Lamé 문제에서는 Korn 상수와 함께 대칭성 위반을 정량화하는 추가 항이 등장하지만, 이는 요소 간 점프(s_h) 형태로 자연스럽게 포함된다.

마지막으로 저자는 수치 실험을 통해 이론적 수렴률과 오차 추정기의 효율성을 검증한다. 특히, 스트레스 재구성 공식 σ_h^* = C ε_h(u_h)+∇_h r_h−(2/d) sym(f_h⊗(id−Π_h id))를 사용해 비대칭 근사에도 불구하고 실제 응력과의 L² 오차가 기대값에 부합함을 보여준다. 전체적으로, 볼록 이중성을 기반으로 한 오류 식이 벡터 문제에 자연스럽게 확장될 수 있음을 입증하고, 최소 정규성 가정 하에서도 데이터 진동 없이 최적 수렴을 달성할 수 있음을 증명한다.