An efficient solver for the high-order accurate time-discontinuous Galerkin (TDG) method for second-order hyperbolic systems

An efficient predictor/multi-corrector algorithm for the high-order accurate time-discontinuous Galerkin (TDG) method is presented. The algorithm overcomes the demanding storage and computational effort of the direct solution of the TDG matrix equations. The success of the block matrix iterative solution method is a result of a specific form of time approximation. An analysis shows that the algorithm needs only a few iteration passes with a single matrix factorization of size equal to the number of spatial degrees-of-freedom to reach the 2p + 1 order of accuracy of the parent TDG solution obtained from a direct solve. For linear (p = 1) and quadratic (p = 2) approximations in time, each iteration pass retains the unconditional and asymptotic stability of the TDG parent solution. Numerical results demonstrate the efficiency and accuracy of the time-integration algorithm over second-order accurate single-step/single-solve (SS/SS) methods.

[1]  R. Alexander.  Diagonally implicit runge-kutta methods for stiff odes , 1977 .

[2]  E. Hairer,et al.  Solving ordinary differential equations I (2nd revised. ed.): nonstiff problems , 1993 .

[3]  Alexander Graham,et al.  Kronecker Products and Matrix Calculus: With Applications , 1981 .

[4]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[5]  G. Hulbert.  Time finite element methods for structural dynamics , 1992 .

[6]  O. C. Zienkiewicz,et al.  An alpha modification of Newmark's method , 1980 .

[7]  G. Dahlquist.  A special stability problem for linear multistep methods , 1963 .

[8]  Thomas J. R. Hughes,et al.  Space-time finite element methods for second-order hyperbolic equations , 1990 .

[9]  B. Owren,et al.  Alternative integration methods for problems in structural dynamics , 1995 .

[10]  T.-Y. Wu,et al.  An improved predictor/multi-corrector algorithm for a time-discontinuous Galerkin finite element method in structural dynamics , 2000 .

[11]  Nils-Erik Wiberg,et al.  Adaptive finite element procedures for linear and non‐linear dynamics , 1999 .

[12]  Nathan M. Newmark,et al.  A Method of Computation for Structural Dynamics , 1959 .

[13]  T. C. Fung,et al.  Weighting parameters for unconditionally stable higher-order accurate time step integration algorithms. Part 2—second-order equations , 1999 .

[14]  T. Hughes,et al.  Space-time finite element methods for elastodynamics: formulations and error estimates , 1988 .

[15]  M. Austin.  High-order integration of smooth dynamical systems - Theory and numerical experiments , 1993 .

[16]  Christopher E. Kees,et al.  Higher order time integration methods for two-phase flow , 2002 .

[17]  Jintai Chung,et al.  A Time Integration Algorithm for Structural Dynamics With Improved Numerical Dissipation: The Generalized-α Method , 1993 .

[18]  G. Fairweather,et al.  A Survey of Higher-order Methods for the Numerical Integration of Semidiscrete Parabolic Problems , 1984 .

[19]  Gene H. Golub,et al.  Matrix computations , 1983 .

[20]  M. Dokainish,et al.  A survey of direct time-integration methods in computational structural dynamics—I. Explicit methods , 1989 .

[21]  Sergio Pissanetzky,et al.  Sparse Matrix Technology , 1984 .

[22]  O. Heinisch,et al.  Burington, R.S.: Handbook of Mathematical Tables and Formulas. 4. Aufl., McGraw‐Hill Book Company, New York, London 1965. XI + 423 S., Preis $ 3,25 , 1966 .

[23]  Thomas J. R. Hughes,et al.  Improved numerical dissipation for time integration algorithms in structural dynamics , 1977 .

[24]  Kumar K. Tamma,et al.  Design, analysis, and synthesis of generalized single step single solve and optimal algorithms for structural dynamics , 2004 .

[25]  Antonio Huerta,et al.  High-order accurate time-stepping schemes for convection-diffusion problems , 2000 .

[26]  An Advanced Time-Discontinuous Galerkin Finite Element Method for Structural Dynamics , 2001 .

[27]  Francesco Ubertini,et al.  An efficient integration procedure for linear dynamics based on a time discontinuous Galerkin formulation , 2003 .

[28]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[29]  F. Ubertini,et al.  The Nørsett time integration methodology for finite element transient analysis , 2002 .

[30]  Lonny L. Thompson,et al.  Adaptive space–time finite element methods for the wave equation on unbounded domains , 2005 .

[31]  Nils-Erik Wiberg,et al.  STRUCTURAL DYNAMIC ANALYSIS BY A TIME‐DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD , 1996 .