Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave prediction”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Discovering new forces with gravitational waves from supermassive black holes

Supermassive black hole binary mergers generate a stochastic gravitational wave background detectable by pulsar timing arrays. While the amplitude of this background is subject to significant uncertainties, the frequency dependence is a robust prediction of general relativity. Here we show that the effects of new forces beyond the Standard Model can modify this prediction and introduce unique features into the spectral shape. In particular, we consider the possibility that black holes in binaries are charged under a new long-range force, and we find that pulsar timing arrays are capable of robustly detecting such forces. Supermassive black holes and their environments can acquire charge due to high-energy particle production or dark sector interactions, making the measurement of the spectral shape a powerful test of fundamental physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Observations and seismoacoustic simulations of earthquake-generated infrasound waves in non-epicentral regions

We analysed infrasound waves associated with the Gyeongju earthquake (ML 5.8) that occurred on 2016 September 12, in the southeastern Korean Peninsula. For infrasound wave detection, the Progressive Multichannel Correlation method was applied to the infrasound data set recorded at seven arrays operating in South Korea at epicentral distances ranging from 178 to 472 km. Based on the back-projection method constrained by array-dependent celerity and azimuth deviation models, the source regions were identified in both the epicentral and non-epicentral regions. Remarkably, the non-epicentral secondary sources of this earthquake were located in regions with shallow water depths: (i) the western coastal area in the Yellow Sea and (ii) the shallow ocean basin and bank in the East Sea. The location results obtained from the earthquake could be corroborated through its foreshock (ML 5.1), yielding location results consistent with those of the main shock. The generation of infrasound waves over shallow water depths was fortuitously validated by direct recordings of dominant single-frequency (~0.3 Hz) infrasound waves at close range via temporary sensors near the ocean basin and bank. We interpreted that low-frequency infrasound signals could be generated from interactions among the ocean floor, shallow seawater and atmosphere. We performed numerical simulations of seismoacoustic fields to predict ground motions on the seafloor and acoustic transmission efficiency between the water and air interface. The simulations quantified the energy transfer through different media and clarified our observational results. We found that because this solid Earth‒water‒atmosphere coupled air wave has a relatively low frequency (~0.3 Hz), it can survive propagation over long distances compared with high-frequency infrasound waves generated in inland and mountain regions. In this study, we extend our understanding of water‒atmosphere coupling and the monitoring framework for earthquake-associated non-epicentral infrasound waves, encompassing not only inland ground shaking but also shallow sea regions located far from the epicentre.

Earthquake monitoring and test-ban treaty verifica↗
DOE OSTI2021

A generalized machine learning workflow to visualize mechanical discontinuity

Accurate detection and mapping of mechanical discontinuity in materials has widespread industrial and research applications. Herein, we developed a generalized machine-learning framework for visualizing single mechanical discontinuity embedded in material of any composition, velocity, density, porosity, and size with limited data. The proposed visualization of discontinuity requires accurate estimations of the length, location, and orientation of the embedded discontinuity by processing multipoint wave-transmission measurements. k-Wave simulator is used to create a large dataset of elastic waveforms recorded during multi-point wave-transmission measurements through materials containing single mechanical discontinuity. k-Wave simulator considers the wave attenuation, dispersion, and mode conversion in wave motion. Discrete wavelet transform (DWT) and statistical feature extraction are essential for data preprocessing prior to the data-driven model development. DWT also minimizes the effect of noise. Using hyper-parameter tuning and cross validation, gradient boosting regression can visualize the mechanical discontinuity with an accuracy of 0.85, in terms of coefficient of determination. A double-layered neural network-based regression has better performance with an accuracy of 0.95. Use of convolutional neural network converts the predictive task from a waveform processing to an image processing problem. Convolutional neural network achieved a generalization performance of 0.91. The proposed generalized workflow requires robust simulation of wave propagation, signal processing, feature engineering, and model evaluation. Sensors closest to the source and those located opposite the source are the most significant for the desired visualization. Notably, the sensors closest to the source capture the non-linear associations, whereas the sensor on the border opposite to the source capture the linear associations between the measured waveforms and the properties of the mechanical discontinuity.

