Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave data”

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,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Porting the WAVEWATCH III (v6.07) wave action source terms to GPU

Abstract. Surface gravity waves play a critical role in several processes, including mixing, coastal inundation, and surface fluxes. Despite the growing literature on the importance of ocean surface waves, wind–wave processes have traditionally been excluded from Earth system models (ESMs) due to the high computational costs of running spectral wave models. The development of the Next Generation Ocean Model for the DOE’s (Department of Energy) E3SM (Energy Exascale Earth System Model) Project partly focuses on the inclusion of a wave model, WAVEWATCH III (WW3), into E3SM. WW3, which was originally developed for operational wave forecasting, needs to be computationally less expensive before it can be integrated into ESMs. To accomplish this, we take advantage of heterogeneous architectures at DOE leadership computing facilities and the increasing computing power of general-purpose graphics processing units (GPUs). This paper identifies the wave action source terms, W3SRCEMD, as the most computationally intensive module in WW3 and then accelerates them via GPU. Our experiments on two computing platforms, Kodiak (P100 GPU and Intel(R) Xeon(R) central processing unit, CPU, E5-2695 v4) and Summit (V100 GPU and IBM POWER9 CPU) show respective average speedups of 2× and 4× when mapping one Message Passing Interface (MPI) per GPU. An average speedup of 1.4× was achieved using all 42 CPU cores and 6 GPUs on a Summit node (with 7 MPI ranks per GPU). However, the GPU speedup over the 42 CPU cores remains relatively unchanged (∼ 1.3×) even when using 4 MPI ranks per GPU (24 ranks in total) and 3 MPI ranks per GPU (18 ranks in total). This corresponds to a 35 %–40 % decrease in both simulation time and usage of resources. Due to too many local scalars and arrays in the W3SRCEMD subroutine and the huge WW3 memory requirement, GPU performance is currently limited by the data transfer bandwidth between the CPU and the GPU. Ideally, OpenACC routine directives could be used to further improve performance. However, W3SRCEMD would require significant code refactoring to make this possible. We also discuss how the trade-off between the occupancy, register, and latency affects the GPU performance of WW3.

58 GEOSCIENCES↗
DOE OSTI2025

Test–retest reliability for a social discounting of personal information task

Increasing cybercrime rates means identifying potential victims is critically important. Social discounting tasks show that individuals share less personally identifying information as social distance increases. However, the test–retest reliability and uniqueness of this measure is unclear. The current study assessed social discounting for personally identifying information (SDPII), delay discounting, risk taking, and personality at two measurement waves 30 days apart for 64 undergraduate students. Test–retest reliability was statistically significant for the SDPII and all other measures, replicating previous studies. SDPII rates were not significantly correlated with other measures during both measurement waves, showing discriminant validity. SDPII rates were lower than those reported in a previous study but were still well described by a hyperbolic discounting function, suggesting replicability across studies. Furthermore, the high test–retest reliability, uniqueness, and replicability of the SDPII suggests that it may quantitatively identify cybercrime victimization. Future research should test which measure or combination of measures can accurately predict scam and cybercrime victimization to inform data-based interventions.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Sound velocities and thermal equation of state of fcc -iron-nickel alloys at high pressure and high temperature: Implications for the cores of Moon and several planets

