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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,531 records · Page 85

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Evidence for AGN-regulated Cooling in Clusters at z ∼ 1.4: A Multiwavelength View of SPT-CL J0607-4448

Abstract We present a multiwavelength analysis of the galaxy cluster SPT-CL J0607-4448 (SPT0607), which is one of the most distant clusters discovered by the South Pole Telescope at z = 1.4010 ± 0.0028. The high-redshift cluster shows clear signs of being relaxed with well-regulated feedback from the active galactic nucleus (AGN) in the brightest cluster galaxy (BCG). Using Chandra X-ray data, we construct thermodynamic profiles and determine the properties of the intracluster medium. The cool-core nature of the cluster is supported by a centrally peaked density profile and low central entropy ( K 0 = 18 − 9 + 11 keV cm 2 ), which we estimate assuming an isothermal temperature profile due to the limited spectral information given the distance to the cluster. Using the density profile and gas cooling time inferred from the X-ray data, we find a mass-cooling rate M ̇ cool = 100 − 60 + 90 M ⊙ yr −1 . From optical spectroscopy and photometry around the [O ii ] emission line, we estimate that the BCG star formation rate is SFR [ O II ] = 1.7 − 0.6 + 1.0 M ⊙ yr −1 , roughly two orders of magnitude lower than the predicted mass-cooling rate. In addition, using ATCA radio data at 2.1 GHz, we measure a radio jet power P cav = 3.2 − 1.3 + 2.1 × 10 44 erg s −1 , which is consistent with the X-ray cooling luminosity ( L cool = 1.9 − 0.5 + 0.2 × 10 44 erg s −1 within r cool = 43 kpc). These findings suggest that SPT0607 is a relaxed, cool-core cluster with AGN-regulated cooling at an epoch shortly after cluster formation, implying that the balance between cooling and feedback can be reached quickly. We discuss the implications for these findings on the evolution of AGN feedback in galaxy clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1994

Hot gas and the origin of the nebular continuum in novae

A featureless blue continuum (with constant f(sub nu)) is a defining feature of declining classical novae. The fact that f(sub nu) is constant into the infrared, and the absence of a Balmer jump, suggests that this continuum originates in hot tenous gas. The electron temperature and density of the hot gas phase of classical nova QU Vul 1984 are estimated from ground-based optical and IUE ultraviolet observations. This region has a temperature of 10(exp 5) - 10(exp 6), and a density in the range of 2 x 10(exp 4) to 3 x 10(exp 5)/cu cm. In addition, a colder, denser shell composed of many 'clumps' or condensations is present, and it is responsible for most of the optical-to-ultraviolet emission lines. The hot has, however, may be much more massive and could dominate the energetics of the ejecta. The total returned mass could be as high as 10(exp -4) and 4 x 10(exp -3) solar mass, depending on conditions in the hot gas. The large mass of the shell, in combination with the high metallicity found in this nova, supports the idea that novae may be an important contributor to the chemical enrichment of the interstellar medium.

Saizar, Pedro↗

NASA NTRS2007

Using MODIS and GLAS Data to Develop Timber Volume Estimates in Central Siberia

