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.

SourceResearch connectionFrom yearThrough year

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

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Understanding Differences in Water Adsorption Isotherms: Structural Variations, Force Fields, and Monte Carlo Simulation Approaches

Accurate prediction of water adsorption in micro- and mesoporous materials with hydrophobic pores is essential for the design and characterization of advanced adsorbent materials for separation and energy applications. Here, we assess the reproducibility and consistency of water adsorption isotherms in two microporous all-silica MFI zeolite structures (MFI-K and MFI-O) using two different zeolite force fields and three simulation approaches: grand canonical Monte Carlo (GCMC), Gibbs ensemble Monte Carlo (GEMC), and transition matrix Monte Carlo (TMMC). We demonstrate that consistent treatment of the bulk fluid phase in GCMC and TMMC simulations is critical for reconciling isotherms across methods, and we construct simulation-based equations of state for the TIP4P water model to enable rigorous fugacity-to-pressure conversions. Large shifts in the isotherms are observed for two zeolite force fields developed using different parametrization strategies, with the GCS force field representing implicitly a defect-containing all-silica zeolite, whereas the TraPPE-zeo force field accurately represents an essentially defect-free all-silica zeolite. While water in the van Koningsveld structure of MFI exhibits a first-order phase transition and condensation-like step for adsorption near room temperature, water in the Olson structure of MFI displays continuous adsorption, attributed to differences in the adsorption free energy landscapes. Structural analysis reveals that small geometric variations, particularly Si–O–Si bond angles near the strongest adsorption sites, lead to these substantial differences in adsorption behavior. Furthermore, our results highlight the sensitivity of simulated water adsorption isotherms in hydrophobic frameworks to seemingly small differences in the framework structures, force field parametrization, and simulation approaches.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Consensus and Disagreement in Atmospheric River Detection: ARTMIP Global Catalogues

Atmospheric rivers (ARs) constitute an important mechanism for water vapor transport, but research on their characteristics and impacts has relied on a diverse assortment of detection methodologies, complicating comparisons. The AR Tracking Method Intercomparison Project (ARTMIP) provides a platform for comparing such methodologies, but analysis of ARTMIP catalogues has heretofore focused primarily on specific regions. Here we investigate ARs as detected by an ensemble of algorithms with global coverage. We find that the frequency of occurrence of the majority–consensus ARs produces a robust distribution, featuring five hot spots over the extratropical oceans, against which we compare individual algorithm results. We further explore the underlying similarities and differences via two case studies of AR evolution. The dominant source of disagreement between detection methodologies globally consists of detections (or lack thereof) of weak features, and the algorithms otherwise tend to agree remarkably well on the footprints of ARs.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Coupled Lake‐Atmosphere‐Land Physics Uncertainties in a Great Lakes Regional Climate Model

Abstract This study develops a surrogate‐based method to assess the uncertainty within a convective permitting integrated modeling system of the Great Lakes region, arising from interacting physics parameterizations across the lake, atmosphere, and land surface. Perturbed physics ensembles of the model during the 2018 summer are used to train a neural network surrogate model to predict lake surface temperature (LST) and near‐surface air temperature (T2m). Average physics uncertainties are determined to be 1.5C for LST and T2m over land, and 1.9C for T2m over lake, but these have significant spatiotemporal variations. We find that atmospheric physics parameterizations alone are the dominant sources of uncertainty (45%–53%), while lake and land parameterizations account for 33% and 38% of the uncertainty of LST and T2m over land respectively. Interactions of atmosphere physics parameterizations with those of the land and lake contribute to an additional 13%–17% of the total variance. LST and T2m over the lake are more uncertain in the deeper northern lakes, particularly during the rapid warming phase that occurs in late spring/early summer. The LST uncertainty increases with sensitivity to the lake model's surface wind stress scheme. T2m over land is more uncertain over forested areas in the north, where it is most sensitive to the land surface model, than the more agricultural land in the south, where it is most sensitive to the atmospheric planetary boundary and surface layer scheme. Uncertainty also increases in the southwest during multiday temperature declines with higher sensitivity to the land surface model.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Estimating viscosity of individual substrate-deposited particles from measurements of their height-to-width ratios

