Abstract and Applied Analysis

Volume 2013 (2013), Article ID 497439, 9 pages

http://dx.doi.org/10.1155/2013/497439

## Parallel Methods and Higher Dimensional NLS Equations

^{1}Department of Mathematics, College of Science, King Abdulaziz University, P.O. Box 80203, Jeddah, Saudi Arabia^{2}Department of Computer Science, University of Georgia, Athens, GA 30602-7404, USA

Received 3 June 2013; Accepted 28 July 2013

Academic Editor: Juan Carlos Cortés López

Copyright © 2013 M. S. Ismail and T. R. Taha. 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

Alternating direction implicit (ADI) schemes are proposed for the solution of the two-dimensional coupled nonlinear Schrödinger equation. These schemes are of second- and fourth-order accuracy in space and second order in time. The resulting schemes in each ADI computation step correspond to a block tridiagonal system which can be solved by using one-dimensional block tridiagonal algorithm with a considerable saving in computational time. These schemes are very well suited for parallel implementation on a high performance system with many processors due to the nature of the computation that involves solving the same block tridiagonal systems with many right hand sides. Numerical experiments on one processor system are conducted to demonstrate the efficiency and accuracy of these schemes by comparing them with the analytic solutions. The results show that the proposed schemes give highly accurate results.

#### 1. Introduction

In this paper, we consider the coupled nonlinear Schrödinger equation with initial conditions and Dirichlet boundary conditions where is a rectangular domain in . We assume that is the boundary of is the time interval, and , , , and are given sufficiently smooth functions. The two functions are representing the amplitudes of the two circularly polarized waves. The values of vary over a wide range; that is, and correspond to kerr type electronic nonlinearity [1]. The physical significance of system (1) can be seen in the transverse effects in nonlinear optics. Since solitons interact like particles, the studies for the interactions of solitons have been of both experimental and theoretical interest. Many numerical methods have been developed for solving the coupled nonlinear Schrödinger [2–6]. Many published works for solving the two-dimensional nonlinear Schrödinger equation are given in [7–11]. In this work, we are going to derive an ADI method for solving the two-dimensional coupled nonlinear Schrödinger equation.

In this paper, we derive two ADI schemes for solving the two-dimensional coupled nonlinear Schrödinger system, one is of second order in space and time directions, and the other one is of fourth order in space and second order in time. Both methods are producing a block nonlinear tridiagonal system; a fixed point method has been developed to solve this system. The proposed schemes are unconditionally stable using Fourier stability analysis.

The ADI method [8–15], which replaces the solution of multidimensional problems by sequences of one-dimensional cases, only needs to solve tridiagonal linear system or block tridiagonal systems, and the resulting schemes are unconditionally stable and received much attention in recent years.

Following Biswas [16], we derive the soliton solution of the system (1) which can be written as where and , , , , and are arbitrary constants (see the appendix). The system has the conserved quantities [15]

To avoid the complex computation, we assume

The system (1) can be written as where

System (10) can be written in matrix vector form as where

The paper is organized as follows. In Section 2 we give two ADI schemes for solving the two-dimensional coupled nonlinear Schrödinger equation, and in Section 3, we present the von Neumann stability analysis for the proposed schemes. Numerical experiments for several problems are presented in Section 4. A parallel algorithm for the proposed ADI schemes is given in Section 5. Conclusions are given in Section 6.

#### 2. Numerical Method

To derive the numerical schemes for solving system (1), we consider the domain of interest , such that . The domain is divided by a uniform mesh in each direction such that where and are the space and time step sizes, respectively. We denote and to be the exact and the numerical solutions at the point , respectively.

##### 2.1. Second-Order ADI Method

To derive the first scheme, we approximate the space derivative using the central difference formulae

We apply (15) to system (12); this will lead us to the following first order differential system in time

Equation (16) is of second order in space.

By applying the Crank-Nicolson method for the temporal discretization, we get the following difference scheme with accuracy : where

This scheme can be written as which can be approximated by using the factored form

By Taylor’s expansion, the last term in (20) can be written as which is of the same order of accuracy of the truncation error, and we ignore it to obtain

By introducing a new intermediate vector , we propose a D’Yakonov [12, 17] ADI-like scheme for the coupled system which is a nonlinear scheme. An iterative algorithm of fixed point nature can be used to solve the system of the nonlinear equations (23)-(24). The fixed point that we propose can be given by where the superscript denotes the th iterate for solving the nonlinear system of equations for each time step. The block tridiagonal matrix equations of (25) can be solved by Crout’s method. The initial iterate is chosen as

The iteration continues until the condition is satisfied.

The boundary value of the intermediate variable can be extracted from (24) and is given by the following formulas:

##### 2.2. Fourth-Order ADI Method

Now, we want to derive a highly accurate fourth order ADI scheme; to do this, we approximate the space derivatives by the following formulas

