Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional modeling”

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

Structure Sensitive Reaction Kinetics of Chiral Molecules on Intrinsically Chiral Surfaces

Enantiospecific heterogeneous catalysis utilizes chiral surfaces to resolve enantiomers via structure sensitive surface chemistry. The catalyst design challenge is the identification of chiral surface structures that maximize enantiospecificity. Herein, we develop data driven models for the enantiospecificity of tartaric acid reactions on chiral Cu(hkl) R&S surfaces. Measurements of enantiospecific rate constants were obtained by using curved Cu(hkl) R&S surfaces that enable kinetic measurements on hundreds of chiral surface orientations. One model uses feature vectors derived from generalized coordination numbers to capture the local structure around Cu atoms exposed by the Cu(hkl) R&S surfaces. The second model introduces the use of chiral cubic harmonic functions to capture the symmetry constraints of the face-centered cubic Cu structure. The model using 58 generalized coordination numbers has a fitting error similar to that of the model using only 5 cubic harmonic functions. The two models predict maxima in the enantiospecificity on surfaces with very similar surface orientations. The models developed in this work are applicable for any enantiospecific reaction happening on any chiral material with a cubic lattice structure, opening the way to understanding the surface structure sensitivity of the enantiospecific reaction kinetics.

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

Fusion of Multiple Models for Improving Gross Primary Production Estimation With Eddy Covariance Data Based on Machine Learning

Abstract Terrestrial gross primary production (GPP) represents the magnitude of CO 2 uptake through vegetation photosynthesis, and is a key variable for carbon cycles between the biosphere and atmosphere. Light use efficiency (LUE) models have been widely used to estimate GPP for its physiological mechanisms and availability of data acquisition and implementation, yet each individual GPP model has exhibited large uncertainties due to input errors and model structure, and further studies of systematic validation, comparison, and fusion of those models with eddy covariance (EC) site data across diverse ecosystem types are still needed in order to further improve GPP estimation. We here compared and fused five GPP models (VPM, EC‐LUE, GOL‐PEM, CHJ, and C‐Fix) across eight ecosystems based on FLUXNET2015 data set using the ensemble methods of Bayesian Model Averaging (BMA), Support Vector Machine (SVM), and Random Forest (RF) separately. Our results showed that for individual models, EC‐LUE gave a better performance to capture interannual variability of GPP than other models, followed by VPM and GLO‐PEM, while CHJ and C‐Fix were more limited in their estimation performance. We found RF and SVM were superior to BMA on merging individual models at various plant functional types (PFTs) and at the scale of individual sites. On the basis of individual models, the fusion methods of BMA, SVM, and RF were examined by a five‐fold cross validation for each ecosystem type, and each method successfully improved the average accuracy of estimation by 8%, 18%, and 19%, respectively.

Environmental Sciences & Ecology↗
DOE OSTI2020

Probabilistic Prediction of Geomagnetic Storms and the K p index

Geomagnetic activity is often described using summary indices to summarize the likelihood of space weather impacts, as well as when parameterizing space weather models. The geomagnetic index K p in particular, is widely used for these purposes. Current state-of-the-art forecast models provide deterministic K p predictions using a variety of methods – including empirically-derived functions, physics-based models, and neural networks – but do not provide uncertainty estimates associated with the forecast. This paper provides a sample methodology to generate a 3-hour-ahead K p prediction with uncertainty bounds and from this provide a probabilistic geomagnetic storm forecast. Specifically, we have used a two-layered architecture to separately predict storm (K p ≥ 5 – ) and non-storm cases. As solar wind-driven models are limited in their ability to predict the onset of transient-driven activity we also introduce a model variant using solar X-ray flux to assess whether simple models including proxies for solar activity can improve the predictions of geomagnetic storm activity with lead times longer than the L1-to-Earth propagation time. By comparing the performance of these models we show that including operationally-available information about solar irradiance enhances the ability of predictive models to capture the onset of geomagnetic storms and that this can be achieved while also enabling probabilistic forecasts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Atomistically resolved hot exciton relaxation dynamics in CdSe quantum dots: Experiment and theory

Semiconductor quantum dots (QDs) are well known to give rise to a quantum confined structure of excitons. Because of this quantum confinement, new physics of hot exciton relaxation dynamics arises. Decades of work using transient absorption (TA) spectroscopy have yielded initial simple observations, such as estimates of the cooling rate from single pump photon energy experiments. More detailed TA experiments employed variable pump photon energies to measure excitonic state-resolved transition rates. These TA measurements, usually the simplest form, have been employed to characterize QDs and their relaxation dynamics to this day. Yet, these TA measurements are fundamentally lacking in their ability to measure energy-resolved hot exciton cooling, which requires observation of the full cooling history through the real excitonic manifold. Here, we employ coherent multi-dimensional spectroscopy (CMDS) to perform an atomistically directed study of hot exciton cooling in CdSe QDs, revealing energy resolved relaxation dynamics. CMDS experiments are compared with simulations and prior TA measurements and simpler theories. Our findings reveal a hot exciton relaxation dynamics landscape. This relaxation dynamics landscape is a linear or sub-linear function of excess energy for different structures of QDs, with a strong size dependence. Our model simulations parameterized by the empirical pseudopotential model reproduces the experimental functional form and the dependence upon QD diameter and shell.

