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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Filling the black hole mass gap: Avoiding pair instability in massive stars through addition of nonnuclear energy

In standard stellar evolution, stars with masses ranging from approximately 150 to 240M ⊙ are expected to evolve to a pair instability supernova with no black hole (BH) remnant. This evolutionary behavior leads to a predicted gap in the black hole mass function from approximately 50 to 140M ⊙ . Yet the LIGO and Virgo Collaborations recently discovered black holes of masses 66M ⊙ and 85M ⊙ in the gravitational wave event GW190521. Here, we propose a new method to populate the BH mass gap. If an energy source is added throughout the star in addition to nuclear fusion, it is possible for the altered evolution to avoid the complete destruction of a pair instability supernova, and instead a BH remnant is left behind. An example of an extra energy source is dark matter annihilation within the star, but our results hold more generally. We show this phenomenon by exploring the effect of adding an energy source independent of temperature and density to a 180M ⊙ star, using the MESA one-dimensional stellar evolution software. If ~50% of the star’s energy is due to this new source, the star is capable of avoiding the pair instability entirely and evolving towards a core-collapse supernova and ultimately a BH remnant with mass ~ 120M ⊙ .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Double-hybrid density-functional theory with meta-generalized-gradient approximations

We extend the previously proposed one-parameter double-hybrid density-functional theory [K. Sharkas, J. Toulouse, and A. Savin, J. Chem. Phys. 134, 064113 (2011)] to meta-generalized-gradient-approximation (meta-GGA) exchange-correlation density functionals. We construct several variants of one-parameter double-hybrid approximations using the Tao-Perdew-Staroverov-Scuseria (TPSS) meta-GGA functional and test them on test sets of atomization energies and reaction barrier heights. The most accurate variant uses the uniform coordinate scaling of the density and of the kinetic energy density in the correlation functional, and improves over both standard Kohn-Sham TPSS and second-order Møller-Plesset calculations.

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

Approximations to the plasma dispersion function

Linear wave propagation in hot collisionless plasmas is described by the linearized Vlasov and Maxwell equations. In uniform media, the utilization of spatial and temporal transforms of those equations leads to the consideration of integrals of the Hilbert transform type. Analysis and comparison of two simple approximations are provided based on the utilization of resonance velocity distributions. Application is then made to the Landau and whistler waves, along with a discussion of the results, and commentary on possible improvements.

Brinca, A. L.↗
DOE OSTI2023

Generalized moving least squares vs. radial basis function finite difference methods for approximating surface derivatives

Approximating differential operators defined on two-dimensional surfaces is an important problem that arises in many areas of science and engineering. Over the past ten years, localized meshfree methods based on generalized moving least squares (GMLS) and radial basis function finite differences (RBF-FD) have been shown to be effective for this task as they can give high orders of accuracy at low computational cost, and they can be applied to surfaces defined only by point clouds. However, there have yet to be any studies that perform a direct comparison of these methods for approximating surface differential operators (SDOs). The first purpose of this work is to fill that gap. For this comparison, we focus on an RBF-FD method based on polyharmonic spline kernels and polynomials (PHS+Poly) since they are most closely related to the GMLS method. Additionally, we use a relatively new technique for approximating SDOs with RBF-FD called the tangent plane method since it is simpler than previous techniques and natural to use with PHS+Poly RBF-FD. Further, the second purpose of this work is to relate the tangent plane formulation of SDOs to the local coordinate formulation used in GMLS and to show that they are equivalent when the tangent space to the surface is known exactly. The final purpose is to use ideas from the GMLS SDO formulation to derive a new RBF-FD method for approximating the tangent space for a point cloud surface when it is unknown. For the numerical comparisons of the methods, we examine their convergence rates for approximating the surface gradient, divergence, and Laplacian as the point clouds are refined for various parameter choices. We also compare their efficiency in terms of accuracy per computational cost, both when including and excluding setup costs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Exact finite elements for conduction and convection

An approach for developing exact one dimensional conduction-convection finite elements is presented. Exact interpolation functions are derived based on solutions to the governing differential equations by employing a nodeless parameter. Exact interpolation functions are presented for combined heat transfer in several solids of different shapes, and for combined heat transfer in a flow passage. Numerical results demonstrate that exact one dimensional elements offer advantages over elements based on approximate interpolation functions.

Thornton, E. A.↗
NASA NTRS1981

Exact finite elements for conduction and convection

An appproach for developing exact one dimensional conduction-convection finite elements is presented. Exact interpolation functions are derived based on solutions to the governing differential equations by employing a nodeless parameter. Exact interpolation functions are presented for combined heat transfer in several solids of different shapes, and for combined heat transfer in a flow passage. Numerical results demonstrate that exact one dimensional elements offer advantages over elements based on approximate interpolation functions. Previously announced in STAR as N81-31507