The boreal forest is the Earth's largest terrestrial biome, covering some 12 million km2 and accounting for about one third of this planet's total forest area. Mapping of boreal forest's type, structure parameters and biomass are critical for understanding the boreal forest's significance in the carbon cycle, its response to and impact on global climate change. Ground based forest inventories, have much uncertainty in the inventory data, particularly in remote areas of Siberia where sampling is sparse and/or lacking. In addition, many of the forest inventories that do exist for Siberia are now a decade or more old. Thus, available forest inventories fail to capture the current conditions. Changes in forest structure in a particular forest-type and region can change significantly due to changing environment conditions, and natural and anthropogenic disturbance. Remote sensing methods can potentially overcome these problems. Multispectral sensors can be used to provide vegetation cover maps that show a timely and accurate geographic distribution of vegetation types rather than decade old ground based maps. Lidar sensors can be used to directly obtain measurements that can be used to derive critical forest structure information (e.g., height, density, and volume). These in turn can used to estimate biomass components using allometric equations without having to use out dated forest inventory. Finally, remote sensing data is ideally suited to provide a sampling basis for a rigorous statistical estimate of the variance and error bound on forest structure measures. In this study, new remote sensing methods were applied to develop estimates timber volume using NASA's MODerate resolution Imaging Spectroradiometer (MODIS) and unique waveform data of the geoscience laser altimeter system (GLAS) for a 10 deg x 10 deg area in central Siberia. Using MODIS and GLAS data, maps were produced for cover type and timber volume for 2003, and a realistic variance (error bound) for timber volume was calculated for the study area. In this 'study we used only GLAS footprints that had a slope value of less than 10 deg. This was done to avoid large errors due to the effect of slope on the GLAS models. The method requires the integration of new remote sensing methods with available ground studies of forest timber volume conducted in Russian forests. The results were compared to traditional ground forest inventory methods reported in the literature and to ground truth collected in the study area.

Ranson, K. Jon↗

DOE OSTI2020

The effects of earth model uncertainty on the inversion of seismic data for seismic source functions

SUMMARY We use Monte Carlo simulations to explore the effects of earth model uncertainty on the estimation of the seismic source time functions that correspond to the six independent components of the point source seismic moment tensor. Specifically, we invert synthetic data using Green’s functions estimated from a suite of earth models that contain stochastic density and seismic wave-speed heterogeneities. We find that the primary effect of earth model uncertainty on the data is that the amplitude of the first-arriving seismic energy is reduced, and that this amplitude reduction is proportional to the magnitude of the stochastic heterogeneities. Also, we find that the amplitude of the estimated seismic source functions can be under- or overestimated, depending on the stochastic earth model used to create the data. This effect is totally unpredictable, meaning that uncertainty in the earth model can lead to unpredictable biases in the amplitude of the estimated seismic source functions.

58 GEOSCIENCES↗

NASA NTRS1994

Regional and seasonal estimates of fractional storm coverage based on station precipitation observations

Simulated climates using numerical atmospheric general circulation models (GCMs) have been shown to be highly sensitive to the fraction of GCM grid area assumed to be wetted during rain events. The model hydrologic cycle and land-surface water and energy balance are influenced by the parameter bar-kappa, which is the dimensionless fractional wetted area for GCM grids. Hourly precipitation records for over 1700 precipitation stations within the contiguous United States are used to obtain observation-based estimates of fractional wetting that exhibit regional and seasonal variations. The spatial parameter bar-kappa is estimated from the temporal raingauge data using conditional probability relations. Monthly bar-kappa values are estimated for rectangular grid areas over the contiguous United States as defined by the Goddard Institute for Space Studies 4 deg x 5 deg GCM. A bias in the estimates is evident due to the unavoidably sparse raingauge network density, which causes some storms to go undetected by the network. This bias is corrected by deriving the probability of a storm escaping detection by the network. A Monte Carlo simulation study is also conducted that consists of synthetically generated storm arrivals over an artificial grid area. It is used to confirm the bar-kappa estimation procedure and to test the nature of the bias and its correction. These monthly fractional wetting estimates, based on the analysis of station precipitation data, provide an observational basis for assigning the influential parameter bar-kappa in GCM land-surface hydrology parameterizations.

Gong, Gavin↗

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↗

NASA NTRSDate not supplied

Battery Cell-to-Pack Scaling Trends for Electric Aircraft

