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

03/23/2021
by   Chongmin Song, et al.
0

A single-step high-order implicit time integration scheme for the solution of transient and wave propagation problems is presented. It is constructed from the Padé expansions of the matrix exponential solution of a system of first-order ordinary differential equations formulated in the state-space. A computationally efficient scheme is developed exploiting the techniques of polynomial factorization and partial fractions of rational functions, and by decoupling the solution for the displacement and velocity vectors. An important feature of the novel algorithm is that no direct inversion of the mass matrix is required. From the diagonal Padé expansion of order M a time-stepping scheme of order 2M is developed. Here, each elevation of the accuracy by two orders results in an additional system of real or complex sparse equations to be solved. These systems are comparable in complexity to the standard Newmark method, i.e., the effective system matrix is a linear combination of the static stiffness, damping, and mass matrices. It is shown that the second-order scheme is equivalent to Newmark's constant average acceleration method, often also referred to as trapezoidal rule. The proposed time integrator has been implemented in MATLAB using the built-in direct linear equation solvers. In this article, numerical examples featuring nearly one million degrees of freedom are presented. High-accuracy and efficiency in comparison with common second-order time integration schemes are observed. The MATLAB-implementation is available from the authors upon request or from the GitHub repository (to be added).

READ FULL TEXT

page 29

page 31

page 34

page 36

research
06/08/2022

High-order implicit time integration scheme with controllable numerical dissipation based on mixed-order Padé expansions

A single-step high-order implicit time integration scheme with controlla...
research
09/28/2022

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

This paper develops an implicit family of sub-step integration algorithm...
research
11/12/2018

An Application of Rubi: Series Expansion of the Quark Mass Renormalization Group Equation

We highlight how Rule-based Integration (Rubi) is an enhanced method of ...
research
09/20/2023

Matrix-based implementation and GPU acceleration of linearized ordinary state-based peridynamic models in MATLAB

Ordinary state-based peridynamic (OSB-PD) models have an unparalleled ca...
research
06/11/2021

High order integrators obtained by linear combinations of symmetric-conjugate compositions

A new family of methods involving complex coefficients for the numerical...
research
12/23/2020

Optimal and Low-Memory Near-Optimal Preconditioning of Fully Implicit Runge-Kutta Schemes for Parabolic PDEs

Runge-Kutta (RK) schemes, especially Gauss-Legendre and some other fully...

Please sign up or login with your details

Forgot password? Click here to reset