Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density estimation”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quasielastic least-squares reverse time migration of PS reflections

Conventional analysis of amplitude variation with offset for elastic PS reflections is based on analytical reflection coefficients in a layered medium, and wave-equation-based PS migration is mainly used to produce a structural image. To overcome this problem, we have developed a least-squares reverse time migration (LSRTM) method for elastic PS reflections based on a quasielastic wave equation. Here, the quasielastic wave equation can accurately model PS reflections with elastic amplitudes under the first-order Born approximation. Our LSRTM method inverts for perturbations of the S-wave velocity and density by minimizing the L 2 norm of the difference between recorded and predicted PS reflections modeled using a quasielastic wave equation. We refer to our new method as quasielastic LSRTM of PS reflections. Numerical tests on synthetic and field data indicate that our method can properly handle the amplitudes of elastic PS reflections and provides an accurate estimate of the perturbations of S-wave velocity and density. Extending the method to the 3D case is not straightforward and might require incorporating certain PS data processing techniques into the inversion itself.

58 GEOSCIENCES↗
DOE OSTI2021

Shear reflectivity compensation in full-waveform inversion using least-squares reverse-time migration

SUMMARY The computational cost of elastic-waveform inversion is too high for inverting PP reflections, while using acoustic full-waveform inversion (FWI) is inaccurate because it does not depend on the shear modulus/velocity/impedance that affects elastic PP wavefield amplitudes. To solve this problem, we develop a waveform inversion method that uses acoustic least-squares reverse-time migration (LSRTM) to compensate the shear reflectivity for acoustic FWI. Our method is based on the quasi-elastic-wave equation developed by Chapman et al. (2014). The quasi-elastic-wave equation uses a linearized acoustic-wave equation with shear modulus μ as a virtual source to correct the acoustic PP wavefield amplitudes toward elastic ones. Our waveform inversion method inverts for elastic parameters by minimizing the L2 norm of the difference between recorded and predicted PP reflections modelled using the quasi-elastic-wave equation. Numerical tests on synthetic and field data show that our method can properly handle the amplitudes of elastic PP reflections and provides an accurate estimate of the P- and S-wave velocities/impedances and, in some cases, the density. The method does not need the computationally expensive numerical solution to the elastic-wave equation. It also gives a better estimate of elastic parameters than a pure LSRTM method for elastic PP reflections.

Feng, Zongcai↗
DOE OSTI2025

Carbon Dioxide Removal Measurement, Reporting, and Verification Simulation toolkit (CDR MRVSim) v0.1

CDR MRVSim is a statistical software toolkit for techno-economic analysis of measurement, reporting, and verification of multiple carbon dioxide removal technologies. The current version applies Monte Carlo simulation to existing datasets estimate the cost and uncertainty of measuring the soil organic carbon (SOC) content of agricultural land, accounting for multiple sources of uncertainty, using data from field measurements. We are planning to incorporate enhanced rock weathering and other technologies into the software. The current model leverages existing data to characterize underlying variability in soil organic carbon and related parameters, including bulk density, in an agricultural field attempting to increase its SOC. We then simulate baseline and post-intervention "measurement campaigns" in which some number of SOC and bulk density measurements are conducted using a selected technology. This enables estimates of the cost and accuracy of measuring changes in SOC in the simulated field. The primary advantages over similar software are: 1) ability to quantify tradeoffs between cost and uncertainty in MRV across a range of possible MRV approaches 2) focus on guiding development of novel sensors by determining desirable sets of characteristics

