Table of Content

Open Access


Computation of Nonlinear Schrödinger Equation on an Open Waveguide Terminated by a PML

Jianxin Zhu1, Zheqi Shen1
Department of Mathematics, Zhejiang University, Hangzhou 310027, China, Email:

Computer Modeling in Engineering & Sciences 2011, 71(4), 347-362.


It is known that the perfectly matched layer (PML) is a powerful tool to truncate the unbounded domain. Recently, the PML technique has been introduced in the computation of nonlinear Schrödinger equations (NSE), in which the nonlinearity is separated by some efficient time-splitting methods. A major task in the study of PML is that the original equation is modified by a factor c which varies fast inside the layer. And a large number of grid points are needed to capture the profile of c in the discretization. In this paper, the possibility is discussed for using some nonuniform finite difference schemes in spatial discretization. It is proved that the uniform refinement inside the PML will cause spurious reflections at the interface, and therefore is invalid. As a remedy, a new method for the discretization is proposed by applying an additional coordinate transformation. This method effectively reduces the error caused by the discretization of PML, and it is quite useful in long time computation of time-dependent problems in open waveguides. These discussions are essential for the further research of the numerical solutions of NSE on unbounded domains.


nonlinear Schrödinger equation, open waveguide, PML, nonuniform discretization.

Cite This Article

Zhu, J., Shen, Z. (2011). Computation of Nonlinear Schrödinger Equation on an Open Waveguide Terminated by a PML. CMES-Computer Modeling in Engineering & Sciences, 71(4), 347–362.

This work is licensed under a Creative Commons Attribution 4.0 International License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • 841


  • 653


  • 0


Share Link

WeChat scan