Related papers: A Two Step, Fourth Order, Nearly-Linear Method wit…
In this paper, we consider the numerical approximations for the fourth order Cahn-Hilliard equation with concentration dependent mobility, and the logarithmic Flory-Huggins potential. One challenge in solving such a diffusive system…
In this work, we consider the error estimates and the long-time conservation or near-conservation of geometric structures, including energy, mass shell and phase-space volume, for four two-step symmetric methods applied to relativistic…
In this paper, we systematically construct two classes of structure-preserving schemes with arbitrary order of accuracy for canonical Hamiltonian systems. The one class is the symplectic scheme, which contains two new families of…
Recently, a whole class of evergy-preserving integrators has been derived for the numerical solution of Hamiltonian problems. In the mainstream of this research, we have defined a new family of symplectic integrators depending on a real…
For Hamiltonian systems with non-canonical structure matrices, a new family of fourth-order energy-preserving integrators is presented. The integrators take a form of a combination of Runge--Kutta methods and continuous-stage Runge--Kutta…
In this paper, a novel sixth order energy-conserved method is proposed for solving the three-dimensional time-domain Maxwell's equations. The new scheme preserves five discrete energy conservation laws, three momentum conservation laws,…
We design a novel, exactly energy-conserving implicit non-symplectic integration method for an eight-dimensional Hamiltonian system with four degrees of freedom. In our algorithm, each partial derivative of the Hamiltonian with respect to…
A rotation-two-component Camassa-Holm (R2CH) system was proposed recently to describe the motion of shallow water waves under the influence of gravity. This is a highly nonlinear and strongly coupled system of partial differential…
In this paper, we consider numerical approximations for the viscous Cahn-Hilliard equation with hyperbolic relaxation. This type of equations processes energy-dissipative structure. The main challenge in solving such a diffusive system…
In this paper, we develop a novel class of arbitrarily high-order energy-preserving schemes for the Camassa-Holm equation. With the aid of the invariant energy quadratization approach, the Camassa-Holm equation is first reformulated into an…
In this paper, we present a novel class of high-order energy-preserving schemes for solving the Zakharov-Rubenchik equations. The main idea of the scheme is first to introduce an quadratic auxiliary variable to transform the Hamiltonian…
We propose an energy-optimized invariant energy quadratization method to solve the gradient flow models in this paper, which requires only one linear energy-optimized step to correct the auxiliary variables on each time step. In addition to…
This work proposes a suite of numerical techniques to facilitate the design of structure-preserving integrators for nonlinear dynamics. The celebrated LaBudde-Greenspan integrator and various energy-momentum schemes adopt a difference…
The two-dimensional n-body problem of classical mechanics is a non-integrable Hamiltonian system for n > 2. Traditional numerical integration algorithms, which are polynomials in the time step, typically lead to systematic drifts in the…
Kahan's method and a two-step generalization of the discrete gradient method are both linearly implicit methods that can preserve a modified energy for Hamiltonian systems with a cubic Hamiltonian. These methods are here investigated and…
In the last few decades, numerical simulation for nonlinear oscillators has received a great deal of attention, and many researchers have been concerned with the design and analysis of numerical methods for solving oscillatory problems. In…
Existing model reduction techniques for high-dimensional models of conservative partial differential equations (PDEs) encounter computational bottlenecks when dealing with systems featuring non-polynomial nonlinearities. This work presents…
In this paper, we introduce a Lagrange multiplier approach to construct linearly implicit energy-preserving schemes of arbitrary order for general Hamiltonian PDEs. Unlike the widely used auxiliary variable methods, this novel approach does…
Structure-preserving algorithms and algorithms with uniform error bound have constituted two interesting classes of numerical methods. In this paper, we blend these two kinds of methods for solving nonlinear Hamiltonian systems with highly…
One main issue, when numerically integrating autonomous Hamiltonian systems, is the long-term conservation of some of its invariants, among which the Hamiltonian function itself. For example, it is well known that classical symplectic…