SEARCH · Search NASA

Results for “WAVE 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

DOE OSTI2021

The research into the propagation law of the shock wave of a gas explosion inside a building

Based on the dissipation rate conservation equations of turbulent kinetic energy in the k‐ ε turbulence model, a complicated three‐dimensional finite element model of a kitchen filled with gas mixture is developed by using the open source field operation and manipulation (OpenFOAM). Two representative kitchens were used to investigate the propagation law of the shock wave of a gas explosion inside a building by considering the key characteristics of the blast shock wave. The influence of some crucial parameters, such as initial conditions and kitchen parameters, on the properties of the blast shock wave is investigated. The basic steps to predict the peak pressure of the blast shock wave are given in consideration of the initial condition and the kitchen whilst the injury effect of the blast shock wave on the humans and animals is evaluated. The research results indicate that the pressure time history and the peak pressure space distribution are greatly influenced by the kitchen design layout. The coupled interaction between the initial temperature and gas volume concentration, especially at the upper and lower explosion limits of the gas, significantly affects the peak pressure. The peak pressure varies significantly with the opening and the buffer; however, it has little relation with the width, length, and height of the kitchen. The proposed method can accurately and effectively predict the peak pressure of the blast shock wave inside buildings. In terms of the peak pressure space distribution of the explosion shock wave, the peak pressure is much higher than the threshold of the killing pressure, which is unsafe for the humans and animals in the building.

Lin, Shu-Chao↗

DOE OSTI2022

Quantifying modeling uncertainties when combining multiple gravitational-wave detections from binary neutron star sources

Here, with the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in the near future. The combined analysis of these gravitational-wave signals offers the possibility to constrain the neutron-star radius and the equation of state of dense nuclear matter with unprecedented accuracy. However, it is crucial to ensure that uncertainties inherent in the gravitational-wave models will not lead to systematic biases when information from multiple detections is combined. To quantify waveform systematics, we perform an extensive simulation campaign of binary neutron-star sources and analyze them with a set of four different waveform models. For our analysis with 38 simulations, we find that statistical uncertainties in the neutron-star radius decrease to ± 250 m (2% at 90% credible interval) but that systematic differences between currently employed waveform models can be twice as large. Hence, it will be essential to ensure that systematic biases will not become dominant in inferences of the neutron-star equation of state when capitalizing on future developments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Computation of high-frequency magnetoelastic waves in layered materials

Here, the direct calculation of magnetoelastic wave dispersion in layered media is presented using an efficient, accurate computational technique. The governing, coupled equations for elasticity and magnetism, the Navier and Landau-Lifshitz equations, respectively, are linearized to form a quadratic eigenvalue problem that determines a complex web of wave-number–frequency dispersion branches and their corresponding mode profiles. Numerical discretization of the eigenvalue problem via a spectral collocation method (SCM) is employed to determine the complete dispersion maps for both a single, finite-thickness magnetic layer and a finite magnetic-nonmagnetic double-layer. The SCM, previously used to study elastic waves in nonmagnetic media, is fast, accurate, and adaptable to a variety of sample configurations and geometries. Emphasis is placed on the extremely high-frequency regimes being accessed in ultrafast magnetism experiments. The dispersion maps and modes provide insight into how energy propagates through the coupled system, including how energy can be transferred between elastic- and magnetic-dominated waves as well as between different layers. The numerical computations for a single layer are further understood by a simplified analytical calculation in the high-frequency, exchange-dominated regime where the resonance condition required for energy exchange (an anticrossing) between quasi-elastic and quasi-magnetic dispersion branches is determined. Nonresonant interactions are shown to be well approximated by the dispersion of uncoupled elastic and magnetic waves. The methods and results provide fundamental theoretical tools to model and understand current and future magnetic devices powering spintronic innovation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Sampling electronic structure quadratic unconstrained binary optimization problems (QUBOs) with Ocean and Mukai solvers

