Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximation methods”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

On the accuracy of density functional theory and wave function methods for calculating vertical ionization energies

The best practice in computational methods for determining vertical ionization energies (VIEs) is assessed, via reference to experimentally determined VIEs that are corroborated by highly accurate coupled-cluster calculations. These reference values are used to benchmark the performance of density functional theory (DFT) and wave function methods: Hartree-Fock theory, second-order Møller-Plesset perturbation theory, and Electron Propagator Theory (EPT). The core test set consists of 147 small molecules. An extended set of six larger molecules, from benzene to hexacene, is also considered to investigate the dependence of the results on molecule size. The closest agreement with experiment is found for ionization energies obtained from total energy difference calculations. In particular, DFT calculations using exchange-correlation functionals with either a large amount of exact exchange or long-range correction perform best. The results from these functionals are also the least sensitive to an increase in molecule size. In general, ionization energies calculated directly from the orbital energies of the neutral species are less accurate and more sensitive to an increase in molecule size. For the single-calculation approach, the EPT calculations are in closest agreement for both sets of molecules. For the orbital energies from DFT functionals, only those with long-range correction give quantitative agreement with dramatic failing for all other functionals considered. The results offer a practical hierarchy of approximations for the calculation of vertical ionization energies. In addition, the experimental and computational reference values can be used as a standardized set of benchmarks, against which other approximate methods can be compared.

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

Hunting ionic liquids with large electrochemical potential windows

Ionic liquids (ILs) with large electrochemical potential windows (EPWs) are ideally suited to improve the energy and power density of electrical devices for energy storage. However, guidelines are scarcely available for the selection or systematic screening of ILs for practical applications. Here, we present theoretical predictions for the EPW of thousands of cation-anion pairs based on a long-range corrected hybrid density functional theory. The cathodic and anodic potential limits were estimated from the highest occupied molecular orbital and lowest unoccupied molecular orbital energies of individual cations and anions. The approximate method predicts EPWs in reasonable agreement with existing experimental data for several commonly used ionic systems and can be used as a screening tool for a large library of ILs. Five most stable anions, with the cathode stability in the order of B[CN] 4 > PF 6 > BF 4 > BOB > NTf 2 , have been identified. These anions provide the largest EPWs when they are paired with conventional organic cations. The cation stability is reduced upon elongation of the alkyl chain length or branching of the molecular backbone.

catalysis (heterogeneous)↗
DOE OSTI2016

Observation of the spectrally invariant properties of clouds in cloudy-to-clear transition zones during the MAGIC field campaign

We use the spectrally invariant method to study the variability of cloud optical thickness τ and droplet effective radius reff in transition zones (between the cloudy and clear sky columns) observed from Solar Spectral Flux Radiometer (SSFR) and Shortwave Array Spectroradiometer-Zenith (SASZe) during the Marine ARM GPCI Investigation of Clouds (MAGIC) field campaign. The measurements from the SSFR and the SASZe are different, however inter-instrument differences of self-normalized measurements (divided by their own spectra at a fixed time) are small. The spectrally invariant method approximates the spectra in the cloud transition zone as a linear combination of definitely clear and cloudy spectra, where the coefficients, slope and intercept, character-ize the spectrally invariant properties of the transition zone. Simulation results from the SBDART (Santa Barbara DISORT Atmospheric Radiative Transfer) model demonstrate that (1) the slope of the visible band is positively correlated with the cloud optical thickness τ while the intercept of the near-infrared band has high negative cor-relation with the cloud drop effective radius reff even without the exact knowledge of τ; (2) the above relations hold for all Solar Zenith Angle (SZA) and for cloud-contaminated skies. In observations using redundant measure-ments from SSFR and SASZe, we find that during cloudy-to-clear transitions, (a) the slopes of the visible band de-crease, and (b) the intercepts of the near-infrared band remain almost constant near cloud edges. The findings in simulations and observations suggest that, while the optical thickness decreases during the cloudy-to-clear transition, the cloud drop effective radius does not change when cloud edges are approached. These results sup-port the hypothesis that inhomogeneous mixing dominates near cloud edges in the studied cases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Understanding the optical properties of ambient sub- and supermicron particulate matter: results from the CARES 2010 field study in northern California

