Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained wave”

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.

358 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Inference of the optical depth to reionization τ from Planck CMB maps with convolutional neural networks

The optical depth to reionization, τ, is the least constrained parameter of the cosmological Λ cold dark matter (ΛCDM) model. To date, its most precise value is inferred from large-scale polarized cosmic microwave background (CMB) power spectra from the High Frequency Instrument (HFI) aboard the Planck satellite. These maps are known to contain significant contamination by residual non-Gaussian systematic effects, which are hard to model analytically. Therefore, robust constraints on τ are currently obtained through an empirical cross-spectrum likelihood built from simulations. In this paper, we present a likelihood-free inference of τ from polarized Planck HFI maps which, for the first time, is fully based on neural networks (NNs). NNs have the advantage of not requiring an analytical description of the data and can be trained on state-of-the-art simulations, combining the information from multiple channels. By using Gaussian sky simulations and Planck SRoll2 simulations, including CMB, noise, and residual instrumental systematic effects, we trained, tested, and validated NN models considering different setups. We inferred the value of τ directly from Stokes Q and U maps at ~4° pixel resolution, without computing angular power spectra. On Planck data, we obtained τ NN = 0.0579 ± 0.0082, which is compatible with current EE cross-spectrum results but with a ~30% larger uncertainty, which can be assigned to the inherent nonoptimality of our estimator and to the retraining procedure applied to avoid biases. While this paper does not improve on current cosmological constraints on τ, our analysis represents a first robust application of NN-based inference on real data, and highlights its potential as a promising tool for complementary analysis of near-future CMB experiments, also in view of the ongoing challenge to achieve the first detection of primordial gravitational waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Temperature-induced densification in compressed basaltic glass revealed by in-situ ultrasonic measurements

Abstract Acoustic velocities of a model basalt glass (64 mol% CaMgSi2O6 + 36 mol% CaAl2Si2O8) were measured along different pressure-temperature (P-T) paths. One set of experiments involved isothermal compression-decompression cycles, performed at temperatures of 300, 641, 823, and 1006 K and pressures up to 12.2 GPa. The other set of experiments involved constant-load heating-cooling cycles at temperatures up to 823 K and pressures up to 7.5 GPa. Both sets of experiments were performed in a multi-anvil apparatus using a synchrotron-based ultrasonic technique. Our results show that the glass compressed isothermally at 300 K (cold-compression) displays anomalously decreasing compressional (VP) and shear (VS) wave velocities with increasing pressure until ~8 GPa. Beyond 8 GPa, both VP and VS start to increase sharply with pressure and irreversible densification of the glass occurred, producing large hysteresis loops of velocities upon decompression. However, for the glass compressed isothermally at increasingly higher temperatures (hot-compression), the velocity minima gradually shift to lower pressures. At temperature close to the glass transition temperature Tg, the velocity minima disappear completely, displaying a monotonic increase of velocities during compression and higher VP and VS during decompression. In addition, constant-load heating-cooling experiments show that velocities generally decrease slightly with increasing temperature, but start to increase once heated above a threshold temperature (~650 K). During cooling the velocities increase almost linearly with decreasing temperature, resulting in higher velocities (~1.5–2.5% higher) when returned to 300 K. This implies that a temperature-induced densification may have occurred in the glass at high pressures. Raman spectra on recovered samples show that the hot-compressed and high-P heated glasses contain distinctly densified and depolymerized structural signatures compared to the initial glass and the cold-compressed glass below the velocity transition pressure PT (~8 GPa). Such densification may be attributed to the breaking of bridging oxygen bonds and compaction in the intermediate-range structure. Our results demonstrate that temperature can facilitate glass densification at high pressures and point out the importance of P-T history in understanding the elastic properties of silicate glasses. Comparison with melt velocity suggests that hot-compressed glasses may better resemble the pressure dependence of velocity of silicate melts than cold-compressed glasses, but still show significantly higher velocities than melts. If the abnormal acoustic behaviors of cold-compressed glasses were used to constrain melt fractions in the mantle low-velocity regions, the melt fractions needed to explain a given velocity reduction would be significantly underestimated at high pressures.

Geochemistry & Geophysics↗
DOE OSTI2023

Uncertainty-quantified phenomenological optical potentials for single-nucleon scattering