42 ENGINEERING↗
DOE OSTI2024

Finite-frequency modeling of regional tropospheric infrasound using realistic atmospheres and terrain

Infrasonic waves have been observed to propagate to regional (greater than 15 km) distances through the troposphere. Infrasound propagation in the geometric acoustics approximation has shown that realistic terrain can scatter acoustic energy from tropospheric ducts; however, ray methods cannot intrinsically capture finite-frequency behavior such as diffraction. A two-dimensional finite-difference time-domain (FDTD) method has been developed to solve linearized equations for infrasound propagation with realistic terrain. Acoustic wave propagation over 100 km with both flat terrain and a Gaussian hill was first simulated in order to compare finite-frequency propagation with ray predictions. The effects of realistic terrain and atmospheres on infrasound signals from a 2012 surface explosion at the Utah Testing and Training Range are then investigated. Propagation through the troposphere is suggested by array processing results, but eigenrays are not predicted due to weak to nonexistent ducting conditions. FDTD modeling suggests that the inclusion of terrain and finite frequency effects helps explain much of the observed signal in a realistic scenario. Furthermore, these results suggest that geometric acoustics may underestimate propagation through the troposphere, and that recorded waveforms at regional distances may be noticeably affected by terrain.

58 GEOSCIENCES↗
DOE OSTI2022

A new woven composite constitutive model validated by shock wave experiments

In this paper, we present results of plate impact simulations of shock compressed woven glass fiber-reinforced plastic (GRP) performed using the Arbitrary Lagrangian–Eulerian three-dimensional finite element code. A hyperelastic large-strain-based empirical Continuum Damage Mechanics (CDM) formulation is employed to describe damage initiation and growth in the shock-compressed GRP. The model parameters calibration scheme utilizes the Velocity Interferometer System for Any Reflector normal particle velocity measurements at the free surface of the GRP target plates. The impact velocity in the experiments ranged from 8.5 to 418 m/s. The finite element model considered planar 0°/90° bidirectional plies with an individual ply thickness of 0.68 mm, stacked to reach a total laminate thickness of 6.8 mm. The anisotropic elastic strains were estimated from the experimentally determined tetragonal symmetry stiffness matrix for the GRP. The strain-based damage model captures several salient features observed in the measured free surface particle wave profiles, including the shock rise time, onset of Elastic—Elastic Cracking, and the shape of the nonlinear portion of the experimental particle velocity profiles. Furthermore, the CDM model predicts the dominant damage mode to be matrix microcracking due to shear and the associated bulk expansion (bulking) under the global compressive loading in the plate impact configuration.

42 ENGINEERING↗
DOE OSTI2021

Focused wave interactions with floating structures: a blind comparative study

The paper presents results from the Collaborative Computational Project in Wave Structure Interaction (CCP-WSI) Blind Test Series 2. Without prior access to the physical data, participants, with numerical methods ranging from low-fidelity linear models to fully non-linear Navier-Stokes (NS) solvers, simulate the interaction between focused wave events and two separate, taut-moored, floating structures: a hemispherical-bottomed cylinder and a cylinder with a moonpool. The ‘blind’ numerical predictions for heave, surge, pitch and mooring load, are compared against physical measurements. Dynamic time warping is used to quantify the predictive capability of participating methods. In general, NS solvers and hybrid methods give more accurate predictions; however, heave amplitude is predicted reasonably well by all methods; and a WEC-Sim implementation, with CFD-informed viscous terms, demonstrates comparable predictive capability to even the stronger NS solvers. Large variations in the solutions are observed (even among similar methods), highlighting a need for standardisation in the numerical modelling of WSI problems.

42 ENGINEERING↗
DOE OSTI2024

Data Imbalance, Uncertainty Quantification, and Transfer Learning in Data‐Driven Parameterizations: Lessons From the Emulation of Gravity Wave Momentum Transport in WACCM

