Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “benchmark correlations”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Thermal mean-field theories

Several closely related ab initio thermal mean-field theories for fermions, both well-established and new ones, are compared with one another at the formalism level and numerically. The theories considered are Fermi–Dirac theory; thermal Hartree–Fock (HF) theory; two modifications of the thermal single-determinant and the first-order finite-temperature many-body perturbation theory based on a zero-temperature or thermal HF reference. Furthermore, thermal full-configuration-interaction theory is used as the benchmark.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2019

A Segmentation Algorithm for Characterizing Rise and Fall Segments in Seasonal Cycles: an Application to Xco2 to Estimate Benchmarks and Assess Model Bias

There is more useful information in the time series of satellite-derived column-averaged carbon dioxide (XCO2) than is typically characterized. Often, the entire time series is treated at once without considering detailed features at shorter timescales, such as nonstationary changes in signal characteristics – amplitude, period and phase. In many instances, signals are visually and analytically differentiable from other portions in a time series. Each rise (increasing) and fall (decreasing) segment in the seasonal cycle is visually discernable in a graph of the time series. The rise and fall segments largely result from seasonal differences in terrestrial ecosystem production, which means that the segment's signal characteristics can be used to establish observational benchmarks because the signal characteristics are driven by similar underlying processes. We developed an analytical segmentation algorithm to characterize the rise and fall segments in XCO2 seasonal cycles. We present the algorithm for general application of the segmentation analysis and emphasize here that the segmentation analysis is more generally applicable to cyclic time series. We demonstrate the utility of the algorithm with specific results related to the comparison between satellite- and model-derived XCO2 seasonal cycles (2009–2012) for large bioregions across the globe. We found a seasonal amplitude gradient of 0.74–0.77 ppm for every 10∘ of latitude in the satellite data, with similar gradients for rise and fall segments. This translates to a south–north seasonal amplitude gradient of 8 ppm for XCO2, about half the gradient in seasonal amplitude based on surface site in situ CO2 data (∼19 ppm). The latitudinal gradients in the period of the satellite-derived seasonal cycles were of opposing sign and magnitude (−9 d per 10∘ latitude for fall segments and 10 d per 10∘ latitude for rise segments) and suggest that a specific latitude (∼2∘ N) exists that defines an inversion point for the period asymmetry. Before (after) the point of asymmetry inversion, the periods of rise segments are lesser (greater) than the periods of fall segments; only a single model could reproduce this emergent pattern. The asymmetry in amplitude and the period between rise and fall segments introduces a novel pattern in seasonal cycle analyses, but, while we show these emergent patterns exist in the data, we are still breaking ground in applying the information for science applications. Maybe the most useful application is that the segmentation analysis allowed us to decompose the model biases into their correlated parts of biases in amplitude, period and phase independently for rise and fall segments. We offer an extended discussion on how such information about model biases and the emergent patterns in satellite-derived seasonal cycles can be used to guide future inquiry and model development.

segmentation algorithm↗
DOE OSTI2025

Desmearing small-angle scattering data by central moment expansions of instrument resolution

This study introduces an enhanced numerical technique tailored specifically for refining 1D small-angle scattering (SAS) intensity profiles affected by smearing. Our primary objective is to address the resolution blurring commonly encountered in SAS data, particularly in systems with clearly defined correlation peaks whose spread aligns with the width of the resolution function at corresponding Q positions. Unlike previous approaches that expanded the SAS intensity using central moment expansion, the new method focuses on expanding the resolution function itself, thus eliminating artificial oscillations observed in smeared spectra due to limitations inherent in our earlier algorithm. This method is straightforward to implement, computationally efficient and consistently performs well in numerical benchmarking. To illustrate its effectiveness, we present a case study of a lamellar phase characterized by distinct peaks in its small-angle neutron scattering intensities.

central moment expansions↗
DOE OSTI2020

Metallic Fuel Benchmark Simulations Based on the X430 Experiments

To gauge the accuracy of current BISON metallic fuel performance models and identify models that require further development, Oak Ridge National Laboratory has developed benchmark problems based on two pins of the X430 experiment. These benchmarks recreate the experimental conditions of the X430 series of experiments with as much accuracy as possible given the available data. The benchmarks were implemented in the BISON nuclear fuel performance code and the simulation results were compared with the experimental results to quantify the simulation accuracy. Findings in this work were consistent with previous benchmark problems. The temperature profiles, burnup, and fission gas behavior were all as accurate as can reasonably be expected. The cladding radial growth was accurate enough to not affect the results of other measurements. The fuel axial growth was accurate for U-Zr fuel but not for U-Pu-Zr fuel, which appears to be caused by a lack of U-Pu-Zr data. Swelling correlations have not been developed for U-Pu-Zr fuel because there were not enough data on this fuel type, so BISON used the U-Zr correlations. Finally, recommendations were made for near-term BISON development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Electronic structure analysis and DFT benchmarking of Rydberg-type alkali-metal-crown ether, -cryptand, and -adamanzane complexes

Density functional theory (DFT) and electron propagator theory (EPT) calculations were performed to study ground and excited electronic structures of alkali-metal (M) coordinated 9-crown-3, 24-crown-8, [2.1.1]cryptand, o-Me 2 -1.1.1, and 3 6 Adamanzane complexes. Each complex bears an expanded electron in the periphery and occupies diffuse 1p-, 1d-, 1f-type molecular orbitals (or superatomic 1P, 1D, 1F orbitals) in excited electronic states. The calculated superatomic shell model of the M(9-crown-3) 2 is 1S, 1P, 1D, 1F, 2S, 2P, 2D, 1G and it is held by all other complexes up to the studied 1F level. Due to the highly diffuse nature of the electron, the ionization energies of these complexes are significantly lower (1.6–2.0 eV) and hence these complexes belong to the superalkali category. Further, the ab initio EPT ionization energy and the excitation energies of the Li(9-crown-3) 2 were used to evaluate DFT errors associated with a series of exchange correlation functionals that span multiple rungs of Jacob's ladder (i.e., GGA, meta-GGA, global GGA hybrid, meta-GGA hybrid, range-separated hybrid, double-hybrid). Among these, the best performing functional is the range-separated hybrid CAM-B3LYP and the errors are within 6% of high-level ab initio EPT results. The accuracy of CAM-B3LYP is indeed transferable to similar complexes and hence the findings are expected to accelerate the progression of studies of Rydberg-type systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Development of machine learning framework for interface force closures based on bubble tracking data

Interfacial force closures in the two-fluid model play a critical role for the predictive capabilities of void fraction distribution. However, the practices of interfacial force modeling have long been challenged by the inherent physical complexity of the two-phase flows. The rapidly expanding computational capabilities in the recent years have made high-fidelity data from the interface-captured direct numerical simulation become more available, and hence potential for data-driven interfacial force modeling has prevailed. In this work, we established a data-driven modeling framework integrated to the HZDR multiphase Eulerian-Eulerian framework for computational fluid dynamics simulations. The data-driven framework is verified in a benchmark problem, where a feedforward neural network managed to capture the non-linear mapping between bubble Reynolds number and drag coefficient and reproduce the void distribution resulting from the baseline model in the test case. The second focus is on utilizing the bubble tracking data set to form a closure for the bubble drag in the turbulent bubbly flow, in which the drag coefficient is set to be correlated with the bubble Reynolds number and the Eötvös number. Pseudo-steady state filtering in the Frenet Frame was carried out to obtain the drag coefficient from the turbulent bubbly flow data. The performance of the data-driven drag model is also examined through a case study, where improvement of model’s prediction near-wall is regarded necessary. In conclusion, discussion and further plans of investigation are provided.

42 ENGINEERING↗
DOE OSTI2023

A data-driven method for modelling dissipation rates in stratified turbulence

We present a deep probabilistic convolutional neural network (PCNN) model for predicting local values of small-scale mixing properties in stratified turbulent flows, namely the dissipation rates of turbulent kinetic energy and density variance, $\varepsilon$ and $\chi$ . Inputs to the PCNN are vertical columns of velocity and density gradients, motivated by data typically available from microstructure profilers in the ocean. The architecture is designed to enable the model to capture several characteristic features of stratified turbulence, in particular the dependence of small-scale isotropy on the buoyancy Reynolds number $Re_b:=\varepsilon /(\nu N^2)$ , where $\nu$ is the kinematic viscosity and $N$ is the background buoyancy frequency, the correlation between suitably locally averaged density gradients and turbulence intensity and the importance of capturing the tails of the probability distribution functions of values of dissipation. Empirically modified versions of commonly used isotropic models for $\varepsilon$ and $\chi$ that depend only on vertical derivatives of density and velocity are proposed based on the asymptotic regimes $Re_b\ll 1$ and $Re_b\gg 1$ , and serve as an instructive benchmark for comparison with the data-driven approach. When trained and tested on a simulation of stratified decaying turbulence which accesses a range of turbulent regimes (associated with differing values of $Re_b$ ), the PCNN outperforms assumptions of isotropy significantly as $Re_b$ decreases, and additionally demonstrates improvements over the fitted empirical models. A differential sensitivity analysis of the PCNN facilitates a comparison with the theoretical models and provides a physical interpretation of the features enabling it to make improved predictions.

42 ENGINEERING↗
DOE OSTI2023

Scalar Breit interaction for molecular calculations

Variational treatment of the Dirac–Coulomb–Gaunt or Dirac–Coulomb–Breit two-electron interaction at the Dirac–Hartree–Fock level is the starting point of high-accuracy four-component calculations of atomic and molecular systems. In this work, we introduce, for the first time, the scalar Hamiltonians derived from the Dirac–Coulomb–Gaunt and Dirac–Coulomb–Breit operators based on spin separation in the Pauli quaternion basis. While the widely used spin-free Dirac–Coulomb Hamiltonian includes only the direct Coulomb and exchange terms that resemble nonrelativistic two-electron interactions, the scalar Gaunt operator adds a scalar spin–spin term. The spin separation of the gauge operator gives rise to an additional scalar orbit-orbit interaction in the scalar Breit Hamiltonian. Benchmark calculations of Au n (n = 2–8) show that the scalar Dirac–Coulomb–Breit Hamiltonian can capture 99.99% of the total energy with only 10% of the computational cost when real-valued arithmetic is used, compared to the full Dirac–Coulomb–Breit Hamiltonian. In conclusion, the scalar relativistic formulation developed in this work lays the theoretical foundation for the development of high-accuracy, low-cost correlated variational relativistic many-body theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Magnetic pair distribution function data using polarized neutrons and ad hoc corrections

Here, we report the first example of magnetic pair distribution function (mPDF) data obtained through the use of neutron polarization analysis. Using the antiferromagnetic semiconductor MnTe as a test case, we present high-quality mPDF data collected on the HYSPEC instrument at the Spallation Neutron Source using longitudinal polarization analysis to isolate the magnetic scattering cross section. Clean mPDF patterns are obtained for MnTe in both the magnetically ordered state and the correlated paramagnet state, where only short-range magnetic order is present. We also demonstrate significant improvement in the quality of high-resolution mPDF data through the application of ad hoc corrections that require only minimal human input, minimizing potential sources of error in the data processing procedure. We briefly discuss the current limitations and future outlook of mPDF analysis using polarized neutrons. Overall, this work provides a useful benchmark for mPDF analysis using polarized neutrons and provides an encouraging picture of the potential for routine collection of high-quality mPDF data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Identifying strengths and weaknesses of methods for computational network inference from single-cell RNA-seq data

Single-cell RNA-sequencing (scRNA-seq) offers unparalleled insight into the transcriptional programs of different cellular states by measuring the transcriptome of thousands of individual cells. An emerging problem in the analysis of scRNA-seq is the inference of transcriptional gene regulatory networks and a number of methods with different learning frameworks have been developed to address this problem. Here, we present an expanded benchmarking study of eleven recent network inference methods on seven published scRNA-seq datasets in human, mouse, and yeast considering different types of gold standard networks and evaluation metrics. We evaluate methods based on their computing requirements as well as on their ability to recover the network structure. We find that, while most methods have a modest recovery of experimentally derived interactions based on global metrics such as Area Under the Precision Recall curve, methods are able to capture targets of regulators that are relevant to the system under study. Among the top performing methods that use only expression were SCENIC, PIDC, MERLIN or Correlation. Addition of prior biological knowledge and the estimation of transcription factor activities resulted in the best overall performance with the Inferelator and MERLIN methods that use prior knowledge outperforming methods that use expression alone. We found that imputation for network inference did not improve network inference accuracy and could be detrimental. Comparisons of inferred networks for comparable bulk conditions showed that the networks inferred from scRNA-seq datasets are often better or at par with the networks inferred from bulk datasets. Our analysis should be beneficial in selecting methods for network inference. At the same time, this highlights the need for improved methods and better gold standards for regulatory network inference from scRNAseq datasets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Hybrid Basis and Multi-Center Grid Method for Strong-Field Processes

We present a time-dependent framework that combines a hybrid basis, consisting of Gaussian-type orbitals (GTOs) and finite-element discrete-variable representation (FEDVR) functions, with a multicenter grid to simulate strong-field and attosecond dynamics in atoms and molecules. The method incorporates the construction of the orthonormal hybrid basis, the evaluation of electronic integrals, a unitary time-propagation scheme, and the extraction of optical and photoelectron observables. Its accuracy and robustness are benchmarked on one-electron systems such as atomic hydrogen and the dihydrogen cation (H$^+_2$) through comparisons with essentially-exact reference results for bound-state energies, high-harmonic generation spectra, photoionization cross sections, and photoelectron momentum distributions. This work establishes the groundwork for its integration with quantum-chemistry methods, which is already operational but will be detailed in future work, thereby enabling ab initio simulations of correlated polyatomic systems in intense ultrafast laser fields.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1986

Benchmark full configuration-interaction calculations on H2O, F, and

Full configuration-interaction calculations are reported, and compared to other methods, for H2O at its equilibrium geometry and at two geometries with the H-O bonds stretched. Since the percentage of the self-consistent field (SCF) reference in the full configuration-interaction (FCI) wave function decreases greatly with the bond elongation, the accuracy of techniques based on a single reference do not compare well with the FCI results. However, the results from a complete active space SCF/multireference configuration-interaction (CASSCF/MRCI) treatment are in good agreement with the FCI. Correlation effects in F compared to Ne are far more similar than for Fcompared to Ne, despite F- and Ne being isoelectronic. Since the importance of higher than double excitations is more important for F- than F, a very high percentage of the correlation must be obtained to accurately compute the electron affinity. In a CASSCF/MRCI treatment the higher than quadruple excitations contribute 0.02 eV to the electron affinity (EA), even for modest basis sets.

Bauschlicher, C. W., Jr.↗
DOE OSTI2021

Chromium evaporation and oxidation characteristics of alumina-forming austenitic stainless steels for balance of plant applications in solid oxide fuel cells

The chromium (Cr) evaporation behavior of several different types of iron (Fe)-based AFA alloys and benchmark Cr 2 O 3 -forming Fe-based 310 and Ni-based 625 alloys was investigated for 500 h exposures at 800 °C to 900 °C in air with 10% H 2 O. The Cr evaporation rates from alumina-forming austenitic (AFA) alloys were ~5 to 35 times lower than that of the Cr 2 O 3 -forming alloys depending on alloy and temperature. The Cr evaporation behavior was correlated with extensive characterization of the chemistry and microstructure of the oxide scales, which also revealed a degree of quartz tube Si contamination during the test. As a result, long-term oxidation kinetics were also assessed at 800 to 1000 °C for up to 10,000 h in air with 10% H 2 O to provide further guidance for SOFC BOP component alloy selection.

08 HYDROGEN↗
DOE OSTI2020

Enabling Efficient and Accurate Computational Studies of MOF Reactivity via QM/MM and QM/QM Methods

