Numerical Solutions of Coupled Klein-Gordon - Semantic Scholar

2 downloads 0 Views 451KB Size Report
The nonlinear coupled Klein-Gordon-Schr dinger equations describes a system of .... [16] G.L. Payne, D.R. Nicholson, R.M. Downie, Numerical solution of the ...
Applied Mathematical Sciences, Vol. 6, 2012, no. 45, 2245 - 2254

Numerical Solutions of Coupled Klein-GordonSchr dinger Equations by Finite Element Method Bouthina S. Ahmed Department of Mathematics, Collage of Girls Ain Shams University, Cairo, Egypt [email protected] Abstract The nonlinear coupled Klein-Gordon-Schr dinger equations describes a system of conserved scalar nucleons interacting with neutral scalar Mesons coupled with Yukawa interaction method . In this paper we derive a finite element scheme to solve these equations, we test this method for stability and accuracy, many numerical tests have been given to show the validity of the scheme. Keywords: Klein-Gordon-Schr dinger equations, finite element method.

1 Introduction During the past decades a wide range of physical phenomena is explained by dynamics of nonlinear waves. One of the most challenging and modern applications of the control of partial differential equations is the control of quantum mechanical system [1, 4, 5, 17-19]. In the present paper we study the standard system of nonlinear coupled Klein-GordonSchr dinger equations:

+ -

+ +

(1) (2)

where is a complex unknown function represent the slowly varying envelope of the highly oscillatory electric field, and ( ) a real scalar function denotes the fluctuation of the ion density about its equilibrium value. We supplement (1) by prescribing the initialboundary value conditions:

2246

B. S. Ahmed

(3) (4) and satisfy the following three conservation law (i) The total energy conservation law (5) (ii) The total momentum conservation law đťś” where Im stands for the imaginary part, and (iii) The charge conservation law

(6) denotes the complex conjugate of

(7) Eq. (1) consists of the classical Schrödinger equation and Klein Gordon equation. Over the last tens of years, there has been a great deal of work on the numerical solving Schrödinger equations. Recently, the multi-simplistic methods have been proposed and investigated for some important Hamiltonian partial differential equations, such as Klein Gordon Schrödinger eq. (1). Eq. (1) is a classical model described the interaction between conservative complex neutron field and neutral meson Yukawa in quantum field theory. Global strong solution in [9, 11] and stability of stationary state in [14] is proved for the Klein Gordon Schrödinger equations. Asymptotic behavior of the solutions is reported [10, 15]. The solitary wave solution of the Klein-Gordon- Schrödinger equations are obtained in [20]. Numerical studies for the KGS equations are few. We constructed an implicit scheme in [21]. However many finite difference schemes have been presented for the Zakharov equations presented [2, 4, 6, 7, 16, 21, 22]. In this paper we derive the finite elements method in section 2. The accuracy and stability of the scheme are given in section 3 and section 4. Numerical tests for single and interaction of two solitons are given in section 5. Conclusions are contained in Section 6.

2 Numerical method For the numerical study we decompose the complex function u in (1) by , also we set

Substituting in (1) we get (8) 2here

Numerical solutions

2247

The finite element method [8, 12, 13] for the resulting system (8) can be described as follows The weak form of (8) is obtained by taking –inner product with elements of the space . After integration by parts, and using the boundary conditions we obtain

where and

–inner product ( . , . ) is interpreted pointwise and denoted by

The interval

is discretized by uniform (N+1) grid points

where the grid spacing h is given by . Let be an N dimensional subspace of . We choose the elements of to be piecewise linear hat functions and denoted by . We assume the approximation of the exact solution is given by (10) where are time dependent coefficients which are to be determined. Then the discrete form of (9) is given by

(11) where , In order to deal with nonlinear term in (11), we use the product approximation [8] by which (12)

Since the trail and the test function is chosen to be piecewise linear functions

2248

B. S. Ahmed

–

We note that the solution in (11) is zero outside the interval and according to this we can write it as

