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

How Useful are Conventional I–Vs for Performance Calibration of Single- and Two-Junction Perovskite Solar Cells? A Statistical Analysis of Performance Data on ≈200 Cells from 30 Global Sources

As perovskite photovoltaics (PV) advance from the laboratory to commercial prototypes, their accurate and reliable performance testing is becoming increasingly important. The well-documented dynamic response of perovskite solar cells to an external applied voltage has led to the development of steady-state performance measurement methods; however, these methods have not been widely adopted by the perovskite PV community. A key reason for this is that steady-state measurement methods take tens of minutes to complete, as opposed to conventional "fast" current-voltage (I-V) measurements usually lasting a few seconds. Fast I-Vs arise from a snapshot, almost always not a steady-state condition of the device; however, given their widespread use, the question arises: how do performance parameters of perovskite PV compare when measured with fast I-V and with a steady-state method? Results compiled from approximately 200 perovskite PV cells, including single junction, and two-terminal perovskite-perovskite and perovskite-Si tandems, show that fast I-Vs can provide a useful measure of the open-circuit voltage of the devices, while the short-circuit current and the overall efficiency can be widely misestimated. Here, the implications of these findings on performance testing protocols are discussed and possible options for fast and accurate testing of perovskite PV are proposed.

14 SOLAR ENERGY↗
DOE OSTI2018

Systems of Semilinear Parabolic Variational Inequalities with Time-Dependent Convex Obstacles

We consider a system of seminlinear parabolic variational inequalities with time-dependent convex obstacles. We prove the existence and uniqueness of its solution. We also provide a stochastic representation of the solution and show that it can be approximated by the penalization method. Our proofs are based upon probabilistic methods from the theory of Markov processes and the theory of backward stochastic differential equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Ab Initio Study of Electronic and Magnetic Properties in ZnO-Doped and Co-doped by Vanadium and Silver

Based upon the ab initio of spin density functional calculation, the electronic and magnetic properties of Zn{sub 1−x}M{sub x}O (M = Ag or V) and Zn{sub 1−x−y}A{sub x}B{sub y}O (A = V; B = Ag) are investigated, using the Korringa-Kohn-Rostoker (KKR) method coupled with the coherent potential approximation (CPA). The total and partial densities of state are calculated. The effect of concentration values in Zn{sub 1−x}V{sub x}O and Zn{sub 1−x−y}V{sub x}Ag{sub y}O are deduced. Moreover, the magnetic disorder local moment (DLM) and the total energy are obtained for different concentration values of doped and co-doped zinc oxide (ZnO).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Strong Correlation Effect in Ferrimagnetic Half-Metallic V{sub 2}CoAl and V{sub 2}CoGa Heusler Compounds

We present first-principles calculations based on the full-potential linearized augmented plane wave (FP-LAPW) method to investigate the structural, elastic, electronic, and magnetic properties of the Heusler compounds V{sub 2}CoAl and V{sub 2}CoGa. We explore different magnetic configurations in both regular and inverse Heusler structures: the Cu{sub 2}MnAl-type and Hg{sub 2}CuTi-type structures respectively. Our compounds are found to be ferrimagnets in the inverse Heusler structure. To take into account the effects of the strong correlation among the localized d-states, we perform a calculation using the GGA+U method (the general gradient approximation with the Hubbard on-site Coulomb correction U). Then we show that both compounds exhibit a half-metallic character and have a magnetic moment of 2.00 μ{sub B} per formula unit. This latter statement is coherent with the generalized Slater-Pauling rule (SP). In terms of the band structure calculations, we explain the formation of the indirect energy band gaps in the spin-majority channels. We also use the mean-field approximation (MFA) to discuss the magnetic exchange couplings and predict the Curie temperature for both compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Effects of Ag Doping and Point Defect on the Magnetism of ZnO

The magnetism sources and magnetic mechanism of Ag doping and point defect, which coexist in the presence of ZnO, are frequently controversial. To solve these problems, the effects of Ag doping and Zn or O vacancy on the magnetism of ZnO were investigated using geometry optimization and energy calculation according to the first-principles generalized gradient approximation + U (GGA+U) method of density functional theory. Results revealed that the system which Ag doping and Zn vacancy coexists in ZnO could achieve room-temperature ferromagnetism. This system had an extremely high spin polarization, which was advantageous for the spin electron injection sources applied in dilute magnetic semiconductors. Result also revealed that the magnetism of Zn{sub 14}AgO{sub 16} was caused by the electron hybrid coupling effects among the O-2p, Ag-4d, Ag-5s, and Zn-4s orbits. And these effects were produced by the hole carriers after complexes were formed by the Ag doping and Zn vacancy. With regard to the most stable structure of the ground state, all doping systems which Ag replacing Zn and O vacancy coexisted in ZnO and those which Ag replacing Zn, interstitial Ag, and O vacancy coexisted in ZnO were all non-magnetic. Thus, these doping systems were unsuitable for dilute magnetic semiconductors. By contrast, the doping systems of interstitial Ag and Zn vacancy were magnetic, although the magnetism was relatively weak, thus rendering such systems inapplicable as well.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Solution of the Problem of Free Vibrations of a Nonthin Orthotropic Shallow Shell of Variable Thickness in the Refined Statement

