Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function 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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Reporting Recommended Patch Density from Vehicle Panel Vibration Convergence Studies using both DAF and TBL Fits of the Spatial Correlation Function

Using the patch method to represent the continuous spatial correlation function of a phased pressure field over a structural surface is an approximation. The approximation approaches the continuous function as patches become smaller. Plotting comparisons of the approximation vs the continuous function may provide insight revealing: (1) For what patch size/density should the approximation be very good? (2) What the approximation looks like when it begins to break down? (3) What the approximation looks like when the patch size is grossly too large. Following these observations with a convergence study using one FEM may allow us to see the importance of patch density. We may develop insights that help us to predict sufficient patch density to provide adequate convergence for the intended purpose frequency range of interest

Smith, Andrew M.↗
DOE OSTI2018

Modified Steklov Functions and Numerical Differentiation Formulas

We consider an approximation method based on Steklov functions of the first and second order. We obtain estimates for the norms in the space C of continuous periodic functions and clarify how they connect with numerical differentiation formulas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Analytical Approximations for Beamstrahlung at Very High Energy Electron-Positron Colliders

Among the many effects that occur in beam-beam electron-positron collisions at TeV energies, emission of hard synchrotron radiation, or beamstrahlung, has special importance. Beamstrahlung determines the energy spectrum of the most energetic electrons, positrons, and photons and supplies the initial condition for the calculation of all other QED processes. In this paper, we show that the description of beamstrahlung simplifies in the limit of large quantum parameter $Υ$, which is realized in 10 TeV collider designs. The beamstrahlung spectra for electrons and photons are given in terms of universal functions. We supply approximations to these functions that will be useful for more general studies of the beam-beam interaction at very high energies.

Accelerator Physics (physics.acc-ph)↗
NASA NTRS1988

Efficient approximation concepts using second order information

The application of second derivative information for solving structural optimization problems is considered. In the present method, rather than building approximate nonlinear forms for the objective function and constraints, only linear approximations are used. A separable quadratic approximation of the Lagrangian function is included in the subproblem statement. The method has been successfully used for simple problems that can be solved in closed form, in addition to the sizing optimization of trusses, and it is shown to converge faster than the convex linearization method or the method of moving asymptotes.

Fleury, Claude↗
DOE OSTI2020

Thermal hybrid exchange-correlation density functional for improving the description of warm dense matter

Finite-temperature density functional theory (FT-DFT) has become a standard tool for first principles calculations of the properties of warm dense matter (WDM) relevant to high-energy density physics (HEDP) applications. In this work we introduce theoretical grounds of thermal hybrid exchange-correlation (XC) functionals within the generalized Mermin-Kohn-Sham scheme for an improved description of WDM. Building on the previously developed KSDT (Karasiev-Sjostrom-Dufty-Trickey) [Karasiev et al., Phys. Rev. Lett. 112, 076403 (2014)] local density approximation (LDA) and the KDT16 (Karasiev-Dufty-Trickey 2016) [Karasiev et al., Phys. Rev. Lett. 120, 076402 (2018)] generalized-gradient approximation (GGA) XC free-energy density functionals, we construct a new thermal hybrid XC functional, referred to here as KDT0. The KDT0 model at low temperature reduces to the popular ground-state PBE0 hybrid due to properties of the used KDT16 density functional approximation. Application to static calculations of electronic band gap and band structure at a wide range of temperatures for various systems of interest to HEDP show that KDT0 provides a significant improvement to the lower LDA and GGA rung XC functionals and to the ground-state PBE0 hybrid.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Study of the Electronic Structure, Magnetic and Elastic Properties and Half-Metallic Stability on Variation of Lattice Constants for CoFeCrZ (Z = P, As, Sb) Heusler Alloys

Structural, elastic, magnetic and electronic properties of CoFeCrZ (Z = P, As, Sb) Heusler alloys in their YI-type structure have been computed by density functional theory-based WIEN2k code within generalized gradient approximation for exchange correlation functions. Values of formation energy and elastic constants verify the stability of these alloys. True half metallic ferromagnetic behaviour with 100% spin polarization and good band gap in the minority spin are obtained for all the alloys. Magnetic moment of these alloys is 4.00 μ{sub B}, which matches well with the value predicted from Slater-Pauling rule. Half-metallicity is maintained over a wide range of lattice constants making these alloys promising for spintronics device applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1993

