ENERGY AND QUADRATIC INVARIANTS PRESERVING METHODS FOR HAMILTONIAN SYSTEMS WITH HOLONOMIC CONSTRAINTS

We introduce a new class of parametrized structure–preserving partitioned RungeKutta(α-PRK)methods for Hamiltonian systems with holonomic constraints.The methods are symplectic for any fixed scalar parameterα,and are reduced to the usual symplectic PRK methods like Shake-Rattle method or PRK schemes based on Lobatto IIIA-IIIB pairs whenα=0.We provide a new variational formulation for symplectic PRK schemes and use it to prove that theα-PRK methods can preserve the quadratic invariants for Hamiltonian systems subject to holonomic constraints.Meanwhile,for any given consistent initial values(p0,q0)and small step size h>0,it is proved that there existsα∗=α(h,p0,q0)such that the Hamiltonian energy can also be exactly preserved at each step.Based on this,we propose some energy and quadratic invariants preservingα-PRK methods.Theseα-PRK methods are shown to have the same convergence rate as the usual PRK methods and perform very well in various numerical experiments.

QUADRATIC INVARIANTS AND SYMPLECTIC STRUCTURE OF GENERAL LINEAR METHODS

Focuses on a study which presented some invariants and conservation laws of general linear methods applied to differential equation systems. Information on the quadratic invariants; Conservation of symplectic structure; Details on the multiple Runge-Kutta methods; Equations of the one-leg methods.

Semi-Implicit Spectral Deferred Correction Method Based on the Invariant Energy Quadratization Approach for Phase Field Problems

This paper presents a high order time discretization method by combining the semi-implicit spectral deferred correction method with energy stable linear schemes to simulate a series of phase field problems.We start with the linear scheme,which is based on the invariant energy quadratization approach and is proved to be linear unconditionally energy stable.The scheme also takes advantage of avoiding nonlinear iteration and the restriction of time step to guarantee the nonlinear system uniquely solvable.Moreover,the scheme leads to linear algebraic system to solve at each iteration,and we employ the multigrid solver to solve it efficiently.Numerical re-sults are given to illustrate that the combination of local discontinuous Galerkin(LDG)spatial discretization and the high order temporal scheme is a practical,accurate and efficient simulation tool when solving phase field problems.Namely,we can obtain high order accuracy in both time and space by solving some simple linear algebraic equations.