Search NASA⌕ Search

SEARCH · Search NASA

Results for “wave modeling”

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.

At least 469 records · Page 26

On the eigenvector bias of Fourier feature networks: From regression to solving multi-scale PDEs with physics-informed neural networks

Physics-informed neural networks (PINNs) are demonstrating remarkable promise in integrating physical models with gappy and noisy observational data, but they still struggle in cases where the target functions to be approximated exhibit high-frequency or multi-scale features. Here in this work we investigate this limitation through the lens of Neural Tangent Kernel (NTK) theory and elucidate how PINNs are biased towards learning functions along the dominant eigen-directions of their limiting NTK. Using this observation, we construct novel architectures that employ spatio-temporal and multi-scale random Fourier features, and justify how such coordinate embedding layers can lead to robust and accurate PINN models. Numerical examples are presented for several challenging cases where conventional PINN models fail, including wave propagation and reaction–diffusion dynamics, illustrating how the proposed methods can be used to effectively tackle both forward and inverse problems involving partial differential equations with multi-scale behavior.

42 ENGINEERING↗

Encoding of linear kinetic plasma problems in quantum circuits via data compression

We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation Aψ = b to be solved by using the quantum signal processing algorithm. The latter requires encoding of matrix A in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode A in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Effects of Field Line Curvature (FLC) Scattering on Ring Current Dynamics and Isotropic Boundary

Abstract In the ring current dynamics, various loss mechanisms contribute to the ring current decay, including losses to the upper atmosphere through particle precipitation. This study implements the field line curvature (FLC) scattering mechanism in a kinetic ring current model and investigates its role in precipitating ions into the ionosphere during the 17 March 2013 storm. Simulation results indicate that (1) the FLC scattering process exerts on energetic ions on the nightside where the magnetospheric configuration is more stretching. It is more effective on heavy ions (e.g., O + ). These ion losses thereafter lead to a faster recovery of the ring current. (2) The FLC‐associated ion precipitation mainly occurs in the outer region (L > 5 for protons and L > 4.5 for oxygen ions) on the nightside. The O + precipitation takes places in a wider region than protons although its intensity is much lower. Comparisons with POES observations suggest that more proton precipitation is needed in the inner region. This is probably caused by the less stretched configuration in the simulation that prevents more precipitation. It may also imply that other loss process is required in the model such as wave‐particle interactions. (3) The storm time precipitating proton flux of tens of keV due to the FLC scattering sometimes becomes comparable to that of electrons on the nightside, although electrons usually dominate the ionospheric energy deposition from the midnight eastward toward the dayside. (4) The FLC scattering process seems to be capable of explaining the formation of isotropic boundary in the ionosphere during the investigated event.

79 ASTRONOMY AND ASTROPHYSICS↗

Classical discrete time crystals

The spontaneous breaking of time-translation symmetry in periodically driven quantum systems leads to a new phase of matter: the discrete time crystal (DTC). This phase exhibits collective subharmonic oscillations that depend upon an interplay of non-equilibrium driving, many-body interactions and the breakdown of ergodicity. However, subharmonic responses are also a well-known feature of classical dynamical systems ranging from predator–prey models to Faraday waves and a.c.-driven charge density waves. This raises the question of whether these classical phenomena display the same rigidity characteristic of a quantum DTC. In this work, we explore this question in the context of periodically driven Hamiltonian dynamics coupled to a finite-temperature bath, which provides both friction and, crucially, noise. Focusing on one-dimensional chains, where in equilibrium any transition would be forbidden at finite temperature, we provide evidence that the combination of noise and interactions drives a sharp, first-order dynamical phase transition between a discrete time-translation invariant phase and an activated classical discrete time crystal (CDTC) in which time-translation symmetry is broken out to exponentially long timescales. Power-law correlations are present along a first-order line, which terminates at a critical point. We analyse the transition by mapping it to the locked-to-sliding transition of a d.c.-driven charge density wave. To conclude, building upon results from the field of probabilistic cellular automata, we conjecture the existence of classical time crystals with true long-range order, where time-translation symmetry is broken out to infinite times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

In-situ monitoring for liquid metal jetting using a millimeter-wave impedance diagnostic

