Abstract:A single-step high-order implicit time integration scheme for the solution of transient as well as wave propagation problems is presented. It is constructed from the Padé expansion of the matrix exponential solution of a system of first-order ordinary differential equations formulated in the state–space. By exploiting the techniques of polynomial factorization and partial fractions of rational functions, and by decoupling the solution for the displacement and velocity vectors a computationally efficient scheme is developed. An important feature of the novel algorithm is that no direct inversion of the mass matrix is required. Based on the diagonal Padé expansion of order M a time-stepping scheme of order 2 M is constructed. When M is odd, one system of equation with a real matrix and ( M − 1 ) / 2 systems of equations with complex matrices are solved recursively. When M is even, there are M / 2 systems of equations with complex matrices. These systems are sparse and 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. The derived second-order scheme is analytically equivalent to Newmark constant average acceleration method. The proposed time integrator has been implemented in MATLAB and FORTRAN using direct linear equation solvers. In this article, numerical examples featuring nearly one million degrees of freedom are presented to highlight the exceptional accuracy and efficiency of the novel time-stepping scheme in comparison with established second-order methods. The MATLAB- and FORTRAN-implementations are available from the authors upon request.

Adaptive Third Order Adams-Bashforth Time Stepping for Extended Boussinesq Equations

An Arbitrary Lagrangian–Eulerian Method Based on the Adaptive Riemann Solvers for General Equations of State

An Adaptive Time-Stepping Procedure Based On The Scaled Boundary Finite Element Method For Elastodynamics

A 2DH fully dispersive and weakly nonlinear Boussinesq-type model based on a finite-volume and finite-difference TVD-type scheme

High-order Accurate Implicit-Explicit Time-Stepping Schemes for Wave Equations on Overset Grids

Local Time-Stepping for the Shallow Water Equations using CFL Optimized Forward-Backward Runge-Kutta Schemes

Implicit Adaptive Mesh Refinement for Dispersive Tsunami Propagation

A Numerical Investigation of Matrix-Free Implicit Time-Stepping Methods for Large CFD Simulations

p-adaptive discontinuous Galerkin method for the shallow water equations on heterogeneous computing architectures

High Order Space-Time Adaptive WENO Finite Volume Schemes for Non-Conservative Hyperbolic Systems

Computationally Efficient Adaptive Time Step Method for the Cahn-Hilliard Equation

Space-time adaptive ADER discontinuous Galerkin finite element schemes with a posteriori sub-cell finite volume limiting

Second-order linear adaptive time-stepping schemes for the fractional Allen–Cahn equation

On an Accurate A-posteriori Error Estimator and Adaptive Time Stepping for the Implicit and Explicit Composite Time Integration Algorithms

Efficient computational hybrid method for the solution of 2D multi-term fractional order advection-diffusion equation

ADAPTIVE TIME SPLITTING METHOD FOR MULTI-SCALE EVOLUTIONARY PARTIAL DIFFERENTIAL EQUATIONS

Adaptive Partitioned Methods for the Time-Accurate Approximation of the Evolutionary Stokes–Darcy System

A High-Order Hybrid-Spectral Incompressible Navier-Stokes Model For Nonlinear Water Waves

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

Two-level dynamic load-balanced p-adaptive discontinuous Galerkin methods for compressible CFD simulations