[논문 리뷰] High order semi-Lagrangian methods for the BGK equation

Maria Groppi, Giovanni Russo외 1명|arXiv (Cornell University)|2014. 11. 28.

Gas Dynamics and Kinetic Theory참고 문헌 15인용 수 4

한 줄 요약

이 논문은 BDF와 DIRK 시간 적분기법을 사용하여 특성선을 따라 BGK 방정식에 대한 고차수 반라그랑주 방법을 제시하며, 희박한 유동 및 유체 영역 모두에서 점근 보존성과 높은 정확도를 달성한다. 특히 보간을 사용하지 않는 BDF3 방법은 수치 실험에서 계산 효율성과 강건성 면에서 뛰어난 성능을 보이며, 비용-성능 균형에서 DIRK 방법을 능가한다.

ABSTRACT

A new class of high-order accuracy numerical methods for the BGK model of the Boltzmann equation is presented. The schemes are based on a semi-lagrangian formulation of the BGK equation; time integration is dealt with DIRK (Diagonally Implicit Runge Kutta) and BDF methods; the latter turn out to be accurate and computationally less expensive than the former. Numerical results and examples show that the schemes are reliable and efficient for the investigation of both rarefied and fluid regimes in gasdynamics.

연구 동기 및 목표

• 희박한 유동 및 유체 영역 모두에서 안정적이고 정확한 고차수 정확도를 갖는 BGK 방정식의 수치적 해법을 개발한다.
• 반라그랑주 설정에서 BGK 방정식에 대한 BDF와 DIRK 시간 적분 기법의 성능을 비교한다.
• 특히 속도 공간에서의 보간 영향을 정확도와 계산 비용 측면에서 조사한다.
• Knudsen 층을 해상도할 필요 없이도 유체 극한(ε → 0)에서 정확도를 유지하는 점근 보존성 스킴을 설계한다.
• 실제 기체역학 응용을 위해 1차원 공간과 3차원 속도 공간으로의 방법 확장을 위해 Chu 감소법을 적용한다.

제안 방법

• 반라그랑주 접근법을 사용하여 특성선을 따라 BGK 방정식을 재구성함으로써 시간 진동을 공간 격자 제약에서 분리한다.
• 시간 적분은 BDF(역차분 공식) 및 DIRK(대각선 음성 룬게쿠타) 방법을 사용하며, BDF3가 더 뛰어난 계산 효율성을 보였다.
• 비균일 분포 함수에 대해 WENO 재구성(WENO23 및 WENO35)을 통해 고차수 공간 정확도를 확보하여 불연속점 근처의 진동을 최소화한다.
• 위상공간 격자에 직접 작용하는 방법의 변종을 제안하여 보간을 회피함으로써 속도 공간에서의 계산 오버헤드를 감소시켰다.
• Chu 감소법을 통해 1차원 공간과 3차원 속도 공간으로 확장하여 차원을 감소시키면서도 정확도를 유지한다.
• 질량, 운동량, 에너지의 보존과 엔트로피 감소를 보장하여 물리적 일관성을 유지한다.

실험 결과

연구 질문

• RQ1 고차수 BDF 및 DIRK 방법을 반라그랑주 방법과 효과적으로 조합하여 BGK 방정식을 고정확도로 해결할 수 있는가?
• RQ2 시간 적분기법 선택(BDF 대비 DIRK)이 희박한 영역과 유체 영역 모두에서 계산 비용과 정확도에 어떤 영향을 미치는가?
• RQ3 속도 공간에서의 공간 보간을 회피할 경우 정확도 손실 없이 계산 효율성이 얼마나 향상되는가?
• RQ4 제안된 스킴이 유체 극한(ε → 0)에서 점근 보존성을 유지하는가? 특히 낮은 Knudsen 수에서 성능이 어떻게 되는가?
• RQ5 WENO 재구성 기법은 분포 함수의 불연속점에서 안정성과 해상도를 어떻게 향상시키는가?

주요 결과

• 보간을 사용하지 않는 BDF3 방법이 모든 시험 케이스에서 정확도와 계산 비용 측면에서 최상의 종합적 성능을 보였다.
• BDF 기반 스킴은 정밀도가 유사한 DIRK 스킴보다 계산 비용이 낮게 나타났다.
• WENO 재구성 기법을 사용한 스킴은 고차수 정확도를 확보하면서도 충격과 불연속점 근처의 비물리적 진동을 억제하였다.
• 스킴은 점근 보존성을 유지하였으며, 작은 ε에서도 정확히 유체 극한(압축성 라우지에르-나비에-스토크스 행동)을 잘 포착하였다.
• 속도 공간에서의 보간이 없어지면서 계산 비용이 크게 감소하였으며, 특히 낮은 속도 격자 해상도 문제에서 두드러졌다.
• 수치 결과는 희박 유동과 근접 평형 상태의 유체역학을 고해상도로 시뮬레이션할 수 있는 스킴의 강건성을 확인하였다.