[논문 리뷰] Methods of arbitrary optimal order with tetrahedral finite-element meshes forming polyhedral approximations of curved domains

Vitoriano Ruas|arXiv (Cornell University)|2017. 06. 24.

Advanced Numerical Methods in Computational Mathematics참고 문헌 28인용 수 9

한 줄 요약

이 논문은 3차원 곡면 영역에서 딜레르흐 조건을 갖는 2계 경계값 문제를 풀기 위한 새로운 유한요소법을 제안한다. 직선 모서리를 갖는 테트라헤드론 메esh를 사용하며, 경계 조건을 실제 경계 점에서 보간하는 다항식 시험 함수를 도입함으로써 에너지 노름에서 최적 수렴 속도와 사전에 알려지지 않은 $L^2$-오차 추정치를 달성한다. 이는 곡면 요소나 유리 함수 적분을 요구하지 않으며, 등가적 방법보다 정확도와 계산의 단순성에서 뛰어나다.

ABSTRACT

In recent papers the author introduced a simple alternative to isoparametric finite elements of the n-simplex type, to enhance the accuracy of approximations of second-order boundary value problems with Dirichlet conditions, posed in smooth curved domains. This technique is based upon trial-functions consisting of piecewise polynomials defined on straight-edged triangular or tetrahedral meshes, interpolating the Dirichlet boundary conditions at points of the true boundary. In contrast the test-functions are defined upon the standard degrees of freedom associated with the underlying method for polytopic domains. While method's mathematical analysis for both second- and fourth-order problems in two-dimensional domains was carried out in arxiv NA-1701.00663 and in a submitted paper, this article is devoted to the study of the three-dimensional case, in which the method is nonconforming. Well-posedness, uniform stability and optimal a priori error estimates in the energy norm are demonstrated for a tetrahedron-based Lagrange family of finite elements. Novel L2-error estimates for the class of problems considered in this work are also proved. A series of numerical examples illustrates the potential of the new technique. In particular its better accuracy at equivalent cost as compared to the isoparametric technique is highlighted. Moreover the great generality of the new approach is exemplified through a method with degrees of freedom other than nodal values.

연구 동기 및 목표

• 3차원 곡면 영역에서 고차 근사에 대해 최적 수렴 속도를 유지하는 유한요소법을 개발한다. 이는 등가적 요소에 의존하지 않는다.
• 3D 기하 구조에서 보간된 경계 조건(Implicit Boundary Conditions, IBC)의 제한 사항, 예를 들어 노드 점 선택의 제약과 구현의 복잡성을 극복한다.
• 삼각형 요소의 차수 >1에 대해 잘 정의된 문제, 균일한 안정성, 최적의 사전 오차 추정치를 갖는 수학적으로 타당한 프레임워크를 제공한다.
• 특히 복잡한 비선형성과 메쉬 생성을 다룰 때 등가적 방법에 비해 뛰어난 정확도와 효율성을 보여준다.
• 노드 자유도를 넘어서서도 더 일반적인 적용을 위해 모서리 또는 면 기반 자유도를 포함할 수 있도록 방법을 확장한다.

제안 방법

• 직선 모서리를 갖는 테트라헤드론 메쉬 위에서 정의된 조각다항식 시험 함수를 사용하며, 이는 메쉬의 정점 뿐 아니라 진짜 곡면 경계 상의 점에서도 딜레르흐 경계 조건을 보간한다.
• 테스트 함수는 다각형 영역에 대한 기초 유한요소 공간과 관련된 표준 자유도를 사용하여 구성되며, 고전적 갈레르킨 공식화와 일致성을 유지한다.
• 이 형식은 페트로프-갈레르킨 방법이며, 대칭 문제일지라도 비대칭 행렬을 유도하지만, 최적 수렴 성질을 유지한다.
• 다항식 적분항 덕분에 정확한 수치 적분이 가능하여, 등가적 방법에서 흔히 사용되는 구분 규칙과 유리 함수 적분을 피할 수 있다.
• 이 방법은 임의의 고차 다항식 요소(차수 >1)에 적용 가능하며, 모서리나 면의 평균과 같은 비노드 자유도로 자연스럽게 확장된다.
• 메쉬 생성은 여전히 단순하고 견고하며, 오직 직선 모서리를 갖는 요소만 사용하므로 등가적 요소에서 흔히 발생하는 음수 자코비안 문제를 피할 수 있다.

실험 결과

연구 질문

• RQ1 직선 모서리를 갖는 테트라헤드론 메쉬를 기반으로 하는 유한요소법이 3차원 곡면 영역에서 2계 경계값 문제에 대해 최적 수렴 속도를 달성할 수 있는가?
• RQ2 경계의 정점이 아닌 점에서 딜레르흐 조건이 지정될 경우, 제안된 방법이 최적의 $L^2$-오차 추정치를 유지하는가?
• RQ3 곡면 영역에서 등가적 유한요소 기법에 비해 정확도와 계산 비용 측면에서 이 방법은 어떻게 비교되는가?
• RQ4 이 방법은 노드 자유도를 넘어서 모서리 또는 면 평균과 같은 비노드 자유도로 확장될 수 있으며, 이 경우에도 최적 수렴 성질을 유지하는가?
• RQ5 이 방법은 굵은 메쉬에서도 안정성과 수렴성을 유지하는가? 특히 미세한 메쉬 가정 없이도 안정적인가?

주요 결과

• 삼각형 라그랑주 유한요소의 차수 >1에 대해 표준 정규성 가정 하에 메쉬 크기와 무관하게 에너지 노름에서 최적 수렴이 달성되며, 오차 한계가 독립적임을 입증하였다.
• 고려된 문제의 범주에 대해 사전에 알려지지 않은 $L^2$-오차 추정치가 증명되었으며, 곡면 영역에서 고차 방법의 정확도 기준을 새롭게 설정하였다.
• 수치 예제에서는 동일한 계산 비용에서 등가적 방법보다 정확도가 뛰어나며, 비선형 문제에서 특히 두각을 나타냈다.
• 경계 모서리나 면의 중점에 딜레르흐 조건이 지정되어도 수렴 차수 감소 없이 최적 수렴 성질을 유지한다.
• 다항식 대수의 활용으로 요소 행렬의 정확한 적분이 가능하여, 구분 오차를 피하고 복잡한 비선형성의 구현을 단순화한다.
• 粗한 메쉬에서도 안정적이며, $h$가 영역 지름의 반일지라도 일관된 수렴 행동이 관찰되어, 점근적 가정을 초월한 실용적 적용 가능성을 시사한다.