Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density functional”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Integrated Experimental and Computational K-Edge X-ray Absorption Near-Edge Structure Analysis of Vanadium Catalysts

X-ray Absorption Near-Edge Structure (XANES) spectroscopy is a powerful tool to reveal key structural and electronic features of isolated catalytic sites, yet insights into molecular structure and more detailed orbital analysis through a combination of experimental and computed XANES analysis are necessary for accurate interpretation of the spectra, especially when significant heterogeneity exits among the catalytic sites. Herein, we present an integrated computational and experimental strategy to determine both primary and secondary bonding interactions within the XANES pre-edge region for organovanadium complexes, which was developed using a series of well-defined molecular vanadium complexes and then applied to the characterization of a supported organovanadium olefin hydrogenation catalyst. Time-dependent density functional theory is used to predict the energy of pre-edge XANES features for a series of vanadium complexes with a variety of oxidation states and local coordination environments. Further, a calibration scheme incorporating different density functionals and basis sets is established, resulting in an optimized scheme that accurately predicts pre-edge energies with a mean absolute error of 0.40 eV. Second-shell coordination (e.g., V---V) effects within XANES are identified through the analysis of the computed dominant orbital contributions for multi-vanadium complexes. Orbital analysis also provided confirmation that the vanadium-hydride formation combined with the heterogeneity of the catalytic active species in ole-fin hydrogenation caused the energy shift and broadening of the pre-edge peak after hydrogen treatment of the silica-supported organovanadium pre-catalyst. This work further elucidates computational XANES simulations and techniques potentially guiding characterization in surface organometallic chemistry.

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

Nonlinear Attitude Filtering Methods

This paper provides a survey of modern nonlinear filtering methods for attitude estimation. Early applications relied mostly on the extended Kalman filter for attitude estimation. Since these applications, several new approaches have been developed that have proven to be superior to the extended Kalman filter. Several of these approaches maintain the basic structure of the extended Kalman filter, but employ various modifications in order to provide better convergence or improve other performance characteristics. Examples of such approaches include: filter QUEST, extended QUEST, the super-iterated extended Kalman filter, the interlaced extended Kalman filter, and the second-order Kalman filter. Filters that propagate and update a discrete set of sigma points rather than using linearized equations for the mean and covariance are also reviewed. A two-step approach is discussed with a first-step state that linearizes the measurement model and an iterative second step to recover the desired attitude states. These approaches are all based on the Gaussian assumption that the probability density function is adequately specified by its mean and covariance. Other approaches that do not require this assumption are reviewed, including particle filters and a Bayesian filter based on a non-Gaussian, finite-parameter probability density function on SO(3). Finally, the predictive filter, nonlinear observers and adaptive approaches are shown. The strengths and weaknesses of the various approaches are discussed.

F Landis Markley↗
DOE OSTI2022

Evaluation of Thermal Neutron Scattering Law and Cross Sections for Calcium Hydride

Presented in this paper are the calculated thermal scattering law (TSL) and thermal neutron scattering cross sections for Calcium Hydride, hereafter referred to by its chemical symbol CaH 2 . The only other such data prior to this evaluation are thermal neutron scattering libraries in the JEFF database, which suffer from nonphysical features and inaccuracies. The data in this evaluation are calculated from first principles; Density Functional Theory (DFT) is used to calculate the phonon density of states (DOS), which is the primary input required to calculate the TSL. The TSL and cross sections have been evaluated for the three non-equivalent atom cites in the CaH 2: Ca, H 1 , and H 2 . Each evaluation has been submitted to the NNDC for consideration in the next ENDF/B release.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

The Vertical Distribution of Thin Features Over the Arctic Analysed from CALIPSO Observations: Optically Thin Clouds - Part 1