The most advanced D-Wave Advantage quantum annealer has 5000+ qubits, however, every qubit is connected to a small number of neighbors. As such, implementation of a fully-connected graph results in an order of magnitude reduction in qubit count. To compensate for the reduced number of qubits, one has to rely on special heuristic software such as qbsolv, the purpose of which is to decompose a large quadratic unconstrained binary optimization (QUBO) problem into smaller pieces that fit onto a quantum annealer. In this work, we compare the performance of the open-source qbsolv which is a part of the D-Wave Ocean tools and a new Mukai QUBO solver from Quantum Computing Inc. (QCI). The comparison is done for solving the electronic structure problem and is implemented in a classical mode (Tabu search techniques). The Quantum Annealer Eigensolver is used to map the electronic structure eigenvalue-eigenvector equation to a QUBO problem, solvable on a D-Wave annealer. We find that the Mukai QUBO solver outperforms the Ocean qbsolv with one to two orders of magnitude more accurate energies for all calculations done in the present work, both the ground and excited state calculations. This work stimulates the further development of software to assist in the utilization of modern quantum annealers.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Verification of a fully kinetic ion model for electromagnetic simulations of high-frequency waves in toroidal geometry

For the research of high-frequency electromagnetic waves in tokamaks, an electromagnetic simulation model, in which the ion dynamics is described by a six-dimensional Vlasov equation and the electron dynamics is described by a drift kinetic equation, is formulated and implemented in the global gyrokinetic toroidal code (GTC). Analytic dispersion relations are derived in reduced systems and compared with various theories to verify the model. Linear simulations of a generalized ion Bernstein wave and ion cyclotron emission are verified by comparing the GTC simulation results with analytic dispersion relation theory and magnetoacoustic cyclotron instability theory, respectively, in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Energy density and spectrum of single-cycle and sub-cycle electromagnetic pulses

Based on the exact solution of the Maxwell equations in the form of a collapsing spherical vector wave specified by an arbitrary function of time, we have calculated the maximum energy density that can be achieved when focusing extremely short pulses of various shapes. It is shown that our earlier formula expressing the maximum energy density in terms of the spectrum parameters for Gaussian quasi-monochromatic pulses is approximately (with an accuracy of ±40 %) valid for the main types of extremely short pulses. (extremely short laser pulses)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Bayesian inference of fine features of the nuclear equation of state from future neutron star radius measurements to 0.1 km accuracy

To more precisely constrain the equation of state (EOS) of supradense neutron-rich nuclear matter, future high-precision x-ray and gravitational wave observatories are proposed to measure the radii of neutron stars (NSs) with an accuracy better than about 0.1 km. However, it remains unclear what particular aspects (other than the stiffness generally spoken of in the literature) of the EOS and to what precision they will be better constrained. In this work, within a Bayesian framework using a metamodel EOS for NSs, we infer the posterior probability distribution functions (PDFs) of incompressibility K 0 and skewness J 0 of symmetric nuclear matter (SNM) as well as the slope L, curvature K sym , and skewness J sym characterizing the density dependence of nuclear symmetry energy E sym ⁡(ρ), respectively, from mean values of NS radii consistent with existing observations and an expected accuracy Δ⁢R ranging from about 1.0 to 0.1 km. Here, we found that (1) the Δ⁢R has little effect on inferring the stiffness of SNM at suprasaturation densities, (2) smaller Δ⁢R reveals more accurately not only the PDFs but also pairwise correlations among parameters characterizing high-density E sym ⁡(ρ), (3) a double-peak feature of the PDF(K sym ) corresponding to the strong K sym – J sym and K sym – L anticorrelations is revealed when Δ⁢R is less than about 0.2 km, and the locations of the two peaks are sensitive to the maximum value of J sym reflecting the stiffness of E sym ⁡(ρ) above about 3 times the saturation density ρ 0 of SNM, and (4) the high-precision radius measurement for canonical NSs is more useful than that for massive ones for constraining the EOS of nucleonic matter around (2–3)⁢ρ 0 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Diverse Manifestations of Electron-Phonon Coupling in a Kagome Superconductor