and the nonlinear term as

By using these definitions, the inner product in (11) can be calculated as follows

= and

= By making use of these integrals in (11), this will lead us to the first order ordinary differential system (13) where

Numerical solutions

2249

2.1 Full discertization In order to solve the system (13), we first discretize the time interval by the grid points where is the time step size. We assume that approximation of the exact solution . Making use of the substitution:

is the

and Where

which is equivalent to the implicit midpoint rule? This will lead us to the implicit nonlinear scheme: (14)

3 Accuracy of the scheme To study the accuracy of the scheme (14) we replace the numerical solution exact solution in (14). Doing this the proposed scheme will be of the form

by the

(15) Taylor expansion of all terms in (15) can be displayed as follows

+

+

+

+

Substituting by this expansion in (15) we get

+

+

2250

B. S. Ahmed

By using the differential eq. (8) all terms inside the brackets are equal zero, then we have

Thus, tends to zero a h,k tends to zero, we deduce that theproposed scheme is a second order in time and it is consistent since the local truncation error tends to zero. 

4 Stability To study the stability of the proposed scheme, the von Neumann stability analysis will be used. This method can only be applied for linear scheme. By freezing all terms which make the scheme nonlinear [12], then eq. (8) has the form

(16) where

The finite difference scheme of eq. (16) is

(17) where

. Assume that (18)

be the test function đť›˝ â„ť â„ť and â„ť be the amplification matrix. The Von Neumann necessary condition for stability of the scheme is

Substituting in (17) we get after manipulation the amplification matrix is given by đťś” where đťś”

đťś”

đťś”

đťś”

đťś”

(19)

Numerical solutions

2251

The eigenvalues of the matrix

are

It is clear that all the modulus of the eigenvalues equal 1, then the scheme is unconditionally stable. the scheme is also consistent, and then according to Lax theorem the scheme is convergent.

5 Numerical results We study the accuracy of the proposed scheme by calculating the infinity norm of the error and we experiment the conservations laws. The trapezoidal rule is used to calculate the conservation quantities and the accuracy is measured by using error norm of u and v which given by    Two cases are considered single soliton and two soliton interactions. We take the parameters and the velocity  Single soliton solution In this test we take

In Table 1 and Table 2, we calculate the errors and the conserved quantities. It is very clear that the method produced accurate results and keeps the conserved quantities almost constants for the charge, momentum and nearly constants for the energy. In Figure 1 shows the evolution of single soliton and v at moving to the right.

Figure 1: Single soliton solution left

right v

2252

B. S. Ahmed Time 10 20 30 40 50 60

0.00538971 0.00943423 0.01168950 0.01493350 0.01904940 0.02281280

0.0539015 0.0941379 0.0116418 0.0148990 0.0189605 0.0227392

Table 1: Single soliton

Two Solitons Interaction In this test we take the interaction of two solitons, where the initial conditions are assumed of the form

Time 10 20 30 40 50 60

Energy Charge -1.184235 1.73641 -1.184230 1.73641 -1.184235 1.73641 -1.184238 1.73641 -1.184238 1.73641 -1.184230 1.73641 Table 2: Single soliton (conserved quantities)

Momentum 0.331929 0.331929 0.331929 0.331929 0.331929 0.331929

which represents the sum of two single solitons located at and respectively. Figure 2 display the interactions of two solitons and at with the parameters and .We noticed that the two waves approaches each other interact and leave the interaction unchanged in shape and velocity.

Figure 2: Interaction of two soliton solutions left

right v

Numerical solutions

2253

6 Conclusions The Klein-Gordon-Schrödinger equations have important role in quantum physics. In this paper, we developed a finite element method for solving the coupled nonlinear Schrödinger equations. The resulting scheme is second order in both directions space and time, and it is unconditionally stable. The numerical experiments are presented to demonstrate the effectiveness of the scheme and keep the conservative laws. The method simulates the interaction picture clearly in the case of two solitons interaction.