Sherwin, Evan [Lawrence Berkeley National Laborato↗
DOE OSTI2022

High Flux Isotope Reactor Low Enriched Uranium U-10Mo Fuel Design Parameters

Activities to convert the HFIR from HEU to LEU are ongoing as part of the US Department of Energy (DOE) National Nuclear Security Administration (NNSA) nuclear nonproliferation mission. Design activities to study the conversion of HFIR from HEU to LEU fuel explored different fuel design features and shapes with a uranium-molybdenum (U-10Mo) monolithic alloy fuel. This high-density alloy contains 90 wt % uranium and 10 wt % molybdenum and has a uranium density of 15.318gU/cm 3 . The goal of these studies is to generate several candidate HFIR LEU fuel designs of varying fuel fabrication complexity that meet the current HEU performance metrics and safety requirements. Recent advancements in modeling and simulation tools and design methods enabled a thorough analysis of the available design space with U-10Mo fuel. A surrogate model used this analysis as training data to quickly determine the performance of a design given specific design parameters. An optimization module used this surrogate model to quickly search this multidimensional search space given specific desired performance characteristics. This approach was made possible by the large available design space with U-10Mo fuel. Shift, a Monte Carlo tool optimized for high-performance computing (HPC) architectures, was used for faster calculation and better data management for reactor physics simulations. Once most of these design studies were complete, a new suite called the Python HFIR Analysis and Measurement Engine (PHAME) was developed to connect all fuel design analysis steps, making design studies more efficient and reproducible. The post-processing capabilities of these new tools are leveraged for the information provided herein. Leveraging these tools, several candidate fuel designs were selected with varying levels of feature complexity and reactor performance. This report provides design feature details for four selected HFIR LEU U-10Mo fuel designs and their corresponding performance and safety metrics. Nominal best-estimate design parameters and irradiation conditions, including fission rate densities, power densities, heat fluxes, and cumulative fission densities, are provided. Simulations show that the high uranium density of U-10Mo fuel provides a large potential design space that enables various LEU designs to meet HEU core performance metrics and safety requirements with a power increase from 85 MW (HEU) to 95 MW or 100 MW (LEU).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Accurate shear estimation with fourth-order moments

ABSTRACT As imaging surveys progress in exploring the large-scale structure of the Universe through the use of weak gravitational lensing, achieving sub-per cent accuracy in estimating shape distortions caused by lensing, or shear, is imperative for precision cosmology. In this paper, we extend the Fourier power function shapelets (FPFS) shear estimator using fourth-order shapelet moments and combine it with the original second-order shear estimator to reduce galaxy shape noise. We calibrate this novel shear estimator analytically to a sub-per cent level-accuracy using the AnaCal framework. This higher order shear estimator is tested with realistic image simulations, and after analytical correction for the detection/selection bias and noise bias, the multiplicative shear bias $|m|$ is below $3\times 10^{-3}$ (99.7 per cent confidence interval) for both isolated and blended galaxies. Once combined with the second-order FPFS shear estimator, the shape noise is reduced by $\sim 35~{{\ \rm per\ cent}}$ for isolated galaxies in simulations with Hyper Suprime-Cam and Vera C. Rubin Observatory Legacy Survey of Space and Time observational conditions. However, for blended galaxies, the effective number density does not significantly improve with the combination of the two estimators. Based on these results, we recommend exploration of how this framework can further reduce the systematic uncertainties in shear due to point spread function leakage and modelling error, and potentially provide improved precision in shear inference in high-resolution space-based images.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Effects of secondary void-initiating particles on the steady-state crack growth resistance of high-strength steel

In this investigation, we focus on the dependence of stable crack growth resistance, as measured by the short-rod fracture toughness test, on secondary void-initiating particles such as AlN and Ti-based grain-refining precipitates in a variety of high-strength steels with tempered martensitic microstructures. The model developed by Ritchie and Thompson is modified to illustrate the significant amounts of toughening that can result from the refinement of secondary particles. Analysis of the data suggests that material strength is a predominant factor in increasing the short-rod fracture toughness relative to linear-elastic measures of initiation fracture toughness, but the extent of toughening is limited by the size and number density of secondary particles in the microstructure. The variation in estimates of secondary microvoid initiation and growth strains with precipitate size reinforce the notion that primary fracture at non-metallic inclusions and secondary fracture at smaller particles occur as sequential processes with a degree of concurrence that is dependent on the state of precipitation in both particle dispersions. Toughening in this connection is maximized by increases in microvoid growth strain that result from decreases in the size and areal number density of secondary void-initiating particles. Finally, the occurrence of transient instabilities during crack extension in short-rod specimens is explained with a phenomenological model that relates crack growth stability to natural variations in the dispersion of secondary void-initiating particles in the microstructure.

42 ENGINEERING↗
DOE OSTI2025

COMSOL Results for the Nominal Steady-State Operation of the Proposed 95-MW LEU Silicide Core for HFIR Conversion

Engineering design studies are being performed to determine the feasibility of converting the High Flux Isotope Reactor (HFIR) from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel at Oak Ridge National Laboratory. This activity is sponsored by the Office of Reactor Conversion and Uranium Supply (ORCUS) under the auspices of the US Department of Energy National Nuclear Security Administration’s Office of Material Management and Minimization. HFIR is a very high flux, pressurized, light water–cooled and moderated, flux trap–type research reactor with a core made of involute shaped U 3 O 8 /Al cermet fuel plates and coolant channels. HFIR currently operates at a thermal power of 85 MW and supports key national and international missions in neutron scattering, isotope production, materials/fuels irradiation, neutron activation analysis, gamma irradiation, and neutrino research. Advanced multiphysics computational fluid dynamics models have been developed in the COMSOL Multiphysics software to simulate the steady-state operating conditions for the proposed low-and high-density LEU U 3 Si 2 -Al (uranium silicide dispersion) fuel designs. The COMSOL models for HFIR inner and outer fuel element models incorporate various essential inputs and physics such as spatially dependent nuclear heat deposition, multilayer heat conduction, conjugate heat transfer, turbulent flows (using Reynolds-averaged Navier Stokes turbulence models), structural mechanics (thermal–structural interactions and fuel swelling), and oxide layer build-up. This report presents the best-estimate thermal hydraulics results for the low- and high-density optimized silicide LEU core designs at 95 MW steady-state nominal operation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

The Challenge of Observing Patchy Reionization with CMB Optical-Depth Fluctuations

Spatial fluctuations in the Thomson optical depth encode information about the inhomogeneous nature of cosmic reionization. We compute the optical-depth angular power spectrum, $C_\ell^{ττ}$, using past lightcones constructed from five Cosmic Reionization on Computers (CROC) radiation-hydrodynamical simulations. By decomposing the electron-density field into patchy and density components, we quantify the separate contributions of ionization-fraction and baryon-density fluctuations to the optical-depth anisotropy. Because the simulations end at $z\approx5$, we supplement the reionization-era signal with an analytic estimate of the fully ionized low-redshift contribution. We find that baryon-density fluctuations dominate the high-redshift signal over most angular scales, while the accumulated low-redshift contribution exceeds the high-redshift signal across the full multipole range considered. Our results demonstrate that a significant fraction of the optical-depth power is not uniquely associated with reionization morphology, implying that future interpretations of $C_\ell^{ττ}$ must account for the density contribution in addition to patchy ionization.

Takoudes, Nick [Chicago U., Astron. Astrophys. Ctr↗
DOE OSTI2020

Mechanical behavior and texture evolution of WE43 magnesium-rare earth alloy in Split-Hopkinson Pressure Bar and Taylor Impact Cylinder Testing

Mechanical behavior and texture evolution of Mg rare-earth alloy WE43 is investigated for strain-rates 10 -3 /s to upwards of 10 5 /s for the two material conditions - as-cast (AC) and T6 age hardened, rolled plate (RT6). The high strain-rate behavior is tested using both Taylor cylinder impact tests (TC) and split Hopkinson pressure bar tests (SHB) and bulk textures are obtained using neutron diffraction. Unlike the quasi-static strained material, AC and RT6 SHB retained high hardening rates throughout the test, even up to 30% true strain. Moreover, the high strain-rate data revealed that the RT6 material has a much higher strength than the AC material, but similar hardening rates despite significantly different initial texture. The flow stress near yield increased up to I 0% for RT6 and up to 30% for AC as the strain-rate increased six orders of magnitude from quasi-static rates 10 -3 /s to 10 3 /s. Neither material exhibited significant plastic anisotropy over the broad range of strain rates, despite the fact that the RT6 material had a moderately strong initial texture. In the TC tests, the geometric cross-sectional changes and texture along the cylinder from the cylindrical sample foot to head are measured and from the neutron diffraction texture analysis, upper-bound estimates of twin volume fraction are obtained as well as dislocation density from analyzing diffraction peak broadening. Recorded geometrical changes along several loading directions show that the material has deformed remarkably homogeneously under impact. Analysis of deformed textures indicates that {101¯2} extension deformation twinning occurred in the RT6 condition over the range of strain rates, with an upper bound estimate of 40% twin volume fraction for approximately 0.10-0.25 true strain. The peak texture components after the impact have their c-axes closely aligned with the impact direction. As a result, these observations are presented and rationalized in the paper.

42 ENGINEERING↗
DOE OSTI2021

A 3D Full Stress Tensor Model for Oklahoma

Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.

58 GEOSCIENCES↗
DOE OSTI2020

Transient solute drag and strain aging of dislocations

The transient drag force exerted by mobile solutes on a moving dislocation is computed using continuum theory. These mobile solutes form so-called Cottrell atmospheres around dislocations during static and dynamic strain aging. We evaluate the evolution of the drag force exerted by the atmosphere under two velocity time-histories: impulsive acceleration to a chosen velocity and a constant acceleration rate. A particular focus is on the conditions under which the stationary limit assumed by theories of dynamic strain aging is obeyed. According to our results, two conditions—one on the dislocation velocity and one on the acceleration rate—must be satisfied for the stationary limit to hold. Using the Orowan relation and a line tension model, we obtain estimates for the temperature, stress, strain rate, and dislocation density regimes where the stationary limit is valid, and compare these results with experiments for a few material systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In-situ two-dimensional temperature measurements using x-ray fluorescence spectroscopy in laminar flames with high silica particle concentrations

X-ray fluorescence spectroscopy (XRF) was used to measure temperatures and study mixing, for the first time in silica particle synthesis flames. Hexamethyldisiloxane (HMDSO) and trimethylsilanol (TMSO) were the particle precursors. A multi-element diffusion burner was used to produce a flat methane flame, and the precursors, dilute in inert gas, were injected via a central jet. Krypton was the fluorescent medium at 3.2 % concentration by volume. Scans with Kr in the central flow and not in the main flow were made to assess the mixing effects between the central and main gas flows. When HMDSO or TMSO were added, a secondary diffusion flame formed between the jet and the main methane flame. The results revealed a dramatic change in the centerline temperature profile of the jet gases when HMDSO or TMSO were added. The main methane flame stoichiometry also affected the temperature profiles. The results show HMDSO and TMSO reactions are initiated in a lowtemperature and low-oxygen concentration region of the jet where thermal decomposition is not expected to be significant. Reaction of the particle precursors is therefore attributed to radical transport from the main methane flame. In the current work, particle number densities of up to 270 g/m 3 locally are estimated. Thus, the study also demonstrates the capability of the XRF technique for high spatial fidelity measurements in flames with high concentrations of condensed-phase particles, leveraging the attribute that the XRF signal is generally not impacted by condensed-phase interferences. In conclusion, the observations and data obtained in this study inform likely reaction pathways for this important class of siloxane compounds.

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

Quantification of Shallow Defect Suppression in Aminosilane-Passivated Lead Halide Perovskites via Optical and Electrochemical Methods

We combine fluence-dependent photoluminescence (PL) and electrochemical characterization to estimate defect-dependent trapping rates, energetic distributions, and densities of defects, in lead halide perovskite films Cs 0.1 FA 0.9 PbI 3 and Cs 0.17 FA 0.83 Pb­(I 0.75 Br 0.25 ) 3 , before and after trap passivation with vapor-deposited (3-aminopropyl)­trimethoxysilane (APTMS). For both compositions, the PL studies show that energetically shallow subgap near-valence defect densities lead to bimolecular nonradiative recombination which is reduced by 92% after APTMS treatment. Electrochemical characterization of these same active layers shows for the first time that APTMS passivation impacts mainly on near-valence (mobile anionic) defects that decrease in density by 4 orders of magnitude, from 10 18 cm –3 to 10 14 cm –3 , as a result of silane modification. Here, the combination of PL and electrochemical characterization promises a unique approach to speciation and quantification of a broader distribution of trap states in perovskites and provides straightforward assessments of the efficacy of defect mitigation strategies.

Defects↗
DOE OSTI2021

Radially resolved active charge exchange measurements of the hydrogenic isotope fraction on DIII-D

Radially resolved hydrogenic isotope fraction measurement capabilities have been developed for DIII-D using the main-ion charge exchange recombination (MICER) spectroscopy system in preparation for mixed hydrogen and deuterium experiments. Constraints on the hydrogenic ion temperatures and velocities based on measurements of the impurity ion properties are required to accurately fit the spectrum. Corrections for cross sectional distortions, spatial smearing due to the halo, and a neoclassical offset between the impurity and hydrogenic toroidal rotation are applied to the constraints prior to fitting the MICER spectrum. Extensive atomic physics calculations have been performed using the FIDASIM code, which has recently been improved to allow simulations using mixtures of hydrogenic species. Overall, these results demonstrate that for the same plasma parameters, the Dα emission is 20%–30% brighter than Hα due to differences in rate coefficients associated with the different ion thermal velocities for the same temperature and therefore must be taken into consideration when calculating absolute densities. However, despite these differences, the absolute error when estimating the hydrogen isotope fraction [n H /(n H + n D )] by using the Hα radiance fraction [L Hα /(L Hα + L Dα )] is typically less than 5% due to the way the fraction is formed, making the radiance fraction a reasonably accurate estimate of the isotope fraction for most cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Probing anisotropies of the Stochastic Gravitational Wave Background with LISA

Here, we investigate the sensitivity of the Laser Interferometer Space Antenna (LISA) to the anisotropies of the Stochastic Gravitational Wave Background (SGWB). We first discuss the main astrophysical and cosmological sources of SGWB which are characterized by anisotropies in the GW energy density, and we build a Signal-to-Noise estimator to quantify the sensitivity of LISA to different multipoles. We then perform a Fisher matrix analysis of the prospects of detectability of anisotropic features with LISA for individual multipoles, focusing on a SGWB with a power-law frequency profile. We compute the noise angular spectrum taking into account the specific scan strategy of the LISA detector. We analyze the case of the kinematic dipole and quadrupole generated by Doppler boosting an isotropic SGWB. We find that β Ω GW ~ 2 × 10 -11 is required to observe a dipolar signal with LISA. The detector response to the quadrupole has a factor ~ 10 3 β relative to that of the dipole. The characterization of the anisotropies, both from a theoretical perspective and from a map-making point of view, allows us to extract information that can be used to understand the origin of the SGWB, and to discriminate among distinct superimposed SGWB sources.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Quantitative comparison of impurity transport in turbulence reduced and enhanced scenarios at Wendelstein 7-X

We assess the turbulent particle transport being responsible for the limitation of the confinement and, thus, the overall performance of the neoclassically optimized stellarator Wendelstein 7-X. The radial particle transport is experimentally inferred from the evaluation of impurity injection into turbulence reduced and enhanced plasma scenarios revealing a completely different confinement behavior. The impact of the density gradient on the turbulent ion transport is theoretically estimated using large-scale non-linear gyro-kinetic simulations enabling, for the first time in Wendelstein 7-X, a quantitative comparison to the experimentally assessed impurity transport properties. We demonstrate that impurity transport in most of the Wendelstein 7-X discharges, up to now impossible to cover only with neoclassical estimations, is dominated by turbulence and can be modelled via gyro-kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Training quantum neural networks using the quantum information bottleneck method

Abstract We provide in this paper a concrete method for training a quantum neural network to maximize the relevant information about a property that is transmitted through the network. This is significant because it gives an operationally well founded quantity to optimize when training autoencoders for problems where the inputs and outputs are fully quantum. We provide a rigorous algorithm for computing the value of the quantum information bottleneck quantity within error ε that requires O ( log 2 ⁡ ( 1 / ϵ ) + 1 / δ 2 ) queries to a purification of the input density operator if its spectrum is supported on { 0 } ⋃ [ δ , 1 − δ ] for δ > 0 and the kernels of the relevant density matrices are disjoint. We further provide algorithms for estimating the derivatives of the QIB function, showing that quantum neural networks can be trained efficiently using the QIB quantity given that the number of gradient steps required is polynomial.

Çatlı, Ahmet Burak (ORCID:0000000152294141)↗
DOE OSTI2022

Cold gas mass measurements for the era of large optical spectroscopic surveys

ABSTRACT Gas plays an important role in many processes in galaxy formation and evolution, but quantifying the importance of gas has been hindered by the challenge to measure gas masses for large samples of galaxies. Data sets of direct atomic and molecular gas measurements are sufficient to establish simple scaling relations, but often not large enough to quantify three-parameter relations, or second-order dependences. As an alternative approach, we derive here indirect cold gas measurements from optical emission lines using photoionization models for galaxies in the Sloan Digital Sky Survey (SDSS) main galaxy sample and the PHANGS-MUSE survey. We calibrate the gas surface density measurements using xCOLD GASS and PHANGS-ALMA molecular gas measurements to ensure that our measurements are reliable. We demonstrate the importance of taking into account the scale dependence of the relation between optical depth (τ V ) and gas surface density (Σ gas ) and provide a general prescription to estimate Σ gas from τ V , metallicity, and the dust-to-metal ratio, at any arbitrary physical resolution. To demonstrate that the indirect cold gas masses are accurate enough to quantify the role of gas in galaxy evolution, we study the mass–metallicity relation of SDSS galaxies and show that as a third parameter gas mass is better than star formation rate at reducing the scatter of the relation, as predicted by models and simulations.

79 ASTRONOMY AND ASTROPHYSICS↗