Abstract This article presents a millimeter-wave diagnostic for the in-situ monitoring of liquid metal jetting additive manufacturing systems. The diagnostic leverages a T-junction waveguide device to monitor impedance changes due to jetted metal droplets in real time. An analytical formulation for the time-domain T-junction operation is presented and supported with a quasi-static full-wave electromagnetic simulation model. The approach is evaluated experimentally with metallic spheres of known diameters ranging from 0.79 to 3.18 mm. It is then demonstrated in a custom drop-on-demand liquid metal jetting system where effective droplet diameters ranging from 0.8 to 1.6 mm are detected. Experimental results demonstrate that this approach can provide information about droplet size, timing, and motion by monitoring a single parameter, the reflection coefficient amplitude at the input port. These results show the promise of the impedance diagnostic as a reliable in-situ characterization method for metal droplets in an advanced manufacturing system.

42 ENGINEERING↗

Ion cyclotron heating at high plasma density in Proto-MPEX

The physics of ion cyclotron heating (ICH) relevant to the steady-state linear machine MPEX (Material Plasma Exposure eXperiment) has been explored in its predecessor, short-pulse device: Proto-MPEX. MPEX will utilize fundamental ICH to increase heat flux at the target and produce ion temperatures and velocity distributions with improved fidelity to those found in a tokamak divertor region, in comparison to those produced by substrate biasing. In the experiments on Proto-MPEX described here, bulk ion temperatures up to ~15 eV have been achieved with 20 kW net ICH power at 6.5 MHz, using ICH heating of a deuterium plasma produced by a helicon plasma source. The heat flux at the target has been observed to increase throughout the plasma cross section, including in the core region. Core $T_i$ and target heat flux are observed to scale linearly with injected ICH power. Further, measurements of plasma loading and target heat flux as a function of the magnetic field strength at the antenna, together with modeling of the wave propagation from the antenna to the ion cyclotron resonance using the ANTENA and COMSOL codes with a warm plasma dielectric tensor, indicate that power is coupled to the core plasma via fast wave excitation of a kinetic Alfvén wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Black hole echoes

In this work, we consider a very simple model for gravitational wave echoes from black hole merger ringdowns which may arise from local Lorentz symmetry violations that modify graviton dispersion relations. If the corrections are sufficiently soft so they do not remove the horizon, the reflection of the infalling waves which trigger the echoes is very weak. As an example, we look at the dispersion relation of a test scalar field corrected by rotonlike operators depending only on spatial momenta, in Gullstrand-Painlevé coordinates. The near-horizon regions of a black hole do become reflective, but only very weakly. The resulting “bounces” of infalling waves can yield repetitive gravity wave emissions but their power is very small. This implies that to see any echoes from black holes we really need an egregious departure from either standard GR or effective field theory, or both. One possibility to realize such strong echoes is the recently proposed classical firewalls which replace black hole horizons with material shells surrounding timelike singularities.

79 ASTRONOMY AND ASTROPHYSICS↗

Interpretable AI forecasting for numerical relativity waveforms of quasicircular, spinning, nonprecessing binary black hole mergers

We present a deep-learning artificial intelligence model (AI) that is capable of learning and forecasting the late-inspiral, merger and ringdown of numerical relativity waveforms that describe quasicircular, spinning, nonprecessing binary black hole mergers. We used the NRHybSur3dq8 surrogate model to produce train, validation and test sets of ℓ = |m| = 2 waveforms that cover the parameter space of binary black hole mergers with mass ratios q ≤ 8 and individual spins |s$^{z}_ {(1,2)}$| ≤ 0.8. These waveforms cover the time range t ∊ [-5000 M, 130 M], where t = 0M marks the merger event, defined as the maximum value of the waveform amplitude. We harnessed the ThetaGPU supercomputer at the Argonne Leadership Computing Facility to train our AI model using a training set of 1.5 million waveforms. We used 16 NVIDIA DGX A100 nodes, each consisting of 8 NVIDIA A100 Tensor Core GPUs and 2 AMD Rome CPUs, to fully train our model within 3.5 h. Our findings show that artificial intelligence can accurately forecast the dynamical evolution of numerical relativity waveforms in the time range t ∊ [-100 M, 130 M]. Sampling a test set of 190,000 waveforms, we find that the average overlap between target and predicted waveforms is ≳99% over the entire parameter space under consideration. We also combined scientific visualization and accelerated computing to identify what components of our model take in knowledge from the early and late-time waveform evolution to accurately forecast the latter part of numerical relativity waveforms. This work aims to accelerate the creation of scalable, computationally efficient and interpretable artificial intelligence models for gravitational wave astrophysics

