[논문 리뷰] An unfitted $hp$-interface penalty finite element method for elliptic interface problems

Haijun Wu, Yuanming Xiao|arXiv (Cornell University)|2010. 07. 17.

Advanced Numerical Methods in Computational MathematicsEngineering참고 문헌 29인용 수 17

한 줄 요약

이 논문은 2차원 및 3차원에서 메esh가 맞지 않는 경우에도 타원형 인터페이스 문제를 위해 $hp$-버전의 맞지 않는 인터페이스 페널티 유한요소법($hp$-IPFEM)을 제안한다. 이는 인터페이스 유량 안정화와 페널티 항을 포함한 대칭 및 비대칭 변분 형식을 사용하여, 파손된 $H^1$ 노름에서 최적의 $h$-수렴과 반차수만큼 열등한 $p$-수렴(하나의 차수만큼 낮음)을 달성하며, 쌍대성 원리를 통해 $L^2$ 오차 추정치를 도출한다.

ABSTRACT

An $hp$ version of interface penalty finite element method ($hp$-IPFEM) is proposed for elliptic interface problems in two and three dimensions on unfitted meshes. Error estimates in broken $H^1$ norm, which are optimal with respect to $h$ and suboptimal with respect to $p$ by half an order of $p$, are derived. Both symmetric and non-symmetric IPFEM are considered. Error estimates in $L^2$ norm are proved by the duality argument.

연구 동기 및 목표

• 메쉬에 맞는 삼각형 메쉬가 필요 없이, 인터페이스가 요소를 가로지르는 메쉬에서 고차수이고 안정적인 유한요소법을 개발하기 위해.
• 메쉬 크기 $h$에 대한 최적의 수렴률과 다항식 차수 $p$에 대한 열등하지만 날카로운 수렴률을 파손된 $H^1$ 노름에서 달성하기 위해.
• 페널티 유한요소법을 $hp$-버전 프레임워크로 확장하여, 대칭 및 비대칭 형식 모두에서 안정성과 수렴성을 보장하기 위해.
• 쌍대성 원리를 사용하여 $H^1$ 및 $L^2$ 노름에서 엄밀한 오차 추정치를 도출하여, 인터페이스 위치에 대해 강건한 성능을 확보하기 위해.

제안 방법

• 해의 값과 유량의 점프를 페널티화하는 변분 형식을 사용하여 대칭 및 비대칭 $hp$-인터페이스 페널티 유한요소법($hp$-IPFEM)을 수립한다.
• 메쉬가 맞지 않는 경우의 정확도와 조건수 개선을 위해, 이차형식에 가중 평균 유량 기법을 적용한다.
• 해의 값과 유량의 점프에 대한 페널티 항을 도입하여 인터페이스 조건을 약한 형태로 강제함으로써 안정성과 수렴성을 확보한다.
• 불연속 갈레르킨 유사 접근법을 사용하며, 인터페이스에 의해 잘린 요소에서 고차수 다항식 근사($p$-refinement)를 허용한다.
• 원추형 기하학과 인터페이스의 매개변수화를 통해 유도된 곡선 인터페이스 패치에서의 추적 부등식을 사용하여 경계항을 제어한다.
• 쌍대성 원리를 적용하여 $L^2$ 오차 추정치를 도출함으로써, $H^1$-노름을 초월한 수렴 분석을 확장한다.

실험 결과

연구 질문

• RQ1 타원형 인터페이스 문제에 대해 메쉬가 맞지 않는 경우에도 최적의 $h$-수렴과 근접한 최적의 $p$-수렴을 유지할 수 있는 $hp$-버전의 유한요소법을 구성할 수 있는가?
• RQ2 대칭 형식과 비대칭 형식의 선택이 $hp$-IPFEM의 안정성과 수렴성에 어떤 영향을 미치는가?
• RQ3 파손된 $H^1$ 노름에서 메쉬가 맞지 않는 $hp$-FEM의 안정성과 최적 수렴을 보장하기 위해 최적의 페널티 파arameter 스케일링은 무엇인가?
• RQ4 메쉬가 맞지 않고 인터페이스가 요소를 가로지르는 상황에서도, 표준 FEM과 유사한 $L^2$ 수렴률을 달성할 수 있는가?
• RQ5 인터페이스의 기하학적 가정(예: $C^2$-연속성)과 메쉬의 형태 정규성은 오차 추정치에 어떤 영향을 미치는가?

주요 결과

• 메쉬의 인터페이스 위치에 관계없이, 파손된 $H^1$ 노름에서 최적의 $h$-수렴 속도를 달성하며, 오차 추정치는 인터페이스의 메쉬에 대한 위치에 독립적이다.
• $p$-수렴 속도는 $p$의 반차수만큼 열등하여, 즉 $O(h^p + h^{p+1/2})$이며, 이는 이와 같은 맥락에서 알려진 바 있는 최고의 결과이다.
• $L^2$ 노름에서의 오차 추정치는 쌍대성 원리를 통해 도출되었으며, 표준 정규성 가정 하에 수렴이 확인되었다.
• 대칭 및 비대칭 형식 모두 안정적이고 수렴하는 해법을 제공하며, 대칭 형식은 더 나은 조건수를 제공한다.
• 분석은 원추 기반 기하학과 매개변수화를 통해 유도된 새로운 곡선 인터페이스 패치에서의 추적 부등식에 의존하며, 소형 요소에 대해서도 강건성을 보장한다.
• 인터페이스가 요소를 임의의 위치에서 가로지르더라도 메쉬 재구성 또는 복잡한 기저 함수 수정 없이도 균일한 안정성과 수렴성을 유지한다.