Local structure-preserving algorithms for general multi-symplectic Hamiltonian PDEs
Introduction
A wide range of -dimensional PDEs can be written into a multi-symplectic Hamiltonian form [1] where and , , are constant skew-symmetric matrices, , , is a state variable, and is a scalar-valued smooth function, such as the Sine-Gordon equation [2], the Korteweg–de Vries (KdV) equation [3], the Camassa–Holm equation [4], the regularized long-wave (RLW) equation [5], the nonlinear Schrödinger (NLS) equation [[6], [7]], Maxwell’s equations [[8], [9]] and so on. It is well-known that the solution of (1) admits a multi-symplectic conservation law (MSCL) with the differential 2-form and where matrices and satisfy A numerical method that satisfies a discrete version of (2) is called multi-symplectic integrator; see [2] for reviews of multi-symplectic integration. Up to now, some multi-symplectic integrators have been proposed for the multi-symplectic Hamiltonian PDEs and exhibit good numerical performances, such as the box scheme obtained by applying simply implicit midpoint rule (a Runge–Kutta method) in both time and space [2], the Euler-box scheme developed by using the non-compact Euler rule in both time and space [10], the pseudospectral scheme obtained by applying the pseudospectral method in space and the implicit midpoint rule in time [11], and the diamond scheme [12]. Although multi-symplecticity by itself does not ensure good performance on traditional criteria like accuracy and stability, multi-symplectic methods do have a number of advantages. On the one hand, they are essentially variational and the standard discrete model in the relevant part of physics. On the other hand, the symplecticity in space is essential for a semidiscretization to be amenable to symplectic time integration, whose advantages in long-time integration are well known, and the symplecticity in space can also preserve periodic, quasi-periodic and heteroclinic solutions. However, multi-symplectic integrators can only preserve some conservation laws, and these do not typically include those of energy and momentum, which can be important for nonlinear stability and the nonlinear convergence analysis of the schemes.
In recent years, there has been an increased emphasis on constructing numerical methods preserving certain invariant quantities such as energy and momentum in the continuous dynamical systems. Nowadays, for computing physical/chemical problem, whether the numerical scheme captures the energy/momentum of the system has been a criterion to judge the success of the numerical simulation. Many PDEs with appropriate boundary conditions admit the global energy/momentum conservation law and there have been some techniques for constructing the global energy-/momentum-preserving methods for these conservative PDEs [[13], [14], [15]]. However, these discrete preserving properties are global and depend on the boundary conditions inevitably, so that the schemes are invalid for the problem without appropriate boundary conditions. As we know, besides possessing the global conservation laws, many PDEs also have local conservation laws which are independent of the boundary conditions and more essential than the global ones. For example, the multi-symplectic Hamiltonian PDE admits a local energy conservation law (LECL) where and represent the energy density and energy fluxes respectively, and also has local momentum conservation laws (LMCLs) in the directions , , given by where , and . It is clear that the local conservation laws (4), (5) are established on any time–space point and independent of boundary conditions. Assume that the PDEs are imposed on some appropriate boundary conditions, these LECL and LMCLs will result in global energy and momentum conservation laws, respectively. And thus the LECL and LMCLs produce richer informations of the PDE system than the corresponding global ones. Naturally, we expect to propose a scheme that satisfies a discrete version of LECL (4) or LMCL (5). In this paper, this kind of schemes will be called local energy-preserving (LEP) integrator or local momentum-preserving (LMP) integrator. On the construction of the LEM and LMP integrators for the PDEs, we have done some works. In [[16], [17]], by using the concatenating method, we have proposed some LEP and LMP schemes for the Sine-Gordon equation and the coupled NLS equations. These schemes have excellent properties and exhibit good numerical performance, but they are not completely systematic either in their derivations or in their applicability to a class of PDEs. Recently, combining the implicit midpoint rule with the averaged vector field (AVF) method, we have developed a LEP integrator and a LMP integrator which are applicable to the entire class (1) in one dimension [18]. The derivation of the LEP and LMP schemes are easy, however, these schemes have some less positive features. The schemes are fully implicit, which makes them expensive; The implicit equations may not have a solution: with periodic boundary conditions, solvability requires that the number of grid points be odd [19]; For a given PDE, it is difficult to establish the error estimates of the obtained LEP/LMP scheme no matter whether there exist positive invariants or not. The reason is that there is no relationship , where constants , , represents the discrete -norm, is an average operator, and represents the error vector whose entries denote the difference between the numerical solution of the LEP/LMP scheme and the exact solution of original problem.
In this paper, instead of discretizing one particular PDE, we wish to develop methods that are applicable to the entire class of multi-symplectic Hamiltonian PDEs (1), specializing to a particular equation or family as late as possible. To this end, we split the skew-symmetric matrix or in the form of (3). And then applying the leap-frog rule to space/time gives a semi-discrete system of ODEs with respect to time/space. Next, we use a discrete gradient (DG) method such as the coordinate increment DG method [20], the midpoint DG method [21], or the AVF method to integrate the obtained ODEs. In this study, an application of the AVF method gives a fully discrete LEP/LMP scheme for the multi-symplectic Hamiltonian PDE (1). The AVF method was first written down in [22] for the system of ODEs , , and identified as an energy-preserving and B-series method in [23]. The second-order AVF method is the map defined by where is the step size. The method (6) is affine-covariant [24] and self-adjoint. When is Hamiltonian with respect to a constant symplectic structure, i.e., with a nonsingular and antisymmetric matrix, the AVF method preserves the Hamiltonian . Besides preserving the LECL or LMCL of the system (1), the present method also has several appealing properties. It is systematic and applicable to a huge of PDEs which has the multi-symplectic Hamiltonian form (1); Applications of the matrix splitting technique and leap-frog rule may result in a fully explicit or linearly implicit scheme for some nonlinear PDEs; Different schemes can be proposed due to different splitting matrices; The method is effective and valid no matter whether the number of the grid points is odd or even; The present method is established on the grid points unlike the method in [18] establishing on the central box, so that the discrete conservation laws are also related to the grid points. For some particular multi-symplectic Hamiltonian PDEs holding the positive discrete conservation laws, the feature allows us to establish the error estimate of the scheme easily.
The present work is organized as follows. In Section 2, we briefly review the multi-symplectic form of the one-dimensional PDE and then develop a class of energy- and momentum-preserving integrators for it. We apply the developed local structure-preserving integrators to some classical one-dimensional PDEs and conduct numerical experiments to show their numerical performance in the same section. In Section 3, the methodology of the localstructure-preserving algorithms for the one-dimensional multi-symplectic Hamiltonian PDE is extended to multi-dimensional multi-symplectic Hamiltonian PDEs. In Section 4, we study the LEP discretization by means of numerical dispersion relation for the linear PDE. Finally, we finish the paper with concluding remarks in Section 5.
Section snippets
Local structure-preserving algorithms for 1D problem
Consider the multi-symplectic Hamiltonian PDE (1) in one dimension (), i.e., where . When dealing with (semi-)discrete systems we use the notation where the index corresponds to increments in space and to increments in time. That is, the point is the discrete equivalent of , where and . Define the finite difference operators , , and the average operators ,
Local structure-preserving methods for multi-dimensional problem and applications
Many PDEs such as two-/three-dimensional (2D/3D) Gross–Pitaevskii equation and Maxwell’s equations can be recast into the form of (1). The present LEP and LMP integrators for the 1D multi-symplectic Hamiltonian PDE (7) can be extended to -dimensional case (1) directly.
Let be a numerical approximation to where , , and . For the simplicity, in the most of the equations we present, one of the indices is held constant, in which case, we
Numerical dispersion
As one would expect, the nonlinearity of (1) makes the study of numerical solutions for the problem more difficult. Hence, the aim of this section is to conduct an analysis of the LEP discretization method for the linear equation where is symmetric. Applying the LEP scheme to (39) yields In the following, we only consider the choice of .
Theorem 4.1 Numerical solutions of the scheme (40) satisfy the numerical dispersion relation
Conclusion
Many PDEs such as the KdV equation, G–P equation, Maxwell’s equations and so on can be reformed into the general multi-symplectic Hamiltonian system which has three local conservation laws, including MSCL, LECL and LMCL. These conservation laws are independent of boundary conditions, so that they are more essential than the global ones. Naturally, it is desired to develop a numerical integrator holding one or more discrete versions of the local conservation laws for the general multi-symplectic
Acknowledgments
The work wassupported by Natural Science Foundation of Jiangsu Province of China (BK20181482), NNSF of China (11771213), the National Key Research and Development Project of China (2016YFC0600310), the Priority Academic Program Development of Jiangsu Higher Education Institutions of China and Jiangsu Overseas Visiting Scholar Program for University Prominent Young & Middle-aged Teachers and Presidents , China. The authors would like to express sincere gratitude to the referee for the
References (35)
- et al.
Phys. Lett. A
(2001) - et al.
Appl. Numer. Math.
(2004) - et al.
J. Comput. Phys.
(2008) - et al.
Comput. Math. Appl.
(2002) - et al.
J. Comput. Phys.
(2018) - et al.
J. Comput. Phys.
(2010) - et al.
J. Comput. Phys.
(2012) - et al.
J. Comput. Phys.
(2013) - et al.
J. Comput. Phys.
(2014) - et al.
Comput. Methods Appl. Mech. Engrg.
(1996)
Appl. Math. Comput.
Appl. Math. Model.
Appl. Math. Comput.
Comput. Phys. Comm.
Comput. Phys. Comm.
J. Comput. Appl. Math.
Cited by (11)
Global energy preserving model reduction for multi-symplectic PDEs
2023, Applied Mathematics and ComputationCitation Excerpt :In the last two decades, many multi-symplectic integrators have been developed; the box/Preissmann scheme [1], the Euler-box scheme [11], the Fourier pseudo-spectral collocation scheme [10], wavelet collocation method [12]. Local and global energy/momentum preserving methods have been developed in [2,4] with the average vector field (AVF) method and in [3,13] with the linearly implicit Kahan’s method. In this paper, a global ROM is investigated for multi-symplectic PDEs that preserve the energy/Hamiltonian globally.
Two classes of linearly implicit local energy-preserving approach for general multi-symplectic Hamiltonian PDEs
2020, Journal of Computational PhysicsCitation Excerpt :In this section, we apply our schemes to solve these equations with different boundary conditions. Finally, we make numerical comparisons between our schemes and the existing second-order finite difference schemes, such as the LEP schemes [12,18], multi-symplectic Preissmann scheme (MPS) [42] and narrow box scheme (NBS) [1]. The results are listed in Table 2.
Efficient mass- and energy-preserving schemes for the coupled nonlinear Schrödinger–Boussinesq system
2019, Applied Mathematics LettersEFFICIENT LINEARIZED LOCAL ENERGY-PRESERVING METHOD FOR THE KADOMTSEV-PETVIASHVILI EQUATION
2022, Discrete and Continuous Dynamical Systems - Series BLinearly implicit local energy-preserving algorithm for a class of multi-symplectic Hamiltonian PDEs
2022, Computational and Applied Mathematics