[논문 리뷰] An asymptotic-preserving 2D-2P relativistic Drift-Kinetic-Equation solver for runaway electron simulations in axisymmetric tokamaks

Luis Chacòn, Don Daniel외 1명|arXiv (Cornell University)|2021. 05. 04.

Magnetic confinement fusion research참고 문헌 25인용 수 3

한 줄 요약

이 논문은 축대칭 토카막에서 상대론적 드리프트-키네틱 방정식을 위한 점근적 보존, 균일 수렴성 있는 2D-2P 준라그랑주 해소법을 제시한다. 연산자 분할을 사용하여 해석적 라그랑주 궤도 적분과 보존적인 충돌 항 단계를 구현하였으며, 굳은 초구형 운반을 반복하지 않고도 모든 충돌성 범위에서 런어웨이 전자 시뮬레이션을 정확하게 수행한다. 시간에 대해 일차 정확도를, 점근 매개변수 ϵ에 대해 이차 정확도를 확보한다.

ABSTRACT

We propose an asymptotic-preserving (AP), uniformly convergent numerical scheme for the relativistic collisional Drift-Kinetic Equation (rDKE) to simulate runaway electrons in axisymmetric toroidal magnetic field geometries typical of tokamak devices. The approach is derived from an exact Green's function solution with numerical approximations of quantifiable impact, and results in a simple, two-step operator-split algorithm, consisting of a collisional Eulerian step, and a Lagrangian orbit-integration step with analytically prescribed kernels. The AP character of the approach is demonstrated by analysis of the dominant numerical errors, as well as by numerical experiments. We demonstrate the ability of the algorithm to provide accurate answers regardless of plasma collisionality on a circular axisymmetric tokamak geometry.

연구 동기 및 목표

• 축대칭 토카막에서 모든 플라즈마 충돌성 범위에서 정확하고 안정한 상대론적 충돌성 드리프트-키네틱 방정식의 수치 해소법을 개발한다.
• 특히 광범위한 충돌성에서 광범위한 운반을 지배하는 광범위한 충돌성에서의 굳은 초구형 운반 문제를 해결한다.
• 해소법이 점증적 보존성, 균일 수렴성, 확장 가능성을 갖추고, 초구형 성분을 반복적으로 해결하지 않도록 보장한다.
• 런어웨이 전자 생성 및 진화의 정확한 시뮬레이션을 가능하게 하며, 상대론적 충돌 및 에너지 전달 메커니즘을 포함한다.
• 향후 3D 자기 기하구조 및 추가 물리 현상(예: 콩-온 충돌, 복사력 포함)으로의 확장을 위한 기초를 마련한다.

제안 방법

• 해소법은 충돌성 유럽형 단계와 라그랑주 궤도 적분 단계로 분할하여 굳은 운반과 충돌 물리학을 분리한다.
• 초구형 운반 성분은 정확한 그린 함수를 사용하여 재구성되며, 정량적으로 근사된 커널을 사용하여 입자 궤도를 자기선 따라 해석적으로 평가할 수 있도록 한다.
• 라그랑주 단계는 자기장 기하구조에 기반한 해석적 커널을 사용하여 메esh 포인트 수에 비례하는 선형 스케일링과 최적의 성능를 보장한다.
• ϵ → 0 근처에서 느린 다각형에 해석해를 투영하기 위해 라그랑주 커널에 타겟된 소산을 통합하여 점증적 보존성을 확보한다.
• 충돌 단계는 이전 연구에서 유래한 완전히 암시적이고 보존적이며 확장 가능한 상대론적 포커-플랑크 알고리즘을 사용하며, 선형화 및 비선형 모드를 모두 지원한다.
• 각 라그랑주 포인트에서 국소적인 이차 스퍼링 보간을 적용하여 리메핑 과정에서 정확도를 유지한다.

실험 결과

연구 질문

• RQ1 축대칭 토카막에서 모든 충돌성 범위에서 상대론적 드리프트-키네틱 방정식에 대해 점증적 보존성을 갖춘 준라그랑주 스킴을 설계할 수 있는가?
• RQ2 강한 초구형 운반을 반복 보정 없이 해결하면서도, 약한 충돌성 및 강한 충돌성 한계 모두에서 정확성을 유지할 수 있는가?
• RQ3 정의된 바와 같이, ϵ(운동론적 시간스케일 대 충돌성 시간스케일 비율)에 대해 균일 수렴성이 확보되는가?
• RQ4 ϵ ≪ 1 범위에서, 스킴이 빈도 평균화된 포커-플랑크 근사에 얼마나 정확하게 수렴하는가?
• RQ5 점증적 보존성과 효율성 특성을 유지하면서 3D 자기 기하구조로 확장할 수 있는가?

주요 결과

• 수치 실험을 통해 ∆t < ϵ 일 때 시간에 대해 일차 정확도, ∆t > ϵ 일 때 ϵ 에 대해 이차 정확도를 확보하였다.
• 다양한 플라즈마 온도와 충돌성 범위에서 축대칭 토카막 기하구조에서 Dreicer 런어웨이 전자 생성에 대한 이전 결과를 정확히 재현하였다.
• 알고리즘이 점증적 보존성이다: ϵ → 0 근처에서 항상 빈도 평균화된 포커-플랑크 방정식으로 수렴한다.
• 초구형 운반 성분에 대한 반복을 피하여 효율성을 유지하면서도 점증적 성질을 유지한다.
• 라그랑주 단계는 미지수의 수에 비례하여 선형 스케일링을 보이며, 대규모 시뮬레이션에서 최적의 확장 가능성을 확보한다.
• 타겟된 소산을 통해 라그랑주 커널이 자동으로 적절한 시점에 해석해를 느린 다각형으로 투영함으로써, 스킴이 스스로 조정된다.