The wave equation and the green's dyadic for bounded magnetoplasmas.
Separation method solution of dyadic-vector Helmholtz equation representing electromagnetic wave propagation in bounded magnetoplasmas
HELMHOLTZ EQUATION↗
SEARCH · Engineering Papers
Results for “Helmholtz equation”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
ARES v1.x - Performance Portable Tool to Simulate Supernovae based on Parthenon Framework
Historically, codes for simulating supernovae (such as Arepo, FLASH or LEAFS) have been at the forefront of scientific high-performance computing to the immense computational resources required for full 3D simulations. However, given the shift towards heterogenous HPC architectures, many current-generation codes are at the risk of losing their competitiveness as they are only designed to run on homogeneous CPU-only systems. There exist several efforts to enable these codes for GPU’s, however, these efforts only consider specific architectures or vendors (e.g., implement only CUDA or HIP), limiting themselves to a small range of exascale computing systems. Frameworks such as Kokkos aim to provide a framework which is agnostic of the targeted architecture, enabling the development of performant and portable code. In the Ares code, we develop a performance portable tool to simulate supernovae based on the Parthenon Framework, which in turn uses Kokkos in the background. Here, the Parthenon Framework provides an interface to the underlying mesh-refinement routines, which form the backbone of our code. In addition, we incorporate the already existing Singularity-EOS toolkit to provide us with various equations of state, primarily the Helmholtz equation of state. We also include the JINA Reaclib as a basis for our nuclear network solver. Finally, we implement a gravity solver to complete the required physics. This setup will provide us with a minimal code base to simulate supernova in a similar style to the tried-and-tested Arepo code, but in a futureproof performance portable framework.
Lim, Hyun↗
Analysis of sound propagation in ducts using the wave envelope concept
A finite difference formulation is presented for sound propagation in a rectangular two-dimensional duct without steady flow for plane wave input. Before the difference equations are formulated, the governing Helmholtz equation is first transformed to a form whose solution does not oscillate along the length of the duct. This transformation reduces the required number of grid points by an order of magnitude, and the number of grid points becomes independent of the sound frequency. Physically, the transformed pressure represents the amplitude of the conventional sound wave. Example solutions are presented for sound propagation in a one-dimensional straight hard-wall duct and in a two-dimensional straight soft-wall duct without steady flow. The numerical solutions show evidence of the existence along the duct wall of a developing acoustic pressure diffusion boundary layer which is similar in nature to the conventional viscous flow boundary layer. In order to better illustrate this concept, the wave equation and boundary conditions are written such that the frequency no longer appears explicitly in them. The frequency effects in duct propagation can be visualized solely as an expansion and stretching of the suppressor duct.
Baumeister, K. J.↗
Generalized wave envelope analysis of sound propagation in ducts with stepped noise source profiles and variable axial impedance
A finite difference formulation is presented for sound propagation in a rectangular two-dimensional duct without steady flow. Before the difference equations are formulated, the governing Helmholtz equation is first transformed to a form whose solution tends not to oscillate along the length of the duct. This transformation reduces the required number of grid points by an order of magnitude. Example solutions indicate that stepped noise source profiles have much higher attenuation than plane waves in a uniform impedance liner. Also, multiple stepped impedance liners are shown to have higher attenuation than uniform ducts if the impedances are chosen properly. For optimum noise reduction with axial variations in impedance, the numerical analysis indicates that for a plane wave input the resistance should be near zero at the entrance of a suppressor duct, while the reactance should be near the optimum value associated with the least-attenuated mode in a uniform duct.
Baumeister, K. J.↗
The Prediction of Noise Scattered by a Wing/Ducted Fan Configuration
In this proof of concept research, a computational method was developed for predicting the sound field created by the scattering of ducted fan engine noise by a blended wing-body (BWB). It was assumed that all acoustic processes were linear and time harmonic with excitation frequency co. Inflow effects were neglected and no penetration boundary conditions were applied to the engine nacelle and BWB surfaces. A scattering approach was adopted in which the total acoustic field is written as the sum of known incident (from the engine duct) and unknown scattered parts. We further assume that the incident field is independent of the scattered field. Application of the above conditions to the equations of linearized acoustics yields the Helmholtz equation (reduced wave equation) for the scattered pressure with Neumann boundary conditions.
Tweed, John↗
Development of a Performance Portable Non-Equilibrium Plasma Fluid Solver on Adaptive Grids
This presentation will describe the numerical techniques, programming paradigms, verification, and performance of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures. Our plasma fluid model solves the conservation equations for self-consistent electrostatic Poisson, electron and heavy species transport, and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive mesh management library, AMReX (Zhang et al., JOSS, 4 (37) 1370, 2019), and can be built and run on widely available vendor specific GPU architectures (NVIDIA/AMD/Intel). We utilize a non-subcycled second order semi-implicit time-stepping method where all adaptive mesh refinement (AMR) levels are advanced with the same time step. The composite multi-level multigrid solver from within AMReX is used for each of the governing equations that are cast into a Helmholtz equation form. We have also developed a python based chemical mechanism parser framework that uses a similar format as CANTERA (Goodwin et al., Zenodo, 2018) yaml files as input. Our custom parser reads the yaml file and provides C++ files with transport and production rate functions that can be executed on both host (CPU) and device (GPU). We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on low-pressure capacitive and high-pressure streamer discharges. Our initial performance studies indicate 10X speed-up using 20 NVIDIA GPUs versus 200 CPUs for an atmospheric streamer discharge problem solved on a 512 x 1024 x 512 grid.
graphics processing units↗
2D Modeling of Plasma Streamer and Glow Phases at Ammonia-Air Flame Conditions
Streamer and glow plasma phases have been modeled at two thermochemical states of an ammonia-air flame: fresh reactants and burnt products. A new AMReX-based solver has been verified against benchmarks in the literature and has been used to perform these simulations. A Helmholtz-equation based photoionization model with parameters accounting for the presence of NH3 in air has been coupled with the solver to accurately model the streamer propagation phase. A detailed plasma kinetics mechanism has been compiled and used to predict the evolution of electrons, excited states, ions, and radicals during streamer propagation and glow formation. The propagation velocity of streamers was found increase by almost two-fold when the mixture was changed from the fresh reactants to the burnt products. Moreover, vibrational excitation was found to be limited to the streamer body, whereas ionization predominantly occurred at the streamer head, as is expected. Finally, the differences in the pathways of O and H radical production during the streamer propagation and glow phases have been briefly discussed.
ammonia-air flame↗
Fast meta-solvers for 3D complex-shape scatterers using neural operators trained on a non-scattering problem
Three-dimensional target identification using scattering techniques requires high accuracy solutions and very fast computations for real-time predictions in some critical applications. We first train a deep neural operator (DeepONet) to solve wave propagation problems described by the Helmholtz equation in a domain without scatterers but at different wavenumbers and with a complex absorbing boundary condition. We then design two classes of fast meta-solvers by combining DeepONet with either relaxation methods, such as Jacobi and Gauss-Seidel, or with Krylov methods, such as GMRES and BiCGStab, using the trunk basis of DeepONet as a coarse-scale preconditioner. We leverage the spectral bias of neural networks to account for the lower part of the spectrum in the error distribution while the upper part is handled inexpensively using relaxation methods or fine-scale preconditioners. The meta-solvers are then applied to solve scattering problems with different shape of scatterers, at no extra training cost. We first demonstrate that the resulting meta-solvers are shape-agnostic, fast, and robust, whereas the standard standalone solvers may even fail to converge without the DeepONet. We then apply both classes of meta-solvers to scattering from a submarine, a complex three-dimensional problem. We achieve very fast solutions, especially with the DeepONet-Krylov methods, which require orders of magnitude fewer iterations than any of the standalone solvers.
97 MATHEMATICS AND COMPUTING↗
Adaptive Activation Functions Accelerate Convergence in Deep and Physics-informed Neural Networks
We employ adaptive activation functions for regression in deep and physics-informed neural networks (PINNs) to approximate smooth and discontinuous functions as well as solutions of linear and nonlinear partial differential equations. In particular, we solve the nonlinear Klein-Gordon equation, which has smooth solutions, the nonlinear Burgers equation, which can admit high gradient solutions, and the Helmholtz equation. We introduce a scalable hyper-parameter in the activation function, which can be optimized to achieve best performance of the network as it changes dynamically the topology of the loss function involved in the optimization process. The adaptive activation function has better learning capabilities than the traditional one (fixed activation) as it improves greatly the convergence rate, especially at early training, as well as the solution accuracy. To better understand the learning process, we plot the neural network solution in the frequency domain to examine how the network captures successively different frequency bands present in the solution. We consider both forward problems, where the approximate solutions are obtained, as well as inverse problems, where parameters involved in the governing equation are identified. Our simulation results show that the proposed method is a very simple and effective approach to increase the efficiency, robustness and accuracy of the neural network approximation of nonlinear functions as well as solutions of partial differential equations, especially for forward problems. We theoretically prove that in the proposed method, gradient descent algorithms are not attracted to suboptimal critical points or local minima.
machine leaning, Bad minima, Inverse problems, Phy↗
Gravity shear waves atop the cirrus layer of intense convective storms
Recent visual satellite photographs of certain intense convective storms have revealed concentric wave patterns. A model for the generation and growth of these waves is proposed. The proposed initial generating mechanism is similar to the effect noticed when a pebble is dropped into a calm pond. The penetration of the tropopause by overshooting convection is analogous to the pebble's penetration of the water's surface. The model for wave growth involves instability due to the wind shear resulting from the cirrus outflow. This model is based on an equation for the waves' phase speed which is similar to the Helmholtz equation. It, however, does not assume an incompressible atmosphere, but rather assumes density is a logarithmic function of height. Finally, the model is evaluated on the two mid-latitude and three tropical cases. The data indicate that shearing instability may be a significant factor in the appearance of these waves.
Stobie, J. G.↗
A Computer Program for the Computation of Running Gear Temperatures Using Green's Function
A new technique has been developed to study two dimensional heat transfer problems in gears. This technique consists of transforming the heat equation into a line integral equation with the use of Green's theorem. The equation is then expressed in terms of eigenfunctions that satisfy the Helmholtz equation, and their corresponding eigenvalues for an arbitrarily shaped region of interest. The eigenfunction are obtalned by solving an intergral equation. Once the eigenfunctions are found, the temperature is expanded in terms of the eigenfunctions with unknown time dependent coefficients that can be solved by using Runge Kutta methods. The time integration is extremely efficient. Therefore, any changes in the time dependent coefficients or source terms in the boundary conditions do not impose a great computational burden on the user. The method is demonstrated by applying it to a sample gear tooth. Temperature histories at representative surface locatons are given.
GREENS FUNCTIONS↗
Interaction of finite-amplitude sound with air-filled porous materials
The propagation of high intensity sound waves through an air-filled porus material was studied. The material is assumed: (1) to be rigid, incompressible, and homogeneous, and (2) to be adequately described by two properties: resistivity r and porosity. The resulting wave equation is still nonlinear, however, because of the u sgn(u) term in the resistivity. The equation is solved in the frequency domain as an infinite set of coupled inhomogeneous Helmholtz equations, one for each harmonic. An approximate but analytical solution leads to predictions of excess attenuation, saturation, and phase speed reduction for the fundamental component. A more general numerical solution is used to calculate the propagation curves for the higher harmonics. The u sgn(u) nonlinearity produces a cubic distortion pattern; when the input signal is a pure tone, only odd harmonic distortion products are generated.
Nelson, D. A.↗
Computation of Nonlinear Backscattering Using a High-Order Numerical Method
The nonlinear Schrodinger equation (NLS) is the standard model for propagation of intense laser beams in Kerr media. The NLS is derived from the nonlinear Helmholtz equation (NLH) by employing the paraxial approximation and neglecting the backscattered waves. In this study we use a fourth-order finite-difference method supplemented by special two-way artificial boundary conditions (ABCs) to solve the NLH as a boundary value problem. Our numerical methodology allows for a direct comparison of the NLH and NLS models and for an accurate quantitative assessment of the backscattered signal.
Fibich, G.↗
Boundary conditions for the numerical solution of elliptic equations in exterior regions
Elliptic equations in exterior regions frequently require a boundary condition at infinity to ensure the well-posedness of the problem. Examples of practical applications include the Helmholtz equation and Laplace's equation. Computational procedures based on a direct discretization of the elliptic problem require the replacement of the condition on a finite artificial surface. Direct imposition of the condition at infinity along the finite boundary results in large errors. A sequence of boundary conditions is developed which provides increasingly accurate approximations to the problem in the infinite domain. Estimates of the error due to the finite boundary are obtained for several cases. Computations are presented which demonstrate the increased accuracy that can be obtained by the use of the higher order boundary conditions. The examples are based on a finite element formulation but finite difference methods can also be used.
Bayliss, A.↗
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.
Belyaev, Mikhail A. (ORCID:0000000224908887)↗
Parameter estimation technique for boundary value problems by spline collocation method
A parameter-estimation technique for boundary-integral equations of the second kind is developed. The output least-squares identification technique using the spline collocation method is considered. The convergence analysis for the numerical method is discussed. The results are applied to boundary parameter estimations for two-dimensional Laplace and Helmholtz equations.
Kojima, Fumio↗
General theory of spherically symmetric boundary-value problems of the linear transport theory.
A general theory of spherically symmetric boundary-value problems of the one-speed neutron transport theory is presented. The formulation is also applicable to the 'gray' problems of radiative transfer. The Green's function for the purely absorbing medium is utilized in obtaining the normal mode expansion of the angular densities for both interior and exterior problems. As the integral equations for unknown coefficients are regular, a general class of reduction operators is introduced to reduce such regular integral equations to singular ones with a Cauchy-type kernel. Such operators then permit one to solve the singular integral equations by the standard techniques due to Muskhelishvili. We discuss several spherically symmetric problems. However, the treatment is kept sufficiently general to deal with problems lacking azimuthal symmetry. In particular the procedure seems to work for regions whose boundary coincides with one of the coordinate surfaces for which the Helmholtz equation is separable.
Kanal, M.↗
Three-dimensional periodic disturbances acting upon airfoils in cascade
A theory for three dimensional periodic gusts interacting with loaded airfoils was developed. The unsteady disturbances are linearized with respect to the mean potential flow of the airfoils. The vorticity transport equation is integrated analytically in a Lagrangian form. The vorticity vector shows strong variations in its magnitude and wavelength as it interacts with the flow field of a lifting airfoil. The streamwise component of vorticity increases significantly with flow acceleration and flow turning. The rectangular fan approximation is used and a single Helmholtz equation is derived that characterizes the unsteady three dimensional flow field.
Atassi, H.↗