[논문 리뷰] Time discretisation and $L_2$ stability of polynomial summation-by-parts schemes with Runge-Kutta methods

Hendrik Ranocha, Jan Glaubitz외 2명|arXiv (Cornell University)|2016. 09. 08.

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

한 줄 요약

이 논문은 다항식 합성부분합(SBP) 스킴과 룬게-쿠타(Runge-Kutta) 시간 적분기를 결합하여 완전 이산화된 형식에서 $L_2$ 안정성을 보장하는 데에 확장한다. 보정 절차 재구성(CPR/FR) 기법을 적응적으로 적용하고 특정 룬게-쿠타 방법의 안정성을 분석함으로써, 선형 문제에서 안정성을 유지하는 시간 이산화 조건을 수립한다. 이는 수치 예제를 통해 검증된다.

ABSTRACT

Splitting the discretisation of a partial differential equation by the method of lines in a semidiscretisation in space and a following integration in time, the coupling of these schemes has to be considered for fully discrete stability. This article focuses at first on the extension of adaptive techniques developed by the authors for semidiscretisations in the correction procedure via reconstruction / flux reconstruction framework. Additionally, general results about stability of some Runge-Kutta methods for linear problems are established, followed by numerical examples.

연구 동기 및 목표

• 준이산화된 공간 스킴에서의 적응형 안정화 기법을 완전 이산화된 시간-공간 형식으로 확장하기 위해.
• 다항식 합성부분합(SBP) 공간 이산화와 결합된 룬게-쿠타 방법의 $L_2$ 안정성 분석하기 위해.
• 특정 룬게-쿠타 스킴이 선형 시간에 의존하는 PDE에 대해 안정성을 유지할 수 있는 일반 조건 수립하기 위해.
• 이론적 결과를 선형 모델 문제에서의 안정성과 수렴성에 대한 수치 예제를 통해 검증하기 위해.

제안 방법

• 시간 이산화를 위해 룬게-쿠타 방법을 포함한 보정 절차 재구성(CPR/FR) 기법 적응하기 위해.
• 공간 준이산화에 대해 합성부분합(SBP) 연산자를 적용하여 안정성 확보하기 위해.
• 단계 값 분석과 SBP 성질과의 일致성 분석을 통해 룬게-쿠타 스킴의 안정성 조건 유도하기 위해.
• 에너지 방법을 사용하여 특정 순서와 안정성 조건 하에서 완전 이산화 시스템의 $L_2$ 안정성 증명하기 위해.
• 선형 모델 문제에서 이론적 안정성 결과를 검증하기 위한 수치 실험 구현하기 위해.

실험 결과

연구 질문

• RQ1 룬게-쿠타 시간 적분기가 SBP 공간 이산화와 결합되었을 때 $L_2$ 안정성을 유지하는 조건은 무엇인가?
• RQ2 CPR/FR 기반 공간 스킴은 어떻게 완전 이산화 설정에서 안정성을 확보할 수 있는가?
• RQ3 SBP 연산자와 함께 $L_2$ 안정성을 유지하기 위해 룬게-쿠타 스킴이 가져야 할 순서와 안정성 특성은 무엇인가?
• RQ4 이론적 안정성 조건은 모델 문제에서의 수치 성능와 어떻게 비교되는가?

주요 결과

• 논문은 선형 문제에 대해 완전 이산화된 SBP-룬게-쿠타 스킴의 $L_2$ 안정성에 충분한 조건을 수립한다.
• 특정 대각형 은밀한 룬게-쿠타 방법이 단계 일致성이 적절할 경우 SBP 연산자와 조합되었을 때 안정성을 유지함을 보여준다.
• 시간 적분에 대한 CPR/FR 기법의 확장은 공간 준이산화의 안정성 특성을 유지한다.
• 수치 결과는 이론적 안정성 조건이 시뮬레이션에서 관측된 수렴성과 유한성과 일치함을 확인한다.