Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,621 records · Page 90
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Singular perturbations in the state regulator problem

Most of the results of singular perturbation theory have been concerned with initial value problems whereas optimal control problems are of two-point boundary value type. The portions of this theory applicable to the open loop state regulator problem are reviewed. For obtaining approximate solutions to the state regulator problem the method of matched asymptotic expansions is employed. This method has been developed in connection with certain fluid mechanics problems and is applicable to nonlinear as well as linear problems. It has been found in the past to be advantageous not to formulate this method generally but to apply it to each individual problem and this approach is adopted here. A general recipe for the method is given and its application is illustrated by using the method to obtain an approximate solution to a simple, specific state regulator problem.

Ardema, M. D.↗
NASA NTRS1994

Relativistic and Correlation Effects in CuH, AgH and AuH: Comparison of Various Relativistic Methods

The effects of relativity on the bond lengths, dissociation energies, and harmonic vibrational frequencies of the 1Epsilon(+) electronic ground states of the group IB hydrides CuH, AgH and AuH have been evaluated with a variety of ab initio methods. These properties were investigated with moderately-sized basis sets at the self-consistent field Hartree Fock (SCF HF) level and with second-order Moller-Plesset (MP2) perturbation theory for electron correlation. Comparisons were made between all-electron results using the nonrelativistic Hamiltonian, perturbation theory (PT) at first-order with only the one-electron non-fine structure terms of the Breit-Pauli Hamiltonian, the spin-free Douglas-Kroll (DK) transformed Dirac Hamiltonian and the untransformed Dirac Hamiltonian, and results using two sets of relativistic effective core potentials (RECPs). The expected trends of bond length decrease, dissociation energy increase and harmonic frequency increase with both relativity and correlation are found. Both sets of RECPs are shown to give good results, if accompanied by a reasonable basis set. The DK method is demonstrated to be an inexpensive, reliable approximation to the DHF method.

Collins, Charlene L.↗
DOE OSTI2018

Multimode Jahn-Teller effect in bulk systems: A case of the N V 0 center in diamond

Here, the multimode Jahn-Teller (JT) effect in a bulk system of a neutral nitrogen-vacancy (NV 0 ) center in diamond is investigated via first-principles density-functional-theory calculations and the intrinsic distortion path (IDP) method. The adiabatic potential energy surface of the electronic ground state of the NV 0 center is calculated based on the local spin-density approximation. Our calculations confirm the presence of the dynamic Jahn-Teller effect in the ground 2 E state of the NV 0 center. Within the harmonic approximation, the IDP method provides the reactive path of JT distortion from unstable high-symmetry geometry to stable low-symmetry energy minimum geometry, and it describes the active normal modes participating in the distortion. We find that there is more than one vibrational mode contributing to the distortion, and their contributions change along the IDP. Several vibrational modes with large contributions to JT distortion, especially those modes close to 44 meV, are clearly observed as the phonon sideband in photoluminescence spectra in a series of experiments, indicating that the dynamic Jahn-Teller effect plays an important role in the optical transition of the NV 0 center.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1996

A Comparison of Neural Networks and Fuzzy Logic Methods for Process Modeling

The goal of this work was to analyze the potential of neural networks and fuzzy logic methods to develop approximate response surfaces as process modeling, that is for mapping of input into output. Structural response was chosen as an example. Each of the many methods surveyed are explained and the results are presented. Future research directions are also discussed.

Cios, Krzysztof J.↗
DOE OSTI2023

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

Logarithmic divergences in the k-inflationary power spectra computed through the uniform approximation

