[논문 리뷰] A Nonconforming Finite Element Method for Fourth Order Curl Equations in R^3

Bin Zheng, Qiya Hu외 1명|arXiv (Cornell University)|2010. 01. 31.

Advanced Numerical Methods in Computational Mathematics참고 문헌 15인용 수 3

한 줄 요약

이 논문은 자기유체역학(MHD) 모델에서 발생하는 3차원에서의 네째단계 컬 방정식을 해결하기 위한 비일치 유한요소법을 제시한다. 이 방법은 컬의 접선 연속성을 강제하는 최소 자유도를 가진 새로운 요소를 사용하여 최적 수렴 속도를 달성하고, 구현을 위한 명시적 기저 함수를 제공한다.

ABSTRACT

In this paper we present a nonconforming finite element method for solving fourth order curl equations in three dimensions arising from magnetohydrodynamics models. We show that the method has an optimal error estimate for a model problem involving both curl^2 and curl^4 operators. The element has a very small number of degrees of freedom and it imposes the inter-element continuity along the tangential direction which is appropriate for the approximation of magnetic fields. We also provide explicit formulae of basis functions for this element.

연구 동기 및 목표

• 3차원 자기유체역학(MHD) 문제에서 네째단계 컬 연산자를 안정적이고 효율적으로 이산화하는 유한요소법을 개발하기 위해.
• 표준 일치 요소로는 고 자유도로 인해 어려운 MHD 모델에서 $(\nabla\times)^4$ 연산자를 이산화하는 데 도전하기 위해.
• 요소 간 컬의 접선 연속성을 보장하는 비일치 유한요소를 구성하여 자석장 근사에 적합하게 하기 위해.
• 실제 구현을 가능하게 하기 위해 새로운 요소에 대한 명시적 기저 함수를 제공하기 위해.
• 기본 소볼레프 정규성 조건 하에 이 방법의 최적 오차 추정을 수립하기 위해.

제안 방법

• 조각별 2차 다항식에 기반한 비일치 유한요소 공간을 사용하며, 컬과 관련된 자유도를 보강한다.
• 유한요소 공간은 컬의 요소 간 접선 성분의 연속성을 강제하여 자석장 물리법칙과 일致하도록 한다.
• 체적 및 면적 적분을 포함하는 $ \nabla\times\mathbf{v}_h $ 를 이용한 이차형식 $ a_h(\cdot, \cdot) $ 를 구성하며, $ \alpha $ 와 $ \beta $ 를 통한 안정화 항을 포함한다.
• 비일치 갈레르킨 설정을 사용하여 분석을 수행하며, 부분적 적분과 경계 추정식을 통해 일치성 및 근사 오차 추정을 유도한다.
• 2차 스트랑의 보조정리를 활용하여 컬과 그 기울기의 $ H^1 $-유사 노름에서 최적 수렴을 확립한다.
• 요소에 대한 명시적 기저 함수를 유도하고 제공하여 유한요소 코드에서 직접 구현 가능하게 한다.

실험 결과

연구 질문

• RQ1 최소 자유도로 3차원에서 네째단계 컬 연산자를 효율적으로 이산화할 수 있는 비일치 유한요소법을 구성할 수 있는가?
• RQ2 요소 경계에서 컬의 접선 연속성을 강제하면 네째단계 컬 문제에서 최적 수렴이 이루어지는가?
• RQ3 MHD 모델의 $(\nabla\times)^4$ 연산자는 비일치 접근을 통해 어떻게 안정적으로 근사할 수 있는가?
• RQ4 기본 정규성 조건 하에 이러한 방법으로 달성 가능한 최적 오차 추정은 무엇인가?
• RQ5 실제 구현을 지원하기 위해 새로운 요소에 대한 명시적 기저 함수를 도출할 수 있는가?

주요 결과

• 제안된 비일치 유한요소법은 에너지 노름에서 최적 수렴 속도를 달성하며, $ \mathbf{u} \in (H^4(\Omega))^3 $ 에 대해 오차가 $ \lesssim h\|\mathbf{u}\|_{4,\Omega} $ 로 유계이다.
• 해의 해, 그 컬, 그리고 컬의 기울기 모두 최적 오차 추정을 보이며 고차 정밀도를 확보한다.
• 각 삼각형 요소당 자유도가 30개 뿐이므로, 기존 일치 방법이 요구하는 220개의 자유도에 비해 계산 비용을 크게 감소시킨다.
• 요소 간면에서 컬의 접선 연속성을 강제함으로써 자석장의 물리적 구조를 유지한다.
• 기저 함수에 대한 명시적 공식을 도출하여 유한요소 소프트웨어에서 직접 구현이 가능하다.
• 일치 오차는 $ \lesssim h(\|\nabla\times\Delta(\nabla\times\mathbf{u})\| + \|\nabla(\nabla\times\mathbf{u})\|_{1} + \|\nabla\times\mathbf{u}\| + \|\nabla\times\nabla\times\mathbf{u}\|) $ 로 유계이며, 최적 근사 성질을 확인한다.