`Journal of Applied MathematicsVolume 2012 (2012), Article ID 705179, 22 pageshttp://dx.doi.org/10.1155/2012/705179`
Research Article

## Nonstandard Finite Difference Variational Integrators for Multisymplectic PDEs

Department of Mathematics, Harbin Institute of Technology, 2 Wenhua West Road, Shandong, Weihai 264209, China

Received 12 August 2012; Accepted 4 October 2012

Copyright © 2012 Cuicui Liao and Xiaohua Ding. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We use the idea of nonstandard finite difference methods to derive the discrete variational integrators for multisymplectic PDEs. We obtain a nonstandard finite difference variational integrator for linear wave equation with a triangle discretization and two nonstandard finite difference variational integrators for the nonlinear Klein-Gordon equation with a triangle discretization and a square discretization, respectively. These methods are naturally multisymplectic. Their discrete multisymplectic structures are presented by the multisymplectic form formulas. The convergence of the discretization schemes is discussed. The effectiveness and efficiency of the proposed methods are verified by the numerical experiments.

#### 1. Introduction

It is a fundamental approach to develop the discrete multisymplectic numerical methods based on the discrete Hamilton’s principle, because it leads in a natural way to multisymplectic integrators [1]. The discrete Euler-Lagrange equation is produced in the discrete variational principle [24]; meanwhile, the discrete multisymplectic structure is also generated [5, 6]. In the other words, the discrete variational integrators are multisymplectic automatically.

##### 1.1. Multisymplectic Structure of Discrete Variational Integrators

By the Hamilton’s principle [69], the discrete multisymplectic structure which is preserved by the discrete variational integrator, is described by Poincaré-Cartan forms, in a differential geometric language. In paper [6], Marsden et al. showed how to obtain this structure directly from the variational principle, on the Lagrangian side. They defined it as the multisymplectic form formula, and they showed that it was conserved by the discrete variational integrator.

Lemma 1.1. If is a solution of discrete Euler-Lagrange equation and , are first variations of , then the following discrete multsisymplectic form formula holds:

The details of this conclusion could be referred to papers [5, 6]. This conclusion states that the discrete variational principles produce discrete variational integrators, and the multisymplecticity of these variational integrators is presented by the discrete multisymplectic form formula (1.1).

Vankerschaver et al. [10] revisited the multisymplectic form formula [6], showing that it could be obtained from the boundary Lagrangian that they defined in their paper. They presented an easy way to derive discrete multisymplectic form formula from discrete variational principle, using the notations of Poincaré-Cartan forms. In this paper, we follow the same way to derive the discrete multisymplectic form formulas of our discrete variational integrators.

When we use the discrete variational principle, we need to make a approximation of the Lagrangian. Here, in our paper, we would use nonstandard finite difference methods, instead of standard finite difference, to approximate the Lagrangian function, and derive the corresponding discrete variational integrators.

##### 1.2. Nonstandard Finite Difference Methods

The nonstandard finite difference schemes are well developed by Mickens [1115] in the past decades. These schemes are developed for compensating the weaknesses that may be caused by standard finite difference methods, for example, the numerical instabilities. Regarding the positivity of solutions, boundedness, and monotonicity of solutions, nonstandard finite difference schemes have a better performance than standard finite difference schemes, due to its flexibility to construct a nonstandard finite difference scheme that can preserve certain properties and structures, which are obeyed by the original equations. Also, the dynamic consistency could be presented well by nonstandard finite difference scheme. These advantages of nonstandard finite difference methods have been shown in many numerical applications. GonzLez-Parra et al. [1618] developed nonstandard finite difference methods to solve population or biological models. The positivity condition and the conservation law of population dynamics are preserved by nonstandard finite difference schemes. Jordan [19] and Malek [20] constructed nonstandard finite difference schemes for heat transfer problems. For the symplectic systems, Mickens [15] derived the nonstandard finite difference variational integrator for symplectic ODEs. Ma et al. [21] developed the nonstandard finite difference variational integrator in stochastic ordinary differential equations.

