High-order accurate multi-sub-step implicit integration algorithms with dissipation control for second-order hyperbolic problems

by   Jinze Li, et al.

This paper develops an implicit family of sub-step integration algorithms, which firstly requires identical effective stiffness matrices and third-order consistency within each sub-step. Consequently, the trapezoidal rule has to be employed in the first sub-step and optimal spectral properties are naturally embedded into the proposed algorithms. The analysis reveals at the first time that the constructed s-sub-step implicit schemes with s≤6 can reach sth-order accuracy when achieving dissipation control and unconditional stability simultaneously. Hence, only four cost-optimal high-order implicit algorithms corresponding to three, four, five, and six sub-steps are developed. Unlike some published high-order algorithms, four novel methods do not suffer from order reduction for solving forced vibrations. Moreover, the novel methods overcome the defect that the authors' previous high-order algorithms require an additional solution to obtain accurate accelerations. Linear and nonlinear examples are solved to confirm the numerical performance and superiority of four novel high-order algorithms.


page 1

page 2

page 3

page 4


Quinpi: integrating conservation laws with CWENO implicit methods

Many interesting applications of hyperbolic systems of equations are sti...

Quinpi: Integrating stiff hyperbolic systems with implicit high order finite volume schemes

Many interesting physical problems described by systems of hyperbolic co...

High-order implicit time integration scheme based on Padé expansions

A single-step high-order implicit time integration scheme for the soluti...

Design of DIRK Schemes with High Weak Stage Order

Runge-Kutta (RK) methods may exhibit order reduction when applied to cer...

Second Order Threshold Dynamics Schemes for Two Phase Motion by Mean Curvature

The threshold dynamics algorithm of Merriman, Bence, and Osher is only f...

Arbitrarily high order implicit ODE integration by correcting a neural network approximation with Newton's method

As a method of universal approximation deep neural networks (DNNs) are c...

A novel Newton-Raphson style root finding algorithm

Many problems in applied mathematics require root finding algorithms. Un...

Please sign up or login with your details

Forgot password? Click here to reset