[논문 리뷰] Enhanced conservation properties of Vlasov codes through coupling with conservative fluid models

Thomas Alexander Trost, S. Lautenbach외 1명|arXiv (Cornell University)|2017. 02. 01.

Gas Dynamics and Kinetic Theory참고 문헌 21인용 수 3

한 줄 요약

이 논문은 운동론적 코드를 순간 기반 피드백을 통해 보정함으로써 비틀림 시뮬레이션에서 에너지 보존을 향상시키는 새로운 결합 방법을 제안한다. 일반화된 10모멘트 피팅 절차를 적용함으로써 수치적 가열을 크게 감소시키면서도 운동론적 특성을 유지하며, 10모멘트 피팅은 등방성 또는 5모멘트 접근 방식과 달리 안정적이고 효과적임을 입증한다.

ABSTRACT

Many phenomena in collisionless plasma physics require a kinetic description. The evolution of the phase space density can be modeled by means of the Vlasov equation, which has to be solved numerically in most of the relevant cases. One of the problems that often arise in such simulations is the violation of important physical conservation laws. Numerical diffusion in phase space translates into unphysical heating, which can increase the overall energy significantly, depending on the time scale and the plasma regime. In this paper, a general and straightforward way of improving conservation properties of Vlasov schemes is presented that can potentially be applied to a variety of different codes. The basic idea is to use fluid models with good conservation properties for correcting kinetic models. The higher moments that are missing in the fluid models are provided by the kinetic codes, so that both kinetic and fluid codes compensate the weaknesses of each other in a closed feedback loop.

연구 동기 및 목표

• 상태공간 이산화와 인위적 가열로 인한 Vlasov 시뮬레이션에서의 수치적 에너지 비보존 문제를 해결하기 위해.
• 운동론적 플라즈마 코드의 보존 성질을 향상시키기 위한 일반적이고 방법에 종속되지 않는 방법을 개발하기 위해.
• 유체 모델의 보존 강점을 운동론적 코드의 고차 모멘트 정밀도와 상호 피드백을 통해 활용하기 위해.
• 선형 및 비선형 플라즈마 영역에서의 순간 기반 보정의 효과성을 테스트하고 검증하기 위해.
• 다중 영역 시뮬레이션에서 보존적 모멘트 피팅을 사용하여 유체 및 운동론적 코드 간의 견고한 결합을 가능하게 하기 위해.

제안 방법

• 유체 모멘트가 운동론적 상태공간 분포를 보정하는 피드백 루프를 갖춘 폐쇄형 시스템이 Vlasov 해석기와 보존적 유체 모델을 연결한다.
• 속도 공간에서의 선형 변환을 사용하여 운동론적 분포를 정해진 유체 모멘트에 맞추며, 이는 이전의 등방성 피팅을 비대칭 케이스로 일반화한 것이다.
• 피팅 행렬은 기울기 하강법을 통해 정확히 계산하거나 계산 비용을 줄이기 위해 근사적으로 계산하며, 두 방법 모두 안정성과 정확도를 테스트하였다.
• 유한차분 Vlasov 해석기에서 CWENO 방법을 사용하여 공간 도함수를 계산하고, 시간 분할 방법을 사용하여 통합하였다.
• 10모멘트 피팅은 밀도, 기체 속도, 에너지 밀도 텐서를 일치시키는 데 사용되었으며, 비교를 위해 5모멘트 피팅도 테스트하였다.
• 시뮬레이션은 GEM 챌린지 벤치마크에서 검증되었으며, 재결합 속도, 에너지 변화, 분포 함수 형태를 비교하였다.

실험 결과

연구 질문

• RQ1 보존적 유체 모델이 인위적 효과를 유발하지 않고 Vlasov 시뮬레이션에서 수치적 가열을 효과적으로 보정할 수 있는가?
• RQ2 등방성 10모멘트 피팅이 등방성 또는 5모멘트 피팅보다 운동론적 특성과 보존 법칙을 얼마나 잘 유지하는가?
• RQ3 피팅 행렬의 정확한 계산과 근사 계산이 에너지 보존과 시뮬레이션 안정성에 미치는 영향은 무엇인가?
• RQ4 이 결합 방법이 다중 영역 플라즈마 시뮬레이션에서의 유체-운동론적 인터페이스 문제에 일반화될 수 있는가?
• RQ5 피드백 루프가 비선형 영역에서 인위적 가열을 크게 줄이면서도 물리적 정확도를 유지하는가?

주요 결과

• 근사 행렬 계산을 사용한 10모멘트 피팅은 보정되지 않은 Vlasov 시뮬레이션 대비 총 에너지 증가를 최대 50% 감소시켜 에너지 보존을 크게 향상시킨다.
• 등방성 피팅 절차는 수치적 불안정성을 유발하므로 실용적 사용에는 부적합하나, 일반 방법의 특수 케이스이기는 하다.
• 5모멘트 피팅은 필수적인 운동론적 특성을 파괴하고 인위적 불안정성을 유도하므로 시뮬레이션은 물리적으로 의미가 없다.
• GEM 챌린지 시뮬레이션 결과, 보정된 Vlasov 코드는 고해상도에서도 기준 솔루션과 유사한 재결합 플럭스와 전류 프로파일을 유지한다.
• 이 방법은 다양한 해상도와 시간 간격에서 견고하며, CWENO 방법의 정규화 파rameter ε의 선택에 대해 거의 민감하지 않다.
• 피드백 루프는 상태공간 확산을 보정함으로써 시뮬레이션을 성공적으로 안정화시키며, 플라즈마의 운동론적 기술을 손상시키지 않는다.