[논문 리뷰] Monte Carlo stochastic Galerkin methods for non-Maxwellian kinetic models of multiagent systems with uncertainties

Andrea Medaglia, Andrea Tosin외 1명|arXiv (Cornell University)|2022. 01. 01.

Probabilistic and Robust Engineering Design참고 문헌 51인용 수 9

한 줄 요약

이 논문은 물리 공간에서의 직접 시뮬레이션 몽테카를로(DSMC)와 무작위 공간에서의 스펙트럼 갈레르킨 방법을 조합한 하이브리드 DSMC-스펙트럼 갈레르킨 방법을 제안한다. 이는 비맥스웰 분포를 따르는 다중 에이전트 시스템의 비결정론적 모델에 적용되며, 물리 공간에서는 직접 시뮬레이션 몽테카를로를, 무작위 공간에서는 스펙트럼 갈레르킨 방법을 사용한다. 이 방법은 양수성과 엔트로피 소산 등의 핵심 물리적 성질을 유지하면서도, 단계 함수의 정규화와 비보존 양을 보완하기 위한 서브스티튜션 Fokker-Planck 모델을 통해 무작위 공간에서 스펙트럼 수렴을 달성한다.

ABSTRACT

In this paper, we focus on the construction of a hybrid scheme for the approximation of non-Maxwellian kinetic models with uncertainties. In the context of multiagent systems, the introduction of a kernel at the kinetic level is useful to avoid unphysical interactions. The methods here proposed, combine a direct simulation Monte Carlo (DSMC) in the phase space together with stochastic Galerkin (sG) methods in the random space. The developed schemes preserve the main physical properties of the solution together with accuracy in the random space. The consistency of the methods is tested with respect to surrogate Fokker-Planck models that can be obtained in the quasi-invariant regime of parameters. Several applications of the schemes to non-Maxwellian models of multiagent systems are reported.

연구 동기 및 목표

• 비결정론적 매개변수를 가진 비맥스웰 운동 모델에 대한 수치적 프레임워크를 개발하는 것.
• 불확실성 정량화된 운동 시뮬레이션에서 물리적 성질(양수성, 엔트로피 소산)을 유지하는 데 도전하는 것.
• 비선형성과 비맥스웰성의 상호작용 커널로 인해 무작위 공간에서의 스펙트럼 수렴을 달성하는 것.
• 근사적 비보존 영역에서의 균형 분포를 정확히 모사하기 위해 서브스티튜션 Fokker-Planck 모델을 구축하는 것.
• 단계 함수의 정규화와 비보존 양의 보정을 통해 방법의 강건성과 일관성을 확보하는 것.

제안 방법

• 물리 공간에서 입자 수준의 동역학을 위해 직접 시뮬레이션 몽테카를로(DSMC)를, 무작위 공간에서의 스펙트럼 갈레르킨(sG) 방법을 조합한다.
• DSMC 샘플의 무작위 의존성을 표현하기 위해 일반화 다항식 확장(gPC)을 사용한다.
• sG에서 스펙트럼 수렴을 복원하기 위해 입자 수준에서 단계 함수의 정규화(mollification)를 도입한다.
• 근사적 비보존 영역에서 유도된 서브스티튜션 Fokker-Planck 모델을 사용하여 평균 속도와 에너지와 같은 비보존 양을 보정한다.
• 두 단계 보정: 충돌 규칙의 정규화와 Fokker-Planck 근사에 기반한 모멘트 재스케일링.
• 정확도와 안정성 검증을 위해 일관성 추정과 수렴 연구를 수행한다.

실험 결과

연구 질문

• RQ1 하이브리드 DSMC-sG 방법이 비맥스웰 운동 모델의 비결정론적 특성에서 물리적 성질(양수성, 엔트로피 소산)을 유지할 수 있는가?
• RQ2 비맥스웰 모델에서 단계 함수가 수렴을 방해할 경우, 어떻게 무작위 공간에서의 스펙트럼 수렴을 복원할 수 있는가?
• RQ3 서브스티튜션 Fokker-Planck 모델이 비보존 영역에서 운동 해의 모멘트를 얼마나 정확히 근사하는가?
• RQ4 정규화와 모멘트 보정이 다중 에이전트 시스템에서의 DSMC-sG 스킴의 정확도와 안정성에 어떤 영향을 미치는가?
• RQ5 이 방법이 불확실한 상호작용을 가진 복잡한 모델(예: 부의 분포, 교통 흐름)에서 높은 정확도와 효율성을 달성할 수 있는가?

주요 결과

• 단계 함수의 정규화와 서브스티튜션 Fokker-Planck 모델을 통한 모멘트 보정을 통한 DSMC-sG 방법은 무작위 공간에서 스펙트럼 수렴을 달성한다.
• 도박 모델에서, M1 = M2 = 50일 때 L2 오차 수렴 차수 ≈5.0을 보이며 스펙트럼 정확도를 확인한다.
• 부의 분포 모델에서, 불확실성 하에 파레토 유사 분포에서 가우시안 유사 분포로의 전이를 정확히 포착한다.
• 교통 흐름 모델에서, 정규화와 재스케일링을 적용한 방법은 평균 속도와 에너지의 기대되는 행동을 재현하며, M1 = M2 = 5일 때 수렴 차수 ≈4.5를 보인다.
• 서브스티튜션 Fokker-Planck 모델은 진정한 평형 분포가 알려져 있지 않은 경우에도 균형 분포와 모멘트 진동을 정확히 근사할 수 있다.
• 모든 테스트 케이스에서 높은 차원의 불확실성 조건 하에서도 물리적 일관성(예: 양수성, 유한한 모멘트)을 유지한다.