[논문 리뷰] Weight-adjusted discontinuous Galerkin methods: wave propagation in heterogeneous media

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

Electromagnetic Simulation and Numerical Methods참고 문헌 20인용 수 5

한 줄 요약

이 논문은 이질성 매질 내 파동 전파를 위한 가중치 조정 불연속 갤러르킨(WADG) 방법을 제안한다. 이 방법은 요소별 가중치 부여 질량 행렬을 등가의 가중치 조정 내적으로 대체하여 저장 비용 증가 없이 에너지 안정성과 고차 정확도를 달성한다. 이 방법은 계산 효율성을 유지하며 기존 DG 프레임워크와 호환되며, 최적 수렴성과 초수렴 보존 오차를 달성한다.

ABSTRACT

Time-domain discontinuous Galerkin (DG) methods for wave propagation require accounting for the inversion of dense elemental mass matrices, where each mass matrix is computed with respect to a parameter-weighted L2 inner product. In applications where the wavespeed varies spatially at a sub-element scale, these matrices are distinct over each element, necessitating additional storage. In this work, we propose a weight-adjusted DG (WADG) method which reduces storage costs by replacing the weighted L2 inner product with a weight-adjusted inner product. This equivalent inner product results in an energy stable method, but does not increase storage costs for locally varying weights. A-priori error estimates are derived, and numerical examples are given illustrating the application of this method to the acoustic wave equation with heterogeneous wavespeed.

연구 동기 및 목표

• 국소적으로 변화하는 파장 속도를 갖는 파동 전파 문제에서 기존 DG 방법의 높은 저장 비용 문제를 해결하기 위해.
• 요소별 질량 행렬을 저장하거나 분해할 필요 없이 에너지 안정성과 고차 정확도를 유지하기 위해.
• 고차 DG 스킴의 메모리 용량을 줄여 GPU 구현을 효율적으로 가능하게 하기 위해.
• 기존 DG 코드에 비침습적인 수정을 통해 계산 구조와 확장성을 유지하기 위해.
• 매끄러운 파장 속도에 대해 점근적 일致성과 초수렴 보존 오차를 보장하기 위해.

제안 방법

• 에너지 안정성을 유지하는 가중치 조정 내적을 사용하여 매개변수 가중 $L^2$ 내적을 대체한다.
• 차수 $2N+1$의 저차수 적분을 사용한 적분 기반 변분 형식을 통해 가중치 조정 질량 행렬을 계산한다.
• 공간적으로 변화하는 파장 속도를 내적 구조에 흡수하는 수정된 기저를 사용하여 가중치 부여 질량 행렬을 근사한다.
• 가중치 조정 내적이 타당한 노름을 유도하도록 보장하여 에너지 안정성을 보장한다.
• 요소 내에서 불연속적이거나 매끄럽게 변화하는 파장 속도를 갖는 음향파 방정식에 방법을 적용한다.
• 비다항식 파장 속도에 대해 정확한 보존성을 복원하기 위해 질량 행렬에 저랭크 보정을 적용한다.

실험 결과

연구 질문

• RQ1 국소적으로 이질적인 파장 속도에 대해 요소별 질량 행렬을 저장하지 않으면서도 에너지 안정성과 고차 정확도를 유지할 수 있는 DG 방법을 제안할 수 있는가?
• RQ2 가중치 부여된 $L^2$ 내적을 가중치 조정 내적으로 대체할 경우 안정성과 수렴 성질이 유지되는가?
• RQ3 WADG 형식에서 적분 차수를 얼마나 낮출 수 있으며, 정확도나 안정성에 영향을 주지 않는가?
• RQ4 WADG 방법에서 보존 오차는 어떻게 행동하는가? 매끄러운 해에 대해 초수렴하는가?
• RQ5 제한된 메모리 자원을 갖는 GPU 아키텍처에서 WADG 방법을 효율적으로 구현할 수 있는가?

주요 결과

• 매끄러운 해에 대해 WADG 방법은 $L^2$-노름에서 최적의 수렴 속도 $O(h^{N+1})$를 달성하며, 표준 DG 성능과 동일하다.
• 모수 $N=4$, $h=1/8$, $T=0.5$일 때 수치 실험에서 WADG와 표준 DG 해는 시각적으로 구분되지 않는다.
• 최적 수렴을 위해 적분 차수 $2N+1$이 충분하며, 이를 초과해도 오차에 유의미한 개선이 없다.
• 매끄러운 파장 속도에 대해 WADG의 보존 오차는 $O(h^{2N+2})$ 속도로 초수렴하며, 보존 성질의 높은 정확도를 나타낸다.
• 수치 실험 결과, WADG의 $L^2$ 오차는 표준 DG와 유사하거나 고차수에서 약간 우수한 편이다.
• 질량 행렬에 저랭크 보정을 적용하면 비다항식 파장 속도에 대해 정확한 보존성을 복원할 수 있으며, 물리적 정확성을 유지한다.