TY - JOUR
T1 - Temporal stability of collocation, Petrov–Galerkin, and other non-symmetric methods in elastodynamics and an energy conserving time integration
AU - Wang, Jiarui
AU - Hillman, Michael C.
N1 - Funding Information:
This work was supported by the National Science Foundation, USA award number 1826221 and the U.S. Army Engineering Research and Development Center through Ordnance Technology Initiative Agreement, USA DOTC-17-01-INIT0880 . Proofing of this manuscript by Jennifer Dougal is also acknowledged and appreciated.
Publisher Copyright:
© 2022 Elsevier B.V.
PY - 2022/4/1
Y1 - 2022/4/1
N2 - Non-symmetric matrices may arise in the discretization of self-adjoint problems when a Petrov–Galerkin, collocation, or finite-volume method is employed. While these methods have been widely applied, in this paper it is shown that the use of these non-symmetric matrices is incompatable with the conservation of energy in elastodynamics. First, the consistency between the continuous forms of the momentum equation and the energy equation is examined. It is shown that the conservation of linear momentum is equivalent to conservation of energy provided the solution is sufficiently smooth. The semi-discrete counterparts are then analyzed, where it is demonstrated that they are also equivalent, but only conditionally: the mass and stiffness matrices must be symmetric. As a result, employing a non-symmetric method in elastodynamics may artificially generate or dissipate energy. The fully discrete forms with Newmark time integration are then examined where it is shown that unconditionally unstable algorithms may arise. An energy-conserving time integration algorithm is then proposed which provides stability in the solutions of non-symmetric systems. The collocation and finite-volume methods are employed in numerical examples to demonstrate stability issues and the effectiveness of the proposed time integration methodology.
AB - Non-symmetric matrices may arise in the discretization of self-adjoint problems when a Petrov–Galerkin, collocation, or finite-volume method is employed. While these methods have been widely applied, in this paper it is shown that the use of these non-symmetric matrices is incompatable with the conservation of energy in elastodynamics. First, the consistency between the continuous forms of the momentum equation and the energy equation is examined. It is shown that the conservation of linear momentum is equivalent to conservation of energy provided the solution is sufficiently smooth. The semi-discrete counterparts are then analyzed, where it is demonstrated that they are also equivalent, but only conditionally: the mass and stiffness matrices must be symmetric. As a result, employing a non-symmetric method in elastodynamics may artificially generate or dissipate energy. The fully discrete forms with Newmark time integration are then examined where it is shown that unconditionally unstable algorithms may arise. An energy-conserving time integration algorithm is then proposed which provides stability in the solutions of non-symmetric systems. The collocation and finite-volume methods are employed in numerical examples to demonstrate stability issues and the effectiveness of the proposed time integration methodology.
UR - http://www.scopus.com/inward/record.url?scp=85126308740&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=85126308740&partnerID=8YFLogxK
U2 - 10.1016/j.cma.2022.114738
DO - 10.1016/j.cma.2022.114738
M3 - Article
AN - SCOPUS:85126308740
SN - 0045-7825
VL - 393
JO - Computer Methods in Applied Mechanics and Engineering
JF - Computer Methods in Applied Mechanics and Engineering
M1 - 114738
ER -