Fcc-Fe-Ni alloy is believed to be the most dominant solid constitute of moderate-sized terrestrial planetary cores. Investigating the physical properties, especially the density and sound velocity of Fe-Ni alloys and comparing them with seismic observations is an indispensable approach to constructing compositional models for planetary interiors. In this study, we conducted sound velocity measurements on Fe-Ni alloys with 10 wt.% and 20 wt.% Ni up to ∼13.5 GPa and 1073 K, using the ultrasonic interferometry technique in a multi-anvil apparatus in conjunction with synchrotron radiation. By fitting the experimental data to finite strain equations, the bulk and shear moduli and their pressure and temperature derivatives are derived, yielding K S0 =145.8(14) GPa, G 0 = 73.2(7) GPa, K S0 ’ = 5.89(24), G 0 ’ = 2.89(8), (∂K S /∂T) P = -0.0181(12) GPa/K and (∂G/∂T) P = -0.0393(10) GPa/K for fcc-Fe 80 Ni 20 . An examination of the density-velocity relationship shows that compressional wave velocity is insensitive to temperature within the current pressure and temperature range, while shear wave velocity exhibits a large reduction with increasing temperature. Here, extrapolation of the sound velocities following the finite strain theories suggests that much slower Vs should be expected at pressure and temperature conditions corresponding to those of the lunar core. Possible core density and velocity profiles for other moderate planets and satellites, such as Mars, Mercury, and Ganymede are also calculated.

Equation of state↗
DOE OSTI2023

Finite-difference time-domain methods

The finite-difference time-domain (FDTD) method is a widespread numerical tool for full-wave analysis of electromagnetic fields in complex media and for detailed geometries. Applications of the FDTD method cover a range of time and spatial scales, extending from subatomic to galactic lengths and from classical to quantum physics. Technology areas that benefit from the FDTD method include biomedicine — bioimaging, biophotonics, bioelectronics and biosensors; geophysics — remote sensing, communications, space weather hazards and geolocation; metamaterials — sub-wavelength focusing lenses, electromagnetic cloaks and continuously scanning leaky-wave antennas; optics — diffractive optical elements, photonic bandgap structures, photonic crystal waveguides and ring-resonator devices; plasmonics — plasmonic waveguides and antennas; and quantum applications — quantum devices and quantum radar. This Primer summarizes the main features of the FDTD method, along with key extensions that enable accurate solutions to be obtained for different research questions. Additionally, hardware considerations are discussed, plus examples of how to extract magnitude and phase data, Brillouin diagrams and scattering parameters from the output of an FDTD model. Furthermore, the Primer ends with a discussion of ongoing challenges and opportunities to further enhance the FDTD method for current and future applications.

42 ENGINEERING↗
DOE OSTI2022

Density reconstruction in convergent high-energy-density systems using x-ray radiography and Bayesian inference

X-ray radiography is a technique frequently used to diagnose convergent high-energy-density (HED) systems, such as inertial confinement fusion implosion, and provides unique information that is not available through self-emission measurements. In this work, we investigate the scope and limits of that information using a radiography simulation combined with the Bayesian inference workflow. The accuracy of density reconstruction from simulated radiographs of spherical implosions driven with 27 kJ laser energy is assessed, including the increase or decrease in accuracy due to the addition of Lagrangian marker layers, Poisson noise, and improved prior information. This work is the first to present the full uncertainty distributions inferred from radiography analysis in HED systems and demonstrates the importance of constructing the full posterior probability density, as opposed to a point estimate, due to the modal structure of the likelihood surface introduced by typical experimental noise sources. This general methodology can be used both for robust analysis of radiographic data and for an improved design of radiography experiments by modeling the full experimental system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Applications of physics informed neural operators

Abstract We present a critical analysis of physics-informed neural operators (PINOs) to solve partial differential equations (PDEs) that are ubiquitous in the study and modeling of physics phenomena using carefully curated datasets. Further, we provide a benchmarking suite which can be used to evaluate PINOs in solving such problems. We first demonstrate that our methods reproduce the accuracy and performance of other neural operators published elsewhere in the literature to learn the 1D wave equation and the 1D Burgers equation. Thereafter, we apply our PINOs to learn new types of equations, including the 2D Burgers equation in the scalar, inviscid and vector types. Finally, we show that our approach is also applicable to learn the physics of the 2D linear and nonlinear shallow water equations, which involve three coupled PDEs. We release our artificial intelligence surrogates and scientific software to produce initial data and boundary conditions to study a broad range of physically motivated scenarios. We provide the source code , an interactive website to visualize the predictions of our PINOs, and a tutorial for their use at the Data and Learning Hub for Science .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

