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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Weighted greedy-optimal design of computer experiments for kernel-based and Gaussian process model emulation and calibration

This article is concerned with the approximation of high-dimensional functions by kernel-based methods. Motivated by uncertainty quantification, which often necessitates the construction of approximations that are accurate with respect to a probability density function of random variables, we aim at minimizing the approximation error with respect to a weighted $L^p$-norm. We present a greedy procedure for designing computer experiments based upon a weighted modification of the pivoted Cholesky factorization. The method successively generates nested samples with the goal of minimizing error in regions of high probability. Numerical experiments validate that this new importance sampling strategy is superior to other sampling approaches, especially when used with non-product probability density functions. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

An Efficient, Globally Convergent Method for Optimization Under Uncertainty Using Adaptive Model Reduction and Sparse Grids

We present a new approach to efficiently solve optimization problems constrained by partial differential equations (PDEs) with uncertain coefficients. The method leverages two sources of inexactness that trade accuracy for speed: (1) stochastic collocation based on dimension-adaptive sparse grids (SGs), which approximates the stochastic objective function with a limited number of quadrature nodes, and (2) projection-based reduced-order models (ROMs), which generate efficient approximations to PDE solutions. These two sources of inexactness lead to inexact objective function and gradient evaluations, which are managed by a trust-region method that guarantees global convergence by adaptively refining the SG and ROM until a proposed error indicator drops below a tolerance specified by trust-region convergence theory. A key feature of the proposed method is that the error indicator---which accounts for errors incurred by both the SG and ROM---must be only an asymptotic error bound, i.e., a bound that holds up to an arbitrary constant that need not be computed. This allows the approach to be applicable to a wide range of problems, including those where sharp, computable error bounds are not available; this distinguishes the proposed method from previous works. Numerical experiments performed on a model problem from optimal flow control under uncertainty verify global convergence of the method and demonstrate the method's ability to outperform previously proposed alternatives.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Full self-consistency versus quasiparticle self-consistency in diagrammatic approaches: Exactly solvable two-site Hubbard model

Self-consistent solutions of Hedin's equations (HE) for the two-site Hubbard model (HM) have been studied. They have been found for three-point vertices of increasing complexity (Γ = 1 (GW approximation), Γ₁ from the first-order perturbation theory, and the exact vertex Γ E ). Comparison is made between the cases when an additional quasiparticle (QP) approximation for Green's functions is applied during the self-consistent iterative solving of HE and when QP approximation is not applied. Results obtained with the exact vertex are directly related to the present open question—which approximation is more advantageous for future implementations, GW + DMFT or QPGW + DMFT. It is shown that in a regime of strong correlations only the originally proposed GW + DMFT scheme is able to provide reliable results. Vertex corrections based on Perturbation Theory systematically improve the GW results when full self-consistency is applied. The application of QP self-consistency combined with PT vertex corrections shows similar problems to the case when the exact vertex is applied combined with QP sc. An analysis of Ward Identity violation is performed for all studied in this work's approximations and its relation to the general accuracy of the schemes used is provided.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2003

Long Focal Length Large Mirror Fabrication System

The goals of this ambitious program are (1) to develop systems to make large superpolished optical mirrors, (2) to develop low scatter polishing techniques using centrifugal elutriation, (3) to develop a means of measuring scatter at any point on the mirror, (4) to polish a Hindle sphere to measure the optical figure of a one meter diameter convex mandrel, and (5) to fabricate low scatter, large adaptive optic graphite filled, cyanate ester replica transfer mirrors using these mandrels. Deliverables are a 30 cm diameter superpolished composite AO mirror. We fabricated a 1/3rd meter superpolished zerodur flat mandrel and with the support of our major subcontractor, Composite Mirror Applications Inc (CMA) we have demonstrated a 30 cm lightweight cyanate ester mirror with an rms microroughness between 0.6 and 0.8 nm and 8 faceplate influence function of 5 cm. The influence function was chosen to be comparable to the atmospheric correlation coefficient r(sub 0) which is about 5 cm at sea level. There was no print-thru of the graphite fibers in the cyanate ester surface (the bane of many previous efforts to use cyanate ester mirrors). Our subcontractor has devised a means for developing a 30-50 nm thick layer of graphite free pure ester resin on the surface of the mirrors. This graphite fiber filled material has a thermal expansion coefficient in the 10(exp -8) centimeter per Kelvin range (the same range of expansion coefficient as Zerodur and ULE glasses) and does not take up water and swell, so it is a nearly ideal mirror material in these areas. Unfortunately for these 0.8mm thick faceplates, the number of plies is not enough to result in isometric coverage. Isolated figure irregularities can appear, making it necessary to go to thicker faceplates. The influence function will then only approximate the length of r(sub 0), at higher altitudes or longer wavelengths. The influence function goes as the cube of the thickness, so we are now making a faceplate optimized for an infrared wavelength, possibly that used for the AO system of the Keck telescope, instead of 0.63 microns. We have polished a 55 cm diameter mandrel to better than 1/20th wave optical figure in the visible using centrifugal elutriation. CMA has just told us that it needs to retool to get optimum mirror faceplate quality in this size, so implementing the 55 cm AO mirror may be delayed somewhat. We expect to complete our 1/3 rd meter AO mirror on time using novel piezoelectric actuators with a throw of one micrometer per volt, as compared to 0.005 micrometers per volt for conventional piezoelectric actuators. We will then demonstrate its AO performance interferometrically.

