[논문 리뷰] High-Order Curvilinear Finite Element Magneto-Hydrodynamics I: A Conservative Lagrangian Scheme

Jan Nikl, Milan Kuchařík외 1명|arXiv (Cornell University)|2021. 10. 22.

Computational Fluid Dynamics and Aerodynamics참고 문헌 53인용 수 7

한 줄 요약

이 논문은 곡선 요소를 사용하는 고차수, 보존적인 라그랑주 유한요소 방법을 제안한다. 이는 저항성 유체역학(MHD) 방정식에 적용되며 질량, 운동량, 자기장 선속, 총 에너지의 정확한 보존을 보장하고, 자기장의 수렴이 없는 조건도 유지한다. 방법은 임의의 차수의 공간 수렴성을 달성하며, 안정성을 확보하기 위해 저항성 확산을 암시적으로 처리하는 고차수 시간 적분기를 사용한다.

ABSTRACT

Magneto-hydrodynamics is one of the foremost models in plasma physics with applications in inertial confinement fusion, astrophysics and elsewhere. Advanced numerical methods are needed to get an insight into the complex physical phenomena. The classical Lagrangian methods are typically limited to the low orders of convergence and suffer from violation of the divergence-free condition for magnetic field or conservation of the invariants. This paper is the first part of a new series about high-order non-ideal magneto-hydrodynamics, where a multi-dimensional conservative Lagrangian method based on curvilinear finite elements is presented. The condition on zero divergence of magnetic field and conservation of mass, momentum, magnetic flux and the total energy are satisfied exactly. The curvilinear elements prevent entangling of the computational mesh and its imprinting into the solution. A high-order conservative time integration is applied, where an arbitrary order of convergence is attained for problems of ideal magneto-hydrodynamics. The resistive magnetic field diffusion is solved by an implicit scheme. Description of the method is given and multiple test problems demonstrating properties of the scheme are performed. The construction of the method and possible future directions of development are discussed.

연구 동기 및 목표

• 기존의 저차수 스킴의 한계를 극복하기 위해, 저항성 MHD에 대한 고차수, 보존적인 라그랑주 방법을 개발한다.
• 다차원 시뮬레이션에서 질량, 운동량, 자기장 선속, 총 에너지의 정확한 보존을 보장한다.
• 적절한 변환 규칙을 사용한 모서리 기반 유한요소를 통해 자기장의 수렴이 없는 조건을 정확히 유지한다.
• 곡선 요소에서 고차수 다항식 기저 함수를 사용하여 임의의 차수의 공간 수렴성을 달성한다.
• 자기장 확산의 암시적 시간 적분 기법을 통해 저항성 영역에서의 안정성을 확보한다.

제안 방법

• 곡선 요소 위에서 약한 갈레르킨 형식을 사용하여 RMHD 방정식을 공간 이산화한다.
• 자기장과 전기장에 대해 모서리 기반 유한요소(레비아르트-톰슨 및 네드레크 유형)를 사용하여 ∇·B = 0 조건을 정확히 강제한다.
• 등가형 매핑을 적용하여 곡선 요소를 처리하고 기하학적 일致성을 유지한다.
• 이상적 및 저항성 MHD에 적합한 보존적인 고차수 시간 적분 기법(예: IMEX 및 RK2-평균 방법)을 구현한다.
• 이동하는 라그랑주 기준에서 약한 형식으로부터 유도된 질량, 강성, 및 유량 행렬을 사용한 반연속형 수식을 적용한다.
• 일관된 수치적 유량을 통해 에너지 균형과 파이팅 벡터의 구조를 유지하는 이산 변분 형식을 적용한다.

실험 결과

연구 질문

• RQ1 저항성 MHD에 대해 고차수, 보존적인 라그랑주 유한요소 방법을 구성할 수 있는가? 이는 주요 물리적 보존량을 유지하는가?
• RQ2 고차수 라그랑주 프레임워크에서 자기장의 수렴이 없는 조건을 어떻게 정확히 강제할 수 있는가?
• RQ3 보존성과 안정성을 유지하면서도 임의의 차수의 공간 수렴성을 어느 정도 달성할 수 있는가?
• RQ4 저항성 확산의 암시적 처리 방식이 수치적 안정성과 정확도에 어떤 영향을 미치는가?
• RQ5 비정규 곡선 메esh에서 보존성과 기하학적 일致성을 유지할 수 있는가?

주요 결과

• 이론적으로 이상적 MHD 문제에서 질량, 운동량, 자기장 선속, 총 에너지의 정확한 보존이 달성된다.
• 적절한 변환 규칙을 사용한 모서리 기반 유한요소의 적용으로 인해 자기장의 수렴이 없는 조건이 정확히 만족된다.
• 이상적 MHD에서 이론적으로 임의의 차수의 공간 수렴성이 확보되며, 수치적 검증을 통해 최적의 수렴 속도가 확인된다.
• 저항성 자기장 확산에 대한 암시적 방법은 강한 압축 조건에서도 장시간 시뮬레이션에서 안정성을 보장한다.
• 테스트 문제(예: 자기장 확산, 굴더레이 유형의 유동)에 대한 수치 벤치마크를 통해 방법의 강건성과 정확도가 확인된다.
• PETE2 코드 내에서의 구현은 다물리학 시뮬레이션과의 호환성을 보이며, 비정규 곡선 메쉬에서 고차수 정확도를 지원한다.