Optical-model potentials (OMPs) continue to play a key role in nuclear reaction calculations. However, the uncertainty of phenomenological OMPs in widespread use—inherent to any parametric model trained on data—has not been fully characterized, and its impact on downstream users of OMPs remains unclear. Here we assign well-calibrated uncertainties for two representative global OMPs, those of Koning-Delaroche and Chapel Hill '89, using Markov-chain Monte Carlo for parameter inference. By comparing the canonical versions of these OMPs against the experimental data originally used to constrain them, we show how a lack of outlier rejection and a systematic underestimation of experimental uncertainties contributes to bias of, and overconfidence in, best-fit parameter values. Our updated, uncertainty-quantified versions of these OMPs address these issues and yield complete covariance information for potential parameters. Scattering predictions generated from our ensembles show improved performance both against the original training corpora of experimental data and against a new “test” corpus comprising many of the experimental single-nucleon scattering data collected over the last twenty years. Finally, we apply our uncertainty-quantified OMPs to two case studies of application-relevant cross sections. We conclude that, for many common applications of OMPs, including OMP uncertainty should become standard practice. Furthermore, to facilitate their immediate use, digital versions of our updated OMPs and related tools for forward uncertainty propagation are included as Supplemental Material.

150 ≤ A ≤ 189↗
DOE OSTI2024

Classification of compact objects and model comparison using EOS knowledge

Nuclear theory and experiments, alongside astrophysical observations, constrain the equation of state (EOS) of supranuclear-dense matter. Conversely, knowledge of the EOS allows an improved interpretation of nuclear or astrophysical data. In this article, we use several established constraints on the EOS and the new NICER measurement of PSR J0437-4715 to comment on the nature of the primary companion in GW230529 and the companion of PSR J0514-4002E. We find that, with a probability of ≳84% and ≳68%, respectively, both objects are black holes. These likelihoods increase to above 95% when one uses GW170817’s remnant as an upper limit on the TOV mass. We also demonstrate that the current knowledge of the EOS substantially disfavors high masses and radii for PSR J⁢0030+0451, inferred recently when combining NICER with XMM-Newton background data and using particular hot-spot models. Lastly, we also use our obtained EOS knowledge to comment on measurements of the nuclear symmetry energy, finding that the large value predicted by the PREX-II measurement displays some mild tension with other constraints on the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Uncertainty quantification of optical models in fission fragment deexcitation

Here, we take the first step towards incorporating compound nuclear observables at astrophysically relevant energies into the experimental evidence used to constrain optical models, by propagating the uncertainty in two global optical potentials, one phenomenological and one microscopic, to correlated fission observables using the Monte Carlo Hauser-Feshbach formalism. We compare to a wide range of historic and recent experimental fission measurements, and discuss in detail regions of disagreement. We find that the parametric optical model uncertainty in neutron-fragment correlated observables involving neutron energy is significant. On the other hand, we observe that other experimental features, particularly neutron-fragment correlations near the 132 Sn shell closure and the high energy component of neutron spectra, are unlikely to be explained by the optical potential, and will require further experimental and theoretical effort to explain.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Dynamical ejecta from binary neutron star mergers: Impact of a small residual eccentricity and of the equation of state implementation

Predicting the properties of the matter ejected during and after a neutron star merger is crucial to our ability to use electromagnetic observations of these mergers to constrain the masses of the neutron stars, the equation of state of dense matter, and the role of neutron star mergers in the enrichment of the Universe in heavy elements. Furthermore, our ability to reliably provide such predictions is however limited by a broad range of factors, including the finite resolution of numerical simulations, their treatment of magnetic fields, neutrinos, and neutrino-matter interactions, and the approximate modeling of the equation of state of dense matter. In this manuscript, we study specifically the role that a small residual eccentricity and different implementations of the same equation of state have on the matter ejected during the merger of a 1.3M ⊙ –1.4M ⊙ binary neutron star system. We find that a residual eccentricity e ~ 0.01, as measured ~ 4–6 orbits before merger, causes O(25%–30%) changes in the amount of ejected mass, mainly due to changes in the amount of matter ejected as a result of core bounces during merger. We note that O(1%) residual eccentricities have regularly been used in binary neutron star merger simulations as proxy for circular binaries, potentially creating an additional source of error in predictions for the mass of the dynamical ejecta.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Renormalized classical theory of quantum magnets