Bennett, H. E.↗
DOE OSTI2021

How well do self-interaction corrections repair the overestimation of static polarizabilities in density functional calculations?

Herein, we examine the effect of removing self-interaction error (SIE) on the calculation of molecular polarizabilities in the local spin density (LSDA) and generalized gradient approximations (GGA). To this end, we utilize a database of 132 molecules taken from a recent benchmark study [Hait and Head-Gordon, Phys. Chem. Chem. Phys. , 2018, 20 , 19800] to assess the influence of SIE on polarizabilities by comparing results with accurate reference data. Our results confirm that the general overestimation of molecular polarizabilities by these density functional approximations can be attributed to SIE. However, removing SIE using the Perdew–Zunger self-interaction-correction (PZ-SIC) method, implemented using the Fermi–Löwdin Orbital SIC approach, leads to an underestimation of molecular polarizabilities, showing that PZ-SIC overcorrects when combined with LSDA or GGA. Application of a recently proposed locally scaled SIC [Zope, et al., J. Chem. Phys. , 2019, 151 , 214108] is found to provide more accurate polarizabilities. We attribute this to the ability of the local scaling scheme to selectively correct for SIE in the regions of space where the correction is needed most.

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

High-energy heavy-ion scattering and the optical phase shift function

Within the framework of the Glauber approximation, an exact expression for the optical phase shift function is derived which takes into account the center-of-mass correlation function. We calculate total and differential cross sections using the first-order term for the shift function and compare them with the usual (first-order) optical limit results and, for the lighter nuclei, with the exact Glauber results. The nucleus-nucleus elastic scattering amplitude does not exhibit the large-q divergence which characterizes the usual optical limit. The obtained results improve the calculated total and differential cross sections dramatically for light nuclei, and significantly for medium and heavy nuclei.

Franco, V.↗
NASA NTRS1972

Methods of computing vocabulary size for the two-parameter rank distribution

A summation method is described for computing the vocabulary size for given parameter values in the 1- and 2-parameter rank distributions. Two methods of determining the asymptotes for the family of 2-parameter rank-distribution curves are also described. Tables are computed and graphs are drawn relating paris of parameter values to the vocabulary size. The partial product formula for the Riemann zeta function is investigated as an approximation to the partial sum formula for the Riemann zeta function. An error bound is established that indicates that the partial product should not be used to approximate the partial sum in calculating the vocabulary size for the 2-parameter rank distribution.

Edmundson, H. P.↗
DOE OSTI2016

A well-scaling natural orbital theory

Here, we introduce an energy functional for ground-state electronic structure calculations. Its variables are the natural spin-orbitals of singlet many-body wave functions and their joint occupation probabilities deriving from controlled approximations to the two-particle density matrix that yield algebraic scaling in general, and Hartree–Fock scaling in its seniority-zero version. Results from the latter version for small molecular systems are compared with those of highly accurate quantum-chemical computations. The energies lie above full configuration interaction calculations, close to doubly occupied configuration interaction calculations. Their accuracy is considerably greater than that obtained from current density-functional theory approximations and from current functionals of the oneparticle density matrix.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2019

Revised M11 Exchange-Correlation Functional for Electronic Excitation Energies and Ground-State Properties

