[논문 리뷰] Adaptive discontinuous Galerkin methods for nonlinear parabolic problems

Stephen Metcalfe|arXiv (Cornell University)|2015. 04. 10.

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

한 줄 요약

이 논문은 비선형 포물형 문제에 적용된 불연속 갈레르킨(dG) 방법을 위한 사후 오차 추정기법과 적응형 알고리즘을 개발한다. 대상 문제로는 대류-확산 방정식, 상미분방정식과 편미분방정식에서의 폭발 현상, 비선형 인터페이스 문제 등이 포함된다. 주요 기여는 잔차 기반 오차 추정기법을 활용한 공간-시간 적응형 알고리즘 설계로, 이는 최적 수렴 속도를 달성하고 폭발 시점과 날카로운 경계층을 정확히 탐지한다.

ABSTRACT

This work is devoted to the study of a posteriori error estimation and adaptivity in parabolic problems with a particular focus on spatial discontinuous Galerkin (dG) discretisations. We begin by deriving an a posteriori error estimator for a linear non-stationary convection-diffusion problem that is discretised with a backward Euler dG method. An adaptive algorithm is then proposed to utilise the error estimator. The effectiveness of both the error estimator and the proposed algorithm is shown through a series of numerical experiments. Moving on to nonlinear problems, we investigate the numerical approximation of blow-up. To begin this study, we first look at the numerical approximation of blow-up in nonlinear ODEs through standard time stepping schemes. We then derive an a posteriori error estimator for an implicit-explicit (IMEX) dG discretisation of a semilinear parabolic PDE with quadratic nonlinearity. An adaptive algorithm is proposed that uses the error estimator to approach the blow-up time. The adaptive algorithm is then applied in a series of test cases to gauge the effectiveness of the error estimator. Finally, we consider the adaptive numerical approximation of a nonlinear interface problem that is used to model the mass transfer of solutes through semi-permiable membranes. An a posteriori error estimator is proposed for the IMEX dG discretisation of the model and its effectiveness tested through a series of numerical experiments.

연구 동기 및 목표

• 선형 비정상 대류-확산 문제에 대해 후행 오차 추정기법을 개발한다. 이는 후행 오차 추정기법을 사용한 역행 오차 dG 이산화를 대상으로 한다.
• 이 추정기법을 활용해 경계층과 내부 경계층을 효율적으로 해결하는 적응형 알고리즘을 설계한다. 이는 대류 지배 문제에서의 해상도를 향상시킨다.
• 이 프레임워크를 비선형 문제로 확장한다. 특히 상미분방정식과 반선형 포물형 편미분방정식에서 유한시간 내 폭발 현상이 나타나는 문제들을 대상으로 한다.
• 반응물이 반투과성 막을 통해 이동하는 과정을 모델링하는 비선형 인터페이스 문제에 대해 사후 오차 추정기법을 구성하고 검증한다.
• 적응형 알고리즘이 날카로운 경계층을 해상하고 폭발 시점을 최적 수렴 속도로 수렴시키는 데 효과적임을 입증한다.

제안 방법

• 선형 비정상 대류-확산 방정식에 대해 후행 오차 dG 이산화를 적용한 후행 오차 추정기법을 유도한다. 이는 dG 공간 이산화와 역행 오차 시간 스텝을 사용한다.
• 오차 추정기법을 활용해 국소 오차 지표에 기반해 공간 및 시간 메쉬를 동적으로 정밀화하는 적응형 알고리즘을 제안한다.
• 이와 동일한 적응형 프레임워크를 이차 비선형성을 가진 반선형 포물형 편미분방정식에 적용한다. 이는 강성 문제 처리를 위해 IMEX dG 이산화를 사용한다.
• 폭발 문제의 맥락에서 IMEX dG 방법에 대한 잔차 기반 사후 오차 추정기법을 개발한다. 오차 한계는 그론발라의 부등식을 통한 추론 방식으로 유도된다.
• 인터페이스 문제에 대해 새로운 적응형 알고리즘을 도입한다. 이는 활성 영역 이외의 영역에서 과다 정밀화되지 않도록 하면서도 인터페이스 경계층을 탐지하고 해상화한다.
• 수치 실험을 통해 공간 및 시간 오차 추정기법의 최적 수렴 속도를 검증하고, 폭발 탐지에서의 강건성을 입증한다.

실험 결과

연구 질문

• RQ1 후행 오차 dG 이산화를 사용한 선형 포물형 문제에 대해 잔차 기반 사후 오차 추정기법을 유도할 수 있으며, 이는 최적 수렴 속도를 달성할 수 있는가?
• RQ2 오차 추정기법을 기반으로 한 메쉬 및 시간 스텝 정밀화가 경계층과 내부 경계층을 효과적으로 해결할 수 있는가?
• RQ3 사후 오차 추정 및 적응형 알고리즘을 비선형 상미분방정식과 반선형 편미분방정식에서 폭발 시점 탐지 및 근사에 확장할 수 있는가?
• RQ4 비선형 인터페이스 문제에서 불연속 계수를 가진 상황에서 오차 추정기법이 날카로운 인터페이스 경계층을 얼마나 잘 포착하는가?
• RQ5 해의 활성도가 낮은 영역에서 자유도를 최소화하면서도 최적 수렴 속도를 달성할 수 있는가?

주요 결과

• 선형 대류-확산 문제에 대한 사후 오차 추정기법은 적응형 메쉬 정밀화 조건 하에서 최적의 공간 및 시간 수렴 속도를 나타낸다.
• 적응형 알고리즘은 경계층과 내부 경계층을 성공적으로 해상하며, 모든 시험 케이스에서 공간 및 시간 효과 지수는 1에 매우 가까운 수준을 유지한다.
• 상미분방정식의 폭발 문제에서 오차 추정기법을 기반으로 한 적응형 알고리즘은 수치 실험에서 λ의 수렴 속도를 확인함으로써 폭발 시점에 최적 순서로 수렴한다.
• 이차 비선형성을 가진 반선형 포물형 편미분방정식에서는 IMEX dG 기반 오차 추정기법이 폭발 시점의 정확한 근사를 가능하게 하며, 적응형 메쉬는 폭발 영역 근처에 집중된다.
• 인터페이스 문제 오차 추정기법은 인터페이스 경계층을 성공적으로 탐지하고 해상화하며, 그림 6.3에서 최적의 공간 및 시간 수렴 속도가 관찰된다. 또한, 인터페이스 반대편 영역에서는 과다 정밀화가 방지된다.
• 오차 추정기법의 상수는 그론발라 부등식을 통해 ε⁻¹에 의존함을 시사하며, 이는 대류 지배 영역에서 이 지수적 의존성을 줄이기 위한 향후 분석 개선이 필요함을 시사한다.