Out-of-distributional risk bounds for neural operators with applications to the Helmholtz equation
Not provided.
Computer Science↗
SEARCH · Search NASA
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.
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↗
Reducing Frequency Bias of Fourier Neural Operators in 3D Seismic Wavefield Simulations Through Multistage Training
The recent development of neural operator (NeurOp) learning for solutions to the elastic wave equation shows promising results and provides the basis for fast large-scale simulations for different seismological applications. In this article, we use the Fourier neural operator (FNO) model to directly solve the 3D Helmholtz wave equation for fast seismic ground-motion simulations on different frequencies and show the frequency bias of the FNO model, that is, it learns the lower frequencies better comparing to the higher frequencies. To reduce the frequency bias, we adopt the multistage FNO training, that is, after training a stage 1 FNO model for estimating the ground motion, we use a second FNO model as the stage 2 to learn from the residual, which greatly reduced the errors on the higher frequencies. By adopting this multistage training, the FNO models show reduced biases on higher frequencies, which enhanced the overall results of the ground-motion simulations. Thus the multistage training FNO improves the accuracy and realism of the ground-motion simulations.
earthquakes↗
Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods
Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.
Axisymmetric system↗
An FFT-based micromechanical model for gradient enhanced brittle fracture
Damage models incorporated within FFT-based micromechanical methods have received much attention recently because of the need to better understand and predict brittle and ductile fracture. An important aspect of a damage model is non-local regularization, which removes the mesh dependence of the predictions that otherwise become physically unacceptable upon grid refinement. In this work, the Helmholtz-type equation for non-local gradient regularization of a damage model on a distorted grid is solved using an FFT-based approach. Further, the resulting system of equations is solved using the Jacobi iterative method. The model is applied to simulate brittle fracture of an intermetallic. The influence of the time and space discretization, the length-scale parameter, and intermetallic crystallographic orientation on crack evolution is studied.
36 MATERIALS SCIENCE↗
Model Development and Analysis of a High-Fidelity Neutron Transport Sensor: The Quadrupole Detector Concept for Measurement of the Neutron Flux Gradient
Accurate reconstruction of the neutron flux distribution within a reactor core is essential for safe and efficient reactor operation. Traditional power shape synthesis in Light Water Reactors relies on hundreds of in-core detectors. However, this approach becomes impractical for Advanced Reactors and Microreactors due to limited space and harsh environments. To address this challenge, we propose a data-driven methodology that combines high-fidelity modeling with real-time ex-core sensor measurements, enabling the reconstruction of core power distribution while minimizing the reliance on intrusive in-core instrumentation. This project began in FY24 and achieved two initial milestones: (1) the definition of a three-year development plan for a Digital Twin framework and (2) the development of high-fidelity neutronics models of the Purdue University Reactor One (PUR-1) using both MCNP6 and OpenMC. The PUR-1 reactor, a zero-power facility, was selected due to its suitability for neutronics-focused modeling and the availability of experimental data for validation. Both models were benchmarked using neutron flux measurements obtained from irradiated gold foils, which were strategically placed within the core during a dedicated campaign in July 2024. This report marks the continuation and completion of those foundational tasks. The OpenMC model has been refined (improved geometric accuracy, expanded cross-section libraries, and refined sampling) and validated using additional experimental data. An updated sensor design—based on quadrupole configuration—was designed to measure both ex-core flux and its spatial gradient. These measurements will serve as inputs to a neural network-based reconstruction algorithm. Finally, the methodology was demonstrated on a two-dimensional test case representative of the heterogeneous material composition of the PUR-1 reactor core. A neural network implementation of the Kirchhoff-Helmholtz integral equation was employed to solve the boundary value problem using peripheral sensor measurements. The preliminary results confirm the strong potential of the proposed approach for accurate and minimally invasive neutron flux reconstruction.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Speed of Sound Measurements of R-1130(E) and an Azeotropic Blend of R-1336mzz(Z)/1130(E)
Sound speed data measured using a dual-path pulse-echo instrument are reported for pure trans-1,2-dichloroethene (R-1130(E)) and an azeotropic blend of cis-1,1,1,4,4,4-hexafluorobutene (R-1336mzz(Z)) and R-1130(E) with a composition of 74.8 mass % R-1336mzz(Z) with the balance being R-1130(E). The azeotropic blend of R-1336mzz(Z)/1130(E) is classified as R-514A in ANSI/ASHRAE standard 34. Liquid phase speed of sound data are reported from just above the saturation pressure of pure R-1130(E) or the bubble point pressure of R-514A to a maximum pressure of 26.7 MPa. The relative combined expanded uncertainty in the speed of sound varies from 0.032 % to 0.148 % with the greatest deviations occurring at the lowest sound speeds. At present, no reference Helmholtz-energy-explicit equation of state (EOS) is available for R-1130(E). Therefore, the reported data for pure R-1130(E) are compared to an extended corresponding states (ECS) model. Deviations between the pure R-1130(E) sound speed data and the ECS model were found to be consistently negative ranging between − 4.1 % and − 3.5 %. The R-514A data are compared to a multifluid model inclusive of the established reference Helmholtz-energy-explicit EOS for R-1336mzz(Z) and ECS model for R-1130(E) with estimated binary interaction parameters. Deviations between the experimental speed of sound data and the multifluid model were also found to be consistently negative. However, deviations from the multifluid model were found to be as great as − 17.1 %. The large deviations from the ECS model and multifluid model underscore the need for a robust Helmholtz-energy-explicit EOS for R-1130(E).
42 ENGINEERING↗
Speed of Sound Measurements of Binary Mixtures of trans -1,2-Difluoroethylene (R-1132(E)) with Difluoromethane (R-32) or 2,3,3,3-Tetrafluoropropene (R-1234yf)
The speed of sound of two binary mixtures containing 0.336 mole fraction trans-1,2- difluoroethylene (R-1132(E)) with difluoromethane (R-32) and 0.435 mole fraction R-1132(E) with 2,3,3,3-tetrafluoropropene (R-1234yf) were measured with a dual-path pulse-echo instrument. The speed of sound was measured along several pseudo-isochores for each blend at temperatures ranging from 230 K to 325 K for blends with R-32 and to 342 K for blends with R- 1234yf. Measurements started at pressures just above each mixtures bubble point pressure and were limited to 8 MPa to avoid potential disproportionation reactions of the R-1132(E). The data were compared to multi-fluid models incorporating Helmholtz-energy-explicit equations of state (EOS) for each pure fluid.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Bubble Point Measurements of cis-1,1,1,4,4,4-Hexafluorobutene [R-1336mzz(Z)] + trans-1,2-Dichloroethene [R-1130(E)] mixtures
Saturation pressures of pure R-1336mzz(Z) and R-1130(E) and bubble point pressures of three R-1336mzz(Z)/1130(E) blends were measured from 265 K to 360 K. For each pure refrigerant or refrigerant blend, a total of twenty unique saturation pressures or bubble points were measured. In total 100 unique state points were obtained. Presently, no Helmholtz-energy-explicit type equation of state (EoS) is available for R-1130(E). While an extended corresponding states EoS for R-1130(E) is available to estimate the properties of the R-1336mzz(Z)/1130(E) blend, this model does not resolve the azeotropic behavior of the mixture. Therefore, the perturbed-chain statistical-associating fluid theory (PC-SAFT) EoS is used to model the vapor–liquid equilibria of the R-1336mzz(Z)/1130(E) blend. PC-SAFT model parameters are reported, and the overall performance of the model is characterized by deviations from the experimental data.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Equation of State for the Thermodynamic Properties of Trans-1,2-dichloroethene [R-1130(E)]
We present an empirical equation of state in terms of the Helmholtz energy for trans-1,2-dichloroethene [R-1130(E)]. The range of validity is from the triple-point temperature, 223.31 K to 525 K with pressures up to 30 MPa. It may be used to calculate all thermodynamic properties in the fluid phase, including liquid, gas, and supercritical regions. Comparisons are given with existing literature data and estimated uncertainties are provided. In addition, checks were made for correct extrapolation behavior so that the equation behaves in a physically realistic manner when used outside of its range of validity, enabling its use in mixture models. The estimated uncertainties (at a k = 2 or 95 % level of confidence) are based on comparisons with critically assessed data and are 0.25 % for vapor pressure for temperatures in the range 300 K < T < 454 K, rising to 1.5 % as the temperature decreases from 300 K to 265 K. For density in the liquid phase the estimated uncertainty is 0.14 % for temperatures 270 K < T < 410 K and for pressures up to 30 MPa. For the vapor phase the estimated uncertainty in density is 3 %. The uncertainty for liquid-phase heat capacity is 1 % at atmospheric pressure over the temperature range 268 K < T < 309 K, and the uncertainty for the speed of sound in the liquid phase is 0.25 % for temperatures 230 K < T < 420 K and for pressures up to 30 MPa. The uncertainties are larger outside of these specified ranges and in the critical region.
1,2-Dichloroethene↗
A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing
A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.
Buoyancy–drag↗
Physics-constrained machine learning for electrodynamics without gauge ambiguity based on Fourier transformed Maxwell’s equations
We utilize a Fourier transformation-based representation of Maxwell’s equations to develop physics-constrained neural networks for electrodynamics without gauge ambiguity, which we label the Fourier–Helmholtz–Maxwell neural operator method. In this approach, both of Gauss’s laws and Faraday’s law are built in as hard constraints, as well as the longitudinal component of Ampère–Maxwell in Fourier space, assuming the continuity equation. An encoder–decoder network acts as a solution operator for the transverse components of the Fourier transformed vector potential, $\hat{A}_⟂(k,t)$, whose two degrees of freedom are used to predict the electromagnetic fields. This method was tested on two electron beam simulations. Among the models investigated, it was found that a U-Net architecture exhibited the best performance as it trained quicker, was more accurate and generalized better than the other architectures examined. We demonstrate that our approach is useful for solving Maxwell’s equations for the electromagnetic fields generated by intense relativistic charged particle beams and that it generalizes well to unseen test data, while being orders of magnitude quicker than conventional simulations. We show that the model can be re-trained to make highly accurate predictions in as few as 20 epochs on a previously unseen data set.
97 MATHEMATICS AND COMPUTING↗
A phase-field fracture formulation for generalized standard materials: The interplay between thermomechanics and damage
Accurately modeling fracture of ductile materials poses open challenges in the field of computational mechanics due to the multiphysics nature of their failure processes. Integrating the interplay between thermodynamics and damage into ductile fracture models is vital for predicting critical failure modes. Here, in this paper, we develop a versatile phase-field (PF) framework for modeling ductile fracture, taking into account finite-strain elasto-plasticity. The framework stems from a variational formulation of constitutive relations for generalized standard materials (GSMs), whose response is described by a Helmholtz free energy and a dissipation pseudo-potential. Its variational structure is based on a minimum principle for a functional that expresses the sum of power densities for reversible and irreversible processes. By minimizing this functional with a constraint on a von Mises yield function, we derive the evolution equation for the equivalent plastic strain and an associative flow rule. This constrained optimization problem is analytically solved for a wide class of thermo-viscoplasticity models. The key innovations of the current work include (i) a cubic plastic degradation function that accounts for a non-vanishing damage-dependent yield stress, (ii) closed-form expressions of the Helmholtz free energy and dissipation pseudo-potential for three thermo-viscoplasticity models, (iii) an extended Johnson–Cook plasticity model with a nonlinear hardening law, and (iv) a plastic work heat source that depends on the plastic degradation function and a variable Taylor–Quinney (TQ) coefficient. The capabilities of the proposed framework are tested with the aid of four ductile fracture problems, including the Sandia Fracture Challenge. In each of these problems, we examine the evolution of relevant field variables such as the PF order parameter, the equivalent plastic strain, the temperature, and the internal power dissipation density, in addition to the overall structural response quantified by the force–displacement curve. These numerical studies demonstrate that the proposed framework effectively represents ductile fracture, yielding computational results that exhibit good agreement with experimental data.
36 MATERIALS SCIENCE↗
Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor mixing induced by power-law accelerations in the small Atwood number limit
Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing turbulent Rayleigh–Taylor mixing driven by a temporal power-law acceleration are derived in the small Atwood number (Boussinesq) limit. The solutions generalize those previously derived for constant acceleration Rayleigh–Taylor mixing for models based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar variance S and its dissipation rate χ [O. Schilling, “Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit,” Phys. Fluids 33, 085129 (2021)]. The turbulent fields are expressed in terms of the model coefficients and power-law exponent, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. Mixing layer growth parameters and other physical observables are obtained explicitly as functions of the model coefficients and parameterized by the exponent of the power-law acceleration. Values for physical observables in the constant acceleration case are used to calibrate the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to a canonical (i.e., constant acceleration) Rayleigh–Taylor turbulent flow. The calibrated four-equation model is then used to numerically reconstruct the mean and turbulent fields, and turbulent equation budgets across the mixing layer for several values of the power-law exponent. Finally, the reference solutions derived here can be used to understand the model predictions for strongly accelerated or decelerated Rayleigh–Taylor mixing in the large Reynolds number limit.
42 ENGINEERING↗
Thermodynamic properties of “near-perfect” gas
When applied thermodynamics requires simple representations of thermodynamic quantities, one often finds quantities such as pressure and the thermal energy density expressed as a monomial linear product of density and temperature power laws. This is a simple generalization of the product of density and temperature used for the perfect gas that admits a broader range of thermodynamic behavior into analytic fluid calculations or provides an analytic form that is readily fit to tabulated equation-of-state data for an arbitrary material. This paper reviews the thermodynamic properties of this generalized perfect-gas model, treating it as a class of “near-perfect-gas” models that are defined and unified here, for the first time, through a shared Helmholtz free energy. This generalization from perfect to near-perfect-gas models preserves some perfect-gas properties (e.g., constant adiabatic exponents) but not others (e.g., constant specific heats) and is constrained by thermodynamic consistency. It is important to be aware of the properties of this class of near-perfect-gas models when using them as surrogates for real nonperfect materials that may have properties that cannot be captured by near-perfect-gas models.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes
We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.
79 ASTRONOMY AND ASTROPHYSICS↗