RATTLE Method for Dissipative Constrained Hamiltonian Systems

摘要

Firstly, the Hamiltonian formulation for the constrained dissipative systems was deduced, it is index-3 DAEs. The index-2 DAEs was obtained based on the GGL stabilized method. Then RATTLE method that proposed for motion equations of conservative constrained Hamiltonian systems was extended to solve dissipative constrained Hamiltonian systems. For the dissipative systems, the symplectic structure of RATTLE method is no longer preserved, but this method can capture the decay of the energy accurately because of no numerical dissipation. Numerical experiment results illustrate the effectiveness of the method.