[논문 리뷰] Spectral Analysis of Continuous FEM for Hyperbolic PDEs: Influence of Approximation, Stabilization, and Time-Stepping

Sixtine Michel, Davide Torlo외 2명|arXiv (Cornell University)|2021. 01. 01.

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

한 줄 요약

이 논문은 1차원 초구속 PDE에 대한 연속 유한요소법의 체계적인 스펙트럼 분석을 제시하며, 근사 공간(베르누이, 등간격, 가우스-로바토), 안정화 기법(SUPG, CIP, LPS), 시간 적분기법(RK, SSPRK, DeC) 간의 상호작용을 평가한다. 주요 기여는 완전 이산화된 푸리에 분석을 통해 최적의 CFL 수치와 안정화 파rameter를 규명한 것으로, 입체 요소에 SSPRK와 CIP 또는 LPS 안정화 기법을 조합할 경우 정확도, 안정성, 계산 효율성의 최적의 균형을 달성함을 결론으로 내린다.

ABSTRACT

We study continuous finite element dicretizations for one dimensional hyperbolic partial differential equations. The main contribution of the paper is to provide a fully discrete spectral analysis, which is used to suggest optimal values of the CFL number and of the stabilization parameters involved in different types of stabilization operators. In particular, we analyze the streamline-upwind Petrov–Galerkin stabilization technique, the continuous interior penalty (CIP) stabilization method and the orthogonal subscale stabilization (OSS). Three different choices for the continuous finite element space are compared: Bernstein polynomials, Lagrangian polynomials on equispaced nodes, and Lagrangian polynomials on Gauss-Lobatto cubature nodes. For the last choice, we only consider inexact quadrature based on the formulas corresponding to the degrees of freedom of the element, which allows to obtain a fully diagonal mass matrix. We also compare different time stepping strategies, namely Runge–Kutta (RK), strong stability preserving RK (SSPRK) and deferred correction time integration methods. The latter allows to alleviate the computational cost as the mass matrix inversion is replaced by the high order correction iterations. To understand the effects of these choices, both time-continuous and fully discrete Fourier analysis are performed. These allow to compare all the different combinations in terms of accuracy and stability, as well as to provide suggestions for optimal values discretization parameters involved. The results are thoroughly verified numerically both on linear and non-linear problems, and error-CPU time curves are provided. Our final conclusions suggest that cubature elements combined with SSPRK and CIP or OSS stabilization are the most promising combinations.

연구 동기 및 목표

• 다양한 이산화 선택에 따라 1차원 초구속 법칙에 대한 연속 갈레르킨 FEM의 안정성과 정확도를 분석하기.
• 스펙트럼 분석을 통해 최적의 CFL 수치와 안정화 파ram터(SUPG, CIP, LPS)를 규명하기.
• 베르누이 다항식, 등간격 라그랑주 다항식, 가우스-로바토 입체 요소를 포함한 다양한 유한요소 공간의 성능을 비교하며, 질량 행렬 구조와 적분 방법에 초점 맞추기.
• 시간 적분 기법인 룬게-쿠타(RK), 강안정성 보존 룬게-쿠타(SSPRK), 연기 보정(DeC)을 평가하며 질량 행렬 처리 방식과 계산 비용에 초점 맞추기.
• 선형 및 비선형 문제, 특히 얕은 수면 방정식에 대해 이론적 결과를 수치적으로 검증하기.

제안 방법

• 모든 기법 조합에 대해 시간 연속 및 완전 이산화된 푸리에 분석(von Neumann)을 수행하여 분산 및 감쇠 오차 평가.
• 세 가지 안정화 기법 분석: SUPG(스트림라인-업윈 페트로프-갈레르킨), CIP(연속 내부 페널티), LPS(국소 투영 안정화)로 각각 질량 행렬과 일致성에 다른 영향을 미침.
• 세 가지 유한요소 공간 비교: 베르누이 다항식, 등간격 노드에 대한 라그랑주 다항식, 가우스-로바토(입체) 노드에 대한 라그랑주 다항식으로, 비정확 적분을 통해 대각질량 행렬 달성.
• 세 가지 시간 적분 방법 평가: 표준 RK, SSPRK(안정성 향상을 위한), DeC(질량 행렬 역행렬 계산을 피하기 위해 반복 보정 방식)
• 상대 오차 ηu를 최소화하기 위한 스펙트럼 분석 프레임워크를 활용하여 각 기법에 최적의 안정화 및 CFL 파ram터 도출.
• 선형 이동, 버거스 방정식, 얕은 수면계에 대해 수치적 검증을 수행하며 오차-CPU 시간 곡선을 통한 성능 비교.

실험 결과

연구 질문

• RQ1 1차원 초구속 PDE에 대한 연속 FEM에서 안정성과 정확도를 확보하기 위해 최적의 CFL 수치와 안정화 파ram터(SUPG, CIP, LPS)는 무엇인가?
• RQ2 특히 비정확 적분을 사용하는 입체 요소와 같은 다양한 유한요소 공간은 스펙트럼 성질과 계산 비용에 어떤 영향을 미치는가?
• RQ3 시간 적분 기법(RK, SSPRK, DeC)의 선택은 질량 행렬 역행렬 계산과 관련하여 안정성, 정확도, 효율성에 어떤 영향을 미치는가?
• RQ4 완전 이산화 기법의 스펙트럼 성질(분산 및 감쇠)은 선형 및 비선형 문제에서 관측된 수치적 수렴 속도와 어느 정도 관련이 있는가?
• RQ5 초구속 문제에서 정확도, 안정성, 계산 비용 간의 최적 균형을 이룰 수 있는 요소 유형, 안정화 기법, 시간 적분기법의 조합은 무엇인가?

주요 결과

• 입체 요소(가우스-로바토 노드에 비정확 적분 적용)는 완전히 대각질량 행렬을 제공하여 역행렬 계산 없이 효율적인 해법이 가능하며, 다른 요소 유형에 비해 스펙트럼 정확도와 안정성에서 뛰어남.
• 입체 요소에 SSPRK 시간 적분과 CIP 또는 LPS 안정화 기법을 조합할 경우 정확도, 안정성, 계산 효율성의 최적 균형을 달성함을 오차-CPU 시간 곡선을 통해 확인.
• DeC 시간 적분 기법은 질량 행렬 역행렬 계산을 피하기 위해 고차 보정 반복을 사용하여 계산 비용을 감소시키지만, 특히 베르누이 다항식과 고차 요소에서는 더 높은 분산 및 감쇠 오차를 유발함.
• 모든 기법에서 최적의 CFL 수치는 SSPRK의 경우 약 2.0, 고차 DeC의 경우 1.5로, 스펙트럼 분 析에서 상대 오차 ηu를 최소화함으로써 유도됨.
• 스펙트럼 분 析 결과, 입체 요소는 럼프드 방법(예: DeC)보다 낮은 분산 오차를 보이며, 전체 질량 행렬 방법보다 더 뛰어난 안정성을 보이며, 특히 고주파 모드에서 유리함.
• 양건 수면 방정식에 대한 수치 결과는 이론적 예측을 확인함: 입체 요소에 LPS 및 CIP 안정화 기법을 적용할 경우 기대되는 수렴 차수를 달성하며, 반면 베르누이 다항식은 특히 P2 및 P3에서 DeC와 함께 최적 수렴을 이룰 수 없음.