By using these approximations together with Crank-Nicolson for the time direction, we get the numerical scheme which can be written as and this can be written in the factored form as

The difference between (31) and (32) is which is of the same order of accuracy of the truncation error, we ignore it, and then we operate on both sides of (32) by

The fourth-order D’Yakonov ADI-like scheme in this case can be displayed as

Now the systems in (35) and (36) are nonlinear. By the similar approach used in the previous scheme, we can derive a fixed point iterative formulas.

The boundary value of the intermediate variable in (35) can be given by the following formulas:

We can easily write a generalized version of D’Yakonov ADI like method for solving the two-dimensional coupled nonlinear Schrödinger system (1) as for arbitrary value of , and for , we recover the second- and fourth-order schemes, respectively.

#### 3. Stability Analysis

To study the stability of the proposed scheme, we consider the von Neumann stability analysis which can be only applied for the linear finite difference scheme, so we consider the linear version of the generalized scheme where is constant and .

Suppose that the numerical solution can be expressed by Fourier series, whose typical term is where is the amplification matrix at time level , and are the wave numbers in and directions. Substituting (40) into (39) will lead to the matrix equation where

The eigenvalues of the matrix are given by

It can be easily shown that the modulus of the maximum eigenvalue of the amplification matrix is less than one; hence, the scheme is unconditionally stable according to the von Neumann stability analysis.

#### 4. Numerical Results

In this section, the efficiency and accuracy of the proposed schemes will be tested by comparing with the exact solutions. We will measure the accuracy of the proposed schemes using the norm. We compute the conserved quantity by using the trapezoidal rule.

##### 4.1. Example 1

In this example, we choose the initial conditions from the exact solutions where at . The following parameters are used

The boundary conditions are extracted from the exact solution (44).

Tables 1, 2, and 3 show the errors (ER) and the conserved quantities (I) for and , respectively. We have noticed that the scheme with produced highly accurate results, and this is due to the fourth order accuracy in space and second-order accuracy in time. The other two methods using and are of second-order accuracy in space and time. Figure 1 shows a single soliton at . Figure 2 shows the numerical solution at .

##### 4.2. Example 2

In this example we will choose the initial condition [1] where which represent two solitons of different amplitudes. The parameters , , , and are complex parameters. In this test, we choose the parameters

In Figures 3 and 4, we display the numerical solution of at and , while in Figures 5 and 6 we display the numerical solution of at and .

##### 4.3. Example 3

To study the interaction of two solitons, many numerical tests have been conducted with different initial conditions, and, among these tests, we select the initial conditions of the form which represent two solitons moving in the opposite direction and centered at and , respectively. In this test we choose the parameters

The interaction scenario is given in Figures 7 and 8. In Figure 7 we show the two solitons with two different amplitudes at . In Figure 8 we display the interaction where the two solitons interact at , and, in Figure 9, we display the two solitons after the interaction at . In Table 4, we display the conserved quantities, and we see that the numerical method we proposed conserves the conserved quantities almost exactly. It is easy to see that the interaction regiem is an inelastic one. See [18].

##### 4.4. Example 4

The system under consideration generates a progressive plane wave solutions [7]: where

Our numerical experiments are conducted in the domain with

Initial and boundary conditions are extracted from the exact solution. In this example we choose , the fourth order ADI method. The numerical results are presented in Table 5. Figure 10 displays the solution at , and Figure 11 displays, the solution at .

Table 4 displays the accuracy of the scheme and preserves the conserved quantities.

#### 5. Parallel Algorithm for the Proposed ADI Schemes

It is to be noted that the implementation of the ADI schemes requires solving the same block-tridiagonal matrix with different right-hand sides. This can be done efficiently by using the fast parallel algorithm given in [19]. This algorithm [19] is a generalization of the parallel dichotomy algorithm for solving tridiagonal liner system of equations [20].

It has been shown that this dichotomy yields almost a linear speedup on a high performance system with many processors [19]. We expect the same speedup for our proposed methods. The parallel implementation of the proposed ADI method will be reported in our future work.

#### 6. Concluding Remarks

In this present work, a generalized alternating direction implicit methods for solving two-dimensional coupled nonlinear Schrödinger equation have been established, the methods are unconditionally stable. The methods produced schemes of second and fourth order accuracy in space and second order in time according to the selected value of . The numerical results have shown that the proposed schemes successfully combine accuracy and efficiency for the two-dimensional coupled nonlinear Schrödinger system. The implementation of the ADI schemes requires solving the same block-tridiagonal matrix with different right-hand sides. This can be done by using the fast parallel algorithm given in [19]. This algorithm yields almost a linear speedup on a high performance system with many processors [19]. We expect the same speedup for our proposed methods. The proposed schemes can be easily extended to solve higher dimensional coupled nonlinear Schrödinger system.

#### Appendix

The exact solution of the nonlinear system can be derived as follows: we assume [19] where