Abstract Neural networks (NNs) are increasingly used for data‐driven subgrid‐scale parameterizations in weather and climate models. While NNs are powerful tools for learning complex non‐linear relationships from data, there are several challenges in using them for parameterizations. Three of these challenges are (a) data imbalance related to learning rare, often large‐amplitude, samples; (b) uncertainty quantification (UQ) of the predictions to provide an accuracy indicator; and (c) generalization to other climates, for example, those with different radiative forcings. Here, we examine the performance of methods for addressing these challenges using NN‐based emulators of the Whole Atmosphere Community Climate Model (WACCM) physics‐based gravity wave (GW) parameterizations as a test case. WACCM has complex, state‐of‐the‐art parameterizations for orography‐, convection‐, and front‐driven GWs. Convection‐ and orography‐driven GWs have significant data imbalance due to the absence of convection or orography in most grid points. We address data imbalance using resampling and/or weighted loss functions, enabling the successful emulation of parameterizations for all three sources. We demonstrate that three UQ methods (Bayesian NNs, variational auto‐encoders, and dropouts) provide ensemble spreads that correspond to accuracy during testing, offering criteria for identifying when an NN gives inaccurate predictions. Finally, we show that the accuracy of these NNs decreases for a warmer climate (4 × CO 2 ). However, their performance is significantly improved by applying transfer learning, for example, re‐training only one layer using ∼1% new data from the warmer climate. The findings of this study offer insights for developing reliable and generalizable data‐driven parameterizations for various processes, including (but not limited to) GWs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Tailoring the Weight of Surface and Intralayer Edge States to Control LUMO Energies

Abstract The energies of the frontier molecular orbitals determine the optoelectronic properties in organic films, which are crucial for their application, and strongly depend on the morphology and supramolecular structure. The impact of the latter two properties on the electronic energy levels relies primarily on nearest‐neighbor interactions, which are difficult to study due to their nanoscale nature and heterogeneity. Here, an automated method is presented for fabricating thin films with a tailored ratio of surface to bulk sites and a controlled extension of domain edges, both of which are used to control nearest‐neighbor interactions. This method uses a Langmuir–Schaefer‐type rolling transfer of Langmuir layers (rtLL) to minimize flow during the deposition of rigid Langmuir layers composed of π‐conjugated molecules. Using UV–vis absorption spectroscopy, atomic force microscopy, and transmission electron microscopy, it is shown that the rtLL method advances the deposition of multi‐Langmuir layers and enables the production of films with defined morphology. The variation in nearest‐neighbor interactions is thus achieved and the resulting systematically tuned lowest unoccupied molecular orbital (LUMO) energies (determined via square‐wave voltammetry) enable the establishment of a model that functionally relates the LUMO energies to a morphological descriptor, allowing for the prediction of the range of accessible LUMO energies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Vaccination compartmental epidemiological models for the delta and omicron SARS-CoV-2 variants

We explore the inclusion of vaccination in compartmental epidemiological models concerning the delta and omicron variants of the SARS-CoV-2 virus that caused the COVID-19 pandemic. We expand on our earlier compartmental-model work by incorporating vaccinated populations. We present two classes of models that differ depending on the immunological properties of the variant. The first one is for the delta variant, where we do not follow the dynamics of the vaccinated individuals since infections of vaccinated individuals were rare. The second one for the far more contagious omicron variant incorporates the evolution of the infections within the vaccinated cohort. We explore comparisons with available data involving two possible classes of counts, fatalities and hospitalizations. We present our results for two regions, Andalusia and Switzerland (including the Principality of Liechtenstein), where the necessary data are available. In the majority of the considered cases, the models are found to yield good agreement with the data and have a reasonable predictive capability beyond their training window, rendering them potentially useful tools for the interpretation of the COVID-19 and further pandemic waves, and for the design of intervention strategies during these waves.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Exact restoration of Galilei invariance in density functional calculations with quantum Monte Carlo

