Arbitrary-order functionally fitted energy-diminishing methods for gradient systems
arXiv:1801.08484
Abstract
It is well known that for gradient systems in Euclidean space or on a Riemannian manifold, the energy decreases monotonically along solutions. In this letter we derive and analyse functionally fitted energy-diminishing methods to preserve this key property of gradient systems. It is proved that the novel methods are unconditionally energy-diminishing and can achieve damping for very stiff gradient systems. We also show that the methods can be of arbitrarily high order and discuss their implementations. A numerical test is reported to illustrate the efficiency of the new methods in comparison with three existing numerical methods in the literature.
References in corpus (3)
- Exponential integrators preserving first integrals or Lyapunov functions for conservative or dissipative systems
- A unified approach to Hamiltonian systems, Poisson systems, gradient systems, and systems with Lyapunov functions and/or first integrals
- Trigonometric collocation methods based on Lagrange basis polynomials for multi-frequency oscillatory second-order differential equations