[논문 리뷰] Trefftz discontinuous Galerkin methods on unstructured meshes for the wave equation

Andrea Moiola|arXiv (Cornell University)|2015. 05. 01.

Advanced Numerical Methods in Computational Mathematics참고 문헌 9인용 수 3

한 줄 요약

이 논문은 공간-시간 티프트스 불연속 갈레르킨 방법을 비정렬 메시와 임의의 면 방향을 가진 파동 방정식에 적용한다. 잘 정의된 문제성, 안정성, 그리고 DG 및 $L^2$ 노름에서 최적의 사전 오차 한계를 증명하며, 파동 방정식을 요소별로 만족하는 티프트스 기저 함수 덕분에 고차 정확도를 달성하고 반추적 방법을 통해 국소적으로 효율적인 해를 계산할 수 있다.

ABSTRACT

We describe and analyse a space-time Trefftz discontinuous Galerkin method for the wave equation. The method is defined for unstructured meshes whose internal faces need not be aligned to the space-time axes. We show that the scheme is well-posed and dissipative, and we prove a priori error bounds for general Trefftz discrete spaces. A concrete discretisation can be obtained using piecewise polynomials that satisfy the wave equation elementwise.

연구 동기 및 목표

• 공간-시간 면이 정렬되어 있지 않은 비정렬 메시에서 파동 방정식에 대해 안정적이고 수렴하는 공간-시간 불연속 갈레르킨 방법을 개발한다.
• 기존의 티프트스 DG 방법을 고차원 및 더 일반적인 메시, 특히 공간적 면과 시간적 면을 포함한 메시로 확장한다.
• 일반적인 티프트스 이산 공간에 대해 DG 및 $L^2$ 노름에서 사전 오차 한계를 증명한다.
• 메시의 조건에 따라 텐트 피치 알고리즘을 사용해 반추적 방법으로 해를 국소적으로 순차적으로 업데이트할 수 있음을 보여준다.

제안 방법

• 파동 방정식을 첫째 단계 시스템으로 표현하며, 이는 이동의 시간 도함수인 $v$와 이동의 기울기인 ${\boldsymbol{\sigma}}$를 변수로 사용한다.
• 모든 함수가 각 요소 $K \in \mathcal{T}_h$ 내부에서 국소적으로 파동 방정식을 만족하도록 티프트스 이산 공간 $\mathbf{V}(\mathcal{T}_h)$를 정의한다.
• 약한 형태는 메시 면에서 수치적 플럭스를 사용하며, 딜리클레, 뉴먼, 루빈 경계에서의 안정화를 위해 펜alty 항을 적용한다.
• 에너지 추정과 쌍대성 논증을 통해 방법이 잘 정의되고 소산적이며, 안정성과 수렴성을 보장한다.
• 공간적 면 조건을 만족하는 메시의 일군에 대해, 텐트 피치 알고리즘을 사용해 반추적이고 요소별로 해를 계산할 수 있는 전략을 제공한다.
• 다항수 티프트스 공간은 $P_{\ell,j}({\mathbf{x}},t) = ({\mathbf{x}} \cdot {\mathbf{d}}_j - ct)^\ell$ 형태의 다항수 파동을 사용해 구성되며, 이는 국소적으로 파동 방정식을 만족함을 보장한다.

실험 결과

연구 질문

• RQ1 비정렬 공간-시간 메시와 임의의 면 방향을 가진 파동 방정식에 대해 티프트스 불연속 갈레르킨 방법을 구성하고 분석할 수 있는가?
• RQ2 공간적 면과 시간적 면을 포함한 일반적인 메시 조건 하에서 제안된 티프트스-DG 방법이 잘 정의되고 안정적인가?
• RQ3 일반적인 티프트스 이산 공간에 대해 DG 및 $L^2$ 노름에서 유도할 수 있는 사전 오차 한계는 무엇인가?
• RQ4 어떤 메시 조건에서 티프트스-DG 방법이 국소적이고 순차적인 업데이트를 허용하는 반추적 방법으로 계산될 수 있는가?
• RQ5 티프트스-DG 방법의 수렴 속도는 자유도와 다항식 차수 측면에서 표준 DG 방법과 비교해 어떻게 되는가?

주요 결과

• 모든 티프트스 이산 공간에 대해 티프트스-DG 방법은 잘 정의되고 소산적이며, 에너지 및 쌍대성 논증을 통해 안정성이 증명된다.
• DG 노름에서 사전 오차 한계가 확립되었으며, 이는 티프트스 공간의 근사 성질과 메시의 정규성에 따라 달라진다.
• 내부 면이 모두 공간적 면인 메시(예: 텐트 피치 알고리즘으로 생성된 메시)에 대해, 반추적이고 국소적으로 계산 가능한 해법 전략을 제공하며, 요소별 순차적 해 업데이트가 가능하다.
• 텐트 피치 메시에서 검증된 충분한 조건 하에, 오차의 $L^2$ 노름은 데이터와 잔여항의 $L^2$ 노름에 비례하는 상수 배로 유계이다.
• 차수 $p$의 다항수 티프트스 공간은 차원이 $\mathcal{O}(p^n)$이며, 전체 다항수 공간($\mathcal{O}(p^{n+1})$)보다 훨씬 낮아 효율성이 높다.
• 해석적 해에 대해 $p$에 대해 지수 수렴이 예상되며, $n=1$일 경우 $h_K$와 $p$에 대해 최적의 대수적 수렴 속도가 증명되었으며, 표준 DG 방법보다 높은 효율성을 보인다.