[논문 리뷰] Gaussian-Moment Relaxation Closures for Verifiable Numerical Simulation of Fast Magnetic Reconnection in Plasma

Evan Johnson|arXiv (Cornell University)|2014. 09. 24.

Gas Dynamics and Kinetic Theory참고 문헌 32인용 수 15

한 줄 요약

이 논문은 빠른 자기재결합의 검증 가능하고 수렴하는 시뮬레이션을 가능하게 하며, 운동학적 시뮬레이션과 일치하는 두상태 자기유체역학(MHD)을 위한 가우시안-모멘트 이완 폐쇄를 제안한다. 압력 텐서의 등방성화와 열확산계수에 대한 가우시안-BGK 폐쇄를 통합함으로써, 모델은 특이성을 피하고 확산 영역 내 핵심 물리 현상을 포착한다. 수치 결과는 최고점에서의 자기재결합과 양호한 일치를 보이며, 명시적 충돌 항목이 없을 경우 불안정성을 보인다.

ABSTRACT

The motivating question for this dissertation was to identify the minimal requirements for fluid models of plasma to allow converged simulations that agree well with converged kinetic simulations of fast magnetic reconnection. We show that truncation closure for the deviatoric pressure or for the heat flux results in singularities. Due to the strong pressure anisotropies that arise, we study magnetic reconnection with a Gaussian-moment two-fluid MHD with isotropization of the pressure tensor. For the GEM magnetic reconnection challenge problem, our deviatoric pressure tensor agrees well with published kinetic simulations at the time of peak reconnection, but sometime thereafter the numerical solution becomes unpredictable and develops near-singularities that crash the simulation unless positivity limiters are applied. To explain these difficulties, we show that steady reconnection requires heat flux. Specifically, for two-dimensional problems invariant under 180-degree rotation about the X-point, entropy production in the vicinity of the X-point is necessary for reconnection and causes solutions to become singular for conservative models that lack heat flux. This prompts the need for a 10-moment gyrotropic heat flux closure. Using a Chapman-Enskog expansion with a Gaussian-BGK collision operator yields a heat flux closure for a magnetized 10-moment charged gas which generalizes the closure of McDonald and Groth.

연구 동기 및 목표

• 빠른 자기재결합의 수렴 가능한 시뮬레이션을 운동학적 결과와 일치시키기 위한 최소한의 유체 모델 요구 조건을 규명하는 것.
• 자기재결합에서 강한 압력 이방성으로 인해 발생하는 잘린 유체 모델에서 특이성이 발생하는 문제를 해결하는 것.
• 이sovotropic 다섯 또는 열모멘트 모델이 명시적 충돌 연산자가 없이도 수렴 가능하고 특이성이 없는 정상상태 자기재결합을 지지할 수 없다는 것을 입증하는 것.
• 확산 영역 내 필수 운동학 물리 현상을 포착하는 최소한의 물리적으로 일관된 모델을 개발하는 것.
• GEM 자기재결합 도전 과제 문제와의 비교를 통해 모델을 검증하고 수치적 안정성 및 수렴성을 분석하는 것.

제안 방법

• 특이성을 방지하기 위해 압력 텐서의 등방성화를 포함한 가우시안-모멘트 두상태 MHD 모델을 도입한다.
• 가우시안-BGK 충돌 연산자를 적용하여 맥Donald과 그로스의 폐쇄를 일반화한 10모멘트 열확산계수 폐쇄를 유도한다.
• 자기장이 있는 플라즈마에서 BGK 충돌 모델로부터 체르노프-엔스코그 전개를 사용해 열확산계수 폐쇄를 유도한다.
• 등온 두상태-맥스웰 시스템과 등방성화를 포함한 GEM 자기재결합 도전 과제의 시뮬레이션에 모델을 구현한다.
• 조정 가능한 등방성화 시간 τ_s를 가진 내성 충돌 항목(R_s 및 Q_s)을 통합하여 수렴성을 연구한다.
• 열확산계수 항목이 생략되었을 경우 시뮬레이션을 안정화시키기 위해 정규성 제한기(positivity limiters)를 적용하여 수치적 불안정성을 드러낸다.

실험 결과

연구 질문

• RQ1 빠른 자기재결합의 수렴 가능하고 특이성이 없는 시뮬레이션을 운동학적 결과와 일치시키기 위해 필요한 최소한의 유체 모델 구조는 무엇인가?
• RQ2 왜 잘린 유체 모델(예: 다섯 또는 열모멘트)이 등온 폐쇄를 사용할 경우 정상상태 자기재결합을 수렴 가능하고 특이성이 없는 상태로 지지할 수 없는가?
• RQ3 열확산계수는 정상 자기재결합을 유지하는 데 어떤 역할을 하는가? 왜 수치적 안정성에 필수적인가?
• RQ4 가우시안-BGK 기반 열확산계수 폐쇄는 엔트로피를 존중하는 폐쇄에 비해 물리적으로 일관되고 수치적으로 안정적인 대안이 될 수 있는가?
• RQ5 명시적 충돌 항목(등방성화)의 포함이 두상태 MHD 시뮬레이션의 수렴성과 정확성에 어떤 영향을 미치는가?

주요 결과

• 편차 압력 또는 열확산계수 텐서에 대한 잘린 폐쇄는 정상상태 자기재결합 시뮬레이션에서 특이성을 유도한다.
• 등온 다섯모멘트 또는 열모멘트 모델은 명시적 충돌 연산자가 없이도 수렴 가능하고 특이성이 없는 정상상태 자기재결합을 지지할 수 없다.
• 압력의 등방성화와 가우시안-BGK 열확산계수 폐쇄를 포함한 제안된 가우시안-모멘트 두상태 MHD 모델은 피크 자기재결합에서 운동학적 시뮬레이션과 양호한 일치를 보이는 편차 압력을 생성한다.
• 명시적 열확산계수 항목이 없는 시뮬레이션은 불안정해지고 근접 특이성을 띠며, 시스템 충돌을 방지하기 위해 정규성 제한기가 필요하다.
• 자기재결합에서 지속적인 열확산계수의 필요성은 동적 자기재결합이 계속적으로 주입되는 과정을 모의하는 것으로 해석되며, 이는 10모멘트 기울기 폐쇄가 필요하다는 것을 의미한다.
• 유도된 가우시안-BGK 열확산계수 폐쇄는 맥Donald과 그로스의 폐쇄를 일반화하며, 이 맥락에서 엔트로피를 존중하는 대안보다 더 물리적으로 일관된 것으로 입증된다.