Abstract. Measurements of the optical properties (absorption, scattering and extinction) of PM1, PM2.5 and PM10 made at two sites around Sacramento, CA, during the June 2010 Carbonaceous Aerosols and Radiative Effects Study (CARES) are reported. These observations are used to establish relationships between various intensive optical properties and to derive information about the dependence of the optical properties on photochemical aging and sources. Supermicron particles contributed substantially to the total light scattering at both sites, about 50 % on average. A strong, linear relationship is observed between the scattering Ångström exponent for PM10 and the fraction of the scattering that is contributed by submicron particles (fsca, PM1) at both sites and with similar slopes and intercepts (for a given pair of wavelengths), suggesting that the derived relationship may be generally applicable for understanding variations in particle size distributions from remote sensing measurements. At the more urban T0 site, the fsca, PM1 increased with photochemical age, whereas at the downwind, more rural T1 site the fsca, PM1 decreased slightly with photochemical age. This difference in behavior reflects differences in transport, local production and local emission of supermicron particles between the sites. Light absorption is dominated by submicron particles, but there is some absorption by supermicron particles ( ∼ 15 % of the total). The supermicron absorption derives from a combination of black carbon that has penetrated into the supermicron mode and from dust, and there is a clear increase in the mass absorption coefficient of just the supermicron particles with increasing average particle size. The mass scattering coefficient (MSC) for the supermicron particles was directly observed to vary inversely with the average particle size, demonstrating that MSC cannot always be treated as a constant in estimating mass concentrations from scattering measurements, or vice versa. The total particle backscatter fraction exhibited some dependence upon the relative abundance of sub- versus supermicron particles; however this was modulated by variations in the median size of particles within a given size range; variations in the submicron size distribution had a particularly large influence on the observed backscatter efficiency and an approximate method to account for this variability is introduced. The relationship between the absorption and scattering Ångström exponents is examined and used to update a previously suggested particle classification scheme. Differences in composition led to differences in the sensitivity of PM2.5 to heating in a thermodenuder to the average particle size, with more extensive evaporation (observed as a larger decrease in the PM2.5 extinction coefficient) corresponding to smaller particles; i.e., submicron particles were generally more susceptible to heating than the supermicron particles. The influence of heating on the particle hygroscopicity varied with the effective particle size, with larger changes observed when the PM2.5 distribution was dominated by smaller particles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Direct Nonlinear Approximation for Security Region Boundary of Integrated Energy Systems: A Polynomial Chaos Expansion Solution

The strong interdependence of electricity, gas, and heating systems can facilitate fault propagation within integrated energy systems (IESs), posing significant challenges to secure operation. This paper proposes a polynomial chaos expansion (PCE)-based approximation method to accurately characterize the IES security region boundary (IES–SRB). By integrating the Karush-Kuhn-Tucker conditions with PCE theory, the IES-SRB approximation problem is reformulated as a set of nonlinear equations concerning the approximation coefficients. Using the Galerkin projection method, these equations are further transformed into a system of projection equations that govern the polynomial approximation coefficients in the IES-SRB approximation. To reduce computational complexity while maintaining high approximation accuracy, a piecewise polynomial approximation method is proposed. Numerical studies on the E39-G20-H6 and E118-G96-H52 IES test systems demonstrate that the proposed method can accurately and effectively construct IES security regions.

Wu, Chenghao [Northeast Electric Power University]↗
DOE OSTI2015

Electronic properties of excess Cr at Fe site in FeCr{sub 0.02}Se alloy

