European Journal of Applied Mathematics

Papers

Numerical simulation of vortex dynamics in Ginzburg-Landau-Schrödinger equation

YANZHI ZHANGa1, WEIZHU BAOa2 and QIANG DUa3

a1 Department of Mathematics, National University of Singapore, Singapore 117543 email: zhyanzhi@gmail.com

a2 Department of Mathematics and Center for Computational Science and Engineering National University of Singapore, Singapore 117543 email: bao@math.nus.edu.sg, URL: http://www.math.nus.edu.sg/~bao/

a3 Department of Mathematics, Pennsylvania State University, University Park, PA 16802, USA email: qdu@math.psu.edu, URL: http://www.math.psu.edu/qdu

Abstract

The rich dynamics of quantized vortices governed by the Ginzburg-Landau-Schrödinger equation (GLSE) is an interesting problem studied in many application fields. Although recent mathematical analysis and numerical simulations have led to a much better understanding of such dynamics, many important questions remain open. In this article, we consider numerical simulations of the GLSE in two dimensions with non-zero far-field conditions. Using two-dimensional polar coordinates, transversely highly oscillating far-field conditions can be efficiently resolved in the phase space, thus giving rise to an unconditionally stable, efficient and accurate time-splitting method for the problem under consideration. This method is also time reversible for the case of the non-linear Schrödinger equation. By applying this numerical method to the GLSE, we obtain some conclusive experimental findings on issues such as the stability of quantized vortex, interaction of two vortices, dynamics of the quantized vortex lattice and the motion of vortex with an inhomogeneous external potential. Discussions on these simulation results and the recent theoretical studies are made to provide further understanding of the vortex stability and vortex dynamics described by the GLSE.

(Received May 08 2007)

(Revised July 23 2007)

(Online publication September 04 2007)