Atomistic simulations↗
DOE OSTI2020

Accelerated Cu2O Reduction by Single Pt Atoms at the Metal-Oxide Interface§

The reducibility of metal oxides, when serving as the catalyst support or the active sites themselves, plays an important role in heterogeneous catalytic reactions. Here we present an integrated experimental and theoretical study that reveals how the addition of small amounts of atomically dispersed Pt at the metal/oxide interface dramatically enhances the reducibility of a Cu2O thin film by H2. X-ray photoelectron spectroscopy (XPS) and temperature-programmed desorption (TPD) results reveal that upon oxidation, a PtCu single-atom alloy (SAA) surface is covered by a thin Cu2O film and is therefore unable to dissociate H2. Despite this, in situ studies using ambient pressure (AP) XPS reveal that the presence of a small amount of Pt under the oxide layer can, at the single-atom limit, promote the reduction of Cu2O by H2 at room temperature. We build two density functional theory based surface models to better understand these experimental findings: A Cu2O/Cu(111)-like surface oxide layer, known as the ‘29’ oxide, in which Pt is alloyed into the Cu(111) surface as well as a PtCu SAA. Our calculations suggest that the increased activity is due to the presence of atomically dispersed Pt under the surface oxide layer, which weakens the Cu-O bonds in its immediate vicinity, thus making the interface between subsurface Pt and the surface oxide a nucleation site for the formation of metallic Cu. This initial step in the reduction process results in the presence of surface Pt atoms surrounded by metallic Cu patches, and the Pt atoms become active in H2 dissociation, which consequently accelerates the reduction of the oxide layer. This work demonstrates how isolated Pt atoms at the metal/oxide interface of a Cu-based catalyst accelerate the reduction of the oxide and, therefore, helps maintain the active, reduced state of the catalyst under reaction conditions, as is required for the synthesis of methanol from CO2 hydrogenation.

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

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

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

Profiling the BLAST bioinformatics application for load balancing on high-performance computing clusters

Abstract Background The Basic Local Alignment Search Tool (BLAST) is a suite of commonly used algorithms for identifying matches between biological sequences. The user supplies a database file and query file of sequences for BLAST to find identical sequences between the two. The typical millions of database and query sequences make BLAST computationally challenging but also well suited for parallelization on high-performance computing clusters. The efficacy of parallelization depends on the data partitioning, where the optimal data partitioning relies on an accurate performance model. In previous studies, a BLAST job was sped up by 27 times by partitioning the database and query among thousands of processor nodes. However, the optimality of the partitioning method was not studied. Unlike BLAST performance models proposed in the literature that usually have problem size and hardware configuration as the only variables, the execution time of a BLAST job is a function of database size, query size, and hardware capability. In this work, the nucleotide BLAST application BLASTN was profiled using three methods: shell-level profiling with the Unix “time” command, code-level profiling with the built-in “profiler” module, and system-level profiling with the Unix “gprof” program. The runtimes were measured for six node types, using six different database files and 15 query files, on a heterogeneous HPC cluster with 500+ nodes. The empirical measurement data were fitted with quadratic functions to develop performance models that were used to guide the data parallelization for BLASTN jobs. Results Profiling results showed that BLASTN contains more than 34,500 different functions, but a single function, RunMTBySplitDB, takes 99.12% of the total runtime. Among its 53 child functions, five core functions were identified to make up 92.12% of the overall BLASTN runtime. Based on the performance models, static load balancing algorithms can be applied to the BLASTN input data to minimize the runtime of the longest job on an HPC cluster. Four test cases being run on homogeneous and heterogeneous clusters were tested. Experiment results showed that the runtime can be reduced by 81% on a homogeneous cluster and by 20% on a heterogeneous cluster by re-distributing the workload. Discussion Optimal data partitioning can improve BLASTN’s overall runtime 5.4-fold in comparison with dividing the database and query into the same number of fragments. The proposed methodology can be used in the other applications in the BLAST+ suite or any other application as long as source code is available.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

MoorPy (Quasi-Static Mooring Analysis in Python)

