# Symplectic and Multi-Symplectic Methods For Coupled Nonlinear Schrödinger Equations With Periodic Solutions | Atılım University Open Archive System

Bu kaydın yasal hükümlere uygun olmadığını düşünüyorsanız lütfen sayfa sonundaki Hata Bildir bağlantısını takip ederek bildirimde bulununuz. Kayıtlar ilgili üniversite yöneticileri tarafından eklenmektedir. Nadiren de olsa kayıtlarla ilgili hatalar oluşabilmektedir. MİTOS internet üzerindeki herhangi bir ödev sitesi değildir!

**
Symplectic and Multi-Symplectic Methods For Coupled Nonlinear Schrödinger Equations With Periodic Solutions**

**BROWSE_DETAIL_CREATION_DATE:** 07-09-2015

### BROWSE_DETAIL_IDENTIFIER_SECTION

**BROWSE_DETAIL_TYPE:**
Article

**BROWSE_DETAIL_PUBLISH_STATE:**
Published

**BROWSE_DETAIL_FORMAT:**
No File

**BROWSE_DETAIL_LANG:**
English

**BROWSE_DETAIL_SUBJECTS:**
SCIENCE, Mathematics,

**BROWSE_DETAIL_CREATORS:**
Aydin , Ayhan (Author), Karasözen, Bülent (Author),

**BROWSE_DETAIL_CONTRIBUTERS:**

**BROWSE_DETAIL_DOI:**
doi:10.1016/j.cpc.2007.05.010

**BROWSE_DETAIL_URL:**
http://www.sciencedirect.com/science/article/pii/S0010465507002500

**BROWSE_DETAIL_PUBLICATION_IDENTIFIERS:** Kaynağa gitmek için URL' ye tıklayınız.

**BROWSE_DETAIL_TAB_FILE:**

BROWSE_DETAIL_NO_FILE

**BROWSE_DETAIL_TAB_KEYWORDS**

**BROWSE_DETAIL_TAB_ABSTRACT**

We consider for the integration of coupled nonlinear Schrödinger equations with periodic plane wave solutions a splitting method from the class of symplectic integrators and the multi-symplectic six-point scheme which is equivalent to the Preissman scheme. The numerical experiments show that both methods preserve very well the mass, energy and momentum in long-time evolution. The local errors in the energy are computed according to the discretizations in time and space for both methods. Due to its local nature, the multi-symplectic six-point scheme preserves the local invariants more accurately than the symplectic splitting method, but the global errors for conservation laws are almost the same.

**BROWSE_DETAIL_TAB_TOC**

## 1. Introduction

The nonlinear Schrödinger equation (NLS) arises as model equation with second-order dispersion and cubic nonlinearity for describing the dynamics of slowly varying wave packets in nonlinear optics and fluid dynamics. If there are two or more modes the coupled nonlinear Schrödinger (CNLS) system would be the relevant model. The two coupled nonlinear Schrödinger (CNLS) equations are given by

where ψ_{1}(x,t) and ψ

_{2}(x,t) are complex amplitudes or ‘envelopes’ of two wave packets,

*i*is the imaginary number,

*x*and

*t*are the space and time variables, respectively. The CNLS system has many applications including nonlinear optics[1] and [2] and geophysical fluid dynamics [3] and [4]. The parameters α

_{j}are the dispersion coefficients, σ

_{j}the Landau constants which describe the self-modulation of the wave packets, and v

_{12}and v

_{21}are the wave–wave interaction coefficients which describe the cross-modulations of the wave packets [5] and [6].

Analytical solutions can be obtained only for a few special integrable cases, like the Manakov model [1] where the self-modulation and wave–wave interaction coefficients are equal, i.e. α_{1}=α_{2}=1, σ_{1}=σ_{2}=v_{12}=v_{21}. Other integrable cases are: α_{1}=α_{2},σ_{1}=σ_{2}=v_{12}=v_{21} and α_{1}=−α_{2}, σ_{1}=σ_{2}=−v_{12}=−v_{21} which can be solved by the inverse scattering method [7]. For nonintegrable cases, where the parameters are different, numerical methods have to be used in order to understand different nonlinear phenomena that arise by the interaction of stable and unstable wave packets in the CNLS system.

There has been a lot work done for the solitary wave solutions of (1). Symplectic and multi-symplectic methods are also used for the numerical solution of soliton collision[8] and [9]. Little attention has been paid to plane wave solutions under periodic boundary conditions. In this work we will consider the plane wave solutions of (1) under periodic boundary conditions with period *L*

_{10}and ψ

_{20}are assumed to be real [5] and [6]. We takev

_{12}=v

_{21}=v because this choice of the parameters allows a multi-symplectic formulation of the CNLS system.

Besides the Manakov's case which corresponds to elliptical polarization we consider two other polarizations [5] and [6]:

- •
linear polarization: v

_{12}=v_{21}=2/3,- •
circular polarization: v

_{12}=v_{21}=2

_{1}=α

_{2}=1, σ

_{1}=σ

_{2}=1.

Several numerical studies have been carried out in recent years in order to understand the behavior of the periodic plane wave solutions of the CNLS system. In [5] a pseudospectral discretization was used in the space variables and the resulting system of ordinary differential equations (ODEs) were integrated using the standard fourth-order Runge–Kutta method. Because the method used in [5] is not conservative, no indication is given about the preservation of the conserved quantities of the CNLS system. In [6] the CNLS system with periodic boundary conditions was discretized in space using second-order finite differences and the resulting system of ordinary differential equations (ODEs) are integrated in time using the so-called Hopscotch method, a mixture of an explicit and an implicit method. The accuracy of the solutions were checked only using the norm conservation of each amplitude of the CNLS system. The choice of particular initial conditions may affect the preservation of some invariants. As mentioned in [10] in the context of multi-symplectic integrators usually symmetric initial conditions are used in the literature. Unsymmetric initial conditions may lead to an increase of errors in preserving some invariants like momentum over time [10].

The numerical solution of nonlinear wave equations using symplectic and multi-symplectic geometric integrators has been the subject of many studies in recent years (see for a recent review [11]). The NLS and CNLS systems represent an infinite-dimensional Hamiltonian system. After semi-discretization in the space variables one usually obtains a system of Hamiltonian ordinary differential equations, for which various symplectic integrators can be applied [8], [12] and [13]. In Section 2 we give a symplectic integrator based on the splitting of the semi-discretized Hamiltonian system in linear and nonlinear parts. Another way is to use the multi-symplectic structure of the NLS and CNLS equations and apply multi-symplectic integrators [9], [14] and [15]. A brief review of multi-symplectic integrators and their application to the CNLS system given in Section3. In Section 4 numerical results obtained with the symplectic splitting method and multi-symplectic six-point Preissman scheme are compared with results found in the literature for symmetric and unsymmetric initial conditions. Section 5 is devoted to concluding remarks.

**BROWSE_DETAIL_TAB_DESCRIPTION**

**BROWSE_DETAIL_TAB_RIGHTS**

**BROWSE_DETAIL_TAB_NOTES**

**BROWSE_DETAIL_TAB_REFERENCES**

**BROWSE_DETAIL_TAB_REFERENCED_BYS**