We investigate some aspects of the connection between mean-field and ab-initio calculations for any-nucleon systems employing the Gogny effective interaction. In particular, we present variational Monte Carlo calculations of 4 He and 16 O, for which we extended the formalism of the linear method to the case of a density-dependent Hamiltonian. Monte Carlo calculations using an uncorrelated basis reduce to Hartree-Fock ones, but Galilean invariance can be preserved, contrary to the usual treatment. An analysis of the results obtained with the Gogny force when used in conjunction with correlated wave functions has also been performed. Here, the calculations show the collapse of the deuteron, leading to a set of unphysical predictions for the structure of heavier nuclei. This behavior remarks the different role of the Gogny effective interaction as compared to bare interactions typically used in fully microscopic calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Two methods to study inelastic neutron-scattering measurements based on ω n (q) versus S(q, ω) applied to the magnetic open honeycomb lattice Tb 2 Ir 3 Ga 9

This work describes two methods to fit the inelastic neutron-scattering spectrum S(q, ω) with wavevector q and frequency ω. The common and well-established method extracts the experimental spin-wave branches ω n (q) from the measured spectra S(q, ω) and then minimizes the difference between the observed and predicted frequencies. When n branches of frequencies are predicted but the measured frequencies overlap to produce only m < n branches, the weighted average of the predicted frequencies must be compared to the observed frequencies. A penalty is then exacted when the width of the predicted frequencies exceeds the width of the observed frequencies. The second method directly compares the measured and predicted intensities S(q, ω) over a grid {q i , ω j } in wavevector and frequency space. After subtracting background noise from the observed intensities, the theoretical intensities are scaled by a simple wavevector-dependent function that reflects the instrumental resolution. Furthermore, the advantages and disadvantages of each approach are demonstrated by studying the open honeycomb material Tb 2 Ir 3 Ga 9 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Millimeter-wave beam scattering and induced broadening by plasma turbulence in the TCV tokamak

The scattering of millimeter-wave beams from electron density fluctuations and the associated beam broadening are experimentally demonstrated. Using a dedicated setup, instantaneous deflection and (de-)focusing of the beam due to density blobs on the beam path are shown to agree with full-wave simulations. The detected time-averaged wave power transmitted through the turbulent plasma is reproduced by the radiative-transfer model implemented in the WKBeam code, which predicts a ~50% turbulence-induced broadening of the beam cross-section. Finally, the role of core turbulence for the considered geometry is highlighted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Chiral interactions and superfluidity in the calcium isotopic chain

We perform ab initio calculations of three-point mass differences in the odd- and even-mass 39−49Ca isotopes to probe nuclear superfluidity via empirical neutron pairing gaps. We also quantify the sensitivity of those gaps to the parameters of the interaction at mean-field level. Recent studies employing accurate chiral nuclear interactions have found these gaps to be too small. We show that experimental values can be reproduced at mean-field level by substantially increasing the attraction of the singlet 𝑆 -wave two-nucleon contact interaction, but doing so induces an unphysical bound state of the dineutron. The sensitivity of these predictions to the full calibration of the nuclear interaction is then studied by performing Bayesian posterior sampling in a delta-full chiral effective field theory at third chiral order. We find that pairing gaps remain largely unaffected, leaving the explanation of nuclear superfluidity as a future task for improved many-body modeling and refined interactions at higher chiral orders.

Hagen, Gaute [ORNL] (ORCID:0000000160191687)↗
DOE OSTI2020

Ab initio calculations of the rate of carrier trapping and release at dopant sites in NaI: Tl beyond the harmonic approximation