We investigate a calculation method for solving the Mukhanov-Sasaki equation in slow-roll k-inflation based on the uniform approximation (UA) in conjunction with an expansion scheme for slow-roll parameters with respect to the number of e-folds about the so-called turning point. Earlier works on this method have so far gained some promising results derived from the approximating expressions for the power spectra among others, up to second order with respect to the Hubble and sound flow parameters, when compared to other semi-analytical approaches (e.g., Green's function and WKB methods). However, a closer inspection is suggestive that there is a problem when higher-order parts of the power spectra are considered; residual logarithmic divergences may come out that can render the prediction physically inconsistent. Looking at this possibility, we map out up to what order with respect to the mentioned parameters several physical quantities can be calculated before hitting a logarithmically divergent result. It turns out that the power spectra are limited up to second order, the tensor-to-scalar ratio up to third order, and the spectral indices and running converge to all orders. This indicates that the expansion scheme is incompatible with the working equations derived from UA for the power spectra but compatible with that of the spectral indices. For those quantities that involve logarithmically divergent terms in the higher-order parts, existing results in the literature for the convergent lower-order parts calculated in the equivalent fashion should be viewed with some caution; they do not rest on solid mathematical ground.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Calculation of the spectrum of quasiparticle electron excitations in organic molecular semiconductors

A quasiparticle electronic spectrum belongs to the characteristics of nanoobjects that are most important for applications. The following methods of calculating the electronic spectrum are analyzed: the Kohn-Sham equations of the density functional theory (DFT), the hybrid functional method, the GW approximation, and the Lehmann approximation used in the spectral representation of one-electron Green’s function. The results of these approaches are compared with the data of photoemission measurements of benzene, PTCDA, and phthalocyanine (CuPc, H{sub 2}Pc, FePc, PtPc) molecules, which are typical representatives of organic molecular semiconductors (OMS). This comparison demonstrates that the Kohn-Sham equations of DFT incorrectly reproduce the electronic spectrum of OMS. The hybrid functional method correctly describes the spectrum of the valence and conduction bands; however, the HOMO-LUMO gap width is significantly underestimated. The correct gap width is obtained in both the GW approximation and the Lehmann approach, and the total energy in this approach can be calculated in the local density approximation of DFT.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Eigenvalue and eigenvector sensitivity and approximate analysis for repeated eigenvalue problems

A set of computationally efficient equations for eigenvalue and eigenvector sensitivity analysis are derived, and a method for eigenvalue and eigenvector approximate analysis in the presence of repeated eigenvalues is presented. The method developed for approximate analysis involves a reparamaterization of the multivariable structural eigenvalue problem in terms of a single positive-valued parameter. The resulting equations yield first-order approximations of changes in both the eigenvalues and eigenvectors associated with the repeated eigenvalue problem. Examples are given to demonstrate the application of such equations for sensitivity and approximate analysis.

Hou, Gene J. W.↗
DOE OSTI2018

Neutrino Oscillations in Supernovae

A (core-collapse) supernova is the final death call of a massive star during which a gigantic amount of a special kind of elementary particles known as neutrinos are released for a brief period of time. An earlier work by the PI and his collaborators has shown that the dense neutrino gas deep inside the supernova can experience "flavor transformation" or "oscillations" in a collective fashion, a macroscopic quantum phenomenon spanning over distances of 100-1000 kilometers. During this process neutrinos of different energies and propagation directions change their identities simultaneously. This phenomenon may be a key to understanding the mechanism of supernova explosion and the origin of many of the chemical elements found on the earth which were produced in the supernovae exploded in the past. The previous understanding of neutrino oscillations in supernovae was based on a stationary 3-dimensional model with 1 spatial dimension and 2 momentum dimensions. In this "neutrino bulb model" neutrinos are emitted from a spherically symmetric surface and they are decoupled from the motion of matter and the production of chemical elements. In order to accurately model the neutrino oscillations in supernovae and determine its physical consequences, sophisticated 7-dimensional (1 temporal dimension + 3 spatial dimensions + 3 momentum dimensions) modelling of neutrino transport must be performed self-consistently together with the hydrodynamic simulation of the flow of energy and matter inside the supernova and with the calculation of the production and destruction of chemical elements through a large nucleosynthesis network. Although it will be long journey to achieve the lofty goal, our project has made important progress in three different fronts towards this goal. In the first front we explored new computing hardware on which the next generation simulation codes may be run. We ported the neutrino bulb code to the platform of Graphics Processing Unit (GPU) and then later developed a new bulb code that can be run on supercomputers equipped with the Intel Xeon CPUs and/or the Xeon Phi coprocessors. In the second front we developed new algorithms that can be used towards computing neutrino oscillations in multi-dimensional models. We have developed a moment expansion method which is approximately 10 times more efficient than the old algorithm. We have also developed numerical codes to solve neutrino oscillations in toy models with two spatial dimensions. In the last front we used analytic methods to study neutrino oscillations in more sophisticated models, e.g. of multiple spatial dimensions and/or with time dependence. By analyzing the neutrino oscillations in toy models with two spacial dimensions we found that the approximate spacial translation symmetry, which is similar to the spherical symmetry in the neutrino bulb model, can be broken spontaneously as neutrino oscillations begin. By analyzing a time-dependent neutrino bulb model we found that certain modes of neutrino oscillations that are suppressed in the stationary model can become active in the time-dependent model. We have also developed theories which give insights to how ordinary matter may affect neutrino oscillations in supernovae and other astrophysical environments.

Duan, Huaiyu↗
DOE OSTI2015

Observations on variational and projector Monte Carlo methods

Variational Monte Carlo and various projector Monte Carlo (PMC) methods are presented in a unified manner. Similarities and differences between the methods and choices made in designing the methods are discussed. Both methods where the Monte Carlo walk is performed in a discrete space and methods where it is performed in a continuous space are considered. It is pointed out that the usual prescription for importance sampling may not be advantageous depending on the particular quantum Monte Carlo method used and the observables of interest, so alternate prescriptions are presented. The nature of the sign problem is discussed for various versions of PMC methods. A prescription for an exact PMC method in real space, i.e., a method that does not make a fixed-node or similar approximation and does not have a finite basis error, is presented. This method is likely to be practical for systems with a small number of electrons. Approximate PMC methods that are applicable to larger systems and go beyond the fixed-node approximation are also discussed.

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

Designing Cladded Components for High Temperature Nuclear Service: Part I—Analysis Methods

Use of corrosion-resistant cladding can greatly extend the design life of structural components in many advanced reactor systems. However, there are currently no ASME design rules for cladded components to guard against creep-fatigue failure and ratcheting strain accumulation in elevated temperature nuclear service. This paper, presented in two parts, addresses this gap by proposing a design strategy for cladded components that do not require long-term testing of clad materials. The proposed approach relies on approximate design analysis methods for two types of clad materials—a soft clad that creeps much faster than and has lower yield stress than the class A base material and a hard clad that creeps much slower than and has higher yield stress than the class A base material. Part I discusses the approximate analysis strategies for the clad materials—treat a soft clad as perfectly compliant and a hard clad as linear elastic—and Part II develops a complete set of design rules for each of the two types of cladded components. In conclusion, finite element analyses of representative high temperature reactor components show that the proposed design analysis methods can bound the design quantities in soft cladded components and approximate the design quantities in hard cladded components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2008

Modeling State-Space Aeroelastic Systems Using a Simple Matrix Polynomial Approach for the Unsteady Aerodynamics

A simple matrix polynomial approach is introduced for approximating unsteady aerodynamics in the s-plane and ultimately, after combining matrix polynomial coefficients with matrices defining the structure, a matrix polynomial of the flutter equations of motion (EOM) is formed. A technique of recasting the matrix-polynomial form of the flutter EOM into a first order form is also presented that can be used to determine the eigenvalues near the origin and everywhere on the complex plane. An aeroservoelastic (ASE) EOM have been generalized to include the gust terms on the right-hand side. The reasons for developing the new matrix polynomial approach are also presented, which are the following: first, the "workhorse" methods such as the NASTRAN flutter analysis lack the capability to consistently find roots near the origin, along the real axis or accurately find roots farther away from the imaginary axis of the complex plane; and, second, the existing s-plane methods, such as the Roger s s-plane approximation method as implemented in ISAC, do not always give suitable fits of some tabular data of the unsteady aerodynamics. A method available in MATLAB is introduced that will accurately fit generalized aerodynamic force (GAF) coefficients in a tabular data form into the coefficients of a matrix polynomial form. The root-locus results from the NASTRAN pknl flutter analysis, the ISAC-Roger's s-plane method and the present matrix polynomial method are presented and compared for accuracy and for the number and locations of roots.

Pototzky, Anthony S.↗
NASA NTRS1986

Implicit multigrid algorithms for the three-dimensional flux split Euler equations

The full approximation scheme multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computations required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. Results comparing pressure distributions with experimental data using both splitting types are shown. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined. Using the multigrid method on both subsonic and transonic wing calculations, the final lift coefficient is obtained to within 0.1 percent of its final value in a few as 15 cycles for a mesh with over 210,000 points. A spectral radius of 0.89 is achieved for both subsonic and transonic flow over the ONERA M6 wing while a spectral radius of 0.83 is obtained for supersonic flow over an analytically defined forebody. Results compared with experiment for all cases show good agreement.

Anderson, W. K.↗
DOE OSTI2020

Hole–hole Tamm–Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation

The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole–hole Tamm–Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole–hole counterpart to the more established particle–particle TDA (pp-TDA) method, both of which are derived from the particle–particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. Furthermore, the treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. Additionally, we show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states—particularly those with both low-lying ππ* and nπ* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.

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

Numerical simulation of scattering of acoustic waves by inelastic bodies using hypersingular boundary integral equation

A numerical method for solving a problem of diffraction of acoustic waves by system of solid and thin objects based on the reduction the problem to a boundary integral equation in which the integral is understood in the sense of finite Hadamard value is proposed. To solve this equation we applied piecewise constant approximations and collocation methods numerical scheme. The difference between the constructed scheme and earlier known is in obtaining approximate analytical expressions to appearing system of linear equations coefficients by separating the main part of the kernel integral operator. The proposed numerical scheme is tested on the solution of the model problem of diffraction of an acoustic wave by inelastic sphere.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Flux tube spectra from approximate integrability at low energies

We provide a detailed introduction to a method we recently proposed for calculating the spectrum of excitations of effective strings such as QCD flux tubes. The method relies on the approximate integrability of the low-energy effective theory describing the flux tube excitations and is based on the thermodynamic Bethe ansatz. The approximate integrability is a consequence of the Lorentz symmetry of QCD. For excited states, the convergence of the thermodynamic Bethe ansatz technique is significantly better than that of the traditional perturbative approach. We apply the new technique to the lattice spectra for fundamental flux tubes in gluodynamics in D = 3 + 1 and D = 2 + 1, and to k-strings in gluodynamics in D = 2 + 1. We identify a massive pseudoscalar resonance on the worldsheet of the confining strings in SU(3) gluodynamics in D = 3 + 1, and massive scalar resonances on the worldsheet of k = 2.3 strings in SU(6) gluodynamics in D = 2 + 1.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Analysis of Electromagnetic Scattering From a Dihedral Cut Into a Sphere Using Vector Spherical Harmonics

Proper calibration or verification of the calibration of fully - polarimetric radar requires more than a simple specular target, such as a conducting sphere or corner reflector. The dihedral reflector is a useful radar - scattering target for calibrating the polarization response of a radar system, because it possesses the ability to change the direction of the scattered electric - field vector in a special way, while a sphere or trihedral target does not. However, the scattered electromagnetic field of a sphere can be computed with arbitrary precision using the well - known Mie series, which is obtained by expanding the scattered field as a series of vector spherical harmonics. In contrast, the computation of the scattered electromagnetic field of dihedral reflectors typically relies on numerical methods or various approximate solutions of Maxwell's equations. Precision calibration of wide - bandwidth polarimetric radar has not been adequately achieved for some applications using this approach. This report describes a method for computing the scattered electromagnetic field of a specially shaped dihedral, derived through the application of a series expansion in vector spherical harmonics, with the goal of achieving the same precision that is available for the sphere.

47 OTHER INSTRUMENTATION↗
NASA NTRS1984

LANDSAT-D investigations in snow hydrology

Two stream methods provide rapid approximate calculations of radiative transfer in scattering and absorbing media. Although they provide information on fluxes only, and not on intensities, their speed makes them attractive to more precise methods. The methods provide a comprehensive, unified review for a homogeneous layer, and solve the equations for reflectance and transmittance for a homogeneous layer over a non reflecting surface. Any of the basic kernels for a single layer can be extended to a vertically inhomogeneous medium over a surface whose reflectance properties vary with illumination angle, as long as the medium can be subdivided into homogeneous layers.

Dozier, J.↗