Space transportation vehicle design evaluation using saturated designs

An important objective in the preliminary design and evaluation of space transportation vehicles is to find the best values of design variables that optimize the performance characteristic (e.g. dry weight). For a given configuration, the vehicle performance can be determined by the use of complex sizing and performance evaluation computer programs. These complex computer programs utilize iterative algorithms and they are generally too expensive and/or difficult to use directly in multidisciplinary design optimization. An alternative is to use response surface methodology (RSM) and obtain quadratic polynomial approximations to the functional relationships between performance characteristics and design variables. In RSM, these approximation models are then used to determine optimum design parameter values and for rapid sensitivity studies. Constructing a second-order model requires that 'n' design parameters be studied at least at 3 levels (values) so that the coefficients in the model can be estimated. There, 3(n) factorial experiments (point designs or observations) may be necessary. For small values of 'n' such as two or three, this design works well. However, when a large number of design parameters are under study, the number of design points required for a full-factorial design may become excessive. Fortunately, these quadratic polynomial approximations can be obtained by selecting an efficient design matrix using central composite designs (CCD) from design of experiments theory. Each unique point design from the CCD matrix is then conducted using computerized analysis tools (e.g. POST, CONSIZ, etc.). In the next step, least squares regression analysis is used to calculate the quadratic polynomial coefficients from the data. However, in some multidisciplinary applications involving a large number of design variables and several disciplines, the computerized performance synthesis programs may get too time consuming and expensive to run even with the use of efficient central composite designs. In such cases, it may be preferable to keep the number of design points to an absolute minimum and trade some model accuracy with cost. For this purpose, another class of experimental designs, called saturated D-optimal designs may be utilized for generating a matrix of vehicle designs. A design is called saturated when the number of design points is exactly equal to the number of terms in the model to be fitted. As a result, saturated designs require the absolute minimum number of design points ((n+1)(n+2)/s) to estimate the quadratic polynomial model coefficients. Saturated designs can be generated using the D-optimality criterion. A good saturated design should give rise to least squares estimates with minimum generalized variance.

Unal, Resit↗
DOE OSTI2023

Introducing the embedded random phase approximation: H 2 dissociative adsorption on Cu(111) as an exemplar

The random phase approximation (RPA) as a means of treating electron correlation recently has been shown to outperform standard density functional theory (DFT) approximations in a variety of cases. However, the computational cost of the RPA is substantially more than DFT, especially when aiming to study extended surfaces. Properly accounting for sufficient surface ensemble size, Brillouin zone sampling, and vacuum separation of periodic images in standard periodic-planewave-based DFT code raises the cost to achieve converged results. Here, we show that sub-system embedding schemes enable use of the RPA for modeling heterogeneous reactions at reduced computational cost. Further, we explore two different embedded RPA (emb-RPA) approaches, periodic emb-RPA and cluster emb-RPA. We use the (experimentally and theoretically) well-studied H 2 dissociative adsorption on Cu(111) as our exemplar, and first perform full periodic RPA calculations as a benchmark. The full RPA results match well the semi-empirical barrier fit to experimental observables and others derived from high-level computations, e.g., from recent embedded n-electron valence second order perturbation theory [Zhao et al., J. Chem. Theory Comput. 16(11), 7078–7088 (2020)] and quantum Monte Carlo [Doblhoff-Dier et al., J. Chem. Theory Comput. 13(7), 3208–3219 (2017)] simulations. Among the two emb-RPA approaches tested, the cluster emb-RPA accurately reproduces the energy profile (maximum error of 50 meV along the reaction pathway) while reducing the computational cost by approximately two orders of magnitude. We therefore expect that the embedded cluster approach will enable wider RPA implementation in heterogeneous catalysis.

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

Three-dimensional analysis of anisotropic spatially reinforced structures

The material-adaptive three-dimensional analysis of inhomogeneous structures based on the meso-volume concept and application of deficient spline functions for displacement approximations is proposed. The general methodology is demonstrated on the example of a brick-type mosaic parallelepiped arbitrarily composed of anisotropic meso-volumes. A partition of each meso-volume into sub-elements, application of deficient spline functions for a local approximation of displacements and, finally, the use of the variational principle allows one to obtain displacements, strains, and stresses at anypoint within the structural part. All of the necessary external and internal boundary conditions (including the conditions of continuity of transverse stresses at interfaces between adjacent meso-volumes) can be satisfied with requisite accuracy by increasing the density of the sub-element mesh. The application of the methodology to textile composite materials is described. Several numerical examples for woven and braided rectangular composite plates and stiffened panels under transverse bending are considered. Some typical effects of stress concentrations due to the material inhomogeneities are demonstrated.