We have studied the effect of substitution of transition-metal chromium (Cr) in excess on Fe sub-lattice in the electronic structure of iron-selenide alloys, FeCr{sub 0.02}Se. In our calculations, we used Korringa-Kohn-Rostoker coherent potential approximation method in the atomic sphere approximation (KKR-ASA-CPA). We obtained different band structure of this alloy with respect to the parent FeSe and this may be reason of changing their superconducting properties. We did unpolarized calculations for FeCr{sub 0.02}Se alloy in terms of density of states (DOS) and Fermi surfaces. The local density approximation (LDA) is used in terms of exchange correlation potential.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Introduction of the Adding and Doubling Method for Solving Bateman Equations for Nuclear Fuel Depletion

This paper introduces and evaluates the Adding and Doubling Method (ADM) for solving the Bateman equations for depletion systems with varying numbers of nuclides and compares it to the Chebyshev Rational Approximation Method (CRAM), both implemented in the reactor physics analysis application Griffin. ADM, when applied to the Crank-Nicolson Finite Difference method, can produce results comparable in accuracy and precision to CRAM with comparable run times for systems with 35 or 297 nuclides. For systems with more than 300 nuclides, the matrix-matrix operations required by ADM are significantly more costly than the matrix-vector operations required by CRAM, making CRAM the more efficient method for systems with large numbers of nuclides. ADM is an accurate method that maintains other advantages over CRAM in that it does not depend on pre-generated coefficients or require complex number operations. ADM also manages to outperform CRAM by a factor of more than 250 in terms of run time for depletion systems that require multiple Bateman solves while the depletion matrix and time step size remain constant over all depletion intervals.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

Calculation of the thermal conductivity of nanostructured Bi{sub 2}Te{sub 3} with the real phonon spectrum taken into account

The results of calculations of lattice thermal conductivity in bismuth telluride are presented and its decrease in a nanostructured material due to boundary scattering is estimated. Calculation is carried out using the lattice-dynamics method with the real phonon spectrum and the phonon–phonon interaction taken into account. Its results well agree with experimental data for the crystalline material. The estimates for the nanostructured material give a decrease in the thermal conductivity by 30% for grain sizes of 20 nm. Calculations by this method are compared with those obtained using approximate methods of description of the spectrum and scattering processes. It is shown that the least difference in the estimates (about 10%) can be obtained when using the approximation of a constant matrix element of phonon–phonon scattering with correction of the frequency dependence for acoustic modes.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Spectral partitioning with blends of eigenvectors

We report many common methods for data analysis rely on linear algebra. We provide new results connecting data analysis error to numerical accuracy in the context of spectral graph partitioning. We provide pointwise convergence guarantees so that spectral blends (linear combinations of eigenvectors) can be employed to solve data analysis problems with confidence in their accuracy. We apply this theory to an accessible model problem, the ring of cliques, by deriving the relevant eigenpairs and finding necessary and sufficient solver tolerances. Analysis of the ring of cliques provides an upper bound on eigensolver tolerances for graph partitioning problems. These results bridge the gap between linear algebra based data analysis methods and the convergence theory of iterative approximation methods. These results explain how the combinatorial structure of a problem can be recovered much faster than numerically accurate solutions to the associated linear algebra problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Stochastic DC optimal power flow with reserve saturation

We propose an optimization framework for stochastic optimal power flow with uncertain loads and renewable generator capacity. Our model follows previous work in assuming that generator outputs respond to load imbalances according to an affine control policy, but introduces a model of saturation of generator reserves by assuming that when a generator’s target level hits its limit, it abandons the affine policy and produces at that limit. This is a particularly interesting feature in models where wind power plants, which have uncertain upper generation limits, are scheduled to provide reserves to balance load fluctuations. Here, the resulting model is a nonsmooth nonconvex two-stage stochastic program, and we use a stochastic approximation method to find stationary solutions to a smooth approximation. Computational results on 6-bus and 118-bus test instances demonstrate that by considering the effects of saturation, our model can yield solutions with lower expected generation costs (at the same target line violation probability level) than those obtained from a model that enforces the affine policy to stay within generator limits with high probability.

17 WIND ENERGY↗
DOE OSTI2025

Exploring Data Set Bias and Decision Support with Predictive Uncertainty Through Bayesian Approximations and Convolutional Neural Networks