Clouds play a crucial role in the Arctic climate system. Therefore, it is essential to accurately and reliably quantify and understand cloud properties over the Arctic. It is also important to monitor and attribute changes in Arctic clouds. Here, we exploit the capability of the CALIPSO-CALIOP instrument and provide comprehensive statistics of tropospheric thin clouds, otherwise extremely difficult to monitor from passive satellite sensors.We use 4 yr of data (June 2006.May 2010) over the circumpolar Arctic, here defined as 67-82 deg. N, and characterize probability density functions of cloud base and top heights, geometrical thickness and zonal distribution of such cloud layers, separately for water and ice phases, and discuss seasonal variability of these properties. When computed for the entire study area, probability density functions of cloud base and top heights and geometrical thickness peak at 200-400, 1000-2000 and 400-800 m, respectively, for thin water clouds, while for ice clouds they peak at 6-8, 7-9 and 400-1000 m, respectively. In general, liquid clouds were often identified below 2 km during all seasons, whereas ice clouds were sensed throughout the majority of the upper troposphere and also, but to a smaller extent, below 2 km for all seasons.

DeVasthale, Abhay↗
NASA NTRS1991

Progress in the development of PDF turbulence models for combustion

A combined Monte Carlo-computational fluid dynamic (CFD) algorithm was developed recently at Lewis Research Center (LeRC) for turbulent reacting flows. In this algorithm, conventional CFD schemes are employed to obtain the velocity field and other velocity related turbulent quantities, and a Monte Carlo scheme is used to solve the evolution equation for the probability density function (pdf) of species mass fraction and temperature. In combustion computations, the predictions of chemical reaction rates (the source terms in the species conservation equation) are poor if conventional turbulence modles are used. The main difficulty lies in the fact that the reaction rate is highly nonlinear, and the use of averaged temperature produces excessively large errors. Moment closure models for the source terms have attained only limited success. The probability density function (pdf) method seems to be the only alternative at the present time that uses local instantaneous values of the temperature, density, etc., in predicting chemical reaction rates, and thus may be the only viable approach for more accurate turbulent combustion calculations. Assumed pdf's are useful in simple problems; however, for more general combustion problems, the solution of an evolution equation for the pdf is necessary.

Hsu, Andrew T.↗
DOE OSTI2023

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Chemistry↗
DOE OSTI2020

Interlayer magnetism in Fe 3-x GeTe 2

Fe 3 - x GeTe 2 is a layered van der Waals magnetic material with a relatively high ordering temperature and large anisotropy. While most studies have concluded the interlayer ordering to be ferromagnetic, there have also been reports of interlayer antiferromagnetism in Fe 3 - x GeTe 2 . In this work, we investigate the interlayer magnetic ordering by neutron diffraction experiments, scanning tunneling microscopy (STM) and spin-polarized STM measurements, density functional theory plus U calculations, and STM simulations. We conclude that the layers of Fe 3 - x GeTe 2 are coupled ferromagnetically and that in order to capture the magnetic and electronic properties of Fe 3 - x GeTe 2 within density functional theory, Hubbard U corrections need to be taken into account.

2-dimensional systems↗
NASA NTRS2003

Probabilistic Aeroelastic Analysis Developed for Turbomachinery Components

Aeroelastic analyses for advanced turbomachines are being developed for use at the NASA Glenn Research Center and industry. However, these analyses at present are used for turbomachinery design with uncertainties accounted for by using safety factors. This approach may lead to overly conservative designs, thereby reducing the potential of designing higher efficiency engines. An integration of the deterministic aeroelastic analysis methods with probabilistic analysis methods offers the potential to design efficient engines with fewer aeroelastic problems and to make a quantum leap toward designing safe reliable engines. In this research, probabilistic analysis is integrated with aeroelastic analysis: (1) to determine the parameters that most affect the aeroelastic characteristics (forced response and stability) of a turbomachine component such as a fan, compressor, or turbine and (2) to give the acceptable standard deviation on the design parameters for an aeroelastically stable system. The approach taken is to combine the aeroelastic analysis of the MISER (MIStuned Engine Response) code with the FPI (fast probability integration) code. The role of MISER is to provide the functional relationships that tie the structural and aerodynamic parameters (the primitive variables) to the forced response amplitudes and stability eigenvalues (the response properties). The role of FPI is to perform probabilistic analyses by utilizing the response properties generated by MISER. The results are a probability density function for the response properties. The probabilistic sensitivities of the response variables to uncertainty in primitive variables are obtained as a byproduct of the FPI technique. The combined analysis of aeroelastic and probabilistic analysis is applied to a 12-bladed cascade vibrating in bending and torsion. Out of the total 11 design parameters, 6 are considered as having probabilistic variation. The six parameters are space-to-chord ratio (SBYC), stagger angle (GAMA), elastic axis (ELAXS), Mach number (MACH), mass ratio (MASSR), and frequency ratio (WHWB). The cascade is considered to be in subsonic flow with Mach 0.7. The results of the probabilistic aeroelastic analysis are the probability density function of predicted aerodynamic damping and frequency for flutter and the response amplitudes for forced response.

Reddy, T. S. R.↗
DOE OSTI2020

Revisiting Competing Paths in Electrochemical CO 2 Reduction on Copper via Embedded Correlated Wavefunction Theory

We re-evaluate two key steps in the mechanism of CO 2 reduction on copper at a higher level of theory capable of correcting inherent errors in density functional theory (DFT) approximations, namely, embedded correlated wavefunction (ECW) theory. Here, we consider the CO reduction step on Cu(111), which is critical to understanding reaction selectivity. Here, we optimize embedding potentials at the periodic plane-wave DFT level using density functional embedding theory (DFET). All possible adsorption sites (adsites) for each adsorbate then are screened with ECW theory at the catalytically active site to refine the local electronic structure. Unsurprisingly, DFT and ECW theory predict different adsite preferences, largely because of DFT’s inability to properly situate the CO 2π* level. Differing preferred adsites suggest that different reaction pathways could emerge from DFT versus ECW theory. Starting from these preferred ECW theory adsites, we then obtain reaction pathways at the plane-wave DFT level using the climbing-image nudged elastic band method to determine minimum energy paths. Thereafter, we perform ECW calculations at the catalytically active site to correct the energetics at each interpolated structure (image) along the reaction pathways. Via this approach, we confirm that the first step in CO reduction via hydrogen transfer on Cu(111) is to form hydroxymethylidyne (*COH) instead of formyl (*CHO). Although the prediction to preferentially form *COH is consistent with that of DFT, the two theories predict quite different structural and mechanistic behaviors, suggesting that verification is needed for other parts of the mechanism of CO 2 reduction, which is the subject of ongoing work.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Examination of probability distribution of mixture fraction in LES/FDF modelling of a turbulent partially premixed jet flame

An accurate prediction of the probability density function (PDF) of the mixture fraction is crucial to the prediction of combustion since mixing plays an important role in turbulent non-premixed and partially premixed flames. This work provides an assessment of the large-eddy simulation (LES)/filtered density function (FDF) method for the prediction of the PDF of the mixture fraction. The advantage of the LES/FDF method is that it provides the full predictions of the statistical distribution of scalars including the mixture fraction. The predictive accuracy of the method for the PDF is yet to be fully validated. Assessing the prediction of the PDF of the mixture fraction, a conserved scalar, is an important starting point. The Sydney/Sandia inhomogeneous inlet jet flame is used as a test case. A quick comparison shows that the LES/FDF predicted PDF shapes of the mixture fraction deviate significantly from the commonly presumed Beta-PDF as well as from the experimental data in the flame. Here, to examine the source of the discrepancy, we clarify the different PDF definitions used in the comparison among the predictions, measurements, and the presumed shape PDFs. The discrepancy observed from the comparison is largely reconciled by clarifying the difference between the PDFs that are examined. The PDF of the resolved mixture fraction is shown to be close to the Beta-PDF in both the measurements and predictions, while the PDF directly deduced from the LES/FDF particles deviates significantly from the Beta-PDF. A multimodal PDF analysis and a pseudo convergence analysis are conducted to provide plausible evidence to support the predicted multimodal PDF shapes. The sub-filter scale FDF is shown to be close to the Beta-PDF too through the construction of a synthesized PDF, which supports the common presumed Beta-PDF assumption used in the presumed PDF methods when combined with LES.

42 ENGINEERING↗
DOE OSTI2020

Operando identification of site-dependent water oxidation activity on ruthenium dioxide single-crystal surfaces

Understanding the nature of active sites is central to controlling (electro)catalytic activity. Here we employed surface X-ray scattering coupled with density functional theory and surface-enhanced infrared absorption spectroscopy to examine the oxygen evolution reaction on RuO 2 surfaces as a function of voltage. At 1.5 V RHE , our results suggest that there is an –OO group on the coordinatively unsaturated ruthenium (Ru CUS ) site of the (100) surface (and similarly for (110)), but adsorbed oxygen on the Ru CUS site of (101). Density functional theory results indicate that the removal of –OO from the Ru CUS site, which is stabilized by a hydrogen bond to a neighbouring –OH (–OO–H), could be the rate-determining step for (100) (similarly for (110)), where its reduced binding on (100) increased activity. A further reduction in binding energy on the Ru CUS site of (101) resulted in a different rate-determining step (–O + H 2 O – (H + + e - ) → –OO–H) and decreased activity. Our study provides molecular details on the active sites, and the influence of their local coordination environment on activity.

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

A numerical method for integrating the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes

An extension of the method of moments is developed for the numerical integration of the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes. The number density function n sub k (x,t) in each separate droplet packet between droplet mass grid points (x sub k, x sub k+1) is represented by an expansion in orthogonal polynomials with a given weighting function. In this way droplet number concentrations, liquid water contents and other moments in each droplet packet are conserved and the problem of solving the kinetic equations is replaced by one of solving a set of coupled differential equations for the number density function moments. The method is tested against analytic solutions of the corresponding kinetic equations. Numerical results are obtained for different coalescence/breakup and condensation/evaporation kernels and for different initial droplet spectra. Also droplet mass grid intervals, weighting functions, and time steps are varied.

Emukashvily, I. M.↗
DOE OSTI2020

GPU acceleration of all-electron electronic structure theory using localized numeric atom-centered basis functions

We present an implementation of all-electron density-functional theory for massively parallel GPU-based platforms, using localized atom-centered basis functions and real-space integration grids. Special attention is paid to domain decomposition of the problem on non-uniform grids, which enables compute- and memory-parallel execution across thousands of nodes for real-space operations, e.g. the update of the electron density, the integration of the real-space Hamiltonian matrix, and calculation of Pulay forces. To assess the performance of our GPU implementation, we performed benchmarks on three different architectures using a 103-material test set. We find that operations which rely on dense serial linear algebra show dramatic speedups from GPU acceleration: in particular, SCF iterations including force and stress calculations exhibit speedups ranging from 4.5 to 6.6. For the architectures and problem types investigated here, this translates to an expected overall speedup between 3–4 for the entire calculation (including non-GPU accelerated parts), for problems featuring several tens to hundreds of atoms. Additional calculations for a 375-atom Bi2Se3 bilayer show that the present GPU strategy scales for large-scale distributed-parallel simulations.

42 ENGINEERING↗
DOE OSTI2023

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

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

From Oxo to Oxyl to Biradical: Systematic Multireference Calculations of Methane Activation at MOF Nodes

Methane C–H activation at transition-metal sites often involves electronic structures that challenge conventional single-reference electronic structure descriptions. Although Kohn–Sham density functional theory (DFT) is widely used to study catalytic trends, its reliability for reactions involving strongly correlated species remains uncertain. Here we present a systematic multireference investigation of methane activation at metal–organic framework (MOF) node catalysts across the 3d transition-metal series. We introduce an automated workflow for active space selection to enable consistent application of multireference methods, including multiconfiguration pair-density functional theory and n-electron valence state perturbation theory, to these catalytic systems. These calculations show substantial static correlation in the C–H activation reaction step and predict activation barriers that differ from DFT by 30–70 kJ mol–1, with DFT often qualitatively disagreeing in barrier height trends across transition metals. Analysis of multireference wave functions shows that reactivity is governed by the electronic structure of the M–O moiety along a continuum from metal–oxo to oxyl radical and O biradical character. Increased oxygen-centered spin density and weakened M–O bonding are identified as descriptors of catalytic activity which correlate with lower activation barriers.

Wardzala, Jacob↗
DOE OSTI2023

Mechanical and Thermal Characterization of Phase-Change Material and High-Density Polyethylene Functional Composites for Thermal Energy Storage

Abstract Phase-change materials (PCMs) can be used to develop thermal energy storage systems as they absorb large amount of latent heat nearly at a constant temperature when changing phase from a solid to a liquid. To prevent leakage when in a liquid state, PCM is shape stabilized in a polymer matrix of high-density polyethylene (HDPE). The present research explores the injection-molded mechanical and thermal properties of different PCM/HDPE composite ratios. The tensile strength and modulus of elasticity at room temperature and with the PCM fully melted within the composite are measured. Additionally, the hardness, latent heat of fusion, phase-change temperature, and thermal conductivity are investigated. An analysis of microstructures of the composite is used to support the findings. The PCM within the PCM/HDPE composite gives it the benefit of thermal storage but causes a decrease in mechanical properties.

Energy & Fuels↗
DOE OSTI2021

From Molecules to Solids: A vdW-DF-C09 Case Study of the Mercury Dihalides

The mercury dihalides show a remarkable diversity in the structural preferences in their minimum energy structure types, spanning molecular to strongly bound ionic solids. A challenge in the development of density functional methods for extended systems is to arrive at strategies that serve equally well such a broad range of bonding modes or structural preferences. The chemical bonding and the stabilities of mercury dihalides and the general utility and reliability of the van der Waals density functional with C09 exchange (vdW-DF-C09) in predicting or describing the energetics and structural preferences in these metal dihalides is examined. Here, we show that, in contrast with the uncorrected generalized gradient approximation of the Perdew-Burke-Erzenhoff (PBE) exchange-correlation functional, qualitative and quantitative patterns in the bonding of the mercury dihalide solids are well reproduced with vdW-DF-C09 for the full series of HgX 2 systems for X = F, Cl, Br, and I. The possible existence of a low-temperature cotunnite polymorph for HgF 2 and PbF 2 is posited.

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

Unveiling the Role of Surface Ir-Oxo Species in O 2 Evolution at IrO 2 Electrocatalysts via Embedded Cluster Multireference Calculations

Understanding the mechanisms driving the oxygen evolution reaction (OER) on iridium oxide (IrO 2 )-based catalysts is essential to improving their performance and enabling an actual scale-up of water-splitting photoelectrochemical cells. The mechanistic pathways at IrO 2 interfaces have been extensively investigated computationally using density functional theory (DFT), which predicts a high-energy barrier for the last step of the OER of molecular oxygen detachment and release from the catalyst surface. Nevertheless, surface O 2 over- and under-binding results by standard generalized gradient approximation and hybrid density functionals, respectively, call for further analysis of this crucial step via multireference methods. Aiming at unveiling the nature of such a barrier, we hereby address the formation of O 2 from the most-stable IrO 2 (110) surface with both periodic DFT and an electrostatic embedded cluster approach at the n-electron valence-state perturbation theory. With this multireference approach, we find a value for the aforementioned energy barrier that is much closer to experimental indications than DFT ones. Here, an in-depth analysis of the involved molecular orbitals suggests that the origin of this barrier is related to the breaking of a π interaction between O 2 and Ir surface atom and to a significant additional O 2 interaction with adjacent electrophilic Ir-oxo species, which is present under experimental operating conditions. Besides shedding light on the mechanism of the OER on IrO 2 , these findings point out the importance of multireference methods for dissecting complex reactions at electrocatalytic interfaces and pave the route for further investigations with effective embedding approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