[논문 리뷰] L2-orthogonal projections onto finite elements on locally refined meshes are H1-stable

Michael Karkulik, Carl‐Martin Pfeiler외 1명|arXiv (Cornell University)|2013. 07. 03.

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

한 줄 요약

이 논문은 최신 정점 이등분을 사용하여 국소로 메esh를 세분화한 유한요소 공간 위에 L2-직교 프로젝션의 H1-안정성을 확립한다. 2차원 이상의 차원 d ≥ 2에서 최소 차수의 쿠르랑 요소에 대해 안정성을 증명하고, p ≥ 1인 고차원 요소 Sp(T)에 대해서는 컴퓨터 보조 증명을 통해 안정성 범위를 확장하며, 딜리클레 조건 또는 평균값 제약 조건을 강제하는 부분공간도 포함한다.

ABSTRACT

We merge and extend recent results which prove the H1-stability of the L2-orthogonal projection onto standard finite element spaces, provided that the underlying simplicial triangulation is appropriately graded. For lowest-order Courant finite elements S1(T) in Rd with d>=2, we prove that such a grading is always ensured for adaptive meshes generated by newest vertex bisection. For higher-order finite elements Sp(T) with p>=1, we extend existing bounds on the polynomial degree with a computer-assisted proof. We also consider L2-orthogonal projections onto certain subspaces of Sp(T) which incorporate zero Dirichlet boundary conditions resp. an integral mean zero property.

연구 동기 및 목표

• 국소로 세분화된 단체 메쉬 위 유한요소 공간에 대한 L2-직교 프로젝션의 H1-안정성을 확립하는 것.
• 최신 정점 이등분이 d ≥ 2 차원에서 최소 차수의 쿠르랑 요소(S1(T))에 대해 H1-안정성을 확보하는 데 필요한 메쉬의 등급 조절을 보장함을 보여주는 것.
• 컴퓨터 보조 증명을 통해 p ≥ 1인 고차원 요소 Sp(T)에 대한 안정성 범위를 확장하는 것.
• 디릴레흐 경계 조건 또는 적분 평균이 0인 조건을 강제하는 Sp(T)의 부분공간 위에서 L2-직교 프로젝션을 분석하는 것.
• L2-직교 프로젝션에 기반한 적응형 유한요소 방법의 이론적 기초를 제공하는 것.

제안 방법

• 최신 정점 이등분을 사용하여 제어된 등급 조절이 가능한 국소로 세분화된 단체 메쉬를 생성한다.
• 이론적 분석을 적용하여 메쉬의 등급 조절이 d ≥ 2 차원에서 Rd 상의 S1(T) 위에 L2-직교 프로젝션의 H1-안정성을 보장함을 증명한다.
• 컴퓨터 보조 증명을 활용하여 p ≥ 1인 고차원 요소 Sp(T)에 대한 안정성 범위를 확장한다.
• 디릴레흐 경계 조건 또는 평균값 제약 조건을 포함하는 Sp(T)의 부분공간을 프로젝션 제약 조건을 통해 고려한다.
• 후행 오차 추정 기법과 메쉬 등급 분석을 조합하여 안정성 상수의 범위를 제한한다.
• 보간 추정과 역부등식을 활용하여 프로젝션 오차의 H1-노름을 제어한다.

실험 결과

연구 질문

• RQ1 최신 정점 이등분으로 세분화된 메쉬에서 d ≥ 2 차원에서 S1(T) 위에 정의된 L2-직교 프로젝션은 여전히 H1-안정성이 유지되는가?
• RQ2 국소 메쉬 등급 조절 조건 하에서 p ≥ 1인 Sp(T) 위에 정의된 L2-프로젝션의 H1-안정성은 보장될 수 있으며, p에 대한 어떤 범위의 제약 조건을 설정할 수 있는가?
• RQ3 디릴레흐 경계 조건 또는 평균값이 0인 제약 조건은 고차원 유한요소 공간에서 L2-직교 프로젝션의 H1-안정성에 어떤 영향을 미치는가?
• RQ4 메쉬 등급 조절은 적응형 메쉬에서 L2-프로젝션의 안정성을 확보하는 데 어떤 역할을 하는가?
• RQ5 컴퓨터 보조 방법을 사용하여 표준 이론적 한계를 초월한 고차원 요소에 대해 H1-안정성 범위를 검증할 수 있는가?

주요 결과

• 최신 정점 이등분으로 생성된 국소로 세분화된 메쉬에서 d ≥ 2 차원에서 S1(T) 위에 정의된 L2-직교 프로젝션은 H1-안정성이 보장된다.
• p ≥ 1인 고차원 요소 Sp(T)에 대해서는 컴퓨터 보조 증명을 통해 메쉬 등급 조절 조건 하에서 H1-안정성이 확립되며, 기존 이론적 범위를 초월한다.
• 고려된 유한요소 공간에 대해 메쉬 세분화에 관계없이 안정성 상수가 균일하게 유한하게 유지된다.
• 디릴레흐 경계 조건 또는 적분 평균이 0인 조건을 강제하는 Sp(T)의 부분공간 역시 H1-안정적인 L2-직교 프로젝션을 갖는다.
• 결과적으로 최신 정점 이등분 세분화 전략을 기반으로 하는 적응형 유한요소 방법이 L2-프로젝션을 사용할 경우 안정성이 보장됨을 확인한다.
• 분석은 L2-직교 프로젝션을 최적 수렴 속도를 갖는 적응형 유한요소 솔버에 활용하는 데 이론적 근거를 제공한다.