AT2017gfo: Bayesian inference and model selection of multicomponent kilonovae and constraints on the neutron star equation of state

The joint detection of the gravitational wave GW170817, of the short γ-ray burst GRB170817A and of the kilonova AT2017gfo, generated by the the binary neutron star (NS) merger observed on 2017 August 17, is a milestone in multimessenger astronomy and provides new constraints on the NS equation of state. We perform Bayesian inference and model selection on AT2017gfo using semi-analytical, multicomponents models that also account for non-spherical ejecta. Observational data favour anisotropic geometries to spherically symmetric profiles, with a log-Bayes’ factor of ~104, and favour multicomponent models against single-component ones. The best-fitting model is an anisotropic three-component composed of dynamical ejecta plus neutrino and viscous winds. Using the dynamical ejecta parameters inferred from the best-fitting model and numerical–relativity relations connecting the ejecta properties to the binary properties, we constrain the binary mass ratio to q < 1.54 and the reduced tidal parameter to $120\lt \tilde{\Lambda }\lt 1110$. Finally, we combine the predictions from AT2017gfo with those from GW170817, constraining the radius of a NS of 1.4 M ⊙ to 12.2 ± 0.5 km (1σ level). This prediction could be further strengthened by improving kilonova models with numerical-relativity information.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Utah FORGE 2439: A Multi-Component Approach to Characterizing In-Situ Stress

Core-based in-situ stress estimation, Triaxial Ultrasonic Velocity (labTUV) data, and Deformation Rate Analysis (DRA) data for Utah FORGE well 16A(78)-32 using triaxial ultrasonic velocity and deformation rate analysis. Report documenting a multi-component approach to characterizing in-situ stress at the U.S. DOE FORGE EGS site: laboratory, modeling and field measurement. Core-based methods for in-situ stress estimation were applied using samples from 5 intervals within the Utah FORGE 16A(78)-32 well. At three of these locations, Triaxial Ultrasonic Velocity (labTUV) tests were performed, resulting in experimentally-determined relationships between wave velocities and stresses. Non-monotonic increase in the velocity-stress relationships are inferred provide evidence of stress history and are therefore used to estimate in-situ stress magnitudes. Additionally, Deformation Rate Analysis (DRA) tests were run on core plugs from various orientations at each of the 5 sampling locations. These, too, provide evidence of stress history based on stress-strain behavior. A novel Weight of Evidence (WoE) method was developed as a means of synthesizing in-situ stress evidence from these two types of tests. Results indicate the minimum horizontal stress gradient ranges from 0.58 psi/ft to 0.69 psi/ft, with 4 of the 5 values between 0.66 psi/ft and 0.69 psi/ft. The vertical stress gradient ranges from 1.05 psi/ft to 1.12 psi/ft, with 4 of the 5 zones given results between 1.09 psi/ft and 1.12 psi/ft. The maximum horizontal stress gradient ranges from 0.98 psi/ft to 1.34 psi/ft, with 4 of the 5 zones falling between 0.98 psi/ft and 1.24 psi/ft. The stress regime thus appears to be on the edge between normal faulting and strike-slip faulting, potentially flipping back and forth between the two regimes due to variability of rock properties, structures such as faults, and/or thermal anomalies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

Dark initial state radiation and the kinetic mixing portal

Data from Planck measurements of the cosmic microwave background (CMB) place important constraints on models with light dark matter (DM) and light mediators especially when both lie in the mass range below ~ 1 GeV. In models involving kinetic mixing where the dark photon acts as the mediator, these constraints are easily satisfied and the appropriate DM relic density achievable if the DM is, e.g., a complex scalar, where p-wave annihilation occurs, or is the lighter component of a split pseudo-Dirac state where co-annihilation dominates. In both of these cases, although higher order in the dark gauge coupling, g D , the corresponding annihilation processes including dark photon initial state radiation (ISR) will be dominantly s-wave with essentially temperature independent cross sections. The rates for these dark ISR associated processes, though not yielding cross sections large enough to contribute to the relic density, can still run into possible conflicts with the bounds arising from the CMB. In this paper we perform a preliminary study of the present and potential future constraints that the CMB imposes on the parameter spaces for both of these scenarios due to the existence of this dark ISR. Further analyses of the effects of dark ISR in DM annihilation is clearly warranted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