MoorPy is a quasi-static mooring model and a suite of associated functions for mooring system analysis. The core model supports quasi-static analysis of moored floating systems including any arrangement of mooring lines and floating platforms. It solves the distributed position and tension of each mooring line segment using standard catenary equations. Floating platforms can be represented with linear hydrostatic characteristics. MoorPy automatically computes a floating system's equilibrium state and can be queried to identify a mooring system's nonlinear force-displacement relationships. Linearized stiffness matrices are efficiently computed using semi-analytic Jacobians. MoorPy also includes plotting functions and a library of mooring component property and cost coefficients. MoorPy can be used directly from Python scripts to perform mooring design and analysis tasks, or it can be coupled with other tools to compute quasi-static mooring reactions as part of a larger simulation.

Hall, Mathew↗
DOE OSTI2024

Modeling Multi-Step Organic Reactions: Can Density Functional Theory Deliver Misleading Chemistry?

Many organic reactions are characterized by a complex mechanism with a variety of transition states and intermediates of different chemical natures. Their correct and accurate theoretical characterization critically depends on the accuracy of the computational method used. In this work, we study a complex ambimodal cycloaddition with five transition states, two intermediates, and three products, and we ask whether density functional theory (DFT) can provide a correct description of this type of complex and multifaceted reaction. Our work fills a gap in that most systematic benchmarks of DFT for chemical reactions have considered much simpler reactions. Our results show that many density functionals not only lead to seriously large errors but also differ from one another in predicting whether the reaction is ambimodal. Only a few of the available functionals provide a balanced description of the complex and multifaceted reactions. The parameters varied in the tested functionals are the ingredients, the treatment of medium-range and nonlocal correlation energy, and the inclusion of Hartree–Fock exchange. Furthermore, these results show a clear need for more benchmarks on the mechanisms of large molecules in complex reactions.

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

Insights into Methodologies and Operational Details of Resource Adequacy Assessment: A Case Study with Application to a Broader Flexibility Framework

Assessing and maintaining resource adequacy (RA) is a core pillar of power systems. However, recent changes in the physical makeup of these systems and the conditions under which these systems must operate have yielded a renewed interest in the methods, metrics, and assumptions that underpin RA assessments. In this paper, we systematically explore a wide range of RA modeling dimensions, including: the objective function and level of operational detail in the underlying model formulation; the quantity (look-ahead) and quality (accuracy) of data that is available for making operational decisions within those models; and the physical configuration of solar photovoltaics (PV) with battery storage hybrid resources. We apply a set of probabilistic RA tools and production cost modeling tools to a realistic test system based loosely on a future Electric Reliability Council of Texas power system dominated by solar PV resources. Under the assumptions of our system and models, we find that multi-stage probabilistic assessments may provide a more robust evaluation of RA by capturing a wider range of operational and system interactions, but this comes at a computational cost of 1-2 orders of magnitude longer run time depending on the specific configuration. In addition, the information on thermal generator availability impacts RA performance by an order of magnitude more than solar resource forecasts, which is driven by the comparatively larger magnitude of thermal outages than solar forecast errors within our test system. Lastly, the flexibility provided by hybrid and other resources can help reduce system load-shedding event frequencies and enable the system to be more robust to inaccurate forecast information, and alternative hybrid inverter sizes can impact RA levels by 1-2 orders of magnitude. Our results point to the importance of a broader flexibility framework to describe the interaction between (1) flexibility "supply" from both physical resource capabilities and operational constraints considered in the modeling, and (2) flexibility "demand" from forecast errors, thermal generator outages, and other sources of uncertainty, as well as their RA impacts. Results are likely sensitive to the system buildout explored; future work could consider additional system configurations and conditions.

ENERGY PLANNING, POLICY, AND ECONOMY,POWER TRANSMI↗
DOE OSTI2020

One-dimensional, transient modeling of a fixed-bed regenerator as a replacement for recuperators in supercritical CO 2 power cycles

Fixed-bed regenerators may be a low cost alternative to recuperators for supercritical CO 2 power cycles. However, the transient nature of regenerators operating in sCO 2 power cycles has not been previously characterized. This study details a one-dimensional, transient model of a fixed-bed regenerator that captures phenomena such as local fluid properties and heat transfer coefficients, entrained fluid heat capacity, and the effects of conservation of mass during process switching. The model is validated with experimental results of a 10 kWth regenerator. This comparison shows that the model predicts temperature as a function of both time and space exceptionally well. The model predicts effectiveness, heat transfer, and discharge temperatures within 5%, and predicts pressure drop and fluid carryover within 10%. This model could be used in future studies to simulate regenerators within a supercritical CO 2 power cycle.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2023

Development of local-power-free, remote α -particle detection using optical fibers