Individual seismic catalogs can contain multiscale observations from fault level to global scales and associated waveforms from discrete events reflect crustal structure across many different scales and locations. Seismic network aperture, geographic location, and observation distance may not provide informative guidance or intuition on how different catalogs will behave across models trained under different conditions. We rely on uncertainty to provide guardrails for when to trust model decisions, but understanding when our uncertainty is trustworthy is an open challenge. Here, in this work, we explore Bayesian approximation methods for assigning predictive uncertainty in seismic event classification problems. We find that computationally expensive Bayesian approximations do not outperform simple ensemble methods. We also find that when exploiting multiple seismic event catalogs, joint training with data from all the catalogs combined with Bayesian approximations and supervised training for classification can obscure bias and result in less robust uncertainty while also not providing substantial performance benefits compared to training individual models for each catalog.

58 GEOSCIENCES↗
DOE OSTI2016

Educating local radial basis functions using the highest gradient of interest in three dimensional geometries: EDUCATING LOCALIZED RBF-SCHEMES

We present a novel methodology to effectively localize radial basis function approximation methods in three dimensions. The local scheme requires shape parameter-dependent functions that can be used to approximate gradients of scattered data and to solve partial differential operators. The optimum shape parameter is obtained from the highest gradient of interest, where a known analytical function, when boundary conditions are not present, or a shape parameter-free global approximation are used to educate the localized scheme. The later option is applicable to problems where the operator needs to be solved multiple times, like in time evolution or stochastic integration. Past shape parameter's optimizations, for two-dimensional domains, based on the condition number of the interpolant matrix, were unable to provide satisfactory approximations. The applicability of our method is illustrated in the context of an analytical expression interpolation and during a Ginzburg-Landau relaxation of a free energy functional. In general, the optimum shape parameter depends on geometry, node distribution, and density, whereas the approximation errors decrease as the node density and the local stencil size increase. The effective localization of radial basis functions motivates its use in moving boundary problems and accelerates solutions through sparse matrix solvers.

free energy functional↗
DOE OSTI2021

Analytic insights into nonlocal energy transport. III. steady state Fokker Planck theory in spherical and planar geometry

This paper develops an approximate steady state Fokker Planck theory for energetic electron transport in a spherical laser fusion target. First, we apply the theory to only a small population of electrons, which is specified at the outset. Then, one determines the nonlocal electron energy flux of these electrons at a particular position. These energetic electrons may either come from the tail of the thermal distribution function or be generated by an instability. This paper develops two approximate methods of solution, which we call the “characteristic method” and “sparse eigenfunction.” The former works only in planar geometry and the latter in both planar and spherical geometry. Comparison of the two methods in planar geometry shows that even though the approximations are very different, they give about the same result, increasing their credibility. For the example we have chosen, spherical effects are not important for electrons from the tail of the distribution function but may well be for instability generated electrons, which have much higher energy. Comparing planar to spherical, one finds an additional spherical barrier protecting the fuel. It turns out that the associated fuel preheat a Fokker Planck model predicts is considerably less than that predicted by the Krook models as developed at both NRL and other places.

Manheimer, Wallace (ORCID:0000000163342591)↗
DOE OSTI2020

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

A posteriori error analysis of fluid–structure interactions: Time dependent error

A posteriori error analysis is a technique to quantify the error in particular simulations of a numerical approximation method. In this article, we use such an approach to analyze how various error components propagate in certain moving boundary problems. We study quasi-steady state simulations where slowly moving boundaries remain in mechanical equilibrium with a surrounding fluid. Such problems can be numerically approximated with the Method of Regularized Stokelets (MRS), a popular method used for studying viscous fluid–structure interactions, especially in biological applications. Our approach to monitoring the regularization error of the MRS is novel, along with the derivation of linearized adjoint equations to the governing equations of the MRS with a elastic elements. Our main numerical results provide a clear illustration of how the error evolves over time in several MRS simulations.

97 MATHEMATICS AND COMPUTING↗