Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble methods”

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 937 records · Page 52

Satellite-Assisted Particulate Matter (SAPM) for the Models, In situ, and Remote Sensing of Aerosols (MIRA) Working Group

Models, In situ, and Remote sensing of Aerosols (MIRA) is an international working group that provides a forum for collaborations amongst these three atmospheric aerosol communities. MIRA consists of a collection of interdisciplinary and independently funded Topic/Project groups with clear goals and generally characterized by requests for additional scientific data. The Satellite-Assisted Particulate Matter (SAPM) Topic group focuses on the study of fine particulate matter (PM2.5), as it is a major contributor to air pollution and negatively impacts human health. Our group aims to provide intercomparisons of various methods and techniques for estimating surface PM2.5 assisted by satellite remote sensors (passive and active), global aerosol models, and in situ aerosol measurements. The overall motivation of our work is to enhance PM2.5 coverage beyond in situ ground stations, which can be limited. In this presentation, we provide an overview of various PM2.5 estimation approaches by current SAPM team members, as each approach has its own strengths and limitations. One approach uses near-surface aerosol extinction retrievals from the spaceborne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and an assumed value of mass extinction efficiency to derive PM2.5 concentrations (Toth et al 2019). A study using this method found a promising agreement (R=0.60; slope=0.89) with U.S. Environmental Protection Agency (EPA) in situ PM2.5 measurements over the contiguous United States (CONUS) for a 12-year (2007-2018) period (Toth et al 2022). In a different approach combining spaceborne lidar and a model, the Cloud Aerosol Transport System (CATS) lidar and Goddard Earth Observing System (GEOS) model are blended in an ensemble variational assimilation scheme to retrieve aerosol extinction and convert speciated mass concentrations to total PM2.5. A recent study (Matus et al in prep.) found agreement within 2 μg/m3 over the CONUS on the 2016 annual mean when compared to EPA ground-based measurements. Another study explored trends in city aerosol optical depth (AOD) from a spaceborne passive remote sensor, Moderate Resolution Imaging Spectroradiometer (MODIS), and found that these trends agree well with trends in city surface PM2.5, which provides confidence in the use of city AOD for city PM2.5 trend studies (Vohra et al 2021). In a follow-on study, steep and significant trends (~2.5 to 7.8% a-1) in AOD are found for several South Asian cities (e.g., Bangalore and Hyderabad), suggesting rapid growth in PM2.5 (Vohra et al 2022). Also in this presentation, we report results on our current collective SAPM study of India, a country that exhibits high levels of PM2.5 pollution. We show the temporal (i.e., seasonal) and spatial variability of PM2.5 over India using the CALIOP and CATS+model approaches, comparisons between the estimated PM2.5 from the retrieval schemes, and validation with in situ PM2.5 ground-based observations from the Central Pollution Control Board in India. This region provides us with an excellent case study to test our PM2.5 retrievals in heavily polluted scenes and over a complex and varying topography. The SAPM Topic group actively seeks international participants/collaborators in our group, including those working with in situ aerosol measurements (e.g., ACTRIS). We are interested in aerosol datasets in order to either improve our PM2.5 estimates from various approaches (e.g., using in situ mass scattering/absorption coefficient and aerosol hygroscopic properties for various aerosol species) and/or validate the PM2.5 estimates (e.g., using in situ ground-based PM2.5 concentrations).

Travis D Toth↗

Calabi-Yau CFTs and random matrices

Using numerical methods for finding Ricci-flat metrics, we explore the spectrum of local operators in two-dimensional conformal field theories defined by sigma models on Calabi-Yau targets at large volume. Focusing on the examples of K3 and the quintic, we show that the spectrum, averaged over a region in complex structure moduli space, possesses the same statistical properties as the Gaussian orthogonal ensemble of random matrix theory.

superstring vacua↗

Inferring the Redshift of More than 150 GRBs with a Machine-learning Ensemble Model