CO 2 induced changes in Mount Simon sandstone: Understanding links to post CO 2 injection monitoring, seismicity, and reservoir integrity

The purpose of this study is to quantify geochemical reactions of CO 2 and brine with subsurface samples taken from the Mt. Simon sandstone and identify any potential alterations of the geomechanical rock properties that could lead to changes observable in seismic monitoring or result in changes of micro seismicity such as those observed at the Illinois Basin-Decatur Project (IBDP) site. Two Mt. Simon sandstone samples from 6919.3 feet (2109.0 m) and 6926.1 feet (2111.1 m) depth were exposed to supercritical CO 2 (scCO 2 ) dissolved in brine at in-situ reservoir conditions for one month. Geochemical, spectral, scanning electron microscopy, and petrophysical methods were used to analyze the samples before and after the one-month exposure. Significant changes were observed. Multiple geomechanical properties were chosen to form the framework against which to interrogate the petrophysical data: Young’s modulus (E), Poisson’s ratio (ν), lambda·rho (λρ), and mu·rho (μρ). In this study we conclude that framework composition, porosity, heterogeneities, effective pressure, and reactive geochemistry are first order controls on trends in the E-μ and λρ-μρ cross plot spaces. Changes in porosity, permeability, dynamic moduli, and brittleness with exposure to these fluids were observed. No change in ultrasonic P-wave attenuation (Q P ) was observed. Geochemical alteration causes a distinct shift in λρ-μρ in both samples as well as changes in E, ν, and P and S wave seismic velocity values. Furthermore, these observations could provide insight into subsurface monitoring using seismic methods including amplitude variation with offset (AVO) classification.

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

The Patterson function as auto-hologram and graph enables the direct solution to the phase problem for coherently illuminated atomistic structures

Abstract The coherent Patterson function, derived from the coherent scattering intensities from N atoms, can be described as an auto-holographic image, i.e., the superposition of N -holographic images in which each atom serves as the source of a different reference wave, and as a mathematical graph of the unknown structure. These insights indicate that the unknown structure is significantly over-determined by the inherent information content in a Patterson function by a factor of ( N − 1)/2 (for non-degenerate structures), independent of dimensionality. However, we also show that the ability to resolve the distinct features in the Patterson function depends strongly on details of the experimental design (data range, sampling frequency, and dimensionality). This re-interpretation suggests that the coherent Patterson function provides a natural context for describing the information content in coherent scattering, reveals that there is no inherent phase problem for coherently illuminated pointwise structures, and enables the design of an algorithm which retrieves the structure directly from the Patterson function, without using error metric minimization.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Molecular Viscosity and Diffusivity Effects in Transitional and Shock-Driven Mixing Flows

Here, this paper investigates the importance of molecular viscosity and diffusivity for the prediction of transitional and shock-driven mixing flows featuring high and low Reynolds and Mach number regions. Two representative problems are computed with implicit large-eddy simulations using the inviscid Euler equations (EE) and viscous Navier-Stokes equations (NSE): the Taylor-Green vortex at Reynolds number Re = 3000 and initial Mach number Ma = 0.28, and an air- SF 6 -air gas curtain subjected to two shock waves at Ma = 1.2. The primary focus is on differences between NSE and EE predictions due to viscous effects. The outcome of the paper illustrates the advantages of utilizing NSE. In contrast to the EE, where the effective viscosity decreases upon grid refinement, NSE predictions can be assessed for simulations of flows with transition to turbulence at prescribed constant Re. The NSE can achieve better agreement between solutions and reference data, and the results converge upon grid refinement. On the other hand, the EE predictions do not converge with grid refinement, and can only exhibit similarities with the NSE results at coarse grid resolutions. We also investigate the effect of viscous effects on the dynamics of the coherent and turbulent fields, as well as on the mechanisms contributing to the production and diffusion of vorticity. The results show that nominally inviscid calculations can exhibit significantly varying flow dynamics driven by changing effective resolution-dependent Reynolds number, and highlight the role of viscous processes affecting the vorticity field. These tendencies become more pronounced upon grid refinement. The discussion of the results concludes with the assessment of the computational cost of inviscid and viscous computations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Knight Shift and Leading Superconducting Instability from Spin Fluctuations in Sr 2 RuO 4

