[논문 리뷰] A globally conservative finite element MHD code and its application to the study of compact torus formation, levitation and magnetic compression

Carl Dunlea, Ivan Khalzov|arXiv (Cornell University)|2019. 07. 31.

Magnetic confinement fusion research참고 문헌 16인용 수 4

한 줄 요약

이 논문은 선형 삼각형 요소와 행렬 기반 이산 미분 연산자를 사용하여 축대칭 자기유체역학(MHD)을 위한 전역적으로 보존적인 유한요소 코드인 DELiTE를 제시한다. 이 코드는 질량, 에너지, 토로이드 릭스, 각운동량의 전역 보존을 보장한다. 코드는 밀도 높은 토루스 형성, 자기 부상, 압축을 성공적으로 모의하여 실험 진단과 강한 일치를 보이며, 특히 압축 중 릭스 보존과 토로이드 자기장 진화에서 뛰어난 일치를 보인다.

ABSTRACT

The DELiTE (Differential Equations on Linear Triangular Elements) framework was developed for spatial discretisation of partial differential equations on an unstructured triangular grid in axisymmetric geometry. The framework is based on discrete differential operators in matrix form, which are derived using linear finite elements and mimic some of the properties of their continuous counterparts. A single-fluid two-temperature MHD code is implemented in this framework. The inherent properties of the operators are used in the code to ensure global conservation of energy, particle count, toroidal flux, and angular momentum. The code was applied to study a novel experiment in which a compact torus (CT), produced with a magnetized Marshall gun, is magnetically levitated off an insulating wall and then magnetically compressed through the action of currents in the levitation/compression coils located outside the wall. We present numerical models for CT formation, levitation, and magnetic compression, and comparisons between simulated and experimental diagnostics.

연구 동기 및 목표

• 축대칭 시스템에서 질량, 에너지, 토로이드 릭스, 각운동량의 내재된 전역 보존을 갖춘 유한요소 MHD 코드를 개발하기 위해.
• 복잡한 MHD 시뮬레이션에서 수치적 불안정성과 물리적으로 타당하지 않은 해를 해결하기 위해 보존 법칙을 이산 연산자 수준에 통합하기 위해.
• General Fusion에서 수행하는 새로운 실험에서 밀도 높은 토루스(CT)의 형성, 부상, 자기적 압축을 모의 및 시뮬레이션하기 위해.
• 절연벽과 진공장 연계를 통합하여 플라즈마-벽 상호작용과 코ils에 의해 유도되는 자기장 영향을 정확하게 시뮬레이션할 수 있도록 하기 위해.
• 특히 수평 및 토로이드 자기장 측정치를 포함한 실험 진단과의 비교를 통해 수치 모델을 검증하기 위해.

제안 방법

• DELiTE 프레임워크는 비구조적 삼각형 메esh에서 선형 유한요소를 사용하여 행렬 형태의 이산 미분 연산자를 유도한다.
• 이산 연산자는 미메틱 성질(예: 이산 곱의 법칙, 발산 정리)을 만족하도록 구성되어 있어 이산 수준에서 보존 법칙이 보장된다.
• 단일 유체 이온-전자 두 온도 MHD 방정식은 연속 도함수를 행렬 연산자로 대체하여 전역적으로 보존성을 유지한다.
• 절연 영역의 수평 진공 자기장은 별도로 해결하고 플라즈마 영역에 연결하여 토로이드 릭스 보존을 유지한다.
• 명시적 시간 적분을 사용하며, 수치적 안정성을 확보하기 위해 인위적 점성항을 추가한다. 이는 시간 간격과 점성계수의 범위를 제한한다.
• 플라즈마 MHD 해와 외부 코ils 전류 및 진공장 해를 결합하여 자기 부상 및 압축을 모의한다.

실험 결과

연구 질문

• RQ1 축대칭 기하학에서 질량, 에너지, 토로이드 릭스, 각운동량을 전역적으로 보존하는 유한요소 MHD 코드를 설계할 수 있는가?
• RQ2 밀도 높은 토루스 형성, 부상, 자기적 압축 동안 DELiTE 코드는 실험 진단을 얼마나 잘 재현하는가?
• RQ3 절연벽과 진공장 연계를 포함함으로써 플라즈마-벽 상호작용과 릭스 보존에 어떤 영향을 미치는가?
• RQ4 압축 중 시뮬레이션은 실험보다 더 좋은 수평 릭스 보존을 보이는 이유는 무엇이며, 그 차이의 원인은 무엇인가?
• RQ5 이 방법에 포함된 인위적 점성항은 밀도 및 자기장 진단의 정확성에 어떤 영향을 미치는가?

주요 결과

• DELiTE 코드는 미메틱 이산 연산자를 통해 총 에너지, 입자 수, 토로이드 릭스, 각운동량을 전역적으로 보존한다.
• 시뮬레이션은 개선된 코일 구성을 통해 플라즈마-벽 상호작용이 감소함을 보여주며, 이는 압축 중 벽에 대한 접촉 감소로 확인된다.
• 내부 프로브에서의 시뮬레이션 수평 자기장은 압축 기간 동안 계속 상승하며, 릭스 보존과 일치하지만, 실험에서는 약 145 μs에서 붕괴되는 것으로 관측된다.
• 시뮬레이션된 토로이드 자기장 ($B_\phi$)은 실험 평균 $B_\phi$ 측정치와 정성적으로 일치하며, 특히 저유도 경로로 전류가 분배되는 동안 잘 일치한다.
• 모델은 코브라드 샤프트 전류 증가에 따라 $B_\phi$의 일반적인 경향을 잘 포착하지만, 축대칭 모델링은 토로이드적으로 비대칭인 불안정성을 재현할 수 없다.
• 코드의 일阶 정밀도와 명시적 시간 적분은 시뮬레이션 속도와 안정성에 제약을 주며, 특히 높은 점성계수에서는 더욱 두드러진다. 이는 향후 작업에서 암시적 시간 적분이 필요함을 시사한다.