Thornton, E. A.↗
NASA NTRS1995

Exponentially accurate approximations to piece-wise smooth periodic functions

A family of simple, periodic basis functions with 'built-in' discontinuities are introduced, and their properties are analyzed and discussed. Some of their potential usefulness is illustrated in conjunction with the Fourier series representations of functions with discontinuities. In particular, it is demonstrated how they can be used to construct a sequence of approximations which converges exponentially in the maximum norm to a piece-wise smooth function. The theory is illustrated with several examples and the results are discussed in the context of other sequences of functions which can be used to approximate discontinuous functions.

Greer, James↗
DOE OSTI2016

Calibration of photoelastic modulator based dichrometers: maintaining constant phase across the spectrum

Photoelastic modulators can alter the polarization state of a beam of ultraviolet, visible or infrared photons by means of periodic stress-induced differences in the refractive index of a transparent material that forms the optical element of the device and is isotropic in the absence of stress. Furthermore, they have found widespread application in instruments that characterize or alter the polarization state of a beam in fields as diverse as astronomy, structural biology, materials science and ultraviolet lithography for the manufacture of nano-scale integrated circuits. Measurement of circular dichroism, the differential absorption of left- and right circularly polarized light, and of strain-induced birefringence of optical components are major applications. Instruments using synchrotron radiation and photoelastic modulators with CaF 2 optical elements have extended circular dichroism measurements down to wavelengths of about 130 nm in the vacuum ultraviolet. Maintaining a constant phase shift between two orthogonal polarization states across a spectrum requires that the amplitude of the modulated stress be changed as a function of wavelength. For commercially available photoelastic modulators, the voltage that controls the amplitude of modulation required to produce a specified phase shift, which is a surrogate for the stress modulation amplitude, has been shown to be an approximately linear function of wavelength in the spectral region where the optical element is transparent. But, extrapolations of such straight lines cross zero voltage at a non-zero wavelength, not at zero-wavelength. For modulators with calcium fluoride and fused silica optical elements, the zero-crossing wavelength is always in the spectral region where the optical element of the modulator strongly absorbs the incident radiation, and at a wavelength less than the longest-wavelength apparent resonance deduced from experimental values of the refractive index fit to the Sellmeier equation. Using a model that relates the refractive indices of a stressed optical element to the refractive index of its unstressed state, an expression for the modulator control voltage was derived that closely fits the experimental data. Our result provides a theoretical rational for the apparently linear constant-phase programming voltage, and thus provides theoretical backing for the calibration procedure frequently used for these modulators. Lastly there are other factors that can influence the calibration of a photoelastic modulator, including temperature and atmospheric pressure, are discussed briefly.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1977

Monte Carlo turbulence simulation

The paper describes turbulence simulation experiments based on the principles of control system theory, that is, the construction of a system characterized by a system function such that upon exciting the system with prescribed noise processes the output of the system is a realization of a random processing the desired statistical attributes of turbulence. An experimental autocorrelation of Jimsphere measurements of wind velocity was approximated to simulate turbulent wind. From the approximate autocorrelation function, the required system function is obtained, and a discrete time system is designed. Another method of simulation is to solve the convolution integral by filter techniques. Other methods include discrete Fourier simulation and self-similar simulation.

Fichtl, G. H.↗
DOE OSTI2016

Tight-binding approximations to time-dependent density functional theory — A fast approach for the calculation of electronically excited states

We propose a new method of calculating electronically excited states that combines a density functional theory based ground state calculation with a linear response treatment that employs approximations used in the time-dependent density functional based tight binding (TD-DFTB) approach. The new method termed time-dependent density functional theory TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive TD-DFT calculations. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB.

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

Analysis of the spin Hall effect in CuIr alloys: Combined approach of density functional theory and Hartree-Fock approximation

We analyze the spin Hall effect in CuIr alloys in theory by the combined approach of the density functional theory (DFT) and Hartree-Fock (HF) approximation. The spin Hall angle (SHA) is obtained to be negative without the local correlation effects. After including the local correlation effects of the 5d orbitals of Ir impurities, the SHA becomes positive with realistic correlation parameters and consistent with experiment [Niimi et al., Phys. Rev. Lett. 106, 126601 (2011)]. Moreover, our analysis shows that the DFT + HF approach is a convenient and general method to study the influence of local correlation effects on the spin Hall effect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Weighted nodal domain averages of eigenstates for quantum Monte Carlo and beyond

