Dune‑Fem에서 hp 적응형 비연속 유한요소 구현

초록

본 논문은 Dune‑Fem 라이브러리에 hp‑적응형 비연속 Galerkin 방법을 도입하기 위한 일반화된 알고리즘과 데이터 구조를 제시한다. 특히 격자 적응 과정에서 사용자 정의 데이터의 제한·보강(prolongation) 작업을 자동화하고, 기존 h‑적응 코드와의 호환성을 유지하면서 p‑적응과 hp‑적응을 지원한다. 구현된 dune‑fem‑hpdg 모듈을 이용해 내부 벌점법을 적용한 elliptic 문제의 hp‑적응 시연 결과를 제시한다.

상세 분석

이 논문은 Dune‑Fem의 기존 h‑적응 프레임워크를 확장하여 p‑적응과 hp‑적응을 동시에 처리할 수 있는 구조를 설계한다. 핵심은 전역 자유도(DOF) 집합 D(m)와 전역 DOF 매핑 μ(m)를 정의하고, 격자 변형 단계마다 새로운 DOF를 삽입·재배열·삭제하는 3단계 알고리즘(정의 3.3)을 제시한 점이다. 단계 1에서는 기존 DOF와 새로 생성된 DOF를 합쳐 중간 집합 D(m+½)와 연속 매핑 μ(m+½)를 만든다. 단계 2에서는 사용자 데이터(예: 해의 계수 벡터)를 임시로 크기를 늘린 뒤, 새 요소에 대해서는 지역 L²‑투영 또는 사용자가 지정한 다른 투영 연산자를 적용해 값을 초기화한다. 특히 p‑적응 시에는 같은 물리적 요소에 대해 다항 차수가 바뀔 때, 기존 DOF와 새로운 DOF 사이의 관계를 정확히 매핑해 주어야 하는데, 이를 전역 매핑 μ(m+½)와 지역 매핑 μ_E를 통해 일관되게 처리한다. 단계 3에서는 사용되지 않는 ‘구멍’(hole) 인덱스를 정리하고, 남은 DOF를 연속적인 인덱스로 재배열해 최종 매핑 μ(m+1)을 만든다. 이 과정은 기존 Dune‑Fem의 적응 단계와 동일한 인터페이스를 유지하면서도, p‑변경에 따른 데이터 손실을 방지한다.

논문은 또한 이 알고리즘을 구현한 dune‑fem‑hpdg 모듈의 설계 원칙을 설명한다. 모듈은 DiscreteFunctionSpace, LocalBasisFunctionSet, 그리고 DOFMapping 클래스를 템플릿화하여 사용자가 任意의 로컬 다항 차수와 형태함수를 정의할 수 있게 한다. 기존 Dune‑Fem 2.4와 호환되도록, 격자 적응 시에 필요한 영속 ID 매핑(persistent id)와 연관 컨테이너를 그대로 활용한다. 메모리 사용량은 단계 2에서 임시 벡터를 할당하더라도 전체 복잡도가 O(|G(m)∪G(m+1)|)에 머무르며, 실제 구현에서는 메모리 오버헤드가 최소화되도록 스마트 포인터와 이동语义를 활용한다.

마지막으로, hp‑적응 interior‑penalty DG 방법을 elliptic 문제에 적용한 실험을 제시한다. 적응 지표로는 에러 추정량과 다항 차수 선택 규칙을 사용했으며, h‑정밀화와 p‑증강을 적절히 조합해 자유도 대비 오류 감소율이 크게 향상됨을 보였다. 이 결과는 제안된 데이터 전이 메커니즘이 복합 적응 환경에서도 정확히 동작함을 실증한다.