We consider the problem of investigation of the spectrum of natural vibrations of a nonthin orthotropic shallow shell variable in two coordinate directions of thickness in the nonclassical statement. The approach to the solution of the obtained two-dimensional boundary-value problem is based on its reduction (by the method of spline-approximation of the unknown functions along one coordinate direction) to the one-dimensional problem with its subsequent solution. We study different cases of boundary conditions imposed on the contours of the shell. We also perform the comparison and analysis of the natural frequencies and modes of vibrations of orthotropic shells of constant and variable thickness.

42 ENGINEERING↗
DOE OSTI2020

Singular value decomposition of adjoint flux distributions for Monte Carlo variance reduction

Monte Carlo (MC) shielding calculations often use weight windows (WWs) and biased sources formed from a deterministic estimate of the adjoint flux to improve the convergence rate of tallies. This requires a significant amount of computer memory, which can limit the memory available for high-resolution tally output. A new method is proposed for reducing these memory requirements by using singular value decomposition (SVD) in linear or logarithmic space to approximate the adjoint flux. This method’s performance is evaluated using the Shift and Denovo codes for streaming and diffusion base case problems, followed by problems using the Westinghouse AP1000 and the Joint European Torus. The log SVD reduced WW memory requirements by an order of magnitude in all cases without a significant performance penalty. Additionally, the linear SVD reduced biased source memory requirements by an order of magnitude, but further investigation is needed to account for observed limitations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

A hybrid approach for quantizing complicated motion of a charged particle in time-varying magnetic field

Quantum characteristics of a charged particle subjected to a singular oscillator potential under an external magnetic field is investigated via SU(1,1) Lie algebraic approach together with the invariant operator and the unitary transformation methods. The system we managed is somewhat complicated since we considered not only the time-variation of the effective mass of the system but also the dependence of the external magnetic field on time in an arbitrary fashion. In this case, the system is a kind of time-dependent Hamiltonian systems which require more delicate treatment when we study it. The complete wave functions are obtained without relying on the methods of perturbation and/or approximation, and the global phases of the system are identified. To promote the understanding of our development, we applied it to a particular case, assuming that the effective mass slowly varies with time under a time-dependent magnetic field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

astroABC : An Approximate Bayesian Computation Sequential Monte Carlo sampler for cosmological parameter estimation

Given the complexity of modern cosmological parameter inference where we arefaced with non-Gaussian data and noise, correlated systematics and multi-probecorrelated data sets, the Approximate Bayesian Computation (ABC) method is apromising alternative to traditional Markov Chain Monte Carlo approaches in thecase where the Likelihood is intractable or unknown. The ABC method is called"Likelihood free" as it avoids explicit evaluation of the Likelihood by using aforward model simulation of the data which can include systematics. Weintroduce astroABC, an open source ABC Sequential Monte Carlo (SMC) sampler forparameter estimation. A key challenge in astrophysics is the efficient use oflarge multi-probe datasets to constrain high dimensional, possibly correlatedparameter spaces. With this in mind astroABC allows for massive parallelizationusing MPI, a framework that handles spawning of jobs across multiple nodes. Akey new feature of astroABC is the ability to create MPI groups with differentcommunicators, one for the sampler and several others for the forward modelsimulation, which speeds up sampling time considerably. For smaller jobs thePython multiprocessing option is also available. Other key features include: aSequential Monte Carlo sampler, a method for iteratively adapting tolerancelevels, local covariance estimate using scikit-learn's KDTree, modules forspecifying optimal covariance matrix for a component-wise or multivariatenormal perturbation kernel, output and restart files are backed up everyiteration, user defined metric and simulation methods, a module for specifyingheterogeneous parameter priors including non-standard prior PDFs, a module forspecifying a constant, linear, log or exponential tolerance level,well-documented examples and sample scripts. This code is hosted online athttps://github.com/EliseJ/astroABC

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

A computationally efficient ductile damage model accounting for nucleation and micro-inertia at high triaxialities

