[논문 리뷰] Strong Stability Preserving Multistep Runge-Kutta Methods

Christopher Bresten, Sigal Gottlieb외 4명|arXiv (Cornell University)|2013. 07. 30.

Computational Fluid Dynamics and Aerodynamics참고 문헌 25인용 수 12

한 줄 요약

이 논문은 단조성과 안정성 특성을 갖는 고차수의 명시적 강안정성 보존(Strong Stability Preserving, SSP) 다단계 룬게쿠타 방법을 개발하여 초구속 편미분방정식(Hyperbolic PDEs)을 해결한다. 순서 조건과 단조성 조건을 유도하고, 2차 이상의 방법에 대해 SSP 계수의 상한을 증명하며, 수치 최적화를 통해 최적화된 방법을 10차까지, 8단계까지, 5단계까지 구성하여 기존 표준 방법보다 훨씬 높은 효과적 SSP 계수를 달성한다.

ABSTRACT

High-order spatial discretizations with strong stability properties (such as monotonicity) are desirable for the solution of hyperbolic PDEs. Methods may be compared in terms of the strong stability preserving (SSP) time-step. We prove an upper bound on the SSP coefficient of explicit multistep Runge--Kutta methods of order two and above. Order conditions and monotonicity conditions for such methods are worked out in terms of the method coefficients. Numerical optimization is used to find optimized explicit methods of up to five steps, eight stages, and tenth order. These methods are tested on the advection and Buckley-Leverett equations, and the results for the observed total variation diminishing and positivity preserving time-step are presented.

연구 동기 및 목표

• 총변동량 감소(Total Variation Diminishing, TVD) 및 양성 보존과 같은 강안정성 특성을 유지하는 고차수의 명시적 다단계 룬게쿠타 방법을 개발하기 위해.
• 다단계 룬게쿠타 방법의 계수에 기반하여 필요한 순서 조건과 단조성 조건을 도출하기 위해.
• 2차 이상의 명시적 다단계 룬게쿠타 방법에 대해 SSP 계수의 상한을 설정하기 위해.
• 수치 알고리즘을 사용하여 순서 2에서 10까지, 최대 5단계 및 8단계까지의 방법을 최적화하여 효과적 SSP 계수를 극대화하기 위해.
• 광고 방정정식과 버클리-레버렛 방정식에 대한 수치 실험을 통해 이론적 SSP 성질을 검증하기 위해.

제안 방법

• 방법 계수를 사용하여 명시적 다단계 룬게쿠타 방법의 순서 조건과 단조성 조건을 유도하기 위해.
• 강안정성 보존(SSP) 성질을 전진 오일러 단계의 볼록 조합으로 공식화하여, 축소된 시간단계 하에서의 안정성을 보장하기 위해.
• SSP 계수 C를 전진 오일러 단계 대비 허용 가능한 최대 시간단계로 정의하고, 효과적 SSP 계수를 C/s로 정의함(여기서 s는 단계 수임).
• 순서 2에서 10까지, 최대 5단계 및 8단계까지의 방법에 대해 수치 최적화를 통해 최대 효과적 SSP 계수를 갖는 방법을 찾기 위해.
• 버클리-레버렛 방정식에 대해 코렌 리미터를 사용하는 보존적 공간 이산화와 광고 방정식에 대해 반-라그랑주 방법을 구현하기 위해.
• 이론적 예측과 비교하여 TVD 및 양성 보존을 위한 관측된 시간단계를 분석하여 SSP 성질을 검증하기 위해.

실험 결과

연구 질문

• RQ1 2차 이상의 명시적 다단계 룬게쿠타 방법에 대해 SSP 계수의 이론적 상한은 무엇인가?
• RQ2 다단계 룬게쿠타 방법에 대해 순서 조건과 단조성 조건을 체계적으로 유도하는 방법은 무엇인가?
• RQ3 제한된 단계 수와 단계 수를 갖는 고차수의 명시적 다단계 룬게쿠타 방법에서 달성 가능한 최대 효과적 SSP 계수는 무엇인가?
• RQ4 수치 실험에서 TVD 및 양성 보존을 위한 관측된 시간단계는 이론적 예측과 어떻게 비교되는가?
• RQ5 최적화된 다단계 룬게쿠타 방법은 초구속 편미분방정식의 안정성과 효율성 측면에서 기존의 룬게쿠타 방법을 능가할 수 있는가?

주요 결과

• 2차 이상의 명시적 다단계 룬게쿠타 방법에 대해 SSP 계수의 상한이 증명되었으며, 이는 볼록 조합 프레임워크 하에서 달성 가능한 최대 안정성에 대한 제약을 둔다.
• 5단계 및 8단계를 갖는 10차 최적화 방법은 효과적 SSP 계수 0.603을 달성하여 비-SSP RK4,4 방법의 0.250보다 훨씬 높은 성능을 보였다.
• 광고 방정식에서 최적화된 (7,3,3) 방법의 관측된 TVD 시간단계는 이론적 값의 1.113배였으며, 이는 이론적 한계를 초월한 강력한 안정성 보존을 나타낸다.
• 버클리-레버렛 문제에서 관측된 양성 보존 시간단계는 항상 TVD 시간단계보다 크게 나타났으며, 이는 안정성 성질의 선택이 관측된 안정성 한계에 영향을 준다는 것을 보여준다.
• (20,3,10) 방법은 관측된 TVD 시간단계 0.356Δx를 기록하여 효과적 SSP 계수 0.110를 달성하였으며, 이는 단계 수를 늘일수록 고차수 방법의 안정성이 향상됨을 보여준다.
• TVD 및 양성 보존을 위한 관측된 시간단계는 단계 수가 증가함에 따라 증가하며, 이는 더 많은 단계가 안정성 특성을 향상시킨다는 이론적 기대를 확인한다.