Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed inference”

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 217 records · Page 12

Role of Surface Hydroxyls and Lattice Oxygen in Governing Selectivity and Restructuring During Acetic Acid Conversion on Fe 3 O 4 (001)

Understanding the reactivity of carboxylic acids on metal oxide surfaces is critical for elucidating ketonization mechanisms relevant to biomass upgrading. Here, we investigate the adsorption and thermal decomposition of acetic acid (CH 3 COOH) on Fe 3 O 4 (001) using scanning tunneling microscopy (STM), temperature-programmed reaction spectroscopy (TPRS), and X-ray photoelectron spectroscopy (XPS). At room temperature, acetic acid adsorbs dissociatively to form ordered bidentate acetate (CH 3 COO) overlayer that lifts the (√2 × √2)R 45° surface reconstruction. TPRS reveals ketene (CH 2 CO) as the dominant product, along with CO, CO 2 , and H 2 O, the latter evolving via a Mars–van Krevelen (MvK) mechanism. Isotopic labeling shows preferential CO 2 formation from the carboxyl carbon and a more balanced CO/CO 2 ratio from the methyl carbon, suggesting distinct oxidation pathways. STM imaging reveals embedded acetate intermediates filling surface oxygen vacancies created in MvK steps. Upon product formation completion (~700 K), extensive surface etching is observed, with pits elongated along the octahedral Fe rows. Approximately 20% of the surface oxygen is removed, consistent with vacancy formation stoichiometry inferred from product distributions. These findings demonstrate that carboxylate-induced restructuring of Fe 3 O 4 (001) involves both surface healing and reduction processes, offering mechanistic insights relevant to ketonization and broader carboxylic acid chemistry on metal oxides.

acetic acid↗

A Framework for Parametric and Predictive Uncertainty Quantification in the E3SM Land Model: Assessing Site and Observable Generalizability

Quantifying parametric uncertainty using observations from individual sites provides a critical foundation for Earth system modeling, serving as a necessary first step before scaling up to regional or global applications. This study introduces a novel computational framework designed to enhance model predictability by reducing parametric uncertainty and assessing site and observable generalizability using various observational constraints. The framework integrates five components: Model Simulation, Statistical Emulation, Global Sensitivity Analysis (GSA), Model Calibration, and Model Prediction. Using the E3SM land model, we simulated site-level land-atmosphere carbon and energy fluxes from 2003 to 2007 across five evergreen needleleaf FLUXNET sites, perturbing 26 vegetation-related model parameters. Gaussian process emulators were employed to expedite GSA and model calibration. Four critical parameters that strongly influence selected land-atmosphere fluxes were identified by GSA. Bayesian approaches were used to infer parameter probability distributions leveraging synthetic data and FLUXNET observations. The results reveal that posterior parameter distributions vary significantly across different sites and observables within the same plant functional type. Probabilistic predictions indicate that parameters calibrated at one site can enhance predictive accuracy at other sites, although site heterogeneity may sometimes outweigh parametric uncertainty. Additionally, the probabilistic predictions demonstrate that calibration for one variable can also improve predictability for other variables, thereby maximizing predictive capabilities with limited observations. This framework provides a powerful approach for reducing parametric uncertainty in Earth system models and deepening our understanding of carbon dynamics and energy cycles. Its adaptability makes it a valuable tool for broader applications in Earth system modeling.

54 ENVIRONMENTAL SCIENCES↗

A forward modeling approach to analyzing galaxy clustering with S IM BIG

