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,765 records · Page 98
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Ab-initio study of electronic structure and elastic properties of ZrC

The electronic and elastic properties of ZrC have been investigated using the linear combination of atomic orbitals method within the framework of density functional theory. Different exchange-correlation functionals are taken into account within generalized gradient approximation. We have computed energy bands, density of states, elastic constants, bulk modulus, shear modulus, Young’s modulus, Poisson’s ratio, lattice parameters and pressure derivative of the bulk modulus by calculating ground state energy of the rock salt structure type ZrC.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Brownian motion of a plasma crystal

The dynamics of a plasma crystal under the action of random external forces is considered. The pair correlation functions of the particle displacements are calculated in the harmonic approximation by using the Langevin equations. The case of a planar hexagonal lattice is analyzed in more detail. Analogues of the Van Hove singularities in the spectral densities are discovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

An examination of the adiabatic approximation in cosmic ray propagation theory

A theory for the pitch-angle scattering of cosmic rays in a turbulent magnetic field can be derived from first principles by applying both the quasi-linear and adiabatic approximations to the master equation for the ensemble averaged distribution function. A proof of the failure of these approximations, taken together, is given. Some predictions of the quasi-linear theory, are given. New wave-like propagation modes have been discovered.

Scudder, J.↗
NASA NTRS1987

Io: Comparison of photometric scans produced by the Minnaert and Hapke functions

Experience has shown that the empirical Minnaert function is a very useful approximation to real photometric behavior near opposition (phase angle = 0), but that in general it cannot accurately model photometric scans across the face of even a homogeneous planet at higher phase angles. Given recent work on fitting the rigorous Hapke function to Voyager data for Io, we can test to what degree the Minnaert function breaks down in the case of Ionian materials by comparing photometric scans produced by the two approaches. At phase angles 2, 20, 30, 60, and 90 deg, we have computed scans of the reflectance along the photometric equator and mirror meridian that would be expected for a homogeneous planet whose surface obeys Hapke's law. The values used for the Hapke parameters omega(0), h, g, and bar-theta were derived for Io by Simonelli and Veverka in both the Voyager narrow angle camera violet filter and orange filter. Each calculated Hapke scan is compared with the corresponding scan predicted by Minnaert's law for various values of the Minnaert limb-darkening parameter k. For a Minnaert scan at a particular k, the value of the reflectance parameter B(0) is arbitrarily chosen so that the Minnaert and Hapke scans coincide at the so-called specular point (the point where the photometric equator and mirror meridian intersect). The violet filter photometric scans that result from this process are shown. Results indicate that the empirical Minnaert function, while not a perfect model of real photometric behavior, does provide a very useful parameterization of limb darkening at phase angles out to 90 deg, and is especially useful near opposition.

Simonelli, Damon P.↗
DOE OSTI2017

Computationally predicted energies and properties of defects in GaN

Recent developments in theoretical techniques have significantly improved the predictive power of density-functional-based calculations. In this article, we discuss how such advancements have enabled improved understanding of native point defects in GaN. We review the methodologies for the calculation of point defects, and discuss how techniques for overcoming the band-gap problem of density functional theory affect native defect calculations. In particular, we examine to what extent calculations performed with semilocal functionals (such as the generalized gradient approximation), combined with correction schemes, can produce accurate results. The properties of vacancy, interstitial, and antisite defects in GaN are described, as well as their interaction with common impurities. We also connect the first-principles results to experimental observations, and discuss how native defects and their complexes impact the performance of nitride devices. Overall, we find that lower-cost functionals, such as the generalized gradient approximation, combined with band-edge correction schemes can produce results that are qualitatively correct. However, important physics may be missed in some important cases, particularly for optical transitions and when carrier localization occurs.

36 MATERIALS SCIENCE↗
DOE OSTI2016

A wave based method to predict the absorption, reflection and transmission coefficient of two-dimensional rigid frame porous structures with periodic inclusions