In this report we study the nodal properties of many-body eigenstates of stationary Schrödinger equation that affect the accuracy of real-space quantum Monte Carlo calculations. In particular, we introduce weighted nodal domain averages that provide a new probe of nodal surfaces beyond the usual expectations. Particular choices for the weight function reveal, for example, that the difference between two arbitrary fermionic eigenvalues is given by the nodal hypersurface integrals normalized by overlaps with the bosonic ground state of the given Hamiltonian. Noninteracting and fully interacting Be atom with corresponding almost exact and approximate wave functions are used to illustrate several aspects of these concepts. Variational formulations that employ different weights are proposed for prospective improvement of nodes in variational and fixed-node diffusion Monte Carlo calculations.

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

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions in the presence of a generalized incident local radiation field computed by Cloudy. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins (0.5-1, 1-4, and 13-16Ry) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors ≳10 times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions computed by Cloudy in the presence of a generalized incident local radiation field. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins ($0.5-1, 1-4$, and $13-16 \, \mathrm{Ry}$) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors $\gtrsim 10$ times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Accurate frozen core approximation for all-electron density-functional theory

We implement and benchmark the frozen core approximation, a technique commonly adopted in electronic structure theory to reduce the computational cost by means of mathematically fixing the chemically inactive core electron states. The accuracy and efficiency of this approach are well controlled by a single parameter, the number of frozen orbitals. Explicit corrections for the frozen core orbitals and the unfrozen valence orbitals are introduced, safeguarding against seemingly minor numerical deviations from the assumed orthonormality conditions of the basis functions. A speedup of over twofold can be achieved for the diagonalization step in all-electron density-functional theory simulations containing heavy elements, without any accuracy degradation in terms of the electron density, total energy, and atomic forces. This is demonstrated in a benchmark study covering 103 materials across the Periodic Table and a large-scale simulation of CsPbBr3 with 2560 atoms. Our study provides a rigorous benchmark of the precision of the frozen core approximation (sub-meV per atom for frozen core orbitals below −200 eV) for a wide range of test cases and for chemical elements ranging from Li to Po. The algorithms discussed here are implemented in the open-source Electronic Structure Infrastructure software package.

Yu, Victor Wen-zhe↗
DOE OSTI2023

Imposing correct jellium response is key to predict the density response by orbital-free DFT

Orbital-free density functional theory constitutes a computationally highly effective tool for modeling electronic structures of systems ranging from room-temperature materials to warm dense matter. Its accuracy critically depends on the employed kinetic energy (KE) density functional, which has to be supplied as an external input. In this work we consider several nonlocal and Laplacian-level KE functionals and use an external harmonic perturbation to compute the static density response at T = 0 K in the linear and beyond-linear response regimes. We test for the satisfaction of exact conditions in the limit of uniform densities and for how approximate KE functionals reproduce the density response of realistic materials (e.g., Al and Si) against the Kohn-Sham DFT reference, which employs the exact KE. The results illustrate that several functionals violate exact conditions in the uniform electron gas (UEG) limit. We find a strong correlation between the accuracy of the KE functionals in the UEG limit and in the strongly inhomogeneous case. This empirically demonstrates the importance of imposing the limit of UEG response for uniform densities and validates the use of the Lindhard function in the formulation of kernels for nonlocal functionals. This conclusion is substantiated by additional calculations for bulk aluminum (Al) with a face-centered cubic (fcc) lattice and silicon (Si) with an fcc lattice, body-centered cubic (bcc) lattice, and semiconducting crystal diamond state. The analysis of fcc Al, and fcc as well as bcc Si data follows closely the conclusions drawn for the UEG, allowing us to extend our conclusions to realistic systems that are subject to density inhomogeneities induced by ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Algebraic grid generation for wing-fuselage bodies

An algebraic procedure for the generation of boundary-fitted grids about wing-fuselage configurations is presented. A wing-fuselage configuration is specified by cross sections and mathematically represented by Coons' patches. A configuration is divided into sections so that several grid blocks that either adjoin each other or partially overlap each other can be generated, and each grid has six surfaces that map into a computational cube. Grids are first determined on the six boundary surfaces and then in the interior. Grid curves that are on the surface of the configuration are derived using plane-patch intersections, and single-valued functions relating approximate arc lengths along the curves to computational coordinates define the distribution of grid points. The two-boundary technique and transfinite interpolation are used to determine the boundary surface grids that are not on the configuration, and transfinite interpolation with linear blending functions is used to determine the interior grids.

Smith, R. E.↗
DOE OSTI2022

Noise effects on Padé approximants and conformal maps

Here, we analyze the properties of Padé and conformal map approximants for functions with branch points, in the situation where the expansion coefficients are only known with finite precision or are subject to noise. We prove that there is a universal scaling relation between the strength of the noise and the expansion order at which Padé or the conformal map breaks down. We illustrate this behavior with some physically relevant model test functions and with two non-trivial physical examples where the relevant Riemann surface has complicated structure.

97 MATHEMATICS AND COMPUTING↗