Domain decomposition algorithms for the two dimensional nonlinear Schrödinger equation and simulation of Bose–Einstein condensates

Besse, Christophe; Xing, Feng

The SMAI Journal of computational mathematics, Tome 2 (2016), pp. 277-300.

In this paper, we apply the optimized Schwarz method to the two dimensional nonlinear Schrödinger equation and extend this method to the simulation of Bose–Einstein condensates (Gross– Pitaevskii equation). We propose an extended version of the Schwartz method by introducing a preconditioned algorithm. The two algorithms are studied numerically. The experiments show that the preconditioned algorithm improves the convergence rate and reduces the computation time. In addition, the classical Robin condition and a newly constructed absorbing condition are used as transmission conditions.

Classification : 35Q55, 65M55, 65Y05, 65M60
Mots clés : nonlinear Schrödinger equation, rotating Bose–Einstein condensate, optimized Schwarz method, preconditioned algorithm, parallel algorithm

[1] Antoine, X.; Bao, W.; Besse, C. Computational methods for the dynamics of the nonlinear Schrödinger/Gross–Pitaevskii equations, Comput. Phys. Commun., Volume 184 (2013) no. 12, pp. 2621-2633 | DOI | Zbl

[2] Antoine, X.; Besse, C.; Klein, P. Absorbing Boundary Conditions for the Two-Dimensional Schrödinger Equation With an Exterior Potential Part I: Construction and a Priori Estimates, Math. Model. Methods Appl. Sci., Volume 22 (2012) no. 10 | DOI | Zbl

[3] Antoine, X.; Besse, C.; Klein, P. Absorbing boundary conditions for the two-dimensional Schrödinger equation with an exterior potential. Part II: Discretization and numerical results, Numer. Math., Volume 125 (2013) no. 2, pp. 191-223 | DOI | Zbl

[4] Antoine, X.; Duboscq, R. Computer Physics cations GPELab , a Matlab Toolbox to solve Gross–Pitaevskii Equations I : computation of stationary solutions, Comput. Phys. Commun., Volume 185 (2014) no. 11, pp. 2969-2991 | DOI | Zbl

[5] Antoine, X.; Lorin, E.; Bandrauk, A. D. Domain Decomposition Method and High-Order Absorbing Boundary Conditions for the Numerical Simulation of the Time Dependent Schrödinger Equation with Ionization and Recombination by Intense Electric Field, J. Sci. Comput. (2014), pp. 1-27 | Zbl

[6] Balay, S.; Adams, M. F.; Brown, J.; Brune, P.; Buschelman, K.; Eijkhout, V.; Gropp, W. D.; Kaushik, D.; Knepley, M. G.; McInnes, L. C.; Rupp, K.; Smith, B. F.; Zhang, H. PETSc Users Manual (2013) no. ANL-95/11 - Revision 3.4 (Technical report)

[7] Bao, W.; Cai, Y. Mathematical theory and numerical methods for Bose-Einstein condensation, Kinet. Relat. Model., Volume 6 (2012) no. 1, pp. 1-135 | DOI | MR

[8] Bao, W.; Chern, I.-L.; Lim, F. Y. Efficient and spectrally accurate numerical methods for computing ground and first excited states in Bose–Einstein condensates, J. Comput. Phys., Volume 219 (2006) no. 2, pp. 836-854 | DOI | MR | Zbl

[9] Bao, W.; Du, Q. Computing the ground state solution of Bose–Einstein condensates by a normalized gradient flow, SIAM J. Sci. Comput., Volume 25 (2004) no. 5, pp. 1674-1697 | DOI | MR | Zbl

[10] Bao, W.; Marahrens, D.; Tang, Q.; Zhang, Y. A Simple and Efficient Numerical Method for Computing the Dynamics of Rotating Bose–Einstein Condensates via Rotating Lagrangian Coordinates, SIAM J. Sci. Comput., Volume 35 (2013) no. 6 | DOI | MR | Zbl

[11] Bao, W.; Markowich, P. A.; Wang, H. Ground, Symmetric and Central Vortex States in Rotating Bose-Einstein Condensates, Commun. Math. Sci., Volume 3 (2005) no. 1, pp. 57-88 | DOI | MR | Zbl

[12] Besse, C.; Xing, F. Domain decomposition algorithms for two dimensional linear Schrödinger equation (2015) (https://arxiv.org/abs/1506.05639)

[13] Besse, C.; Xing, F. Schwarz waveform relaxation method for one dimensional Schrödinger equation with general potential, Numerical Algorithms, accepted (2016) | DOI

[14] Boubendir, Y.; Antoine, X.; Geuzaine, C. A quasi-optimal non-overlapping domain decomposition algorithm for the Helmholtz equation, J. Comput. Phys., Volume 231 (2012) no. 2, pp. 262-280 | DOI | MR | Zbl

[15] Durán, A.; Sanz-Serna, J.-M. The numerical integration of relative equilibrium solutions. The nonlinear Schrodinger equation, IMA J. Numer. Anal., Volume 20 (2000) no. 2, pp. 235-261 | DOI | MR | Zbl

[16] Gander, M. J. Optimized Schwarz Methods, SIAM J. Numer. Anal., Volume 44 (2006) no. 2, pp. 699-731 | DOI | MR | Zbl

[17] Gander, M. J. Schwarz methods over the course of time, Electron. Trans. Numer. Anal., Volume 31 (2008), pp. 228-255 | MR | Zbl

[18] Gander, M. J.; Halpern, L. Optimized Schwarz Waveform Relaxation Methods for Advection Reaction Diffusion Problems, SIAM J. Numer. Anal., Volume 45 (2007) no. 2, pp. 666-697 | DOI | MR | Zbl

[19] Gander, M. J.; Halpern, L. Méthodes de décomposition de domaine, Encyclopédie électronique pour les ingénieurs, 2012, 81 pages

[20] Gander, M. J.; Magoules, F.; Nataf, F. Optimized Schwarz methods without overlap for the Helmholtz equation, SIAM J. Sci. Comput., Volume 24 (2002) no. 1, pp. 38-60 | DOI | MR | Zbl

[21] Halpern, L.; Szeftel, J. Optimized and quasi-optimal Schwarz waveform relaxation for the one dimensional Schrödinger equation, Math. Model. Methods Appl. Sci., Volume 20 (2010) no. 12, pp. 2167-2199 | DOI | Zbl

[22] Lions, P.-L. On the Schwarz alternating method. III: a variant for nonoverlapping subdomains, Third Int. Symp. domain Decompos. methods Partial Differ. equations, Volume 6 (1990), pp. 202-223 | Zbl

[23] Loisel, S. Condition Number Estimates for the Nonoverlapping Optimized Schwarz Method and the 2-Lagrange Multiplier Method for General Domains and Cross Points, SIAM J. Numer. Anal., Volume 51 (2013) no. 6, pp. 3062-3083 | DOI | MR | Zbl

[24] Message Passing Interface Forum MPI : A Message-Passing Interface Standard Version 3.0, 2012 (https://spcl.inf.ethz.ch/Publications/index.php?pub=160)

[25] Nataf, F.; Rogier, F. Factorization of the convection-diffusion operator and a (possibly) non overlapping Schwarz method, Contemp. Math. (1994) | DOI | Zbl