Here, we derive a renormalized classical spin (RCS) theory for 𝑆 >1/2 quantum magnets by constraining a generalized classical theory that includes all multipolar fluctuations to a reduced CP 1 phase space of dipolar SU(2) coherent states. When the spin Hamiltonian $\hat{ℋ}$(𝑆) is linear in the spin operators $\hat{𝑺}$ 𝑗 for each lattice site 𝑗, the RCS Hamiltonian $\tilde{ℋ}$ cl coincides with the usual classical model ℋ cl = lim 𝑆→∞⁡ $\hat{ℋ}$(𝑆). In the presence of nonlinear terms, however, the RCS theory is more accurate than ℋ cl . For the many materials modeled by spin Hamiltonians with (nonlinear) single-ion anisotropy terms, the use of the RCS theory is essential to accurately model phase diagrams and to extract the correct Hamiltonian parameters from neutron-scattering data.

magnetic anisotropy↗
DOE OSTI2026

Asymptotic-preserving semi-implicit finite volume scheme for extended magnetohydrodynamics

A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge–Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Direct nonparametric multimessenger constraints on the equation of state of cold dense nuclear matter

We utilize the now substantial amount of astrophysical observations of neutron stars (NSs), along with perturbative quantum chromodynamics (pQCD) calculations at high density, to directly constrain the NS equation of state (EOS). To this end, we construct nonparametric EOS priors by using Gaussian processes trained on 75 EOSs, which include models with either hadrons, hyperons, or quarks at high densities. We create a prior using the full EOS sample (model agnostic), and one prior for each EOS family to test model discrimination. We introduce a novel inference approach, which allows the simultaneous sampling of intrinsic and extrinsic parameters of binary NS mergers, as well as a nonparametric equation of state. We showcase this method in a Bayesian updating scheme by first performing a complete analysis of the binary NS merger event GW170817 with minimal assumptions, and sequentially adding information from x-ray and radio NS observations, along with pQCD calculations. Besides providing standard constraints, such as the pressure at twice nuclear saturation density 𝑝⁡(2⁢𝜌 sat ) = 4.3$^{+0.6}_{−0.6}$ × 10 34 dyne/cm 2 , at 95% confidence level, for the model agnostic prior, our methodology shows how the choice of EOS families used in conditioning changes the inferred astrophysical properties of the EOS, namely tidal deformability and maximum supported NS mass. We find hyperonic priors predicting higher tidal deformabilities for a 1.4⁢𝑀 ⊙ NS, and hadronic priors being preferred by the considered astrophysical data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Spurious solar-wind effects on acceleration noise in LISA Pathfinder

Spurious solar-wind effects are a potential noise source in future Laser Interferometer Space Antenna (LISA) measurements. One noise coupling mechanism is constrained by estimating solar-wind effects on acceleration noise in LISA Pathfinder (LPF). While LISA is designed for drag-free differential measurement, predicting the realistic impact both bounds the operational environment and assesses whether LISA could provide serendipitous space-weather observations. Data from NASA's Advanced Composition Explorer (ACE), situated at the L1 Lagrange point, serves as a reliable source of solar-wind data. The data sets are compared over the 114 d time period from 1 March 2016 to 23 June 2016. This period gives the longest readily-available open data set, without interference from other commissioning activities. To evaluate space weather effects, the data from both satellites are formatted, gap-filled/interpolated, and fast-Fourier transformed for amplitude spectral density and coherence comparisons. Solar wind effects are not seen in a coherence plot between LPF and ACE; modest coherence in the planned LISA observational frequency band can be attributed to chance. This result indicates that measurable correlation due to solar-wind acceleration noise over 3 month timescales will be a negligible noise source. LISA is unlikely to inform solar wind measurements routinely. Another source of noise from the Sun, solar radiation pressure, is estimated to impart greater acceleration noise, but has yet to be analyzed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Modeling dust mineralogical composition: sensitivity to soil mineralogy atlases and their expected climate impacts

