Saved in:
Bibliographic Details
Main Authors: Akramov, Ikrom, Götschel, Sebastian, Minion, Michael, Ruprecht, Daniel, Speck, Robert
Format: Preprint
Published: 2023
Subjects:
Online Access:https://arxiv.org/abs/2310.08352
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866914795048927232
author Akramov, Ikrom
Götschel, Sebastian
Minion, Michael
Ruprecht, Daniel
Speck, Robert
author_facet Akramov, Ikrom
Götschel, Sebastian
Minion, Michael
Ruprecht, Daniel
Speck, Robert
contents Spectral deferred corrections (SDC) are a class of iterative methods for the numerical solution of ordinary differential equations. SDC can be interpreted as a Picard iteration to solve a fully implicit collocation problem, preconditioned with a low-order method. It has been widely studied for first-order problems, using explicit, implicit or implicit-explicit Euler and other low-order methods as preconditioner. For first-order problems, SDC achieves arbitrary order of accuracy and possesses good stability properties. While numerical results for SDC applied to the second-order Lorentz equations exist, no theoretical results are available for SDC applied to second-order problems. We present an analysis of the convergence and stability properties of SDC using velocity-Verlet as the base method for general second-order initial value problems. Our analysis proves that the order of convergence depends on whether the force in the system depends on the velocity. We also demonstrate that the SDC iteration is stable under certain conditions. Finally, we show that SDC can be computationally more efficient than a simple Picard iteration or a fourth-order Runge-Kutta-Nyström method.
format Preprint
id arxiv_https___arxiv_org_abs_2310_08352
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle Spectral deferred correction methods for second-order problems
Akramov, Ikrom
Götschel, Sebastian
Minion, Michael
Ruprecht, Daniel
Speck, Robert
Numerical Analysis
Spectral deferred corrections (SDC) are a class of iterative methods for the numerical solution of ordinary differential equations. SDC can be interpreted as a Picard iteration to solve a fully implicit collocation problem, preconditioned with a low-order method. It has been widely studied for first-order problems, using explicit, implicit or implicit-explicit Euler and other low-order methods as preconditioner. For first-order problems, SDC achieves arbitrary order of accuracy and possesses good stability properties. While numerical results for SDC applied to the second-order Lorentz equations exist, no theoretical results are available for SDC applied to second-order problems. We present an analysis of the convergence and stability properties of SDC using velocity-Verlet as the base method for general second-order initial value problems. Our analysis proves that the order of convergence depends on whether the force in the system depends on the velocity. We also demonstrate that the SDC iteration is stable under certain conditions. Finally, we show that SDC can be computationally more efficient than a simple Picard iteration or a fourth-order Runge-Kutta-Nyström method.
title Spectral deferred correction methods for second-order problems
topic Numerical Analysis
url https://arxiv.org/abs/2310.08352