Airborne particles alter the radiative forcing of climate and have further consequences on air visibility, atmospheric chemistry, and human health. Recent studies reported the existence of highly viscous semisolid and even solid amorphous organic aerosol (OA) particles. Particle viscosity has an impact on the heterogeneous chemistry, gas-particle partitioning, and ice nucleation properties. Consequently, variations in particle viscosity must be considered when predicting the atmospheric impact of OA. Here, in this work, we use scanning electron microscopy (SEM) and scanning transmission X-ray microscopy (STXM) to estimate the viscosity of individual particles deposited on substrates based on their characteristic height-to-width ratios, which are affected by changes in morphology upon deposition. The height-to-width ratios obtained from SEM and STXM exhibit a strong correlation, demonstrating that both imaging approaches can be applied separately for viscosity assessment of the substrate-deposited particles. While these metrics are largely qualitative, this method enables rapid assessment of particle viscosity ranges, distinguishing between semisolid (>10 10 Pa·s), viscous (10 4– 10 8 Pa·s), and liquid (10°–10 1 Pa·s) particles within ensembles of ambient particles collected for microscopy studies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

TDEP: Temperature Dependent Effective Potentials

The Temperature Dependent Effective Potential (TDEP) method is a versatile and efficient approach to include temperature in ab initio materials simulations based on phonon theory. TDEP can be used to describe thermodynamic properties in classical and quantum ensembles, and several response properties ranging from thermal transport to Neutron and Raman spectroscopy. A stable and fast reference implementation is given in the software package of the same name described here. The underlying theoretical framework and foundation is briefly sketched with an emphasis on discerning the conceptual difference between bare and effective phonon theory, in both self-consistent and non-self-consistent formulations. References to numerous applications and more in-depth discussions of the theory are given.

36 MATERIALS SCIENCE↗

NASA NTRS1984

Direct numerical simulation of reacting flows

The objectives of this work are: (1) to extend the technique of direct numerical simulations to turbulent, chemically reacting flows, (2) to test the validity of the method by comparing computational results with laboratory data, and (3) to use the simulations to gain a better understanding of the effects of turbulence on chemical reactions. The effects of both the large scale structure and the smaller scale turbulence on the overall reaction rates are addressed. The relationship between infinite reaction rate and finite reaction rate chemistry is compared with some of the results of calculations with existing theories and laboratory data. The direct numerical simulation method involves the numerical solution of the detailed evolution of the complex turbulent velocity and concentration fields. Using very efficient numerical methods (e.g., pseudospectral methods), the fully nonlinear (possibly low pass filtered) equations of motion are solved and no closure assumptions or turbulence models are used. Statistical data are obtained by performing spatial, temporal, and/or ensemble averages over the computed flow fields.

Riley, J. J.↗

NASA NTRSDate not supplied

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗

DOE OSTI2020

Transient flow modeling in fractured media using graphs

In this work, we describe a method to simulate transient fluid flows in fractured media using an approach based on graph theory. Our approach builds on past work where the graph-based approach was successfully used to simulate steady-state fluid flows in fractured media. We find a mean computational speedup of the order of 1400 from an ensemble of a 100 discrete fracture networks in contrast to the O(10 4 ) speedup that was obtained for steady-state flows earlier. However, the transient flows considered here involve an additional degree of complexity that was not present in the steady-state flows considered previously with a graph-based approach, that of time marching and solution of the flow equations within a time-stepping scheme. We verify our method with an analytical test case and demonstrate its use on a practical problem related to fluid flows in hydraulically fractured reservoirs. By enabling the study of transient flows, we create an opportunity for a wide set of possibilities where a steady-state approximation is not sufficient, such as the example motivated by hydraulic fracturing that we present here. This work validates the concept that graphs are able to reliably capture the topological properties of the fracture network and serve as effective surrogates in an uncertainty-quantification framework.

58 GEOSCIENCES↗

DOE OSTI2026

An uncertainty visualization framework for large-scale cardiovascular flow simulations: A case study on aortic stenosis

