[논문 리뷰] Partially implicit Runge-Kutta methods for wave-like equations in spherical-type coordinates

I. Cordero-Carrión, P. Cerdá–Durán|arXiv (Cornell University)|2012. 11. 26.

Numerical methods for differential equations참고 문헌 13인용 수 4

한 줄 요약

이 논문은 구형 좌표계와 유사한 좌표계에서 파동 유사 편미분방정식의 시간 진화를 안정화하기 위해 부분적으로 암시적 룬게쿠타 방법을 제안한다. 이는 강성 항이나 좌표계에 기인한 요인으로 인한 수치적 불안정성을 해결한다. 이 방법은 명시적 룬게쿠타 방법보다 안정성을 향상시키면서도 장기 시뮬레이션에서 고차 정확도를 유지한다.

ABSTRACT

Partially implicit Runge-Kutta methods are presented in this work in order to numerically evolve in time a set of partial differential equations. These methods are designed to overcome numerical instabilities appearing during the evolution of a system of equations due to potential numerical unstable terms in the sources, such as stiff terms or the presence of factors as a result of a partic- ular chosen system of coordinates. In this article, partially implicit Runge-Kutta methods for several convergence orders have been derived and stability properties have been analyzed. These methods are shown to be appropriated to avoid the de- velopment of numerical instabilities in the evolution in time of wave-like equations in spherical-type coordinates, in contrast to the explicit Runge-Kutta methods.

연구 동기 및 목표

• 강성 있는 소스 항이나 좌표계 선택으로 인해 발생하는 시간 진화 과정에서의 수치적 불안정성을 해결한다.
• 구형 유사 좌표계에 특화된 고차 부분적으로 암시적 룬게쿠타 스킴을 개발한다.
• 명시적 방법이 실패하는 장기 시뮬레이션에서 안정성과 정확도를 확보한다.
• 제안된 방법의 안정성 특성을 분석하여 그 견고성을 검증한다.

제안 방법

• 파동 유사 PDE의 시간 적분을 위한 다중 수렴 차수를 가진 부분적으로 암시적 룬게쿠타 방법을 유도한다.
• 소스 내 강성 있거나 불안정한 항을 암시적 구성요소로 분리하여 시간 진화를 안정화시킨다.
• 좌표 특이점이나 요인으로 인해 불안정성이 유발될 수 있는 구형 유사 좌표계에서의 편미분방정식 시스템에 방법을 적용한다.
• 룬게쿠타 스킴의 계수를 구조화하기 위해 부터 테이블 표현을 사용한다.
• von Neumann 안정성 분석을 수행하여 유도된 방법의 안정성 특성을 평가한다.
• 표준 명시적 룬게쿠타 방법과 비교하여 부분적으로 암시적 스킴의 성능과 안정성을 평가한다.

실험 결과

연구 질문

• RQ1 강성 있는 항이나 좌표계에 기인한 요인으로 인한 파동 유사 방정식의 시간 적분에서 발생하는 수치적 불안정성은 어떻게 완화할 수 있는가?
• RQ2 구형 유사 좌표계에서 파동 유사 PDE에 부분적으로 암시적 룬게쿠타 방법을 적용했을 때의 안정성 행동은 어떠한가?
• RQ3 명시적 룬게쿠타 방법보다 개선된 안정성을 확보하면서도 고차 정확도를 유지할 수 있는가?
• RQ4 강성 있는 시스템에서 안정성을 보장하기 위해 부분적으로 암시적 룬게쿠타 스킴의 최적 계수는 무엇인가?
• RQ5 구형 좌표계의 맥락에서 제안된 방법의 안정성 영역은 명시적 룬게쿠타 방법의 안정성 영역과 어떻게 비교되는가?

주요 결과

• 제안된 부분적으로 암시적 룬게쿠타 방법은 강성 있거나 좌표계에 기인한 항으로 인해 발생하는 시간 진화 과정에서의 수치적 불안정성을 효과적으로 억제한다.
• 안정성 분석 결과, 명시적 룬게쿠타 방법이 유사 조건에서 발산하는 것과는 달리, 이 방법들은 장기적인 통합 시간 동안 안정성을 유지함을 확인하였다.
• 유도된 스킴에서 고차 수렴성이 유지되어 정확한 장기 시뮬레이션을 가능하게 한다.
• 이 방법들은 구형 유사 좌표계에서 발생하는 특수한 수학적 과제, 예를 들어 특이점이나 반경 성분 등에 특별히 대응하도록 설계되었다.
• 부분적으로 암시적 형식은 명시적 방법보다 더 큰 안정적인 시간 스텝을 허용하여 계산 효율성을 향상시킨다.
• 제안된 스킴의 안정성 영역은 특히 소스 항의 강성 성분에 대해 명시적 룬게쿠타 방법보다 상당히 넓다.