An explicit and practically invariants-preserving method for conservative systems

09/15/2020
βˆ™
by   Wenjun Cai, et al.
βˆ™
0
βˆ™

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.

READ FULL TEXT

page 13

page 19

page 21

research
βˆ™ 03/02/2022

High-order linearly implicit schemes conserving quadratic invariants

In this paper, we propose linearly implicit and arbitrary high-order con...
research
βˆ™ 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...
research
βˆ™ 11/01/2021

Conservative Integrators for Vortex Blob Methods

Conservative symmetric second-order one-step integrators are derived usi...
research
βˆ™ 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...
research
βˆ™ 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...
research
βˆ™ 07/03/2019

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

We present linearly implicit methods that preserve discrete approximatio...
research
βˆ™ 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