Numerical solution of the nonlinear Helmholtz equation using nonorthogonal expansions
Abstract
In [J. Comput. Phys. 171 (2001) 632677] we developed a fourthorder numerical method for solving the nonlinear Helmholtz equation which governs the propagation of timeharmonic laser beams in media with a Kerrtype nonlinearity. A key element of the algorithm was a new nonlocal twoway artificial boundary condition (ABC), set in the direction of beam propagation. This twoway ABC provided for reflectionless propagation of the outgoing waves while also fully transmitting the given incoming beam at the boundaries of the computational domain. Altogether, the algorithm of [J. Comput. Phys. 171 (2001) 632677] has allowed for a direct simulation of nonlinear selffocusing without neglecting nonparaxial effects and backscattering. To the best of our knowledge, this capacity has never been achieved previously in nonlinear optics. In the current paper, we propose an improved version of the algorithm. The principal innovation is that instead of using the Dirichlet boundary conditions in the direction orthogonal to that of the laser beam propagation, we now introduce Sommerfeldtype local radiation boundary conditions, which are constructed directly in the discrete framework. Numerically, implementation of the Sommerfeld conditions requires evaluation of eigenvalues and eigenvectors for a nonHermitian matrix. Subsequently, the separation of variables, which is a key building blockmore »
 Authors:
 Department of Applied Mathematics, School of Mathematical Sciences, Tel Aviv University, Ramat Aviv, Tel Aviv 69978 (Israel). Email: fibich@math.tau.ac.il
 Department of Mathematics and Center for Research in Scientific Computation (CRSC), North Carolina State University, Box 8205, Raleigh, NC 27695 (United States). Email: tsynkov@math.ncsu.edu
 Publication Date:
 OSTI Identifier:
 20687265
 Resource Type:
 Journal Article
 Resource Relation:
 Journal Name: Journal of Computational Physics; Journal Volume: 210; Journal Issue: 1; Other Information: DOI: 10.1016/j.jcp.2005.04.015; PII: S00219991(05)002068; Copyright (c) 2005 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved; Country of input: International Atomic Energy Agency (IAEA)
 Country of Publication:
 United States
 Language:
 English
 Subject:
 71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; ALGORITHMS; BACKSCATTERING; BOUNDARY CONDITIONS; DIFFRACTION; EIGENVALUES; EIGENVECTORS; EQUATIONS; HELMHOLTZ THEOREM; HERMITIAN MATRIX; IMPLEMENTATION; ITERATIVE METHODS; LASERS; NONLINEAR OPTICS; NONLINEAR PROBLEMS; NUMERICAL SOLUTION; PERFORMANCE; SIMULATION; SOMMERFELDWATSON THEORY
Citation Formats
Fibich, G., and Tsynkov, S. Numerical solution of the nonlinear Helmholtz equation using nonorthogonal expansions. United States: N. p., 2005.
Web. doi:10.1016/j.jcp.2005.04.015.
Fibich, G., & Tsynkov, S. Numerical solution of the nonlinear Helmholtz equation using nonorthogonal expansions. United States. doi:10.1016/j.jcp.2005.04.015.
Fibich, G., and Tsynkov, S. Sun .
"Numerical solution of the nonlinear Helmholtz equation using nonorthogonal expansions". United States.
doi:10.1016/j.jcp.2005.04.015.
@article{osti_20687265,
title = {Numerical solution of the nonlinear Helmholtz equation using nonorthogonal expansions},
author = {Fibich, G. and Tsynkov, S.},
abstractNote = {In [J. Comput. Phys. 171 (2001) 632677] we developed a fourthorder numerical method for solving the nonlinear Helmholtz equation which governs the propagation of timeharmonic laser beams in media with a Kerrtype nonlinearity. A key element of the algorithm was a new nonlocal twoway artificial boundary condition (ABC), set in the direction of beam propagation. This twoway ABC provided for reflectionless propagation of the outgoing waves while also fully transmitting the given incoming beam at the boundaries of the computational domain. Altogether, the algorithm of [J. Comput. Phys. 171 (2001) 632677] has allowed for a direct simulation of nonlinear selffocusing without neglecting nonparaxial effects and backscattering. To the best of our knowledge, this capacity has never been achieved previously in nonlinear optics. In the current paper, we propose an improved version of the algorithm. The principal innovation is that instead of using the Dirichlet boundary conditions in the direction orthogonal to that of the laser beam propagation, we now introduce Sommerfeldtype local radiation boundary conditions, which are constructed directly in the discrete framework. Numerically, implementation of the Sommerfeld conditions requires evaluation of eigenvalues and eigenvectors for a nonHermitian matrix. Subsequently, the separation of variables, which is a key building block of the aforementioned nonlocal ABC, is implemented through an expansion with respect to the nonorthogonal basis of the eigenvectors. Numerical simulations show that the new algorithm offers a considerable improvement in its numerical performance, as well as in the range of physical phenomena that it is capable of simulating.},
doi = {10.1016/j.jcp.2005.04.015},
journal = {Journal of Computational Physics},
number = 1,
volume = 210,
place = {United States},
year = {Sun Nov 20 00:00:00 EST 2005},
month = {Sun Nov 20 00:00:00 EST 2005}
}

Numerical simulation of reflecting structures by solution of the twodimensional Helmholtz equation
A method is described for modeling twodimensional reflecting structures based on a solution of the scalar Helmholtz equation. The equation is solved by use of an alternatingdirectionimplicit iterative method together with a semioptimum sequence of acceleration parameters similar to those introduced decades ago for the solution of elliptic equations with positivedefinite operators. The resulting technique is efficient and simple to program, permits the simulation of complex structures with modest storage requirements, and is of very general applicability. 
Accurate numerical solution of the Helmholtz equation by iterative Lanczos reduction
The Lanczos recursion algorithm is used to determine forwardpropagating solutions for both the paraxial and Helmholtz wave equations for longitudinally invariant refractive indices. By eigenvalue analysis it is demonstrated that the method gives extremely accurate solutions to both equations. 
Numerical solution of the Helmholtz equation for twodimensional polygonal regions
The Helmholtz equation together with associated boundary conditions can be solved using a dipole distribution on the boundary of any region of interest. If the region has corners, the distribution satisfies a singular integral equation. In this paper numerical techniques for the solution of this equation which take into account the analytic properties of the solution are discussed. Although a complete error analysis has not been developed, some indicators of the errors generated are considered. The technique is illustrated by comparing the numerical solution of the eigenvalue problem associated with various twodimensional polygonal regions with exact solutions. Numerical solutions ofmore »