References [1] C. Bardeen, Feedback quantum control of molecular electronic population transfer, Chem. Phys. Lett. 280 (1997) 151.158. [2] Q. Chang, H. Jiang, A conservative scheme for the Zakharov equations, J. Comput. Phys. 113 (1994) 309-319. [3] Q. Chang, L. Zhang, Convergence and stability of a conservative finite difference scheme for a class of equation system in interaction of complex Schrödinger field and real Klein-Gordon field, Num. Math. J. Chinese Universities 22 (3) (2000) 362-370. [4] M. Dehghan, A. Mohebbi, Z. Asghari, Fourth-order compact solution of the nonlinear KleinGordon equation, Numer. Algor. 52 (2009) 523.540. [5] M. Dehghan, A. Shokri, A numerical method for solution of the two dimensional sine Gordon equation using the radial basis functions, Math. Comput. Simul. 79 (2008) 700.715. [6] R.T. Glassey, Convergence of an energy-preserving scheme for the Zakharov equations in one space dimension, Math. Comput. 58 (1992) 83-103. [7] R.T. Glassey, Approximate solutions to the Zakharov equations via finite difference, J. Comput. Phys. 100 (1992) 377-383. [8] D.F. Gri¢ ths, A.R. Mitchell, J. Li Morris, A numerical study of the nonlinear Schrödinger equation, NA/52, University of Dundee, 1982. [9] B. Guo, global solution for some problem of class of equations in interaction of complex Schrödinger field and real Klien-Goddon field, Science in China, series A 2 (1982) 97-107. [10] B. Guto, C. Miao, Global existence and asymptotic behavior of the solution for the Klein Gordon Schrödinger, Sinice in China, series A 25(7) (1995) 705-714. [11] N. Hayashi, W. Von Wahl, On the global strong solution of coupled Klein-Gordon-Schrödinger equations, J. Marh. Soc. Jpn. 39 (1987) 489-497. [12] M.S. Ismail, T. Taha, A finite element solution for the coupled Schrödinger equation, in: Proceedings of the 16th IMACS World Congress on Scientific Computation, Lausanne, 2000 [13] M.S. Ismail, Numerical solution of coupled nonlinear Schrödinger equation by Galerkin method, Math and Comp in Simulation 78 (2008) 532.547. [14] M. Ohto, Stability of stationary states for the coupled Klein-Gordon-Schrödinger equations, Non. Anal. 27 (4) (1996) 455-461. [15] T. Ozawa, Y. Tsutsumi, Asymptotic behavior of solutions for the coupled Klein-GordonSchrödinger equations, Adv. Stud. Pure Math. 23(1994) 295-305. [16] G.L. Payne, D.R. Nicholson, R.M. Downie, Numerical solution of the Zakharov equation, J. Comput. Phys. 50 (1983) 482-498. [17] A.P. Peirce, M.A. Dahleh, H. Rabitz, Optimal control of quantum mechanical systems: existence, numerical approximation, and applications, Phys. Rev.A 37 (1988), 4950.4964. [18] F. Shakeri, M. Dehghan, Numerical solution of the Klein-Gordon equation via Hes variational iteration method, Nonlinear Dynam. 51 (2008) 89.97. [19] S. Weinberg, The Quantum Theory of Fields, Foundations, vol. 1, Cambridge University Press, 1995.

2254

B. S. Ahmed

[20] H. Zhang, Solitary wave solutions for the coupled Klein-Gordon-Schrödinger equations, Acta Mathematica Scientia. 22A (3) (2002) 332-335. [21] L. Zhang, Q. Chang, Convergence and stability of a conservative finite difference scheme for a class of equation system in interaction of complex Schrödinger field and real Klein-Gordon field, Num. Math. J. Chinese Universities 22 (3) (2000) 362-370. [22] R.T. Glassey, Convergence of an energy-preserving scheme for the Zakharov equations in one space dimension, Math. Comput. 58 (1992) 83-103.