[논문 리뷰] Construction and application of algebraic dual polynomial representations for finite element methods

Varun Jain, Yi Zhang외 2명|arXiv (Cornell University)|2017. 12. 27.

Advanced Numerical Methods in Computational Mathematics인용 수 4

한 줄 요약

이 논문은 유한요소법에서 대수적 쌍대 다항식 기저를 체계적으로 구성하는 방법을 제안하며, 딜리클레 문제와 뉴먼 문제에 대한 정확한 혼합형식을 가능하게 한다. 이는 유한요소 근사가 어떤 격자에서도 항등식 $\phi^h = \text{div}\, \mathbf{q}^h$ 를 만족하며, 제약 조건과 라그랑주 승수의 이산화가 격자 크기, 형태, 다항식 차수와 무관하게 이루어진다는 것을 보여준다.

ABSTRACT

Given a polynomial basis $\Psi_i$ which spans the polynomial vector space $\mathcal{P}$, this paper addresses the construction and use of the algebraic dual space $\mathcal{P}'$ and its canonical basis. Differentiation of dual variables will be discussed. The method will be applied to a Dirichlet and Neumann problem presented in \cite{CarstensenDemkowiczGopalakrishnan} and it is shown that the finite dimensional approximations satisfy $\phi^h = \mbox{div}\, \mathbf{q}^h$ on any grid. The dual method is also applied to a constrained minimization problem, which leads to a mixed finite element formulation. The discretization of the constraint and the Lagrange multiplier will be independent of the grid size, grid shape and the polynomial degree of the basis functions.

연구 동기 및 목표

• 유한요소법에서 주어진 다항식 기저 $\Psi_i$ 에 대응하는 대수적 쌍대 공간 $\mathcal{P}'$ 의 체계적 구성 방법을 개발하는 것.
• 모든 메esh에서 이산 해가 $\phi^h = \text{div}\, \mathbf{q}^h$ 를 만족하도록 보장하여 정확한 혼합 유한요소 형식을 가능하게 하는 것.
• 제약 조건과 라그랑주 승수의 이산화를 격자 기하학성과 다항식 차수에서 분리하여 강건성을 향상시키는 것.
• 쌍대 변수의 미분을 위한 프레임워크를 제공하여 혼합 방법에서 고차 정확도를 지원하는 것.

제안 방법

• 다항식 공간 $\mathcal{P}$ 에 대해 $\Psi_i$ 가 생성하는 기저 $\mathcal{P}'$ 에 대응하는 대수적 쌍대 기저를 구성하고, 이중수직성 조건을 보장하는 것.
• 특히 딜리클레 및 뉴먼 경계 조건을 가진 문제에 대해 쌍대 기저를 혼합 유한요소 형식에 적용하는 것.
• 약약한 형식에서 일致성을 유지하기 위해 쌍대 변수의 미분 규칙을 유도하고 적용하는 것.
• 쌍대 공간을 사용하여 제약 조건이 포함된 최소화 문제를 구현하고, 제약 조건 이산화에 독립적인 혼합 형식을 도출하는 것.
• 제약 조건과 관련된 라그랑주 승수의 이산화가 격자 크기와 형태의 변화에 영향을 받지 않도록 보장하는 것.
• 유한요소 근사가 메쉬의 구조와 관계없이 어떤 격자에서도 항등적으로 $\phi^h = \text{div}\, \mathbf{q}^h$ 를 만족함을 보여주는 것.

실험 결과

연구 질문

• RQ1 유한요소법에서 주어진 다항식 기저 $\Psi_i$ 에 대해 대수적 쌍대 다항식 기저를 어떻게 체계적으로 구성할 수 있는가?
• RQ2 쌍대 기저 형식이 임의의 비정규 격자에서 $\phi^h = \text{div}\, \mathbf{q}^h$ 가 정확히 성립하도록 보장할 수 있는가?
• RQ3 제약 조건과 라그랑주 승수의 이산화가 격자 크기와 형태로부터 얼마나 독립적으로 만들 수 있는가?
• RQ4 혼합 유한요소 체계에서 쌍대 변수의 미분을 어떻게 일致성 있게 정의하고 적용할 수 있는가?
• RQ5 혼합 유한요소법에서 제약 최소화 문제에 대해 쌍대 다항식 기저를 사용할 경우 어떤 영향을 미치는가?

주요 결과

• 구성된 대수적 쌍대 기저는 메쉬의 구조와 관계없이 어떤 격자에서도 유한요소 해가 $\phi^h = \text{div}\, \mathbf{q}^h$ 를 항등적으로 만족함을 보장한다.
• 이 방법은 제약 조건과 라그랑주 승수의 이산화가 격자 크기, 형태, 다항식 차수와 무관하게 불변함을 달성한다.
• 쌍대 변수의 미분은 일致성 있게 정의되어 변분 형식에 통합되어 정확도를 유지한다.
• 쌍대 기저에서 유도된 혼합 유한요소 형식은 주어진 조건 하에서 최적 수렴 성질을 유지한다.
• 이 방법은 동일한 쌍대 프레임워크를 사용하여 딜리클레 문제와 뉴먼 문제 모두에 대해 강건하고 일致한 혼합 형식을 가능하게 한다.
• 이 프레임워크는 안정적이고 독립적인 제약 조건 이산화를 지원하는 제약 최소화 문제를 처리할 수 있어 수치적 신뢰성을 향상시킨다.