Gamma-ray bursts (GRBs), due to their high luminosities, are detected up to a redshift of 10, and thus have the potential to be vital cosmological probes of early processes in the Universe. Fulfilling this potential requires a large sample of GRBs with known redshifts, but due to observational limitations, only 11% have known redshifts (z). There have been numerous attempts to estimate redshifts via correlation studies, most of which have led to inaccurate predictions. To overcome this, we estimated GRB redshift via an ensemble-supervised machine-learning (ML) model that uses X-ray afterglows of long-duration GRBs observed by the Neil Gehrels Swift Observatory. The estimated redshifts are strongly correlated (a Pearson coefficient of 0.93) and have an rms error, namely, the square root of the average squared error $\langle$Δz 2 $\rangle$, of 0.46 with the observed redshifts showing the reliability of this method. The addition of GRB afterglow parameters improves the predictions considerably by 63% compared to previous results in peer-reviewed literature. Finally, we use our ML model to infer the redshifts of 154 GRBs, which increase the known redshifts of long GRBs with plateaus by 94%, a significant milestone for enhancing GRB population studies that require large samples with redshift.

79 ASTRONOMY AND ASTROPHYSICS↗

Multiclass Flight Anomaly Detection Using Sensor Fusion Based on Dempster-Shafer Theory

As aviation systems in commercial operations continue to grow in complexity, the anomalies exhibited by these systems become more elaborate and difficult to detect. To address the challenge of detecting these complex anomalies, deep learning models have been used extensively in aviation anomaly detection studies, at the expense of end-user interpretability. Aiming to maintain the same level of interpretability as traditional threshold-exceedance methods, we continue our development of prediction models using ordinal patterns and their distributions throughout the flight. Specifically, this study extends our work into multiclass anomaly detection using sensor fusion based on Dempster-Shafer theory (DST), a second-order probability theory used to combine information from different sources of evidence. Our approach uses DST toreduce the uncertainty in the class predictions of an ensemble of classifiers. These classifiers rely on the similarity between flight data and class templates to make a prediction of the state of the aircraft. Our approach aims to take advantage of simple models trained on interpretable features (ordinal patterns) to correctly predict an anomaly and identify the flight dynamics linked to the anomaly. Our results show an improvement when using DST-based sensor fusion over simple majority voting. Additionally, our results provide insight into aircraft states linked to rare high-risk anomalies.

Risk detection↗

Semiclassical shell-structure micro-macroscopic approach for the level density

Level density ρ(E,A) is derived for a one-component nucleon system with a given energy E and particle number A within the mean-field semiclassical periodic-orbit theory beyond the saddle-point method of the Fermi gas model. We obtain ρ∝I ν (S)/S ν , with I ν (S) being the modified Bessel function of the entropy S. Within the micro-macro-canonical approximation (MMA), for a small thermal excitation energy U, with respect to rotational excitations E rot , one obtains ν = 3/2 for ρ(E,A). In the case of excitation energy U larger than E rot but smaller than the neutron separation energy, one finds a larger value of ν = 5/2. A role of the fixed spin variables for rotating nuclei is discussed. The MMA level density ρ reaches the well-known grand-canonical ensemble limit (Fermi gas asymptote) for large S related to large excitation energies, and also reaches the finite micro-canonical limit for small combinatorial entropy S at low excitation energies (the constant “temperature” model). Fitting the ρ(E,A) of the MMA to the experimental data for low excitation energies, taking into account shell and, qualitatively, pairing effects, one obtains for the inverse level density parameter K a value which differs essentially from that parameter derived from data on neutron resonances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Alchemical Free Energy Estimators and Molecular Dynamics Engines: Accuracy, Precision, and Reproducibility

