[논문 리뷰] On Local Super-Penalization of Interior Penalty Discontinuous Galerkin Methods

Andrea Cangiani, John Chapman외 2명|arXiv (Cornell University)|2012. 05. 25.

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

한 줄 요약

이 논문은 내부 페널티 계수의 무한으로 수렴함에 따라 불연속 갈레르킨(dG) 방법이 연속 갈레르킨(cG) 방법으로 수렴함을 이론적으로 및 수치적으로 입증한다. 비선형 및 시간에 의존하는 설정에서도 성립한다. 또한 局부 초페널티화를 도입하여, 안정성에 손상을 주지 않고 자유도를 줄일 수 있는 하이브리드 연속-불연속 갈레르킨 방법을 개발하였으며, 이는 적응 메쉬 정밀화와 부동산 확산 및 혼합 이동 문제에 대한 수치 실험을 통해 검증되었다.

ABSTRACT

We prove in an abstract setting that standard (continuous) Galerkin finite element approximations are the limit of interior penalty discontinuous Galerkin approximations as the penalty parameter tends to infinity. We apply this result to equations of non-negative characteristic form and the non-linear, time dependent system of incompressible miscible displacement. Moreover, we investigate varying the penalty parameter on only a subset of a triangulation and the effects of local super-penalization on the stability of the method, resulting in a partly continuous, partly discontinuous method in the limit. An iterative automatic procedure is also proposed for the determination of the continuous region of the domain without loss of stability of the method.

연구 동기 및 목표

• 내부 페널티 dG 방법이 페널티 계수의 무한으로 수렴함에 따라 cG 방법으로 수렴하는 일반적인 이론적 프레임워크를 수립하는 것.
• 특정 요소나 면에만 고 페널티를 적용하는 국소적 초페널티화의 영향을 조사하여, 하이브리드 연속-불연속 갈레르킨(cdG) 방법을 만드는 것.
• 해당 영역 내 안정적인 연속 영역을 자동으로 식별하는 반복적이고 자동화된 절차를 개발하고 테스트하는 것.
• 특히 경계층이나 급격한 프론트가 존재하는 부동산 지배 및 비선형 시간에 의존하는 문제에서 cdG 방법의 안정성과 정확도를 평가하는 것.
• dG 방법의 안정성과 정확도를 유지하면서 국소적으로 연속적인 요소를 사용하여 자유도를 줄이는 것.

제안 방법

• 무한 페널티에서 dG에서 cG로의 수렴을 증명하기 위해, inf-sup 조건과 에너지 추정에 기반한 추상적 프레임워크를 사용한다.
• 비음성 특성 형태의 방정식과 비선형 불압축성 혼합 이동 시스템에 이 프레임워크를 적용하며, 혼합 Raviart-Thomas-dG 포지션을 사용한다.
• 메쉬 면의 일부에만 페널티 계수를 변화시켜 국소적 초페널티화를 도입함으로써 하이브리드 cdG 방법을 생성한다.
• 요소 간 농도 필드의 점프 크기를 기반으로 연속 근사에 적합한 요소를 자동으로 식별하고 표시하기 위한 반복 알고리즘을 제안한다.
• 요소 간 해의 점프에 대한 $L^2$-노름에 대한 허용 오차 기준을 설정하여 연속 영역에 포함할 요소를 표시한다.
• 시간에 따라 $L^2$-노름 차이와 자유도 수를 비교하여 cdG 근사와 완전한 dG 및 cG 해를 평가한다.

실험 결과

연구 질문

• RQ1 내부 페널티 dG 방법이 페널티 계수의 무한으로 수렴함에 따라 비선형 및 시간에 의존하는 문제에서도 표준 cG 방법으로 수렴하는가?
• RQ2 특정 요소 면에만 고 페널티를 적용하는 국소적 초페널티화가 안정적이고 정확한 하이브리드 연속-불연속 갈레르킨 방법을 생성할 수 있는가?
• RQ3 국소적으로 연속적인 요소를 사용함으로써 자유도를 얼마나 줄일 수 있으며, 수치적 진동이나 정확도 저하 없이 가능한가?
• RQ4 자동이고 반복적인 절차로 안정적인 연속 영역을 식별할 수 있으며, 해의 안정성에 영향을 주지 않는가?
• RQ5 cdG와 dG 근사 간 오차는 시간에 따라 어떻게 변화하며, 자유도 감소와 정확도 사이의 상충 관계는 어떠한가?

주요 결과

• dG 방법이 페널티 계수의 무한으로 수렴함에 따라 cG 방법으로 수렴함을 확인하였으며, 이는 비선형 및 시간에 의존하는 문제에 적용 가능한 일반적인 추상적 프레임워크 내에서 이론적 수렴을 입증한다.
• 국소적 초페널티화를 통해 자유도를 줄이면서도 dG 수준의 정확도를 유지하는 안정적인 하이브리드 cdG 방법을 구축할 수 있다.
• 불압축성 혼합 이동 문제에서 자동 표시 절차의 허용 오차를 $10^{-5}$로 설정함으로써 자유도를 dG 방법의 74.6%로 줄일 수 있었으며, dG 해와의 $L^2(L^2)$ 노름 오차는 단지 $1.0455 \times 10^{-4}$에 불과했다.
• $L^2(L^2)$ 노름에서 cdG와 dG 근사 간의 차이가 $\mathtt{tol} = 10^{-3}$일 때 $1.2073 \times 10^{-2}$에서 $\mathtt{tol} = 10^{-5}$일 때 $1.0455 \times 10^{-4}$로 감소하여, 더 엄격한 표시 조건에서 정확도 향상을 보였다.
• 수치 결과에 따르면 $\mathtt{tol} = 10^{-4}$ 및 $10^{-3}$일 때 dG와 cdG 근사 간에 눈에 띄는 차이가 없었으며, 완전한 cG 해는 경계층에서 심한 진동을 보였다.
• 연속 영역을 식별하기 위한 반복 절차는 안정성과 정확도를 유지하였으며, 연속 메쉬를 다수 번 재초기화한 후에도 오차가 낮게 유지되었다.