We demonstrate the application of fluorescence optical fiber coupled to a telecom grade fiber as a sensor for alpha particles using alpha-specific ZnS(Ag) scintillation materials whose wavelength is down-shifted into a low-loss region of the telecom grade fiber transmission band. Telecom-grade fiber optics offer a solution for sensing alpha radiation in deep repositories and cask storage for radioactive materials due to the stability of SiO 2 under normal environmental conditions and its relative radiation hardness at low radiation doses. Long-term nuclear waste storage facilities require sensors for the detection of leakage of radioactive materials that are maintenance-free, do not require power and can survive with no ‘wear out’ mechanisms for decades. By accomplishing the wavelength transformation, we maximize efficiencies in the detection of α-particles and signal transport and can detect alpha scintillation at distances on the order of >1 km with a sensor that is ~3% efficient and can be easily scaled as a sensor array. This paper describes the construction and testing of the sensor including manufacture of the controlled thickness films, verification of the wavelength shift from 450 to 620 nm and optimization of the sensitivity as a function of thickness. We also model the relative sensitivity of the film as a function of film thickness, and we demonstrate a signal-to-noise ratio of 10 at a range of greater than 1 km.

alpha particles↗
DOE OSTI2023

Desmearing small-angle scattering data by central moment expansions

This article proposes a robust method for desmearing 1D small-angle scattering (SAS) intensity profiles. Using the central-moment-expansion technique, smearing in SAS intensity profiles is achieved in a model-free manner, eliminating the need for presumptive input during the desmearing. This innovative approach enables data analysis of SAS data without the resolution being convoluted with the structural model, thereby facilitating faster implementation and performance of fitting functions directly from the scattering models.

42 ENGINEERING↗
DOE OSTI2021

Sodium Fire Collaborative Study Progress CNWG Fiscal Year 2021

This report discusses the progress on the collaboration between Sandia National Laboratories (Sandia) and Japan Atomic Energy Agency (JAEA) on the sodium fire research in fiscal year (FY) 2021 and is a continuation of the FY2020 progress report. In this report, we only report the changes made to the current sodium pool fire model in MELCOR. We modified and corrected many control functions to enhance the current model from the FY2020 report. This year’s enhancements relate to better agreement of the suspended aerosol measurement from JAEA’s F7 series tests. Staff from Sandia and JAEA conducted the validation studies of the sodium pool fire model in MELCOR. To validate this pool fire model with the latest enhancement, JAEA sodium pool fire experiments (F7-1 and F7-2) were used. The results of the calculation, including the code-to-code comparisons are discussed as well as suggestions for further model improvement. Finally, recommendations are made for new MELCOR simulations for FY2022.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A self-adjusting platinum surface for acetone hydrogenation

We show that platinum displays a self-adjusting surface that is active for the hydrogenation of acetone over a wide range of reaction conditions. Reaction kinetics measurements under steady-state and transient conditions at temperatures near 350 K, electronic structure calculations employing density-functional theory, and microkinetic modeling were employed to study this behavior over supported platinum catalysts. The importance of surface coverage effects was highlighted by evaluating the transient response of isopropanol formation following either removal of the reactant ketone from the feed, or its substitution with a similarly structured species. The extent to which adsorbed intermediates that lead to the formation of isopropanol were removed from the catalytic surface was observed to be higher following ketone substitution in comparison to its removal, indicating that surface species leading to isopropanol become more strongly adsorbed on the surface as the coverage decreases during the desorption experiment. This phenomenon occurs as a result of adsorbate–adsorbate repulsive interactions on the catalyst surface which adjust with respect to the reaction conditions. Reaction kinetics parameters obtained experimentally were in agreement with those predicted by microkinetic modeling when the binding energies, activation energies, and entropies of adsorbed species and transition states were expressed as a function of surface coverage of the most abundant surface intermediate (MASI, C 3 H 6 OH*). It is important that these effects of surface coverage be incorporated dynamically in the microkinetic model (e.g., using the Bragg–Williams approximation) to describe the experimental data over a wide range of acetone partial pressures.

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

Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics

Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/$\mathrm{cm^3}$ to 16 g/$\mathrm{cm^3}$ and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method. Here, we find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the sesame database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT-MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Superposed periodic kink and pulse solutions of coupled nonlinear equations

Here, we present novel previously unexplored periodic solutions, expressed in terms of Jacobi elliptic functions, for both a coupled Φ 4 model and a coupled nonlinear Schrödinger equation (NLS) model. Remarkably, these solutions can be elegantly reformulated as a linear combination of periodic kinks and antikinks, or as a combination of two periodic kinks or two periodic pulse solutions. However, we also find that for $m=0$ and a specific value of the periodicity (or at a nonzero value of the elliptic modulus $m$) this superposition does not hold. These results demonstrate that the notion of superposed solutions extends to the coupled nonlinear equations as well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