Hindawi Publishing Corporation Journal of Applied Mathematics Volume 2012, Article ID 705179,22 pages doi:10.1155/2012/705179

Research Article

Nonstandard Finite Difference Variational Integrators for Multisymplectic PDEs

Cuicui Liao and Xiaohua Ding

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

Correspondence should be addressed to Xiaohua Ding, mathdxh@hit.edu.cn Received 12 August 2012; Accepted 4 October 2012 Academic Editor: Renat Zhdanov

Copyright © 2012 C. Liao and X. 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.

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 [2-4]; 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 [6-9], the discrete multisymplectic structure which is preserved by the discrete variational integrator, is described by Poincare-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 u is a solution of discrete Euler-Lagrange equation and V, W are first variations of u, 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 mul-tisymplecticity 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 Poincare-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 [11-15] 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. [16-18] 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

du Mj+1 - y(At)u.i

dt <p(At)

where the numerator functions f (At) and the denominator functions $(At) satisfy ^(At) = 1 + O((At)), $(At) = At + 0( (At)2). (1.3)

(3) Both linear and nonlinear terms should be represented by nonlocal discrete representations on the discrete computational lattice.

Note. For example,

u —> 2ui - ui+1, 2

u —> uiui+1,

2 ( ui-1 + ui + ui+1 \ u2 ( -3- Wi, (1.4)

u3 —> 2ui3 - ui2ui+1,

u —> ui-1uiui+1.

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,

utt = uxx, (2.1)

where u(x, t) is a scalar field function with two independent variables, x and t.

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 [23-28] and Lagrangian viewpoints [3,5, 6,10]. Based on the Lagrangian viewpoint, we could also obtain (2.1) from the Euler-Lagrange equation

dL d dL + d dL (2 2)

du dt dut dx dux' '

with the Lagrangian function L(u, ut, ux),

L(u, ut,ux) = 2ut2 - 2Ux2. (2.3)

Assume that we have a uniform quadrangular mesh in the base space, with mesh lengths Ax and At. The nodes in this mesh are denoted by (i,j) e Z x Z, corresponding to the points (xi,tj) := (iAx,jAt) in R2. We denote the value of the field u at the node (i,j) by uj. We label the triangle at (i,j) with three ordered triple ((i,j), (i + 1,j), (i,j + 1)) as Aij, and we define XA to be the set of all such triangles. Then the discrete jet bundle [6,10] is defined as follows:

Ji r := { («W+j e R3 : ((i,j), (i + 1, j), (i,j + 1)) e Xt], (2.4)

which is equal to XA x R3.

Now we use nonstandard finite difference to define the discrete Lagrangian Ld on Jir, which is the discrete version of Lagrangian density [10, 29],

j j j+1 j+1 j j j := 1 A'AxL( "u^.^), (2*5)

where denominator functions $(At) and f (Ax) are defined according the exact solution of wave equation [5, 6,12]

0(At) = ±sm( f (Ax) = isi^ Ax)- (2.6)

We have followed the rules of constructing nonstandard finite difference schemes in Mickens'papers [11-15].

(1) The discrete first-derivative is represented by

. (2.7)

du ui - ui du u'i+1 - ui

dt $(At) dx f(Ax)

where denominator functions $(At), f (Ax) are defined in (2.6). Using Taylor series expansion,

At\ At 1

n-;=-- 48(At)3+•••. (2.8)

Then the denominator functions satisfy

#(At) = At + (At)3), f (Ax) - Ax + (Ax)3). (2.9)

(2) Nonlocal representation on the discrete computational lattice are used here by

j j j+1 ui + ui+1 + ui (2.10)

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

„ (".'"■'+1,"j+1) - 2 A'Ax( Y - i( fAffy (2.11)

By the discrete Hamilton's principle [6, 10], we have the discrete Euler-Lagrange equation,

D1Ld(uj,uj+1,uj+1) + D2L^uj-1,uj,uj+1) + D3Ld(u{-1,uj-1,uQ - 0, (2.12) where Ld(uj_ 1,uj,uj+1) and Ld(u ^uj+^ul) are defined similarly as (2.11), which are

jj) = 2 A'Axi 2 ( jAh )-1 ( jé,

(2.13)

Ld (u;1, u-1, uÙ ) = 1 AtAxf U { -j -

dV i ' i-1 ' V 2 ( 2 I $(At) ) 2 I f(Ax)

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

uj+1 - 2uj + uj 1 uj+ - 2uj + uj , i i i i+1 i i-1 = 0. (2.14)

(<£(Af))2 (f (Ax))2

We could find that this scheme is symmetric in (i,j + 1) and (i,j - 1), (i + 1,j) and (i - 1,j). 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 multis-ymplectic 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 uj, denote this area by U. Following the idea in [10], the discrete boundary Lagrangian is given by

Ldu(udu) := extu1[L^unl,unl+l,unl+l) + Ld^.^.u^1) + Ld(<-1,un+;11,<)], (2.15)

n+1 ,,n+1

Udu := («¿+1,«r 'ui-1,uï-1'uï

-1,«n+-i1

(2.16)

Taking twice exterior derivative of both sides, we have, by the fact that d2LdU = 0, the discrete multisymplectic form formula with following form [10]:

X E A®)= 0, (2.17)

k=1 l=1;l= k

where ^L = -d©L (for k = 1,2,3). The discrete Poincare-Cartan forms ©[, and ©L are defined by

©L(un,un+1,un+1) := D1L^un,un+1,un+^dunlf (2.18)

and similarly for ©L and ©3L. 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:

duj+1 A duj + duj 1 A duj duj+1 A dù + duj_ 1 A duj

№(Af))2 (f(Ax))2

= 0. (2.19)

Now we have the first conclusion.

Theorem 2.1. The nonstandard finite difference variational integrator (2.14),

M+1 - 2u\ + Wi 1 uL - 2W + W , 1 1 1 '+1 1 i-1 = 0, (2.20)

(<£(Af))2 (f (Ax))2

/or //near wave equation (2.1) is multisymplectic, and the discrete multisymplectic structure is

du{+l h duj

+ du', A du

du\+1 A dui

■ dui_ 1 A dui

(0(Af))

(f (Ax)) 2

(2.21)

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

W+1 - 2ui + uf1 _ w+1 - 2ui + uf1 _ 1 -O(Af)0 '+1 „ j '-1N

($(Af))2 (Af + o( (Af)3))2 (Af)2

(W+1 - 2u{ + W-1)

1 -O Af)2), . (2.22)

-^^ Ç(Af)2utt(xi,f j) + O(Af)4)

utt (Xi,f) + o( (Af)2).

Similarly,

uj+1 - 2uj + uj-1 (f (Ax))2

_ Uxx(xl,fj) + o((Ax)2). (2.23)

The above two equations show that the scheme is consistent and the truncation error for the integrator (2.14) is O((Af)2 + (Ax)2).

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

j j-1 uj uj vj _ Ui - Ui wj _ ui -Ui-1 (2.24)

Vi $(Af) , Wi-1/2 f(Af) . ' '

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

J+1 _ vj wj - wj Vi_VL = Wi+1/2 t-1/2

<(At) = w (Ax) .

•+i • • (2.25)

Wj-1/2 - wj-1/2 = j -<(At) w(Ax) .

By using the Von Neumann method [30], we could get the amplification matrix of the above scheme,

At) = ^2ir sinAx) 1 - 4r2sSd5Ax)). (2.26)