The computational formulation of a micro-mechanical material model for the dynamic failure of ductile metals is presented in this paper. The statistical nature of porosity initiation is accounted for by introducing an arbitrary probability density function which describes the pores nucleation pressures. Each micropore within the representative volume element is modeled as a thick spherical shell made of plastically incompressible material. The treatment of porosity by a distribution of thick-walled spheres also allows for the inclusion of micro-inertia effects under conditions of shock and dynamic loading. The second order ordinary differential equation governing the microscopic porosity evolution is solved with a robust implicit procedure. A new Chebyshev collocation method is employed to approximate the porosity distribution and remapping is used to optimize memory usage. The adaptive approximation of the porosity distribution leads to a reduction of computational time and memory usage of up to two orders of magnitude. Moreover, the proposed model affords consistent performance: changing the nucleation pressure probability density function and/or the applied strain rate does not reduce accuracy or computational efficiency of the material model. The numerical performance of the model and algorithms presented is tested against three problems for high density tantalum: single void, one-dimensional uniaxial strain, and two-dimensional plate impact. Here, the results using the integration and algorithmic advances suggest a significant improvement in computational efficiency and accuracy over previous treatments for dynamic loading conditions.

42 ENGINEERING↗
DOE OSTI2021

Two-color pyrometry based flame to fuel surface radiative heat flux diagnostic using flamelets

In this study, a non-intrusive two-color pyrometry based radiative heat flux diagnostic is presented that accounts for the self-absorption effects. This diagnostic relies on looking up the flame temperature and soot volume fraction using the ratio of color intensities obtained from digital cameras and numerical solutions of 1D steady-state diffusion flames. Virtual two color pyrometry of computed flame show a unique one-to-one mapping of strained flames with intensity ratio; thereby allowing flames to be indexed and creation of a two-color flamelet manifold (TCFM). The TCFM is a function of two variables; the ratio of red to green intensities (flame index) and physical distance to the stoichiometric surface defined using a level-set function. The developed diagnostic is used in upward flame spread experiments where an approximate flame hull reconstruction method is applied to provide the 3D level-set function and intensity ratio for use of TCFM to provide radiative absorption properties for 3D ray tracing. Accounting for self- absorption effects provide improved soot volume fraction estimates. In conclusion, the radiative heat flux estimates are found to be in very good agreement with literature and improved results are obtained near the flame base.

42 ENGINEERING↗
DOE OSTI2017

Efficient block preconditioned eigensolvers for linear response time-dependent density functional theory

We present two efficient iterative algorithms for solving the linear response eigen- value problem arising from the time dependent density functional theory. Although the matrix to be diagonalized is nonsymmetric, it has a special structure that can be exploited to save both memory and floating point operations. In particular, the nonsymmetric eigenvalue problem can be transformed into a product eigenvalue problem that is self-adjoint with respect to a K-inner product. This product eigenvalue problem can be solved efficiently by a modified Davidson algorithm and a modified locally optimal block preconditioned conjugate gradient (LOBPCG) algorithm that make use of the K-inner product. The solution of the product eigenvalue problem yields one component of the eigenvector associated with the original eigenvalue problem. However, the other component of the eigenvector can be easily recovered in a postprocessing procedure. Therefore, the algorithms we present here are more efficient than existing algorithms that try to approximate both components of the eigenvectors simultaneously. The efficiency of the new algorithms is demonstrated by numerical examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Performance of multifunctional piezoelectric energy harvesting microgyroscopes with material degradation

Most microelectromechanical systems (MEMS) are made from brittle materials that are prone to cracking. Damage from the formation of cracks can occur anytime during the manufacturing process or due to mechanical loading of MEMS devices. Cracks can increase the structures flexibility and lead to changes in performance or total failure. Understanding the potential influences cracks in MEMS is necessary for the longevity of these systems. The energy harvesting microgyroscope studied in this work is susceptible to crack propagation in its silicon substrate and lead zirconate titanate piezoelectric patch. To this end, an electromechanical model is developed to capture the impacts of cracks on the performance of a damaged energy harvesting microgyroscope. Numerical methods are used to approximate the damage associated with distributed crack networks. Then, the static and dynamic responses of the damaged microgyroscope are analyzed for two possible multifunctional setups. Here, the current study revealed that cracks can significantly affect the static and dynamic behaviors of the microgyroscope energy harvester. Combinations of cracking and piezoelectric property degradation cause a decrease in the static pull-in voltage and shift the resonance frequencies of the system. Additionally, it is shown that the systems sensitivity to changes in the base rotation and energy harvesting capabilities vary widely depending on the severity of material degradation.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Single-ion hydration thermodynamics from clusters to bulk solutions: Recent insights from molecular modeling

The importance of single-ion hydration thermodynamic properties for understanding the driving forces of aqueous electrolyte processes, along with the impossibility of their direct experimental measurement, have prompted a large number of experimental, theoretical, and computational studies aimed at separating the cation and anion contributions. Here we provide an overview of historical approaches based on extrathermodynamic assumptions and more recent computational studies of single-ion hydration in order to evaluate the approximations involved in these methods, quantify their accuracy, reliability, and limitations in the light of the latest developments. Finally, we also offer new insights into the factors that influence the accuracy of ion–water interaction models and our views on possible ways to fill this substantial knowledge gap in aqueous physical chemistry.

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