The initial foundation of nonstandard finite difference methods is formed by the exact finite difference schemes [22]. After generalizing these results, Mickens summarizes the following three basic rules to construct nonstandard finite difference schemes. (1) The orders of the discrete derivatives should be equal to the orders of the corresponding derivatives appearing in the differential equations.Note. If the orders of the discrete derivatives are larger than those occurring in the differential equations, then numerical instabilities will in general occur.(2) Discrete representations for derivatives must, in general, have nontrivial denominator functions.Note. For example, the discrete first-derivative is generally represented by where the numerator functions and the denominator functions satisfy (3) Both linear and nonlinear terms should be represented by nonlocal discrete representations on the discrete computational lattice.Note. For example,

In our paper, we combine the advantages of nonstandard finite difference methods and discrete variational principles to construct nonstandard finite difference variational integrators, for two multisymplectic PDEs. These integrators are multysimplectic and their multysimplecticity are presented by their discrete multisymplectic form formulas, respectively.

In Section 2, we consider a simple linear wave equation. With the triangle discretization, we define the discrete Lagrangian using the idea of nonstandard finite difference and derive discrete variational integrator and the corresponding multisymplectic form formula, by discrete variational principle. The convergence of this method is analyzed. In Section 3, for the nonlinear Klein-Gordon equation, triangle discretization and square discretization are considered to obtain the nonstandard finite difference variational integrators. The discrete multisymplectic structures are presented, respectively. The convergence orders of the two methods are also discussed, and the convergence orders are shown in error tables in the numerical experiment section. Section 4 is devoted to showing the numerical behaviors of the developed nonstandard finite difference variational integrators.

#### 2. Nonstandard Finite Difference Variational Integrator for Linear Wave Equation

We first consider a simple linear wave equation, where is a scalar field function with two independent variables, and .

This linear wave equation is actually a multisymplectic PDE. As a classical and simple multisymplectic example, wave equation and its multisymplectic structure have been studied from both Hamiltonian [2328] and Lagrangian viewpoints [3, 5, 6, 10]. Based on the Lagrangian viewpoint, we could also obtain (2.1) from the Euler-Lagrange equation with the Lagrangian function ,

Assume that we have a uniform quadrangular mesh in the base space, with mesh lengths and . The nodes in this mesh are denoted by , corresponding to the points in . We denote the value of the field at the node by . We label the triangle at with three ordered triple as , and we define to be the set of all such triangles. Then the discrete jet bundle [6, 10] is defined as follows: which is equal to .

Now we use nonstandard finite difference to define the discrete Lagrangian on , which is the discrete version of Lagrangian density [10, 29], where denominator functions and are defined according the exact solution of wave equation [5, 6, 12]

We have followed the rules of constructing nonstandard finite difference schemes in Mickens’ papers [1115]. (1)The discrete first-derivative is represented by where denominator functions , are defined in (2.6). Using Taylor series expansion, Then the denominator functions satisfy (2)Nonlocal representation on the discrete computational lattice are used here by

So, for the linear wave equation (2.1) with the Lagrangian (2.3), the discrete Lagrangian becomes

By the discrete Hamilton’s principle [6, 10], we have the discrete Euler-Lagrange equation, where and are defined similarly as (2.11), which are

After some simple calculations, the discrete Euler-Lagrange equation (2.12) becomes

We could find that this scheme is symmetric in and , and . This is the nonstandard finite difference variational integrator, for the linear wave equation.

As we mentioned in Section 1 and Lemma 1.1, the advantages of deriving the multisymplectic numerical schemes from discrete variational principle are that they are naturally multisymplectic and the discrete multisymplectic structures are also generated in the variational principle. Now it is meaningful to show the multisymplectic structure of this discrete variational integrator (2.14) based on nonstandard finite difference method.

Since we consider triangulation discretization here, we focus on three adjacent triangles around , denote this area by . Following the idea in [10], the discrete boundary Lagrangian is given by where

Taking twice exterior derivative of both sides, we have, by the fact that , the discrete multisymplectic form formula with following form [10]: where (for ). The discrete Poincaré-Cartan forms , , and are defined by and similarly for and . Thus, for the linear wave equation (2.1), the multisymplectic form formula of this scheme (2.14) based on nonstandard finite difference method can be obtained as follows:

Now we have the first conclusion.

Theorem 2.1. The nonstandard finite difference variational integrator (2.14), for linear wave equation (2.1) is multisymplectic, and the discrete multisymplectic structure is

We now discuss the convergence of this variational integrator (2.14) based on the nonstandard finite difference method. From the Lax equivalence theorem we know that, for a well-posed linear initial value problem, the consistent finite difference method is convergent if and only if it is stable.

By Taylor series expansion, we have

Similarly,

The above two equations show that the scheme is consistent and the truncation error for the integrator (2.14) is .

To explore the stability of the nonstandard finite difference variational integrator (2.14), we introduce the following notations:

Then the three-level explicit integrator (2.14) is equivalent to the following two-level scheme:

By using the Von Neumann method [30], we could get the amplification matrix of the above scheme, where . Note that, in the above matrix, . Let . We have the characteristic equation and the eigenvalues

When ,, the scheme (2.14) satisfies the Von Neumann conditions, which is a necessary condition of the stability of the scheme (2.14). If , , where is an integer, then has two different eigenvalues. If , , then is an identity matrix, but has two different eigenvalues [30]. So is the sufficient condition of the stability for integrator (2.14). Note that, if , there is an unbounded solution , . So the scheme (2.14) is not stable when . Now, we find the necessary and sufficient condition of the stability for integrator (2.14), which is

With the consistence and stability conditions, we have following conclusion.

Theorem 2.2. The nonstandard finite difference variational integrator (2.14) is convergent, when the step sizes and satisfy .

We have shown the idea of using the nonstandard finite difference method to get the discrete variational integrator and the corresponding discrete multisymplectic form formula. In the next section, we will consider the discrete variational integrators for a more complicated example, the nonlinear Klein-Gordon equation.

#### 3. Nonstandard Finite Difference Variational Integrators for Nonlinear Klein-Gordon Equation

In this section, we consider the nonlinear Klein-Gordon equation [31],

As known, this equation can be obtained by Euler-Lagrange equation (2.2) with the Lagrangian function

Now we consider the triangle discretization and square discretization, respectively, to get the nonstandard finite difference variational integrators.

##### 3.1. Triangle Discretization

Following the steps in last section and using the idea of nonstandard finite difference, we define the discrete Lagrangian as based on the following constructing rules, (1)The discrete first derivative is represented by where the denominator functions are defined by (2.6), and (2)Nonlocal representations for and are given by where , , , and are positive parameters. Such discretizations for and guarantee the symmetric property of the discrete Lagrangian function [15].

Similarly, we define discrete Lagrangians on other two adjoint triangles,

Now, the discrete variational integrator with nonstandard finite difference methods could be obtained by discrete Euler-Lagrange equation (2.12):

Substituting , , and into above equation, we arrive at

Using the definition of discrete Lagrangian functions, one can find that this scheme is symmetric with respect to and , and ; that it is multisymplectic, and that it preserves the multisymplectic structure of the original equation.

Its corresponding discrete multisymplectic form formula can be obtained from (2.17), that is, where

It shows the multisymplectic structure of scheme (3.9), and the relations between the field values on the three adjoint triangles are around .

We now analyze the truncation error of integrator (3.9). By Taylor series expansion [32, 33], we have

Combining the above two equations and (2.22), (2.23), we can observe that the nonstandard finite difference variational integrator (3.9) has the truncation error .

The above results are summarized in the following theorem.

Theorem 3.1. The nonstandard finite difference variational integrator (3.9) for the nonlinear Klein-Gordon equation (3.1) is multisymplectic, and its truncation error is . The discrete multisymplectic structure of this scheme is presented by (3.10).

##### 3.2. Square Discretization

In this case, we denote a square at with four ordered quaternion by and define to be the set of all such squares. Then the discrete jet bundle [6, 10] is defined as which is equal to .

Following the philosophy of the nonstandard finite difference method, we define the discrete Lagrangian on as