The binding free energy between a ligand and its target protein is an essential quantity to know at all stages of the drug discovery pipeline. Assessing this value computationally can offer insight into where efforts should be focused in the pursuit of effective therapeutics to treat a myriad of diseases. In this work, we examine the computation of alchemical relative binding free energies with an eye for assessing reproducibility across popular molecular dynamics packages and free energy estimators. The focus of this work is on 54 ligand transformations from a diverse set of protein targets: MCL1, PTP1B, TYK2, CDK2, and thrombin. These targets are studied with three popular molecular dynamics packages: OpenMM, NAMD2, and NAMD3 alpha. Trajectories collected with these packages are used to compare relative binding free energies calculated with thermodynamic integration and free energy perturbation methods. The resulting binding free energies show good agreement between molecular dynamics packages with an average mean unsigned error between them of 0.50 kcal/mol. The correlation between packages is very good, with the lowest Spearman’s, Pearson’s and Kendall’s tau correlation coefficients being 0.92, 0.91, and 0.76, respectively. Agreement between thermodynamic integration and free energy perturbation is shown to be very good when using ensemble averaging.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Roadmap on nanoscale magnetic resonance imaging

Abstract The field of nanoscale magnetic resonance imaging (NanoMRI) was started 30 years ago. It was motivated by the desire to image single molecules and molecular assemblies, such as proteins and virus particles, with near-atomic spatial resolution and on a length scale of 100 nm. Over the years, the NanoMRI field has also expanded to include the goal of useful high-resolution nuclear magnetic resonance (NMR) spectroscopy of molecules under ambient conditions, including samples up to the micron-scale. The realization of these goals requires the development of spin detection techniques that are many orders of magnitude more sensitive than conventional NMR and MRI, capable of detecting and controlling nanoscale ensembles of spins. Over the years, a number of different technical approaches to NanoMRI have emerged, each possessing a distinct set of capabilities for basic and applied areas of science. The goal of this roadmap article is to report the current state of the art in NanoMRI technologies, outline the areas where they are poised to have impact, identify the challenges that lie ahead, and propose methods to meet these challenges. This roadmap also shows how developments in NanoMRI techniques can lead to breakthroughs in emerging quantum science and technology applications.

Budakian, Raffi↗

Multiclass Flight Anomaly Detection Using Sensor Fusion Based on Dempster-Shafer Theory

As aviation systems in commercial operations continue to grow in complexity, the anomalies exhibited by these systems become more elaborate and difficult to detect. To address the challenge of detecting these complex anomalies, deep learning models have been used extensively in aviation anomaly detection studies, at the expense of end-user interpretability. Aiming to maintain the same level of interpretability as traditional threshold-exceedance methods, we continue our development of prediction models using ordinal patterns and their distributions throughout the flight. Specifically, this study extends our work into multiclass anomaly detection using sensor fusion based on Dempster-Shafer theory (DST), a second-order probability theory used to combine information from different sources of evidence. Our approach uses DST to reduce the uncertainty in the class predictions of an ensemble of classifiers. These classifiers rely on the similarity between flight data and class templates to make a prediction of the state of the aircraft. Our approach aims to take advantage of simple models trained on interpretable features (ordinal patterns) to correctly predict an anomaly and identify the flight dynamics linked to the anomaly. Our results show an improvement when using DST-based sensor fusion over simple majority voting. Additionally, our results provide insight into aircraft states linked to rare high-risk anomalies.

Risk detection↗

Single-molecule epitranscriptomic analysis of full-length HIV-1 RNAs reveals functional roles of site-specific m6As