The ability of Kohn-Sham density functional theory (KS-DFT) to accurately predict various types of electronic excitation energies with (necessarily approximate) exchange-correlation functionals faces several challenges. Chief among these is that valence excitations are usually inherently multiconfigurational and therefore best treated by functionals with local exchange, whereas Rydberg and charge transfer excitations are often better treated with nonlocal exchange. The question arises of whether one can optimize a functional such that all three kinds of excitations (valence, Rydberg, and charge transfer – including long-range charge transfer) are treated in a balanced and accurate way. The goal of the present work is to try to answer that question and then to optimize a functional with the best possible balanced behavior. Of the variety of functional types available, we select range-separated hybrid meta functionals because (i) range separation allows the percentage of Hartree–Fock (HF) exchange to change with interelectronic separation, and therefore, one can have 100% HF exchange at large interelectronic separations, which gives good performance for long-range charge-transfer excitations, while the range separation allows one to simultaneously have smaller values of HF exchange at small and intermediate inter-electronic separations, which give good performance for valence and Rydberg excitations and (ii) meta functionals allow one to obtain better accuracy with high HF exchange than is possible with functionals whose local part depends only on spin densities and their gradients. Furthermore, this work starts with the range-separated hybrid meta functional, M11, and re-optimizes it (with stronger smoothness restraints) against electronic excitation energies and ground-state properties to obtain a new functional called revM11 that gives good performance for all three types of electronic excitations and at the same time gives very good predictions across-the-board for ground-state properties.

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

Electronic energy states

One-electron wave functions are reviewed and approximate solutions of two-electron systems are given in terms of these one-electron functions. The symmetry effects associated with electron spin are reviewed and the effects of electron exchange on energy levels of the two-electron system are given. The coupling of electronic orbital and spin angular momentum is considered next and the Lande interval rule for Russell-Saunders or LS coupling is derived. The configurations possible for various multi-electron LS couplings are enumerated (examples from the first two rows of the periodic table are given), and the meaning of the spectroscopic nomenclature is discussed, particularly with respect to the degeneracies of the electron states involved. Next the nomenclature, symmetries, and degeneracies for electron states of diatomic molecules are discussed, and some examples for N2, O2, and NO are presented. The electronic partition functions and derivative thermodynamic properties are expressed in terms of these energies and degeneracies, and examples are given for some of the simple gas species encountered in the earth's atmosphere.

Source record↗
DOE OSTI2023

Two-pion femtoscopic correlations in Be+Be collisions at $\sqrt{s_{\text {NN}}} = 16.84$ GeV measured by the NA61/SHINE at CERN

This paper reports measurements of two-pion femtoscopic correlations in Be+Be collisions at a beam momentum of 150$A\,\hbox {GeV}\!/\!c$ (energy available in the center-of-mass system for nucleon pair $\sqrt{s_{\text {NN}}} = 16.84$ GeV) by the NA61/SHINE experiment at the CERN SPS accelerator. The obtained momentum space correlation functions can be well described by a Lévy distributed source model. The transverse mass dependence of the Lévy source parameters is presented, and their possible theoretical interpretations are discussed. The results show that the Lévy exponent $\alpha $ is approximately constant as a function of $m_{\text {T}}$ , and far from both the Gaussian case of $\alpha = 2$ or the conjectured value at the critical endpoint, $\alpha = 0.5$. The radius scale parameter R shows a slight decrease in $m_{\text {T}}$ , which can be explained as a signature of transverse flow. Finally, an approximately constant trend of the intercept parameter $\lambda $ as a function of $m_{\text {T}}$ was observed, similar to previous NA44 S + Pb results (obtained with a Gaussian approximation, but unlike RHIC results).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Exchange interaction between the triplet exciton and the localized spin in copper-phthalocyanine

Triplet excitonic state in the organic molecule may arise from a singlet excitation and the following inter-system crossing. Especially for a spin-bearing molecule, an exchange interaction between the triplet exciton and the original spin on the molecule can be expected. In this paper, such exchange interaction in copper-phthalocyanine (CuPc, spin-1/2 ) was investigated from first-principles by using density-functional theory within a variety of approximations to the exchange correlation, ranging from local-density approximation to long-range corrected hybrid-exchange functional. The magnitude of the computed exchange interaction is in the order of meV with the minimum value (1.5 meV, ferromagnetic) given by the long-range corrected hybrid-exchange functional CAM-B3LYP. This exchange interaction can therefore give rise to a spin coherence with an oscillation period in the order of picoseconds, which is much shorter than the triplet lifetime in CuPc (typically tens of nanoseconds). This implies that it might be possible to manipulate the localized spin on Cu experimentally using optical excitation and inter-system crossing well before the triplet state disappears.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