We present ab initio calculations of electron capture and release coefficients as a function of temperature for Tl dopants in the widely used scintillator NaI, which is a soft ionic crystal. The modelled capture and release events occur by transitions mediated by multiple phonon absorption and emission between states with significantly different local geometries around the trapping site. We demonstrate that such transitions are not well-described by the normal harmonic approximation to the nuclear dynamics. We avoid the harmonic approximation by numerically solving the vibrational Schrödinger equation for the motion along/of a single effective phonon coordinate. The localized trapped state is not correctly described by semilocal density functionals, so we employ a hybrid functional tuned to reproduce the band gap and lattice constant of the host material. Our calculations, which combine plane-wave and cluster electronic structure methods in a novel way, are in reasonable agreement with available experimental data for detrapping and predict an unusual temperature dependence for trapping rates.

Prange, Micah P.↗
DOE OSTI2024

Impact of dark matter on tidal signatures in neutron star mergers with the Einstein Telescope

If dark matter (DM) accumulates inside neutron stars (NS), it changes their internal structure and causes a shift of the tidal deformability from the value predicted by the dense-matter equation of state (EOS). In principle, this shift could be observable in the gravitational-wave (GW) signal of binary neutron star (BNS) mergers. We investigate the effect of fermionic, noninteracting DM when observing a large number of GW events from DM-admixed BNSs with the precision of the proposed Einstein telescope (ET). Specifically, we study the impact on the recovery of the baryonic EOS and whether DM properties can be constrained. For this purpose, we create event catalogs of BNS mock events with DM fraction up to 1%, from which we reconstruct the posterior uncertainties with the Fisher matrix approach. Using this data, we perform joint Bayesian inference on the baryonic EOS, DM particle mass, and DM particle fraction in each event. Here, our results reveal that when falsely ignoring DM effects, the EOS posterior is biased toward softer EOSs, though the offset is rather small. Further, we find that within our assumptions of our DM model and population, ET will likely not be able to test the presence of DM in BNSs, even when combining many events and adding Cosmic Explorer (CE) to the next-generation detector network. Likewise, the potential constraints on the DM particle mass will remain weak because of degeneracies with the fraction and EOS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Dynamical Generation of Higher-Order Spin-Orbit Coupling, Topology, and Persistent Spin Texture in Light-Irradiated Altermagnets

Altermagnets have been identified as the third category of magnetic materials, exhibiting momentum-dependent spin splitting characterized by even powers of momentum. Here, in this Letter, we show that when subjected to elliptically polarized light, these materials serve as an exemplary framework for the dynamic generation of topological bands featuring higher-order spin-orbit coupling (SOC). Notably, while the generated Zeeman field remains invariant to the particular altermagnetic ordering, the induced higher-order SOCs are related to the magnitude and symmetry of the altermagnetic order. Specifically, we show that an altermagnet exhibiting 𝐤 2⁢𝑛 -spin splitting can generate odd-in-𝐤 spin-orbit couplings up to 𝐤 2⁢𝑛−1 . In the limit of circularly polarized light, the only correction is 𝑘 2⁢𝑛−1 , with all lower-order contributions being nullified. Interestingly, light-induced SOCs significantly impact the low-energy band topology, where their Chern numbers change by Δ⁢𝐶 = ±1, 2, 3 for 𝑑, 𝑔, 𝑖-wave altermagnets. Finally, we find a critical field in which a persistent spin texture is realized, a highly desirable state with predicted infinite spin lifetime. Our Letter showcases light as a powerful, controllable tool for engineering complex and exciting phenomena in altermagnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Quantitative Correction of Ohmic Effects on Square Wave Voltammetry for High-Concentration Soluble-Soluble Redox Reactions in Molten Salts

Molten salts are of particular interest for a variety of industrial applications. As such, accurate characterization of species’ concentrations within the molten salt media is critical to ensure a well-controlled unit operation. Although electroanalytical tools are properly suited for precise in situ monitoring in these systems, uncompensated ohmic resistance R Ω can lead to erroneous, technique-dependent results. Using numerical simulations applied to a model system, this work will first illustrate then quantify the extent to which R Ω attenuates square wave voltammograms. This approach allows for post-experiment correction that leads to converging results across voltammetric techniques and facilitates accurate predictions of species’ concentrations.

Shaheen, Nora A. [Argonne National Laboratory (ANL↗