79 ASTRONOMY AND ASTROPHYSICS↗

Spatiotemporal mapping of mesoscopic liquid dynamics

The study of liquid dynamics at mesoscopic scales is still strewn with difficulty due to limitations in theory and experiment. Historically, significant attention has been given to the analysis of space-time correlation functions and their frequency-Fourier transforms at a few discrete wave numbers. The massive computing power afforded by modern high performance computing clusters and the advent of a wide-angle neutron spin-echo spectrometer, however, have unlocked a more intuitive and fruitful approach to this problem. Using molecular dynamics simulations, here we demonstrate the benefits of spatiotemporally mapping intermediate scattering functions on a dense grid of correlation times and wave numbers. Four model systems are investigated: a Lennard-Jones liquid, a coarse-grained bead-spring polymer, a molten sodium chloride, and a poly(ethylene oxide) melt. We show that the spatiotemporal mapping approach is particularly useful for elucidating the mesoscopic dynamics in these liquids, where several underlying mechanisms, such as molecular relaxations, hydrodynamic modes, and nonhydrodynamic excitations, are potentially at play. Therefore, compared to the traditional method, direct visualization of density space-time correlation functions on two-dimensional color maps permits appraisals of complicated dynamical behavior at mesoscales in a global manner. For example, the scaling relations between space and time for different types of molecular motions can be straightforwardly identified on these plots, without any model-dependent analysis. Additionally, we show how theoretical ideas regarding collective mesoscopic dynamics, such as the classical hydrodynamic theory, the convolution approximation, and a recently proposed phenomenological model, can be discussed in terms of the global features of spatiotemporal maps of intermediate scattering functions. The new perspective offered by the spatiotemporal mapping method should prove useful for the study of liquid dynamics in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Constraining Axion Inflation with Gravitational Waves across 29 Decades in Frequency

We demonstrate that gravitational waves generated by efficient gauge preheating after axion inflation generically contribute significantly to the effective number of relativistic degrees of freedom $N_{\text{eff}}$. In this work we show that, with existing Planck limits, gravitational waves from preheating already place the strongest constraints on the inflaton’s possible axial coupling to Abelian gauge fields. We demonstrate that gauge preheating can completely reheat the Universe regardless of the inflationary potential. Further, we quantify the variation of the efficiency of gravitational wave production from model to model and show that it is correlated with the tensor-to-scalar ratio. In particular, when combined with constraints on models whose tensor-to-scalar ratios would be detected by next-generation cosmic microwave background experiments, $\textit{r}$ ≳ 10 –3 , constraints from $N_{\text{eff}}$ will probe or rule out the entire coupling regime for which gauge preheating is efficient.

79 ASTRONOMY AND ASTROPHYSICS↗

Gemma

Gemma is a linear electromagnetic code developed at Sandia National Laboratories based on the method of moments in the frequency domain, with the objective of supporting next-generation high-performance computing (HPC) platforms. Gemma provides full-wave computational electromagnetic modeling. As an example, Gemma models electromagnetic radiation (EMR) environments, including those resulting from radar (both hostile and friendly) and high-power-electromagnetic (HPEM) sources. Gemma refers to a suite of source code and scripts that aid in compiling, running, and testing the primary Gemma executable, which can be compiled for diverse platforms including CPUs, GPUs, and MICs by leveraging Sandia's Kokkos library.SAND2019-15221 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Zinser, Brian↗

Energy density functional and sensitivity of energies of giant resonances to bulk nuclear matter properties

