[논문 리뷰] Coarse-grained analysis of a lattice Boltzmann model for planar streamer fronts

Wim Vanroose, Giovanni Samaey외 1명|ArXiv.org|2007. 01. 02.

Lattice Boltzmann Simulation Studies참고 문헌 3인용 수 4

한 줄 요약

이 논문은 평면 스크리머 프론트를 위한 격자 붕소 모델에서 속도에 의존하는 이온화 반응을 포함하여 이동파 해를 계산하기 위해 군집화된 시간 스텝퍼 방법을 제안한다. 짧은 시뮬레이션과 리프팅 및 제한 단계를 조합함으로써, 최소 물리적 파동 속도를 제공하는 固定점 해를 계산하는 방법으로, 기울기가 급격한 모델에 대해 전통적인 찰프먼-엔스코그 분석의 강력한 대안을 제공한다.

ABSTRACT

We study the traveling wave solutions of a lattice Boltzmann model for the planar streamer fronts that appear in the transport of electrons through a gas in a strong electrical field. To mimic the physical properties of the impact ionization reaction, we introduce a reaction matrix containing reaction rates that depend on the electron velocities. Via a Chapman--Enskog expansion, one is able to find only a rough approximation for a macroscopic evolution law that describes the traveling wave solution. We propose to compute these solutions with the help of a coarse-grained time-stepper, which is an effective evolution law for the macroscopic fields that only uses appropriately initialized simulations of the lattice Boltzmann model over short time intervals. The traveling wave solution is found as a fixed point of the sequential application of the coarse-grained time-stepper and a shift-back operator. The fixed point is then computed with a Newton-Krylov Solver. We compare the resulting solutions with those of the approximate PDE model, and propose a method to find the minimal physical wave speed.

연구 동기 및 목표

• 미세구조적 전자 충격 이온화 물리법을 속도 해상도를 갖춘 격자 붕소 모델을 통해 통합함으로써 최소 스크리머 모델을 확장하기 위해.
• 전자 밀도 및 전기장의 급격한 기울기를 포착하는 데에 전통적인 찰프먼-엔스코그 분석이 제한됨을 해결하기 위해.
• 정확한 매크로스코픽 PDE를 유도하기 어려운 비이상적인 운동학 모델에서 이동파 해를 계산하기 위한 계산 방법을 개발하고 검증하기 위해.
• 군집화된 시간 스텝퍼와 뉴턴-크로이로프 솔버를 사용하여 스크리머 프론트의 최소 물리적 파동 속도를 수치적으로 결정하기 위해.
• 군집화된 시간 스텝퍼와 찰프먼-엔스코그 전개를 통해 유도된 근사 PDE 모델 간의 비교를 위해.

제안 방법

• 다중 속도를 갖는 격자 붕소 모델을 구성하여 전자 동역학을 시뮬레이션하고, 미세구조적 이온화 단면적에 기반한 속도에 의존하는 반응 속도를 통합한다.
• 군집화된 시간 스텝퍼(CGTS)는 매크로스코픽 장을 기반으로 분포를 초기화하는 리프팅 단계, 군집화된 시간 간격 ΔT 동안 LBM을 시뮬레이션하는 단계, 매크로스코픽 장을 추출하는 제한 단계로 구성된 순서이다.
• 이동파 해는 파동 전파를 고려하기 위해 시프트-백 연산자를 결합한 CGTS의 고정점으로서 발견된다.
• 고정점을 계산하기 위해 뉴턴-크로이로프 솔버를 사용하여 이동파 프로파일에 수렴할 수 있다.
• 이 방법은 매크로스코픽 PDE의 명시적 유도를 피함으로써 급격한 기울기가 있는 경우에도 정확한 시뮬레이션을 가능하게 한다.
• 찰프먼-엔스코그 전개를 통해 유도된 근사 PDE 모델의 결과와의 비교를 통해 방법의 유효성을 검증한다.

실험 결과

연구 질문

• RQ1 군집화된 시간 스텝퍼는 속도에 의존하는 이온화 반응을 포함한 격자 붕소 모델에서 이동파 해를 정확하게 계산할 수 있는가?
• RQ2 미세구조적 모델에서 이온화 단면적의 강도가 최소 물리적 파동 속도에 어떻게 영향을 미치는가?
• RQ3 군집화된 시간 스텝퍼는 전통적인 찰프먼-엔스코그 전개에 비해 임계 파동 속도를 예측하는 데 얼마나 우수한가?
• RQ4 이 방법은 공진향 및 반대진향 스크리머 프론트, 특히 양전하 프론트를 모두 탐지할 수 있는가?
• RQ5 고차수 리프팅 절차는 군집화된 시간 스텝퍼의 정확도에 어떤 영향을 미치는가?

주요 결과

• 군집화된 시간 스텝퍼는 속도에 의존하는 이온화 반응을 포함한 격자 붕소 모델에서 평면 스크리머 프론트의 이동파 해를 성공적으로 계산한다.
• CGTS와 시프트-백 연산자를 결합한 고정점 해를 찾음으로써 최소 물리적 파동 속도를 수치적으로 계산할 수 있다.
• 계산된 파동 속도는 이온화 단면적의 강도에 의존하며, 이 의존성은 군집화된 접근법을 통해 정량적으로 분석할 수 있다.
• 군집화된 방법의 결과는 찰프먼-엔스코그 전개를 통해 유도된 근사 PDE 모델의 결과와 일치하여 이 방법의 타당성을 검증한다.
• 이 방법은 공진향 및 반대진향 스크리머 프론트, 양전하 프론트를 포함한 모든 종류의 프론트를 탐지할 수 있을 정도로 강력하다.
• 고차수 리프팅 절차는 초기 분포 함수를 더 잘 근사함으로써 정확도를 향상시키고, 시뮬레이션 단계를 줄이며 수렴성을 향상시킬 것으로 기대된다.