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 433 records · Page 24

Estimation and Modeling of Enceladus Plume Jet Density Using Reaction Wheel Control Data

The Cassini spacecraft was launched on October 15, 1997 by a Titan 4B launch vehicle. After an interplanetary cruise of almost seven years, it arrived at Saturn on June 30, 2004. In 2005, Cassini completed three flybys of Enceladus, a small, icy satellite of Saturn. Observations made during these flybys confirmed the existence of a water vapor plume in the south polar region of Enceladus. Five additional low-altitude flybys of Enceladus were successfully executed in 2008-9 to better characterize these watery plumes. The first of these flybys was the 50-km Enceladus-3 (E3) flyby executed on March 12, 2008. During the E3 flyby, the spacecraft attitude was controlled by a set of three reaction wheels. During the flyby, multiple plume jets imparted disturbance torque on the spacecraft resulting in small but visible attitude control errors. Using the known and unique transfer function between the disturbance torque and the attitude control error, the collected attitude control error telemetry could be used to estimate the disturbance torque. The effectiveness of this methodology is confirmed using the E3 telemetry data. Given good estimates of spacecraft's projected area, center of pressure location, and spacecraft velocity, the time history of the Enceladus plume density is reconstructed accordingly. The 1-sigma uncertainty of the estimated density is 7.7%. Next, we modeled the density due to each plume jet as a function of both the radial and angular distances of the spacecraft from the plume source. We also conjecture that the total plume density experienced by the spacecraft is the sum of the component plume densities. By comparing the time history of the reconstructed E3 plume density with that predicted by the plume model, values of the plume model parameters are determined. Results obtained are compared with those determined by other Cassini science instruments.

Lee, Allan Y.↗

System-Independent X-Ray Characterization

We have developed a method for image reconstruction and analysis of dual-energy X-ray scans that has been demonstrated to yield a physically intuitive characterization of materials (by electron density and effective atomic number) in a system-independent way. Our methodology stands in contrast to common practice of reconstructing images to yield linear attenuation coefficients, which depend strongly on system spectral response. At the core of our method, a technique called dual-energy decomposition is used to generate spectral energy-independent characterizations of the attenuation properties of a material. Inputs to this decomposition include system spectral response estimates and dual-energy X-ray projection data. X-ray cross-section tables are used to convert decomposition results to physical feature estimates (electron density, effective atomic number). Our method currently works best for materials with atomic numbers below 22 (titanium). In prior work for DHS, we demonstrated in a controlled, no-clutter environment that the decomposition technique could be applied to scans of a well-characterized materials on multiple systems with different X-ray spectral responses. In cases where LAC values showed a standard deviation of as much as 20% of mean value, application of our method yielded results with standard deviations as low as 1–3% of mean value.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Variability Timescale and Spectral Index of Sgr A* in the Near Infrared: Approximate Bayesian Computation Analysis of the Variability of the Closest Supermassive Black Hole

Sagittarius A* (Sgr A* ) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We present an analysis of the most comprehensive NIR variability data set of Sgr A* to date: eight 24 hr epochs of continuous monitoring of Sgr A* at 4.5 μm with the IRAC instrument on the Spitzer Space Telescope, 93 epochs of 2.18 μm data from Naos Conica at the Very Large Telescope, and 30 epochs of 2.12 μm data from the NIRC2 camera at the Keck Observatory, in total 94,929 measurements. A new approximate Bayesian computation method for fitting the first-order structure function extracts information beyond current fast Fourier transformation (FFT) methods of power spectral density (PSD) estimation. With a combined fit of the data of all three observatories, the characteristic coherence timescale of Sgr A* is b 243 57 82 t = - + minutes (90% credible interval). The PSD has no detectable features on timescales down to 8.5 minutes (95% credible level), which is the ISCO orbital frequency for a dimensionless spin parameter a = 0.92. One light curve measured simultaneously at 2.12 and 4.5 μm during a low flux-density phase gave a spectral index αs = 1.6 ± 0.1 (F µ n s n -a ). This value implies that the Sgr A* NIR color becomes bluer during higher flux-density phases. The probability densities of flux densities of the combined data sets are best fit by log-normal distributions. Based on these distributions, the Sgr A* spectral energy distribution is consistent with synchrotron radiation from a nonthermal electron population from below 20 GHz through the NIR.

accretion↗

Converting sWeights to probabilities with density ratios