This paper presents an extension to the Wave Based Method to predict the absorption, reflection and transmission coefficients of a porous material with an embedded periodic set of inclusions. The porous unit cell is described using the Multi-Level methodology and by embedding Bloch–Floquet periodicity conditions in the weighted residual scheme. The dynamic pressure field in the semi-infinite acoustic domains is approximated using a novel wave function set that fulfils the Helmholtz equation, the Bloch–Floquet periodicity conditions and the Sommerfeld radiation condition. The method is meshless and computationally efficient, which makes it well suited for optimisation studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Generalized plasma dispersion function: One-solve-all treatment, visualizations, and application to Landau damping

A unified, fast, and effective approach is developed for numerical calculation of the well-known plasma dispersion function with extensions from Maxwellian distribution to almost arbitrary distribution functions, such as the δ, flat top, triangular, κ or Lorentzian, slowing down, and incomplete Maxwellian distributions. The singularity and analytic continuation problems are also solved generally. Given that the usual conclusion γ∝∂f{sub 0}/∂v is only a rough approximation when discussing the distribution function effects on Landau damping, this approach provides a useful tool for rigorous calculations of the linear wave and instability properties of plasma for general distribution functions. The results are also verified via a linear initial value simulation approach. Intuitive visualizations of the generalized plasma dispersion function are also provided.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Excitation energies of molecules within time-independent density functional theory

Recently proposed exchange energy functional for excited-states is tested for obtaining excitation energies of diatomic molecules. The functional is the ground-state counterpart of the local-density approximation, the modified local spin density (MLSD). The MLSD functional is tested for the N{sub 2} and CO diatomic molecules. The excitation energy obtained with the MLSD functional for the N{sub 2} molecule is in close vicinity to that obtained from the exact exchange orbital functional, Krieger, Li and Iafrate (KLI). For the CO molecule, the departure in excitation energy is observed and is due to the overcorrection of self-interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Effect of mechanical strain on electronic properties of bulk MoS{sub 2}

Ab-initio density functional theory based calculations of electronic properties of bulk and monolayer Molybdenum di-Sulfide (MoS{sub 2}) have been performed using all electron Full Potential Linearised Augmentad Plane Wave (FPLAPW) method using Elk code. We have used Generalised Gradient Approximation (GGA) for exchange and correlation functionals and performed calculaitons of Lattice parameters, Density Of States (DOS) and Band Structure (BS). Band structure calculations revealed that bulk MoS{sub 2} has indirect band gap of 0.97 eV and mono-layer MoS{sub 2} has direct band gap which has increased to 1.71 eV. These are in better agreement with experimental values as compared with the other calculations using pseudo-potential code. The effect of mechanical strain on the electronic properties of bulk MoS{sub 2} has also been studied. For the different values of compressive strain (varying from 2% to 8% in steps of 2%) along the c-axis, the corresponding DOS and BS are obtained. We observed that the band gap decreases by about 15% for every 2% increase in strain along the c-axis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Electron scattering on A = 3 nuclei from quantum Monte Carlo based approaches

We perform first-principle calculations of electron-nucleus scattering on 3 He and 3 H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of 3 He, and the inclusive cross sections of both 3 He and 3 H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. Here, we also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Electron scattering on ${\mathbf{A=3}}$ nuclei from quantum Monte Carlo based approaches

We perform first-principle calculations of electron-nucleus scattering on $^3$He and $^3$H using the Green's function Monte Carlo method and two approaches based on the factorization of the final hadronic state: the spectral-function formalism and the short-time approximation. These three methods are benchmarked among each other and compared to the experimental data for the longitudinal and transverse electromagnetic response functions of $^3$He, and the inclusive cross sections of both $^3$He and $^3$H. Since these three approaches are based on the same description of nuclear dynamics of the initial target state, comparing their results enables a precise quantification of the uncertainties inherent to factorization schemes. At sufficiently large values of the momentum transfer, we find an excellent agreement of the Green's function Monte Carlo calculation with experimental data and with both the spectral-function formalism and the short-time approximation. We also analyze the relevance of relativistic effects, whose inclusion becomes crucial to explain data at high momentum and energy transfer.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1991

Test functions for the moment method which yield the minimum mean square error