Abstract Although the significance of chemical modifications on RNA is acknowledged, the evolutionary benefits and specific roles in human immunodeficiency virus (HIV-1) replication remain elusive. Most studies have provided only population-averaged values of modifications for fragmented RNAs at low resolution and have relied on indirect analyses of phenotypic effects by perturbing host effectors. Here we analysed chemical modifications on HIV-1 RNAs at the full-length, single RNA level and nucleotide resolution using direct RNA sequencing methods. Our data reveal an unexpectedly simple HIV-1 modification landscape, highlighting three predominant N 6 -methyladenosine (m 6 A) modifications near the 3′ end. More densely installed in spliced viral messenger RNAs than in genomic RNAs, these m 6 As play a crucial role in maintaining normal levels of HIV-1 RNA splicing and translation. HIV-1 generates diverse RNA subspecies with distinct m 6 A ensembles, and maintaining multiple of these m 6 As on its RNAs provides additional stability and resilience to HIV-1 replication, suggesting an unexplored viral RNA-level evolutionary strategy.

60 APPLIED LIFE SCIENCES↗

Data and code for Daily and Multi-Day Extreme Rainfall Analysis Under Future Climates Using Stochastic Storm Transposition and NEX-GDDP-CMIP6 Over CONUS

This data package provides inputs, codes, and outputs for a comprehensive analysis of projected changes in extreme precipitation across 10 regions of the continental United States, using 34 downscaled Earth System Models (ESMs) from the NASA Earth Exchange Global Daily Downscaled Projections, Coupled Model Intercomparison Project Phase 6 (NEX-GDDP-CMIP6) dataset. These models are part of the Coupled Model Intercomparison Project Phase 6 (CMIP6), a coordinated climate modeling framework widely used to assess climate change impacts. The analysis applies a stochastic storm transposition method to quantify changes in extreme rainfall under two Shared Socioeconomic Pathway (SSP) climate scenarios—SSP2-4.5 (moderate emissions) and SSP5-8.5 (high emissions)—compared to historical conditions (1995–2014 vs. 2081–2100). The dataset includes rainfall depth estimates for extreme events with return periods from 2 to 500 years across multiple storm durations (1, 3, and 5 days) for each of the 10 U.S. regions. Weighted ensemble statistics are derived from individual ESM performance against historical precipitation patterns, enabling robust uncertainty quantification through both sign-based and permutation-test-based model agreement assessments. Key analyses address: (1) relative changes in extreme precipitation for each climate scenario, (2) differences between SSP scenarios (SSP5-8.5 vs. SSP2-4.5), (3) contrasts between rare and frequent events, and (4) variations between multi-day and daily storm durations. The workflow produces ensemble statistics—median, 5th, 25th, 75th, and 95th percentiles—along with model agreement metrics that identify regions and event types with robust climate change signals. The dataset includes: processed rainfall depth outputs (netCDF format) from the RainyDay Python package, ESM weights from historical performance evaluation using DayMet observations, ensemble statistics across all storm dimensions, and figures summarizing key findings.

54 ENVIRONMENTAL SCIENCES↗

A machine learning approach to emulation and biophysical parameter estimation with the Community Land Model, version 5

Abstract. Land models are essential tools for understanding and predicting terrestrial processes and climate–carbon feedbacks in the Earth system, but uncertainties in their future projections are poorly understood. Improvements in physical process realism and the representation of human influence arguably make models more comparable to reality but also increase the degrees of freedom in model configuration, leading to increased parametric uncertainty in projections. In this work we design and implement a machine learning approach to globally calibrate a subset of the parameters of the Community Land Model, version 5 (CLM5) to observations of carbon and water fluxes. We focus on parameters controlling biophysical features such as surface energy balance, hydrology, and carbon uptake. We first use parameter sensitivity simulations and a combination of objective metrics including ranked global mean sensitivity to multiple output variables and non-overlapping spatial pattern responses between parameters to narrow the parameter space and determine a subset of important CLM5 biophysical parameters for further analysis. Using a perturbed parameter ensemble, we then train a series of artificial feed-forward neural networks to emulate CLM5 output given parameter values as input. We use annual mean globally aggregated spatial variability in carbon and water fluxes as our emulation and calibration targets. Validation and out-of-sample tests are used to assess the predictive skill of the networks, and we utilize permutation feature importance and partial dependence methods to better interpret the results. The trained networks are then used to estimate global optimal parameter values with greater computational efficiency than achieved by hand tuning efforts and increased spatial scale relative to previous studies optimizing at a single site. By developing this methodology, our framework can help quantify the contribution of parameter uncertainty to overall uncertainty in land model projections.