In this case,(1)the discrete first-derivative is represented by where the denominator functions are defined by (2.6), and (2)nonlocal representations for and are

Similarly, we have the definitions of on the other three squares adjoint to :

Taking derivate of action functional with respect to , we have the discrete Euler-Lagrange equation in this square discretization [5, 6, 10, 34], which is

Substituting the discrete Lagrangian , , , and into the previous equation, we arrive at

After simple calculations, it becomes

It is multisymplectic and symmetric in and , and . Similarly, we have the discrete multisymplectic form formula:

To study the truncation error of the integrator, we do the Taylor expansion which leads to

Combing these equations, we can readily observe that the nonstandard finite difference variational integrator (3.21) has truncation error . To verify this conclusion, we investigate the numerical convergence order in our numerical experiments. See Section 4.

We summarize our conclusion in the following theorem.

Theorem 3.2. The nonstandard finite difference variational integrator (3.21) for the nonlinear Klein-Gordon equation (3.1) is multisymplectic, and its truncation error is . The discrete multisymplectic structure of this scheme is presented by (3.22).

#### 4. Numerical Simulations

In this section, we report the performance of the nonstandard finite difference variational integrator (2.14) for solving linear wave equation (2.1) and the nonstandard finite difference variational integrators (3.21) and (3.9) for the nonlinear Klein-Gordon equation (3.1).

##### 4.1. Linear Wave Equation

For linear wave equation (2.1), we consider the initial conditions and the periodic boundary conditions

The nonstandard finite difference variational integrator (2.14) is an explicit five points scheme. We choose the denominator functions and in as and .

From Figure 1, we can see that the nonstandard finite difference variational integrator (2.14) for the linear wave equation performs very well and the periodicity of the linear wave equation is preserved accurately.

Figure 1: The waveforms of linear wave equation (2.1) by the nonstandard finite difference variational integrator (2.14) , .
##### 4.2. Nonlinear Klein-Gordon Equation

We now consider the nonlinear Klein-Gordon equation (3.1) with the initial condition and periodic boundary conditions.

We use nonstandard finite difference variational integrator (3.21) to simulate this problem with amplitude . The nonstandard finite difference variational integrator (2.14) is an implicit nine-points nonstandard finite difference scheme. The denominator functions and are defined the same as before.

As depicted in Figure 2, the nonstandard finite difference variational integrator (3.21) simulates the wave propagation perfectly at the beginning. After a long time simulation, the integrator still performs very accurate and stable, without showing any blowup. With periodic boundary condition, the wave going out the computational domain shows up in the other direction periodically.

Figure 2: Waveforms of nonlinear Klein-Gordon equation by the nonstandard finite difference variational integrator (3.21) , .
##### 4.3. Convergence Order of the Nonlinear Integrators (3.9) and (3.21)

To further investigate the numerical convergence of the proposed schemes, we conduct a series of numerical tests of our nonlinear integrators. In this example, we consider the nonlinear Klein-Gordon equation (3.1) with the initial boundary conditions as follows: where . The exact solution of the problem is

The nonstandard finite difference integrators (3.9) and (3.21) are applied to simulate the Klein-Gordon equation. In the integrator (3.9), we choose here. The -norm errors at , , and are listed in Tables 1 and 2. The orders in the tables are calculated with the formula [35, 36]

Table 1: The error and convergence orders of (3.9) for Klein-Gordon problem with .
Table 2: The error and convergence orders of (3.21) for Klein-Gordon problem with .

Overall, it is clear that the error decreases as the mesh size goes to zero, indicating the convergence of our nonlinear integrators (3.9) and (3.21). Moreover, numerical orders clearly exhibit second order convergence when the mesh size decreases with , which further confirms our theoretical derivation of the truncation errors of the numerical schemes.

In all, the numerical tests verify that the nonstandard finite difference variational integrators that we developed are capable of preserving the characteristics of original equations. They are all accurate, effective, and suitable for solving multisymplectic systems.

#### 5. Conclusion

