[논문 리뷰] Structure-preserving Finite Element Methods for Stationary MHD Models

Kaibo Hu, Jinchao Xu|arXiv (Cornell University)|2015. 03. 20.

Advanced Numerical Methods in Computational Mathematics참고 문헌 13인용 수 5

한 줄 요약

이 논문은 정적인 자기유체역학(MHD)에 대해 자기장 $\bm{B}$와 전류 밀도 $\bm{j}$를 주요 변수로 사용하는 구조 보존 혼합 유한요소 방법을 제안한다. 이 방법은 $H(\mathrm{div})$-적합한 유한요소와 라그랑주 승수를 사용하여 발산이 없는 조건 $\nabla\cdot\bm{B}=0$와 시스템의 에너지 법칙을 정확히 유지함으로써, 레이놀즈 수에 제한 없이 잘 정의되고 수렴하는 결과를 보장한다.

ABSTRACT

In this paper, we develop a class of mixed finite element scheme for stationary magnetohydrodynamics (MHD) models, using magnetic field $\bm B$ and current density $\bm j$ as the discretization variables. We show that the Gauss's law for the magnetic field, namely $ abla\cdot\bm{B}=0$, and the energy law for the entire system are exactly preserved in the finite element schemes. Based on some new basic estimates for $H^{h}(\mathrm{div})$, we show that the new finite element scheme is well-posed. Furthermore, we show the existence of solutions to the nonlinear problems and the convergence of Picard iterations and finite element methods under some conditions.

연구 동기 및 목표

• 정적인 MHD를 위한 유한요소 스킴을 개발하여 이산 수준에서 자기장에 대한 가우스 법칙 $\nabla\cdot\bm{B}=0$를 정확히 유지한다.
• 전체 MHD 시스템의 이산 에너지 법칙을 정확히 유지하여 수치 시뮬레이션에서 물리적 일관성을 유지한다.
• 소규모 레이놀즈 수 조건이 필요 없이 비선형 이산 문제의 잘 정의됨을 확립한다. 이는 이전의 수식화 방식의 한계를 극복한다.
• 시간에 의존하는 MHD에서의 구조 보존 기법을 정적인 경우로 확장한다. 여기서 패러데이 법칙이 자동으로 발산이 없는 $\bm{B}$를 보장하지 않기 때문이다.
• 적절한 정규성 및 유한요소 근사 조건 하에서 피카르 반복과 유한요소 근사의 수렴성을 보여준다.

제안 방법

• 이 방법은 $\bm{B}$와 $\bm{j}$를 주요 미지수로 사용하며, $\|\bm{j}\|$로 표현된 에너지 표현을 고려하여 $\bm{E}$보다 이들을 선택하는 데 동기가 된다.
• 이산 수준에서 $\nabla\cdot\bm{B}=0$를 강제로 유지하기 위해 라그랑주 승수를 도입하여 가우스 법칙의 정확한 유지 보장을 한다.
• 이산 수식은 $\bm{B}$에 대해 $H(\mathrm{div})$-적합한 유한요소와 관련 변수에 대해 $H(\mathrm{curl})$-적합한 유한요소를 사용하여, 드 라함 복합체와의 호환성을 확보한다.
• 로렌츠 힘 항은 $(\bm{E} + \mathbb{P}(\bm{u}\times\bm{B}), \mathbb{P}(\bm{v}\times\bm{B}))$로 이산화되며, 여기서 $\mathbb{P}$는 사영 연산자이다. 이는 시간에 의존하는 방법과 다름을 보이며, 소규모 $R_e$ 조건 없이 잘 정의됨을 가능하게 한다.
• 이산 인프-스터프 조건과 방법의 안정성 분석을 위해 $H^h(\mathrm{div})$ 공간에 대한 새로운 추정식을 개발한다.
• 속도, 압력, 자기장, 전류 밀도, 라그랑주 승수에 적절한 함수 공간을 사용한 혼합 변분 수식을 통해 방법을 분석한다.

실험 결과

연구 질문

• RQ1 정적인 MHD에 대해 이산 수준에서 $\nabla\cdot\bm{B}=0$를 정확히 유지할 수 있는 유한요소 방법을 구성할 수 있는가?
• RQ2 정적인 유한요소 이산화에서 MHD 시스템의 에너지 법칙을 어떻게 정확히 유지할 수 있는가?
• RQ3 레이놀즈 수 조건을 소규모로 제한하지 않고도 이산 문제의 잘 정의됨을 확립할 수 있는가?
• RQ4 변수 선택($\bm{B}$와 $\bm{j}$)이 구조 보존과 제한적인 가정을 피하는 데 어떤 역할을 하는가?
• RQ5 특히 로렌츠 힘 항의 이산 비선형 항은 방법의 안정성과 수렴성에 어떤 영향을 미치는가?

주요 결과

• 제안된 유한요소 스킴은 $H(\mathrm{div})$-적합한 요소와 라그랑주 승수를 통해 이산 수준에서 가우스 법칙 $\nabla\cdot\bm{B}=0$를 정확히 유지한다.
• 전체 MHD 시스템의 이산 에너지 법칙이 정확히 유지되어 장기적 안정성과 물리적 일관성을 확보한다.
• 레이놀즈 수 $R_e$에 대한 제한 없이 이산 문제의 잘 정의됨이 입증되었으며, $\bm{E}$를 기반으로 한 수식화와는 달리 이는 이전의 제약을 피한다.
• 적절한 정규성 및 근사 조건 하에서 피카르 반복과 유한요소 근사의 수렴성이 증명되었다.
• 유한요소 공간이 $r_1 = r_2 = r_3 + 1 = s_u = s_p + 1$ 조건을 만족할 경우 최적 수렴률을 보여주는 오차 추정식을 유도하였다. 여기서 BDM와 네데일 요소는 타당한 선택이다.
• 분석의 기술적 한계로 인해 최저 순서의 라비아르트-톰슨 요소를 사용하지 않으며, 이 제한은 본질적이지 않다고 여겨진다.