[논문 리뷰] Robust A Posteriori Error Estimates for Stabilized Finite Element Methods

Lutz Tobiska, R. Verfürth|arXiv (Cornell University)|2014. 02. 24.

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

한 줄 요약

이 논문은 대류-확산 방정식에 적용된 광범위한 안정화된 유한요소법에 대해 강건한 사후 오차 추정을 수립한다. 스트림라인 확산, 국소 투영, 미세구조 스케일, 연속 내부 페널티 방법이 모두 점도, 대류 또는 반응 항에 관계없이 일정한 상수를 갖는 동일한 오차 지표를 도출함을 보여주며, 모든 방법에 대해 균일한 신뢰성과 효율성을 보장한다.

ABSTRACT

There is a wide range of stabilized finite element methods for stationary and non-stationary convection-diffusion equations such as streamline diffusion methods, local projection schemes, subgrid-scale techniques, and continuous interior penalty methods to name only a few. We show that all these schemes give rise to the same robust a posteriori error estimates, i.e. the multiplicative constants in the upper and lower bounds for the error are independent of the size of the convection or reaction relative to the diffusion. Thus, the same error indicator can be used modulo higher order terms caused by data approximation.

연구 동기 및 목표

• 대류-확산 문제에서 다수의 안정화된 유한요소법에 대한 사후 오차 추정을 통합하는 것.
• 다양한 안정화된 방법에서 유도된 오차 지표가 동일한 강건한 오차 경계를 갖는다는 것을 증명하는 것.
• 대류 지배 영역에서도 점도 ε 및 반응 매개변수 β에 관계없이 오차 추정이 균일하게 유지되도록 보장하는 것.
• 대류 또는 반응 항의 상대 크기에 영향을 받지 않는 계산 가능한 상한 및 하한을 잔차에 대해 도출하는 것.
• 잔차 기반 오차 추정기법을 정적 및 비정적 문제 모두에 확장하여 잔차와 일致성 오차 통제를 기반으로 하는 통합 프레임워크를 사용하는 것.

제안 방법

• [39]의 일반적 프레임워크를 사용하여 오차와 잔차 사이의 일반적 강건한 등가성을 수립한다.
• 레마 2.3–2.6를 통해 각 안정화된 방법의 일치 오차에 대한 명시적이고 계산 가능한 상한을 유도한다.
• 대칭 부분에 맞는 에너지 노름과 쌍대 노름을 적용하여 강건성을 확보한다.
• 다양한 방법에 대해 국소 기여도 η, θ, Θ를 포함하는 잔차 기반 오차 추정기법을 적용하며, 연속 내부 페널티 방법의 경우 파라미터 σ_cip를 사용한다.
• 시간 이산화된 오차 추정기법을 비정적 문제에 적용하여 시간 평균 자료와 이산 시간 단계를 활용한다.
• 공간과 시간 모두에서 오차의 전역 상한과 국소 하한을 수립하며, 이 상한과 하한의 상수는 ε 및 β에 의존하지 않는다.

실험 결과

연구 질문

• RQ1 다양한 안정화된 유한요소법에 대해 하나의 사후 오차 추정기법이 대류-확산 문제에 대해 강건하게 적용될 수 있는가?
• RQ2 다른 안정화된 방법에서 유도된 오차 추정이 대류 또는 반응 항의 상대 크기에 관계없이 균일하게 신뢰성 있고 효율적인가?
• RQ3 각 안정화 방법에서 유도되는 일치 오차가 메esh 파라미터에 대해 균일하고 명시적으로 유계일 수 있는가?
• RQ4 동일한 강건성 특성을 갖는 한에서, 동일한 오차 지표가 정적 및 비정적 대류-확산 방정식에 모두 적용 가능한가?
• RQ5 모든 안정화된 방법에 대해 점도 ε 및 반응 매개변수 β에 대해 강건한 잔차 기반 오차 추정기법을 만들 수 있는가?

주요 결과

• 검토된 모든 안정화된 유한요소법—스트림라인 확산, 국소 투영, 미세구조 스케일, 연속 내부 페널티 방법—이 동일한 강건한 사후 오차 추정을 도출한다.
• 오차의 상한 및 하한은 ε 및 β에 대해 균일하며, 이는 경계 상수들이 대류 또는 반응 강도에 의존하지 않음을 의미한다.
• 정적 경우, 잔차 노름은 오차 추정기 항인 η, θ, Θ에 의해 상하한으로 유계지며, 상수 c♭ 및 c♭는 ε 및 β에 의존하지 않는다.
• 비정적 경우, 전역 오차는 시간 단계에 걸친 추정기 합으로 상한으로 유계지며, 상수 c* 및 c*는 T, ε, β에 의존하지 않는다.
• 시간 이산화된 오차 추정기에는 초기 자료, 시간 진화, 자료 근사 기여도가 포함되며, 이 모든 것이 균일하게 유계이다.
• 연속 내부 페널티 방법의 경우 파라미터 σ_cip는 1이며, 다른 방법의 경우 0으로 설정되며, 이는 방법별로 특화된 일치 오차 기여도를 반영한다.