[논문 리뷰] Analysis of an asymptotic preserving scheme for stochastic linear kinetic equations in the diffusion limit

Nathalie Ayi, Erwan Faou|arXiv (Cornell University)|2018. 03. 16.

Gas Dynamics and Kinetic Theory참고 문헌 18인용 수 9

한 줄 요약

이 논문은 다중성 노이즈를 가진 스토크래틱 선형 운동 방정식에 대해 새로운 점근 보존(AP) 스킴을 제안한다. 이 스킴은 분포 함수를 매크로스코픽 및 마이크로스코픽 성분으로 분리하는 마이크로-매크로 분해를 사용한다. 스킴은 텔레그래프 모델과 일반 케이스에서 명시적 CFL 조건 하에 ε → 0일 때도 운동론적 및 확산 영역에서 정확하고 효율적인 시뮬레이션을 보장하는 바, 균일하게 안정함을 증명하였다.

ABSTRACT

We present an asymptotic preserving scheme based on a micro-macro decomposition for stochastic linear transport equations in kinetic and diffusive regimes. We perfom a mathematical analysis and prove that the scheme is uniformly stable with respect to the mean free path of the particles in the simple telegraph model and in the general case. We present several numerical tests which validate our scheme.

연구 동기 및 목표

• 다중성 노이즈를 가진 스토크래틱 선형 운동 방정식에 대해 운동론적 및 확산 영역에서 모두 안정적이고 정확한 수치 스킴을 개발하는 것.
• 점근 보존(AP) 프레임워크를 다중성 위너 노이즈를 가진 스토크래틱 편미분 방정식에 확장하여, AP 문헌에서의 새로운 기여를 이룩하는 것.
• 스케일링 파라미터 ε(평균 자유 경로)에 대해 균일한 안정성을 증명하여 ε → 0일 때의 강건성을 확보하는 것.
• 다양한 ε 값에서 명시적 및 크랭크-니콜슨 스킴과의 비교를 통해 수치적으로 스킴을 검증하는 것.

제안 방법

• 충돌 연산자의 보존 및 평형 성질에 기반하여 분포 함수 f를 매크로스코픽 성분 ρ와 마이크로스코픽 성분 g로 분해하는 마이크로-매크로 분해를 사용한다.
• 고정된 영역에서 안정성을 확보하기 위해 연립 마이크로-매크로 시스템에 대해 반명시적 시간 이산화를 적용한다.
• 자기수반 연산자 Q가 합산 조건을 만족하도록, 원통형 위너 과정 dWt를 통해 다중성 노이즈를 통합한다.
• 이토 미적분과 스토크래틱 그론발라 유형의 추론을 사용하여 에너지 추정을 유도하고, CFL 조건 하에서 안정성을 확립한다.
• ε → 0일 때 확산 해법으로 수렴하는 스토크래틱 AP 스킴(SMM)을 구현하여 점근적 행동을 유지한다.
• 주기적 경계 조건을 가진 유한 차분 공간 이산화를 사용하고, 100회의 스토크래틱 실현을 통해 집합 평균을 계산한다.

실험 결과

연구 질문

• RQ1 다중성 노이즈를 가진 스토크래틱 선형 운동 방정식에 대해 점근 보존 스킴을 구성할 수 있는가?
• RQ2 제안된 스킴은 텔레그래프 모델과 일반 케이스에서 평균 자유 경로 ε에 대해 균일하게 안정한가?
• RQ3 스킴은 운동론적(ε = 1) 및 확산적(ε = 10⁻²) 영역 모두에서 정확성과 안정성을 유지하는가?
• RQ4 일반 케이스에서 안정성을 보장하는 CFL 조건은 무엇이며, 이는 ε 및 공간 해상도에 어떻게 의존하는가?
• RQ5 SMM 스킴의 성능은 명시적 및 크랭크-니콜슨 스킴과 비교해 볼 때 평균 밀도 진화 측면에서 어떻게 다른가?

주요 결과

• 텔레그래프 모델에서 스킴은 ε에 대해 균일하게 안정하며, 에너지 추정과 이토 미적분을 통해 증명되었다.
• 일반 케이스에서 안정성은 명시적 CFL 조건 하에 확립되었다: Δt ≤ (σ̃Δx²)/(2(2+ε)) + (εΔx)/(2+ε), 이는 에너지 노름의 성장이 유계임을 보장한다.
• 에너지 추정은 매크로스코픽 밀도 ρ와 마이크로스코픽 변동성 g의 기대 제곱 L² 노름이 시간에 따라 유계로 유지됨을 보여주며, 이는 노이즈 계수에 의존하는 상수에 의해 제어된다.
• 수치 실험은 SMM 스킴이 운동론적(ε=1) 및 확산적(ε=10⁻²) 영역 모두에서 평균 밀도 ρ를 정확히 캡처함을 확인하였으며, 기준 솔루션과 일치함을 보였다.
• 스킴은 ε → 0의 극한에서 확산 해법으로 수렴하여, 스토크래틱 설정에서의 점근 보존 성질을 확인하였다.
• 스킴은 유한 시간 간격 동안 안정성을 유지하지만, 스토크래틱 에너지 추정의 성격상 증명이 무한 시간으로 확장되지 않는다.