[논문 리뷰] Numerical study of goal-oriented error control for stabilized finite element methods

Marius Paul Bruchhäuser, Kristina Schwegler외 1명|arXiv (Cornell University)|2018. 03. 27.

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

한 줄 요약

이 논문은 이중 가중 잔여법(DWR)을 안정화된 유한요소 형식과 결합하여 대류지배 문제에 대한 목표 지향적 적응형 유한요소법을 제안한다. 이는 사용자가 정의한 관심량에 대해 강건하고 계산 가능한 오차 통제를 달성하기 위해 이중 문제의 고차원 이산화를 사용한다. 이로 인해 2차원 및 3차원 기준 테스트에서 허위 진동이 크게 감소하고 효과성 지수(effectivity index)가 1에 매우 가까워진다.

ABSTRACT

The efficient and reliable approximation of convection-dominated problems continues to remain a challenging task. To overcome the difficulties associated with the discretization of convection-dominated equations, stabilization techniques and a posteriori error control mechanisms with mesh adaptivity were developed and studied in the past. Nevertheless, the derivation of robust a posteriori error estimates for standard quantities and in computable norms is still an unresolved problem and under investigation. Here we combine the Dual Weighted Residual (DWR) method for goal-oriented error control with stabilized finite element methods. By a duality argument an error representation is derived on that an adaptive strategy is built. The key ingredient of this work is the application of a higher order discretization of the dual problem in order to make a robust error control for user-chosen quantities of interest feasible. By numerical experiments in 2D and 3D we illustrate that this interpretation of the DWR methodology is capable to resolve layers and sharp fronts with high accuracy and to further reduce spurious oscillations.

연구 동기 및 목표

• 기존의 후행 추정기법이 종종 실패하는 대류지배 문제에서 신뢰성 있고 효율적인 오차 통제 문제를 해결한다.
• 기존의 안정화된 유한요소 방법에 대한 후행 오차 추정에서의 강건하지 않거나 계산 불가능한 오차 노름의 한계를 극복한다.
• 전역 노름이 아닌 사용자가 정의한 관심량에 초점을 맞춘 목표 지향적 적응 전략을 개발한다.
• 이중 문제의 근사 향상을 통해 층과 급격한 전 fronts 를 정확하게 포착함으로써 강건성과 정확성을 확보한다.
• 복잡한 속도장이 존재하는 2차원 및 3차원 대류지배 이동 문제에서 이 방법의 효과성을 입증한다.

제안 방법

• 이중 가중 잔여법(DWR)을 적용하여 안정화된 유한요소 해에 대한 목표 지향적 후행 오차 표현식을 유도한다.
• 오차 추정기의 정확도를 향상시키기 위해 이중 문제에 고차원 유한요소 공간을 사용하여 저차원 이중 근사의 함정을 피한다.
• 주 문제에 잔여기반 안정화(SUPG)를 구현하여 대류지배 영역에서 허위 진동을 억제한다.
• 이중 해의 값으로 국소 잔여치를 가중하는 방식으로 계산 가능한 오차 지표를 구성한다.
• 오차 지표에 기반해 메쉬를 적응적으로 세분화하고 굴절시켜 관심량에 가장 큰 영향을 미치는 곳에 자유도를 집중시킨다.
• 주 문제보다 더 높은 차수의 요소를 사용해 이중 문제를 해결함으로써, 층이 풍부한 영역에서도 신뢰성 있고 강건한 오차 추정을 확보한다.

실험 결과

연구 질문

• RQ1 고차원 이중 이산화가 대류지배 문제에 대한 안정화된 유한요소 방법에서 목표 지향적 오차 추정의 강건성과 신뢰성 향상에 기여하는가?
• RQ2 DWR 기반 적응 전략은 허위 진동을 최소화하면서 날카로운 층과 전 fronts 를 얼마나 효과적으로 해결하는가?
• RQ3 제안된 방법은 2차원 및 3차원 설정에서 관심량에 대해 효과성 지수가 1에 매우 가까운지를 확인할 수 있는가?
• RQ4 SUPG 안정화와 목표 지향적 적응성의 조합이 전역 오차 통제보다 더 효율적인 메쉬를 생성하는가?
• RQ5 변동하는 대류장과 복잡한 기하구조를 가진 문제에서 이 방법은 3차원에서 어떻게 성능을 발휘하는가?

주요 결과

• 모든 테스트된 목표 양에 대해 효과성 지수가 1에 매우 가까운 것을 확인하여 매우 강건한 오차 통제가 이루어짐을 입증하였다.
• 특히 경계층 및 내부층에서 수치 해의 허위 진동이 크게 감소하였으며, 이는 개선된 이중 문제 근사 덕분이었다.
• 2차원 및 3차원 테스트 케이스에서 적응형 메쉬 세분화가 흐름 방향에 따라 정렬된 층과 같은 관심량에 중요한 영역에 요소를 집중시키는 데 성공하였다.
• 점값 및 평균값 목표 함수의 경우, 조밀한 초기 메쉬 조건에서도 목표 영역에서 오차가 최소화된 매우 정확한 근사를 제공하였다.
• 변동하는 대류장이 존재하는 3차원 문제에서도 이 방법은 강건성을 유지하여 복잡한 실제 이동 문제에의 적용 가능성을 입증하였다.
• 고차원 이중 이산화를 사용함으로써 저차원 이중 근사보다 더 정확한 오차 추정이 가능해졌으며, 오차 표현식 내의 원천 오차가 감소하였다.