[논문 리뷰] Curl-curl conforming elements on tetrahedra

Qian Zhang, Zhimin Zhang|arXiv (Cornell University)|2020. 07. 20.

Advanced Numerical Methods in Computational Mathematics참고 문헌 26인용 수 4

한 줄 요약

이 논문은 3차원 큐드커울 문제를 해결하기 위해 단순체계 메esh에서 H(curl²)-적합한 유한요소를 구축한다. 이는 전자기학 및 자기유체역학에서 나타나는 네째계수 컬 시스템이다. 제안된 요소들은 최적의 보간 오차 추정치를 달성하며, 일반 리프시츠 도메인에서 적합한 유한요소 이산화를 가능하게 하며, 수치적 검증을 통해 정확성과 수렴성을 확인한다.

ABSTRACT

In [24], we proposed H(curl^2)-conforming elements on both a triangle and a rectangle. This family of elements provides a brand new method to solve the quad-curl problem in 2 dimensions. In this paper, we turn our focus to 3 dimensions and construct an H(curl^2)-conforming tetrahedral finite element. The newly proposed element has been proved to have the optimal interpolation error estimate. Having tetrahedral elements, we can solve the quad-curl problem in any Lipschitz domain by conforming finite element method. We also provide several numerical examples of using our element to solve the quad-curl problem. The results of the numerical experiments show the effectiveness and correctness of our element.

연구 동기 및 목표

• 3차원 큐드커울 문제를 해결하기 위해 단순체계 메쉬에서 H(curl²)-적합한 유한요소를 개발한다.
• H¹(curl) 요소보다 더 낮은 정규성 요구 조건을 가진 3차원에서 컬-컬-적합 요소를 구성하는 데 도전하는 문제를 해결한다.
• 최저차수 요소에서 315개의 자유도를 포함하여 매우 큰 자유도 수에도 불구하고 최적의 보간 오차 추정치를 확보한다.
• H²-적합한 방법에서 발생하는 투영 오차를 피하기 위해 정확한 해에 대한 H(curl²) 공간에서 수렴하는 적합한 유한요소 방법을 제공한다.
• 임의의 리프시츠 도메인에서 단순체계 메쉬를 사용하여 큐드커울 문제의 정확한 수치해를 얻을 수 있도록 한다.

제안 방법

• 벡터장의 점 평가, 컬, 변형된 법선 도함수를 모서리와 면에서 포함하는 자유도(DOF)를 정의하여 단순체계에서 H(curl²)-적합한 유한요소를 구성한다.
• 기준 요소와 변환 행렬 D를 사용하여 자유도를 기준 요소에서 물리적 단순체계로 매핑함으로써 유일성과 적합성을 보장한다.
• 최적의 보간 오차 추정치를 증명하기 위해 기준 요소와 관련된 중간 요소를 도입한다.
• 전통적인 라그랑주 유형 기저 함수 계산이 불가능하므로, [8]에서 제시한 방법을 사용하여 고자유도 요소(315개 자유도)의 기저 함수를 계산한다.
• 자기행렬과 모서리/면의 법선 벡터를 포함하는 변환 행렬을 통해 물리적 요소의 모든 자유도를 기준 요소의 기능 값의 선형 조합으로 표현한다.
• 자유도를 기반으로 보간 연산자를 표현하고 일致성과 안정성을 검증함으로써, 보간 오차가 H(curl²) 노름에서 최적으로 수렴함을 증명한다.

실험 결과

연구 질문

• RQ1 3차원 단순체계 메쉬에서 H(curl²)-적합한 유한요소를 구성하여 큐드커울 문제에 대해 최적의 수렴률을 달성할 수 있는가?
• RQ2 자유도가 모서리와 면에서 법선 도함수를 포함할 경우, 이들이 자연스럽게 기준 요소와 물리적 요소 사이로 전달되지 않기 때문에, 컬-컬-적합 요소의 유일성과 적합성을 어떻게 확보할 수 있는가?
• RQ3 제안된 H(curl²)-적합한 유한요소 공간이 단순체계에서의 최적 보간 오차 추정치는 무엇인가?
• RQ4 제안된 방법은 H²-적합한 요소에서 발생하는 투영 오차를 피하면서도 정확한 해에 대해 H(curl²) 공간에서 수렴할 수 있는가?
• RQ5 자유도가 315개에 이르는 고차수 H(curl²)-적합한 요소의 기저 함수를 효율적으로 계산하는 방법은 무엇인가?

주요 결과

• 제안된 단순체계에서의 H(curl²)-적합한 유한요소는 H(curl²) 노름에서 최적의 보간 오차 추정치를 달성하여 이론적 신뢰성을 입증한다.
• 요소의 구성은 수학적으로 엄밀하며, 자유도와 그 변환에 대한 세부 분석을 통해 유일성과 적합성이 검증된다.
• 최저차수 요소는 315개의 자유도를 가지며, 표준 네데lec 요소보다 크게 높아, 컬과 두 번째 도함수를 포괄하기 위한 필요성을 반영한다.
• 수치 실험을 통해 방법의 정확성과 수렴성이 확인되었으며, 3차원 큐드커울 문제에서 최적의 수렴률을 보였다.
• 이 방법은 H²-적합한 방법과 비적합한 방법의 한계를 극복하고, 임의의 리프시츠 도메인에서 큐드커울 문제의 적합한 유한요소 이산화를 가능하게 한다.
• https://github.com/QianZhangMath/3D-curl-curl-conforming-FE에 공개된 코드 기반의 계산 프레임워크를 통해 고자유도 수치 계산이 효율적으로 수행된다.