Recent angle-resolved photoemission spectroscopy (ARPES) experiments on the kagome metal CsV_{3}Sb_{5} revealed distinct multimodal dispersion kinks and nodeless superconducting gaps across multiple electron bands. The prominent photoemission kinks suggest a definitive coupling between electrons and certain collective modes, yet the precise nature of this interaction and its connection to superconductivity remain to be established. Here, employing the state-of-the-art ab initio many-body perturbation theory computation, we present direct evidence that electron-phonon (e-ph) coupling induces the multimodal photoemission kinks in CsV_{3}Sb_{5}, and profoundly, drives the nodeless s-wave superconductivity, showcasing the diverse manifestations of the e-ph coupling. Our calculations well capture the experimentally measured kinks and their fine structures, and reveal that vibrations from different atomic species dictate the multimodal behavior. Results from anisotropic GW-Eliashberg equations predict a phonon-mediated superconductivity with nodeless s-wave gaps, in excellent agreement with various ARPES and scanning tunneling spectroscopy measurements. Despite the universal origin of the e-ph coupling, the contributions of several characteristic phonon vibrations vary in different phenomena, highlighting a versatile role of e-ph coupling in shaping the low-energy excitations of kagome metals.

Electron-phonon coupling↗

DOE OSTI2022

Constraints on the Maximum Densities of Neutron Stars from Postmerger Gravitational Waves with Third-Generation Observations

Using data from 289 numerical relativity simulations of binary neutron star mergers, we identify, for the first time, a robust quasiuniversal relation connecting the postmerger peak gravitational-wave frequency and the value of the density at the center of the maximum mass nonrotating neutron star. This relation offers a new possibility for precision equation-of-state constraints with next-generation ground-based gravitational-wave interferometers. Mock Einstein Telescope observations of fiducial events indicate that Bayesian inferences can constrain the maximum density to ~15% (90% credibility level) for a single signal at the minimum sensitivity threshold for a detection. If the postmerger signal is included in a full-spectrum (inspiral-merger-postmerger) analysis of such a signal, the pressure-density function can be tightly constrained up to the maximum density, and the maximum neutron star mass can be measured with an accuracy better than 12% (90% credibility level).

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Parametrically driven nonlinear Dirac equation with arbitrary nonlinearity

Here, the damped and parametrically driven nonlinear Dirac equation with arbitrary nonlinearity parameter $\kappa$ is analyzed, when the external force is periodic in space and given by $\mathcal{f(x) = r}$ cos(Kx), both numerically and in a variational approximation using five collective coordinates (time dependent shape parameters of the wave function). Our variational approximation satisfies exactly the low-order moment equations. Because of competition between the spatial period of the external force λ = 2π/K , and the soliton width $l_s$, which is a function of the nonlinearity $\kappa$ as well as the initial frequency ω 0 of the solitary wave, there is a transition (at fixed ω 0 ) from trapped to unbound behavior of the soliton, which depends on the parameters $\mathcal{r}$ and K of the external force and the nonlinearity parameter $\kappa$. We previously studied this phenomena when $\kappa$ = 1 where we showed that for λ $\gg l_s$ the soliton oscillates in an effective potential, while forλ $\ll l_s$ it moves uniformly as a free particle. In this paper we focus on the $\kappa$ dependence of the transition from oscillatory to particle behavior and explicitly compare the curves of the transition regime found in the collective coordinate approximation as a function of $\mathcal{r}$ and K when $\kappa$ = 1/2,1,2 at fixed value of the frequency ω 0 . Since the solitary wave gets narrower for fixed ω 0 as a function of $\kappa$, we expect and indeed find that the regime where the solitary wave is trapped is extended as we increase $\kappa$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Shock Hugoniot of Forged and Additively Manufactured 304L Stainless Steel

