An official website of the United States government
In this work, we develop numerical methods to solve forward and inverse wave problems for a nonlinear Helmholtz equation defined in a spherical shell between two concentric spheres centred at the origin. A spectral method is developed to solve the forward problem while a combination of a finite difference approximation and the least squares method are derived for the inverse problem. Numerical examples are given to verify the method. ReferencesR. Askey. Orthogonal polynomials and special functions. CBMS-NSF Regional Conference Series in Applied Mathematics. SIAM, 1975. doi: 10.1137/1.9781611970470G. Baruch, G. Fibich, and S. Tsynkov. High-order numerical method for the nonlinear Helmholtz equation with material discontinuities in one space dimension. Nonlinear Photonics. Optica Publishing Group, 2007. doi: 10.1364/np.2007.ntha6G. Fibich and S. Tsynkov. High-Order Two-Way Artificial Boundary Conditions for Nonlinear Wave Propagation with Backscattering. J. Comput. Phys. 171 (2001), pp. 632–677. doi: 10.1006/jcph.2001.6800G. Fibich and S. Tsynkov. Numerical solution of the nonlinear Helmholtz equation using nonorthogonal expansions. J. Comput. Phys. 210 (2005), pp. 183–224. doi: 10.1016/j.jcp.2005.04.015P. M. Morse and K. U. Ingard. Theoretical Acoustics. International Series in Pure and Applied Physics. McGraw-Hill Book Company, 1968G. N. Watson. A treatise on the theory of Bessel functions. International Series in Pure and Applied Physics. Cambridge Mathematical Library, 1996. url: https://www.cambridge.org/au/universitypress/subjects/mathematics/real-and-complex-analysis/treatise-theory-bessel-functions-2nd-edition-1?format=PB&isbn=9780521483919
more »
« less
Exact wave solver for nonparaxial laser beam propagation
Simulations of inertial confinement fusion (ICF) experiments require high-fidelity models for laser beam propagation in a nonuniform plasma with varying index of refraction. We describe a new numerical wave solver that is applicable to centimeter-scale length plasmas encountered in indirect drive ICF applications. The one-way Helmholtz equation (OHE) generalizes the time-harmonic paraxial wave equation to large angles. Here, we present a methodology to numerically evaluate the exact solution to the OHE. This solution is computed by analytically advancing eigenfunctions of the one-way Helmholtz operator along a propagation direction and is applicable to any given index of a refraction profile. We compare our exact method with a commonly used approximate split-step technique for solving the OHE. As a test problem, we consider nonparaxial propagation of Gaussian and speckled beams in a plasma density channel with internal reflection. We find that the split-step approach incurs significant errors compared to the exact solution computed using the novel algorithm.
more »
« less
- Award ID(s):
- 2209472
- PAR ID:
- 10512242
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- Physics of Plasmas
- Volume:
- 31
- Issue:
- 5
- ISSN:
- 1070-664X
- Subject(s) / Keyword(s):
- Wave propagation, Operator theory, Partial differential equations, Numerical integration, Function approximation, Geometrical optics, Gaussian beam, Plasma confinement, Plasma properties and parameters
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This study explores the generation of electrostatic (ES) electron Kelvin–Helmholtz instability (EKHI) in collisionless plasma with a step-function electron velocity shear akin to that developed in the electron diffusion region in magnetic reconnection. In incompressible plasma, ES EKHI does not arise in any velocity shear profile due to the decoupling of the electric potential from the electron momentum equation. Instead, a fluid-like Kelvin–Helmholtz instability (KHI) can arise. However, in compressible plasma, the compressibility couples the electric potential with the electron dynamics, leading to the emergence of a new ES mode EKHI on Debye length λDe, accompanied by the co-generation of an electron acoustic-like wave. The minimum threshold of ES EKHI is ΔU>2cse, i.e., the electron velocity shear is larger than twice the electron acoustic speed cse. The corresponding growth rate is Im(ω)=((ΔU/cse)2−4)1/2ωpe, where ωpe is the electron plasma frequency.more » « less
-
This paper concerns the analysis of a passive, broadband approximate cloaking scheme for the Helmholtz equation in Rd for d = 2 or d = 3. Using ideas from transformation optics, we construct an approximate cloak by “blowing up” a small ball of radius ϵ > 0 to one of radius 1. In the anisotropic cloaking layer resulting from the “blow-up” change of variables, we incorporate a Drude-Lorentz- type model for the index of refraction, and we assume that the cloaked object is a soft (perfectly conducting) obstacle. We first show that (for any fixed ϵ) there are no real transmission eigenvalues associated with the inhomogeneity representing the cloak, which implies that the cloaking devices we have created will not yield perfect cloaking at any frequency, even for a single incident time harmonic wave. Secondly, we establish estimates on the scattered field due to an arbitrary time harmonic incident wave. These estimates show that, as ϵ approaches 0, the L2 -norm of the scattered field outside the cloak, and its far field pattern, approach 0 uniformly over any bounded band of frequencies. In other words: our scheme leads to broadband approximate cloaking for arbitrary incident time harmonic waves.more » « less
-
The YBJ equation (Young & Ben Jelloul, J. Marine Res. , vol. 55, 1997, pp. 735–766) provides a phase-averaged description of the propagation of near-inertial waves (NIWs) through a geostrophic flow. YBJ is obtained via an asymptotic expansion based on the limit $$\mathit{Bu}\rightarrow 0$$ , where $$\mathit{Bu}$$ is the Burger number of the NIWs. Here we develop an improved version, the YBJ + equation. In common with an earlier improvement proposed by Thomas, Smith & Bühler ( J. Fluid Mech. , vol. 817, 2017, pp. 406–438), YBJ + has a dispersion relation that is second-order accurate in $$\mathit{Bu}$$ . (YBJ is first-order accurate.) Thus both improvements have the same formal justification. But the dispersion relation of YBJ + is a Padé approximant to the exact dispersion relation and with $$\mathit{Bu}$$ of order unity this is significantly more accurate than the power-series approximation of Thomas et al. (2017). Moreover, in the limit of high horizontal wavenumber $$k\rightarrow \infty$$ , the wave frequency of YBJ + asymptotes to twice the inertial frequency $2f$ . This enables solution of YBJ + with explicit time-stepping schemes using a time step determined by stable integration of oscillations with frequency $2f$ . Other phase-averaged equations have dispersion relations with frequency increasing as $$k^{2}$$ (YBJ) or $$k^{4}$$ (Thomas et al. 2017): in these cases stable integration with an explicit scheme becomes impractical with increasing horizontal resolution. The YBJ + equation is tested by comparing its numerical solutions with those of the Boussinesq and YBJ equations. In virtually all cases, YBJ + is more accurate than YBJ. The error, however, does not go rapidly to zero as the Burger number characterizing the initial condition is reduced: advection and refraction by geostrophic eddies reduces in the initial length scale of NIWs so that $$\mathit{Bu}$$ increases with time. This increase, if unchecked, would destroy the approximation. We show, however, that dispersion limits the damage by confining most of the wave energy to low $$\mathit{Bu}$$ . In other words, advection and refraction by geostrophic flows does not result in a strong transfer of initially near-inertial energy out of the near-inertial frequency band.more » « less
-
We study the paraxial wave equation with a randomly perturbed index of refraction, which can model the propagation of a wave beam in a turbulent medium. The random perturbation is a stationary and isotropic process with a general form of the covariance that may or may not be integrable. We focus attention mostly on the nonintegrable case, which corresponds to a random perturbation with long-range correlations, that is, relevant for propagation through a cloudy turbulent atmosphere. The analysis is carried out in a high-frequency regime where the forward scattering approximation holds. It reveals that the randomization of the wave field is multiscale: The travel time of the wave front is randomized at short distances of propagation, and it can be described by a fractional Brownian motion. The wave field observed in the random travel time frame is affected by the random perturbations at long distances, and it is described by a Schr\"odinger-type equation driven by a standard Brownian field. We use these results to quantify how scattering leads to decorrelation of the spatial and spectral components of the wave field and to a deformation of the pulse emitted by the source. These are important questions for applications, such as imaging and free space communications with pulsed laser beams through a turbulent atmosphere. We also compare the results with those used in the optics literature, which are based on the Kolmogorov model of turbulence. We show explicitly that the commonly used approximations for the decorrelation of spatial and spectral components are appropriate for the Kolmogorov model but fail for models with long-range correlations.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0198523
