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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Ares I-X In-Flight Modal Identification

Operational modal analysis is a procedure that allows the extraction of modal parameters of a structure in its operating environment. It is based on the idealized premise that input to the structure is white noise. In some cases, when free decay responses are corrupted by unmeasured random disturbances, the response data can be processed into cross-correlation functions that approximate free decay responses. Modal parameters can be computed from these functions by time domain identification methods such as the Eigenvalue Realization Algorithm (ERA). The extracted modal parameters have the same characteristics as impulse response functions of the original system. Operational modal analysis is performed on Ares I-X in-flight data. Since the dynamic system is not stationary due to propellant mass loss, modal identification is only possible by analyzing the system as a series of linearized models over short periods of time via a sliding time-window of short time intervals. A time-domain zooming technique was also employed to enhance the modal parameter extraction. Results of this study demonstrate that free-decay time domain modal identification methods can be successfully employed for in-flight launch vehicle modal extraction.

Bartkowicz, Theodore J.↗
DOE OSTI2018

Collapse of the electron gas from three to two dimensions in Kohn-Sham density functional theory

Here under pressure, a quasi-2D electron gas can collapse toward the true 2D limit. In this limit, the exact exchange-correlation energy per electron has a known finite limit, but general-purpose semilocal approximate density functionals, like the local density approximation (LDA) and Perdew-Burke-Ernzerhof generalized gradient approximation (PBE GGA), are known to diverge to minus infinity. Here we consider a model density for a non-interacting electron gas confined to a thickness $L$ by infinite-barrier walls, with a fixed 2D density $1/(\pi (r_s^{\text{2D}})^2)$, and $r_s^{\text{2D}}=4$ Bohr. We estimate that LDA, PBE, and the strongly constrained and appropriately normed (SCAN) meta-GGA are accurate for the exchange-correlation energy over a wide quasi-2D range, $1.5< L/r_s^{\text{2D}} < 3.85$, but not for smaller $L$. Of these functionals, only SCAN tends to a finite limit when $L$ tends to 0. Since the non-interacting kinetic energy, treated exactly in Kohn-Sham theory, dominates in this limit within a deformable jellium model, all of the general-purpose functionals can estimate the pressure required to achieve any thickness (with SCAN and LDA better than PBE). This pressure vanishes around $L/r_s^{\text{2D}} = 3.85$, where the 3D electron density is roughly that of the valence electrons in metallic potassium, and reaches about 20 GPa at $L/r_s^{\text{2D}} = 1.5$ and 400 GPa at $L/r_s^{\text{2D}} = 0.6$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Exploiting Kronecker structure in exponential integrators: Fast approximation of the action of φ $-$functions of matrices via quadrature

Here, in this article, we propose an algorithm for approximating the action of φ $-$ functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the φ $-$ functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive a priori bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1976

On the use of Pade approximants to represent unsteady aerodynamic loads for arbitrarily small motions of wings

The general behavior of unsteady airloads in the frequency domain is explained. Based on this, a systematic procedure is described whereby the airloads, produced by completely arbitrary, small, time-dependent motions of a thin lifting surface in an airstream, can be predicted. This scheme employs as raw materials any of the unsteady linearized theories that have been mechanized for simple harmonic oscillations. Each desired aerodynamic transfer function is approximated by means of an appropriate Pade approximant, that is, a rational function of finite degree polynomials in the Laplace transform variable. Although these approximations have many uses, they are proving especially valuable in the design of automatic control systems intended to modify aeroelastic behavior.

Vepa, R.↗
DOE OSTI2019

A Semi-Lagrangian Spectral Method for the Vlasov–Poisson System Based on Fourier, Legendre and Hermite Polynomials

In this study, we apply a semi-Lagrangian spectral method for the Vlasov–Poisson system, previously designed for periodic Fourier discretizations, by implementing Legendre polynomials and Hermite functions in the approximation of the distribution function with respect to the velocity variable. We discuss second-order accurate-in-time schemes, obtained by coupling spectral techniques in the space–velocity domain with a BDF time-stepping scheme. The resulting method possesses good conservation properties, which have been assessed by a series of numerical tests conducted on some standard benchmark problems including the two-stream instability and the Landau damping test cases. In the Hermite case, we also investigate the numerical behavior in dependence of a scaling parameter in the Gaussian weight. Confirming previous results from the literature, our experiments for different representative values of this parameter, indicate that a proper choice may significantly impact on accuracy, thus suggesting that suitable strategies should be developed to automatically update the parameter during the time-advancing procedure.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Phosphorus allotropes: Stability of black versus red phosphorus re-examined by means of the van der Waals inclusive density functional method

