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

09/28/2022
by   Jinze Li, et al.
0

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.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
02/01/2021

Quinpi: integrating conservation laws with CWENO implicit methods

Many interesting applications of hyperbolic systems of equations are sti...
research
07/27/2023

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

Many interesting physical problems described by systems of hyperbolic co...
research
03/23/2021

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

A single-step high-order implicit time integration scheme for the soluti...
research
04/24/2022

Design of DIRK Schemes with High Weak Stage Order

Runge-Kutta (RK) methods may exhibit order reduction when applied to cer...
research
11/12/2019

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

The threshold dynamics algorithm of Merriman, Bence, and Osher is only f...
research
02/28/2022

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...
research
09/01/2023

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