Search NASA⌕ Search

SEARCH · Search NASA

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 505 records · Page 28

Thinking outside the box: effects of modes larger than the survey on matter power spectrum covariance

Accurate power spectrum (or correlation function) covariance matrices are a crucial requirement for cosmological parameter estimation from large scale structure surveys. In order to minimize reliance on computationally expensive mock catalogs, it is important to have a solid analytic understanding of the different components that make up a covariance matrix. Considering the matter power spectrum covariance matrix, it has recently been found that there is a potentially dominant effect on mildly non-linear scales due to power in modes of size equal to and larger than the survey volume. This beat coupling effect has been derived analytically in perturbation theory and while it has been tested with simulations, some questions remain unanswered. Moreover, there is an additional effect of these large modes, which has so far not been included in analytic studies, namely the effect on the estimated average density which enters the power spectrum estimate. In this article, we work out analytic, perturbation theory based expressions including both the beat coupling and this local average effect and we show that while, when isolated, beat coupling indeed causes large excess covariance in agreement with the literature, in a realistic scenario this is compensated almost entirely by the local average effect, leaving only ∼ 10% of the excess. We test our analytic expressions by comparison to a suite of large N-body simulations, using both full simulation boxes and subboxes thereof to study cases without beat coupling, with beat coupling and with both beat coupling and the local average effect. For the variances, we find excellent agreement with the analytic expressions for k < 0.2 hMpc{sup −1} at z = 0.5, while the correlation coefficients agree to beyond k = 0.4 hMpc{sup −1}. As expected, the range of agreement increases towards higher redshift and decreases slightly towards z = 0. We finish by including the large-mode effects in a full covariance matrix description for arbitrary survey geometry and confirming its validity using simulations. This may be useful as a stepping stone towards building an actual galaxy (or other tracer's) power spectrum covariance matrix.

79 ASTRONOMY AND ASTROPHYSICS↗

DETERMINATION OF TRANSVERSE DENSITY STRUCTURING FROM PROPAGATING MAGNETOHYDRODYNAMIC WAVES IN THE SOLAR ATMOSPHERE

We present a Bayesian seismology inversion technique for propagating magnetohydrodynamic transverse waves observed in coronal waveguides. The technique uses theoretical predictions for the spatial damping of propagating kink waves in transversely inhomogeneous coronal waveguides. It combines wave amplitude damping length scales along the waveguide with theoretical results for resonantly damped propagating kink waves to infer the plasma density variation across the oscillating structures. Provided that the spatial dependence of the velocity amplitude along the propagation direction is measured and the existence of two different damping regimes is identified, the technique would enable us to fully constrain the transverse density structuring, providing estimates for the density contrast and its transverse inhomogeneity length scale.

79 ASTRONOMY AND ASTROPHYSICS↗

Variable Blue Straggler Stars in the Open Cluster NGC 6819 Observed in the Kepler “Superstamp” Field

Abstract NGC 6819 is an open cluster of age 2.4 Gyr that was in the NASA Kepler spacecraft’s field of view from 2009 to 2013. The central part of the cluster was observed in a 200 × 200 pixel “superstamp” during these four years in 30 minute cadence photometry, providing a unique, high-precision, long time-series data set. The cluster contains “blue straggler” stars, i.e., stars on the main sequence above the cluster turnoff that should have left the main sequence to become red giants. We present light curves and pulsation frequency analyses derived from custom photometric reductions for five confirmed cluster members—four blue stragglers and one star near the main-sequence turnoff. Two of these stars show a rich spectrum of δ Scuti pulsation modes, with 236 and 124 significant frequencies identified, respectively, while two stars show mainly low-frequency modes, characteristic of γ Doradus variable stars. The fifth star, a known active X-ray binary, shows only several harmonics of two main frequencies. For the two δ Scuti stars, we use a frequency separation–mean density relation to estimate their mean densities, and then use these values along with their effective temperature to derive their stellar masses and radii. For the two stars showing low frequencies, we searched for period-spacing sequences that may be representative of gravity-mode or Rossby-mode sequences, but found no clear sequences. The common age for the cluster members, considered along with the frequencies, will provide valuable constraints for asteroseismic analyses, and may shed light on the origin of the blue stragglers.

79 ASTRONOMY AND ASTROPHYSICS↗

An interstellar line coincident with the P(2,1) transition of hydronium (H3O+)

Double-sideband observations of the Orion KL nebular region with the NRAO 12 m radiotelescope in March 1986 detected a weak emission spectral line that coincided with the P(2,1) rotation-inversion transition of H3O+. Spectral maps are provided for the observed line at 307,192 MHz with 2 and 1 MHz resolution. The line parameters, along with a calculated transition dipole moment of 1.44 D, are used to estimate column densities of 7.7 and 8.7 x 10 to the 13/sq cm for excitation temperatures of 50 and 90 K, respectively. A minimum H2 density of 1 million/cu cm is estimated for the region if the thermalization of the transition to the excitation temperatures is to occur. Finally, fractional abundances of H2O in the Orion cloud are projected to be 1/100,000 and 1/1 million for the nonthermal and thermal transition cases, respectively.

Hollis, J. M.↗

Spatio-temporal and weather characterization of road loads of electrified heavy-duty commercial vehicles across U.S. interstate roads

Adoption of battery electric vehicles (BEV) in heavy duty (HD) commercial freight transportation is difficult due to technological and economic hurdles. Beyond safety and compliance, fleet and operational logistics necessitate both high uptime and parity with diesel system productivity/Total Cost of Ownership to support widespread deployment of electric powertrains. However, relatively high energy storage costs, along with the higher weight of BEV systems, limit the viability of HD commercial freight transport to shorter-range applications where smaller batteries will serve for mission energy requirements (single operational shift). Knowing the energy consumption and operating variations of these commercial vehicle systems is crucial for effectively sizing the energy storage systems. This paper is the first in a series of studies to understand the regional specific operating design domain variations of commercial Class 8 HD trucks and the associated impact to their energy requirements. In particular, the local weather conditions are shown to influence the total vehicle energy usage. Further, the impact of temperature, pressure, and humidity changes are shown to impact the local air density. This is needed to calculate aerodynamic drag on vehicles and has been found to play a significant role in the overall performance of a vehicle. Although the methodology of calculating air density varies only slightly throughout the literature, the application to transportation has still been limited. This study provides a means by which air density can be estimated for the contiguous U.S. using the NOAA MADIS dataset. These air density estimates were then used to determine vehicle performance in varying regions of the country to highlight the importance of consideration of weather variables when monitoring vehicle performance, but also to provide recommendations based on these locales upon fleet conversion to BEVs.

Moore, Amy↗

A NEW RECIPE FOR OBTAINING CENTRAL VOLUME DENSITIES OF PRESTELLAR CORES FROM SIZE MEASUREMENTS

We propose a simple analytical method for estimating the central volume density of prestellar molecular cloud cores from their column density profiles. Prestellar cores feature a flat central part of the column density and volume density profiles of the same size indicating the existence of a uniform-density inner region. The size of this region is set by the thermal pressure force which depends only on the central volume density and temperature of the core, and can provide a direct measurement of the central volume density. Thus, a simple length measurement can immediately yield a central density estimate independent of any dynamical model for the core and without the need for fitting. Using the radius at which the column density is 90% of the central value as an estimate of the size of the flat inner part of the column density profile yields an estimate of the central volume density within a factor of two for well-resolved cores.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal Estimation Inversion of Ionospheric Electron Density from GNSS-POD Limb Measurements: Part I-Algorithm and Morphology

GNSS-LEO radio links from Precise Orbital Determination (POD) and Radio Occultation (RO) antennas have been used increasingly in characterizing the global 3D distribution and variability of ionospheric electron density (N e ). In this study, we developed an optimal estimation (OE) method to retrieve N e profiles from the slant total electron content (hTEC) measurements acquired by the GNSS-POD links at negative elevation angles (ε < 0°). Although both OE and onion-peeling (OP) methods use the Abel weighting function in the N e inversion, they are significantly different in terms of performance in the lower ionosphere. The new OE results can overcome the large N e oscillations, sometimes negative values, seen in the OP retrievals in the E-region ionosphere. In the companion paper in this Special Issue, 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 N e morphology. The unprecedented spatiotemporal sampling of the ionosphere from these constellations now allows a detailed analysis of the frequency–wavenumber spectra for the N e 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-existence of eastward- and westward-propagating wave4 components implies the presence of a stationary wave4 (SPW4), as suggested by other earlier studies. Further improvements to the OE method are proposed, including a tomographic inversion technique that leverages the asymmetric sampling about the tangent point associated with GNSS-LEO links.

GNSS constellation↗

Effects of levitated dust on astronomical observations from the lunar surface

It is believed that a substantial population of levitated dust is present in the terminator region of the moon. Stray light scattered by this dust layer may contaminate astronomical observations made from the lunar surface using infrared, visible, and ultraviolet light. The evidence for dust levitation stems from: Surveyor vidicon images of horizon glow; anomalous brightness in photographs of the solar corona taken by Apollo astronauts while the spacecraft was just inside the moon's shadow; and observations by Apollo astronauts of streamers just prior to lunar orbital sunrise or just after lunar orbital sunset. It has been proposed that the differential charging of the lunar surface in the terminator region due to photoemission and the consequent strong local electric fields comprise the mechanism responsible for this levitation. Although quantitative data on the levitated lunar dust distribution are meager, it is possible to estimate column densities and sizes. In this paper we summarize the estimates of particulate sizes and number densities of previous authors, and construct a nominal terminator dust distribution, as a function of particulate radius and altitude above the lunar surface. Using the model we estimate the brightness of scattered sunshine for three wavelength bands. For the results in the visible wavelengths, we compare the estimated brightness with the known brightness of selected astronomical objects and discuss the implications for lunar-based astronomy.

Murphy, D. L.↗

Modeled Salt Density for Nuclear Material Estimation in the Treatment of Spent Nuclear Fuel

Spent metallic nuclear fuel is being treated in a pyrometallurgical process that includes electrorefining the uranium metal in molten eutectic LiCl-KCl as the supporting electrolyte. We report a model for determining the density of the molten salt. Inventory operations account for the net mass of salt and for the mass of actinides present. It was necessary to know the molten salt density but difficult to measure, and it was decided to model the salt density for the initial treatment operations. The model assumes that volumes are additive for the ideal molten salt solution as a starting point; subsequently a correction factor for the lanthanides and actinides was developed. After applying the correction factor, the percent difference between the net salt mass in the electrorefiner and the resulting modeled salt mass decreased from more than 4.0% to approximately 0.1%. As a result, there is no need to measure the salt density at 500 C for inventory operations; the model for the salt density is found to be accurate.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Parameter estimation with maximal updated densities

A recently developed measure-theoretic framework solves a stochastic inverse problem (SIP) for models where uncertainties in model output data are predominantly due to aleatoric (i.e., irreducible) uncertainties in model inputs (i.e., parameters). The subsequent inferential target is a distribution on parameters. Another type of inverse problem is to quantify uncertainties in estimates of “true” parameter values under the assumption that such uncertainties should be reduced as more data are incorporated into the problem, i.e., the uncertainty is considered epistemic. A major contribution of this work is the formulation and solution of such a parameter identification problem (PIP) within the measure-theoretic framework developed for the SIP. The approach is novel in that it utilizes a solution to a stochastic forward problem (SFP) to update an initial density only in the parameter directions informed by the model output data. In other words, this method performs “selective regularization” only in the parameter directions not informed by data. Further, the solution is defined by a maximal updated density (MUD) point where the updated density defines the measure-theoretic solution to the PIP. Another significant contribution of this work is the full theory of existence and uniqueness of MUD points for linear maps with Gaussian distributions. Data-constructed Quantity of Interest (QoI) maps are also presented and analyzed for solving the PIP within this measure-theoretic framework as a means of reducing uncertainties in the MUD estimate. We conclude with a demonstration of the general applicability of the method on two problems involving either spatial or temporal data for estimating uncertain model parameters. The first problem utilizes spatial data from a stationary partial differential equation to produce a MUD estimate of an uncertain boundary condition. The second problem utilizes temporal data obtained from the state-of-the-art ADvanced CIRCulation (ADCIRC) model to obtain a MUD estimate of uncertain wind drag coefficients for a simulated extreme weather event near the Shinnecock Inlet located in the Outer Barrier of Long Island, NY, USA.

42 ENGINEERING↗

Detection of Missing Assemblies and Estimation of the Scattering Densities in a VSC-24 Dry Storage Cask with Cosmic-Ray-Muon-Based Computed Tomography

We report that highly energetic, cosmic-ray muons can penetrate a dry storage cask and yield information about the material inside it by making use of the physics of multiple Coulomb scattering. Work by others has shown this information may be used for verification of dry storage cask contents after continuity of knowledge has been lost. In our modeling and simulation approach, we use ideal planar radiation detectors to record the trajectories and momentum of both incident and exiting cosmic ray muons; this choice allows us to demonstrate the fundamental limit of the technology for a particular measurement and reconstruction method. In a method analogous to computed tomography with the attenuation coefficient replaced by scattering density, we apply a filtered back projection algorithm in order to reconstruct the geometry in modeled scenarios for a VSC-24 concrete-walled cask. We also report on our attempt to estimate material-specific information. A scenario where one of the middle four spent nuclear fuel assemblies is missing—undetectable with a simple PoCA-based approach—is expected to be detectable with a CT-based approach. Moreover, a trickier scenario where one or more assemblies is replaced by a dummy assembly is put forward. Lastly, in this case, we expect that this dry storage cask should be found to be not as declared based on our simulation and reconstruction results.

36 MATERIALS SCIENCE↗

System-independent characterization of materials using dual-energy computed tomography

A system for characterizing the material of an object scanned via a dual-energy computed tomography scanner is provided. The system generates photoelectric and Compton sinograms based on a photoelectric-Compton decomposition of low-energy and high-energy sinograms generated from the scan and based on a scanner spectral response model. The system generates a Compton volume with Compton attenuation coefficients from the Compton sinogram and a photoelectric volume with photoelectric attenuation coefficients from the photoelectric sinogram. The system generates an estimated effective atomic number for a voxel and an estimated electron density for the voxel from the Compton attenuation coefficient and photoelectric coefficient for the voxel and scanner-specific parameters. The system then characterizes the material within the voxel based on the estimated effective atomic number and estimated electron density for the voxel. This information can be used to provide a mapping of known effective atomic numbers and known electron densities to known materials.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Hydrometeor Size Distribution Measurements by Imaging the Attenuation of a Laser Spot

The optical extinction of a laser due to scattering of particles is a well-known phenomenon. In a laboratory environment, this physical principle is known as the Beer-Lambert law, and is often used to measure the concentration of scattering particles in a fluid or gas. This method has been experimentally shown to be a usable means to measure the dust density from a rocket plume interaction with the lunar surface. Using the same principles and experimental arrangement, this technique can be applied to hydrometeor size distributions, and for launch-pad operations, specifically as a passive hail detection and measurement system. Calibration of a hail monitoring system is a difficult process. In the past, it has required comparison to another means of measuring hydrometeor size and density. Using a technique recently developed for estimating the density of surface dust dispersed during a rocket landing, measuring the extinction of a laser passing through hail (or dust in the rocket case) yields an estimate of the second moment of the particle cloud, and hydrometeor size distribution in the terrestrial meteorological case. With the exception of disdrometers, instruments that measure rain and hail fall make indirect measurements of the drop-size distribution. Instruments that scatter microwaves off of hydrometeors, such as the WSR-88D (Weather Surveillance Radar 88 Doppler), vertical wind profilers, and microwave disdrometers, measure the sixth moment of the drop size distribution (DSD). By projecting a laser onto a target, changes in brightness of the laser spot against the target background during rain and hail yield a measurement of the DSD's second moment by way of the Beer-Lambert law. In order to detect the laser attenuation within the 8-bit resolution of most camera image arrays, a minimum path length is required. Depending on the intensity of the hail fall rate for moderate to heavy rainfall, a laser path length of 100 m is sufficient to measure variations in optical extinction using a digital camera. For hail fall only, the laser path may be shorter because of greater scattering due to the properties of hailstones versus raindrops. A photodetector may replace the camera in automated installations. Laser-based rain and hail measurement systems are available, but they are based on measuring the interruption of a thin laser beam, thus counting individual hydrometeors. These systems are true disdrometers since they also measure size and velocity. The method reported here is a simple method, requiring far less processing, but it is not a disdrometer.

Lane, John↗

Using X-ray spectrum imaging to quantify fission gas in high-burnup UO2 fuel

Analysis of fission gas bubbles (FGBs) in light-water reactor (LWR) fuel is needed to improve the understanding and predictive capability of fuel evolution under normal- and transient conditions. One of the most important parameters of a FGB is the pressure of the gas, primarily Xe, inside it. Bubble volume and location (inter- vs. intragranular) are important considerations as well. However, experimental analysis of such bubbles is challenging due to their small size and embedded nature, and usefulness of the data requires large numbers be analyzed. This paper proposes a method to measure the pressure of Xe bubbles using X-ray spectrum imaging (XSI) in scanning transmission electron microscopy (STEM). From X-ray generation and instrumental parameters, the number of Xe atoms yielding a given number of Xe L-series X-ray counts is estimated, and then from the estimated volume of the bubble the Xe density is estimated, which can then be converted to pressure via an equation-of-state. We apply the method to XSI datasets from high burnup (HBu) fuel from North Anna 1 reactor in as-irradiated and post-LOCA-test conditions and find Xe pressures in the FGBs clustered around 1 GPa.

Parish, Chad [ORNL] (ORCID:0000000312097439)↗

Upper limits to (O) in the lower thermosphere from airglow

A method for estimating the upper limit to oxygen atom density in the lower thermosphere is presented. The method is based on an analytical equation for the energy density of recombination process, and the radiation of O2 atoms at excited levels. As an example of the method, upper limits of oxygen atom density are estimated for the CIRA (O2) reference atmosphere of 1972. The estimates are given in a table.

Sharp, W. E.↗

Galactic H I and the interstellar medium in Ursa Major

Sensitive neutral hydrogen measurements have been made over a 100 square degree field in Ursa Major. The observations cover the area which has the least neutral hydrogen and probably the least interstellar matter of any direction in the sky. It is likely that this data set contains every line of sight that intersects a total H I column density of less than 6.5 x 10 to the 19th/sq cm. The smallest colum density is estimated to be 4.4 + or - 0.5 x 10 to the 19th/sq cm. Observations at 10 arcmin angular resolution reveal little structure that is not seen at 21 arcmin resolution, and it is believed that the estimates of the total H I column density are accurate on yet smaller angular scales. The lowest N(H I) area is also extremely deficient in H-alpha emission, and shows no evidence for significant amounts of H2.

Jahoda, Keith↗

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

Analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed.

transient↗

Investigation of Transient Gas Phase Column Density Due to Droplet Evaporation

Analytical expressions are developed for estimating column density near a rapidly evaporating droplet along general paths. The influence for instantaneous evaporation is created first as a limiting case, where the peak value occurs at the time it takes a wave of vapor to reach the closest point along the optical path traveling at its most probable thermal speed. Next the case for finite-period evaporation is evaluated for constant conditions. Compared to the instantaneous case, peak column density occurs shortly after droplet extinction at a lower intensity. A new mathematical function is discovered that solves the integrals associated with this case. Finally, ways to account for droplet motion and changes in evaporation rate with size and temperature are discussed.

Column Density↗