We present cosmological constraints from a simulation-based inference (SBI) analysis of galaxy clustering from the SimBIG forward modeling framework. SimBIG leverages the predictive power of high-fidelity simulations and provides an inference framework that can extract cosmological information on small nonlinear scales. In this work, we apply SimBIG to the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy sample and analyze the power spectrum, P ℓ (k), to k max = 0.5 h/Mpc. We construct 20,000 simulated galaxy samples using our forward model, which is based on 2,000 high-resolution Quijote N -body simulations and includes detailed survey realism for a more complete treatment of observational systematics. We then conduct SBI by training normalizing flows using the simulated samples and infer the posterior distribution of ΛCDM cosmological parameters: Ω m , Ω b , h, n s , σ 8 . We derive significant constraints on Ω m and σ 8 , which are consistent with previous works. Our constraint on σ 8 is 27% more precise than standard P ℓ analyses because we exploit additional cosmological information on nonlinear scales beyond the limit of current analytic models, k > 0.25 h/Mpc. This improvement is equivalent to the statistical gain expected from a standard P ℓ analysis of galaxy sample ~ 60% larger than CMASS. While we focus on P ℓ in this work for validation and comparison to the literature, SimBIG provides a framework for analyzing galaxy clustering using any summary statistic. We expect further improvements on cosmological constraints from subsequent SimBIG analyses of summary statistics beyond P ℓ .

79 ASTRONOMY AND ASTROPHYSICS↗

Characterization of core neutrals using passive measurements of the D α spectrum near the X -point in the DIII-D tokamak

New spectroscopic measurements of deuterium Balmer-α emission are used to infer the spatial distribution of neutral particles near the X-point in a diverted high confinement mode plasma. The charge exchange neutral spectroscopy (CENS) diagnostic—recently installed on the DIII-D tokamak—uses 15 lines of sight extending from the edge of the confined region (ρ ≈ 0.7) to the X-point in lower single-null plasmas. Each CENS line of sight is spectrally resolved around the deuterium Balmer-α wavelength (6561 Å) to measure the Doppler shift, line broadening, and intensity of Dα emission. Thermal line broadening is used to identify emission from neutrals undergoing charge-exchange (CX) with high temperature ions in the confined plasma. This spectral information allows a more accurate determination of the neutral density deeper inside the confined plasma compared to traditional filter-based diagnostics. We present two methods of analyzing CENS measurements: (1) fitting the spectrum from each view-chord separately, and (2) a tomographic method for inverting the neutral density over a 2D region of space using all CENS views collectively. The neutral density profile is found to decay exponentially in the radial direction with two scale-lengths. In the pedestal region the neutral density decays at a rate approximately equal to the local mean-free-path for CX collisions, $L_{n_\mathrm{D0}}\approx \lambda_\text{CX}$. Further inside the plasma the neutral density decays at a rate equal to the mean-free-path for ionization, $L_{n_\mathrm{D0}}\approx \lambda_\text{inz.}$. The separatrix $n_\mathrm{D0}$ value is found to be approximately $2\times10^{15}$ m−3. Based on results from the 2D inversion, the density of neutrals is found to increase along the separatrix approaching the X-point by a factor of $\approx 5\times$ over the region covered by the CENS diagnostic.

X-point neutrals↗

Halo-independent analysis of direct dark matter detection through electron scattering

Sub-GeV mass dark matter particles whose collisions with nuclei would not deposit sufficient energy to be detected, could instead be revealed through their interaction with electrons. Analyses of data from direct detection experiments usually require assuming a local dark matter halo velocity distribution. In the halo-independent analysis method, properties of this distribution are instead inferred from direct dark matter detection data, which allows then to compare different data without making any assumption on the uncertain local dark halo characteristics. Furthermore, this method has so far been developed for and applied to dark matter scattering off nuclei. Here we demonstrate how this analysis can be applied to scattering off electrons.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimal Bayesian supervised domain adaptation for RNA sequencing data

Abstract Motivation When learning to subtype complex disease based on next-generation sequencing data, the amount of available data is often limited. Recent works have tried to leverage data from other domains to design better predictors in the target domain of interest with varying degrees of success. But they are either limited to the cases requiring the outcome label correspondence across domains or cannot leverage the label information at all. Moreover, the existing methods cannot usually benefit from other information available a priori such as gene interaction networks. Results In this article, we develop a generative optimal Bayesian supervised domain adaptation (OBSDA) model that can integrate RNA sequencing (RNA-Seq) data from different domains along with their labels for improving prediction accuracy in the target domain. Our model can be applied in cases where different domains share the same labels or have different ones. OBSDA is based on a hierarchical Bayesian negative binomial model with parameter factorization, for which the optimal predictor can be derived by marginalization of likelihood over the posterior of the parameters. We first provide an efficient Gibbs sampler for parameter inference in OBSDA. Then, we leverage the gene-gene network prior information and construct an informed and flexible variational family to infer the posterior distributions of model parameters. Comprehensive experiments on real-world RNA-Seq data demonstrate the superior performance of OBSDA, in terms of accuracy in identifying cancer subtypes by utilizing data from different domains. Moreover, we show that by taking advantage of the prior network information we can further improve the performance. Availability and implementation The source code for implementations of OBSDA and SI-OBSDA are available at the following link. https://github.com/SHBLK/BSDA. Supplementary information Supplementary data are available at Bioinformatics online.

