[논문 리뷰] An efficient dynamical low-rank algorithm for the Boltzmann-BGK equation close to the compressible viscous flow regime

Lukas Einkemmer, Jingwei Hu외 1명|arXiv (Cornell University)|2021. 01. 18.

Model Reduction and Neural Networks참고 문헌 36인용 수 5

한 줄 요약

이 논문은 볼츠만-비거키(Boltzmann-BGK) 방정식에 대한 새로운 동적 저질서 적분기법을 제안하며, $f = Mg$로 표현된 해를 통해 압축성 레이놀즈 수 나비에-스토크스 근사에 효율적으로 도달한다. 여기서 $M$은 맥클라우린 분포이고, 저질서 근사법은 오직 $g$에만 적용된다. 이 방법은 압축성 영역에서도 정확하고 저질서의 시뮬레이션을 가능하게 하며, 상당히 감소된 질서 요구 조건과 충격파 및 급격한 기울기를 잘 다룬다.

ABSTRACT

It has recently been demonstrated that dynamical low-rank algorithms can provide robust and efficient approximation to a range of kinetic equations. This is true especially if the solution is close to some asymptotic limit where it is known that the solution is low-rank. A particularly interesting case is the fluid dynamic limit that is commonly obtained in the limit of small Knudsen number. However, in this case the Maxwellian which describes the corresponding equilibrium distribution is not necessarily low-rank; because of this, the methods known in the literature are only applicable to the weakly compressible case. In this paper, we propose an efficient dynamical low-rank integrator that can capture the fluid limit -- the Navier-Stokes equations -- of the Boltzmann-BGK model even in the compressible regime. This is accomplished by writing the solution as $f=Mg$, where $M$ is the Maxwellian and the low-rank approximation is only applied to $g$. To efficiently implement this decomposition within a low-rank framework requires, in the isothermal case, that certain coefficients are evaluated using convolutions, for which fast algorithms are known. Using the proposed decomposition also has the advantage that the rank required to obtain accurate results is significantly reduced compared to the previous state of the art. We demonstrate this by performing a number of numerical experiments and also show that our method is able to capture sharp gradients/shock waves.

연구 동기 및 목표

• 기존의 저질서 방법이 맥클라우린 분포의 비저질서 성질로 인해 실패하는, 압축성 점성 유동 영역에서 볼츠만-비거키 방정식을 효율적으로 시뮬레이션하는 데 도전하는 것.
• 작은 쿤델 수($\varepsilon \ll 1$)의 유체 영역과 전이 영역을 포함한 광범위한 쿤델 수 범위에서 정확성과 효율성을 유지하는 강력한 수치적 프레임워크를 개발하는 것.
• 저질서 근사를 통해 물리적 거동을 유지하면서 급격한 기울기와 충격파를 정확하게 해결할 수 있도록 하는 것.
• 해를 $f = Mg$로 분해함으로써 기존 방법보다 훨씬 낮은 질서로도 정확한 시뮬레이션을 가능하게 하여, 필요한 질서가 상당히 감소함을 입증하는 것.

제안 방법

• 해를 $f = Mg$로 분해하며, $M$은 평형 맥클라우린 분포이고 $g$는 비평형 성분으로, 저질서 근사는 오직 $g$에만 적용된다.
• 프로젝터-스플릿팅 방법에 기반한 동적 저질서 적분기를 사용하여 $g$를 시간에 따라 진화시키며, 저질서 구조를 유지한다.
• 이소othermal 경우, 빠른 푸리에 변환(FFT) 기반 스펙트럼 방법을 이용해 적분형 계수를 효율적으로 평가한다.
• 공간 및 속도 이산화에 스펙트럼 콜로케이션 유한차분(SCFD) 방법을 적용하고, 이동항의 안정성을 확보하기 위해 업윈딩 기법을 사용한다.
• 체이프먼-엔스코 변환을 활용하여, 소규모 $\varepsilon$ 근사에서 나비에-스토크스 방정식과의 점근적 일致성을 정당화한다.
• 적응형 시간 스텝을 구현하고, 고질서($r=60$) 또는 고해상도 스펙트럼 방법으로 계산된 기준 솔루션과의 비교를 통해 방법을 검증한다.

실험 결과

연구 질문

• RQ1 맥클라우린 분포가 저질서가 아니더라도, 볼츠만-비거키 방정식의 압축성 나비에-스토크스 근사에 정확하게 도달할 수 있는 동적 저질서 적분기를 설계할 수 있는가?
• RQ2 $f = Mg$ 분해가 직접 $f$에 대한 저질서 근사와 비교해 필요한 질서를 상당히 감소시키는가?
• RQ3 이 방법은 저질서 근사를 유지하면서도 급격한 기울기와 충격파를 정확하게 해결할 수 있는가?
• RQ4 공간적으로 비균일한 $\varepsilon$ 분포를 포함한 다양한 쿤델 수에서 이 방법의 성능은 어떠한가?
• RQ5 유체 영역($\varepsilon \to 0$)과 전이 영역($\varepsilon \sim 0.01$)에서 이 방법의 수렴 특성과 정확도는 어떠한가?

주요 결과

• $\varepsilon = 10^{-5}$일 때, 질서 $r=3$에서도 기준 솔루션과 구분되지 않을 정도로 정확한 결과를 도출하여, 매우 낮은 질서로도 유체 영역에서 충분히 정확한 시뮬레이션이 가능함을 보여준다.
• $\varepsilon = 0.1$ 및 $\varepsilon = 0.5$일 때, 비균형 빔 편포의 지수 감쇠 특성(예측된 $\sim \exp(-t/\varepsilon)$)을 정확히 포착하며, 이는 이론적 예측과 일치한다.
• 공간적으로 변화하는 $\varepsilon$를 가진 빔 안정화 문제에서, 열화 속도가 $\varepsilon$에 따라 적절히 변화함을 확인하여, 다양한 영역 간의 물리적 일致성을 입증한다.
• 공간적으로 변화하는 $\varepsilon$의 경우, 질서 $r=10$일 때 기준 솔루션($r=60$)과 잘 일치하는 반면, $r=3$일 경우 상당한 오차가 발생함을 확인하여, 정확도를 확보하기 위한 실용적 질서 한계를 제시한다.
• 폭발 문제에서 반경 방향으로 확장되는 충격파를 성공적으로 해결하였으며, 낮은 질서에서도 날카운 충격파 전면을 유지함을 보였다.
• 이소othermal 경우에 복합형 계수 평가를 FFT 기반으로 구현함으로써 $f = Mg$ 분해의 효율적 구현이 가능해졌으며, 계산 비용이 감소하였다.