We present a generalizable uncertainty quantification (UQ) and visualization framework for lattice Boltzmann method simulations of high Reynolds number vascular flows, demonstrated on a patient-specific stenosed aorta. The framework combines EasyVVUQ for parameter sampling with large-eddy simulation turbulence modeling in HemeLB, and executes ensembles on the Frontier exascale supercomputer. Spatially resolved metrics, including entropy and isosurface-crossing probability, are used to map uncertainty in pressure and wall shear stress fields directly onto vascular geometries. Two sources of model variability are examined: inlet peak velocity and the Smagorinsky constant. Inlet velocity variation produces high uncertainty downstream of the stenosis where turbulence develops, while upstream regions remain stable. Smagorinsky constant variation has little effect on the large-scale pressure field but increases WSS uncertainty in localized high-shear regions. In both cases, the stenotic throat manifests low entropy, indicative of robust identification of elevated WSS. By linking quantitative UQ measures to three-dimensional anatomy, the framework improves interpretability over conventional 1D UQ plots and supports clinically relevant decision-making, with broad applicability to vascular flow problems requiring both accuracy and spatial insight.

Hemodynamics↗

NASA NTRS2020

Automated Cardiovascular Pathology Assessment using Semantic Segmentation and Ensemble Learning

Cardiac magnetic resonance imaging provides high spatial resolution, enabling improved extraction of important functional and morphological features for cardiovascular disease staging. Segmentation of ventricular cavities and myocardium in cardiac cine sequencing provides a basis to quantify cardiac measures such as ejection fraction. A method is presented that curtails the expense and observer bias of manual cardiac evaluation by combining semantic segmentation and disease classification into a fully automatic processing pipeline. The initial processing element consists of a robust dilated convolutional neural network architecture for voxel-wise segmentation of the myocardium and ventricular cavities. The resulting comprehensive volumetric feature matrix captures diagnostic clinical procedure data and is utilized by the final processing element to model a cardiac pathology classifier. Our approach evaluated anonymized cardiac images from a training data set of 100 patients (4 pathology groups, 1 healthy group, 20 patients per group) examined at the University Hospital of Dijon. The top average Dice index scores achieved were 0.940, 0.886, 0.849 for structure segmentation of the left ventricle (LV), myocardium and right ventricle (RV) respectively. A 5-ary pathology classification accuracy of 90% was recorded on an independent test set using the trained model. Performance results demonstrate potential for advanced machine learning methods to deliver accurate, efficient and reproducible cardiac pathological assessment.

Semantic Segmentation↗

DOE OSTI2021

Comparing Global Spray Combustion Characteristics and Local Shot-to-Shot Variations in a Reacting n -Heptane Spray

In support of efforts to develop improved models of turbulent spray behavior and combustion in diesel engines, experimental data and analysis must be obtained for guidance and validation. For Reynolds-averaged Navier–Stokes (RANS)-based Computational fluid dynamics (CFD) modeling approaches, representative ensemble average experimental results are important. For high-fidelity models such as large eddy simulations (LES)-based CFD, precise individual experimental results are desirable. However, making comparisons between a given experiment and LES is a challenge since local parameters cannot be directly compared. In this work, an optically accessible constant pressure flow rig (CPFR) is utilized to acquire diesel-like fuel injection and reaction behavior simultaneously with three optical diagnostic techniques: rainbow Schlieren deflectometry (RSD), OH* chemiluminescence (OH*), and two-color pyrometry (2CP). The CPFR allows a large number of repeated injection experiments to be performed for statistical analysis and convergence using ensemble-averaging techniques, while maintaining highly repeatable test conditions. Even for stable test conditions, variations in local turbulent fuel–air mixing introduce variability, which manifests as significant differences in OH* and 2CP results. Experimental measurements of characteristic parameters including liquid and vapor jet penetration, liftoff length, soot temperature and concentration, and turbulent flame speed, along with the shot-to-shot variability of each dataset, are presented and discussed. Finally, a statistical method is utilized to analyze the extent of this variability, and to identify superlative injections within the dataset for discussion and analysis of shot-to-shot variations.

42 ENGINEERING↗

DOE OSTI2022

Development of an open-source regional data assimilation system in PEcAn v. 1.7.2: application to carbon cycle reanalysis across the contiguous US using SIPNET

