Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximate computing”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Application of the probabilistic approximate analysis method to a turbopump blade analysis

An eigenvalue analysis of a typical space propulsion system turbopump blade is presented using an approximate probabilistic analysis methodology. The methodology was developed originally to investigate the feasibility of computing probabilistic structural response using closed-form approximate models. This paper extends the methodology to structures for which simple closed-form solutions do not exist. The finite element method will be used for this demonstration, but the concepts apply to any numerical method. The results agree with detailed analysis results and indicate the usefulness of using a probabilistic approximate analysis in determining efficient solution strategies.

Thacker, B. H.↗
NASA NTRS1992

More About Gravitational Bumps And Repeating Orbits

Improved algorithm reduces computation needed to estimate dynamics affecting trajectory of satellite in orbit around Earth with approximately repeating ground track. Underlying concepts presented in "Computing Gravitational Bumps From Repeating-Orbit Data", (NPO-17925). Satellite presumed to undergo same gravitational perturbation upon each near-repetition. Tracking data from multiple nearly repeating arcs collected into segments called "bins" and processed together to reduce effects of random errors, and thereby obtain more-accurate estimates of dynamical parameters. Bumps in gravitational field of Earth estimated with less computation.

Wu, Jiun-Tsong↗
DOE OSTI2018

Emulator Generation Gadget

This software automates the production of standalone emulators for complex computer models. The end result of running the software is a piece of code called an emulator, a statistical approximation which predicts the output of complex computer simulation, but at a much faster speed. The emulator is built using a training set of inputs and outputs from the simulations. The input-output pairs are used in a Bayesian scheme to estimate the parameters of a Gaussian process based regression model that learns the input-output relationship. The estimated parameters are combined with a generic C skeleton to produce a compile-able instantiation of the emulator.

Sweeney, Christine↗
NASA NTRS1993

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. The Roe approximate Riemann solution scheme or the computationally less expensive Advection Upstream Splitting Method (AUSM) flux-splitting scheme are used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passage and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗
NASA NTRS1994

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. Roe's approximate Riemann solution scheme or the computationally less expensive advection upstream splitting method (AUSM) flux-splitting scheme is used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passages and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗
DOE OSTI2018

A new approach to the generalization of Planck’s law of black-body radiation

Highlights: A new approximation approach to generalized Planck’s radiation law has been developed. Average energy of radiation derived using the nonextensive partition function. The spectral energy density and spectral radiance have been computed. The present approach has been compared with approximate schemes and the exact method. Deviations of present approach and AA and FA schemes from exact result computed. -- Abstract: In this study, Planck’s law of black-body radiation has been modified within the framework of nonextensive statistical mechanics. The average energy of radiation has been derived by introducing the nonextensive partition function in the statistical relation of internal energy. The spectral energy density and spectral radiance have also been computed. The derived expression has been compared with the earlier developed approximate schemes (i.e. asymptotic approach and factorization approach) and with the recently obtained exact result. We utilize the exact and approximate Stefan–Boltzmann laws in order to compare with the new approach introduced here.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1976

Lamination residual stresses in hybrid composites, part 1

An experimental investigation was conducted to study lamination residual stresses for various material and loading parameters. The effects of hybridization on residual stresses and residual properties after thermal cycling under load were determined in angle-ply graphite/Kevlar/epoxy and graphite/S-glass/epoxy laminates. Residual strains in the graphite plies are not appreciably affected by the type and number of hybridizing plies. Computed residual stresses at room temperature in the S-glass plies reach values up to seventy-five percent of the transverse strength of the material. Computed residual stresses in the graphite plies exceed the static strength by approximately ten percent. In the case of Kevlar plies, computed residual stresses far exceed the static strength indicating possible early failure of these plies. Static testing of the hybrids above indicates that failure is governed by the ultimate strain of the graphite plies. In thermally cycled hybrids, in general, residual moduli were somewhat lower and residual strengths were higher than initial values.

Daniel, I. M.↗
NASA NTRS2005

Efficient Computational Model of Hysteresis

