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.
Type de document :
Pré-publication, Document de travail
2016
Liste complète des métadonnées

https://hal.archives-ouvertes.fr/hal-01016684
Contributeur : Max Duarte <>
Soumis le : vendredi 1 avril 2016 - 20:08:40
Dernière modification le : samedi 2 avril 2016 - 01:01:25
Document(s) archivé(s) le : lundi 14 novembre 2016 - 13:15:13

Fichiers

paper_ol.pdf
Fichiers produits par l'(les) auteur(s)

Identifiants

  • 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>

Partager

Métriques

Consultations de
la notice

129

Téléchargements du document

47