Soil dust aerosols are a key component of the climate system, as they interact with short- and long-wave radiation, alter cloud formation processes, affect atmospheric chemistry and play a role in biogeochemical cycles by providing nutrient inputs such as iron and phosphorus. The influence of dust on these processes depends on its physicochemical properties, which, far from being homogeneous, are shaped by its regionally varying mineral composition. The relative amount of minerals in dust depends on the source region and shows a large geographical variability. However, many state-of-the-art Earth system models (ESMs), upon which climate analyses and projections rely, still consider dust mineralogy to be invariant. The explicit representation of minerals in ESMs is more hindered by our limited knowledge of the global soil composition along with the resulting size-resolved airborne mineralogy than by computational constraints. In this work we introduce an explicit mineralogy representation within the state-of-the-art Multiscale Online Nonhydrostatic AtmospheRe CHemistry (MONARCH) model. We review and compare two existing soil mineralogy datasets, which remain a source of uncertainty for dust mineralogy modeling and provide an evaluation of multiannual simulations against available mineralogy observations. Soil mineralogy datasets are based on measurements performed after wet sieving, which breaks the aggregates found in the parent soil. Our model predicts the emitted particle size distribution (PSD) in terms of its constituent minerals based on brittle fragmentation theory (BFT), which reconstructs the emitted mineral aggregates destroyed by wet sieving. Our simulations broadly reproduce the most abundant mineral fractions independently of the soil composition data used. Feldspars and calcite are highly sensitive to the soil mineralogy map, mainly due to the different assumptions made in each soil dataset to extrapolate a handful of soil measurements to arid and semi-arid regions worldwide. For the least abundant or more difficult-to-determine minerals, such as iron oxides, uncertainties in soil mineralogy yield differences in annual mean aerosol mass fractions of up to ~ 100 %. Although BFT restores coarse aggregates including phyllosilicates that usually break during soil analysis, we still identify an overestimation of coarse quartz mass fractions (above 2 µm in diameter). In a dedicated experiment, we estimate the fraction of dust with undetermined composition as given by a soil map, which makes up ~ 10 % of the emitted dust mass at the global scale and can be regionally larger. Changes in the underlying soil mineralogy impact our estimates of climate-relevant variables, particularly affecting the regional variability of the single-scattering albedo at solar wavelengths or the total iron deposited over oceans. All in all, this assessment represents a baseline for future model experiments including new mineralogical maps constrained by high-quality spaceborne hyperspectral measurements, such as those arising from the NASA Earth Surface Mineral Dust Source Investigation (EMIT) mission.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Radiation and heat transport in divergent shock–bubble interactions