What Is the Fractional Laplacian? A Comparative Review with New Results

The fractional Laplacian in $\mathbb{R}^{d}$, which we write as (–Δ) α/2 with, α ϵ (0, 2) has multiple equivalent characterizations. Furthermore, in bounded domains, boundary conditions must be incorporated in these characterizations in mathematically distinct ways, and there is currently no consensus in the literature as to which definition of the fractional Laplacian in bounded domains is most appropriate for a given application. The Riesz (or integral) definition, for example, admits a nonlocal boundary condition, where the value of a function must be prescribed on the entire exterior of the domain in order to compute its fractional Laplacian. Yet, the spectral definition requires only the standard local boundary condition. These differences, among others, lead us to ask the question: “What is the fractional Laplacian?” Beginning from first principles, we compare several commonly used definitions of the fractional Laplacian theoretically, through their stochastic interpretations as well as their analytical properties. Next, we present quantitative comparisons using a sample of state-of-the-art methods. We discuss recent advances on nonzero boundary conditions and present new methods to discretize such boundary value problems: radial basis function collocation (for the Riesz fractional Laplacian) and nonharmonic lifting (for the spectral fractional Laplacian). In our numerical studies, we aim to compare different definitions on bounded domains using a collection of benchmark problems. We consider the fractional Poisson equation with both zero and nonzero boundary conditions, where the fractional Laplacian is defined according to the Riesz definition, the spectral definition, the directional definition, and the horizon-based nonlocal definition. We verify the accuracy of the numerical methods used in the approximations for each operator, and we focus on identifying differences in the boundary behaviors of solutions to equations posed with these different definitions. Through our efforts, we aim to further engage the research community in open problems and assist practitioners in identifying the most appropriate definition and computational approach to use for their mathematical models in addressing anomalous transport in diverse applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Geostatistical interpolation of streambed hydrologic attributes with addition of left censored data and anisotropy

Spatial geostatistical interpolation of point measurements of streambed attributes in the hyporheic zone may be constrained by the streambed anisotropy, and data density and spatial distribution may significantly impact the results. Spatial clustering and low spatial data density can be caused by bedrock outcropping at the streambed limiting installation of in-stream piezometers. This study examines parameter error variability of the geostatistical interpolation using anisotropic interpolation methods and increasing the data density by adding left censored values (i.e., data below measurement limit) to locations where measurements were limited by exposed bedrock lining the streambed. The reduction in relative standard error of the interpolation was determined for the spatial distributions of streambed attributes including hydraulic conductivity, seepage flux, and mercury solute flux measured in two different years along a study reach in East Fork Poplar Creek, Tennessee, USA. Here, two methods to impute the left censored values were compared including the conventional half the detection limit substitution method, and the Stochastic Approximation of Expectation-Maximization (SAEM) algorithm, which both had comparable results. Imputing left censored data increased the data density to recommended ranges, reduced data clustering, increased the spatial dependence for some attributes, and reduced the standard error for each of the three attributes. For the reach considered herein, addition of the left censored values resulted in a larger error reduction than the consideration of anisotropy within the interpolation, which confirms the benefit of data addition to increase data density within data-limited river corridors.

58 GEOSCIENCES↗
DOE OSTI2015

Thermoelectric, electronic, optical and chemical bonding properties of Ba{sub 2}PrRuO{sub 6}: At temperature 7 K and 150 K

Highlights: • DFT-FPLAPW method used for calculating the electronic structure. • The Fermi surface of BPRO (7 K and 150 K) is also calculated. • The complex dielectric function has been calculated. • Thermoelectric properties were also calculated using BoltzTraP code. • Power factor shows that both compounds are good thermoelectric materials at 600 K. - Abstract: We present first principles calculations of the band structure, density of states, electronic charge density, Fermi surface and optical properties of Ba{sub 2}PrRuO{sub 6} single crystals at two different temperatures. The atomic positions were optimized by minimizing the forces acting on the atoms. We have employed the full potential linear augmented plane wave method within local density approximation, generalized gradient approximation and Engel–Vosko generalized gradient approximation to treat the exchange correlation potential. The calculation shows that the compound is superconductor with strong hybridization near the Fermi energy level. Fermi surface is composed of two sheets. The calculated electronic specific heat capacities indicate, very close agreement with the experimental one. The bonding features of the compounds are analyzed using the electronic charge density in the (1 0 0) and (0–10) crystallographic planes. The dispersion of the optical constants was calculated and discussed. The thermoelectric properties are also calculated using the BoltzTrap code.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