We evaluate the energetic stabilities of white, red, and black allotropes of phosphorus using density functional theory (DFT) and hybrid functional methods, van der Waals (vdW) corrections (DFT+vdW and hybrid+vdW), vdW density functionals, and random phase approximation (RPA). In this work, we find that stability of black phosphorus over red-V (i.e., the violet form) is not ubiquitous among these methods, and the calculated enthalpies for the reaction phosphorus (red-V)→ phosphorus (black) are scattered between –20 and 40 meV/atom. With local density and generalized gradient approximations, and hybrid functionals, mean absolute errors (MAEs) in densities of P allotropes relative to experiments are found to be around 10%–25%, whereas with vdW-inclusive methods, MAEs in densities drop below ~5%. While the inconsistency among the density functional methods could not shed light on the stability puzzle of black versus red phosphorus, comparison of their accuracy in predicting densities and the supplementary RPA results on relative stabilities indicate that opposite to the common belief, black and red phosphorus are almost degenerate, or the red-V (violet) form of phosphorus might even be the ground state.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Quasiparticle and hybrid density functional methods in defect studies: An application to the nitrogen vacancy in GaN

Defects in semiconductors can play a vital role in the performance of electronic devices, with native defects often dominating the electronic properties of the semiconductor. Understanding the relationship between structural defects and electronic function will be central to the design of new high-performance materials. In particular, it is necessary to quantitatively understand the energy and lifetime of electronic states associated with the defect. In this work, we apply first-principles density functional theory (DFT) and many-body perturbation theory within the GW approximation to understand the nature and energy of the defect states associated with a charged nitrogen vacancy on the electronic properties of gallium nitride (GaN), as a model of a well-studied and important wide gap semiconductor grown with defects. We systematically investigate the sources of error associated with the GW approximation and the role of the underlying atomic structure on the predicted defect state energies. Additionally, analysis of the computed electronic density of states (DOS) reveals that there is one occupied defect state 0.2 eV below the valence band maximum and three unoccupied defect states at energy of 0.2–0.4 eV above the conduction band minimum, suggesting that this defect in the +1 charge state will not behave as a carrier trap. Furthermore, we compare the character and energy of the defect state obtained from GW and DFT using the HSE approximate density functional and find excellent agreement. This systematic study provides a more complete understanding of how to obtain quantitative defect energy states in bulk semiconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1975

Compression of ephemerides

An algorithm is proposed for generating sequences of Chebyshev series which are the best approximations of an astronomical ephemeris in the sense of Chebyshev over large intervals of time. The criterion for a polynomial approximation of a function to be the best polynomial approximation of the function is that the error function present certain rippling characteristics as described by Remez (1957). General features of the program in PL/1 are described.

Deprit, A.↗
NASA NTRS1984

Rational-spline approximation with automatic tension adjustment

An algorithm for weighted least-squares approximation with rational splines is presented. A rational spline is a cubic function containing a distinct tension parameter for each interval defined by two consecutive knots. For zero tension, the rational spline is identical to a cubic spline; for very large tension, the rational spline is a linear function. The approximation algorithm incorporates an algorithm which automatically adjusts the tension on each interval to fulfill a user-specified criterion. Finally, an example is presented comparing results of the rational spline with those of the cubic spline.

Schiess, J. R.↗
NASA NTRS1990

Linear and non-linear dynamic models of a geared rotor-bearing system

A three degree of freedom non-linear model of a geared rotor-bearing system with gear backlash and radial clearances in rolling element bearings is proposed here. This reduced order model can be used to describe the transverse-torsional motion of the system. It is justified by comparing the eigen solutions yielded by corresponding linear model with the finite element method results. Nature of nonlinearities in bearings is examined and two approximate nonlinear stiffness functions are proposed. These approximate bearing models are verified by comparing their frequency responses with the results given by the exact form of nonlinearity. The proposed nonlinear dynamic model of the geared rotor-bearing system can be used to investigate the dynamic behavior and chaos.

Kahraman, Ahmet↗
NASA NTRS2001

Prediction of Unsteady Aerodynamic Coefficients at High Angles of Attack

The nonlinear indicial response method is used to model the unsteady aerodynamic coefficients in the low speed longitudinal oscillatory wind tunnel test data of the 0.1 scale model of the F-16XL aircraft. Exponential functions are used to approximate the deficiency function in the indicial response. Using one set of oscillatory wind tunnel data and parameter identification method, the unknown parameters in the exponential functions are estimated. The genetic algorithm is used as a least square minimizing algorithm. The assumed model structures and parameter estimates are validated by comparing the predictions with other sets of available oscillatory wind tunnel test data.

Pamadi, Bandu N.↗
DOE OSTI2018

Convergence of Probability Densities Using Approximate Models for Forward and Inverse Problems in Uncertainty Quantification

