[논문 리뷰] Enforcing strong stability of explicit Runge--Kutta methods with superviscosity

Zheng Sun, Chi‐Wang Shu|arXiv (Cornell University)|2019. 12. 25.

Numerical methods for differential equations참고 문헌 50인용 수 8

한 줄 요약

이 논문은 선형 및 비선형 문제에 대해 명시적 룬게-쿠타(RK) 방법에서 강한 안정성을 확보하기 위해 두 가지 초점성(viscosity) 기반 방법—수정된 방법과 필터링 방법—을 제안한다. 고차수의 확산성 또는 분산성 초점성 항을 추가함으로써, 정확도를 저하시키지 않고 에너지 감쇠를 유지할 수 있으며, 강한 안정성을 확보하기 위한 초점성 강도에 대한 날카운 경계를 유도하였다. 이 경계는 선형 및 비선형 사례(보존법칙의 DG 공간 이산화 포함)에 대해 수치적으로 검증되었다.

ABSTRACT

A time discretization method is called strongly stable, if the norm of its numerical solution is nonincreasing. It is known that, even for linear semi-negative problems, many explicit Runge--Kutta (RK) methods fail to preserve this property. In this paper, we enforce strong stability by modifying the method with superviscosity, which is a numerical technique commonly used in spectral methods. We propose two approaches, the modified method and the filtering method for stabilization. The modified method is achieved by modifying the semi-negative operator with a high order superviscosity term; the filtering method is to post-process the solution by solving a diffusive or dispersive problem with small superviscosity. For linear problems, most explicit RK methods can be stabilized with either approach without accuracy degeneration. Furthermore, we prove a sharp bound (up to an equal sign) on diffusive superviscosity for ensuring strong stability. The bound we derived for general dispersive-diffusive superviscosity is also verified to be sharp numerically. For nonlinear problems, a filtering method is investigated for stabilization. Numerical examples with linear non-normal ordinary differential equation systems and for discontinuous Galerkin approximation of conservation laws are performed to validate our analysis and to test the performance.

연구 동기 및 목표

• 비정규 또는 비강한 강제 조건에서의 명시적 RK 방법의 강한 안정성 부족 문제를 해결하기 위해.
• 원래의 정확도 순서를 유지하면서 에너지 감쇠가 단조롭게 유지되는 안정화 기법을 개발하기 위해.
• 선형 문제에서 강한 안정성을 확보하기 위해 필요한 초점성 강도에 대한 날카운 이론적 경계를 제공하기 위해.
• 지역 해의 행동에 기반한 적응형 초점성을 사용하는 필터링 전략을 통해 비선형 문제로 접근을 확장하기 위해.
• 선형 ODE 및 보존법칙에 대한 DG 스킴을 대상으로 수치 실험을 통해 방법을 검증하기 위해.

제안 방법

• ODE 시스템 내 반음성 연산자에 고차수 초점성 항을 추가하는 수정된 방법을 도입한다.
• 작은 초점성을 가진 확산성 또는 분산성 문제를 통해 해를 후처리하는 필터링 방법을 제안한다.
• 초점성에 의해 유도되는 오차를 제어하기 위해 안정화된 해와 원래 해 간의 에너지 추정을 유도한다.
• 선형 문제에서 강한 안정성을 확보하기 위한 확산성 초점성에 대한 날카운 이론적 경계를 확립하며, 이 경계는 수치적으로 최적임을 입증한다.
• 비선형 문제에 대해 지역 해의 행동에 기반한 적응형 초점성을 사용하는 필터링 방법을 적용한다.
• 초기값 문제에 대한 하이퍼볼릭 보존법칙에 대해 엔트로피-보존적 플럭스를 사용하는 불연속 갈레르킨 공간 이산화를 적용하여 성능을 테스트한다.

실험 결과

연구 질문

• RQ1 초점성을 사용하여 선형 ODE 시스템에 대해 명시적 RK 방법을 안정화시킬 수 있으며, 정확도를 유지하고 강한 안정성을 확보할 수 있는가?
• RQ2 선형 문제에서 강한 안정성을 보장하기 위해 필요한 최소한의 초점성 강도는 얼마이며, 이 경계는 날카로운가?
• RQ3 에너지 감쇠가 보장되지 않는 비선형 문제로 초점성 안정화 기법을 어떻게 확장할 수 있는가?
• RQ4 필터링 방법은 비선형 보존법칙에서 정확도를 유지하면서 비물리적 진동을 억제하는가?
• RQ5 초점성 파rameter의 선택은 장기 시뮬레이션에서 에너지 변화와 해의 품질에 어떤 영향을 미치는가?

주요 결과

• 확산성 초점성에 대한 날카운 경계가 도출되었으며, 이는 선형 문제에서 강한 안정성을 확보하기 위해 필수적이며 충분함을 의미하며, 이 경계는 수치적으로 최적임을 입증하였다.
• 수정된 방법과 필터링 방법은 선형 문제에서 원래의 정확도 순서를 유지하며, 수렴률은 기반 RK 스킴과 일치한다 (예: $p=2$일 때 2차, $p=5$일 때 5차).
• P^2-DG와 2차 RK를 사용한 버거스 방정식에서 $L^2$ 오차는 수렴률 2.01을 보였으며, $P^4$-DG와 5차 펠러버그 방법을 사용한 경우 $L^ atural$ 노름에서 수렴률은 2.01이었다.
• 수치 결과는 원래 방법이 해의 노름이 증가하는 반면, 필터링 방법은 에너지 감쇠를 유지하며 노름이 시간이 지남에 따라 감소함을 보여주었다.
• 적응형 필터링을 통한 초점성은 충격 포착에서 진동을 억제하지만, 초점성 강도가 너무 클 경우 해를 과도하게 부드럽게 만들 수 있다.
• 필터링 방법은 표준 RK 방법이 강한 안정성을 유지하지 못하는 비선형 문제(예: 버거스 방정식)에서도 성공적으로 안정화시켰으며, 엔트로피-보존적 플럭스를 사용하는 경우에도 효과적이었다.