A Mathematical Framework Towards a Unified Set of Discontinuous State-Phase Hierarchical Time Operators for Computational Dynamics
 Source CMES: Computer Modeling in Engineering & Sciences, Vol. 4, No. 1, pp. 103-118, 2003 Download Full length paper in PDF format. Size = 263,228 bytes Keywords Time discretized discontinuous operators; Generalized weighted residual approach; Hierarchical Time Formulations; Computational Dynamics Abstract Of general interest here is the time dimension aspect wherein discretized operators in time may be continuous or discontinuous; and of particular interest and focus here is the design of time discretized operators in the context of discontinuous state-phase for computational dynamics applications. Based on a generalized bi-discontinuous time weighted residual formulation, the design leading to a new unified set of hierarchical energy conserving and energy dissipating time discretized operators are developed for the first time that are fundamentally useful for time adaptive computations for dynamic problems. Unlike time discontinuous Galerkin approaches, the design is based upon a time discontinuous Petrov-Galerkin-like approach employing an asymptotic series type approximations for the state variables involving derivatives at the beginning of the time step. As a consequence, this enables to design methods that have spectral properties corresponding to the diagonal, first sub-diagonal and second sub-diagonal Pad\'{e} entries. Thus, \emph {A}-stable schemes of order$2q$, \emph {L}-stable schemes of order$2q-1$ and$2q-2$ are obtained. The spectral equivalent algorithms to the diagonal Pad\'{e} entry are energy conserving algorithms. The spectral equivalent algorithms to the first and second sub-diagonal Pad\'{e} entries are energy dissipating algorithms with the property of asymptotic annihilation of the high frequency response. Additionally, these time operators naturally inherit a hierarchical structure that are extremely useful for time adaptive computations. Moreover, since Pad\'{e} entries have the lowest relative error the developed schemes are optimal in terms of order of accuracy in time, dissipation, dispersion and zero-order displacement and velocity overshoot characteristics.