where r = <(At)/y(Ax). Note that, in the above matrix, i = V-1. Let n = 4r2sin2(^Ax). We have the characteristic equation

I2 - (2 - n)X + 1 = 0 (2.27)

and the eigenvalues

1 = 1 - 1 n ±(4n2 - ^1/2- (2-28)

When |1| < 1, that is,r < 1, the scheme (2.14) satisfies the Von Neumann conditions, which is a necessary condition of the stability of the scheme (2.14). If r < 1, f5h = nn, where n is an integer, then G has two different eigenvalues. If r < 1, f5h = nn, then G is an identity matrix, but (d/dph)G has two different eigenvalues [30]. So r < 1 is the sufficient condition of the stability for integrator (2.14). Note that, if r = 1, there is an unbounded solution vj = (-1)i+j(1 - 2j), wj+1/2 = (-1)i+j2j. So the scheme (2.14) is not stable when r = 1. Now, we find the necessary and sufficient condition of the stability for integrator (2.14), which is

r = < 1. (2.29)

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 At and Ax satisfy <(At) < w(Ax).

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],

Utt = uxx - u3 + u. (3.1)

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

L(u,ut,ux) = 2ut2 - 1 ux2 - 4u4 - 2u2. (3.2)

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 Ld as

Ld^w^wf1)

