[논문 리뷰] Efficient discontinuous Galerkin finite element methods via Bernstein polynomials

Robert C. Kirby|arXiv (Cornell University)|2015. 04. 15.

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

한 줄 요약

이 논문은 단순체 요소에서 국소 기저 함수로 베르슈타인 다항식을 사용하여 초월 보존법칙에 대한 효율적인 비연속 갈레르킨(DG) 유한요소 방법을 제시한다. 요소 질량행렬의 재귀적 블록 구조를 활용하여 질량행렬을 포함한 선형계를 해결하는 최적 복잡도 알고리즘을 개발하였으며, O(n^{d+1})의 복잡도를 달성하면서도 다항식 차수 10 이하에서 높은 정확도를 유지한다.

ABSTRACT

We consider the discontinuous Galerkin method for hyperbolic conservation laws, with some particular attention to the linear acoustic equation, using Bernstein polynomials as local bases. Adapting existing techniques leads to optimal-complexity computation of the element and boundary flux terms. The element mass matrix, however, requires special care. In particular, we give an explicit formula for its eigenvalues and exact characterization of the eigenspaces in terms of the Bernstein representation of orthogonal polynomials. We also show a fast algorithm for solving linear systems involving the element mass matrix to preserve the overall complexity of the DG method. Finally, we present numerical results investigating the accuracy of the mass inversion algorithms and the scaling of total run-time for the function evaluation needed in DG time-stepping.

연구 동기 및 목표

• 초월 보존법칙에 대한 고차수 비연속 갈레르킨(DG) 방법에서 밀도 높은 요소 질량행렬의 역행렬 계산에서 발생하는 계산적 병목 현상을 해결한다.
• 베르슈타인 다항식의 기하학적 및 대수적 구조를 활용하여 DG 항의 최적 복잡도 계산을 가능하게 한다.
• 질량행렬을 포함한 선형계를 해결하는 빠르고 정확한 알고리즘을 개발하면서 전체 O(n^{d+1}) 복잡도를 유지한다.
• 실제 DG 시간적분 시나리오에서 질량행렬 역행렬화 기법의 정확도와 성능을 조사한다.
• 구조적 텐서 기반 방법과 노드 기반 DG 접근법 사이의 격차를, 최적 알고리즘 복잡도를 갖는 이원화된 이산화 기저를 통해 메우고자 한다.

제안 방법

• 단순체 요소에서 국소 유한요소 기저로 베르슈타인 다항식을 사용하고, 그 차원별 재귀적 블록 구조를 활용한다.
• Kirby(2011)의 재귀적 블록 구조를 적용하여 일정 계수 질량행렬을 포함한 선형계를 해결하는 O(n^{d+1}) 알고리즘을 설계한다.
• 직교 다항식의 베르슈타인 표현을 사용하여 질량행렬의 고유값과 고유공간을 특성화한다.
• 2차 및 3차 차원에서 하향적으로 낮은 차원의 콜레스키 분해를 재귀적으로 사용하는 블록별 콜레스키 분해 전략을 구현한다.
• 고정된 반복 횟수(n+1)를 가진 공액 경사법 솔버를 사용하고 정확도 및 수렴성을 콜레스키 방법과 비교한다.
• 고성능 구현 및 구조적 메esh에서의 시간 측정을 위해 Cython과 DOLFIN을 활용한다.

실험 결과

연구 질문

• RQ1 베르슈타인 다항식을 통해 초월 보존법칙에 대한 DG 항의 계산, 특히 질량행렬 역행렬화의 최적 복잡도 계산이 가능할 수 있는가?
• RQ2 블록별 분해 알고리즘이 콜레스키 분해에 비해 DG 질량행렬을 포함한 선형계를 해결할 때 얼마나 정확한가?
• RQ3 실제 메쉬와 다항식 차수에서, 베르슈타인 다항식 기반 솔버를 사용한 DG 함수 평가의 실제 런타임 스케일링은 어떻게 되는가?
• RQ4 고정된 반복 횟수(n+1)를 가진 공액 경사법이 고차수 DG 시뮬레이션에서 충분한 정확도를 달성할 수 있는가?
• RQ5 질량행렬의 재귀적 블록 구조는 정확도를 희생시키지 않고 고성능·저복잡도 DG 솔버를 가능하게 할 정도로 얼마나 유용한가?

주요 결과

• 제안된 블록별 분해 알고리즘은 콜레스키 분해와 유사한 상대 정확도를 확보하며, 다항식 차수 10 이하에서 최소 10자리 정밀도를 유지한다.
• 상대 잔차 허용오차 10^{-12}를 가진 공액 경사법 솔버는 1차원, 2차원, 3차원에서 거의 10자리 정밀도를 달성한다.
• 정확히 n+1회의 반복만 허용할 경우, 공액 경사법의 정확도는 다항식 차수가 높아질수록 심각하게 떨어진다.
• 시간 측정 결과, 32×32 메쉬에서 DG 함수 평가의 스케일링 기울기는 약 1.7로 나타나, 2차원에서 이론적 O(n^3) 복잡도에 비해 예상보다 우수한 성능을 보였다.
• 블록 알고리즘은 높은 정확도를 유지하면서도 구현 복잡도를 줄여주며, 전체 콜레스키 분해에 비해 최소한의 정확도 손실만을 초래한다.
• 현재 알고리즘은 복잡도 측면에서 최적화되어 있으나, 고성능 하드웨어에 최적화되지 않아 구현이 복잡하고 아직 완전히 최적화되지 않은 상태이며, 향후 최적화 여지가 있다.