[논문 리뷰] Efficient IMEX Runge-Kutta methods for nonhydrostatic dynamics

Andrew Steyer, Christopher J. Vogl외 2명|arXiv (Cornell University)|2019. 06. 17.

Meteorological Phenomena and Simulations참고 문헌 35인용 수 4

한 줄 요약

이 논문은 수평적으로는 명시적, 수직적으로는 암시적인(HEVI) 분할 방식을 사용하는 비수압 대기 모델의 효율적인 시간 적분을 위해 설계된 새로운 암시-명시적 룬게쿠타(IMEX RK) 방법의 가족인 IMKG2 및 IMKG3를 소개한다. 이 방법들은 넓은 안정성 영역을 갖추고 있어 제3차 정확도를 달성하며, 기존의 IMEX 방법들과 비교해 HOMME-NH 모델에서 뛰어난 효율성과 안정성을 보였다.

ABSTRACT

We analyze the stability and accuracy (up to third order) of a new family of implicit-explicit Runge-Kutta (IMEX RK) methods. This analysis expedites development of methods with various balances in the number of explicit stages and implicit solves. We emphasize deriving methods with large stability regions for horizontally explicit vertically implicit (HEVI) partitionings of nonhydrostatic atmosphere models. The IMKG2 and IMKG3 families of IMEX RK methods are formulated within this framework. The HOMME-NH model with a HEVI partitioning is used for testing the accuracy and stability of various IMKG2-3 methods. The efficiency of several IMKG2-3 methods is demonstrated in HOMME-NH and compared to other IMEX RK methods in the literature.

연구 동기 및 목표

• 강한 수직 동역학으로 인해 발생하는 다중 시간 스케일을 다룰 수 있는 비수압 대기 모델에 대한 효율적인 시간 적분 방법을 개발하기 위해.
• 반응성 있는 시간 적분에서의 안정성과 정확도 문제를 해결하기 위해, 넓은 안정성 영역을 갖춘 IMEX 룬게쿠타 스킴을 설계하기 위해.
• IMEX 방법에서 명시적 단계와 암시적 단계의 균형을 최적화하여 계산 비용을 줄이면서도 정확도를 유지하기 위해.
• 실제 비수압 대기 모델(HOMME-NH)에서 HEVI 분할 방식을 적용한 신규 IMEX 방법의 성능을 평가하기 위해.
• 단계 수와 안정성 특성을 사용자 정의할 수 있는 IMEX RK 방법을 유도하는 체계적인 프레임워크를 제공하기 위해.

제안 방법

• 저자는 계수 벡터 α, β, α̂, β̂, δ̂로 매개변수화된 임의의 내부 단계 수를 가진 새로운 IMEX 룬게쿠타 방법의 가족(IMKG2 및 IMKG3)을 제안한다.
• 이 방법들은 수평 방향의 이동은 명시적으로, 수직 동역학은 암시적으로 처리하는 HEVI 분할 방식을 위해 설계되었으며, 이는 계산 비용을 줄이면서도 안정성을 유지한다.
• 정확도와 안정성 분석은 제3차까지 수행되었으며, A-안정성, 강성 감쇠(SD), 강한 안정성 보존(SSP) 특성에 중점을 두었다.
• 계수는 순서 조건과 유리한 안정성 영역을 확보하기 위해 유도된 부처 표기법 프레임워크를 사용하여 구성되었다.
• IMKG2 및 IMKG3 가족은 다양한 변형(예: 232a, 242b, 354a)을 포함하며, 이는 수치 시험 결과를 바탕으로 안정성과 정확도 간의 트레이드오프를 고려해 선택되었다.
• HEVI 분할 방식을 적용한 HOMME-NH 비수압 대기 모델이 시간 적분 성능 평가를 위한 시험 플랫폼으로 사용되었다.

실험 결과

연구 질문

• RQ1 비수압 대기 모델에 대해 제3차 정확도를 달성하면서도 넓은 안정성 영역을 유지하는 IMEX 룬게쿠타 방법을 체계적으로 설계할 수 있는가?
• RQ2 특히 HEVI 방식인 명시-암시적 분할 전략이 비수압 모델의 시간 적분 안정성과 효율성에 어떤 영향을 미치는가?
• RQ3 명시적 단계와 암시적 단계의 수를 변화시킬 경우, 비수압 동역학에서 IMEX 방법의 안정성과 정확도에 어떤 영향을 미치는가?
• RQ4 HOMME-NH 모델에 적용했을 때, 새로운 IMKG2 및 IMKG3 방법은 기존 문헌에 있는 IMEX RK 방법들과 비교해 효율성과 안정성 면에서 어떻게 다른가?
• RQ5 비수압 대기 시뮬레이션에서 안정성, 정확도, 계산 비용 간 최적의 트레이드오프를 보이는 IMKG2-3 변종은 무엇인가?

주요 결과

• IMKG2 및 IMKG3 가족의 IMEX 룬게쿠타 방법은 유리한 안정성 특성(예: A-안정성, 일부 변종에 대한 강성 감쇠(SD))을 갖추고 있으며 제3차 정확도를 달성한다.
• HOMME-NH 모델에서 IMKG2의 254c 변종가 가장 뛰어난 안정성과 효율성의 균형을 보였으며, 정확도 대비 계산 비용 측면에서 다른 IMEX 방법들을 능가했다.
• HEVI 분할을 적용한 IMKG2-3 방법은 표준 IMEX 방법보다 뛰어난 안정성을 보였으며, 특히 단계 크기 제약 없이 빠른 수직 중력파를 해석하는 데 유리했다.
• IMKG3의 354a 변종는 넓은 안정성 영역과 순서 감소 없이 뛰어난 안정성과 정확도를 확보하여 장시간 적분에 적합했다.
• 243a 및 253a 변종는 뛰어난 SSP 성질을 보였으며, 급격한 기울기나 불연속성이 있는 시뮬레이션에서 강건성을 높였다.
• 종합적으로, IMKG2-3 방법은 HOMME-NH 시뮬레이션에서 계산 비용을 줄이면서도 고차 정확도와 안정성을 유지함으로써 뚜렷한 효율성 향상을 보였다.