[논문 리뷰] Divergence-free H(div)-conforming hierarchical bases for magnetohydrodynamics (MHD)

Wei Cai, Jian Wu외 1명|arXiv (Cornell University)|2012. 10. 20.

Electromagnetic Scattering and Analysis참고 문헌 21인용 수 7

한 줄 요약

이 논문은 삼각형 요소의 경우 최고 차수 4, 사면체 요소의 경우 최고 차수 3까지 잘 조절된 질량 행렬과 강성 행렬을 보장하는, 자기유체역학(MHD)을 위한 새로운 계층적 H(div)-일치 유한요소 기반을 제안한다. 이 방법은 각 요소에 고차수 내부 버블 함수를 추가하여 정확한 수렴 자유 자기장 조건을 강제한다. 정규 연속성과 버블 함수의 영이 정규 성분을 통해 국소 및 전역 수렴 자유 조건을 확보하며, 수치 결과는 삼각형의 경우 최고 차수 4, 사면체의 경우 최고 차수 3까지 잘 조절된 질량 행렬과 강성 행렬을 확인한다.

ABSTRACT

In order to solve the magnetohydrodynamics (MHD) equations with a $\mathbf{\mathcal{H}}(\mathbf{div})$-conforming element, a novel approach is proposed to ensure the exact divergence-free condition on the magnetic field. The idea is to add on each element an extra interior bubble function from higher order hierarchical $\mathbf{\mathcal{H}}(\mathbf{div})$-conforming basis. Four such hierarchical bases for the $\mathbf{\mathcal{H}}% (\mathbf{div})$-conforming quadrilateral, triangular, hexahedral and tetrahedral elements are either proposed (in the case of tetrahedral) or reviewed. Numerical results have been presented to show the linear independence of the basis functions for the two simplicial elements. Good matrix conditioning has been confirmed numerically up to the fourth order for the triangular element and up to the third order for the tetrahedral element.

연구 동기 및 목표

• 정확한 수렴 자유 자기장을 갖는 MHD 방정식을 해결하기 위한 강력한 H(div)-일치 유한요소 방법을 개발한다.
• 이전의 계층적 H(div) 기반에서 2차 근사에서 선형 종속성 문제로 인해 실패하는 문제를 해결한다.
• 요소 경계에서 사라지는 내부 버블 함수를 통해 국소 수렴 자유 조건을 강제한다.
• 단순체 요소에 대해 제안된 기반 함수의 선형 독립성과 양호한 행렬 조건을 수치적으로 입증한다.
• H(div)-일치 요소를 사용한 고차수, 안정적이고 정확한 MHD 시뮬레이션을 위한 기초를 제공한다.

제안 방법

• 기준 요소 프레임워크를 사용하여 사각형, 삼각형, 육각형, 사면체 요소에 대한 새로운 계층적 H(div)-일치 기반을 도입한다.
• 모서리 및 면 기반 자유도를 사용하고, 일반 메esh에서 정규 연속성을 보장하기 위해 Piola 변환을 적용한다.
• p차수 H(div) 기반에 (p+1)차수 내부 버블 함수를 추가하여 각 요소에서 국소 수렴 자유 조건을 강제한다.
• 버블 함수는 모든 요소 경계에서 정규 성분이 0이므로 정규 연속성을 유지하고 전역 수렴 자유 자기장을 가능하게 한다.
• 노드, 모서리, 면 기반 함수를 포함한 계층적 구성 방식을 사용하며, 탄성 성분을 위한 고차수 모멘트도 포함한다.
• 질량 행렬과 강성 행렬을 사용하여 조건 수를 수치적으로 평가하며, 조건 수는 최대 고유값과 최소 고유값의 비율로 계산된다.

실험 결과

연구 질문

• RQ1 이전 수식에서 발생하는 선형 종속성 문제를 피하면서도 사면체 요소에 대해 계층적 H(div)-일치 기반을 구성할 수 있는가?
• RQ2 요소 수준에서 자기장의 수렴 자유 조건을 정확히 강제하면서도 정규 연속성을 유지할 수 있는가?
• RQ3 고차수 근사에 대해 결과로 얻어진 질량 행렬과 강성 행렬의 조건 수 행동은 어떠한가?
• RQ4 고차수 내부 버블 함수의 포함이 기반 함수의 선형 독립성과 안정성을 유지하는가?
• RQ5 제안된 방법이 실질적인 MHD 시뮬레이션에서 삼각형의 경우 p=4, 사면체의 경우 p=3까지 잘 조절된 행렬을 달성할 수 있는가?

주요 결과

• 제안된 사면체 요소에 대한 계층적 H(div) 기반은 질량 행렬의 유한 조건 수(최고 차수 p=4까지)로 인해 선형 독립임이 확인되었다.
• 삼각형 요소의 경우, p=4일 때 질량 행렬의 조건 수는 1.286×10⁴, 강성 행렬의 조건 수는 8.843×10³로 유지되어 양호한 조건 수를 보였다.
• 사면체 요소의 경우, 질량 행렬의 조건 수는 3.412×10⁶, 강성 행렬의 조건 수는 1.094×10⁶로 p=3까지 잘 조절된 상태였다.
• p=2일 때, 첫 번째 종류의 모서리 기반 면 기반 기반은 질량 조건 수 6.987×10³으로 두 번째 종류의 7.733×10⁴보다 더 양호한 조건 수를 보였다.
• p=4일 때, 두 번째 종류의 모서리 기반 면 기반 기반은 첫 번째 종류의 5.972×10⁹보다 우수한 조건 수(질량 조건 수 2.717×10⁷)를 보였다.
• 정규 연속성의 p차수 기반 함수와 요소 경계에서 정규 흐름이 0인 (p+1)차수 내부 버블 함수를 조합함으로써 전역 수렴 자유 자기장을 확보한다.