[논문 리뷰] Implicit-explicit Runge-Kutta schemes and applications to hyperbolic systems with relaxation

L. Pareschi, G. Russo|arXiv (Cornell University)|2010. 09. 14.

Computational Fluid Dynamics and Aerodynamics인용 수 6

한 줄 요약

이 논문은 고차수, 점근적 보존성(Asymptotic-Preserving)인 음이 아닌-양성(implicit-explicit, IMEX) 룬게쿠타 스킴을 초래하는 하이퍼볼릭 시스템에 제안한다. 이 스킴은 비스티프한 유동항에 대해 강건한 안정성 보존성(Strong-Stability-Preserving, SSP)을 갖는 명시적 스킴과, 스티프한 리콜리케이션 항에 대해 L-안정성 있는 대각형 은닉 룬게쿠타(DIRK) 방법을 조합하여, 비스티프 및 스티프 영역 모두에서 고차수 정확도를 유지한다. 또한 WENO 공간 재구성 기법을 통해 고해상도 충격 포착 기능을 실현하며, 얕은 수면, 교통 흐름, 미세입자 기체 등 다양한 응용 분야에서 저차수 방법보다 뛰어난 성능을 보여준다.

ABSTRACT

We consider new implicit-explicit (IMEX) Runge-Kutta methods for hyperbolic systems of conservation laws with stiff relaxation terms. The explicit part is treated by a strong-stability-preserving (SSP) scheme, and the implicit part is treated by an L-stable diagonally implicit Runge-Kutta methods (DIRK). The schemes proposed are asymptotic preserving (AP) in the zero relaxation limit. High accuracy in space is obtained by Weighted Essentially Non Oscillatory (WENO) reconstruction. After a description of the mathematical properties of the schemes, several applications will be presented.

연구 동기 및 목표

• 스티프한 리콜리케이션 항을 갖는 하이퍼볼릭 시스템에 대해 고차수, 점근적 보존성(Asymptotic-Preserving, AP) 시간 통합 스킴을 개발하는 것.
• 비스티프한 유동항에 대해 강건한 안정성 보존성(SSP)을 갖는 명시적 스킴과, 스티프한 소스항에 대해 L-안정성 있는 대각형 은닉 룬게쿠타(DIRK) 방법을 조합하는 것.
• 스킴이 리콜리케이션 항이 0에 수렴할 때(즉, 점근적 정확도)에도 여전히 정확도의 순서를 유지하도록 보장하는 것.
• 충격 포착을 위해 가중치 기반 본질적으로 진동이 없는(WENO) 재구성 기법을 사용하여 고해상도 공간 이산화를 달성하는 것.
• 얕은 수면 유동, 교통 역학, 미세입자 기체 역학 등 실제 물리적 응용 사례에서 스킴을 검증하는 것.

제안 방법

• IMEX 룬게쿠타 프레임워크는 시스템을 명시적 및 은닉적 부분으로 분리한다: 유동항 ∂xF(U)는 명시적으로 처리하고, 스티프한 리콜리케이션 항 R(U)/ε는 은닉적으로 처리한다.
• 명시적 단계 값은 유한한 변동성과 불연속성 근처의 진동 방지를 위해 강건한 안정성 보존성(SSP) 스킴을 사용하여 계산한다.
• 은닉적 단계 값은 스티프 영역에서 안정성을 확보하기 위해 L-안정성 있는 대각형 은닉 룬게쿠타(DIRK) 방법을 사용하여 계산한다.
• 고차수 공간 정확도는 WENO 재구성을 통해 최소한의 진동을 유도하는 보존형 유한체적 또는 유한차분 이산화 기법을 통해 달성된다.
• 스킴이 점근적 보존성(AP)임을 증명하였으며, 이는 ε → 0일 때 수치해가 올바른 평형 시스템으로 수렴함을 의미한다.
• 미세입자 기체 응용 사례에서 발생하는 비선형 대수계를 효율적으로 해결하기 위해 뉴턴 방법을 사용한다.

실험 결과

연구 질문

• RQ1 스티프한 리콜리케이션 극한(ε → 0)에서도 정확도 순서를 유지하는 고차수 IMEX 룬게쿠타 스킴을 구성할 수 있는가?
• RQ2 명시적 부분에서는 강건한 안정성 보존성(SSP) 성질을 유지하면서도 은닉적 부분에서는 L-안정성을 확보할 수 있는가?
• RQ3 WENO 재구성 기법은 저차수 공간 스킴에 비해 하이퍼볼릭 시스템에 대한 충격 해상도를 얼마나 향상시키는가?
• RQ4 다양한 영역에서 전통적인 분리 방법에 비해 제안된 IMEX 스킴은 정확도와 안정성 측면에서 어떻게 성능을 발휘하는가?
• RQ5 미세입자 기체 모델과 같은 복잡한 비선형 소스항을 갖는 시스템에서도 과도한 수치적 확산 없이 스킴이 처리할 수 있는가?

주요 결과

• 제안된 IMEX-SSP2-WENO 및 IMEX-SSP3-WENO 스킴은 각각 비스티프 및 스티프 영역에서 제2차 및 제3차 정확도를 달성하며, 제3차 스킴은 충격의 블러링을 현저히 감소시킨다.
• ε = 0.2인 교통 흐름 모델에서 제3차 스킴은 제2차 스킴보다 충격을 더 날카롭게 포착했으며, 진동은 최소화되고 밀도 및 유량 프로파일이 일관되게 유지되었다.
• ε = 0.01인 미세입자 기체 모델에서 제3차 스킴은 제2차 스킴보다 더 낮은 수치적 확산을 유지하며 충격 구조를 잘 보존했지만, 보존 변수 설정으로 인해 압력 프로파일에 소규모 비물리적 진동이 나타났다.
• 스킴은 점근적 보존성 행동을 보였으며, ε → 0일 때 수치해가 평형 시스템으로 수렴하여 감소된 모델과의 일致성을 확인했다.
• 뉴턴 방법을 사용함으로써 미세입자 기체 사례에서 비선형 은닉계를 효율적으로 해결할 수 있었으며, 소스항의 비선형성에도 불구하고 안정성과 신뢰성을 확보했다.
• 수치적 결과는 스킴이 전통적인 스트랑 분리 방법보다 스티프 영역에서 뛰어난 성능을 보였음을 확인하였으며, 특히 분리 방법이 정확도 저하 문제를 겪는 상황에서 두드러진다.