Electronic structure calculations can provide unique insight into metal–organic framework (MOF) reactivity and defect formation. Such calculations can be broadly categorized as utilizing either periodic or cluster models, each with their respective advantages and disadvantages. In the present study, we demonstrate how multiscale methods can leverage the advantages of both approaches to enable high levels of accuracy and computational efficiency in studies of MOF reactivity. Using defect formation in a zeolitic imidazolate framework (ZIF) as a prototypical example, we benchmark a quantum mechanics/molecular mechanics (QM/MM) scheme that enables highly efficient cluster-based calculations on MOFs. We demonstrate the importance of correctly accounting for the influence of both dative bond cleavage in the QM cluster and long-range mechanical coupling to the bulk to achieving accurate QM/MM studies of MOFs. We subsequently leverage these cluster models in a QM/QM scheme that goes beyond standard DFT to yield “gold-standard” correlated wave function results on MOFs at modest computational cost. Crucially, we find several cases in which the incorporation of these correlated corrections yields qualitatively important corrections over conventional DFT values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1984

Benchmark data on the separability among crops in the southern San Joaquin Valley of California

Landsat MSS data were input to a discriminant analysis of 21 crops on each of eight dates in 1979 using a total of 4,142 fields in southern Fresno County, California. The 21 crops, which together account for over 70 percent of the agricultural acreage in the southern San Joaquin Valley, were analyzed to quantify the spectral separability, defined as omission error, between all pairs of crops. On each date the fields were segregated into six groups based on the mean value of the MSS7/MSS5 ratio, which is correlated with green biomass. Discriminant analysis was run on each group on each date. The resulting contingency tables offer information that can be profitably used in conjunction with crop calendars to pick the best dates for a classification. The tables show expected percent correct classification and error rates for all the crops. The patterns in the contingency tables show that the percent correct classification for crops generally increases with the amount of greenness in the fields being classified. However, there are exceptions to this general rule, notably grain.

Morse, A.↗
NASA NTRS1986

Benchmark full configuration-interaction calculations on H2O, F- and F

Full configuration-interaction calculations are reported, and compared to other methods, for H2O at its equilibrium geometry and at two geometries with the H-O bonds stretched. Since the percentage of the self-consistent field (SCF) reference in the full configuration-interaction (FCI) wave function decreases greatly with the bond elongation, the accuracy of techniques based on a single reference do not compare well with the FCI results. However, the results from a complete active space SCF/multireference configuration-interaction (CASSCF/MRCI) treatment are in good agreement with the FCI. Correlation effects in F compared to Ne are far more similar than for F- compared to Ne, despite F- and Ne being isoelectronic. Since the importance of higher than double excitations is more important for F- than F, a very high percentage of the correlation must be obtained to accurately compute the electron affinity. In a CASSCF/MRCI treatment the higher than quadruple excitations contribute 0.02 eV to the electron affinity (EA), even for modest basis sets.

Bauschlicher, C. W., Jr.↗
DOE OSTI2020

Hetero-site Double Core Ionization Energies with Sub-electronvolt Accuracy from Delta-Coupled-Cluster Calculations

Benchmark scalar-relativistic delta-coupled-cluster calculations of hetero-site double core ionization energies of small molecules containing second-row elements are reported. The present study has focused on the high-spin triplet components of two-site double core-ionized states, which are single reference in character and consistent with the use of standard coupledcluster methods. Contributions to computed double core ionization energies from electron-correlation and basis-set effects as well as corrections to the core-valence separation approximation have been analyzed. Furthermore, on the basis of systematic convergence of computational results with respect to these effects, delta-coupled-cluster calculations have been shown to be capable of providing accurate double core ionization energies with remaining errors estimated to be below 0.3 eV, and thus they are recommended for use to facilitate experimental studies of two-site double core-ionized states that are involved in X-ray pump/X-ray probe studies of electronic and molecular dynamics following inner shell ionization or excitation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Towards excitations and dynamical quantities in correlated lattices with density matrix embedding theory

Density matrix embedding theory (DMET) provides a framework to describe ground-state expectation values in strongly correlated systems, but its extension to dynamical quantities is still an open problem. We show one route to obtaining excitations and dynamical spectral functions by using the techniques of DMET to approximate the matrix elements that arise in a single-mode inspired excitation ansatz. We demonstrate this approach in the one-dimensional Hubbard model, comparing the neutral excitations, single-particle density of states, charge, and spin dynamical structure factors to benchmarks from the Bethe ansatz and density matrix renormalization group. Finally, our work highlights the potential of these ideas in building computationally efficient approaches for dynamical quantities.

1-dimensional systems↗