3차원 Quaddiv 문제를 위한 GradDiv 적합 가상 요소법

초록

본 논문은 3차원에서 발생하는 Quad‑div 연산자를 위한 새로운 변분 형식을 제시하고, 이를 기반으로 H(grad‑div) 적합 가상 요소법(VEM)을 다항형 격자 위에 고차 정확도로 구축한다. 세 가지 VEM 계열을 정의하고, 저차원 원소는 정점·면당 하나의 자유도만을 갖는다. 해석적 측면에서 보간 오차, 이산 이중형식의 안정성, 이산 문제의 잘 정의됨, 최적 수렴율을 엄밀히 증명하고, 수치 실험을 통해 이론을 검증한다.

상세 분석

이 연구는 기존의 Quad‑div 연산자에 대한 유한요소 접근이 제한적이었던 점을 극복하기 위해, 먼저 (∇ div)² 연산자를 위한 새로운 변분 문제를 정의한다. 핵심은 V₀(Ω)= {v∈L²(Ω) | div v∈H¹(Ω), v·n=0, div v=0 on ∂Ω} 라는 H(grad‑div) 적합 공간을 정확히 기술하고, 그 쌍대공간 V′(Ω)를 H^{‑2}(curl;Ω) 로 동등시킨 점이다. 이를 통해 curl‑free 소스 f∈V′(Ω)∩ker(curl) 를 허용하면서도, 라그랑주 승수 φ와 압력 p 를 도입한 혼합 변분식 (3.3)을 구성한다. 프리드리히스와 바우어‑스키–브레즈키 조건을 이용해 A와 B 이중형식이 각각 coercive와 inf‑sup 조건을 만족함을 증명함으로써, 연속 문제의 존재·유일성을 확보한다.

가상 요소법 설계 단계에서는 de Rham 복합체를 강화된 매끄러움과 함께 확장한 (1.5) 형태의 이산 복합체를 구축한다. 여기서 U₁(Ω)와 W_k(Ω) 를 H¹‑공간의 서브스페이스로 정의하고, Σ_{0,r}(Ω) 를 H(curl)‑적합 VEM으로, V_{r‑1,k+1}(Ω) 를 H(grad‑div)‑적합 VEM으로 만든다. r=k, k+1, k+2 세 경우에 대해 각각 다른 자유도 배치를 제공하며, 특히 r=k=1 일 때는 정점·면당 하나의 자유도만을 갖는 가장 단순한 원소가 얻어진다.

보간 연산자는 각 원소 내부에서 경계값 문제를 푸는 방식으로 정의되며, 연속 복합체와 이산 복합체 사이의 교환도표(commutative diagram)를 만족하도록 설계한다. 이를 통해 보간 오차 ‖v−I_h v‖{H(div)} ≤ C h^{k}‖v‖{H^{k+1}} 와 같은 최적 차수의 추정식을 얻는다.

이산 이중형식의 안정성은 요소 내부의 “stabilization term”을 적절히 선택함으로써 확보한다. 특히 Lemma 5.1 의 이산 프리드리히스 부등식은 H(grad‑div)‑공간에서의 norm equivalence을 제공하여, 연산자 A_h 가 coercive 하고 B_h 가 inf‑sup 조건을 만족함을 보인다. 따라서 이산 문제는 연속 문제와 동일한 해 존재·유일성을 가진다.

오차 분석에서는 연속 해 u∈V₀(Ω) 가 H^{s}(Ω) (s>½) 에 속함을 이용해, ‖u−u_h‖{V} ≤ C h^{k}‖u‖{H^{k+1}} 와 같은 최적 수렴률을 증명한다. 또한 φ와 p 에 대해서도 동일 차수의 수렴을 얻는다.

마지막으로, 다각형·다면체 격자에서의 구현 가능성을 강조하며, 정규 격자와 비정규 격자 모두에서 수치 실험을 수행한다. 실험 결과는 이론적 수렴률과 일치하고, 특히 저차원 원소가 정점·면당 하나의 자유도만으로도 높은 정확도를 보임을 확인한다. 전체적으로, 이 논문은 Quad‑div 연산자를 위한 H(grad‑div) 적합 VEM을 최초로 체계화하고, 수학적 엄밀성과 실용적 구현 양면에서 중요한 기여를 한다.