In this paper, we have considered a linear wave equation and a nonlinear Klein-gordon equation. We have derived the nonstandard finite difference variational integrators and the corresponding multisymplectic form formulas from these two multisymplectic PDEs. We have shown that the nonstandard finite difference methods are flexible in constructing numerical schemes and can be employed to derive multisymiplectic schemes for multisymplectic systems. The convergence of our methods has been discussed. The numerical experiments have shown effectiveness and efficiency of these nonstandard finite difference variational integrators.

#### Acknowledgments

The authors are grateful to the editor and anonymous reviewers for their careful reading and many constructive suggestions which lead to a great improvement of this paper. This work is supported by the NNSF of China (no. 11271101) and the NNSF of Shandong Province (no. ZR2010AQ021).

#### References

1. R. Abraham and J. E. Marsden, Foundations of Mechanics, Benjamin/Cummings Publishing, Reading, Mass, USA, 2nd edition, 1978.
2. J.-B. Chen, “Variational integrators and the finite element method,” Applied Mathematics and Computation, vol. 196, no. 2, pp. 941–958, 2008.
3. J.-B. Chen, M.-Z. Qin, and R. Scherer, “Multisymplectic and variational integrators,” International Journal of Pure and Applied Mathematics, vol. 44, no. 4, pp. 509–536, 2008.
4. P. G. Ciarlet, A. Iserles, R. V. Kohn, and M. H. Wright, Simulating Hamltonian Dynamics, Cambridge Monographs on Applied and Computational Mathematics, Cambridge University Press, New York, NY, USA, 2004.
5. J. E. Marsden and M. West, “Discrete mechanics and variational integrators,” Acta Numerica, vol. 10, pp. 357–514, 2001.
6. J. E. Marsden, G. W. Patrick, and S. Shkoller, “Multisymplectic geometry, variational integrators, and nonlinear PDEs,” Communications in Mathematical Physics, vol. 199, no. 2, pp. 351–395, 1998.
7. M. Leok and J. Zhang, “Discrete Hamiltonian variational integrators,” IMA Journal of Numerical Analysis, vol. 31, no. 4, pp. 1497–1532, 2011.
8. K. I. Nawafleh, “Canonical quantization of higher-order Lagrangians,” Journal of Applied Mathematics, vol. 2011, Article ID 375838, 11 pages, 2011.
9. S. Pekarsky and J. E. Marsden, “Abstract mechanical connection and abelian reconstruction for almost Kähler manifolds,” Journal of Applied Mathematics, vol. 1, no. 1, pp. 1–28, 2001.
10. J. Vankerschaver, C. Liao, and M. Leok, “Generating functionals and Lagrangian PDEs,” Journal of Nonlinear Science. In press.
11. R. E. Mickens, Applications of Nonstandard Finite Difference Schemes, World Scientific Publishing, Singapore, 2000.
12. R. E. Mickens, “Nonstandard finite difference schemes for differential equations,” Journal of Difference Equations and Applications, vol. 8, no. 9, pp. 823–847, 2002.
13. R. E. Mickens, “A nonstandard finite difference scheme for the diffusionless Burgers equation with logistic reaction,” Mathematics and Computers in Simulation, vol. 62, no. 1-2, pp. 117–124, 2003.
14. R. E. Mickens, “Dynamic consistency: a fundamental principle for constructing nonstandard finite difference schemes for differential equations,” Journal of Difference Equations and Applications, vol. 11, no. 7, pp. 645–653, 2005.
15. R. E. Mickens, “A numerical integration technique for conservative oscillators combining nonstandard finite-difference methods with a Hamilton's principle,” Journal of Sound and Vibration, vol. 285, no. 1-2, pp. 477–482, 2005.
16. A. J. Arenas, G. González-Parra, and B. M. Chen-Charpentier, “A nonstandard numerical scheme of predictor-corrector type for epidemic models,” Computers & Mathematics with Applications, vol. 59, no. 12, pp. 3740–3749, 2010.
17. G. González-Parra, A. J. Arenas, and B. M. Chen-Charpentier, “Combination of nonstandard schemes and Richardson's extrapolation to improve the numerical solution of population models,” Mathematical and Computer Modelling, vol. 52, no. 7-8, pp. 1030–1036, 2010.
18. R. J. Villanueva, A. J. Arenas, and G. Gonzlez-Parra, “A nonstandard dynamically consistent numerical scheme applied to obesity dynamics,” Journal of Applied Mathematics, vol. 2008, Article ID 640154, 14 pages, 2008.
19. P. M. Jordan, “A nonstandard finite difference scheme for nonlinear heat transfer in a thin finite rod,” Journal of Difference Equations and Applications, vol. 9, no. 11, pp. 1015–1021, 2003.
20. A. Malek, “Applications of nonstandard finite difference methods to nonlinear heat transfer problems heat transfer-mathematical modelling,” Numerical Methods and Information Technology In press.
21. Q. Ma, D. Ding, and X. Ding, “A nonstandard fnite-difference method for a linear oscillator with additive noise,” Applied Mathematics and Information Science. In press.
22. P. M. Manning and G. F. Margrave, “Introduction to non-standard fnite-difference modelling,” CREWES Research Report 18, 2006.
23. T. J. Bridges, “Multi-symplectic structures and wave propagation,” Mathematical Proceedings of the Cambridge Philosophical Society, vol. 121, no. 1, pp. 147–190, 1997.
24. T. J. Bridges and S. Reich, “Numerical methods for Hamiltonian PDEs,” Journal of Physics A, vol. 39, no. 19, pp. 5287–5320, 2006.
25. T. J. Bridges and S. Reich, “Multi-symplectic integrators: numerical schemes for Hamiltonian PDEs that conserve symplecticity,” Physics Letters A, vol. 284, no. 4-5, pp. 184–193, 2001.
26. J. Hong and Y. Liu, “A novel numerical approach to simulating nonlinear Schrödinger equations with varying coefficients,” Applied Mathematics Letters, vol. 16, no. 5, pp. 759–765, 2003.
27. J. Hong, Y. Liu, H. Munthe-Kaas, and A. Zanna, “Globally conservative properties and error estimation of a multi-symplectic scheme for Schrödinger equations with variable coefficients,” Applied Numerical Mathematics, vol. 56, no. 6, pp. 814–843, 2006.
28. X.-C. Li and W.-A. Yao, “Symplectic analytical solutions for the magnetoelectroelastic solids plane problem in rectangular domain,” Journal of Applied Mathematics, vol. 2011, Article ID 165160, 15 pages, 2011.
29. J. Vankerschaver and M. Leok, “A variational, second-order integrator for point vortices on the sphere,” . In press.
30. T. Meis and U. Marcowitz, Numerical Solution of Partial Differential Equations, Springer, Berlin, Germany, 1981.
31. A. M. A. El-Sayed, A. Elsaid, and D. Hammad, “A reliable treatment of homotopy perturbation method for solving the nonlinear klein-gordon equation of arbitrary (fractional) orders,” Journal of Applied Mathematics, vol. 2012, Article ID 581481, 13 pages, 2012.
32. J.-B. Chen, “A multisymplectic integrator for the periodic nonlinear Schrödinger equation,” Applied Mathematics and Computation, vol. 170, no. 2, pp. 1394–1417, 2005.
33. J.-B. Chen, M.-Z. Qin, and Y.-F. Tang, “Symplectic and multi-symplectic methods for the nonlinear Schrödinger equation,” Computers & Mathematics with Applications, vol. 43, no. 8-9, pp. 1095–1106, 2002.
34. M. Leok and T. Shingel, “General techniques for constructing variational integrators,” Frontiers of Mathematics in China, vol. 7, no. 2, pp. 273–303, 2012.
35. S. Zhou and X. Cheng, “Numerical solution to coupled nonlinear Schrödinger equations on unbounded domains,” Mathematics and Computers in Simulation, vol. 80, no. 12, pp. 2362–2373, 2010.
36. S. Zhou and X. Cheng, “A linearly semi-implicit compact scheme for the Burgers-Huxley equation,” International Journal of Computer Mathematics, vol. 88, no. 4, pp. 795–804, 2011.