The method of moments has found extensive applications in the solution of a wide class of electromagnetic problems. Additionally, there are a number of basis functions which have found favor and with them different types of testing methods. This work examines the selection of test functions, given a selected set of basis functions, which serve to minimize the mean square error. These test functions will guarantee equal or lesser error compared with any other test function. Additionally, these test functions provide a monotonic improvement of the approximation as the increase in the order of the system. Hence while these test functions require the application of the differential or integral operator one additional time, they provide an absolute lower bound to the mean square error in the approximation and a means to systematically improve accuracy with each increase in the order of the unknowns.

Shamansky, H.↗
NASA NTRS1987

Efficient probabilistic fracture mechanics analysis

A systematic and efficient method for probabilistic fracture mechanics analysis is proposed. The method is based on a most-probable-point-locus concept. The locus is obtained iteratively where the initial locus is determined using the linear approximation of the service life N(X) function about the mean values of the random variables X. Linear and quadratic approximations of N(X) are established locally at the most probable points, and the reliability analysis methods are used to compute the cumulative probabilities. By using two examples, the proposed method is demonstrated to be efficient and accurate. One example involved a random loading and N(X) was computed using cycle-by-cycle integration. The method is general and can be applied to other performance functions. It is particularly suitable when the computation of the performance function is time consuming such that Monte Carlo simulation is prohibitively costly.

Wu, Y.-T.↗
DOE OSTI2013

LIKELIHOOD-FREE COSMOLOGICAL INFERENCE WITH TYPE Ia SUPERNOVAE: APPROXIMATE BAYESIAN COMPUTATION FOR A COMPLETE TREATMENT OF UNCERTAINTY

Cosmological inference becomes increasingly difficult when complex data-generating processes cannot be modeled by simple probability distributions. With the ever-increasing size of data sets in cosmology, there is an increasing burden placed on adequate modeling; systematic errors in the model will dominate where previously these were swamped by statistical errors. For example, Gaussian distributions are an insufficient representation for errors in quantities like photometric redshifts. Likewise, it can be difficult to quantify analytically the distribution of errors that are introduced in complex fitting codes. Without a simple form for these distributions, it becomes difficult to accurately construct a likelihood function for the data as a function of parameters of interest. Approximate Bayesian computation (ABC) provides a means of probing the posterior distribution when direct calculation of a sufficiently accurate likelihood is intractable. ABC allows one to bypass direct calculation of the likelihood but instead relies upon the ability to simulate the forward process that generated the data. These simulations can naturally incorporate priors placed on nuisance parameters, and hence these can be marginalized in a natural way. We present and discuss ABC methods in the context of supernova cosmology using data from the SDSS-II Supernova Survey. Assuming a flat cosmology and constant dark energy equation of state, we demonstrate that ABC can recover an accurate posterior distribution. Finally, we show that ABC can still produce an accurate posterior distribution when we contaminate the sample with Type IIP supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Gravitation- And Conduction-Driven Melting In A Sphere

Simplifying assumptions lead to approximate closed-form solution. Theoretical paper discusses melting of solid sphere in spherical container. Develops mathematical model of melting process, based in part on simplifying assumptions like those used in theories of lubrication and film condensation. Resulting equation for melting speed as function of melting distance solved approximately in closed form.

Bahrami, Parviz A.↗
NASA NTRS1991

A thermal/nonthermal approach to solar flares

An approach for modeling solar flare high-energy emissions is developed in which both thermal and nonthermal particles coexist and contribute to the radiation. The thermal/nonthermal distribution function is interpreted physically by postulating the existence of DC sheets in the flare region. The currents then provide both primary plasma heating through Joule dissipation, and runaway electron acceleration. The physics of runaway acceleration is discussed. Several methods are presented for obtaining approximations to the thermal/nonthermal distribution function, both within the current sheets and outside of them. Theoretical hard x ray spectra are calculated, allowing for thermal bremsstrahlung from the heated plasma electrons impinging on the chromosphere. A simple model for hard x ray images of two-ribbon flares is presented. Theoretical microwave gyrosynchrotron spectra are calculated and analyzed, uncovering important new effects caused by the interplay of thermal and nonthermal particles. The theoretical spectra are compared with observed high resolution spectra of solar flares, and excellent agreement is found, in both hard x rays and microwaves. The future detailed application of this approach to solar flares is discussed, as are possible refinements to this theory.

