[논문 리뷰] Construction and Analysis of an HDG Method for Incompressible Magnetohydrodynamics

Jeonghun J. Lee, Stephen Shannon외 2명|arXiv (Cornell University)|2017. 02. 16.

Advanced Numerical Methods in Computational Mathematics참고 문헌 27인용 수 2

한 줄 요약

이 논문은 선형화된 비압축성 자기유체역학(MHD) 방정식을 위한 하이브리드 불연속 갈레르킨(HDG) 방법을 제안한다. 이 방법은 이중 사다리꼴 형태에서 유도된 새로운 HDG 흐름을 사용하며, 유체 속도 및 자석 변수에 대해 최적 수렴성을 확보하고, 다른 변수에 대해서는 준최적 수렴성을 달성한다. 이는 2차원 및 3차원의 단체 메esh에서 사전 오차 분석과 수치 실험을 통해 검증되었다.

ABSTRACT

We present a hybridized discontinuous Galerkin (HDG) method for stationary linearized incompressible magnetohydrodynamics (MHD) equations. At the heart of the paper is the introduction of an HDG flux of the dual saddle-point form of the MHD equations that facilitates the hybridization of discontinuous Galerkin (DG) method. We carry out the $ extit{a priori}$ error estimates for the proposed HDG method on simplicial meshes in both two- and three-dimensions. The analysis provides optimal convergence for the fluid velocity and the magnetic variables, and quasi-optimal convergence for the remaining quantities. Numerical examples are presented to verify the theoretical findings.

연구 동기 및 목표

• 정적인 선형화된 비압축성 자기유체역학(MHD) 방정식을 해결하기 위한 안정적이고 정확한 수치 방법을 개발한다.
• 이중 사다리꼴 형태의 설정에 새로운 HDG 흐름을 도입하여 하이브리드 불연속 갈레르킨(HDG) 프레임워크를 MHD 시스템으로 확장한다.
• 이중 및 삼중 차원의 단체 메esh에서 제안된 HDG 방법에 대한 엄밀한 사전 오차 추정을 수립한다.
• 2차원 및 3차원 영역에서의 수치 실험을 통해 이론적 수렴 속도를 검증한다.

제안 방법

• 방법은 선형화된 MHD 방정식의 이중 사다리꼴 형태에서 유도된 하이브리드 불연속 갈레르킨(HDG) 설정을 사용한다.
• 하이브리드화를 가능하게 하기 위해 새로운 HDG 수치적 흐름을 도입하여 전역 자유도를 감소시키면서도 안정성과 정확성을 유지한다.
• 이 설정은 이중 및 삼중 차원의 단체 메쉬에 적용되어 유연하고 국소 질량 보존 성질을 갖는 근사치를 가능하게 한다.
• 요소 경계에서의 추적 불확실성에 의해 전역 시스템이 분리되어 효율적인 해법 전략을 구현할 수 있다.
• 함수 해석학 및 일치하는 단체 메쉬 삼각분할에 대한 근사 이론을 사용하여 사전 오차 추정을 유도한다.

실험 결과

연구 질문

• RQ1 선형화된 비압축성 MHD 방정식에 대해 하이브리드 불연속 갈레르킨 방법을 성공적으로 설정할 수 있는가?
• RQ2 HDG 방법이 이중 및 삼중 차원 설정에서 유체 속도, 자기장 및 기타 변수에 대해 수렴 속도는 어떻게 되는가?
• RQ3 제안된 HDG 흐름이 모든 해 성분에 대해 최적 또는 준최적 수렴성을 보장하는가?
• RQ4 이론적 오차 추정치는 단체 메쉬에서의 수치 결과와 어떻게 비교되는가?

주요 결과

• HDG 방법은 이중 및 삼중 차원에서 모두 유체 속도 및 자기장 변수에 대해 최적 수렴 속도를 달성한다.
• 압력 및 자기포텐셜 변수에 대해 이론적 분석에서 예측한 바와 같이 준최적 수렴이 확보된다.
• 사전 오차 추정치는 엄밀히 유도되었으며, 단체 메쉬에서의 수치 실험을 통해 검증되었다.
• 수치 예제를 통해 다양한 다항식 차수와 메쉬 유형에서 이론적 수렴 행동이 확인되었다.