The purpose of this research was to measure the equation of state for additively manufactured (AM) and forged 304L stainless steel using a novel experimental technique. An understanding of the dynamic behavior of AM metals is integral to their timely adoption into various applications. The Hugoniot of the AM 304L was compared to that of the forged 304L at particle velocities where the material retains a two-wave structure. This comparison enabled us to determine the sensitivity of the equation of state to microstructure as varied due to processing. Our results showed that there was a measurable difference in the measured shock velocity between the AM and forged 304L. The shock wave velocities for the AM 304L were found to be ~3% slower than those for the forged 304L at similar particle velocities. To understand these differences, properties such as densities, sound speeds, and texture were measured and compared between the forged and AM materials. Our results showed that no measurable difference was found in these properties. Additionally, it is possible that differing elastic wave amplitudes may influence shock velocity.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Variational Monte Carlo Calculations of A ≤ 4 Nuclei with an Artificial Neural-Network Correlator Ansatz

Here, the complexity of many-body quantum wave functions is a central aspect of several fields of physics and chemistry where nonperturbative interactions are prominent. Artificial neural networks (ANNs) have proven to be a flexible tool to approximate quantum many-body states in condensed matter and chemistry problems. In this work we introduce a neural-network quantum state ansatz to model the ground-state wave function of light nuclei, and approximately solve the nuclear many-body Schrodinger equation. Using efficient stochastic sampling and optimization schemes, our approach extends pioneering applications of ANNs in the field, which present exponentially scaling algorithmic complexity. We compute the binding energies and point-nucleon densities of A ≤ 4 nuclei as emerging from a leading-order pionless effective field theory Hamiltonian. We successfully benchmark the ANN wave function against more conventional parametrizations based on two- and three-body Jastrow functions, and virtually exact Green's function Monte Carlo results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Adaptive Variational Quantum Imaginary Time Evolution Approach for Ground State Preparation

Abstract An adaptive variational quantum imaginary time evolution (AVQITE) approach is introduced that yields efficient representations of ground states for interacting Hamiltonians on near‐term quantum computers. It is based on McLachlan's variational principle applied to imaginary time evolution of variational wave functions. The variational parameters evolve deterministically according to equations of motions that minimize the difference to the exact imaginary time evolution, which is quantified by the McLachlan distance. Rather than working with a fixed variational ansatz, where the McLachlan distance is constrained by the quality of the ansatz, the AVQITE method iteratively expands the ansatz along the dynamical path to keep the McLachlan distance below a chosen threshold. This ensures the state is able to follow the quantum imaginary time evolution path in the system Hilbert space rather than in a restricted variational manifold set by a predefined fixed ansatz. AVQITE is used to prepare ground states of H 4 , H 2 O, and BeH 2 molecules, where it yields compact variational ansätze and ground state energies within chemical accuracy. Polynomial scaling of circuit depth with system size is shown through a set of AVQITE calculations of quantum spin models. Finally, quantum Lanczos calculations are demonstrated alongside AVQITE without additional quantum resource costs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Incoherent diffractive dijet production and gluon Bose enhancement in the nuclear wave function

We investigate the effect of gluon Bose enhancement in the nuclear wave function on the dijet production in incoherent diffractive processes in DIS and ultraperipheral collisions. We demonstrate that Bose enhancement leads to an enhancement of diffractive dijet production cross section when the transverse momenta of the two jets are aligned at zero relative angle. This enhancement is maximal when the magnitude of the transverse momenta of the two jets are equal, and disappears rather quickly as a function of the ratio of the two momenta. We study both the dilute limit and fully nonlinear dense regime where the nuclear wave function is evolved with the leading order JIMWLK equation. In both cases we observe a visible effect, with it being enhanced by the evolution due to the dynamical generation of the color neutralization scale.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

SchrödingerNet: A Universal Neural Network Solver for the Schrödinger Equation

Recent advances in machine learning have facilitated numerically accurate solution of the electronic Schrödinger equation (SE) by integrating various neural network (NN)-based wave function ansatzes with variational Monte Carlo methods. Nevertheless, such NN-based methods are all based on the Born–Oppenheimer approximation (BOA) and require computationally expensive training for each nuclear configuration. In this work, we propose a novel NN architecture, SchrödingerNet, to solve the full electronic-nuclear SE by defining a loss function designed to equalize local energies across the system. This approach is based on a translationally, rotationally and permutationally symmetry-adapted total wave function ansatz that includes both nuclear and electronic coordinates. Furthermore, this strategy not only allows for an efficient and accurate generation of a continuous potential energy surface at any geometry within the well-sampled nuclear configuration space, but also incorporates non-BOA corrections, through a single training process. Comparison with benchmarks of atomic and small molecular systems demonstrates its accuracy and efficiency.