Battery pack gravimetric energy density is one of the most important, yet often miss- estimated design parameters for sizing all-electric aircraft. Proper accounting for thermal, structural, and operational safety margins are frequently lost when extrapolating performance from the cell level to the aircraft level. This paper summarizes the relevant engineering and certification details needed to better account for the penalties associated when assembling battery packs. The relationship between the cell and pack energy density is not linear, as is often assumed. Furthermore, the relationship varies depending on pack requirements, cell chemistry, and architecture. Parametric, high-fidelity models are used to determine optimal battery pack sizes over a range of conditions to better quantify technology scaling effects.

Battery Electric Aircraft↗

NASA NTRS1994

The study of coronal plasma structures and fluctuations with Faraday rotation measurements

We report dual-frequency, polarimetric measurements of Faraday rotation of extragalactic radio sources viewed through the solar corona. The observations were made at the Very Large Array in 1990 during solar maximum. Of the nine observed, an excess rotation measure of -12.6 rad/sq m was detected for one source (0010+005), which was observed at an elongation of about 9 solar radii. This measurement is in fair agreement with an a priori model rotation measure of -8.6 rad/sq m estimated from coronal potential field models and the electron density model of Paetzold et al. (1992). Our measurement provides a value for the coronal magnetic field strength at 9 solar radii given a knowledge of the magnetic field sector structure, of 12.5 +/- 2.3 mG. Rotation measurements of 0010+005 were made approximately once per hour over an 11 hr period. During this interval, a slow change of about 1 rad/sq m/hr in rotation measure was detected. Although we are not absolutely certain that this drift is not unremoved ionospheric Faraday rotation, extensive analysis of data from the other sources suggests that this is not the case (Sakurai & Spangler 1994). The very long timescale for this variation argues against the agency of magnetohydrodynamics (MHD) waves, and we suggest occultation of 0010+005 by relatively static plasma structures in the corona. We filtered our rotation measure time series to search for variations on an hourly timescale, such as those reported by Hollweg et al. (1992), which could be attributed to coronal MHD waves. We were unable to detect such fluctuations and can report only an upper limit to the rms variation of 1.6 rad/sq m. This upper limit is of the same order, but slightly larger than the values typically reported by Hollweg et al. (1982). This upper limit to the rotation measure fluctuations limits the dimensionless wave amplitude (delta B)/B in the corona to be less than 0.7. Using the number, we estimate the MHD wave flux at the coronal base to be less than 1.6 x 10(exp 5) ergs/sq cm/s. This is less than the amount of wave energy flux required by wave-driven models of the solar wind. Finally, we discuss a number of ways in which such observations could be improved in the future.

Sakurai, Takayuki↗

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↗

NASA NTRS1976

Lossy radial diffusion of relativistic Jovian electrons

The radial diffusion equation with synchrotron losses is solved by the Laplace-transform method for near equatorially mirroring relativistic electrons. The evolution of a power-law distribution function is found, and the characteristics of synchrotron burnoff are stated in terms of explicit parameters for an arbitrary diffusion coefficient of a specific form. The peaking of the 10.4-cm volume emissivity from Jupiter at an L shell of about 1.8 provides an estimate of the diffusion coefficient in the radiation belts; one value is suggested as the appropriate modification, for an equatorial field strength of 4.2 G, of the Birmingham et al. (1974) result. Nonsynchrotron losses are included phenomenologically; from the phase-space densities reported by McIlwain and Fillius (1975), the particle lifetime is estimated. Asymptotic forms for the distribution in the strong synchrotron loss regime are provided.

Barbosa, D. D.↗

NASA NTRS1992

Novel multireceiver communication systems configurations based on optimal estimation theory

A novel multireceiver configuration for carrier arraying and/or signal arraying is presented. The proposed configuration is obtained by formulating the carrier and/or signal arraying problem as an optimal estimation problem, and it consists of two stages. The first stage optimally estimates various phase processes received at different receivers with coupled phase-locked loops wherein the individual loops acquire and track their respective receivers' phase processes but are aided by each other in an optimal manner via LF error signals. The proposed configuration results in the minimization of the the effective radio loss at the combiner output, and thus maximization of energy per bit to noise power spectral density ratio is achieved. A novel adaptive algorithm for the estimator of the signal model parameters when these are not known a priori is also presented.