The development of a modern and more realistic nuclear energy density functional (EDF) for accurate predictions of properties of nuclei is the subject of enhanced activity, since it is very important for the study of properties of nuclear matter (NM), giant resonances and, in particular, of properties of rare nuclei with unusual neutron-to-proton ratios. Here, we provide a short review of the current status of the nuclear EDF and the theoretical results obtained for properties of nuclei and nuclear matter. We will first describe a method for determining the parameters of the EDF, associated with the Skyrme type effective interaction, by carrying out a Hartree-Fock based fit to extensive set of data of ground state properties and constraints. Next we will describe the fully self-consistent Hartree-Fock (HF) based random-phase-approximation (RPA) theory for calculating the strength functions S(E) and centroid energies E CEN of giant resonances and the folding model (FM) distorted wave Вorn approximation (DWBA) to calculate the excitation cross section of giant resonances by α scattering. Then we will provide results for: (i) the Skyrme parameters of the KDE0v1 EDF; (ii) consequences of violation of self-consistently in HF-based RPA; (iii) FM-DWBA calculation of excitation cross section; (iv) values of the E CEN of isoscalar and isovector giant resonances of multipolarities L=0–3 for a wide range of spherical nuclei, using 33 EDFs associated with standard form of the Skyrme type interactions, commonly employed in the literature; and (v) the sensitivities E CEN of the giant resonances to bulk properties of NM. We also determine constraints on NM properties, such as the incompressibility coefficient and effective mass, by comparing with experimental data on E CEN of giant resonances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Nonlocal Operator Learning with Uncertainty Quantification

The goal of this work is to develop a Bayesian framework to characterize the uncertainty of material response when using a nonlocal, homogenized model to describe wave propagation through heterogeneous, disordered materials. Our approach is based on an operator regression technique combined with Bayesian optimization, through which the nonlocal kernel for a specific disordered microstructure is investigated.

36 MATERIALS SCIENCE↗

Multimode Metastructures: Novel Hybrid 3D Lattice Topologies

With the rapid proliferation of additive manufacturing and 3D printing technologies, architected cellular solids including truss-like 3D lattice topologies offer the opportunity to program the effective material response through topological design at the mesoscale. The present report summarizes several of the key findings from a 3-year Laboratory Directed Research and Development Program. The program set out to explore novel lattice topologies that can be designed to control, redirect, or dissipate energy from one or multiple insult environments relevant to Sandia missions, including crush, shock/impact, vibration, thermal, etc. In the first 4 sections, we document four novel lattice topologies stemming from this study: coulombic lattices, multi-morphology lattices, interpenetrating lattices, and pore-modified gyroid cellular solids, each with unique properties that had not been achieved by existing cellular/lattice metamaterials. The fifth section explores how unintentional lattice imperfections stemming from the manufacturing process, primarily sur face roughness in the case of laser powder bed fusion, serve to cause stochastic response but that in some cases such as elastic response the stochastic behavior is homogenized through the adoption of lattices. In the sixth section we explore a novel neural network screening process that allows such stocastic variability to be predicted. In the last three sections, we explore considerations of computational design of lattices. Specifically, in section 7 using a novel generative optimization scheme to design novel pareto-optimal lattices for multi-objective environments. In section 8, we use computational design to optimize a metallic lattice structure to absorb impact energy for a 1000 ft/s impact. And in section 9, we develop a modified micromorphic continuum model to solve wave propagation problems in lattices efficiently.

36 MATERIALS SCIENCE↗

LLNL FESP Theory Highlights: October 2024

I. Novikau, I. Y. Dodin, E. A. Startsev, I. Joseph, Quantum algorithms for simulating dissipative linear and nonlinear dynamics of plasmas. Invited talk at the 66th Annual Meeting of the APS Division of Plasma Physics, Atlanta, Georgia. Novikau I., Dodin I.Y., Startsev E.A., Encoding of linear kinetic plasma problems in quantum circuits via data compression, Journal of Plasma Physics. 2024;90(4):805900401, doi:10.1017/S0022377824000795. We propose an algorithm for encoding linear kinetic plasma problems in quantum circuits. The focus is on modelling electrostatic linear waves in a one-dimensional Maxwellian electron plasma. The waves are described by the linearized Vlasov–Ampère system with a spatially localized external current that drives plasma oscillations. This system is formulated as a boundary-value problem and cast in the form of a linear vector equation to be solved by using the quantum signal processing algorithm. The latter requires encoding of a matrix in a quantum circuit as a sub-block of a unitary matrix. We propose how to encode in a circuit in a compressed form and discuss how the resulting circuit scales with the problem size and the desired precision.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

