임베디드‑하이브리드 불연속 갈루아 방법을 이용한 Stokes 방정식 고성능 해법

초록

본 논문은 Stokes 방정식에 대해 임베디드‑하이브리드(EDG‑HDG) 불연속 갈루아(DG) 방법을 제안한다. 기존 HDG 방법의 압력 강건성·H(div) 연속성을 유지하면서, EDG 방식처럼 연속적인 면 속도 자유도를 사용해 전역 자유도를 크게 감소시킨다. 이론적 안정성·오차 분석을 수행하고, 압력에 독립적인 속도 오차 추정(pressure‑robustness)을 증명한다. 수치 실험을 통해 수렴율을 확인하고, 전처리된 Krylov 솔버에서 EDG‑HDG가 HDG 대비 반복 횟수·시간 모두에서 우수함을 보인다.

상세 분석

논문은 먼저 기존 하이브리드 불연속 갈루아(HDG)와 임베디드 불연속 갈루아(EDG)의 장단점을 정리한다. HDG는 면 압력 라그랑주 승수를 불연속적으로 배치해 H(div)‑연속성을 확보하고, 압력 강건성을 얻지만 면 자유도가 많아 전역 시스템이 커진다. 반면 EDG는 면 속도와 압력을 연속적으로 사용해 자유도를 연속 Galerkin 수준으로 낮추지만, H(div) 연속성이 사라져 압력 강건성이 손실된다.
이 두 접근법을 결합한 EDG‑HDG는 면 속도는 연속, 면 압력은 불연속으로 두어 전역 자유도를 크게 줄이면서도 속도에 대한 H(div) 연속성을 유지한다. 핵심 아이디어는 면 압력 라그랑주 승수를 HDG와 동일하게 두어 정상 성분의 연속성을 강제하고, 연속적인 면 속도 베이스를 통해 전역 행렬의 스펙트럼을 개선한다는 점이다.
수학적 분석에서는 다음을 입증한다.

1. 일관성: 정확 해가 충분히 정규화된 Sobolev 공간에 있으면, 제안된 변분 형식 B_h는 연속적인 선형 형태와 동일하게 동작한다.
2. 안정성: a_h(·,·)에 대한 coercivity와 bounded성은 기존 HDG 결과를 그대로 차용하고, b_1, b_2(압력‑속도 결합) 각각에 대해 Brezzi‑Douglas‑Marini(BDM) 보간 연산자를 이용해 inf‑sup 조건을 증명한다. 특히 b_1에 대해서는 ∇·(Π_BDM v)=q_h를 만족하는 보간을 구성해 압력 자유도에 대한 하한을 확보한다. b_2는 면 압력과 연속 면 속도 사이의 결합으로, 기존 전처리 논문