[논문 리뷰] Weight-adjusted discontinuous Galerkin methods: curvilinear meshes

Jesse Chan, Russell J. Hewett외 1명|arXiv (Cornell University)|2016. 08. 12.

Advanced Numerical Methods in Computational Mathematics참고 문헌 29인용 수 4

한 줄 요약

이 논문은 곡선 메esh에서 고차 정확도와 에너지 안정성을 유지하면서 저장 비용을 줄이기 위해 기하학적 가중치를 적용한 질량 행렬의 역행렬을 가중치 조정된 적분을 통해 근사하는 가중치 조정 불연속 갈레르킨(WADG) 방법을 제안한다. 이 방법은 최적의 수렴 속도와 GPU에서의 높은 효율성을 확보하며, 메쉬 변형에 대해서도 강건함을 입증하였다.

ABSTRACT

Traditional time-domain discontinuous Galerkin (DG) methods result in large storage costs at high orders of approximation due to the storage of dense elemental matrices. In this work, we propose a weight-adjusted DG (WADG) methods for curvilinear meshes which reduce storage costs while retaining energy stability. A priori error estimates show that high order accuracy is preserved under sufficient conditions on the mesh, which are illustrated through convergence tests with different sequences of meshes. Numerical and computational experiments verify the accuracy and performance of WADG for a model problem on curved domains.

연구 동기 및 목표

• 밀도 높은 요소 질량 행렬로 인해 기존의 불연속 갈레르킨(DG) 방법이 곡선 메쉬에서 높은 저장 비용을 유발하는 문제를 해결하기 위해.
• 곡선 영역에서 시간 도메인 파동 시뮬레이션에서 에너지 안정성과 고차 정확도를 유지하기 위해.
• 비용이 많이 드는 행렬 분해나 실시간 해법을 피하는 저저장량의 기존 곡선 DG 방법의 대안을 개발하기 위해.
• 메모리 대역폭과 저장 요구량을 최소화하여 GPU에서의 효율적 구현을 가능하게 하기 위해.
• 추가 계산 비용 없이 공간적으로 변화하는 파동 속도를 갖는 이종 매질로의 확장을 위해.

제안 방법

• WADG 방법은 기하학적으로 가중진 질량 행렬의 역행렬을 가중치 조정된 질량 행렬을 통해 근사함으로써 저장량을 대각형 행렬 수준으로 줄인다.
• L² 프로젝션을 근사하기 위해 가중치 조정된 적분점 기반의 가짜 프로젝션을 사용하며, 오차 한계는 자코비안의 변동성에 따라 달라진다.
• 체적 및 표면 적분을 분리하기 위해 강건-약한 공식화를 사용하여 GPU 커널 실행을 효율적으로 구현한다.
• 수치적 표면 플럭스는 적분점에서의 노드 값으로 계산되며, 업데이트 커널은 공간적으로 변화하는 파동 속도 $c^2$를 수정된 대각형 가중치를 통해 통합한다.
• 이 방법은 자코비안과 파동 속도의 변동을 질량 행렬 근사에 통합함으로써 DG 변분 공식과의 일致성을 유지한다.
• 표면 적분을 위해 분할 커널 전략을 사용하여 고차 요소에서 특히 메모리 대역폭 활용도를 향상시킨다.

실험 결과

연구 질문

• RQ1 저저장량 DG 방법은 공간적으로 변화하는 기하학적 요소를 갖는 곡선 메쉬에서 고차 정확도와 에너지 안정성을 유지할 수 있는가?
• RQ2 WADG 방법이 질량 행렬의 역행렬을 근사함으로써 이산 해의 조건수와 수렴 특성에 어떤 영향을 미치는가?
• RQ3 메쉬 곡률은 WADG 방법의 안정성과 정확도에 표준 DG 및 LSC-DG와 비교해 어떤 영향을 미치는가?
• RQ4 WADG 방법은 추가 계산 비용 없이 공간적으로 변화하는 파동 속도를 갖는 이종 매질을 효과적으로 처리할 수 있는가?
• RQ5 WADG 방법의 계산 성능는 다양한 다항식 차수에서 GPU 아키텍처에서 어떻게 스케일링되는가?

주요 결과

• WADG 방법은 곡선 영역에서 파동 전파 문제에 대해 최적의 $h$- 및 $p$-수렴 속도를 달성하며, 오차 한계는 자코비안 변동성에 명시적으로 의존한다.
• 수치 실험을 통해 WADG 방법은 심한 메쉬 곡률 조건에서도 안정성과 고차 정확도를 유지함을 확인하였으며, 메쉬 변형에 대한 강건성에서 표준 DG 및 LSC-DG를 뛰어넘었다.
• Nvidia GTX 980 GPU에서 $N=5$일 때 WADG 형식은 최대 1500 GFLOPS/s 및 200 GB/s 대역폭을 달성하였으며, 분할 표면 커널이 메모리 대역폭 활용도를 향상시켰다.
• 업데이트 커널의 대각형 가중치를 수정함으로써 공간적으로 변화하는 파동 속도를 갖는 이종 매질을 지원하며, 계산 비용 증가 없이도 가능하다.
• WADG 가짜 프로젝션 오차는 $L^2$ 노름에서 유계이며, 자코비안의 정규성과 적분 가중치 함수에 따라 명시적인 상수가 존재한다.
• 비약선형 요소 또는 특이 매핑을 갖는 전이 피라미드와 같은 경우, 수렴이 $h$- 및 $p$-재분할 조건에서 멈추므로 이 방법은 적용 불가능하다.