10,731 research outputs found

    NLSEmagic: Nonlinear Schr\"odinger Equation Multidimensional Matlab-based GPU-accelerated Integrators using Compact High-order Schemes

    Full text link
    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

    Full text link
    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

    Full text link
    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
    • …
    corecore