[논문 리뷰] Optimized Runge-Kutta Methods with Automatic Step Size Control for Compressible Computational Fluid Dynamics

Hendrik Ranocha, Lisandro Dalcín외 2명|arXiv (Cornell University)|2021. 04. 14.

Numerical methods for differential equations참고 문헌 84인용 수 51

한 줄 요약

이 논문은 압축성 유체역학(CFD) 시뮬레이션을 위한 적응형 시간스텝 제어 기반의 최적화된 명시적 룬게쿠타 방법을 제시한다. 기존의 고정된 CFL 수치가 아닌 오차 기반 제어기를 사용하며, 안정성에 의해 제한되는 영역에서는 거의 최적의 효율을 달성하고, 정확도에 의해 제한되는 영역에서는 정확한 시간 오차 제어를 제공한다. 이는 산업용 CFD 벤치마크와 도전적인 점성 충격 파동 시뮬레이션 모두에서 기존의 CFL 기반 방법을 능가한다.

ABSTRACT

We develop error-control based time integration algorithms for compressible fluid dynamics (CFD) applications and show that they are efficient and robust in both the accuracy-limited and stability-limited regime. Focusing on discontinuous spectral element semidiscretizations, we design new controllers for existing methods and for some new embedded Runge-Kutta pairs. We demonstrate the importance of choosing adequate controller parameters and provide a means to obtain these in practice. We compare a wide range of error-control-based methods, along with the common approach in which step size control is based on the Courant-Friedrichs-Lewy (CFL) number. The optimized methods give improved performance and naturally adopt a step size close to the maximum stable CFL number at loose tolerances, while additionally providing control of the temporal error at tighter tolerances. The numerical examples include challenging industrial CFD applications.

연구 동기 및 목표

• 압축성 CFD 시뮬레이션에서 CFL 기반 시간스텝 제어의 비효율성과 수작업 캘리브레이션의 부담을 해결한다.
• 안정성에 의해 제한되는 영역과 정확도에 의해 제한되는 영역 모두에서 높은 효율성을 유지하는 오차 기반 시간 적분 알고리즘을 개발한다.
• 불연속 스펙트럼 요소 방법에 특화된 최적화된 임bedded 룬게쿠타 쌍과 PID 기반 시간스텝 제어기를 설계한다.
• 오차 기반 제어가 느슨한 정밀도에서 최대 안정성 CFL 수치에 자연스럽게 수렴할 수 있고, 엄격한 정밀도에서 시간 오차 제어를 보장할 수 있음을 입증한다.
• 실제 CFD 응용에서 안정성과 효율성을 보장하기 위해 제어기 파라미터를 선택하는 데 실용적인 지침을 제공한다.

제안 방법

• h-p 적응형 스펙트럼 요소 방법에서 사용할 수 있도록 안정성 영역이 최적화된 새로운 임bedded 룬게쿠타 쌍과 내장된 오차 추정기를 설계한다.
• 국소 절단 오차 추정에 기반해 시간스텝을 조정하는 PID 기반 시간스텝 제어기를 구현하여 안정성과 정확도를 모두 확보한다.
• Butcher 표기법을 사용해 저메모리, 첫 번째 값이 마지막 값과 동일한(FSAL) 룬게쿠타 쌍과 내장된 오차 추정을 정의한다.
• 압축성 오일러 및 라우저-스토크스 방정식의 불연속 갈레르킨 반연속화와 시간 적분기를 통합한다.
• PETSc 기반의 SSDC 솔버 프레임워크 내에서 제어기를 적용하여 확장성 있고, 비정렬, 곡선 메esh 기반 시뮬레이션을 가능하게 한다.
• 수치 실험을 통해 제어기 파라미터(예: 비례, 적분, 미분 이득)를 캘리브레이션하여 강건성과 안정성을 확보한다.

실험 결과

연구 질문

• RQ1 오차 기반 시간스텝 제어는 수작업 CFL 캘리브레이션 없이 안정성에 의해 제한되는 CFD 영역에서 CFL 기반 제어의 효율성을 달성할 수 있는가?
• RQ2 다양한 정밀도에서 최적화된 임bedded 룬게쿠타 쌍과 적응형 제어기를 사용한 방법이 기존 방법에 비해 정확도와 계산 비용 측면에서 어떻게 비교되는가?
• RQ3 고차수, 불연속 스펙트럼 요소 이산화에서 시간스텝 제어의 안정성과 강건성을 보장하기 위해 어떤 제어기 파라미터 설정이 필요한가?
• RQ4 오차 기반 제어는 정밀도가 높은 경우에 안정성 한계 근처에서 거의 최적의 시간스텝을 유지하면서도 시간 오차를 효과적으로 제어할 수 있는가?
• RQ5 새로운 방법은 점성 충격 및 난류 흐름과 같은 도전적인 산업용 CFD 문제에서 어떻게 성능을 발휘하는가?

주요 결과

• 최적화된 오차 기반 방법은 안정성에 의해 제한되는 영역에서 거의 최적의 효율을 달성하며, 느슨한 정밀도에서 최대 안정성 CFL 수치에 가까운 시간스텝을 자동으로 선택한다.
• 엄격한 정밀도에서 방법은 신뢰할 수 있는 시간 오차 제어를 제공하며, 오차가 유한하고 예상되는 수렴 차수를 유지한다. 점성 충격 시험에서 L2 오차 노름을 통해 이를 입증하였다.
• p=4 다항식에서 RK5(4)10F[3S*+] 방법은 정밀도 10−6에서 오직 135,081회의 함수 평가로 L2 오차 5.42×10−5를 달성하여 높은 효율성을 보였다.
• RK4(3)9F[3S*+] 방법은 정밀도 10−8에서 상대적 L2 오차 1.00×10−7를 기록했고, 함수 평가 수는 108,879회로, 뛰어난 정확도와 안정성을 나타낸다.
• 제안된 제어기는 특히 정확도에 의해 제한되는 영역에서 CFL 기반 제어에 비해 효율성과 강건성 측면에서 뛰어나며, 오차 감소가 효과적으로 이루어지지 않는 상황에서의 성능을 뛰어나게 한다.
• 제어기 파라미터(예: β = (0.45, -0.13, 0.00))는 안정성에 매우 중요하며, 부적절한 설정은 시간스텝 거부나 발산을 유도한다. 이는 체계적인 캘리브레이션의 필요성을 강조한다.