THz Radiation Generation to Enable Internal Magnetic Field Measurement of Burning Plasmas (Final Scientific Report)

A one-year experimental program was performed using the Multi-TeraWatt (MTW) laser at the Laboratory for Laser Energetics (LLE) to investigate the efficiency and physical process of laser-driven—wire THz generation. This process was evaluated as a potential high-repetition-rate source of powerful THz pulses for use in THz pulsed polarimetry of a moderate-density magnetic-confinement fusion concept developed by TAE. The goal was to develop a source for short pulses in the THz regime (150 μm wavelength) containing 10’s of mJ. The campaign tested the impact of Cu wire diameter (25, 75, 160 μm), length (250 – 750 μm), and laser irradiation orientation on THz production efficiency in comparison to 20-μm thick Cu foils. Of these targets, the 160 µm diameter, 750 µm long wires with sharpened tips and irradiated at the base were found to produce the highest energy and shortest wavelength THz pulses, exceeding 2.5 Joule/sr and meeting the requirements for the polarimetry diagnostic. The brightest THz signal was radiated in the direction of the sharpened wire tip. A combination of two models – a traveling wave antenna combined with coherent transition radiation – describes the emission pattern, pulse energy, and spectrum observed in the experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetohydrodynamic Wave Mode Conversion in a Pseudostreamer Structure

The interaction of a fast magnetoacoustic wave with a magnetic null point in a pseudostreamer magnetic configuration is modeled. The incident fast wave is excited by an impulsive velocity driver located in the corona outside the Alfvén-acoustic equipartition layer. The fast wave front experiences refraction and wraps around the null point, while the front segment approaching the null point along the magnetic bisector passes through it. In the equipartition layer surrounding the null point, part of the fast wave energy is converted into slow magnetoacoustic waves. Slow waves propagate outward from the null point along the magnetic separatrices. The fronts of slow waves exhibit a characteristic chevron shape. The fast-to-slow mode conversion efficiency increases with the width of the initial perturbation, following a power-law dependence with an index of approximately 1.06. For initial wave amplitudes below about 10% of the local Alfvén speed, the dependence of the conversion efficiency on the initial amplitude is weak. After interacting with the null point, the amplitudes of the fast waves decrease significantly owing to energy conversion and the spreading of wave energy over a broad spatial domain. In contrast, the energy of slow waves remains concentrated in the vicinity of the magnetic separatrices, resulting in less pronounced geometric attenuation of these waves. These findings are consistent with observational reports of the interaction of a coronal fast wave with a pseudostreamer.

Zhong, Yu [Univ. of Warwick, Coventry (United King↗

Snowmass2021 Cosmic Frontier White Paper: Cosmological Simulations for Dark Matter Physics

Over the past several decades, unexpected astronomical discoveries have been fueling a new wave of particle model building and are inspiring the next generation of ever-more-sophisticated simulations to reveal the nature of Dark Matter (DM). This coincides with the advent of new observing facilities coming online, including JWST, the Rubin Observatory, the Nancy Grace Roman Space Telescope, and CMB-S4. The time is now to build a novel simulation program to interpret observations so that we can identify novel signatures of DM microphysics across a large dynamic range of length scales and cosmic time. This white paper identifies the key elements that are needed for such a simulation program. We identify areas of growth on both the particle theory side as well as the simulation algorithm and implementation side, so that we can robustly simulate the cosmic evolution of DM for well-motivated models. We recommend that simulations include a fully calibrated and well-tested treatment of baryonic physics, and that outputs should connect with observations in the space of observables. We identify the tools and methods currently available to make predictions and the path forward for building more of these tools. A strong cosmic DM simulation program is key to translating cosmological observations to robust constraints on DM fundamental physics, and provides a connection to lab-based probes of DM physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