Biochemistry & Molecular Biology↗

DAmodel: hierarchical Bayesian modelling of DA white dwarfs for spectrophotometric calibration

We use hierarchical Bayesian modelling to calibrate a network of 32 all-sky faint DA white dwarf (DA WD) spectrophotometric standards (⁠16.5 < V , 19.5⁠) alongside three CALSPEC standards, from 912 Å to 32 μm. The framework is the first of its kind to jointly infer photometric zero points and WD parameters (surface gravity log g⁠, effective temperature T eff ⁠, extinction A V ⁠, dust relation parameter R V ) by simultaneously modelling both photometric and spectroscopic data. We model panchromatic Hubble Space Telescope Wide Field Camera 3 (HST/WFC3) UVIS and IR photometry, HST/STIS UV spectroscopy, and ground-based optical spectroscopy to sub-per cent precision. Photometric residuals for the sample are the lowest yet yielding < 0.004 mag RMS on average from the UV to the NIR, achieved by jointly inferring time-dependent changes in system sensitivity and WFC3/IR count-rate nonlinearity. Our GPU-accelerated implementation enables efficient sampling via Hamiltonian Monte Carlo, critical for exploring the high-dimensional posterior space. The hierarchical nature of the model enables population analysis of intrinsic WD and dust parameters. Inferred spectral energy distributions from this model will be essential for calibrating the James Webb Space Telescope as well as next-generation surveys, including Vera Rubin Observatory’s Legacy Survey of Space and Time and the Nancy Grace Roman Space Telescope.

methods: statistical↗

FiXR : a framework to reconstruct fiber cross-sections from X-ray fiber diffraction experiments

Ab initio reconstruction methods have revolutionized the capabilities of small-angle X-ray scattering (SAXS), allowing the data-driven discovery of previously unknown molecular conformations, exploiting optimization heuristics and assumptions behind the composition of globular molecules. While these methods have been successful for the analysis of small particles, their impact on fibrillar assemblies has been more limited. The micrometre-range size of these assemblies and the complex interaction of their periodicities in their scattering profiles indicate that the discovery of fibril structures from SAXS measurements requires novel approaches beyond extending existing tools for molecular discovery. In this work, it is proposed to use SAXS measurements, together with diffraction theory, to infer the electron distribution of the average cross-section of a fiber. This cross-section is modeled as a discrete electron density with continuous support, allowing representations beyond binary distributions. Additional constraints, such as non-negativity or smoothness/connectedness, can also be added to the framework. The proposed approach is tested using simulated SAXS data from amyloid β fibril models and using measured data of Tobacco mosaic virus from SAXS experiments, recovering the geometry and density of the cross-sections in all cases. The approach is further tested by analyzing SAXS data from different amyloid β fibril assemblies, with results that are in agreement with previously proposed models from cryo-EM measurements. The limitations of the proposed method, together with an analysis of the robustness of the method and the combination with different experimental sources, are also discussed.

Biochemistry & Molecular Biology↗

Data-Driven Probabilistic Voltage Risk Assessment of MiniWECC System with Uncertain PVs and Wind Generations using Realistic Data

Here, it is found from actual data that due to generation dispatch and uncertain renewable generations and loads with complicated correlations, inferring the probabilistic distributions for uncertain inputs is challenging. Many probabilistic power flow approaches have been developed in the literature but their validations using realistic systems and data are lacking. This paper proposes a data-driven probabilistic analysis approach for system risk assessment of the miniWECC system using actual data. The sparse Gaussian process (SGP) is advocated to quantify the impacts of uncertain inputs on voltage security. SGP does not need the probability distribution function of uncertain inputs, can handle correlations and is highly computationally efficient. Results on the miniWECC system using realistic data show that SGP outperforms existing approaches and is able to quantify the voltage violation risks.

