Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Wave Function”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

code for "Improving the QBO forcing by resolved waves with vertical grid refinement in E3SMv2"

This is all code that used in the journal article "Improving the QBO forcing by resolved waves with vertical grid refinement in E3SMv2", including: Running and post processing E3SM: Run E3SM: run_E3SM.2023.scidac.MMF.amip.py Regrid: regrid.E3SM.py post processing, diagnostics, archive: run_post.2023.scidac.MMF.py TEM calculation, including filter TEM, E3SM: calculate.TEM.v2.py calculate.TEM.monthly_mean.v1.py TEM, ERA5: ERA5_calculate_TEM_90.py merge TEM data: merge_TEM_E3SM_ERA5.py TEM, filter, E3SM: filter_TEM_E3SM.py TEM, filter, ERA5: filter_TEM_ERA5.py WK wave analysis code from (https://github.com/brianpm/wavenumber_frequency): wavenumber_frequency_functions.py calculate the wave spectra so easy to plot: wk_calculation.py Make figures some functions for plot: plot.py code for drawing paper plots: figures_for_paper.ipynb

Hannah, WalterM↗
DOE OSTI2021

Multi-Scenario Extreme Weather Simulator

The Multi-scenario Extreme Weather Simulator (MEWS) takes existing historical weather data and inserts trends and extreme events such as heat waves, hurricanes, extreme precipitation, or winter storms. The code includes functionality for projecting increasing frequency and duration of such events over a large range of statistically sampled realizations of the future. These weather files can then be used in building energy simulations through the Energy Plus software created by DOE. 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 SAND2021-8658 O

Villa, Daniel↗
DOE OSTI2021

Ab initio electronic density in solids by many-body plane-wave auxiliary-field quantum Monte Carlo calculations

We present accurate many-body results of the electronic densities in several solid materials, including Si, NaCl, and Cu. These results are obtained using the ab initio auxiliary-field quantum Monte Carlo (AFQMC) method working in a plane-wave basis with norm-conserving, multiple-projector pseudopotentials. AFQMC has been shown to be an excellent many-body total energy method. Computation of observables and correlation functions other than the ground-state energy requires back-propagation, whose adaption and implementation in the plane-wave basis AFQMC framework are discussed in the present paper. This development allows us to compute correlation functions, electronic densities, and interatomic forces, paving the way for geometry optimizations and calculations of thermodynamic properties in solids. Finite supercell size effects are considerably more subtle in the many-body framework than in independent-electron calculations. We analyze the convergence of the electronic density, and obtain best estimates for the thermodynamic limit. The densities from several typical density functionals are benchmarked against our near-exact results. The electronic densities we obtained can also be used to help construct improved density functionals.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Unified wavefront singularity characterization of three-dimensional elastodynamic time-domain half-space Green's function under impulsive boundary and internal loads

Founded on a novel analytical formulation that led to a rigorous yet compact path-integral representation of the time-domain elastodynamic half-space Green's function, a unified analysis of the possible occurrence of different singular wavefront behaviour in the response under arbitrary impulsive internal or surface point loads at arbitrary source-receiver locations is presented. With the decomposition of the general solution into distinct initiating and reflected wave group integrals that share a common factored format and simple contour definitions, the mathematical framework is shown to allow a straightforward identification of the specific conditions and the particular wave groups that are responsible for the singular wavefront phenomena without resorting to advanced analytic function theories or asymptotic methods. Analytic characterizations of the nature, strength and direction of all intrinsic singular wavefront behaviours of the three-dimensional Green's function in three canonical cases of source-receiver configurations are given in a dual integral-closed form format to facilitate their theoretical understanding as well as computational applications. Graphical illustrations of their variation with the source-receiver configuration and the medium's Poisson's ratio together with relevant comparison and clarifications of some classical treatments are included.

Green’s functions↗
DOE OSTI2022

A Realistic Theory of Quantum Measurement

Abstract We propose that the ontic understanding of quantum mechanics can be extended to a fully realistic theory that describes the evolution of the wavefunction at all times, including during a measurement. In such an approach the wave equation should reduce to the standard wave equation when there is no measurement, and describe state reduction when the system is measured. The general wave equation must be nonlinear and nonlocal, and we require it to be time-symmetric; consequently, this approach is not a new interpretation but a new theory. The wave equation is an integrodifferential equation (IDE). The time symmetry requirement leads to a retrocausal approach, in which the wave equation is solved subject to initial and final conditions to determine history at intermediate times. We propose that different outcomes from (apparently) identically prepared experiments may result from uncontrolled parameters; both the nonlocality and the retrocausality of the theory imply that Bell’s Theorem cannot rule out such “hidden variables.” Beginning with Hamilton’s principle, we demonstrate the construction of such a theory by replacing the action with a functional designed to give rise to a nonlinear, nonlocal IDE as the wave equation. This IDE reduces to the standard wave equation (a differential equation) in the absence of a measurement, but exhibits state reduction to a single eigenvalue when the system interacts with another system with the properties of a measurement apparatus. We demonstrate several desirable features of this theory; for other properties we indicate their plausibility and possible avenues to a proof.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Nonlinear saturation of whistler modes driven by runaway electrons

Runaway electrons exhibit kinetic instabilities with potentially beneficial consequences. The anomalous Doppler resonance between the electrons and whistler modes is a primary underlying mechanism. These instabilities require a first-principle nonlinear theoretical analysis, which would ultimately enable assessment of their impact on disruption mitigation and diagnostics, especially in ITER-relevant conditions. This paper presents recent progress in developing the required theoretical framework for isolated nonlinear resonances. By employing the generic bump-on-tail model, we predict the wave saturation levels in both near-threshold and strongly driven regimes. Remarkably, for the waves of interest, the parallel component of the wave vector dictates the parallel momentum of the resonant particles. This feature, together with the expected wave saturation level, provides a complete description of the resonance impact on the runaway electron distribution function. We show that the parametric dependence of the nonlinearly saturated mode provides a way to recover certain features of the runaway momentum distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The influence of physical and algorithmic factors on simulated far-field waveforms and source–time functions of underground explosions using unsupervised machine learning

SUMMARY Characterizing explosion sources and differentiating between earthquake and underground explosions using distributed seismic networks becomes non-trivial when explosions are detonated in cavities or heterogeneous ground material. Moreover, there is little understanding of how changes in subsurface physical properties affect the far-field waveforms we record and use to infer information about the source. Simulations of underground explosions and the resultant ground motions can be a powerful tool to systematically explore how different subsurface properties affect far-field waveform features, but there are added variables that arise from how we choose to model the explosions that can confound interpretation. To assess how both subsurface properties and algorithmic choices affect the seismic wavefield and the estimated source functions, we ran a series of 2-D axisymmetric non-linear numerical explosion experiments and wave propagation simulations that explore a wide array of parameters. We then inverted the synthetic far-field waveform data using a linear inversion scheme to estimate source–time functions (STFs) for each simulation case. We applied principal component analysis (PCA), an unsupervised machine learning method, to both the far-field waveforms and STFs to identify the most important factors that control variance in the waveform data and differences between cases. For the far-field waveforms, the largest variance occurs in the shallower radial receiver channels in the 0–50 Hz frequency band. For the STFs, both peak amplitude and rise times across different frequencies contribute to the variance. We find that the ground equation of state (i.e. lithology and rheology) and the explosion emplacement conditions (i.e. tamped versus cavity) have the greatest effect on the variance of the far-field waveforms and STFs, with the ground yield strength and fracture pressure being secondary factors. Differences in the PCA results between the far-field waveforms and STFs could possibly be due to near-field non-linearities of the source that are not accounted for in the estimation of STFs and could be associated with yield strength, fracture pressure, cavity radius and cavity shape parameters. Other algorithmic parameters are found to be less important and cause less variance in both the far-field waveforms and STFs, meaning algorithmic choices in how we model explosions are less important, which is encouraging for the further use of explosion simulations to study how physical Earth properties affect seismic waveform features and estimated STFs.

58 GEOSCIENCES↗
DOE OSTI2021

Detecting the stochastic gravitational wave background from massive gravity with pulsar timing arrays

We explore the potential of pulsar timing arrays (PTAs) such as NANOGrav, EPTA, and PPTA to detect the stochastic gravitational wave background in theories of massive gravity. In general relativity, the function describing the dependence of the correlation between the arrival times of signals from two pulsars on the angle between them is known as the Hellings-Downs curve. We compute the analogous overlap reduction function for massive gravity, including the additional polarization states and the correction due to the mass of the graviton, and compare the result with the Hellings-Downs curve. The primary result is a complete analytical form for the analog Hellings-Downs curve, providing a starting point for future numerical studies aimed at a detailed comparison between PTA data and the predictions of massive gravity. We study both the massless limit and the stationary limit as checks on our calculation, and discuss how our formalism also allows us to study the impact of massive spin-2 dark matter candidates on data from PTAs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

First-principles study of the low-temperature charge density wave phase in the quasi-one-dimensional Weyl chiral compound ( TaSe 4 ) 2 I

Using ab initio density functional theory, we study herein the lattice phase transition of quasi-one-dimensional (TaSe 4 ) 2 I. In the undistorted state, the strongly anisotropic semimetal band structure presents two nonequivalent Weyl points. In previous efforts, two possible Ta-tetramerization patterns were proposed to be associated with the low-temperature structure. Our phonon calculations indicate that the orthorhombic F222 CDW-I phase is the most likely ground state for this quasi-one-dimensional system. In addition, the monoclinic C2 CDW-II phase may also be stable according to the phonon dispersion spectrum. Since these two phases have very similar energies in our density functional theory calculations, both these Ta-tetramerization distortions likely compete or coexist at low temperatures. The semimetal-to-insulator transition is induced by a Fermi-surface-driven instability that supports the Peierls scenario, which affects the Weyl physics developed above T CDW . Furthermore, the spin-orbit coupling generates Rashba-like band splittings in the insulating charge density wave phases.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Tracking Seismic Velocity Perturbations at Ridgecrest Using Ballistic Correlation Functions

Abstract We present results based on data of a dense nodal array composed of 147 stations, deployed in 2022 near the epicenter of the 2019 Mw 7.1 Ridgecrest earthquake to investigate characteristics of the seismic wavefields. Through array analyses, we identified two primary components. First, we observed far-field P waves dominating the 0.5–1.2 Hz frequency range, which are likely primarily generated by wind-driven oceanic swell activity. Second, we detected near-field body waves resulting from anthropogenic activities in the frequency range 2–8 Hz. We examined noise correlation functions derived from data of the dense deployment and regional stations to explore fault-zone seismic velocity changes using ballistic arrivals, with a focus on velocity perturbation shortly before and after the Ridgecrest earthquake sequence. Our findings exhibit distinct behavior compared to results obtained through standard coda-wave interferometry. Particularly, we observed a decrease in P-wave travel time on certain station pairs prior to the 2019 earthquake sequence. Supported by detailed investigation of the local seismic wavefields, we interpret the decreasing P-wave travel time as likely caused by a velocity increase away from the fault, possibly related to fluid migration. However, additional information is necessary to verify this hypothesis.

Geochemistry & Geophysics↗
DOE OSTI2021

Microscopic mechanism of unusual lattice thermal transport in TlInTe 2

We investigate the microscopic mechanism of ultralow lattice thermal conductivity (κ l ) of TlInTe 2 and its weak temperature dependence using a unified theory of lattice heat transport, that considers contributions arising from the particle-like propagation as well as wave-like tunneling of phonons. While we use the Peierls–Boltzmann transport equation (PBTE) to calculate the particle-like contributions (κ l (PBTE)), we explicitly calculate the off-diagonal (OD) components of the heat-flux operator within a first-principles density functional theory framework to determine the contributions (κ l (OD)) arising from the wave-like tunneling of phonons. At each temperature, T, we anharmonically renormalize the phonon frequencies using the self-consistent phonon theory including quartic anharmonicity, and utilize them to calculate κ l (PBTE) and κ l (OD). With the combined inclusion of κ l (PBTE), κ l (OD), and additional grain-boundary scatterings, our calculations successfully reproduce the experimental results. Our analysis shows that large quartic anharmonicity of TlInTe 2 (a) strongly hardens the low-energy phonon branches, (b) diminishes the three-phonon scattering processes at finite T, and (c) recovers the weaker than T –1 decay of the measured κ l .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Wave Dark Matter

Here we review the physics and phenomenology of wave dark matter: a bosonic dark matter candidate lighter than about 30 eV. Such particles have a de Broglie wavelength exceeding the average interparticle separation in a galaxy like the Milky Way and are, thus, well described as a set of classical waves. We outline the particle physics motivations for such particles, including the quantum chromodynamics axion as well as ultralight axion-like particles such as fuzzy dark matter. The wave nature of the dark matter implies a rich phenomenology: Wave interference gives rise to order unity density fluctuations on de Broglie scale in halos. One manifestation is vortices where the density vanishes and around which the velocity circulates. There is one vortex ring per de Broglie volume on average. For sufficiently low masses, soliton condensation occurs at centers of halos. The soliton oscillates and undergoes random walks, which is another manifestation of wave interference. The halo and subhalo abundance is expected to be suppressed at small masses, but the precise prediction from numerical wave simulations remains to be determined. For ultralight ∼10 −22 eV dark matter, the wave interference substructures can be probed by tidal streams or gravitational lensing. The signal can be distinguished from that due to subhalos by the dependence on stream orbital radius or image separation. Axion detection experiments are sensitive to interference substructures for wave dark matter that is moderately light. The stochastic nature of the waves affects the interpretation of experimental constraints and motivates the measurement of correlation functions. Current constraints and open questions, covering detection experiments and cosmological, galactic, and black hole observations, are discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

An efficient method to estimate the probability density of seismic Green's functions

We present a computationally efficient method to approximate the probability distribution of seismic Green's functions given the uncertainty of an Earth model. The method is based on the Karhunen-Loève (KL) theorem and an approximation of the Green's function (or seismogram) covariance. Using Monte Carlo (MC) simulations as a control case, we demonstrate that our KL-based method can accurately reproduce a probability distribution of seismograms that results from an uncertain Earth model for a MC-derived seismogram covariance. We then describe a method to estimate the covariance of the seismograms resulting from those Earth models that is not based on MC simulations. We use the estimated Green's function covariance in conjunction with our KL-based method to produce a Green's function probability distribution, and compare that distribution to a Green's function probability distribution produced using a MC finite difference method. We find that the Green's function probability distribution approximated using our KL-based method generally mimics that produced using the MC simulations, especially for direct-arriving body waves. However the accuracy of the KL-based method generally decreases for later times in the simulated Green's function distribution.

58 GEOSCIENCES↗
DOE OSTI2022

Reconfigurable training and reservoir computing in an artificial spin-vortex ice via spin-wave fingerprinting

Strongly interacting artificial spin systems are moving beyond mimicking naturally occurring materials to emerge as versatile functional platforms, from reconfigurable magnonics to neuromorphic computing. Typically, artificial spin systems comprise nanomagnets with a single magnetization texture: collinear macrospins or chiral vortices. Here, by tuning nanoarray dimensions we have achieved macrospin–vortex bistability and demonstrated a four-state metamaterial spin system, the ‘artificial spin-vortex ice’ (ASVI). ASVI can host Ising-like macrospins with strong ice-like vertex interactions and weakly coupled vortices with low stray dipolar field. Vortices and macrospins exhibit starkly differing spin-wave spectra with analogue mode amplitude control and mode frequency shifts of Δf = 3.8 GHz. The enhanced bitextural microstate space gives rise to emergent physical memory phenomena, with ratchet-like vortex injection and history-dependent non-linear fading memory when driven through global magnetic field cycles. We employed spin-wave microstate fingerprinting for rapid, scalable readout of vortex and macrospin populations, and leveraged this for spin-wave reservoir computation. ASVI performs non-linear mapping transformations of diverse input and target signals in addition to chaotic time-series forecasting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Self-sustaining sound in collisionless, high- β plasma

Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $δn/n_0≳2/β$ (where $β\gg1$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves’ otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more ‘fluid-like’ than kinetic. A nonlinear fluctuation–dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels and weakly collisional thermodynamics when IAWs with $δn/n_0≳2/β$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Influence of Polymorphs and Local Defect Structures on NMR Parameters of Graphite Fluorides

In this study, the role of local molecular structure on calculated 13 C and 19 F NMR chemical shifts for graphite fluoride materials was explored by using gauge-including projector augmented wave (GIPAW) computational methods for different periodic crystal polymorphs and density functional theory (DFT) gauge-including atomic orbital (GIAO) computational methods for individual graphite fluoride platelets, i.e., fluorinated graphene (FG). The impact of stacking sequences, d -spacing, and ring conformations on fully fluorinated graphite fluoride structures was investigated. A range of different defects including Stone–Wales, F and C vacancies, void formation, and F inversion were also evaluated using FG structures. These calculations show that distinct chemical shift signatures exist for many of these polymorphs and defects, therefore providing a basis for spectral assignment and development of models describing the mean local CF structure in disordered graphite fluoride materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Tunable magnons of an antiferromagnetic Mott insulator via interfacial metal-insulator transitions

Antiferromagnetic insulators present a promising alternative to ferromagnets due to their ultrafast spin dynamics essential for low-energy terahertz spintronic device applications. Magnons, i.e., quantized spin waves capable of transmitting information through excitations, serve as a key functional element in this paradigm. However, identifying external mechanisms to effectively tune magnon properties has remained a major challenge. Here we demonstrate that interfacial metal-insulator transitions offer an effective method for controlling the magnons of Sr 2 IrO 4 , a strongly spin-orbit coupled antiferromagnetic Mott insulator. Resonant inelastic x-ray scattering experiments reveal a significant softening of zone-boundary magnon energies in Sr 2 IrO 4 films epitaxially interfaced with metallic 4 d transition-metal oxides. Therefore, the magnon dispersion of Sr 2 IrO 4 can be tuned by metal-insulator transitions of the 4 d transition-metal oxides. We tentatively attribute this non-trivial behavior to a long-range phenomenon mediated by magnon-acoustic phonon interactions. Our experimental findings introduce a strategy for controlling magnons and underscore the need for further theoretical studies to better understand the underlying microscopic interactions between magnons and phonons.

magnetic properties and materials↗