10,731 research outputs found
NLSEmagic: Nonlinear Schr\"odinger Equation Multidimensional Matlab-based GPU-accelerated Integrators using Compact High-order Schemes
We present a simple to use, yet powerful code package called NLSEmagic to
numerically integrate the nonlinear Schr\"odinger equation in one, two, and
three dimensions. NLSEmagic is a high-order finite-difference code package
which utilizes graphic processing unit (GPU) parallel architectures. The codes
running on the GPU are many times faster than their serial counterparts, and
are much cheaper to run than on standard parallel clusters. The codes are
developed with usability and portability in mind, and therefore are written to
interface with MATLAB utilizing custom GPU-enabled C codes with the
MEX-compiler interface. The packages are freely distributed, including user
manuals and set-up files.Comment: 37 pages, 13 figure
A High-Order Numerical Method for the Nonlinear Helmholtz Equation in Multidimensional Layered Media
We present a novel computational methodology for solving the scalar nonlinear
Helmholtz equation (NLH) that governs the propagation of laser light in Kerr
dielectrics. The methodology addresses two well-known challenges in nonlinear
optics: Singular behavior of solutions when the scattering in the medium is
assumed predominantly forward (paraxial regime), and the presence of
discontinuities in the % linear and nonlinear optical properties of the medium.
Specifically, we consider a slab of nonlinear material which may be grated in
the direction of propagation and which is immersed in a linear medium as a
whole. The key components of the methodology are a semi-compact high-order
finite-difference scheme that maintains accuracy across the discontinuities and
enables sub-wavelength resolution on large domains at a tolerable cost, a
nonlocal two-way artificial boundary condition (ABC) that simultaneously
facilitates the reflectionless propagation of the outgoing waves and forward
propagation of the given incoming waves, and a nonlinear solver based on
Newton's method.
The proposed methodology combines and substantially extends the capabilities
of our previous techniques built for 1Dand for multi-D. It facilitates a direct
numerical study of nonparaxial propagation and goes well beyond the approaches
in the literature based on the "augmented" paraxial models. In particular, it
provides the first ever evidence that the singularity of the solution indeed
disappears in the scalar NLH model that includes the nonparaxial effects. It
also enables simulation of the wavelength-width spatial solitons, as well as of
the counter-propagating solitons.Comment: 40 pages, 10 figure
The Vanishing Moment Method for Fully Nonlinear Second Order Partial Differential Equations: Formulation, Theory, and Numerical Analysis
The vanishing moment method was introduced by the authors in [37] as a
reliable methodology for computing viscosity solutions of fully nonlinear
second order partial differential equations (PDEs), in particular, using
Galerkin-type numerical methods such as finite element methods, spectral
methods, and discontinuous Galerkin methods, a task which has not been
practicable in the past. The crux of the vanishing moment method is the simple
idea of approximating a fully nonlinear second order PDE by a family
(parametrized by a small parameter \vepsi) of quasilinear higher order (in
particular, fourth order) PDEs. The primary objectives of this book are to
present a detailed convergent analysis for the method in the radial symmetric
case and to carry out a comprehensive finite element numerical analysis for the
vanishing moment equations (i.e., the regularized fourth order PDEs). Abstract
methodological and convergence analysis frameworks of conforming finite element
methods and mixed finite element methods are first developed for fully
nonlinear second order PDEs in general settings. The abstract frameworks are
then applied to three prototypical nonlinear equations, namely, the
Monge-Amp\`ere equation, the equation of prescribed Gauss curvature, and the
infinity-Laplacian equation. Numerical experiments are also presented for each
problem to validate the theoretical error estimate results and to gauge the
efficiency of the proposed numerical methods and the vanishing moment
methodology.Comment: 141 pages, 16 figure
- …