Chemical calculations↗

DOE OSTI2022

Coupled electronic and magnetic excitations in the cuprates and their role in the superconducting transition

The formation of Cooper pairs, a bound state of two electrons of opposite spin and momenta by exchange of a phonon, is a defining feature of conventional superconductivity. In the cuprate high temperature superconductors, even though the superconducting state also consists of Cooper pairs, the pairing mechanism remains intensely debated. Here, we investigate superconducting pairing in the Bi 2 Sr 2 CaCu 2 O 8+δ (Bi2212) cuprate by employing spectral functions obtained from angle-resolved photoemission as input to the Bethe-Salpeter equation. Assuming Cooper pairing is driven by spin fluctuations, we construct the spin-fluctuation-mediated pairing interaction and use it to compute the eigenfunctions and eigenvalues of the Bethe-Salpeter equation for multiple Bi2212 samples. The leading d-wave eigenvalue increases as the temperature is decreased toward $T_c$, reaching a value of approximately 1 at the $T_c$ corresponding to each doping value. This suggests that spin fluctuations can approximately account for $T_c$ and mediate pairing in the cuprate superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Distinguishing prompt-collapse binary neutron star mergers from binary black Holes: Tidal effects and remnant properties

We study the properties of remnants formed in prompt-collapse binary neutron star mergers. We consider nonspinning neutron star binaries over a range of total masses and mass ratios across a set of 22 equations of state, totaling 107 numerical relativity simulations. We report the final mass and spin of the systems (including the accretion disk and ejecta) to be constrained in a narrow range—0.98 ≲ 𝑀 𝑓 /𝑀 ≲ 0.99 for the mass and 0.85 ≲ 𝑎 𝑓 ≲ 0.95 for the dimensionless spin—regardless of the binary configuration and matter effects. This sets them apart from binary black hole merger remnants. We assess the detectability of the postmerger signal in a future 40 km Cosmic Explorer observatory and find that the signal-to-noise ratio in the postmerger of an optimally located and oriented binary at a distance of 100 Mpc can range from <1 to 8, depending on the binary configuration and equation of state, with a majority of them greater than 4 in the set of simulations that we consider. We also consider the distinguishability between prompt-collapse binary neutron star and binary black hole mergers with the same masses and spins. We find that Cosmic Explorer will be able to distinguish such systems primarily via the measurement of tidal effects in the late inspiral. Neutron star binaries with reduced tidal deformability $\tilde{Λ}$ as small as ∼ 3.5 can be identified up to a distance of 100 Mpc, while neutron star binaries with $\tilde{Λ}$ ∼ 22 can be identified to distances greater than 250 Mpc. This is larger than the distance up to which the postmerger will be visible. Finally, we discuss the possible implications of our findings for the equation of state of neutron stars from the gravitational wave event GW230529.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Landau-level composition of bound exciton states in magnetic field

Here, we present a theory that studies the state composition of a bound exciton in magnetic field. Using a basis set made of products of free electron and hole wave functions in Landau gauge, we derive a secular equation which shows the relation between Landau levels (LLs) of the electron and hole when a bound exciton is formed. Focusing on excitons in the light cone, we establish a scattering selection rule for the interaction of an electron in LL 𝑛 e with a hole in LL 𝑛 h . We solve the resulting secular equation and identify a simple pairing law, 𝑛 e =𝑛 h +𝑙, which informs us on the construction of a bound exciton state with magnetic quantum number 𝑙, and on the interaction of the exciton magnetic moment with magnetic field. We obtain good agreement between theory results and recent measurements of the diamagnetic shifts of exciton states in WSe 2 monolayers.

2-dimensional systems↗