Gianluca Frasca-Caccia

Luigi Brugnano, Gianluca Frasca-Caccia, Felice Iavernaro

In this paper, we report about recent findings in the numerical solution of Hamiltonian Partial Differential Equations (PDEs), by using energy-conserving line integral methods in the Hamiltonian Boundary Value Methods (HBVMs) class. In particular, we consider the semilinear wave equation, the nonlinear Schrödinger equation, and the Korteweg-de Vries equation, to illustrate the main features of this novel approach.

Gianluca Frasca-Caccia, Peter E. Hydon

Conservation laws are among the most fundamental geometric properties of a partial differential equation (PDE), but few known finite difference methods preserve more than one conservation law. All conservation laws belong to the kernel of the Euler operator, an observation that was first used recently to construct approximations symbolically that preserve two conservation laws of a given PDE. However, the complexity of the symbolic computations has limited the effectiveness of this approach. The current paper introduces some key simplifications that make the symbolic-numeric approach feasible. To illustrate the simplified approach, we derive bespoke finite difference schemes that preserve two discrete conservation laws for the Korteweg-de Vries (KdV) equation and for a nonlinear heat equation. Numerical tests show that these schemes are robust and highly accurate compared to others in the literature.

Gianluca Frasca-Caccia

By exploiting the fact that conservation laws form the kernel of a discrete Euler operator, we use a recently introduced symbolic-numeric approach to construct a new class of finite difference methods for the modified Korteweg-de Vries (mKdV) equation, that preserve the local conservation laws of mass and energy.