안정화 없는 가상요소법의 잔여 사후오차 추정

초록

본 논문은 2차원 포아송 방정식에 대해 안정화 항을 사용하지 않는 가상요소법(VEM)의 사후오차 추정량을 제시한다. 안정화 항이 없으므로 오차 측정값과 전통적인 잔여 기반 추정량 사이의 등가성을 엄밀히 증명하고, 다양한 다각형 메쉬와 확산계수 불연속성에 대한 수치 실험을 통해 추정량의 신뢰성을 확인한다.

상세 분석

이 연구는 기존 VEM에서 필수적으로 도입되는 안정화 연산자(stabilization operator)의 존재가 사후오차 분석에 미치는 영향을 근본적으로 재검토한다. 저자들은 먼저 다각형 메쉬에 대한 전형적인 규칙(별형성, 최소 면 길이 비율 등)을 가정하고, 각 원소 E 에 대해 고차 다항식 공간 P_k(E) 와 그에 대한 H¹ 직교 투영 Π_{∇,E}^k 를 정의한다. 이후, 가상요소 공간 V_{E}^{h,k} 을 “Δv∈P_k(E)”, “경계값은 P_k(e)와 일치” 등으로 구성하고, 자유도는 정점값, 내부 가우스‑롭탈 점값, 내부 모멘트 등으로 설정한다.

핵심은 원소별 그래디언트 투영 연산자 Π_{0,E}^{P∇} 를 도입해, 실제 그래디언트를 직접 계산하지 않고도 K_E Π_{0,E}^{P∇} u_h 만으로 이산 이중형 a_h(u_h,v_h) =∑E K_E(Π{0,E}^{P∇}u_h, Π_{0,E}^{P∇}v_h)E 을 정의한다는 점이다. 이때 안정화 항이 전혀 포함되지 않으므로, 연산자 Π{0,E}^{P∇} 의 안정성(정리 3.1)만으로도 연속 문제와 이산 문제 사이의 유일해 존재성을 확보한다.

사후오차 분석에서는 새로운 오차 측정값 |||u−u_h|||{K,1,Ω} 을 정의하고, 이를 전통적인 잔여 r_E = f_h+∇·(K_E Π{0,E}^{P∇}u_h) 와 면 흐름 j_e =