[논문 리뷰] Some numerical aspects of the conservative PSM scheme in a 4D drift-kinetic code

Jean‐Philippe Braeunig, Nicolas Crouseilles외 4명|arXiv (Cornell University)|2013. 03. 09.

Magnetic confinement fusion research참고 문헌 14인용 수 7

한 줄 요약

이 논문은 fusion 플라즈마 시뮬레이션을 위한 4D 드리프트-키네틱 코드에서 보존형 조각 다항식 재구성 기법(PSM)의 수치적 안정성과 보존 성질을 조사한다. 유한체적 형태의 PSM이 방향 분할 변형보다 최대 원리(maximum principle)를 더 잘 유지하고 물리적으로 의미 없는 진동을 억제함을 보이며, 특히 장시간의 난류 진화 상황에서 그러한 우월성이 두드러진다. 기울기 제한 스플라인(SLS) 제한기 또한 격자 스케일 이하에서 제어된 확산을 통해 안정성을 더욱 향상시킨다.

ABSTRACT

The purpose of this work is simulation of magnetised plasmas in the ITER project framework. In this context, kinetic Vlasov-Poisson like models are used to simulate core turbulence in the tokamak in a toroidal geometry. This leads to heavy simulations because a 6D dimensional problem has to be solved, even if reduced to a 5D in so called gyrokinetic models. Accurate schemes, parallel algorithms need to be designed to bear these simulations. This paper describes the numerical studies to improve robustness of the conservative PSM scheme in the context of its development in the GYSELA code. In this paper, we only consider the 4D drift-kinetic model which is the backbone of the 5D gyrokinetic models and relevant to build a robust and accurate numerical method.

연구 동기 및 목표

• GYSELA 5D 기립킨티크 코드에서 보존형 PSM 기법의 수치적 안정성을 향상시키기 위해 감소된 4D 드리프트-키네틱 모델에서의 행동을 연구하는 것.
• 이산 수준에서 최대 원리가 유지되어 물리적으로 의미 없는 돌출 또는 과도한 낙폭이 발생하지 않도록 보장하는 것.
• 방향 분할 대비 유한체적 형태의 PSM 변형 간 비교를 통해 보존 성질과 장시간 통합에서의 안정성 평가.
• SLS 제한기가 보존 법칙을 위반하지 않으면서 비물리적 진동을 억제하는 데 얼마나 효과적인지 평가하는 것.
• 질량 및 에너지 보존을 유지하면서도 난류 구조를 정확히 포착할 수 있는지에 대한 성능 검증.

제안 방법

• 이온 난류를 위한 전체 5D 기립킨티크 문제를 단순화하기 위해 원통좌표계(r, θ, z, v∥)에서 4D 드리프트-키네틱 모델을 사용한다.
• 보존형 PSM 기법을 방향 분할(DS) 및 유한체적(FV) 형태로 구현하며, 후자는 발산이 없는 유량을 보장한다.
• 격자 스케일 이하에서 제어된 확산을 추가하기 위해 PSM에 기울기 제한 재구성(SLS)을 적용한다.
• 분포 함수의 극값을 시간에 따라 추적하여 최대 원리 유지 여부를 평가하기 위해 초기 가우시안 펄스를 기준 시험으로 사용한다.
• 비선형 포화 초기기와 완전히 발달한 난류 상태를 각각 반영하기 위해 t = 1800 및 t = 4400에서 시뮬레이션을 실행한다.
• 정확도, 진동 수준, 극값 유지 여부를 평가하기 위해 PSM-DS, PSM-FV, SLS-DS, SLS-FV 및 기준 BSL 기법 간 결과를 비교한다.

실험 결과

연구 질문

• RQ1 장시간 시뮬레이션 중에 PSM 기법의 유한체적 형태가 방향 분할 대비 최대 원리를 얼마나 잘 유지하는가?
• RQ2 SLS 제한기가 보존 성질을 위반하지 않으면서 PSM 기법에서 비물리적 진동을 어느 정도 억제하는가?
• RQ3 보존형 PSM 기법이 난류 영역에서 장시간 통합 동안 분포 함수의 유한한 극값을 유지할 수 있는가?
• RQ4 유한체적 형태의 발산 없는 성질이 불안정성에 대비하여 해를 안정화시키는 데 어떤 역할을 하는가?
• RQ5 다양한 수치적 구현 방식(PGM-DS, PSM-FV, SLS-FV)이 비물리적 잔상 없이 소규모 난류 구조를 얼마나 잘 해석하는가?

주요 결과

• t = 1800에서 PSM 유한체적(FV) 형태는 분포 함수의 초기 극값(min = 0.331, max = 0.4187)을 유지하지만, PSM 방향 분할은 이미 편차를 보이며(min = 0.3086, max = 0.4427).
• t = 4400에서 PSM-DS 기법은 심각한 비물리적 진동을 유발하고 시뮬레이션 충돌 위험이 있으며, 반면 SLS-FV 기법은 허용 가능한 범위 내에서 극값을 유지한다(min = 0.3263, max = 0.4221).
• SLS 제한기는 제어된 확산을 추가함으로써 진동을 효과적으로 억제하여 안정성을 향상시키며, 전체 해 구조에 영향을 주지 않는다.
• 유한체적 형태는 특히 고기울기 및 복잡한 유동 구조가 존재하는 영역에서 방향 분할 대비 더 뛰어난 수치적 안정성을 확보한다.
• BSL 기준 기법 역시 재구성 단계로 인해 최대 원리를 위반하며, 극값 이동을 보인다(min = 0.3124, max = 0.4430 at t = 1800), 그러나 PSM-DS보다 진동 제어에서 더 우수하다.
• 유한체적 기법의 발산 없는 성질은 불안정성을 방지하는 데 핵심적이며, 특히 r-방향 중간 반경 부근에서 SLS-DS는 비물리적 리bbing이 나타나지만 SLS-FV에서는 그러한 현상이 관찰되지 않는다.