Benka, Stephen G.↗
NASA NTRS1992

Electromagnetic scattering by coated convex surfaces and wedges simulated by approximate boundary conditions

Asymptotic/high-frequency solutions are developed for analyzing the non-specular scattering mechanisms associated with coated convex surfaces and edges simulated by approximate boundary conditions. In particular, the standard impedance boundary conditions (SIBC's) and the second order generalized impedance boundary conditions (GIBC's) are employed for a characterization of the edge diffraction, creeping wave, and surface diffracted wave contributions. To study the creeping wave and surface diffracted wave mechanisms, rigorous UTD (uniform geometrical theory of diffraction) diffraction coefficients are developed for a convex coated cylinder simulated with SIBC's and GIBC's. The ray solutions obtained remain valid in the transition region and reduce uniformly to those in the deep lit and shadow regions. A uniform asymptotic solution is also presented for observations in the close vicinity of the cylinder. The diffraction coefficient for a convex cylinder are obtained via a generalization of the corresponding ones of the circular cylinder. To validate the asymptotic/high-frequency solution, integral equations are derived for both E and H-polarization and solved numerically using the method of moments. Results are presented for a single and three layered coated convex cylinder. Some insights are also provided on the accuracy of the employed GIBC's versus SIBC's for application to curved surfaces. To characterize the scattering by impedance wedges illuminated at skew incidence, diffraction coefficients are derived from an approximate solution of the governing functional difference equations. This solution exactly recovers the known ones for an impedance half plane or an arbitrary wedge at normal incidence, and to validate it for other wedge angles, a moment method code was used. Finally, to test the usefulness of the approximate skew incidence impedance wedge diffraction coefficient for three dimensional structures, equivalent currents are derived in the context of the physical theory of diffraction (PTD) for a finite length impedance wedge of arbitrary internal angles. These are incorporated in a standard general purpose PTD code and results are presented for a number of different impedance structures.

Syed, H. H.↗
DOE OSTI2018

Predicting the thermodynamic stability of double-perovskite halides from density functional theory

Recently, a series of double-perovskite halide compounds such as Cs 2 AgBiCl 6 and Cs 2 AgBiBr 6 have attracted intensive interest as promising alternatives to the solar absorber material CH 3 NH 3 PbI 3 because they are Pb-free and may exhibit enhanced stability. The thermodynamic stability of a number of double-perovskite halides has been predicted based on density functional theory (DFT) calculations of compound formation energies. In this paper, we found that the stability prediction can be dependent on the approximations used for the exchange-correlation functionals, e.g., the DFT calculations using the widely used Perdew, Burke, Ernzerhof (PBE) functional predict that Cs 2 AgBiBr 6 is thermodynamically unstable against phase-separation into the competing phases such as AgBr, Cs 2 AgBr 3 , Cs 3 Bi 2 Br 9 , etc., obviously inconsistent with the good stability observed experimentally. The incorrect prediction by the PBE calculation results from its failure to predict the correct ground-state structures of AgBr, AgCl, and CsCl. By contrast, the DFT calculations based on local density approximation, optB86b-vdW, and optB88-vdW functionals predict the ground-state structures of these binary halides correctly. Furthermore, the optB88-vdW functional is found to give the most accurate description of the lattice constants of the double-perovskite halides and their competing phases. Given these two aspects, we suggest that the optB88-vdW functional should be used for predicting thermodynamic stability in the future high-throughput computational material design or the construction of the Materials Genome database for new double-perovskite halides. As a result, using different exchange-correlation functionals has little influence on the dispersion of the conduction and the valence bands near the electronic bandgap; however, the calculated bandgap can be affected indirectly by the optimized lattice constant, which varies for different functionals.

36 MATERIALS SCIENCE↗