17 WIND ENERGY↗

Bayesian High-Rank Hankel Matrix Completion for Nonlinear Synchrophasor Data Recovery

Phasor measurement units (PMUs) provide high temporal-resolution synchrophasor measurements for power system monitoring and control. The frequent data quality issues, such as missing and bad data, prevent the incorporation of synchrophasor data in real-time operations. Most existing data-driven data recovery methods assume the power system dynamics can be approximated by a linear dynamical system, and the recovery performance degrades significantly when the power system is experiencing nonlinear dynamics during significant events. Here, this paper proposes a data-driven Bayesian nonlinear synchrophasor data recovery method (Ba-NSDR) that can recover a consecutive time period of simultaneous data losses or errors across all channels, even when the underlying system is highly nonlinear. The idea is to lift the Hankel matrix of the spatial-temporal synchrophasor data to a higher dimension such that the lifted Hankel matrix is low-rank in that space and can be processed with the kernel trick. Our proposed Bayesian method then infers the probabilistic distributions of synchrophasor from the partial observations. Some distinctive features of Ba-NSDR include an uncertainty index to measure the accuracy of the recovery result and the robustness to parameter selections. Our method is verified on both synthetic and recorded event datasets.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Shell structure in 32P

The nuclear structure of 32P was investigated by detecting multi-γ coincidences along with charged light ions following the bombardment of an 18O target with a 30-MeV 16O beam from the Florida State University accelerator facility using the Clarion2-Trinity array of Clover γ spectrometers and GAGG (Gd-Al-Ga-garnet) scintillators for charged particles. A number of new states up to 12 MeV excitation and spins up to 9+ were observed by their γ decay patterns. Spins and parities were assigned by comparing measured γ angular distributions with polarizations inferred from Compton-scattering asymmetries between the Ge crystals in the Clover spectrometers. The new level scheme of 32P compares well with shell model calculations using the FSU cross-shell interaction and reasonably well with those using the sdpf-m interaction. In particular, structures with higher spin arise from promotion of more nucleons up to the f-p shell.

Allmond, James [ORNL] (ORCID:0000000165338721)↗

Neutrons produced by known energies of ions abundant in space

Particle accelerator radiation measurements are applied to the problem of calculating biological dose from radiation produced in the walls of a spacecraft by various ions in space. Neutrons, one of the products of the interactions of energetic ions with matter, are usually quite penetrating and have large values of Q.F. or R.B.E. Ions of helium, boron, carbon, nitrogen, and oxygen were accelerated and directed onto target materials of copper or tantalum. The secondary neutron production was determined. Studies were made of the angular distribution and an inferred neutron spectrum was calculated from activities of threshold reaction detectors.

Wadman, W. W., III↗

Cloud structure and heating rates in the atmosphere of Venus

Ground-based observations and Venera 8 entry-probe measurements are used to infer the vertical distribution of cloud particles in the atmosphere of Venus. In the cloud-top region, from a few millibars to a few hundred millibars pressure, the mixing ratio of cloud particles to gas increases with depth. The visible clouds are diffuse with a scale height about one-half of the gaseous atmosphere. Although the presence of significant vertical structure could escape detection by available observations, the diffuse haze appears to extend over at least 20 km in altitude. The Venera 8 measurements suggest considerable vertical structure in the deep atmosphere. A unique solution for the cloud structure is not possible, but if it is assumed that the cloud optical properties are independent of height, then some characteristics of the relative cloud structure can be deduced. Under this assumption, the results show a maximum cloud density near 40 km, a nearly homogeneous particle mixing in the region from about 40 to 50 km, and a fairly sharp cloud bottom near 30 km. Relative maxima in the cloud density are also implied near 55 and 10 km, but with much greater uncertainty.

Lacis, A. A.↗