A recently developed mathematical model of the output (displacement) versus the input (applied voltage) of a piezoelectric transducer accounts for hysteresis. For the sake of computational speed, the model is kept simple by neglecting the dynamic behavior of the transducer. Hence, the model applies to static and quasistatic displacements only. A piezoelectric transducer of the type to which the model applies is used as an actuator in a computer-based control system to effect fine position adjustments. Because the response time of the rest of such a system is usually much greater than that of a piezoelectric transducer, the model remains an acceptably close approximation for the purpose of control computations, even though the dynamics are neglected. The model (see Figure 1) represents an electrically parallel, mechanically series combination of backlash elements, each having a unique deadband width and output gain. The zeroth element in the parallel combination has zero deadband width and, hence, represents a linear component of the input/output relationship. The other elements, which have nonzero deadband widths, are used to model the nonlinear components of the hysteresis loop. The deadband widths and output gains of the elements are computed from experimental displacement-versus-voltage data. The hysteresis curve calculated by use of this model is piecewise linear beyond deadband limits.

Shields, Joel↗
DOE OSTI2018

Effect of energy dependence of the density of states on pressure-dependent rate constants

The F E integral for the normalized Boltzmann-weighted number of unimolecular states above the threshold energy is the key quantity for computing the collision efficiency in the pressure-dependent unimolecular rate theory developed by Troe, who calls this the energy dependence factor of the density of states. By using the Whitten–Rabinovitch approximation and assuming that the Whitten–Rabinovitch a(E) function is independent of energy, F E can be approximated by an analytical formula; this approximate formula is widely used because of its convenience and computational efficiency. Here we test its validity by comparing the rate constants computed by using the approximate F E to the ones determined by using the numerically integrated F E . For small-sized molecules and for reactions with high threshold energies E 0 , the differences are negligible at all temperatures, but in other cases, the approximate formula tends to underestimate F E and thus overestimates the collision efficiency, and this leads to smaller pressure falloff. When a(E) at high energies differs appreciably from a(E 0 ), we find that the underestimation of pressure-dependent rate constants by using the approximate formula can be greater than a factor of 5 at high temperatures. The physical insight we draw from this study is that, for reactions with threshold energies below about 30 kcal mol -1 , the rate of collisional energy transfer can be appreciably slowed down by the increase in the density of states at higher energies, and this increases the falloff effect by which finite-pressure rate constants become lower than the high-pressure limit, especially at higher temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

CREW CHIEF: A computer graphics simulation of an aircraft maintenance technician

Approximately 35 percent of the lifetime cost of a military system is spent for maintenance. Excessive repair time is caused by not considering maintenance during design. Problems are usually discovered only after a mock-up has been constructed, when it is too late to make changes. CREW CHIEF will reduce the incidence of such problems by catching design defects in the early design stages. CREW CHIEF is a computer graphic human factors evaluation system interfaced to commercial computer aided design (CAD) systems. It creates a three dimensional man model, either male or female, large or small, with various types of clothing and in several postures. It can perform analyses for physical accessibility, strength capability with tools, visual access, and strength capability for manual materials handling. The designer would produce a drawing on his CAD system and introduce CREW CHIEF in it. CREW CHIEF's analyses would then indicate places where problems could be foreseen and corrected before the design is frozen.

Aume, Nilss M.↗
DOE OSTI2018

A Galerkin method with two-dimensional Haar basis functions for the computation of the Karhunen–Loève expansion

We study the numerical approximation of a homogeneous Fredholm integral equation of second kind associated with the Karhunen–Loève expansion of Gaussian random fields. We develop, validate, and discuss an algorithm based on the Galerkin method with two-dimensional Haar wavelets as basis functions. The shape functions are constructed from the orthogonal decomposition of tensor-product spaces of one-dimensional Haar functions, and a recursive algorithm is employed to compute the matrix of the discrete eigenvalue system without the explicit calculation of integrals, allowing the implementation of a fast and efficient algorithm that provides considerable reduction in CPU time, when compared with classical Galerkin methods. Numerical experiments confirm the convergence rate of the method and assess the approximation error and the sparsity of the eigenvalue system when the wavelet expansion is truncated. We illustrate the numerical solution of a diffusion problem with random input data with the present method. In this problem, accuracy was retained after dropping the coefficients below a threshold value that was numerically determined. A similar method with scaling functions rather than wavelet functions does not need a discrete wavelet transform and leads to eigenvalue systems with better conditioning but lower sparsity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Computational aspects of sensitivity calculations in linear transient structural analysis