54 ENVIRONMENTAL SCIENCES↗

Assessing Global and Local Radiative Feedbacks Based on AGCM Simulations for 1980–2014/2017

In order to avoid contamination of unobservable and uncertain effective radiative forcing (ERF) on the diagnosis of radiative feedbacks based on short-term climate variability, we verify the Kernel-Gregory feedback calculation method using atmospheric model experiments with prescribed ERF. We show that both clear-sky radiative fluxes and all-sky radiative feedbacks have a closure between model simulations and kernel derivations. A near-zero global mean net cloud feedback is found in the simulation with prescribed ERF, which results in a more negative global net climate feedback, -2 W m -2 K -1 . Consistent with AMIP6 ensemble mean results, the lapse rate feedback is the largest contributor among all feedbacks to temperature amplification over the three poles (Arctic, Antarctic and Tibetan Plateau), followed by surface albedo feedback and Planck feedback deviation from its global mean. Except for the higher surface albedo feedback in the Antarctic, other feedbacks are almost same between the Arctic and Antarctic.

54 ENVIRONMENTAL SCIENCES↗

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg (1997). This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple- ζ basis set predicts the spin splittings of [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ to be –35.7 and –17.1 kcal mol –1 , respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A procedure for combining acoustically induced and mechanically induced loads (first passage failure design criterion)

The combined load statistics are developed by taking the acoustically induced load to be a random population, assumed to be stationary. Each element of this ensemble of acoustically induced loads is assumed to have the same power spectral density (PSD), obtained previously from a random response analysis employing the given acoustic field in the STS cargo bay as a stationary random excitation. The mechanically induced load is treated as either (1) a known deterministic transient, or (2) a nonstationary random variable of known first and second statistical moments which vary with time. A method is then shown for determining the probability that the combined load would, at any time, have a value equal to or less than a certain level. Having obtained a statistical representation of how the acoustic and mechanical loads are expected to combine, an analytical approximation for defining design levels for these loads is presented using the First Passage failure criterion.

Crowe, D. R.↗

Ensemble Kalman filter for data assimilation coupled with low-resolution computations techniques applied in fluid dynamics

This paper presents an innovative Reduced-order model (ROM) for merging experimental and simulation data using data assimilation (DA) to estimate the "True" state of a fluid dynamics system, leading to more accurate predictions. Our methodology introduces a novel approach by implementing the ensemble Kalman filter (EnKF) within a reduced-dimensional framework, grounded in a robust theoretical foundation and applied to fluid dynamics. To address the substantial computational demands of DA, the proposed ROM employs low-resolution (LR) techniques to drastically reduce computational costs. This innovative approach involves downsampling datasets for DA computations, followed by an advanced reconstruction technique based on low-cost singular value decomposition (lcSVD). The lcSVD method, a key innovation in this paper, has never been applied to DA before and offers a highly efficient way to enhance resolution with minimal computational resources. Our results demonstrate significant reductions in both computation time and RAM usage through these LR techniques without compromising the accuracy of the estimations. For instance, in a turbulent test case, for a data compression rate of 15.9, the LR approach can achieve a speed-up of 13.7 and a RAM compression of 90.9% while maintaining a low relative root mean square error (RRMSE) of 2.6%, compared to 0.8% in the high-resolution (HR) reference. Furthermore, we highlight the effectiveness of the EnKF in estimating and predicting the state of fluid flow systems based on limited observations and given low-fidelity numerical data. This paper highlights the potential of the proposed DA method in fluid dynamics applications, particularly for improving computational efficiency in CFD and related fields. Its ability to balance accuracy with low computational and memory costs makes it especially suitable for large-scale and real-time applications, such as environmental monitoring or engineering design. This method will be incorporated into ModelFLOWs-app.

Data Assimilation↗

National CO 2 budgets (2015–2020) inferred from atmospheric CO 2 observations in support of the global stocktake

Accurate accounting of emissions and removals of CO 2 is critical for the planning and verification of emission reduction targets in support of the Paris Agreement. Here, we present a pilot dataset of country-specific net carbon exchange (NCE; fossil plus terrestrial ecosystem fluxes) and terrestrial carbon stock changes aimed at informing countries' carbon budgets. These estimates are based on “top-down” NCE outputs from the v10 Orbiting Carbon Observatory (OCO-2) modeling intercomparison project (MIP), wherein an ensemble of inverse modeling groups conducted standardized experiments assimilating OCO-2 column-averaged dry-air mole fraction (X CO 2 ) retrievals (ACOS v10), in situ CO 2 measurements or combinations of these data. The v10 OCO-2 MIP NCE estimates are combined with “bottom-up” estimates of fossil fuel emissions and lateral carbon fluxes to estimate changes in terrestrial carbon stocks, which are impacted by anthropogenic and natural drivers. These flux and stock change estimates are reported annually (2015–2020) as both a global 1° x 1° gridded dataset and a country-level dataset and are available for download from the Committee on Earth Observation Satellites' (CEOS) website. Across the v10 OCO-2 MIP experiments, we obtain increases in the ensemble median terrestrial carbon stocks of 3.29–4.58 Pg CO2 yr -1 (0.90–1.25 Pg C yr -1 ). This is a result of broad increases in terrestrial carbon stocks across the northern extratropics, while the tropics generally have stock losses but with considerable regional variability and differences between v10 OCO-2 MIP experiments. We discuss the state of the science for tracking emissions and removals using top-down methods, including current limitations and future developments towards top-down monitoring and verification systems.

54 ENVIRONMENTAL SCIENCES↗

A benchmark dataset for Hydrogen Combustion

The generation of reference data for deep learning models is challenging for reactive systems, and more so for combustion reactions due to the extreme conditions that create radical species and alternative spin states during the combustion process. Here, we extend intrinsic reaction coordinate (IRC) calculations with ab initio MD simulations and normal mode displacement calculations to more extensively cover the potential energy surface for 19 reaction channels for hydrogen combustion. A total of ~290,000 potential energies and ~1,270,000 nuclear force vectors are evaluated with a high quality range-separated hybrid density functional, ωB97X-V, to construct the reference data set, including transition state ensembles, for the deep learning models to study hydrogen combustion reaction.

08 HYDROGEN↗

Reversible Surface Energy Storage in Molecular-Scale Porous Materials

Forcible wetting of hydrophobic pores represents a viable method for energy storage in the form of interfacial energy. The energy used to fill the pores can be recovered as pressure–volume work upon decompression. For efficient recovery, the expulsion pressure should not be significantly lower than the pressure required for infiltration. Hysteresis of the wetting/drying cycle associated with the kinetic barrier to liquid expulsion results in energy dissipation and reduced storage efficiency. In the present work, we use open ensemble (Grand Canonical) Monte Carlo simulations to study the improvement of energy recovery with decreasing diameters of planar pores. Near-complete reversibility is achieved at pore widths barely accommodating a monolayer of the liquid, thus minimizing the area of the liquid/gas interface during the cavitation process. At the same time, these conditions lead to a steep increase in the infiltration pressure required to overcome steric wall/water repulsion in a tight confinement and a considerable reduction in the translational entropy of confined molecules. In principle, similar effects can be expected when increasing the size of the liquid particles without altering the absorbent porosity. While the latter approach is easier to follow in laboratory work, we discuss the advantages of reducing the pore diameter, which reduces the cycling hysteresis while simultaneously improving the stored-energy density in the material.

36 MATERIALS SCIENCE↗