[논문 리뷰] Legendre-Gauss-Lobatto grids and associated nested dyadic grids

Kolja Brix, Claudio Canuto외 1명|arXiv (Cornell University)|2013. 10. 31.

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

한 줄 요약

이 논문은 레지스터-갈-로블레르토(Legendre-Gauss-Lobatto, LGL) 격자와 연관된 새로운 종류의 중첩된 이진격자(family of nested dyadic grids)를 제안하며, 고차수 스펙트럴 및 비연속 갈레르킨 방법에서 비중첩되는 LGL 노드 문제를 해결한다. LGL 격자의 정밀한 분석적 성질을 도출하고, 재귀적 정밀화 알고리즘을 통해 이진격자를 구성함으로써, 결과적으로 생성된 격자가 LGL 격자와 국소적으로 균일하게 등가되며 최적의 기수(cardinality)를 유지함을 입증한다. 이는 고차수 PDE 해법기의 강력한 조건부 조절(preconditioning)을 가능하게 한다.

ABSTRACT

Legendre-Gauss-Lobatto (LGL) grids play a pivotal role in nodal spectral methods for the numerical solution of partial differential equations. They not only provide efficient high-order quadrature rules, but give also rise to norm equivalences that could eventually lead to efficient preconditioning techniques in high-order methods. Unfortunately, a serious obstruction to fully exploiting the potential of such concepts is the fact that LGL grids of different degree are not nested. This affects, on the one hand, the choice and analysis of suitable auxiliary spaces, when applying the auxiliary space method as a principal preconditioning paradigm, and, on the other hand, the efficient solution of the auxiliary problems. As a central remedy, we consider certain nested hierarchies of dyadic grids of locally comparable mesh size, that are in a certain sense properly associated with the LGL grids. Their actual suitability requires a subtle analysis of such grids which, in turn, relies on a number of refined properties of LGL grids. The central objective of this paper is to derive just these properties. This requires first revisiting properties of close relatives to LGL grids which are subsequently used to develop a refined analysis of LGL grids. These results allow us then to derive the relevant properties of the associated dyadic grids.

연구 동기 및 목표

• 레지스터-갈-로블레르토(Legendre-Gauss-Lobatto, LGL) 격자가 중첩되어 있지 않다는 근본적인 한계를 해결함으로써, 고차수 비연속 갈레르킨 방법에서 효율적인 조건부 조절을 가능하게 하기 위해.
• LGL 격자와 국소적으로 균일하게 등가되는 중첩된 이진격자 계층을 개발하여 보조공간 방법((auxiliary space method))에서 보조 문제의 안정적이고 효율적인 해법을 가능하게 하기 위해.
• LGL 격자의 정밀한 분석적 성질, 특히 노드 간격과 간격 크기의 단조성(monotonicity)을 확립하여, 관련된 이진격자를 구성하는 기초를 마련하기 위해.
• 결과로 생성된 이진격자가 LGL 격자에 대해 최적의 기수를 유지하고 강력한 국소 등가성 경계를 만족함을 증명하기 위해.
• 중첩성을 유지하면서 메쉬 품질과 국소 균일성을 보존하는 이진격자를 생성하는 구조적 알고리즘을 제공하기 위해.

제안 방법

• LGL 격자의 정밀한 성질, 특히 간격 크기의 단조성과 경계 근처에서의 노드 집중 현상 등을 재검토하고 도출함. 이를 위해 체비셰프-갈-로블레르토 노드와의 연결 고리를 활용함.
• 국소적으로 LGL 격자와의 등가성을 유지하기 위해 필요한 경우에만 간격을 분할하는 재귀적 격자 정밀화 알고리즘(알고리즘 2)을 제안함.
• 임의의 이진격자 D와 LGL 격자 Δ가 겹치는 경우, 비율 |Δ|/|D|가 α와 LGL 간격 크기 추정에서의 상수 Cg에 의존하는 상수로 유계임을 통해 국소 균일성을 확립함.
• 정밀화 과정에서 LGL 간격 크기의 단조성(정리 4.2)을 이용하여, 더粗한 격자에서의 간격 크기가 정밀화 과정에서 증가하지 않음을 증명함으로써 안정성을 확보함.
• 유도된 LGL 노드 성질에 기반하여, 제안된 이진격자가 요구되는 국소 등가성 및 기수 경계를 만족함을 보여주기 위해, 정리 3을 적용함.
• 정밀화 임계값을 제어하기 위해 매개변수 α를 도입함으로써, 중첩성을 유지하기 위해 추가로 필요한 최소한의 점만을 추가함.

실험 결과

연구 질문

• RQ1 국소적으로 균일하게 등가된 LGL 격자와 중첩된 이진격자 계층을 구성할 수 있는가?
• RQ2 이진격자가 국소 균일성과 최적 기수를 유지하기 위해 필요한 LGL 격자의 정밀한 분석적 성질은 무엇인가?
• RQ3 고차수 비연속 갈레르킨 방법과 보조공간 조건부 조절의 맥락에서 LGL 격자의 비중첩성 문제는 어떻게 해결할 수 있는가?
• RQ4 동일한 영역 내에서 이진격자 간격과 해당 LGL 간격 사이의 관계는 무엇이며, 이를 균일하게 유계로 둘 수 있는가?
• RQ5 결과로 생성된 이진격자에 포함된 점의 수가 원래 LGL 격자와 비교하여 유사한가를 입증할 수 있는가?

주요 결과

• 알고리즘 2에 의해 생성된 이진격자는 α와 LGL 간격 크기 추정에서의 상수 Cg에 의존하는 상수로 유계된 비율 |Δ|/|D|를 만족하는 국소 (A,B)-균일 등가성을 가짐.
• 이진격자의 기수는 #D_N ≃ #G^LGL_N를 만족함으로써, 점의 수가 LGL 격자와 동일한 渐近적 비율로 증가함을 의미함.
• 수치적 도표(그림 7)를 통해 중첩성을 확보하기 위해 추가로 추가된 점의 수가 최소한임을 확인함.
• 증명은 LGL 간격 크기의 정밀화 과정에서의 단조성에 기반하며, 이를 체비셰프-갈-로블레르토 노드와의 비교를 통해 확립함.
• 이진격자는 보조공간 방법에 사용하기 위한 핵심 조건을 충족함: 중첩되어 있으며, LGL 격자와 국소적으로 메쉬 크기 비교가 가능함.
• 결과는 고차수 스펙트럴 및 비연속 갈레르킨 방법에서 강력한 조건부 조절기를 구성하는 이론적 기초를 확인함.