Bogdanovich, Alexander E.↗
DOE OSTI2022

Statistics of Green's functions on a disordered Cayley tree and the validity of forward scattering approximation

The accuracy of the forward scattering approximation for two-point Green's functions of the Anderson localization model on the Cayley tree is studied. A relationship between the moments of the Green's function and the largest eigenvalue of the linearized transfer-matrix equation is proved in the framework of the supersymmetric functional-integral method. The new large-disorder approximation for this eigenvalue is derived and its accuracy is established. Using this approximation the probability distribution of the two-point Green's function is found and compared with that in the forward scattering approximation (FSA). It is shown that FSA overestimates the role of resonances and thus the probability for the Green's function to be significantly larger than its typical value. The error of FSA increases with increasing the distance between points in a two-point Green's function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

The Sobolev approximation for line formation with partial frequency redistribution

Attention is given to the formation of a spectral line in a uniformly expanding infinite medium in the Sobolev approximation, with emphasis on the various mechanisms for frequency redistribution. Numerical and analytic solutions of the transfer equation are presented of a number of redistribution functions and their approximations, including type I and type II partial redistribution, coherent scattering and complete redistribution, and the Fokker-Planck and uncorrelated approximation to the R sub II function. The solutions for the mean intensity are shown to depend very much on the type of redistribution mechanism, while for the frequency-weighted mean intensity, which enters the rate equations, this dependence is weak. It is inferred that use of Sobolev escape probabilities based on complete redistribution can be an adequate approximation for many calculations for which only the radiative excitation rates are needed.

Hummer, D. G.↗
DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Mappings and accuracy for Chebyshev pseudo-spectral approximations

The effect of mappings on the approximation, by Chebyshev collocation, of functions which exhibit localized regions of rapid variation is studied. A general strategy is introduced whereby mappings are adaptively constructed which map specified classes of rapidly varying functions into low order polynomials which can be accurately approximated by Chebyshev polynomial expansions. A particular family of mappings constructed in this way is tested on a variety of rapidly varying functions similar to those occurring in approximations. It is shown that the mapped function can be approximated much more accurately by Chebyshev polynomial approximations than in physical space or where mappings constructed from other strategies are employed.

Bayliss, Alvin↗
NASA NTRS1992

Mappings and accuracy for Chebyshev pseudo-spectral approximations

The effect of mappings on the approximation, by Chebyshev collocation, of functions which exhibit localized regions of rapid variation is studied. A general strategy is introduced whereby mappings are adaptively constructed which map specified classes of rapidly varying functions into low order polynomials which can be accurately approximated by Chebyshev polynomial expansions. A particular family of mappings constructed in this way is tested on a variety of rapidly varying functions similar to those occurring in approximations. It is shown that the mapped function can be approximated much more accurately by Chebyshev polynomial approximations than in physical space or where mappings constructed from other strategies are employed.

Bayliss, Alvin↗
DOE OSTI2019

Can Convolution and Deconvolution Be Used as Tools for Modeling Multicomponent, Mixing-Limited Reaction Networks?

Transfer functions (i.e., convolutions) have been powerful tools for both forward modeling and reconstructing past states of transport systems in a wide range of subsurface flow problems. The advantage of the transfer function is that it is a simple alternative to complicated, distributed parameter models of flow and transport, but the majority of applications of transfer functions in hydrology have been limited to relatively simple cases, like passive tracers or first-order decay. The central question evaluated in this note is whether or not multicomponent mixing-limited reactive transport can be represented within a transfer function framework. Our examples consider forward-in-time predictions and backward-in-time reconstructions of a carbonate system that represents the intrusion of seawater into a freshwater aquifer. The main result is that accurate forward-in-time and backward-in-time models are developed by posing the problem in terms of conservative components. As with all convolution-based methods, the results are sensitive to errors and/or noise in the input functions, but we show that smoothed approximations of the requisite functions provide good representations of transport. Given the vast unknowns in any subsurface transport problem, such generalized, reactive transfer function models may have yet unexplored advantages when the trade-offs between overall computational cost, accuracy, and uncertainty are explored in more detail.

42 ENGINEERING↗