Diterpenoid compounds and other lipids in deep-sea sediments and their geochemical significance

Cyclic diterpenoid compounds have been found by various investigators in the geosphere (e.g., fossil resins, coals, soil, shale, and deep-sea sediments). These compounds occur in significant amounts only in higher plants and are therefore potential markers of terrigenous plant lipids. Diterpenoids with the abietane skeleton (mainly dehydroabietic acid) have been identified in the lipids of sediment samples from the northeast Pacific Ocean, Black Sea, and North Atlantic Ocean. The presence of these resin-derived compounds was correlated with the terrigenous clay components and with the presence of pollen. The presence of polycyclic diterpenoids was also correlated with the distribution patterns and inferred sources of other sediment lipid constituents. Potamic transport, followed by turbidite redistribution, is the probable input mechanism of these resin-derived compounds to the deep-sea sediments. These diterpenoids appear to be excellent biological markers of resinous higher plants.

Simoneit, B. R. T.↗

An extension to the Chahine method of inverting the radiative transfer equation

An extension of the Chahine relaxation method (1970) for inverting the radiative transfer equation is presented. This method is superior to the original method in that it takes into account in a realistic manner the shape of the kernel function, and its extension to nonlinear systems is much more straightforward. A comparison of the new method with a matrix method due to Twomey (1965), in a problem involving inference of vertical distribution of ozone from spectroscopic measurements in the near ultraviolet, indicates that in this situation this method is stable with errors in the input data up to 4%, whereas the matrix method breaks down at these levels. The problem of non-uniqueness of the solution, which is a property of the system of equations rather than of any particular algorithm for solving them, remains, although it takes on slightly different forms for the two algorithms.

Twomey, S.↗

Topside sounder observations of equatorial bubbles

Large scale regions of depleted equatorial ionospheric plasma, called equatorial bubbles, are investigated using topside sounder data. The sounder's unique remote measuring capability enables the magnetic field-aligned nature of the bubbles to be investigated. A search of all available Alouette 2 and ISIS 1 ionograms during nighttime perigee passes near the magnetic equator has revealed a variety of echo signatures associated with bubbles. In addition to a sudden drop in electron density, these signatures usually include in situ spread F and ducted traces. The ducted traces have been used to determine the electron density distribution and to infer changes in ion composition along the magnetic field line within the duct associated with the bubble. In some cases it can be determined that the bubble is asymmetric with respect to the magnetic equator. Even though such features require 3 dimensional models for their explanation, the great field-aligned extent of the bubbles (relative to their cross section) suggests that current theories, which ignore variations along the magnetic field, are still applicable.

Dyson, P. L.↗

Atmospheric ammonia - Measurements and modeling

Ammonia possesses a unique position in the terrestrial atmosphere in that it is the only gaseous basic constituent. Ammonia readily forms aerosols, and by virtue of its high solubility controls the pH of cloud droplets and precipitation. Over the past year a ground-based solar viewing Infrared Heterodyne Radiometer has been used at Langley Research Center to infer the vertical distribution of ammonia. Ground level in situ measurements of ammonia have also been obtained to supplement the profile data. The ammonia profiles have been analyzed and interpreted with a one-dimensional photochemical model of the troposphere to assess the sources and sinks of NH3.

Hoell, J. M., Jr.↗

The geology and petrology of the Apollo 11 landing site

Geochemical and petrologic data indicate that the 73 Apollo 11 basalts thus far identified can be divided into five petrologic groups (A, B1, B2, B3, D) which must represent at least five separate igneous cooling units. These five igneous bodies range in age from 3.90 b.y. to 3.60 b.y. Photogeologic studies indicate that three mare units are present, and that the lunar module set down on the oldest of the three. The exposure age data suggest that the high-K flow is the surficial rock type at the landing area, and is, therefore, probably the oldest of the three mare units. By examining the size frequency distribution and the inferred cooling rates of the individual samples, it is possible to calculate the formation thicknesses within the 30-m-deep West Crater. This suggests that A = 9 m, B1 = 2 m (and may be an ejecta blanket), B2 is equal to or greater than 8 m, and B3 = 6 m.

Beaty, D. W.↗