We report that recent nuclear magnetic resonance studies [A. Pustogow et al. , Nature 574, 72 (2019)] have challenged the prevalent chiral triplet pairing scenario proposed for Sr 2 RuO 4 . To provide guidance from microscopic theory as to which other pair states might be compatible with the new data, we perform a detailed theoretical study of spin fluctuation mediated pairing for this compound. We map out the phase diagram as a function of spin-orbit coupling, interaction parameters, and band structure properties over physically reasonable ranges, comparing when possible with photoemission and inelastic neutron scattering data information. We find that even-parity pseudospin singlet solutions dominate large regions of the phase diagram, but in certain regimes spin-orbit coupling favors a near-nodal odd-parity triplet superconducting state, which is either helical or chiral depending on the proximity of the γ band to the van Hove points. A surprising near degeneracy of the nodal s ' and d x 2 - y 2 wave solutions leads to the possibility of a near-nodal time-reversal symmetry broken s ' + i d x 2 - y 2 pair state. Lastly, predictions for the temperature dependence of the Knight shift for fields in and out of plane are presented for all states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Data-Driven Transition Path Analysis Yields a Statistical Understanding of Sudden Stratospheric Warming Events in an Idealized Model

Abstract Atmospheric regime transitions are highly impactful as drivers of extreme weather events, but pose two formidable modeling challenges: predicting the next event (weather forecasting) and characterizing the statistics of events of a given severity (the risk climatology). Each event has a different duration and spatial structure, making it hard to define an objective “average event.” We argue here that transition path theory (TPT), a stochastic process framework, is an appropriate tool for the task. We demonstrate TPT’s capacities on a wave–mean flow model of sudden stratospheric warmings (SSWs) developed by Holton and Mass, which is idealized enough for transparent TPT analysis but complex enough to demonstrate computational scalability. Whereas a recent article (Finkel et al. 2021) studied near-term SSW predictability, the present article uses TPT to link predictability to long-term SSW frequency. This requires not only forecasting forward in time from an initial condition, but also backward in time to assess the probability of the initial conditions themselves. TPT enables one to condition the dynamics on the regime transition occurring, and thus visualize its physical drivers with a vector field called the reactive current . The reactive current shows that before an SSW, dissipation and stochastic forcing drive a slow decay of vortex strength at lower altitudes. The response of upper-level winds is late and sudden, occurring only after the transition is almost complete from a probabilistic point of view. This case study demonstrates that TPT quantities, visualized in a space of physically meaningful variables, can help one understand the dynamics of regime transitions.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

The propagation of seismic waves, misinformation, and disinformation from the 2024-10-05 M 4.5 Iran earthquake

The 2024-10-05 Iran M 4.5 earthquake took place at a time of heightened tensions in the Middle East. We perform a discrimination and moment tensor analysis and identify a shallow-dipping, reverse fault source commensurate with the compressional setting of the Iranian interior. Nonetheless, the event's aftermath saw widespread dissemination of misinformation, and potentially active disinformation, concluding that it was in fact a test of an Iranian nuclear weapon. The 'evidence' for many of these claims was based on inaccurate interpretation of seismic data. In this paper, we analyze how geophysical 'fake news' propagated through social media (mainly Twitter/X) following this event, eventually gaining traction in mainstream, earned media. This event is an illustrative warning of how seismic data can be misinterpreted and/or manipulated in public discourse.