Kumar, Rajendra↗

NASA NTRS1980

Geoid anomalies in the vicinity of subduction zones

The regional geoid of the southwest Pacific is matched reasonably well by results from a model of the upper mantle density structure (including slabs) associated with subduction zones of the region. Estimates of the geoid are obtained from Geos-3 and Seasat radar altimeter data. These data are very well suited to the task of detecting intermediate wavelength (600-4000 km) geopotential variations. Actually, subducting slabs can be expected to produce primarily intermediate and longer wavelength variations. Gravimetric profiles across trench/island arc complexes resolve primarily short wavelengths. The model represents subducting slabs as thin surfaces of anomalous mass per unit area. These surfaces are positioned using published seismicity results which detail the configuration of the Benioff zones. Crustal effects are ignored. Effects due to the contrast between the young thermal lithosphere of the behind-arc regions (marginal basins) and the older lithosphere seaward of the trench are modelled. Results indicate that the New Hebrides slab possesses an average areal density anomaly of about 300,000 gm/sq cm. This is about three times that which is estimated for the Tonga-Kermadec slab. Additional modelling suggests that slabs worldwide may be an important source of large, long wavelength gravity highs; i.e., they may contribute substantially to geopotential power of harmonic degree as low as three or four up to twenty or more.

Mcadoo, D. C.↗

NASA NTRS1983

The spatial distribution of fixed mutations within genes coding for proteins

An examination has been conducted of the extensive amino acid sequence data now available for five protein families - the alpha crystallin A chain, myoglobin, alpha and beta hemoglobin, and the cytochromes c - with the goal of estimating the true spatial distribution of base substitutions within genes that code for proteins. In every case the commonly used Poisson density failed to even approximate the experimental pattern of base substitution. For the 87 species of beta hemoglobin examined, for example, the probability that the observed results were from a Poisson process was the minuscule 10 to the -44th. Analogous results were obtained for the other functional families. All the data were reasonably, but not perfectly, described by the negative binomial density. In particular, most of the data were described by one of the very simple limiting forms of this density, the geometric density. The implications of this for evolutionary inference are discussed. It is evident that most estimates of total base substitutions between genes are badly in need of revision.

Holmquist, R.↗

NASA NTRSDate not supplied

Global Ionospheric Electron Density from GNSS-POD Limb Measurements

GNSS-LEO radio links from Precise Orbital Determination (POD) and Radio Occultation (RO) antennas have been used increasingly in studying and monitoring the global ionospheric electron density (Ne). In this study we developed an optimal estimation (OE) method to retrieve Ne profiles from the slant total electron content (hTEC) measurements acquired by the GNSS-POD links at negative elevation angles. The hmF2 and NmF2 from the OE retrieval are validated against ground-based ionosondes and radar observations, showing generally good agreements in NmF2 from all sites. Nighttime hmF2 measurements tend to agree better than the daytime when the ionosonde heights tend to be slightly lower. The OE algorithm has been applied to all GNSS-POD data acquired from the COSMIC-1 (2006-2019), COSMIC-2 (2019-present), and Spire (2019-present) constellations, showing a consistent ionospheric Ne morphology. A detailed analysis of the frequency-wavenumber spectra is made for the Ne variability at different heights. In the lower ionosphere (~150 km) we found significant spectral power in DE1, DW6, DW4, SW5, and SE4 wave components, in addition to well-known DW1, SW2 and DE3 waves. In the upper ionosphere (~450 km), additional wave components are still present, including DE4, DW4, DW6, SE4 and SW4. The co-exist of eastward and westward propagating WN4 components implies the presence of a stationary planetary wave4 (sPW4), as suggested by other earlier studies.

GNSS-POD↗

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↗