An explicit and practically invariants-preserving method for conservative systems

by   Wenjun Cai, et al.

An explicit numerical strategy that practically preserves invariants is derived for conservative systems by combining an explicit high-order Runge-Kutta (RK) scheme with a simple modification of the standard projection approach, which is named the explicit invariants-preserving (EIP) method. The proposed approach is shown to have the same order as the underlying RK method, while the error of invariants is analyzed in the order of π’ͺ(h^2(p+1)), where h is the time step and p represents the order of the method. When p is appropriately large, the EIP method is practically invariants-conserving because the error of invariants can reach the machine accuracy. The method is illustrated for the cases of single and multiple invariants, with regard to both ODEs and high-dimensional PDEs. Extensive numerical experiments are presented to verify our theoretical results and demonstrate the superior behaviors of the proposed method in a long time numerical simulation. Numerical results suggest that the fourth-order EIP method preserves much better the qualitative properties of the flow than the standard fourth-order RK method and it is more efficient in practice than the fully implicit integrators.


page 13

page 19

page 21

βˆ™ 03/02/2022

High-order linearly implicit schemes conserving quadratic invariants

In this paper, we propose linearly implicit and arbitrary high-order con...
βˆ™ 04/25/2021

A new symmetric linearly implicit exponential integrator preserving polynomial invariants or Lyapunov functions for conservative or dissipative systems

We present a new linearly implicit exponential integrator that preserves...
βˆ™ 11/01/2021

Conservative Integrators for Vortex Blob Methods

Conservative symmetric second-order one-step integrators are derived usi...
βˆ™ 04/12/2023

Fully Conservative Difference Schemes for the Rotation-Two-Component Camassa-Holm System with Smooth/Nonsmooth Initial Data

The rotation-two-component Camassa–Holm system, which possesses strongly...
βˆ™ 07/27/2020

A staggered-projection Godunov-type method for the Baer-Nunziato two-phase model

When describing the deflagration-to-detonation transition in solid granu...
βˆ™ 07/03/2019

Linearly implicit local and global energy-preserving methods for Hamiltonian PDEs

We present linearly implicit methods that preserve discrete approximatio...
βˆ™ 02/10/2023

Multiple-Relaxation Runge Kutta Methods for Conservative Dynamical Systems

We generalize the idea of relaxation time stepping methods in order to p...

Please sign up or login with your details

Forgot password? Click here to reset