The use of machine learning approaches continues to have many benefits in experimental nuclear and particle physics. One common issue is generating training data which is sufficiently realistic to give reliable results. Here we advocate using real experimental data as the source of training data and demonstrate how one might subtract background contributions through the use of probabilistic weights which can be readily applied to training data. The sPlot formalism is a common tool used to isolate distributions from different sources. However, the negative sWeights produced by the sPlot technique can cause training problems and poor predictive power. This article demonstrates how density ratio estimation can be applied to convert sWeights to event probabilities, which we call drWeights. The drWeights can then be applied to produce the distributions of interest and are consistent with direct use of the sWeights. This article will also show how decision trees are particularly well suited to convert sWeights, with the benefit of fast prediction rates and adaptability to aspects of experimental data such as the data sample size and proportions of different event sources. We also show that a density ratio product approach in which the initial drWeights are reweighted by an additional converter gives substantially better results.

Glazier, Derek I. [Univ. of Glasgow, Scotland (Uni↗

Tentative Identification of Interstellar CS by UV Absorption in the C-X Transition

Current chemical models of translucent clouds imply that sulfur is depleted with respect to its average abundance in diffuse clouds as determined by Morton (1975 ApJ,197, 85) by factors exceeding 10. Existing gigahertz emission observations yield sizable uncertainties in the column densities to background stars, and attempts to measure column densities of CS in the lines-of-sight to various stars through the absorption in the A-X (0,0) electronic transition near 257.6 nm have been unsuccessful. CS dissociates primarily through discrete absorptions into predissociating states. In analogy with CO the strongest of these is the C-X band found near 140.0 nm. We report the tentative identification of the C-X band of CS in the interstellar spectra of HD 24398 (zeta Per), HD 148184 (chi Oph), and HD 154368. The estimated column densities will be compared to the predictions of prevailing chemical models.

Smith, Andrew↗

Infrared and X-Ray Spectroscopy of the Kes 75 Supernova Shell Characterizing the Dust and Gas Properties

We present deep Chandra observations and Spitzer Space Telescope infrared (IR) spectroscopy of the shell in the composite supernova remnant (SNR) Kes 75 (G29.7-0.3). The remnant is composed of a central pulsar wind nebula and a bright partial shell in the south that is visible at radio, IR, and X-ray wavelengths. The X-ray emission can be modeled by either a single thermal component with a temperature of approx 1.5 keV, or with two thermal components with temperatures of 1.5 and 0.2 keY. Previous studies suggest that the hot component may originate from reverse-shocked SN ejecta. However, our new analysis shows no definitive evidence for enhanced abundances of Si, S, Ar, Mg, and Fe, as expected from supernova (SN) ejecta, or for the IR spectral signatures characteristic of confirmed SN condensed dust, thus favoring a circumstellar or interstellar origin for the X-ray and IR emission. The X-ray and IR emission in the shell are spatially correlated, suggesting that the dust particles are collisionally heated by the X-ray emitting gas. The IR spectrum of the shell is dominated by continuum emission from dust with little, or no line emission. Modeling the IR spectrum shows that the dust is heated to a temperature of approx 140 K by a relatively dense, hot plasma, that also gives rise to the hot X-ray emission component. The density inferred from the IR emission is significantly higher than the density inferred from the X-ray models, suggesting a low filling factor for this X-ray emitting gas. The total mass of the warm dust component is at least 1.3 x 10(exp -2) Solar Mass, assuming no significant dust destruction has occurred in the shell. The IR data also reveal the presence of an additional plasma component with a cooler temperature, consistent with the 0.2 keV gas component. Our IR analysis therefore provides an independent verification of the cooler component of the X-ray emission. The complementary analyses of the X-ray and IR emission provide quantitative estimates of density and filling factors of the clumpy medium swept up by the SNR.

Temim, Tea↗

Infrared and X-Ray Spectroscopy of the KES 75 Supernova Remnant Shell: Characterizing the Dust and Gas Properties

We present deep Chandra observations and Spitzer Space Telescope infrared (IR) spectroscopy of the shell in the composite supernova remnant (SNR) Kes 75 (G29.7-0.3). The remnant is composed of a central pulsar wind nebula and a bright partial shell in the south that is visible at radio, IR, and X-ray wavelengths. The X-ray emission can be modeled by either a single thermal component with a temperature of ~ 1.5 keV, or with two thermal components with temperatures of 1.5 and 0.2 keV. Previous studies suggest that the hot component may originate from reverse-shocked SN ejecta. However, our new analysis shows no definitive evidence for enhanced abundances of Si, S, Ar, Mg, and Fe, as expected from supernova (SN) ejecta, or for the IR spectral signatures characteristic of confirmed SN condensed dust, thus favoring a circumstellar or interstellar origin for the X-ray and IR emission. The X-ray and ill emission in the shell are spatially correlated, suggesting that the dust particles are collisionally heated by the X-ray emitting gas. The IR spectrum of the shell is dominated by continuum emission from dust with little, or no line emission. Modeling the IR spectrum shows that the dust is heated to a temperature of ~ 140 K by a relatively dense, hot plasma, that also gives rise to the hot X-ray emission component. The density inferred from the IR emission is significantly higher than the density inferred from the X-ray models, suggesting a low filling factor for this X-ray emitting gas. The total mass of the warm dust component is at least 1.3 x 10(exp -2) solar mass, assuming no significant dust destruction has occurred in the shell. The IR data also reveal the presence of an additional plasma component with a cooler temperature, consistent with the 0.2 keV gas component. Our IR analysis therefore provides an independent verification of the cooler component of the X-ray emission. The complementary analyses of the X-ray and IR emission provide quantitative estimates of density and filling factors of the clumpy medium swept up by the SNR.

Temim, Tea↗

Infrared and X-Ray Spectroscopy of the Kes 75 Supernova Remnant Shell: Characterizing the Dust and Gas Properties

We present deep Chandra observations and Spitzer Space Telescope infrared (IR) spectroscopy of the shell in the composite supernova remnant (SNR) Kes 75 (G29.7-0.3). The remnant is composed of a central pulsar wind nebula and a bright partial shell in the south that is visible at radio, IR, and X-ray wavelengths. The X-ray emission can be modeled by either a single thermal component with a temperature of approximately 1.5 keY, or with two thermal components with temperatures of 1.5 and 0.2 keY. Previous studies suggest that the hot component may originate from reverse-shocked supernova (SN) ejecta. However, our new analysis shows no definitive evidence for enhanced abundances of Si, S, Ar, Mg, and Fe, as expected from SN ejecta, or for the IR spectral signatures characteristic of confirmed SN condensed dust, thus favoring a circumstellar or interstellar origin for the X-ray and IR emission. The X-ray and IR emission in the shell are spatially correlated, suggesting that the dust particles are collisionally heated by the X-ray emitting gas. The IR spectrum of the shell is dominated by continuum emission from dust with little, or no line emission. Modeling the IR spectrum shows that the dust is heated to a temperature of approximately 140 K by a relatively dense, hot plasma that also gives rise to the hot X-my emission component. The density inferred from the IR emission is significantly higher than the density inferred from the X-ray models, suggesting a low filling factor for this X-my emitting gas. The total mass of the warm dust component is at least 1.3 x 10(exp -2) x solar mass, assuming no significant dust destruction has occurred in the shell. The IR data also reveal the presence of an additional plasma component with a cooler temperature, consistent with the 0.2 keV gas component. Our IR analysis therefore provides an independent verification of the cooler component of the X-ray emission. The complementary analyses of the X-ray and IR emission provide quantitative estimates of density and filling factors of the clumpy medium swept up by the SNR.

Temim, Tea↗

On the flux and the energy spectrum of interstellar ions in the solar system

The flux density of ions created by ionization of interstellar neutral particles in the solar system and picked up by the solar wind is calculated as a function of the neutral particles. For atomic hydrogen the flux density is estimated to exceed 10,000/sq cm/sec over the distance range from a few to nearly 100 AU. The velocity space distribution of the interstellar ions is calculated under the assumption of no significant energy diffusion but with inclusion of adiabatic effects as well as a possible strong pitch angle diffusion. The energy spectrum is highly nonthermal and much broader than that of the solar wind ions; interstellar protons are easily distinguishable from solar wind protons by their location in velocity space. If charge exchange is an important contributor to the ionization of hydrogen, the observed local intensity of interstellar protons should exhibit time variations correlated with the density changes of the solar wind stream structure.

Vasyliunas, V. M.↗

Multiple Asteroid Systems: Dimensions and Thermal Properties from Spitzer Space Telescope and Ground-based Observations

We collected mid-IR spectra from 5.2 to 38 microns using the Spitzer Space Telescope Infrared Spectrograph of 28 asteroids representative of all established types of binary groups. Photometric light curves were also obtained for 14 of them during the Spitzer observations to provide the context of the observations and reliable estimates of their absolute magnitudes. The extracted mid-IR spectra were analyzed using a modified standard thermal model (STM) and a thermophysical model (TPM) that takes into account the shape and geometry of the large primary at the time of the Spitzer observation. We derived a reliable estimate of the size, albedo, and beaming factor for each of these asteroids, representing three main taxonomic groups: C, S, and X. For large (volume-equivalent system diameter Deq > 130 km) binary asteroids, the TPM analysis indicates a low thermal inertia (Lambda < or = approx.100 J/1/2 s/K/sq m2) and their emissivity spectra display strong mineral features, implying that they are covered with a thick layer of thermally insulating regolith. The smaller (surface-equivalent system diameter Deff < 17 km) asteroids also show some emission lines of minerals, but they are significantly weaker, consistent with regoliths with coarser grains, than those of the large binary asteroids. The average bulk densities of these multiple asteroids vary from 0.7-1.7 g/cu cm (P-, C-type) to approx. 2 g/cu cm (S-type). The highest density is estimated for the M-type (22) Kalliope (3.2 +/- 0.9 g/cu cm). The spectral energy distributions (SEDs) and emissivity spectra, made available as a supplement document, could help to constrain the surface compositions of these asteroids.

Marchis, F.↗

Improving Ocean Reanalyses and ENSO Forecasts By Assimilation of Rain-Corrected Satellite Sea Surface Salinity Using the GMOA S2s Forecast System

The ENSO phenomenon has a significant global socio-economic impact and has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of satellite altimetry and subsurface temperature and salinity from (mostly) Argo help improve the initialization of the thermocline, while satellite SST aids in constraining surface heat-fluxes, leading to improved coupled system sub-seasonal to seasonal forecasts. However, few studies have focused on improving the near-surface density and mixing through satellite sea surface salinity (SSS) assimilation. The few ocean models that assimilate satellite SSS, bias correct to normalize towards the near-surface Argo data for expediency. This assumption is likely inadequate in rainy regions, where buoyant water forms a fresh surface lens. In previous work, we showed that adjusting SSS to bulk salinity (Sb) using the Rain Impact Model (RIM) of Santos-Garcia et al., 2014 improves the near-surface density and mixed layer depth, leading to deeper thermocline and improved NINO3.4 SST forecasts. We now utilize the Soil Moisture and Ocean Salinity rain-corrected (SMOS_RC) SSS, available in SMOS-CATDS products, to represent Sb more accurately at the first model layer (e.g., 5 m). Rather than a diffusivity model as RIM, SMOS_RC uses a statistical correction dependent on Integrated Multi-satellitE Retrievals for GPM (IMERG) rain rates, established on observed SMOS SSS decreases related to Sb in the presence of rain (Supply et al., 2020). For all experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF scheme (Penny et al., 2013). One reanalysis additionally assimilates SMOS SSS data as is, and a separate reanalysis assimilates SMOS_RC. We assess the impact on near-surface and subsurface dynamics by validating against observations and explore how SSS assimilation (SMOS vs SMOS_RC) impacts ENSO forecasts using the NASA GMAO Sub-seasonal to Seasonal coupled forecast system (S2S-v3, Molod et al. 2020). We show that improved estimates of density and near-surface mixing led to more accurate coupled air/sea interaction and better ENSO forecasts. The increased SSS, resulting from the removal of the instantaneous rain effect, modifies the ocean state by enhancing mixing and deepening the thermocline.

Veronica Ruiz Xomchuk↗

Minimum entropy filtering for a single output non-Gaussian stochastic system using state transformation

This paper presents a novel filter design for the single-output stochastic non-linear systems subjected to non-Gaussian noises and the proposed assumptions. Based on a state transformation, the unmeasurable states of the systems can be estimated where non-linear terms in the systems have been eliminated. It has been shown that the estimation error is linearly dynamical regarding to the presented vector-valued filter gain which can be optimised by minimising the entropy-based performance criterion. In addition, the convergence of the presented algorithm is analysed in mean-square sense and a numerical example is given to verify the effectiveness of the presented filtering algorithm. Meanwhile, the extended Kalman filter, unscented particle filter and minimum entropy filter are given for the comparisons of the filtering performance. Following the presented framework, some extensions of the presented filtering algorithm are discussed to indicate the flexibility of the filter design. The contribution of this paper can be summarised as establishing a novel minimum entropy filtering framework which consists of model transformation, entropy optimisation and convergence analysis.

42 ENGINEERING↗

Estimation of distributions via multilevel Monte Carlo with stratified sampling

We design and implement a novel algorithm for computing a multilevel Monte Carlo (MLMC) estimator of the joint cumulative distribution function (CDF) of a vector-valued quantity of interest in problems with random input parameters and initial conditions. Our approach combines MLMC with stratified sampling of the input sample space by replacing standard Monte Carlo at each level with stratified Monte Carlo initialized with proportionally allocated samples. We show that the resulting stratified MLMC (sMLMC) algorithm is more efficient than its standard MLMC counterpart due to the additional variance reduction provided by the stratification of the random parameter's domain, especially at the coarsest levels. Additional computational cost savings are obtained by smoothing the indicator function with a Gaussian kernel, which proves to be an efficient and robust alternative to recently developed polynomial-based techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Probabilistic learning and updating of a digital twin for composite material systems

This paper presents an approach for characterizing and estimating statistical dependence between a large number of observables in a composite material system. Conditional regression is carried out using the estimated joint density function, permitting a systematic exploration of interdependence between fine scale and coarse observables that can be used for both prognosis and design of complex material systems. An example demonstrates the integration of experimental data with a computational database. The statistical approach is based on the probabilistic learning on manifolds recently developed by the authors. This approach leverages intrinsic structure detected through diffusion on graphs with projected stochastic differential equations to generate samples constrained to that structure.

42 ENGINEERING↗

Uncovering acoustic signatures of pore formation in laser powder bed fusion

Abstract We present a machine learning workflow to discover signatures in acoustic measurements that can be utilized to create a low-dimensional model to accurately predict the location of keyhole pores formed during additive manufacturing processes. Acoustic measurements were sampled at 100 kHz during single-layer laser powder bed fusion (LPBF) experiments, and spatio-temporal registration of pore locations was obtained from post-build radiography. Power spectral density (PSD) estimates of the acoustic data were then decomposed using non-negative matrix factorization with custom $$\varvec{k}$$ k -means clustering (NMF $$\varvec{k}$$ k ) to learn the underlying spectral patterns associated with pore formation. NMF $$\varvec{k}$$ k returned a library of basis signals and matching coefficients to blindly construct a feature space based on the PSD estimates in an optimized fashion. Moreover, the NMF $$\varvec{k}$$ k decomposition led to the development of computationally inexpensive machine learning models which are capable of quickly and accurately identifying pore formation with classification accuracy of supervised and unsupervised label learning greater than 95% and 90%, respectively. The intrinsic data compression of NMF k , the relatively light computational cost of the machine learning workflow, and the high classification accuracy makes the proposed workflow an attractive candidate for edge computing toward in-situ keyhole pore prediction in LPBF.

36 MATERIALS SCIENCE↗

A Multi-Region SEIR Model Incorporating Inter-County Mobility and Time-Dependent Transmission Dynamics: Application to COVID-19 Disease Outbreak Data in North Carolina.

Classical infectious disease compartmental models typically do not incorporate spatial heterogeneity or mobility. We develop a multi-region susceptible-exposed-infected-recovered (SEIR) model in which disease dynamics are coupled via inter-region mobility and the transmission rate is both region and time dependent. We calibrate the model using rolling averages of daily COVID-19 data in all 100 North Carolina counties. Mobility parameters are prescribed using daily inter-county commuter data. The number of transmission rate parameters is substantially reduced by hypothesizing that the dynamics correlate with county-level population density. Parameter estimation is carried out using several objective functions with error terms at different scales. An additive combination of least squares error at the county-level and the state-level, along with a quadratic transmission rate polynomial, yields the lowest overall error at both spatial scales. The calibrated model is used to simulate regional effects of perturbing disease transmission rates in adjacent counties and to illustrate effects of the state’s mobility infrastructure on disease dynamics and spread for a new disease outbreak.

COVID-19 modeling↗

Thermal diffusivity recovery and defect annealing kinetics of self-ion implanted tungsten probed by insitu transient grating spectroscopy

Tungsten is a promising candidate material for plasma-facing armour components in future fusion reactors. A key concern is irradiation-induced degradation of its normally excellent thermal transport properties. In this comprehensive study, thermal diffusivity degradation in ion-implanted tungsten and its evolution from room temperature (RT) to 1073 K is considered. Five samples were exposed to 20 MeV self-ions at RT to achieve damage levels ranging from 3.2 x 10 -4 to 3.2 displacements per atom (dpa). Transient grating spectroscopy with in-situ heating was then used to study thermal diffusivity evolution as a function of temperature. Using a kinetic theory model, an equivalent point defect density is estimated from the measured thermal diffusivity. The results showed a prominent recovery of thermal diffusivity between 450 K and 650 K, which coincides with the onset of mono-vacancy mobility. After 1073 K annealing samples with initial damage of 3.2 x 10 -3 dpa or less recover close to the pristine value of thermal diffusivity. For doses of 3.2 x 10 -2 dpa or higher, on the other hand, a residual reduction in thermal diffusivity remains even after 1073 K annealing. Transmission electron microscopy reveals that this is associated with extended, irradiation-induced dislocation structures that are retained after annealing. A sensitivity analysis shows that thermal diffusivity provides an efficient tool for the assessing total defect content in tungsten up to 1000 K.

36 MATERIALS SCIENCE↗