Abstract. The ability to monitor, understand, and predict the dynamics of the terrestrial carbon cycle requires the capacity to robustly and coherently synthesize multiple streams of information that each provide partial information about different pools and fluxes. In this study, we introduce a new terrestrial carbon cycle data assimilation system, built on the PEcAn model–data eco-informatics system, and its application for the development of a proof-of-concept carbon “reanalysis” product that harmonizes carbon pools (leaf, wood, soil) and fluxes (GPP, Ra, Rh, NEE) across the contiguous United States from 1986–2019. We first calibrated this system against plant trait and flux tower net ecosystem exchange (NEE) using a novel emulated hierarchical Bayesian approach. Next, we extended the Tobit–Wishart ensemble filter (TWEnF) state data assimilation (SDA) framework, a generalization of the common ensemble Kalman filter which accounts for censored data and provides a fully Bayesian estimate of model process error, to a regional-scale system with a calibrated localization. Combined with additional workflows for propagating parameter, initial condition, and driver uncertainty, this represents the most complete and robust uncertainty accounting available for terrestrial carbon models. Our initial reanalysis was run on an irregular grid of ∼ 500 points selected using a stratified sampling method to efficiently capture environmental heterogeneity. Remotely sensed observations of aboveground biomass (Landsat LandTrendr) and leaf area index (LAI) (MODIS MOD15) were sequentially assimilated into the SIPNET model. Reanalysis soil carbon, which was indirectly constrained based on modeled covariances, showed general agreement with SoilGrids, an independent soil carbon data product. Reanalysis NEE, which was constrained based on posterior ensemble weights, also showed good agreement with eddy flux tower NEE and reduced root mean square error (RMSE) compared to the calibrated forecast. Ultimately, PEcAn's new open-source regional data assimilation framework provides a scalable workflow for harmonizing multiple data constraints and providing a uniform synthetic platform for carbon monitoring, reporting, and verification (MRV) as well as accelerating terrestrial carbon cycle research.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

A theoretical calculation of stacking fault energy of Ni alloys: The effects of temperature and composition

Combining cluster expansion (CE) method with one dimensional axial Ising model, this work investigated the effects of alloying elements and configurational variations due to temperature on the stacking fault energy (SFE) of FCC Ni binary alloys. Ensembles of large numbers of atomistic structures, each with more than ~400 atoms, were generated to consider sufficient long-range chemical disorder and temperature effects due to configurational entropy. Here, a Monte Carlo Metropolis algorithm was used to generate these structures, whose energies were then evaluated based on the effective cluster interactions obtained from CE. As a baseline, this work had shown the SFE of pure Ni and Al to be 127 mJ/m 2 and 137 mJ/m 2 , respectively, which agreed with the experimental values of 125 mJ/m 2 and 150 mJ/m 2 reported in the literature. Additions of Al, Ti, Cr and Co to pure Ni were found to decrease the SFE to different extents. Although temperature does not strongly influence the SFE of the FCC Ni-Al and Ni-Cr binary alloys, it can lead to significant changes to the SFE of the FCC Ni-Ti and Ni-Co alloys. While effects of temperature and composition on SFE observed in this work were calculated from binary Ni alloys, the general trends are nonetheless expected to be valid in the γ phase of multicomponent Ni superalloys.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Speciation of Nanocatalysts Using X-ray Absorption Spectroscopy Assisted by Machine Learning

Structure and morphology of supported nanoparticle catalysts play important roles in many industrial reactions. Recent progress has identified key aspects of structure-activity relationships at nanoscale and novel methods to study the local environment of the active sites. X-ray absorption fine structure (XAFS) spectroscopy, despite being a leading technique for this purpose, is hampered significantly by its ensemble-averaging nature which often leads to a bias towards a single "representative" structure. Learning heterogeneous distributions of nanostructures at the inter- and intra-particle level from the average XAFS spectrum is a formidable challenge that can be overcome in some cases, described in this Perspective. Here we also discuss emerging machine learning techniques for extracting the information about heterogeneity of metal species from XAFS data.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Probing boron vacancy defects in hBN via single spin relaxometry

