[논문 리뷰] Field conserving adaptive mesh refinement (AMR) scheme on massively parallel adaptive octree meshes

Kumar Saurabh, Makrand A. Khanwale외 3명|arXiv (Cornell University)|2026. 02. 08.

Solidification and crystal growth phenomena인용 수 0

한 줄 요약

논문은 연속 Galerkin 이산화에서 병렬 옥트리 기반 AMR을 위한 필드 보존 코어싱 연산자를 도입하여 코아상에서의 이산적 글로벌 보존을 보장하고, CH–유형 시스템에서 이를 평가한다.

ABSTRACT

Adaptive mesh refinement (AMR) is widely used to efficiently resolve localized features in time-dependent partial differential equations (PDEs) by selectively refining and coarsening the mesh. However, in long-horizon simulations, repeated intergrid interpolations can introduce systematic drift in conserved quantities, especially for variational discretizations with continuous basis functions. While interpolation from parent-to-child during refinement in continuous Galerkin (CG) discretizations is naturally conservative, the standard injection-based child-to-parent coarsening interpolation is generally not. We propose a simple, scalable field-conserving coarsening operator for parallel, octree-based AMR. The method enforces discrete global conservation during coarsening by first computing field conserving coarse-element values at quadrature points and then recovering coarse nodal degrees of freedom via an $L^2$ projection (mass-matrix solve), which simultaneously controls the $L_2$ error. We evaluate the approach on mass-conserving phase-field models, including the Cahn--Hilliard and Cahn--Hilliard--Navier--Stokes systems, and compare against injection in terms of conservation error, solution quality, and computational cost.

연구 동기 및 목표

• AMR을 사용한 장기 다물리 시뮬레이션에서 질량 보존을 동기화.
• AMR 주기 동안 이산적 글로벌 보존을 보장하는 보존적 코아싱 연산자 개발.
• 오크트리 메쉬에서 CG 이산화와 호환되는 확장 가능한 알고리즘 제공.
• CH 및 CHNS와 같은 상-필드 모델에서 보존 및 정확성 개선 시연.

제안 방법

• 먼저 대 coarse-element 적분점에서 보존을 강제한 다음, 질량-행렬 풀을 통한 L2 투영으로 대조점의 노드 DOF를 회복하는 필드 보존 코아싱 방법을 사용한다.
• 정제(부모에서 자식으로의) 보존성을 필요로 하는 반면, 코아싱(자식에서 부모로의)은 각 대 요소에 걸친 글로벌 적분을 보존하는 대 coarse 격자에 대한 국소 L2 투영을 사용한다.
• 로컬 L2 투영을 통해 미세-요소 가우스 포인트 값을 대-요소 가우스 포인트 값으로 매핑하는 제한 연산자를 구성하여 표준 가우스-루지 근처에서 질량 행렬이 대각선형(또는 쉽게 가역)으로 얻어지도록 한다.
• 대-포인트 데이터를 대-유한 요소 공간으로 매핑하는 글로벌 L2 투영(M G = b)을 풀어 대-노드 값을 조립하고 질량 보존을 보장한다.
• 1D 제한 연산자의 텐서 곱을 사용하여 임의의 다항 차수(Pp/Qp)와 다차원에 확장하고, 2D/3D에 대해 크로네커-구조로 확장한다.

Figure 1 : Refinement: Figure representing the refinement case. The red marker represents the extra points that are added on the refined mesh. The value of $\bigg[\int_{\Omega}g(x)\;d\Omega\bigg]_{\mathcal{M}_{O}}$ = $\bigg[\int_{\Omega}g(x)\;d\Omega\bigg]_{\mathcal{M}_{R}}$ = 10.6284 , where $\math

실험 결과

연구 질문

• RQ1 CG의 옥트리 AMR에서 축소(coarsening)를 정확도를 해치지 않으면서 질량을 보존하게 만들 수 있는가?
• RQ2 제안된 보존적 코아싱 방식이 주입(injection) 기반 코아싱에 비해 질량 드리프트를 줄이고 장시간 적분에서 보존 불변량을 유지하는가?
• RQ3 상-필드 형태의 PDE에서 보존적 코아싱이 수렴 속도와 계산 비용에 미치는 영향은 무엇인가?
• RQ4 단일 레벨 AMR 업데이트를 가진 대규모 병렬 고해상도 시뮬레이션에 이 방법이 확장 가능한가?

주요 결과

• 보존적 코아싱은 AMR 사이클 전반에 걸쳐 질량의 전역 적분을 수치적 정밀도로 보존한다.
• MMS 확산 테스트에서 보존적 방식은 선형(≈2) 및 이차(≈3)에서 최적 수렴 속도를 유지하고, 선형의 경우 주입 대비 약간 더 작은 오차를 보인다.
• 질량 드리프트 테스트에서 보존적 방식은 매 코아싱 단계에서 드리프트를 없애고, 주입은 더 작은 메쉬에서만 개선되는 비제로 드리프트를 보인다.
• 보존적 코아싱 하에서 이차 및 선형 기저가 비슷하거나 향상된 정확도를 보이며, 오차의 크기에 큰 차이는 없다.
• CH 및 CHNS-유형 모델에서 이 방법을 시연하여 장시간 상-필드 시뮬레이션에서 질량 보존의 이점을 강조한다.

Figure 2 : Coarsening by Injection : Figure demonstrating the coarsening procedure by injection. The gray shaded region shows the area corresponding to the difference between two mesh. $\bigg[\int_{\Omega}g(x)\;d\Omega\bigg]_{\mathcal{M}_{O}}$ = 10.6284, whereas $\bigg[\int_{\Omega}g(x)\;d\Omega\big