5,325 research outputs found
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
High-order numerical method for the nonlinear Helmholtz equation with material discontinuities in one space dimension
The nonlinear Helmholtz equation (NLH) models the propagation of
electromagnetic waves in Kerr media, and describes a range of important
phenomena in nonlinear optics and in other areas. In our previous work, we
developed a fourth order method for its numerical solution that involved an
iterative solver based on freezing the nonlinearity. The method enabled a
direct simulation of nonlinear self-focusing in the nonparaxial regime, and a
quantitative prediction of backscattering. However, our simulations showed that
there is a threshold value for the magnitude of the nonlinearity, above which
the iterations diverge. In this study, we numerically solve the one-dimensional
NLH using a Newton-type nonlinear solver. Because the Kerr nonlinearity
contains absolute values of the field, the NLH has to be recast as a system of
two real equations in order to apply Newton's method. Our numerical simulations
show that Newton's method converges rapidly and, in contradistinction with the
iterations based on freezing the nonlinearity, enables computations for very
high levels of nonlinearity. In addition, we introduce a novel compact
finite-volume fourth order discretization for the NLH with material
discontinuities.The one-dimensional results of the current paper create a
foundation for the analysis of multi-dimensional problems in the future.Comment: 47 pages, 8 figure
An efficient high-order algorithm for acoustic scattering from penetrable thin structures in three dimensions
This paper presents a high-order accelerated algorithm for the solution of the integral-equation formulation of volumetric scattering problems. The scheme is particularly well suited to the analysis of âthinâ structures as they arise in certain applications (e.g., material coatings); in addition, it is also designed to be used in conjunction with existing low-order FFT-based codes to upgrade their order of accuracy through a suitable treatment of material interfaces. The high-order convergence of the new procedure is attained through a combination of changes of parametric variables (to resolve the singularities of the Green function) and âpartitions of unityâ (to allow for a simple implementation of spectrally accurate quadratures away from singular points). Accelerated evaluations of the interaction between degrees of freedom, on the other hand, are accomplished by incorporating (two-face) equivalent source approximations on Cartesian grids. A detailed account of the main algorithmic components of the scheme are presented, together with a brief review of the corresponding error and performance analyses which are exemplified with a variety of numerical results
A hybridizable discontinuous Galerkin method for electromagnetics with a view on subsurface applications
Two Hybridizable Discontinuous Galerkin (HDG) schemes for the solution of
Maxwell's equations in the time domain are presented. The first method is based
on an electromagnetic diffusion equation, while the second is based on
Faraday's and Maxwell--Amp\`ere's laws. Both formulations include the diffusive
term depending on the conductivity of the medium. The three-dimensional
formulation of the electromagnetic diffusion equation in the framework of HDG
methods, the introduction of the conduction current term and the choice of the
electric field as hybrid variable in a mixed formulation are the key points of
the current study. Numerical results are provided for validation purposes and
convergence studies of spatial and temporal discretizations are carried out.
The test cases include both simulation in dielectric and conductive media
Benchmark of FEM, Waveguide and FDTD Algorithms for Rigorous Mask Simulation
An extremely fast time-harmonic finite element solver developed for the
transmission analysis of photonic crystals was applied to mask simulation
problems. The applicability was proven by examining a set of typical problems
and by a benchmarking against two established methods (FDTD and a differential
method) and an analytical example. The new finite element approach was up to
100 times faster than the competing approaches for moderate target accuracies,
and it was the only method which allowed to reach high target accuracies.Comment: 12 pages, 8 figures (see original publication for images with a
better resolution
- âŠ