CrI 3 revisited with a many-body ab initio theoretical approach

CrI 3 has recently been shown to exhibit low-dimensional, long-range magnetic ordering from few layers to single layers of CrI 3 . The properties of CrI 3 bulk and few layered systems are uniquely defined by a combination of short-range intralayer and long-range interlayer interactions; including strong correlations, exchange and spin orbit coupling. Unfortunately, both the long-range van der Waals interactions, which are driven by dynamic, many-body electronic correlations and the competing strong intralayer correlations, present a formidable challenging for the local or semi local mean-field approximations employed in workhorse electronic structure approaches like density functional theory. In this work, we employ a sophisticated many-body approach that can simultaneously describe long- and short-range correlations. Herein, we establish that the fixed-node diffusion Monte Carlo (FNDMC) method reproduces the experimental interlayer separation distance of bulk CrI 3 for the high-temperature, monoclinic phase, with a reliable prediction of the interlayer binding energy. We subsequently employed, the FNDMC results to benchmark the accuracy of several density functional theory exchange-correlation approximations.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Numerical multi-scaling method to solve the linear stochastic partial differential equations

In this paper, we propose a multi-scaling method to solve a class of the stochastic partial differential equations (SPDEs) of Itô-type. By this method, we approximate the solution of the SPDEs with strongly elliptic operator on a finite domain. The method is based on B-spline wavelet approximations; some of these functions are reshaped to satisfy on boundary conditions exactly. We prove the consistency and stability of the method. Illustrative examples are included to demonstrate the validity and applicability of our proposed technique.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

The delta-Eddington approximation for radiative flux transfer

Simple approximations, like the Eddington, are often incapable of coping with the highly asymmetric phase functions typical of particulate scattering. A simple yet accurate method called the delta-Eddington approximation is proposed for determining monochromatic radiative fluxes in an absorbing-scattering atmosphere. In this method, the governing phase function is approximated by a Dirac delta function forward scatter peak and a two-term expansion of the phase function. The fraction of scattering into the truncated forward peak is taken proportional to the square of the phase function asymmetry factor, which distinguishes the delta-Eddington approximation from others of similar nature. The transmission, reflection, and absorption predicted by the delta-Eddington approximation are compared with doubling method calculations for realistic ranges of optical depth, single-scattering albedo, surface albedo, sun angle and asymmetry factor. The approximation is shown to provide an accurate and analytically simple parameterization of radiation to replace the empirism currently encountered in many general circulation and climate models.

Joseph, J. H.↗
DOE OSTI2018

Supersymmetric Dynamics and Zeta-Functions

Boson, fermion, and super oscillators and (statistical) mechanism of cosmological constant; finite approximation of the zeta-function and fermion factorization of the bosonic statistical sum concidered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Scenario generation for stochastic optimization problems via the sparse grid method

We study the use of sparse grids in the scenario generation (or discretization) problem in stochastic programming problems where the uncertainty is modeled using a continuous multivariate distribution. We show that, under a regularity assumption on the random function involved, the sequence of optimal objective function values of the sparse grid approximations converges to the true optimal objective function values as the number of scenarios increases. The rate of convergence is also established. We treat separately the special case when the underlying distribution is an affine transform of a product of univariate distributions, and show how the sparse grid method can be adapted to the distribution by the use of quadrature formulas tailored to the distribution. We numerically compare the performance of the sparse grid method using different quadrature rules with classic quasi-Monte Carlo (QMC) methods, optimal rank-one lattice rules, and Monte Carlo (MC) scenario generation, using a series of utility maximization problems with up to 160 random variables. The results show that the sparse grid method is very efficient, especially if the integrand is sufficiently smooth. In such problems the sparse grid scenario generation method is found to need several orders of magnitude fewer scenarios than MC and QMC scenario generation to achieve the same accuracy. As a result, it is indicated that the method scales well with the dimension of the distribution--especially when the underlying distribution is an affine transform of a product of univariate distributions, in which case the method appears scalable to thousands of random variables.

97 MATHEMATICS AND COMPUTING↗