In this paper, we begin with the nonlinear Schrödinger/Gross-Pitaevskii equation (NLSE/GPE) for modeling Bose-Einstein condensation (BEC) and nonlinear optics as well as other applications, and discuss their dynamical properties ranging from time reversible, time transverse invariant, mass and energy conservation, dispersion relation to soliton solutions. Then, we review and compare different numerical methods for solving the NLSE/GPE including finite difference time domain methods and time-splitting spectral method, and discuss different absorbing boundary conditions. In addition, these numerical methods are extended to the NLSE/GPE with damping terms and/or an angular momentum rotation term as well as coupled NLSEs/GPEs. Finally, applications to simulate a quantized vortex lattice dynamics in a rotating BEC are reported.

No abstract

Abstract. In this paper, we consider the nonlinear Schrödinger equation ut + i∆u − F (u) = 0 in two dimensions. We show, by an operator-theoretic proof, that the well-known Lie and Strang formulae (which are splitting methods) are approximations of the exact solution of order 1 and 2 in time.

No abstract

No abstract

This paper addresses the construction of nonlinear integro-differential artificial boundary conditions for one-dimensional nonlinear cubic Schrödinger equations. Several ways of designing such conditions are provided and a theoretical classification of their accuracy is given. Semi-discrete time schemes based on the method developed by Durán and Sanz-Serna [IMA J. Numer. Anal. 20 (2) (2000), pp. 235-261] are derived for these unusual boundary conditions. Stability results are stated and several numerical tests are performed to analyze the capacity of the proposed approach.

This paper deals with the modeling of the ionospheric plasma. Starting from the two-fluid Euler–Maxwell equations, we present two hierarchies of models. The MHD hierarchy deals with large plasma density situations while the dynamo hierarchy is adapted to lower density situations. Most of the models encompassed by the dynamo hierarchy are classical ones, but we shall give a unified presentation of them which brings a new insight into their interrelations. By contrast, the MHD hierarchy involves a new (at least to the authors) model, the massless-MHD model. This is a diffusion system for the density and magnetic field which could be of great practical interest. Both hierarchies terminate with the "classical" Striation model, which we shall investigate in detail.

Abstract. This paper adresses the construction and study of a Crank-Nicolson-type discretization of the two-dimensional linear Schrödinger equation in a bounded domain Ω with artificial boundary conditions set on the arbitrarily shaped boundary of Ω. These conditions present the features of being differential in space and nonlocal in time since their definition involves some time fractional operators. After having proved the well-posedness of the continuous truncated initial boundary value problem, a semi-discrete Crank-Nicolson-type scheme for the bounded problem is introduced and its stability is provided. Next, the full discretization is realized by way of a standard finite-element method to preserve the stability of the scheme. Some numerical simulations are given to illustrate the effectiveness and flexibility of the method.