Now from , we deduce and similar expressions can be obtained for the function . By substituting these expressions into the given system, and equating the real and imaginary parts, this will produce the following relations: By equating the coefficient of in , we obtain and by equating the coefficients of , we will get the relation From the second equation of the proposed system , we will get the relation and from and , we deduce that , and then we get the relation with , , , and as arbitrary complex parameters.

#### Acknowledgments

This work was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under Grant no. (130-074-D1433). The authors, therefore, acknowledge with thanks the DSR technical and financial support.

#### References

- H.-Q. Zhang, X.-H. Meng, T. Xu, L.-L. Li, and B. Tian, “Interactions of bright solitons for the ($2+1$)-dimensional coupled nonlinear Schrödinger equations from optical fibres with symbolic computation,”
*Physica Scripta*, vol. 75, no. 4, pp. 537–542, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - M. S. Ismail and T. R. Taha, “Numerical simulation of coupled nonlinear Schrödinger equation,”
*Mathematics and Computers in Simulation B*, vol. 56, no. 6, pp. 547–562, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. S. Ismail, “A fourth-order explicit schemes for the coupled nonlinear Schrödinger equation,”
*Applied Mathematics and Computation*, vol. 196, no. 1, pp. 273–284, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. S. Ismail, “Numerical solution of coupled nonlinear Schrödinger equation by Galerkin method,”
*Mathematics and Computers in Simulation*, vol. 78, no. 4, pp. 532–547, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. S. Ismail and S. Z. Alamri, “Highly accurate finite difference method for coupled nonlinear Schrödinger equation,”
*International Journal of Computer Mathematics*, vol. 81, no. 3, pp. 333–351, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. R. Taha and X. Xu, “Parallel split-step fourier methods for the coupled nonlinear Schrödinger type equations,”
*The Journal of Supercomputing*, vol. 32, no. 1, pp. 5–23, 2005. View at Publisher · View at Google Scholar · View at Scopus - Z. Gao and S. Xie, “Fourth-order alternating direction implicit compact finite difference schemes for two-dimensional Schrödinger equations,”
*Applied Numerical Mathematics*, vol. 61, no. 4, pp. 593–614, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - L. Kong, Y. Duan, L. Wang, X. Yin, and Y. Ma, “Spectral-like resolution compact ADI finite difference method for the multi-dimensional Schrödinger equations,”
*Mathematical and Computer Modelling*, vol. 55, no. 5-6, pp. 1798–1812, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. F. Tian and P. X. Yu, “High-order compact ADI (HOC-ADI) method for solving unsteady 2D Schrödinger equation,”
*Computer Physics Communications*, vol. 181, no. 5, pp. 861–868, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. Wang, B. Guo, and Q. Xu, “Fourth order compact and energy conservative difference schemes for the nonlinear Schrödinger equation in two dimensions,”
*Journal of Computational Physics*, vol. 243, pp. 382–399, 2013. View at Google Scholar - Y. Xu and L. Zhang, “Alternating direction implicit method for solving two-dimensional cubic nonlinear Schrödinger equation,”
*Computer Physics Communications*, vol. 183, no. 5, pp. 1082–1093, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - S. Karaa and J. Zhang, “High order ADI method for solving unsteady convection-diffusion problems,”
*Journal of Computational Physics*, vol. 198, no. 1, pp. 1–9, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Sapagovas and K. Jakubeliene, “Alternating direction method for two-dimensional parabolic equation with nonlocal integral condition,”
*Nonlinear Analysis*, vol. 17, no. 1, pp. 91–98, 2012. View at Google Scholar · View at MathSciNet - M. Subasi, “On the finite difference schemes for the numerical solution of two dimensional Schrödinger equation,”
*Numerical Methods for Partial Differential Equations*, vol. 18, pp. 124–134, 2002. View at Google Scholar - D. You, “A high-order Padé ADI method for unsteady convection-diffusion equations,”
*Journal of Computational Physics*, vol. 214, no. 1, pp. 1–11, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Biswas, “1-soliton solution of $1+2$ dimensional nonlinear Schrödinger equations in power law media,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 14, pp. 1830–1833, 2009. View at Google Scholar - A. R. Mitchell and D. F. Griffiths,
*The Finite Difference Method in Partial Differential Equations*, John Wiley & Sons, Chichester, UK, 1980. View at MathSciNet - C. Besse, N. J. Mauser, and H. P. Stimming, “Numerical studies for nonlinear Schrodinger equations: the Schrodinger-Poisson-X
*α*model and Davey Stewartson systems,” In press. - A. V. Terekhov, “A fast parallel algorithm for solving block-tridiagonal systems of liner equations including the domain decomposition method,”
*Parallel Computing*, vol. 39, pp. 245–258, 2013. View at Google Scholar - A. V. Terekhov, “Parallel dichotomy algorithm for solving tridiagonal system of linear equations with multiple right-hand sides,”
*Parallel Computing*, vol. 36, no. 8, pp. 423–438, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet