Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional approximation”

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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Derivative discontinuity and exchange-correlation potential of meta-GGAs in density-functional theory

We investigate fundamental properties of meta-generalized-gradient approximations (meta-GGAs) to the exchange-correlation energy functional, which have an implicit density dependence via the Kohn-Sham kinetic-energy density. To this purpose, we construct the most simple meta-GGA by expressing the local exchange-correlation energy per particle as a function of a fictitious density, which is obtained by inverting the Thomas-Fermi kinetic-energy functional. This simple functional considerably improves the total energy of atoms as compared to the standard local density approximation. The corresponding exchange-correlation potentials are then determined exactly through a solution of the optimized effective potential equation. These potentials support an additional bound state and exhibit a derivative discontinuity at integer particle numbers. We further demonstrate that through the kinetic-energy density any meta-GGA incorporates a derivative discontinuity. However, we also find that for commonly used meta-GGAs the discontinuity is largely underestimated and in some cases even negative.

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

Simple Models of Lateral Heterostructures

General analytical expressions for densities of states of a lateral heterostructure, formed by the contact of two square semi-infinite lattices with single-band and two-band spectra, were obtained in the tight-binding approximation by the Green’s function method. The semi-elliptical density of states was used for numerical estimates, and the model of two interacting dimers was proposed to estimate the charge transfer. Application of this approach to description of lateral epitaxial and graphene-like heterostructures is discussed.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Approximating and incorporating model uncertainty in an inversion for seismic source functions: Preliminary results

We present preliminary work on propagating model uncertainty into the estimation of the time domain source time functions of the seismic source. Our method is based on an estimated model covariance function, which we estimate from the data. The model covariance function is then used to construct a suite of surrogate Greens functions which we use in a Monte Carlo type inversion scheme. The result is a probability density function of the six independent source time functions, each of which corresponds to an individual component of the seismic moment tensor. We compare the results of our method with those obtained using a computationally expensive finite difference Monte Carlo method and find that our new method produces results that are deficient in low frequencies. The advantage of our new method, which we term the Karhunen-Loeve Monte Carlo (KLMC) method, is that is several orders of magnitude faster than our current method, which uses a finite difference scheme to produce the suite of forward models.

42 ENGINEERING↗
DOE OSTI2020

A Multi-Function AAA Algorithm Applied to Frequency-Dependent Line Modeling: Preprint

We present a new multi-input rational approximation algorithm (MIAAA) and illustrate its advantages respect to the popular Vector Fitting (VF) algorithm using examples of approximations of admittance matrices discussed in a very recent paper and by the creators of VF. We show that not only MIAAA outperforms VF in terms of achieving better accuracy using a smaller number of poles, but also has no numerical issues to achieve convergence. In contrast to VF, MIAA is not sensitive to the location of input sample points and it does not require good estimates for the location of the desired approximation poles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Assessing density functionals for describing methane dissociative chemisorption on Pt(110)-(2×1) surface

In this work, we explore the suitability of several density functionals with the generalized gradient approximation (GGA) and beyond for describing the dissociative chemisorption of methane on the reconstructed Pt(110)-(2×1) surface. The bulk and surface structures of the metal, methane adsorption energy, and dissociation barrier are used to assess the functionals. A van der Waals corrected GGA functional (optPBE-vdW) and a metaGGA functional with van der Waals correction (MS PBEl-rVV10) are selected for ab initio molecular dynamics calculations of the sticking probability. Our results suggest that the use of these two functionals may lead to a better agreement with existing experimental results, thus serving as a good starting point for future development of reliable machine-learned potential energy surfaces for the dissociation of methane on the Pt(110)-(2×1) surface.

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

The Problem of Torsion in Prismatic Members of Circular Segmental Cross Section

The problem is solved by approximation, by setting up a function complying with the differential equation of the stress function, and determining the coefficients appearing in it in such a way that the boundary condition is fulfilled as nearly as possible. For the semicircle, for which the solution is known, the method yields very accurate values; the approximated stress distribution is in good agreement with the accurately computed distribution. Stress and strain measurements indicate that the approximate solution is in sufficiently exact agreement with reality for segmental cross sections.

Weigand, A.↗
DOE OSTI2016

Methane dissociation on Pt(111): Searching for a specific reaction parameter density functional

The theoretical description of methane dissociating on metal surfaces is a current frontier in the field of gas-surface dynamics. Dynamical models that aim at achieving a highly accurate description of this reaction rely on potential energy surfaces based on density functional theory calculations at the generalized gradient approximation. We focus here on the effect that the exchange-correlation functional has on the reactivity of methane on a metal surface, using CHD{sub 3} + Pt(111) as a test case. We present new ab initio molecular dynamics calculations performed with various density functionals, looking also at functionals that account for the van der Waals (vdW) interaction. While searching for a semi-empirical specific reaction parameter density functional for this system, we find that the use of a weighted average of the PBE and the RPBE exchange functionals together with a vdW-corrected correlation functional leads to an improved agreement with quantum state-resolved experimental data for the sticking probability, compared to previous PBE calculations. With this semi-empirical density functional, we have also investigated the surface temperature dependence of the methane dissociation reaction and the influence of the rotational alignment on the reactivity, and compared our results with experiments.

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

Results of APL rain gauge network measurements in mid-Atlantic coast region and comparisons of distributions with CCIR models

In this effort are described cumulative rain rate distributions for a network of nine tipping bucket rain gauge systems located in the mid-Atlantic coast region in the vicinity of the NASA Wallops Flight Facility, Wallops Island, Virginia. The rain gauges are situated within a gridded region of dimensions of 47 km east-west by 70 km north-south. Distributions are presented for the individual site measurements and the network average for the year period June 1, 1986 through May 31, 1987. A previous six year average distribution derived from measurements at one of the site locations is also presented. Comparisons are given of the network average, the CCIR (International Radio Consultative Committee) climatic zone, and the CCIR functional model distributions, the latter of which approximates a log normal at the lower rain rate and a gamma function at the higher rates.

Goldhirsh, Julius↗
DOE OSTI2019

Status of free-energy representations for the homogeneous electron gas

Renewed interest in the homogeneous electron gas (HEG) has been stimulated by recent accurate simulations of it over a wide range of densities and temperatures. Those data, combined with known theoretical limits, have led to analytical representations of the free energy. Such a representation is, at least in principle, the complete HEG equation of state. The initial objective here is to establish that the two best representations of the simulation data and constraints are effectively the same in both functional form and accuracy of representation. The second objective is to disclose and delineate a significant difficulty. Despite their expected accuracy for the free energy, the underlying functional form is not adequate for derived thermodynamic properties. As an example, the specific heats obtained from the representations exhibit anomalies suggesting the need first for additional simulation data in critical regimes, then for refined fitting functions. Here, the existing representations are, however, almost certainly adequate for applications based on the free energy alone (e.g., density functional theory for warm dense matter). The third objective is to show that, despite their inability to provide a complete thermodynamic description of the HEG, the best analytical representations do provide a fully adequate exchange-correlation local density approximation for free energy density functional calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1982

Difference equation state approximations for nonlinear hereditary control problems

Discrete approximation schemes for the solution of nonlinear hereditary control problems are constructed. The methods involve approximation by a sequence of optimal control problems in which the original infinite dimensional state equation has been approximated by a finite dimensional discrete difference equation. Convergence of the state approximations is argued using linear semigroup theory and is then used to demonstrate that solutions to the approximating optimal control problems in some sense approximate solutions to the original control problem. Two schemes, one based upon piecewise constant approximation, and the other involving spline functions are discussed. Numerical results are presented, analyzed and used to compare the schemes to other available approximation methods for the solution of hereditary control problems.

Rosen, I. G.↗
NASA NTRS1984

Difference equation state approximations for nonlinear hereditary control problems

Discrete approximation schemes for the solution of nonlinear hereditary control problems are constructed. The methods involve approximation by a sequence of optimal control problems in which the original infinite dimensional state equation has been approximated by a finite dimensional discrete difference equation. Convergence of the state approximations is argued using linear semigroup theory and is then used to demonstrate that solutions to the approximating optimal control problems in some sense approximate solutions to the original control problem. Two schemes, one based upon piecewise constant approximation, and the other involving spline functions are discussed. Numerical results are presented, analyzed and used to compare the schemes to other available approximation methods for the solution of hereditary control problems. Previously announced in STAR as N83-33589

Rosen, I. G.↗
DOE OSTI2014

Low-rank separated representation surrogates of high-dimensional stochastic functions: Application in Bayesian inference

This study introduces a non-intrusive approach in the context of low-rank separated representation to construct a surrogate of high-dimensional stochastic functions, e.g., PDEs/ODEs, in order to decrease the computational cost of Markov Chain Monte Carlo simulations in Bayesian inference. The surrogate model is constructed via a regularized alternative least-square regression with Tikhonov regularization using a roughening matrix computing the gradient of the solution, in conjunction with a perturbation-based error indicator to detect optimal model complexities. The model approximates a vector of a continuous solution at discrete values of a physical variable. The required number of random realizations to achieve a successful approximation linearly depends on the function dimensionality. The computational cost of the model construction is quadratic in the number of random inputs, which potentially tackles the curse of dimensionality in high-dimensional stochastic functions. Furthermore, this vector-valued separated representation-based model, in comparison to the available scalar-valued case, leads to a significant reduction in the cost of approximation by an order of magnitude equal to the vector size. The performance of the method is studied through its application to three numerical examples including a 41-dimensional elliptic PDE and a 21-dimensional cavity flow.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2013

A Monte Carlo algorithm for degenerate plasmas

A procedure for performing Monte Carlo calculations of plasmas with an arbitrary level of degeneracy is outlined. It has possible applications in inertial confinement fusion and astrophysics. Degenerate particles are initialised according to the Fermi–Dirac distribution function, and scattering is via a Pauli blocked binary collision approximation. The algorithm is tested against degenerate electron–ion equilibration, and the degenerate resistivity transport coefficient from unmagnetised first order transport theory. The code is applied to the cold fuel shell and alpha particle equilibration problem of inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

First-principles investigation of uranium mononitride (UN): Effect of magnetic ordering, spin-orbit interactions and exchange correlation functional

Uranium mononitride (UN) is a promising nuclear fuel that combines the advantageous properties of readily used UO 2 and uranium alloys. Various properties of UN have been previously studied using different density functional theory (DFT) methodologies; however, there are still inconsistencies when it comes to the dynamical stability and defect properties of UN. We address these inconsistencies by studying the UN phonons and defect properties using DFT calculations employing two generalized gradient approximation (GGA) exchange-correlation functionals: PBE and AM05, with and without an added on-site Coulomb repulsion term (+U). Furthermore, we investigate the importance of spin-orbit coupling (SOC) when calculating the properties of UN. We use the different methodologies to determine the preference of UN to have antiferromagnetic (AFM) ordering, as seen in experiments, or ferromagnetic (FM) ordering of the uranium spins. We compare the crystallographic properties, density of states, the DFT X-ray photoelectron spectra and phonon dispersions calculated using the different methodologies. We demonstrate that GGA + U reproduces the AFM ordering in UN, but the crystal structure is dynamically unstable. We also show that magnetic ordering is important in finding the lowest energy defective structure, and that SOC has a distinct influence on the energy of the different uranium interstitial defects. Lastly, we discuss the point defect formation energies under U-rich and N-rich conditions, and the stoichiometric formation energies calculated with the different methodologies, providing insight into the observed tendency for forming hypostoichiometric UN.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Introducing constricted variational density functional theory in its relaxed self-consistent formulation (RSCF-CV-DFT) as an alternative to adiabatic time dependent density functional theory for studies of charge transfer transitions

We have applied the relaxed and self-consistent extension of constricted variational density functional theory (RSCF-CV-DFT) for the calculation of the lowest charge transfer transitions in the molecular complex X-TCNE between X = benzene and TCNE = tetracyanoethylene. Use was made of functionals with a fixed fraction (α) of Hartree-Fock exchange ranging from α = 0 to α = 0.5 as well as functionals with a long range correction (LC) that introduces Hartree-Fock exchange for longer inter-electronic distances. A detailed comparison and analysis is given for each functional between the performance of RSCF-CV-DFT and adiabatic time-dependent density functional theory (TDDFT) within the Tamm-Dancoff approximation. It is shown that in this particular case, all functionals afford the same reasonable agreement with experiment for RSCF-CV-DFT whereas only the LC-functionals afford a fair agreement with experiment using TDDFT. We have in addition calculated the CT transition energy for X-TCNE with X = toluene, o-xylene, and naphthalene employing the same functionals as for X = benzene. It is shown that the calculated charge transfer excitation energies are in as good agreement with experiment as those obtained from highly optimized LC-functionals using adiabatic TDDFT. We finally discuss the relation between the optimization of length separation parameters and orbital relaxation in the RSCF-CV-DFT scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