[논문 리뷰] On deterministic approximation of the Boltzmann equation in a bounded domain

Francis Filbet|arXiv (Cornell University)|2011. 06. 06.

Gas Dynamics and Kinetic Theory참고 문헌 25인용 수 27

한 줄 요약

이 논문은 유한 도메인에서 주기적, 반사적, 확산적 경계 조건을 갖는 볼츠만 방정식을 해결하기 위한 완전히 결정론적이고 고차수의 수치 방법을 제시한다. 이는 공간에 대해 두 번째 차수의 유한 체적 방법과 충돌 항에 대해 푸리에 스펙트럼 방법을 조합하여 속도 공간에서 스펙트럼 정밀도를 달성하고, 다양한 쿠드센 수 범위에서 효율적인 시뮬레이션을 가능하게 하며, 유체역학적 극한과 희박한 영역 모두를 포함한다.

ABSTRACT

In this paper we present a fully deterministic method for the numerical solution to the Boltzmann equation of rarefied gas dynamics in a bounded domain for multi-scale problems. Periodic, specular reflection and diffusive boundary conditions are discussed and investigated numerically. The collision operator is treated by a Fourier approximation of the collision integral, which guarantees spectral accuracy in velocity with a computational cost of $N\\,\\log(N)$, where $N$ is the number of degree of freedom in velocity space. This algorithm is coupled with a second order finite volume scheme in space and a time discretization allowing to deal for rarefied regimes as well as their hydrodynamic limit. Finally, several numerical tests illustrate the efficiency and accuracy of the method for unsteady flows (Poiseuille flows, ghost effects, trend to equilibrium).

연구 동기 및 목표

• 유계 영역에서 다스케일 행동을 보이는 공간-시간에 의존하는 볼츠만 방정식을 위한 완전히 결정론적인 수치 방법을 개발하기 위해.
• 주기적, 반사 반사, 확산(맥스웰 유형) 조건을 포함한 다양한 경계 조건을 정확히 처리하기 위해.
• 두 번째 차수의 유한 체적 방법과 스펙트럼 방법을 조합하여 공간과 속도 모두에서 고차수 정밀도를 달성하기 위해.
• 시간 분리 기반의 분리 기법을 통해 희박한 영역에서부터 유체역학적 영역에 이르기까지 전체 쿠드센 수 범위에서의 시뮬레이션을 가능하게 하기 위해.
• 연속체 모델이 실패하는 영역, 특히 고스트 효과와 평형으로의 경향성에 대해 알려진 점근적 극한과의 검증을 위해.

제안 방법

• 볼츠만 방정식의 운반 항에 대해 두 번째 차수의 유한 체적 방법을 사용하여 공간 이산화를 수행한다.
• 충돌 적분에 대해 푸리에 스펙트럼 방법을 적용하여 속도 공간에서 스펙트럼 정밀도를 확보하고, $N\log N$의 계산 비용을 유도한다.
• 운반과 충돌 단계를 연산자 분리 기법을 통해 분리함으로써 효율적인 시간 적분과 병렬 처리를 가능하게 한다.
• 다중 척도, 특히 유체역학적 극한을 포함한 안정성과 정밀도를 유지하는 시간 이산화를 구현한다.
• 경계 조건은 하이브리드 맥스웰 유형 모델을 통해 처리한다: 비율 $\alpha$는 맥스웰 분포로 재방출되고, $1-\alpha$는 반사된다.
• 벽면의 에너지 흡수 계수 $\mu(t,x)$는 질량 보존 조건을 통해 구하고, 물리적 일관성을 확보한다.

실험 결과

연구 질문

• RQ1 유계 영역에서 볼츠만 방정식에 대해 완전히 결정론적인 방법이 공간과 속도 양쪽에서 고차수 정밀도를 달성할 수 있는가?
• RQ2 작은 쿠드센 수를 갖는 주기적 온도 구동 흐름에서 고스트 효과를 얼마나 정확히 포착할 수 있는가?
• RQ3 유체역학적 극한에서 수치 해가 고전적 열전도 방정식이 아닌 점근적 이론으로 수렴하는가?
• RQ4 하나의 알고리즘과 최소한의 격자 해상도로 희박한 영역과 유체역학적 영역을 모두 효율적으로 시뮬레이션할 수 있는가?
• RQ5 표준 적분 기반 충돌 해법 대비 속도 공간에서의 스펙트럼 방법은 정밀도와 계산 비용 측면에서 어떻게 비교되는가?

주요 결과

• 이 방법은 공간에서 두 번째 차수 정밀도와 속도에서 스펙트럼 정밀도를 달성하여, 비교적 적은 수의 속도 격자 포인트로도 고해상도 결과를 도출할 수 있다.
• 포아죄유 흐름의 수치 시뮬레이션 결과는 이론적 예측과 벤치마크 솔루션과 뛰어난 일치를 보였다.
• 고스트 효과 문제에서는 쿠드센 수 $\varepsilon$가 감소함에 따라 온도장이 증가하며, Sone 등이 제시한 점근적 이론으로 수렴함을 확인하였고, 열전도 방정식의 해로 수렴하지는 않았다.
• 한편, 평균 속도장은 $\varepsilon \to 0$ 극한에서 0이 되며, 점근적 이론과 일치하여 임의의 흐름이 존재하지 않음을 확인하였다.
• 주기적 벽면 온도 문제에서 비평형 온도 프로파일을 성공적으로 포착하였고, 열전도 방정식이 예측한 등온선과는 다름을 확인하였다.
• 충돌 단계의 계산 비용은 $N\log N$ 비례로 증가하여 고해상도 속도 격자에 대해 효율적이며, 분리된 구조 덕분에 병렬 구현이 가능하다.