Spin defects in solids offer promising platforms for quantum sensing and memory due to their long coherence times and optical addressability. Here, we integrate a single nitrogen-vacancy (NV) center in diamond with scanning probe microscopy to detect, read out, and spatially map spin-based quantum sensors at the nanoscale. Using the boron vacancy ($V$$^{–}_{B}$) center in hexagonal boron nitride—an emerging two-dimensional spin system—as a model, we detect its electron spin resonance indirectly via changes in the spin relaxation time (T 1 ) of a nearby NV center, eliminating the need for optical excitation or fluorescence detection of the $V$$^{–}_{B}$. Cross-relaxation between NV and $V$$^{–}_{B}$ ensembles significantly reduces NV T1, enabling quantitative nanoscale mapping of defect densities beyond the optical diffraction limit and clear resolution of hyperfine splitting in isotopically enriched h 10 B 15 N. Our method demonstrates interactions between spin sensors in 3D and 2D materials, establishing NV centers as versatile probes for characterizing otherwise inaccessible spin defects.

Quantum metrology↗

DOE OSTI2020

A global ensemble of ocean wave climate projections from CMIP5-driven models

This dataset, produced through the Coordinated Ocean Wave Climate Project (COWCLIP) phase 2, represents the first coordinated multivariate ensemble of 21st Century global wind-wave climate projections available (henceforth COWCLIP2.0). COWCLIP2.0 comprises general and extreme statistics of significant wave height (H S ), mean wave period (T m ), and mean wave direction (θ m ) computed over time-slices 1979–2004 and 2081–2100, at different frequency resolutions (monthly, seasonally and annually). The full ensemble comprising 155 global wave climate simulations is obtained from ten CMIP5-based state-of-the-art wave climate studies and provides data derived from alternative wind-wave downscaling methods, and different climate-model forcing and future emissions scenarios. The data has been produced, and processed, under a specific framework for consistency and quality, and follows CMIP5 Data Reference Syntax, Directory structures, and Metadata requirements. Technical comparison of model skill against 26 years of global satellite measurements of significant wave height has been undertaken at global and regional scales. This new dataset provides support for future broad scale coastal hazard and vulnerability assessments and climate adaptation studies in many offshore and coastal engineering applications.

Science & Technology - Other Topics↗

NASA NTRS1996

The Calculation of Accurate Harmonic Frequencies of Large Molecules: The Polycyclic Aromatic Hydrocarbons, a Case Study

The vibrational frequencies and infrared intensities of naphthalene neutral and cation are studied at the self-consistent-field (SCF), second-order Moller-Plesset (MP2), and density functional theory (DFT) levels using a variety of one-particle basis sets. Very accurate frequencies can be obtained at the DFT level in conjunction with large basis sets if they are scaled with two factors, one for the C-H stretches and a second for all other modes. We also find remarkably good agreement at the B3LYP/4-31G level using only one scale factor. Unlike the neutral PAHs where all methods do reasonably well for the intensities, only the DFT results are accurate for the PAH cations. The failure of the SCF and MP2 methods is caused by symmetry breaking and an inability to describe charge delocalization. We present several interesting cases of symmetry breaking in this study. An assessment is made as to whether an ensemble of PAH neutrals or cations could account for the unidentified infrared bands observed in many astronomical sources.

Bauschlicher, Charles W., Jr.↗

NASA NTRS2009

Loss Factor Estimation Using the Impulse Response Decay Method on a Stiffened Structure

High-frequency vibroacoustic modeling is typically performed using energy-based techniques such as Statistical Energy Analysis (SEA). Energy models require an estimate of the internal damping loss factor. Unfortunately, the loss factor is difficult to estimate analytically, and experimental methods such as the power injection method can require extensive measurements over the structure of interest. This paper discusses the implications of estimating damping loss factors using the impulse response decay method (IRDM) from a limited set of response measurements. An automated procedure for implementing IRDM is described and then evaluated using data from a finite element model of a stiffened, curved panel. Estimated loss factors are compared with loss factors computed using a power injection method and a manual curve fit. The paper discusses the sensitivity of the IRDM loss factor estimates to damping of connected subsystems and the number and location of points in the measurement ensemble.

Cabell, Randolph↗