Shock–bubble interactions (SBIs) are important across a wide range of physical systems. In inertial confinement fusion, interactions between laser-driven shocks and micro-voids in both ablators and foam targets generate instabilities that are a major obstacle in achieving ignition. Experiments imaging the collapse of such voids at high energy densities (HED) are constrained by spatial and temporal resolution, making simulations a vital tool in understanding these systems. In this study, we benchmark several radiation and thermal transport models in the xRAGE hydrodynamic code against experimental images of a collapsing mesoscale void during the passage of a 300 GPa shock. We also quantitatively examine the role of transport physics in the evolution of the SBI. This allows us to understand the dynamics of the interaction at timescales shorter than experimental imaging framerates. We find that all radiation models examined reproduce empirical shock velocities within experimental error. Radiation transport is found to reduce shock pressures by providing an additional energy pathway in the ablation region, but this effect is small (∼1% of total shock pressure). Employing a flux-limited Spitzer model for heat conduction, we find that flux limiters between 0.03 and 0.10 produce agreement with experimental velocities, suggesting that the system is well-within the Spitzer regime. Higher heat conduction is found to lower temperatures in the ablated plasma and to prevent secondary shocks at the ablation front, resulting in weaker primary shocks. Finally, we confirm that the SBI-driven instabilities observed in the HED regime are baroclinically driven, as in the low energy case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Measurement of U 238 ( n , n ' γ ) cross section data and their impact on reaction models

A better knowledge of (n,xn) reaction cross sections is important for both reaction modeling and energy applications. This article focuses on inelastic scattering of neutrons off 238 U for which improvements are needed to better constrain evaluations and solve inconsistencies in nuclear power reactor calculations. A new precise measurement of (n,xnγ) reaction cross sections on 238 U has been performed at the GELINA (Geel Electron LINear Accelerator) neutron facility operated by EC-JRC-Geel (Belgium) with the GRAPhEME (GeRmanium array for Actinides PrEcise MEasurements) setup. The prompt γ-ray spectroscopy method coupled to time-of-flight measurements is used to extract (n,xnγ) cross section values which can be further combined to infer the total neutron inelastic scattering cross section. Cross section data for 18 γ transitions (five never measured before) are presented and compared to the data in the literature. Here, emphasis is especially given to the uncertainty determination to produce partial cross section data as accurate as possible. Due to intrinsic limitations of the experimental method, the use of additional nuclear structure information coupled with theoretical modeling is required to determine the total (n, n') cross section over the whole neutron energy range. We have investigated modeling aspects of the 238 U(n,n'γ) cross sections related to the description of compound nucleus and preequilibirum mechanisms as well as the discrete part of nuclear structure. Through comparison between experimental and calculated (n,n'γ) cross sections, we pinpoint inaccuracies in the description of specific reaction mechanisms and challenge recently implemented models. This helps improving the whole modeling of the (n, n') reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Observations of the Marine Atmospheric Boundary Layer’s Response to a Solar Eclipse

The atmospheric response to the solar eclipse of 8 April 2024 in North America is investigated with a specific focus on the marine atmospheric boundary layer (MABL). We leverage measurements collected during the Third Wind Forecast Improvement Project (WFIP3), including Doppler lidars, sonic anemometers, and thermodynamic profiler data to investigate the atmospheric response across sites that experienced partial eclipse conditions with nearly 90% obscuration. Using these measurements, we examine eclipse-induced changes in key meteorological parameters, such as temperature, wind speed, and turbulent fluxes. Most previous eclipse studies have been conducted over land, whereas this study provides new observations for both coastal and marine environments, offering additional insight into eclipse-driven variability in the MABL. The findings confirm a notable decrease in downwelling shortwave radiation during the eclipse, which results in rapid cooling of surface air. The temperature reduction ranges from $1.2^\circ \text {C}$ to $1.4^\circ \text {C}$ in coastal regions and from $0.3^\circ \text {C}$ to $0.5^\circ \text {C}$ over the ocean. This analysis suggests that the MABL’s higher thermal inertia compared to coastal regions moderates the temperature decrease during the eclipse. Wind speed exhibits a more complex behavior, as it is influenced by both the MABL and preexisting synoptic conditions. Although a reduction in wind speed is observable up to approximately 140 m above ground level (AGL) at more inland sites, at other locations closer to the coast, this reduction is constrained to the lowest 100 m AGL. Turbulence parameters retrieved from sonic anemometers, such as turbulence kinetic energy, turbulent heat flux, and friction velocity, decrease during the eclipse at coastal sites, accompanied by a brief transition of atmospheric stability from unstable to neutral or weakly stable conditions. For the open-ocean sites, the variability in turbulence statistics and atmospheric stability is minimal during the occurrence of the eclipse.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Direct experimental constraints on the spatial extent of a neutrino wavepacket

Despite their high relative abundance in our Universe, neutrinos are the least understood fundamental particles of nature. In fact, the quantum properties of neutrinos emitted in experimentally relevant sources are theoretically contested and the spatial extent of the neutrino wavepacket is only loosely constrained by reactor neutrino oscillation data with a spread of 13 orders of magnitude. Here we present a method to directly access this quantity by precisely measuring the energy width of the recoil daughter nucleus emitted in the radioactive decay of beryllium-7. The final state in the decay process contains a recoiling lithium-7 nucleus, which is entangled with an electron neutrino at creation. The lithium-7 energy spectrum is measured to high precision by directly embedding beryllium-7 radioisotopes into a high-resolution superconducting tunnel junction that is operated as a cryogenic sensor. Under this approach, we set a lower limit on the Heisenberg spatial uncertainty of the recoil daughter of 6.2 pm, which implies that the final-state system is localized at a scale more than a thousand times larger than the nucleus itself. From this measurement, the first, to our knowledge, direct lower limit on the spatial extent of a neutrino wavepacket is extracted. These results may have implications in several areas including the theoretical understanding of neutrino properties, the nature of localization in weak nuclear decays and the interpretation of neutrino physics data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