High order schemes based on operator splitting and deferred corrections for stiff time dependent PDEs

Abstract : We consider quadrature formulas of high order in time based on Radau-type, L-stable implicit Runge-Kutta schemes to solve time dependent stiff PDEs. Instead of solving a large nonlinear system of equations, we develop a method that performs iterative deferred corrections to compute the solution at the collocation nodes of the quadrature formulas. The numerical stability is guaranteed by a dedicated operator splitting technique that efficiently handles the stiffness of the PDEs and provides initial and intermediate solutions to the iterative scheme. In this way the low order approximations computed by a tailored splitting solver of low algorithmic complexity are iteratively corrected to obtain a high order solution based on a quadrature formula. The mathematical analysis of the numerical errors and local order of the method is carried out in a finite dimensional framework for a general semi-discrete problem, and a time-stepping strategy is conceived to control numerical errors related to the time integration. Numerical evidence confirms the theoretical findings and assesses the performance of the method in the case of a stiff reaction-diffusion equation.
Document type :
Preprints, Working Papers, ...
Complete list of metadatas

https://hal.archives-ouvertes.fr/hal-01016684
Contributor : Max Duarte <>
Submitted on : Friday, April 1, 2016 - 8:08:40 PM
Last modification on : Saturday, April 2, 2016 - 1:01:25 AM
Long-term archiving on : Monday, November 14, 2016 - 1:15:13 PM

Files

paper_ol.pdf
Files produced by the author(s)

Identifiers

  • HAL Id : hal-01016684, version 2
  • ARXIV : 1407.0195

Collections

Citation

Max Duarte, Matthew Emmett. High order schemes based on operator splitting and deferred corrections for stiff time dependent PDEs. 2016. ⟨hal-01016684v2⟩

Share

Metrics

Record views

187

Files downloads

315