A study was performed focusing on the calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, transient response problems. One significant goal of the study was to develop and evaluate sensitivity calculation techniques suitable for large-order finite element analyses. Accordingly, approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite element model. Much of the research focused on the accuracy of both response quantities and sensitivities as a function of number of vectors used. Two types of sensitivity calculation techniques were developed and evaluated. The first type of technique is an overall finite difference method where the analysis is repeated for perturbed designs. The second type of technique is termed semi-analytical because it involves direct, analytical differentiation of the equations of motion with finite difference approximation of the coefficient matrices. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models. In several cases this fixed mode approach resulted in very poor approximations of the stress sensitivities. Almost all of the original modes were required for an accurate sensitivity and for small numbers of modes, the accuracy was extremely poor. To overcome this poor accuracy, two semi-analytical techniques were developed. The first technique accounts for the change in eigenvectors through approximate eigenvector derivatives. The second technique applies the mode acceleration method of transient analysis to the sensitivity calculations. Both result in accurate values of the stress sensitivities with a small number of modes and much lower computational costs than if the vibration modes were recalculated and then used in an overall finite difference method.

Greene, William H.↗
DOE OSTI2025

Computation of Auger Electron Spectra in Organic Molecules with Multiconfiguration Pair-Density Functional Theory

Efficient and accurate computation of molecular Auger electron spectra for larger systems is limited by the rapid increase in the number of doubly ionized final states as the system size grows. Here, in this work, we benchmark the application of multiconfiguration pair-density functional theory with a restricted active space (RAS) reference wave function for computing the carbon K-edge decay spectra of 20 organic molecules. Decay rates are computed within the one-center approximation. We evaluate the performance of different basis sets and on-top functionals and find that multiconfiguration pair-density functional theory achieves accuracy comparable to RAS followed by second-order perturbation theory, but at significantly lower computational cost.

Fouda, Adam E. A. [Argonne National Laboratory (AN↗
DOE OSTI2024

Large-scale harmonic balance simulations with Krylov subspace and preconditioner recycling

The multi-harmonic balance method combined with numerical continuation provides an efficient framework to compute a family of time-periodic solutions, or response curves, for large-scale, nonlinear mechanical systems. The predictor and corrector steps repeatedly solve a sequence of linear systems that scale by the model size and number of harmonics in the assumed Fourier series approximation. In this paper, a novel Newton–Krylov iterative method is embedded within the multi-harmonic balance and continuation algorithm to efficiently compute the approximate solutions from the sequence of linear systems that arise during the prediction and correction steps. Further, the method recycles, or reuses, both the preconditioner and the Krylov subspace generated by previous linear systems in the solution sequence. A delayed frequency preconditioner refactorizes the preconditioner only when the performance of the iterative solver deteriorates. The GCRO-DR iterative solver recycles a subset of harmonic Ritz vectors to initialize the solution subspace for the next linear system in the sequence. The performance of the iterative solver is demonstrated on two exemplars with contact-type nonlinearities and benchmarked against a direct solver with traditional Newton–Raphson iterations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1940

The frequencies of cantilever wings in beam and torsional vibrations

Methods are described for calculating the period and frequency of vibration of cantilever wings and similar structures in which the weight and moment of inertia vary along the span. Both the beam and torsional frequencies may be calculated by these methods. The procedure is illustrated by examples. It is shown that a surprisingly close approximation to the beam frequency may be obtained by a very brief calculation in which the curvature of the wing in vibration is assumed to be constant. A somewhat longer computation permits taking account of the true curvature of the beam by a series of successive approximations which are shown to be strongly convergent. Analogous methods are applied to calculations of the torsional frequency. For the first approximation it is assumed that the angle of twist varies linearly alone the semispan. True variation of the twist is computed by successive approximations which are strongly convergent, as in the case of beam vibrations.

Burgess, C P↗
DOE OSTI2018

Eccentric Form of the Loss of Stability of a Rotating Elastoplastic Radially Inhomogeneous Stepped Circular Disk

We propose a procedure for the investigation of the possible loss of stability of a rotating elastoplastic radially inhomogeneous stepped circular disk by the method of small parameter. The characteristic equation for the critical radius of the plastic zone is obtained in the first approximation. The values of the critical angular rotational velocity are numerically found for various parameters of the disk.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Exact sum rules with approximate ground states

Electromagnetic and weak transitions tell us a great deal about the structure of atomic nuclei. Yet modeling transitions can be difficult: it is often easier to compute the ground state, if only as an approximation, than excited states. One alternative is through transition sum rules, in particular the non-energy-weighted and energy-weighted sum rules, which can be computed as expectation values of operators. We investigate by computing sum rules for a variety of nuclei, comparing the numerically exact full configuration-interaction shell model, as a reference, to Hartree-Fock, projected Hartree-Fock, and the nucleon pair approximation. Furthermore, these approximations yield reasonable agreement, which we explain by prior work on the systematics of transition moments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