= 1 AtAxf Wj+1 - ^2 ^ j - ^2

2 \2\ $(At) J 2 V f(Ax)

( j4 j 4 j+14 \ u ( j2 j j+1 j 2 j j+1 j+12 j j an ui + ui+1 + ui 1+ bu u\ ui+1ui + ui+1 uiui + ui uiui+1

Il 3a1 + 3b1

y a^uj2 + uj+12 + uj+12) + b^ujt

■j2 + ^j 2 + ^ + b?(ujuj+1 + uj+1uj+1 + u{u{+1

2 ^ 3a2 + 3b2 y j

based on the following constructing rules,

(1) The discrete first derivative is represented by

j+1 j j j

. (3.4)

du uj - uj du uj+1 - u

dt $(At) ' dx f(Ax)

where the denominator functions are defined by (2.6), and

<(At) = At + ©( (At)3), w (Ax) = Ax + ©( (Ax)3). (3.5)

(2) Nonlocal representations for u4 and u2 are given by

I j4 j 4 j + 14\ 7 ( j2 j j + 1 j 2 j j + 1 j + 12 j j

aW u\ + ui+1 + u\ 1+ bu uiui+1ui + ui+1 uiui + u\ uiui+1

j2 j 2 j+12 a2 ui + ui+1 + ui

3a1 + 3b1

+ b2 uj

j j j+1 j j+1

u ui+1 + ui+1ui + u u

3a2 + 3b2

where a1, b1, a2, and b2 are positive parameters. Such discretizations for u4 and u2 guarantee the symmetric property of the discrete Lagrangian function [15].

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

Ld(uj_1,uj,j1)

=1 AtAxi \(uS_zukY_1/uj-uj-A2

= 2AtAxl ^ <(At) ) 2\ w(Ax) )

j4 j 4 j+14 j 2 j j+1 j2 j j+1 j+12 j j a1 ui + ui-1 + ui-1 + b1 ui-1 uiui-1 + ui ui-1ui-1 + ui-1 ui-1ui

3a1 + 3b1

/ ( j 2 j2 j + 12\ 7 / j j j j + 1 j + 1 j \

1 ( ui-1 + uj + uj-1 j + b2\ui-1ui + ui<1 + ui-1 ui-1)

3a2 + 3b2

i... 11fM-uj-1\ 1/i1 -uj-1V

= 2 AtAx

^ ty(At) J 2 y y(Ax) y

/ j-14 j-14 j4\ / j-12 j-1 j j-12 j j-1 j2 j-1 j-1

^ uj + u'i+1 + uj J + Mi Mi+1 Mi + ui+1 u\u\ + ui+1

3a1+ 3b1

j-12 j-12 j2 j-1 j-1 j-1 j j j-1 ^ uj + uj+1 + uj j + b2{uj uj+1 + uj+1 u\ + «X J

3a2 + 3b2

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

D1Ld(u{, uj+1, «J+1) + D2Ld(u4-1,uÙi, uj-1) + DsLd (uj-1, uj-J, j = 0. (3.8)

Substituting Ld (uj,uj+1,uj+1), Ld (uj_ pu^u^1), and Ld (uj 1,uj+11,uj) into above equation, we arrive at

uj+1 - 2uj + uj 1 uj+1 - 2uj + uj_ 1

(^(Ai))2 + (y (Ax))2

4 3a1 + 3b1

12a1uj + 2b^ujuj+1uj+1 + ujuj-1uj_+1 + ujuj 1uj+1)

, / i 2 j+1 j+12 j j 2 j+1 j+12 j j-12 j-1 j-12 j-1Y + ui+1 ui + uj uj+1 + ui-1 uh + uj-1 uh + uj uj+1 + ui+1 ui )

\3ai + 3b2 [^A2uj +b2 ( uj+1+«T1+uj-1+um + +=0.

Using the definition of discrete Lagrangian functions, one can find that this scheme is symmetric with respect to (i,j +1) and (i,j -1), (i + 1,j) and (i- 1,j); 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,

d^+1 A dui - dui A d^du\+l A dui - dui A du\_1 ($m)2 + (f(Ax))2

C1 uiui+1 + C1u\+1u\+1 + 2C1 u[+1 - C^jdu\+1 A dui + ^C1u[u[+1 + 1 C1uj+1 + C1uj+1uj+1 - C^j du\+l A duj + ^C1u\u\++l + C1u(-1u(+11 + 1 C1u{+11 - C2^ du(-1 A du[ (3.10)

+ ^C1ujuj-1 + 2C1uj-1 + Quj-^uj^ - C2^duj-1 A duj + ^C1ujuj+11 + C1uj 1uj+11 + 1 C1uj+11 - C2^duj 1 A duj + ^ C1ujuj 1 + 2 C1uj 1 + C1uj+]1 uj 1 - C^j dul+1 A duj^ = 0,

Cl 2(3ai + 3ki) ' C 3a2 + 3fr2. (3.11)

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

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

4 3a + 3b

12a1uj + 2b1(u[u\+l uj+1 + ujuj-1uj-11 + ujuj 1uj+l)

, ( j 2 j+1 j+12 j j 2 j+1 j+12 j j-12 j-1 j-12 j-A + b1{ uj+1 uj + u uj+1 + uj- 1 uj- 1 + uj- 1 uj- 1 + uj uj+1 + uj+1 uj )

u3 (xi,tj ) + ©( (Ax)2 + AxAt + (At)2),

l6fl2uj + uj+1 + uj+1 + uj-1 + uj-+i + uj-i + uj-1)l

u(xi,tj ) + (Ax)2 + (At)2).

(3.12)

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 O((Ax)2 + AxAt + (At)2).

The above results are summarized in the following theorem.

Theorem 3.1. The nonstandard finite difference variational integrator (3.9) for the nonlinear KleinGordon equation (3.1) is multisymplectic, and its truncation error is O((Ax)2 + AxAt + (At)2). The discrete multisymplectic structure of this scheme is presented by (3.10).

3.2. Square Discretization

In this case, we denote a square at (i,j) with four ordered quaternion ((i,j), (i + 1,j), (i + 1,j + 1), and (i,j + 1)) by Dj and define XQ to be the set of all such squares. Then the discrete jet bundle [6,10] is defined as

jQ Y := { («! ,Mj;i1,Mj+^ e R4 : ((i,j), (i + 1,j), (i + 1,j + 1), (i,j + 1)) e Xn}, (3.13)

which is equal to XQ x R4.

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

Ldfauh ,u{++\,wj+1)

+1 -ui + uâzUiiA A(-ui+1 + Ui+1-_ 1 ujuj u+u+i

2\24>(At) ty(At) J A 2f(Ax) 2y(AxW 4uiui+iui+1 u

- / i j i j+i i j+i i j+i j+i j+i j+i j+i \ \

1 / u • u. 1 + u• u. 1 + u • u • + u. 1 u. 1 + u • 1 u • + u • ,u' \ \

± I l l+i l l+i ^ ^t+i^ l + i l+i l l+i l » ll,!.,

2 y-6-y yAtAx-

(3.i4)

In this case,

(i) the discrete first-derivative is represented by

j+i j j+i j

du u\ - u\ uj+i - uj+i

dt 2è(At) 2è(At) '

. 1 . . (3.i5)

, j + i j+i j j

du ^ uj+i - uj + uj+i - uj

dx 2y(Ax) 2y(Ax) '

where the denominator functions are defined by (2.6), and

#(At) = At + 0( (At)3), y (Ax) = Ax + o( (Ax)3);

(3.i6)

14 Journal of Applied Mathematics

(2) nonlocal representations for u4 and u2 are

j j j+1 j+1 u u u u ,

I 1+1 1+1 I '

j uj+1 li+1ui+1

j+1 j+1 j+1 j+1 uj+1 u\ + uj+1 ui

(3.17)

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

Ld(uj-1,uj',uj'+1

- i j+1 j j+1 j 1 / ui - ui ui-1 - ui-1

- / j j j+1 j+1' 1 / ui - ui-1 ui - ui-1

2 V ty(At) fy(At) J 2 V 2f(Ax) 2f(Ax)

1 j j j+1 j+1 - 4ui-1u\uj uj-1

„ / j j j j+1 j j+1 j j+1 j j+1 j+1 j+1'

1 / ui_ 1ui + ui_ 1ui + ui_ 1ui_ 1 + uiui + uiui_ 1 + ui ui_ 1

A tA x,

Ld(u\-\,u\ \uj,uj-1)

j j-1 j j-1 1 ui 1 - ui 1 ui - ui

If uj - u^ 1 ui - ui_ 1

j-1 j-1

2 V 20(At) 2fy(At) 2 V 2^(Ax) 2^(Ax)

1 j-1 j-1 j j - 4 uj-1uj ujuj-1

j-1 j-1 j-1 j j-1 j j-1 j j-1 j j j 1 / ui_ 1 ui + u\_ 1 ui + u\_ 1 ui_ 1 + ui ui + ui ui_ 1 + uiui_ 1

I AtAx,

, / j j-1 j j-1 \ 2 , / j j j-1 j-1 \ 2 1 uj - uj . uj+1 - uj+1 \ 1 / uj+1 - uj . uj+1 - uj \

2 V 2#(At) ty(At) J 2 V 2f(Ax) 2f(Ax)

1 j-1 j-1 j j —u- ui ,ui 4 i i+1 i+1 i

j-1 j-1 j-1 j j-1 j j-1 j j-1 j j j 1 / uj uj+1 + u\ uj+1 + u\ uj + uj+1 uj+1 + uj+1 u\ + uj+1 uj

(3.18)

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

D1Ld(u1i,u1i+1, u^l wf1) +

(3.19)

+ D3Ld(ul-l, uiu\, + D4Ld(u{ ~1,u'i-L ,u\+1, uQ = 0.

Substituting the discrete Lagrangian Ld(uj,uj+1,uj+1,uj+1), Ld(uj_1,uj,uj+1,uj+1), Ld(uj_1,uj l,u\,u{_ 1), and Ld(uj 1,uj+11,uj+1,uj) into the previous equation, we arrive at

("- + 6) (uj+1 + ^ + (+ + 6) (uj+1 + ^ / 1 1 \ j / 1 1 1 \ / j+1 j+1 j-1 ;-1\

+ -7--7 K +--T +-T + 77 ui+1 + ui 1 + ui 1 + u' )

V^(Ai)2 y(Ax)2/ i V 4^(Ai)2 4y(Ax)2 17 J V i+1 i-1 i-1 i+1^

1 j j+1 j+1 j j+1 j+1 j-1 j-1 j j-1 j-1 j - 4 (uj+1uj+1uj + u'i-1u'i uj-1 + uj-1uj uh + uj uj+1 uj+^ = 0-

(3.20)

After simple calculations, it becomes

1 / uj++1 ~ 2uj+1 + + 2 uj+1 - 2uj + uj-1 + ui+1 ~ 2uj-1 + «j-1 4 (ф(А0)2 (ф(А0)2 (ф(А0)2

\( uj++11 - 2uj+1 + ui-+1 + 2 uj+1 - 2uj + uj-1 + ui+-1 - 2uj-1 + ui-11 \ 4 V (<KA*))2 (f(Ax))2 (f (Ax))2 J (3.21)

1 / i i+1 i+1 i i+1 i+1 i-1 i-1 i i-1 i-1 i \ 4[u'i+1 ui+1 ui + ui- 1ui ui-1+ ui- 1 ui ui- 1 + ui ui+1 ui+J

Ï2 (ui++i + + uri + ui+-1 + 2ui+1 + 2ui-1 + 2ui+1 + 2uM) = 0.

It is multisymplectic and symmetric in (i,j + 1) and (i,j — 1), (i + 1,j) and (i - 1,j). Similarly, we have the discrete multisymplectic form formula:

I 1 1 11 j j + 1 1 j j + M J j + 1 л J j

(--г--т + 7 — 7 ui+1ui+1 — — ui 1 )ам; л au;

V 2(^(Af))2 2(y(Ax))2 6 4 i+1 i+1 4 i—1 i—7 i i

/1 1 11 j—1 j+1 1 j—1 j \ j—1 j

+--т--т + 7 — 7 M 1 m; — 7 ui+i ам л dM

V 2(ф(А0)2 2(y(Ax))2 6 4 i—1 ; 4 ^ / ; ;

1 1 1 1 j+1 j+1 1 j—1 j—1 j j + -т +-n + 7 — 7 ui+1 ui — - ui ui+1 )dui+1 л dui

\2(фЩ))2 2(y(Ax))2 6 4 ^ ; 4 ; ;+7 ;

i 111 j+i j+i 1 j-i j-A j-i j

--1--+ - - —u, U i - -u. i U )aU: A au

-o n o T T T U U 1 T U i u

2(0(At))2 2(^(Ax))2 6 4 1 1-1 4 1-1 1

1 1 11 j j+A , j+1 , j

-2 +12 - 4 uj+1uj jduj+1 A duj

4( <KAt))2 1 9 4( <KAx))2

4(^(At))2 1 0 4( <KAx))2

4(^(At))2 1 0 4( <KAx))2

1 1 1- -r

4(0(Af))2 4(^(Ax))

1 1 j j+1 j+1 j + 12 - 4 U-1u. \aui-1 A du

+ ï^ - I uj 1uj-1 )auj-]1 A duj

1 1 j-1 j+1 j-1 j 2 + 12 - 4U Uj+1 ' j+1 A j =

(3.22)

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

uj+_2uj + ^U 1

^-p 2 p = Utt(xp,fj) + g((Af)2), p = i - 1, i, i + 1,

($(At))2 P j V 1

uh - 2uqq + uqq-1 (f(Ax))2 1

Uxx(xi,fq) + o((Ax)2), q = j - 1, j, j + 1,

1 j j+1 j+1 j j+1 j+1 j-1 j-1 j j-1 j-1 j 4 uj + uj-1uj uj-1 + UU U uh + U uj+1 <1)

= u3(xi,fj) + (Ax)2 (2^xi,fj)ux^xi,fj) + u2x{xi,tj)) (3.23)

+ (At)2 (2u(xi,tj)utt(xi,tj) + u2(xi,fj))

= u3(xi,tj) + o( (Ax)2 + (At)2),

1 j+1 j+1 j-1 j-1 j+1 j-1 j j — (j + j + j + j + 2uj + 2uj + 2uj+1 + 2uj-1 j

= u(xi,tj) + o( (Ax)2 + (At)2).

Combing these equations, we can readily observe that the nonstandard finite difference variational integrator (3.21) has truncation error O((Ax)2 + (At)2). 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.

0.4 0.2 0 ■0.2 60

■20 -10

Figure 1: The waveforms of linear wave equation (2.1) by the nonstandard finite difference variational integrator (2.14) (Ai = 0.1, Ax = 0.2).

Theorem 3.2. The nonstandard finite difference variational integrator (3.21) for the nonlinear KleinGordon equation (3.1) is multisymplectic, and its truncation error is O((Ax)2 + (Ai)2). 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

The nonstandard finite difference variational integrator (2.14) is an explicit five points scheme. We choose the denominator functions $ and f in as $(At) = 2sin(At/2) and f(At) = 2sin(Ax/2).

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.

u(x,0) = sechx, -10 <x< 10, ut(x, 0) = 0, -10 < x < 10,

and the periodic boundary conditions

u(-10, t) = u(10, t), ux(-10,t) = ux (10, t).

4.2. Nonlinear Klein-Gordon Equation

We now consider the nonlinear Klein-Gordon equation (3.1) with the initial condition

u(x, 0) = a(i 1 + cos^)), ut(x, 0)= 0, (4.3)

and periodic boundary conditions.

We use nonstandard finite difference variational integrator (3.21) to simulate this problem with amplitude A = 5. The nonstandard finite difference variational integrator (2.14) is an implicit nine-points nonstandard finite difference scheme. The denominator functions p and f 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.

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:

ut (x, 0) = c

i(x, 0) = v^sech

1 - c2sech (x -ct) ytanh (x - ct)j,

i(-10,t) = V2sech

u(10, t) = V2sech(

(-10 - ct) ,

(10 - ct)

where c = 1.2. The exact solution of the problem is

i(x, t) = v^sech

(10 - ct)

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 a\ = b = a2 = b2 = 1 here.

-0.5 -1

(a) Waveforms at the beginning from t = 0 to t = 2

-5 -10 -15

(b) Waveforms from t = 5 to t = 8

Figure 2: Waveforms of nonlinear Klein-Gordon equation by the nonstandard finite difference variational integrator (3.21) (At = 0.01, Ax = 0.01).

The Z^-norm errors at t = 1, t = 1.5, and t = 2 are listed in Tables 1 and 2. The orders in the tables are calculated with the formula [35, 36]

ln(Error(Ax1)/Error(Ax2)) ,.

Order «-——-——--. (4.6)

ln(Ax1/Ax2)

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 At = Ax, 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.

Journal of Applied Mathematics Table 1: The error and convergence orders of (3.9) for Klein-Gordon problem with At = Ax.

Mesh size t = 1 t = 1.5 t = 2

Ax = 0.1

Error 6.3e - 3 8.9e - 3 1.07e - 2

Order — — —

Ax = 0.05

Error 1.6e - 3 2.3e - 3 2.8e - 3

Order 1.9773 1.9522 1.9341

Ax = 0.025

Error 4.1645e - 4 5.8258e - 4 7.0059e - 4

Order 1.9419 1.9811 1.9988

Ax = 0.0125

Error 1.0480e - 4 1.4622e - 4 1.7565e - 4

Order 1.9905 1.9943 1.9959

Table 2: The error and convergence orders of (3.21) for Klein-Gordon problem with At = Ax.

Mesh size t = 1 t = 1.5 t = 2

Ax = 0.1

Error 1.28e - 2 1.81e - 2 2.25e - 2

Order — — —

Ax = 0.05

Error 3.3e - 3 4.6e - 3 5.7e - 3

Order 1.9556 1.9763 1.9809

Ax = 0.025

Error 8.3445e - 4 1.2e - 3 1.4e - 3

Order 1.9836 1.9386 2.0255

Ax = 0.0125

Error 2.0990e - 4 2.9458e - 4 3.6187e - 4

Order 1.9911 2.0263 1.9984

5. Conclusion

In this paper, we have considered a linear wave equation and a nonlinear Kleingordon 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

R. Abraham and J. E. Marsden, Foundations of Mechanics, Benjamin/Cummings Publishing, Reading, Mass, USA, 2nd edition, 1978.

J.-B. Chen, "Variational integrators and the finite element method," Applied Mathematics and Computation, vol. 196, no. 2, pp. 941-958, 2008.

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.

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.

J. E. Marsden and M. West, "Discrete mechanics and variational integrators," Acta Numerica, vol. 10, pp. 357-514, 2001.

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. M. Leok and J. Zhang, "Discrete Hamiltonian variational integrators," IMA Journal of Numerical Analysis, vol. 31, no. 4, pp. 1497-1532, 2011.

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 Kahler 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. Gonzalez-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. Gonzalez-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 Schrodinger 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 Schrodinger 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 Schrodinger 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 Schrodinger 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 Schrodinger 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.

Copyright of Journal of Applied Mathematics is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.