[논문 리뷰] A DGFEM for Uniformly Elliptic Two Dimensional Oblique Boundary Value Problems

Ellya L. Kawecki|arXiv (Cornell University)|2017. 11. 06.

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

한 줄 요약

이 논문은 곡선 영역에서 비발산형 타입의 2차원 균일 타원형 편미분방정식에 대해 경계 조건이 기울기 있는 경우를 고려한 새로운 비연속 갈레르킨 유한요소법(DGFEM)을 제안한다. 이 방법은 최적 수렴률을 달성하며 경계 조건에 존재하는 일致성 상수를 일致하게 근사함으로써 이전의 딜레르흐 경계 조건에 대한 연구를 더 일반적인 기울기 경계 조건으로 확장하며, 엄밀한 안정성 및 오차 분석을 제공한다.

ABSTRACT

In this paper we present and analyse a discontinuous Galerkin finite element method (DGFEM) for the approximation of solutions to elliptic partial differential equations in nondivergence form, with oblique boundary conditions, on curved domains. In "E. Kawecki, \emph{A DGFEM for Nondivergence Form Elliptic Equations with Cordes Coefficients on Curved Domains}", the author introduced a DGFEM for the approximation of solutions to elliptic partial differential equations in nondivergence form, with Dirichlet boundary conditions. In this paper, we extend the framework further, allowing for the oblique boundary condition. The method also provides an approximation for the constant occurring in the compatibility condition for the elliptic problems that we consider.

연구 동기 및 목표

• 비발산형 타입의 타원형 편미분방정식에서 딜레르흐 경계 조건에서 기울기 경계 조건으로의 비연속 갈레르킨 유한요소법(DGFEM)의 확장을 위해
• 곡선 2차원 영역에서 균일 타원형 방정식을 위한 안정적이고 수렴하는 수치적 스킴을 개발하기 위해
• 비동차 뉴먼형 제약 조건을 수반하는 문제에서 발생하는 기울기 경계 조건 내의 상수를 근사하기 위해
• 다항식 차수와 메쉬 세분화에 대해 최적의 오차 추정치를 확립하기 위해, 심지어 제한된 소볼레프 정규성 조건 하에서도
• 선형화를 통해 몽주–암페르 방정식 및 하미르톤–자코비–벨만 방정식과 같은 비선형 편미분방정식을 해결하기 위한 기초를 마련하기 위해

제안 방법

• 비발산형 타입의 편미분방정식에 특화된 비연속 갈레르킨 수식을 사용하며, 파손된 소볼레프 공간과 요소 간면에서의 수치적 플럭스를 적용한다.
• 곡선 경계를 다항식 매핑을 통해 정확히 표현하기 위해 곡선 삼각형 요소를 사용하여 기하학적 일致성을 확보한다.
• 기울기 벡터장 $\beta$의 방향 도함수를 포함하는 수치적 플럭스를 통해 기울기 경계 조건을 약하게 강제한다.
• 요소 경계에서 기울기의 점프를 제어하고 이산 안정성을 확보하기 위해 안정화 항을 도입한다.
• 기울기 조건을 약하게 강제하는 변분 수식을 통합하여 경계 조건 내 상수의 근사가 가능하도록 한다.
• 이론적 분석은 코르데스 유형 조건과 미르라도–탈렌티 유형 추정치를 기반으로 하여 파손된 $H^2$-준노름에서 최적의 오차 한계를 유도한다.

실험 결과

연구 질문

• RQ1 비발산형 타입의 타원형 편미분방정식에 대해 기울기 경계 조건이 있는 곡선 영역에서 DGFEM을 구성하고 분석할 수 있는가?
• RQ2 제안된 DGFEM이 다항식 차수와 메쉬 세분화에 대해 최적 수렴률을 달성하는가?
• RQ3 유한요소 프레임워크 내에서 기울기 경계 조건 내의 일치성 상수를 일致하게 근사할 수 있는가?
• RQ4 적응형 메쉬 세분화를 사용할 경우, 제한된 소볼레프 정규성(예: $u \in H^{5/2 - \delta}$)이 수렴에 미치는 영향은 무엇인가?
• RQ5 이 방법은 고차원 기울기 경계값 문제로 확장될 수 있으며, 이 경우 이론적 장벽은 무엇인가?

주요 결과

• DGFEM은 다항식 차수 $p$에 대해 파손된 $H^2$-준노름에서 최적 수렴률 $\mathcal{O}(h^p)$을 달성하며, $p=2,3,4$에 대해 수치적으로 확인되었다.
• 해의 정규성이 떨어지는 원점 쪽으로 메쉬를 세분화하는 적응형 메쉬에서, $p=3$ 및 $p=4$일 경우 수렴률이 $\mathcal{O}(\text{ndofs}^{-1/4})$에서 $\mathcal{O}(\text{ndofs}^{-1/2})$로 향상된다.
• $p=2$일 경우 수렴률은 $\mathcal{O}(\text{ndofs}^{-1/4})$로 유지되며, 이는 해의 정규성이 제한되어 있음을 반영한다.
• 이 방법은 기울기 경계 조건 내의 일치성 상수를 일致하게 근사하며, 이는 문제의 잘 정의됨을 보장하는 데 필수적이다.
• 이산 인프-스터프 조건을 통한 안정성 증명을 통해 수치해의 존재성과 유일성이 보장된다.
• 수치 실험을 통해 조건부 정규성 가정 하에서도 $H^1$-준노름과 $H^2$-준노름에서 최적 수렴률이 확인되었다.