Here, we analyze the convergence of probability density functions utilizing approximate models for both forward and inverse problems. We consider the standard forward uncertainty quantification problem where an assumed probability density on parameters is propagated through the approximate model to produce a probability density, often called a push-forward probability density, on a set of quantities of interest (QoI). The inverse problem considered in this paper seeks to update an initial probability density assumed on model input parameters such that the subsequent push-forward of this updated density through the parameter-to-QoI map matches a given probability density on the QoI. We prove that the densities obtained from solving the forward and inverse problems, using approximate models, converge to the true densities as the approximate models converge to the true models. Numerical results are presented to demonstrate convergence rates of densities for sparse grid approximations of parameter-to-QoI maps and standard spatial and temporal discretizations of PDEs and ODEs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1977

Elastic scattering of deuterons by carbon

Deuteron-carbon elastic scattering is studied within the framework of the Glauber approximation. The full Glauber multiple scattering series is evaluated with simplified nuclear wave functions. An approximation in which deuteron-nucleus scattering is expressed in terms of nucleon-nucleus scattering amplitudes is shown to be accurate for small momentum transfers. The effects of the Coulomb field and the deuteron D state are investigated. The impressive fit to the 650 MeV d-(C-12) elastic scattering data, obtained in a recent calculation, is shown to be due to additional approximations made in that analysis.

Varma, G. K.↗
DOE OSTI2022

Electronic structure of intrinsic defects in c -gallium nitride: Density functional theory study without the jellium approximation

Here, we report the first nonjellium, systematic, density functional theory (DFT) study of intrinsic and extrinsic defects and defect levels in zinc-blende (cubic) gallium nitride. We use the local moment counter charge (LMCC) method, the standard Perdew-Becke-Ernzerhoff (PBE) exchange-correlation potential, and two pseudopotentials, where the Ga 3$\textit{d}$ orbitals are either in the core ($d^0$) or explicitly in the valence set ($d^{10}$). We studied 64, 216, 512, and 1000 atom supercells, and demonstrated convergence to the infinite limit, crucial for delineating deep from shallow states near band edges, and for demonstrating the elimination of finite cell-size errors. Contrary to common claims, we find that exact exchange is not required to obtain defect levels across the experimental band gap. As was true in silicon, silicon carbide, and gallium arsenide, the extremal LMCC defect levels of the aggregate of defects yield an effective LMCC defect band gap that is within 10% of the experimental gap (3.3 eV) for both pseudopotentials. We demonstrate that the gallium vacancy is more complicated than previously reported. There is dramatic metastability–a nearest-neighbor nitrogen atom shifts into the gallium site, forming an antisite, nitrogen vacancy pair, which is more stable than the simple vacancy for positive charge states. Our assessment of the $d^0$ and $d^{10}$ pseudopotentials yields minimal differences in defect structures and defect levels. The better agreement of the $d^0$ lattice constant with experiment suggests that the more computationally economical $d^0$ pseudopotentials are sufficient to achieve the fidelity possible within the physical accuracy of DFT, and thereby enable calculations in larger supercells necessary to demonstrate convergence with respect to finite size supercell errors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Discrete Approximations of Determinantal Point Processes on Continuous Spaces: Tree Representations and Tail Triviality

We prove tail triviality of determinantal point processes μ on continuous spaces. Tail triviality has been proved for such processes only on discrete spaces, and hence we have generalized the result to continuous spaces. To do this, we construct tree representations, that is, discrete approximations of determinantal point processes enjoying a determinantal structure. There are many interesting examples of determinantal point processes on continuous spaces such as zero points of the hyperbolic Gaussian analytic function with Bergman kernel, and the thermodynamic limit of eigenvalues of Gaussian random matrices for Sine{sub 2},Airy{sub 2},Bessel{sub 2}, and Ginibre point processes. Our main theorem proves all these point processes are tail trivial.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1986

Free vibration analysis by BEM using particular integrals

A new method for the free-vibration analysis using the boundary element technique is presented. The method utilizes a fictitious vector function to approximate the inertia forces and then uses the well-known concept of complementary functions and particular integrals to solve the resulting governing differential equations. The necessary particular integrals are defined for the two and three-dimensional analyses, and the present formulation is applied to a number of two-dimensional problems to show its accuracy and efficiency in the solution of realistic engineering problems.

Ahmad, S.↗
NASA NTRS1991

The effects of compressor seventh-stage bleed air extraction on performance of the F100-PW-220 afterburning turbofan engine

A study was conducted to determine the effects of seventh-stage compressor bleed on the performance of the F100 afterburning turbofan engine. The effects of bleed on thrust, specific fuel consumption, fan turbine inlet temperature, bleed total pressure, and bleed total temperature were obtained from the engine manufacturer's status deck computer simulation. These effects were determined for power settings of intermediate, partial afterburning, and maximum afterburning for Mach numbers between 0.6 and 2.2 and for altitudes of 30,000, 40,000, and 50,000 ft. It was found that thrust loss and specific fuel consumption increase were approximately linear functions of bleed flow and, based on a percent-thrust change basis, were approximately independent of power setting.

Evans, Alison B.↗