58 GEOSCIENCES↗
DOE OSTI2022

ENSO Impact on Winter Precipitation in the Southeast United States through a Synoptic Climate Approach

The ENSO impact on winter precipitation in the Southeast United States was analyzed from the perspective of daily weather types (WTs). We calculated the dynamic contribution associated with the change in frequency of the WTs and the thermodynamic contribution due to changes in the spatial patterns of the environmental fields of the WTs. Six WTs were obtained using a k-means clustering analysis of 850 hPa winds in reanalysis data from November to February of 1948–2022. All the WTs can only persist for a few days. The most frequent winter weather type is WT1 (shallow trough in Eastern U.S.), which can persist or likely transfer to WT4 (Mississippi River Valley ridge). WT1 becomes less frequent in El Niño years, while the frequency of WT4 does not change much. WTs 2–6 correspond to a loop of eastward propagating waves with troughs and ridges in the mid-latitude westerlies. Three WTs with a deep trough in the Southeast U.S., which are WT2 (east coast trough), WT3 (off east coast trough) and WT6 (plains trough), become more frequent in El Niño years. The more frequent deep troughs (WTs 2, 3 and 6) and less frequent shallow trough (WT1) result in above-normal precipitation in the coastal Southeast U.S. in the winter of El Niño years. WT5 (off coast Carolina High), with maximum precipitation extending from Mississippi Valley to the Great Lakes, becomes less frequent in El Niño years, which corresponds to the below-normal precipitation from the Great Lakes to Upper Mississippi and Ohio River Valley in El Niño years, and vice versa in La Niña years. The relative contribution of the thermodynamic and dynamic contribution is location dependent. On the east coast, the two contributions are similar in magnitude.

Qian, Jian-Hua↗
DOE OSTI2021

COVID19 Disease Map, a computational knowledge repository of virus–host interaction mechanisms

We need to effectively combine the knowledge from surging literature with complex datasets to propose mechanistic models of SARS-CoV-2 infection, improving data interpretation and predicting key targets of intervention. Here, we describe a large-scale community effort to build an open access, interoperable and computable repository of COVID-19 molecular mechanisms. The COVID-19 Disease Map (C19DMap) is a graphical, interactive representation of disease-relevant molecular mechanisms linking many knowledge sources. Notably, it is a computational resource for graph-based analyses and disease modelling. To this end, we established a framework of tools, platforms and guidelines necessary for a multifaceted community of biocurators, domain experts, bioinformaticians and computational biologists. The diagrams of the C19DMap, curated from the literature, are integrated with relevant interaction and text mining databases. We demonstrate the application of network analysis and modelling approaches by concrete examples to highlight new testable hypotheses. This framework helps to find signatures of SARS-CoV-2 predisposition, treatment response or prioritisation of drug candidates. Such an approach may help deal with new waves of COVID-19 or similar pandemics in the long-term perspective.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Unitarity methods in AdS/CFT

We develop a systematic unitarity method for loop-level AdS scattering amplitudes, dual to non-planar CFT correlators, from both bulk and boundary perspectives. We identify cut operators acting on bulk amplitudes that put virtual lines on shell, and show how the conformal partial wave decomposition of the amplitudes may be efficiently computed by gluing lower-loop amplitudes. A central role is played by the double discontinuity of the amplitude, which has a direct relation to these cuts. We then exhibit a precise, intuitive map between the diagrammatic approach in the bulk using cutting and gluing, and the algebraic, holographic unitarity method of [1] that constructs the non-planar correlator from planar CFT data. Our analysis focuses mostly on four-point, one-loop diagrams — we compute cuts of the scalar bubble, triangle and box, as well as some one-particle reducible diagrams — in addition to the five-point tree and four-point double-ladder. Analogies with S-matrix unitarity methods are drawn throughout.

79 ASTRONOMY AND ASTROPHYSICS↗