[논문 리뷰] Well-balanced and asymptotic preserving schemes for kinetic models

Casimir Emako, Min Tang|arXiv (Cornell University)|2016. 03. 10.

Mathematical Biology Tumor Growth참고 문헌 17인용 수 4

한 줄 요약

이 논문은 ε → 0일 때의 확산 극한과 ε = O(1)일 때의 정 steady 상태 행동이라는 두 영역에서 모두 잘 균형 잡힌(WB) 및 점점 수렴하는(AP) 수치적 스킴을 설계하기 위한 일반적인 프레임워크를 제안한다. 방법은 AP 시간 진동 단계와 국소 정 steady 상태 해법을 조합하여 WB 성질을 강제함으로써, ε에 관계없이 균일한 정확도로 두 점점 수렴 극한과 정 steady 상태를 성공적으로 포착한다.

ABSTRACT

In this paper, we propose a general framework for designing numerical schemes that have both well-balanced (WB) and asymptotic preserving (AP) properties, for various kinds of kinetic models. We are interested in two different parameter regimes, 1) When the ratio between the mean free path and the characteristic macroscopic length $ε$ tends to zero, the density can be described by (advection) diffusion type (linear or nonlinear) macroscopic models; 2) When $ε$ = O(1), the models behave like hyperbolic equations with source terms and we are interested in their steady states. We apply the framework to three different kinetic models: neutron transport equation and its diffusion limit, the transport equation for chemotaxis and its Keller-Segel limit, and grey radiative transfer equation and its nonlinear diffusion limit. Numerical examples are given to demonstrate the properties of the schemes.

연구 동기 및 목표

• 운동 방정식에 대해 동시에 잘 균형 잡힌(WB) 및 점점 수렴하는(AP) 성질을 만족하는 통합 수치적 프레임워크를 개발하는 것.
• ε ≪ 1일 때의 확산 극한과 ε = O(1)일 때의 강한 소스 항 영역이라는 두 개별 매개변수 영역에서 운동 모델을 정확하게 시뮬레이션하는 도전 과제를 해결하는 것.
• 특히 표준 스킴이 실패하는 작은 ε 영역에서도 수치적 스킴이 안정적이고 정확하게 유지되도록 보장하는 것.
• ε → 0일 때 매크로스코픽 극한(예: Keller-Segel 또는 비선형 확산)이 이산 수준에서 정확히 유지되도록 보장하는 것.
• 강한 소스 항을 가진 초구속 시스템의 정 steady 상태를 정확히 포착함으로써 장시간 시뮬레이션에서도 높은 정확도를 유지하는 것.

제안 방법

• 프레임워크는 두 단계 절차를 사용한다: 첫째, AP 스킴이 시간에 의존하는 운동 방정식을 풀어 예측값을 생성한다.
• 둘째, 예측된 해를 각 계산 셀 내의 국소 정 steady 상태 문제에 대한 입구 경계 조건으로 사용한다.
• 정 steady 상태 문제의 해법은 매크로스코픽 극한과 일치하도록 특화된 AP 스킴을 사용하여 풀며, 일致성을 보장한다.
• 최종적으로 각 셀에서의 외부 유량을 계산하여 수치적 유량을 보정함으로써, 이산 수준에서 잘 균형 잡힌 성질을 강제한다.
• 이 방법은 중성자 전달, 화학적 이동(Othmer-Dunbar-Alt), 그레이 복사 전달 세 가지 모델에 적용되었으며, 구체적인 이산화 스킴이 제공된다.
• 이 접근법은 기존의 AP 스킴을 시간 진동 및 정 steady 상태 문제에 대해 유연하게 활용할 수 있어, 효율적인 구현이 가능하다.

실험 결과

연구 질문

• RQ1 일관된 수치적 프레임워크가 다양한 매개변수 영역에서 운동 모델에 대해 동시에 잘 균형 잡힌(WB) 및 점점 수렴하는(AP) 성질을 달성할 수 있는가?
• RQ2 ε의 크기와 관계없이 ε → 0일 때 이산 스킴이 정확한 매크로스코픽 극한(예: 확산 또는 비선형 확산)을 유지할 수 있는가?
• RQ3 ε = O(1)일 때, 강한 소스 항을 가진 초구속 시스템의 정 steady 상태 해를 정확히 포착할 수 있는가?
• RQ4 스킴의 수렴 행동은 메쉬 크기 Δx와 매개변수 ε에 대해 어떻게 되며, ε에 관계없이 균일한가?
• RQ5 비선형 시스템에 대한 반복적 해법에 의존하지 않고도 이산 수준에서 잘 균형 잡힌 성질을 강제할 수 있는가?

주요 결과

• 제안된 스킴은 모든 ε 값에서 메쉬 크기 Δx에 대해 균일한 수렴성을 보이며, 그림 7에서 일관된 수렴 차수를 보여줌으로써 점점 수렴(AP) 성질이 확인된다.
• ε → 0일 때의 올바른 극한 행동을 포착한다: 복사 강도 Iε는 ψ⁰/2에 수렴하며, 여기서 ψ⁰ = ac(T⁰)⁴이며, ε에 대한 수렴 차수는 약 0.6에 가까운 것으로 나타났다(그림 8 참조).
• ε = 1일 때, 스킴은 정확한 정 steady 상태 해 I(x,v) = |v|x 및 ρ = x를 성공적으로 포착하며, T = 1일 때 잔여항이 감소하고 안정화됨을 확인하여 잘 균형 잡힌(WB) 성질을 입증한다(그림 9 및 10 참조).
• 수치 결과는 스킴이 ε ≪ 1일 때의 확산 극한과 ε = O(1)일 때의 강한 소스 영역에서 모두 정확성과 안정성을 유지하며, ε이 감소함에 따라 성능 저하 없이 작동함을 확인한다.
• 이 프레임워크는 계산적으로 효율적이며, 예측 단계에서 명시적 또는 단순한 암시적 AP 스킴을 사용할 수 있고, 비선형 시스템에 대한 고비용의 뉴턴 유형 해법을 피할 수 있다.
• 이 방법은 일반적이며 다른 운동 모델, 특히 회복을 가진 초구속 시스템으로의 확장이 가능하여, 다스케일 시뮬레이션 분야에서 광범위한 응용 가능성이 있다.