Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 685 records · Page 38

Manufacturing in space: Fluid dynamics numerical analysis

Natural convection in a spherical container with cooling at the center was numerically simulated using the Lockheed-developed General Interpolants Method (GIM) numerical fluid dynamic computer program. The numerical analysis was simplified by assuming axisymmetric flow in the spherical container, with the symmetry axis being a sphere diagonal parallel to the gravity vector. This axisymmetric spherical geometry was intended as an idealization of the proposed Lal/Kroes growing experiments to be performed on board Spacelab. Results were obtained for a range of Rayleigh numbers from 25 to 10,000. For a temperature difference of 10 C from the cooling sting at the center to the container surface, and a gravitional loading of 0.000001 g a computed maximum fluid velocity of about 2.4 x 0.00001 cm/sec was reached after about 250 sec. The computed velocities were found to be approximately proportional to the Rayleigh number over the range of Rayleigh numbers investigated.

Robertson, S. J.↗

Behavior of bubbles in glassmelts. II - Dissolution of a stationary bubble containing a diffusing and a nondiffusing gas

The effect of a foreign nondiffusing gas on the rate of shrinkage of an oxygen bubble in a soda-lime-silica melt was studied. The rate of change of bubble radius with time was computed using the quasi-stationary approximation. The effects of melt undersaturation and initial fraction of foreign gas in the bubble are considered and compared with those calculated using previously derived expressions.

Weinberg, M. C.↗

Kinematics of velocity and vorticity correlations in turbulent flow

The kinematic problem of calculating second-order velocity moments from given values of the vorticity covariance is examined. Integral representation formulas for second-order velocity moments in terms of the two-point vorticity correlation tensor are derived. The special relationships existing between velocity moments in isotropic turbulence are expressed in terms of the integral formulas yielding several kinematic constraints on the two-point vorticity correlation tensor in isotropic turbulence. Numerical evaluation of these constraints suggests that a Gaussian curve may be the only form of the longitudinal velocity correlation coefficient which is consistent with the requirement of isotropy. It is shown that if this is the case, then a family of exact solutions to the decay of isotropic turbulence may be obtained which contains Batchelor's final period solution as a special case. In addition, the computed results suggest a method of approximating the integral representation formulas in general turbulent shear flows.

Bernard, P. S.↗

A tentative unified sea model for scattering and emission

This work is an attempt to relate active and passive observations of the sea surface to each other, considering that more information can be obtained from a set of instruments when they are not just considered as a juxtaposition of separate and independent boxes, but as a single complex instrument. A careful review of several approaches to the solution of the scattering problem for a random rough surface has been performed. Inherent simplifications, and inconsistencies, are pointed out. To reduce the computation time to reasonable values, an approximate model is selected and some numerical results will be presented in order to illustrate the benefit that can be gained by a unified approach.

Guissard, A.↗

Composite laminate optimization program suitable for microcomputers

The design of a composite panel requires some way of finding the minimum thickness laminate which will withstand the load requirements without failure. The mathematical complexity of this problem dictates the use of nonlinear optimization techniques. Specialized laminate optimization programs were developed which are compact and efficient enough to run on microcomputers. Only stresses at a point and inplane loads and deflections are considered. The programs are simple to use and require no knowledge of optimization. Techniques are developed which find minimum thickness laminates with either ply ratios or ply angles as design variables. A method is presented for finding the optimum orientation for the axis of symmetry of an orthotropic laminate. The orthotropic laminate program uses an approximate failure theory, which speed up computations dramatically.

Flanagan, G. V.↗

Studies of the aurorally-induced ultraviolet emissions on the nightside of Venus

The effect of a monoenergetic flux of electrons on a model atmosphere of the nightside thermosphere of Venus was examined. The neutral model chosen is that of Hedin for high solar activity and l65 degrees solar zenith angle. The model is based on measurements made by the Pioneer Venus Orbiter Neutral Mass Spectrometer. Four species were included in the calculation: CO2, O, CO, and N2. The numerical method that was chosen for energy deposition of the primary electrons is the continuous slowing down approximation. The secondary electron distribution was computed using the empirically determined shape of the differential cross section.

Fox, J. L.↗

Performance evaluation of fault-tolerant systems with application to the IUS

A new method for quantitatively evaluating the performance of fault-tolerant systems is developed and applied to an example. The method assumes that the random failure and diagnostic decision behavior of the system can be modeled by a finite state Markov process. A performance value must be assigned to each of the states of the model. The method then generates the moments of the probability mass function of the cumulative performance and uses these to generate a maximum entropy approximation to the performance PMF. Some computational considerations are discussed. The method is applied to a typical mission for the Inertial Upper Stage to examine the attitude accuracy performance of the inertial system.

Missana, J.-O.↗

Unsteady Navier-Stokes computations past oscillating delta wing at high incidence

The unsteady, thin-layer, compressible Navier-Stokes equations, written in the moving frame of reference for the flow relative motion, is solved for the steady and unsteady supersonic flow around a round-edged delta wing. For supersonic flow, local conical flow solution has been obtained from the three-dimensional equations. Pseudotime stepping is used for the steady flow, while time-accurate stepping is used for the unsteady flow. The computational scheme is an implicit approximately-factored finite volume scheme which uses explicit and implicit dissipation terms. The scheme is verified for the steady flow solution. The scheme is then applied to a delta wing undergoing rolling oscillation at a reduced frequency of 1.137 with 15- deg maximum amplitude about a mean angle of attack of 10 deg for a Mach number of 2 and a Reynolds number of 500,000.

Kandil, Osama A.↗

Spike: Artificial intelligence scheduling for Hubble space telescope

Efficient utilization of spacecraft resources is essential, but the accompanying scheduling problems are often computationally intractable and are difficult to approximate because of the presence of numerous interacting constraints. Artificial intelligence techniques were applied to the scheduling of the NASA/ESA Hubble Space Telescope (HST). This presents a particularly challenging problem since a yearlong observing program can contain some tens of thousands of exposures which are subject to a large number of scientific, operational, spacecraft, and environmental constraints. New techniques were developed for machine reasoning about scheduling constraints and goals, especially in cases where uncertainty is an important scheduling consideration and where resolving conflicts among conflicting preferences is essential. These technique were utilized in a set of workstation based scheduling tools (Spike) for HST. Graphical displays of activities, constraints, and schedules are an important feature of the system. High level scheduling strategies using both rule based and neural network approaches were developed. While the specific constraints implemented are those most relevant to HST, the framework developed is far more general and could easily handle other kinds of scheduling problems. The concept and implementation of the Spike system are described along with some experiments in adapting Spike to other spacecraft scheduling domains.

Johnston, Mark↗

Finite-volume application of high order ENO schemes to multi-dimensional boundary-value problems

The finite volume approach in developing multi-dimensional, high-order accurate essentially non-oscillatory (ENO) schemes is considered. In particular, a two dimensional extension is proposed for the Euler equation of gas dynamics. This requires a spatial reconstruction operator that attains formal high order of accuracy in two dimensions by taking account of cross gradients. Given a set of cell averages in two spatial variables, polynomial interpolation of a two dimensional primitive function is employed in order to extract high-order pointwise values on cell interfaces. These points are appropriately chosen so that correspondingly high-order flux integrals are obtained through each interface by quadrature, at each point having calculated a flux contribution in an upwind fashion. The solution-in-the-small of Riemann's initial value problem (IVP) that is required for this pointwise flux computation is achieved using Roe's approximate Riemann solver. Issues to be considered in this two dimensional extension include the implementation of boundary conditions and application to general curvilinear coordinates. Results of numerical experiments are presented for qualitative and quantitative examination. These results contain the first successful application of ENO schemes to boundary value problems with solid walls.

Casper, Jay↗

Polynomial Transformations For Discrete-Time Linear Systems

Transformations based on polynomial matrices of finite degree developed for use in computing functions for compensation, inversion, and approximation of discrete-time, multivariable, linear systems. Method derived from z-transform transfer-function form of matrices. Applicable to cascade-compensation problems in design of control systems.

Baram, Yoram↗

Spatial adaption procedures on unstructured meshes for accurate unsteady aerodynamic flow computation

Spatial adaption procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaption procedures were developed and implemented within a two-dimensional unstructured-grid upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in a high gradient region or the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational costs. A detailed description is given of the enrichment and coarsening procedures and comparisons with alternative results and experimental data are presented to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady transonic results, obtained using spatial adaption for the NACA 0012 airfoil, are shown to be of high spatial accuracy, primarily in that the shock waves are very sharply captured. The results were obtained with a computational savings of a factor of approximately fifty-three for a steady case and as much as twenty-five for the unsteady cases.

Rausch, Russ D.↗

Spatial adaption procedures on unstructured meshes for accurate unsteady aerodynamic flow computation

Spatial adaption procedures for the accurate and efficient solution of steady and unsteady inviscid flow problems are described. The adaption procedures were developed and implemented within a two-dimensional unstructured-grid upwind-type Euler code. These procedures involve mesh enrichment and mesh coarsening to either add points in high gradient regions of the flow or remove points where they are not needed, respectively, to produce solutions of high spatial accuracy at minimal computational cost. The paper gives a detailed description of the enrichment and coarsening procedures and presents comparisons with alternative results and experimental data to provide an assessment of the accuracy and efficiency of the capability. Steady and unsteady transonic results, obtained using spatial adaption for the NACA 0012 airfoil, are shown to be of high spatial accuracy, primarily in that the shock waves are very sharply captured. The results were obtained with a computational savings of a factor of approximately fifty-three for a steady case and as much as twenty-five for the unsteady cases.

Rausch, Russ D.↗

Domain decomposition methods in aerodynamics

Compressible Euler equations are solved for two-dimensional problems by a preconditioned conjugate gradient-like technique. An approximate Riemann solver is used to compute the numerical fluxes to second order accuracy in space. Two ways to achieve parallelism are tested, one which makes use of parallelism inherent in triangular solves and the other which employs domain decomposition techniques. The vectorization/parallelism in triangular solves is realized by the use of a recording technique called wavefront ordering. This process involves the interpretation of the triangular matrix as a directed graph and the analysis of the data dependencies. It is noted that the factorization can also be done in parallel with the wave front ordering. The performances of two ways of partitioning the domain, strips and slabs, are compared. Results on Cray YMP are reported for an inviscid transonic test case. The performances of linear algebra kernels are also reported.

Venkatakrishnan, V.↗

Implicit solvers for unstructured meshes

Implicit methods for unstructured mesh computations are developed and tested. The approximate system which arises from the Newton-linearization of the nonlinear evolution operator is solved by using the preconditioned generalized minimum residual technique. These different preconditioners are investigated: the incomplete LU factorization (ILU), block diagonal factorization, and the symmetric successive over-relaxation (SSOR). The preconditioners have been optimized to have good vectorization properties. The various methods are compared over a wide range of problems. Ordering of the unknowns, which affects the convergence of these sparse matrix iterative methods, is also investigated. Results are presented for inviscid and turbulent viscous calculations on single and multielement airfoil configurations using globally and adaptively generated meshes.

Venkatakrishnan, V.↗

Optimal Chebyshev polynomials on ellipses in the complex plane

The design of iterative schemes for sparse matrix computations often leads to constrained polynomial approximation problems on sets in the complex plane. For the case of ellipses, we introduce a new class of complex polynomials which are in general very good approximations to the best polynomials and even optimal in most cases.

Fischer, Bernd↗

Resonance charge transfer, transport cross sections, and collision integrals for N(+)(3P)-N(4S0) and O(+)(4S0)-O(3P) interactions

N2(+) and O2(+) potential energy curves have been constructed by combining measured data with the results from electronic structure calculations. These potential curves have been employed to determine accurate charge exchange cross sections, transport cross sections, and collision integrals for ground state N(+)-N and O(+)-O interactions. The cross sections have been calculated from a semiclassical approximation to the scattering using a computer code that fits a spline curve through the discrete potential data and incorporates the proper long-range behavior of the interactions forces. The collision integrals are tabulated for a broad range of temperatures 250-100,000 K and are intended to reduce the uncertainty in the values of the transport properties of nonequilibrium air, particularly at high temperatures.

Stallcop, James R.↗

An engineering aerodynamic heating method for hypersonic flow

A capability to calculate surface heating rates has been incorporated in an approximate three-dimensional inviscid technique. Surface streamlines are calculated from the inviscid solution, and the axisymmetric analog is then used along with a set of approximate convective-heating equations to compute the surface heat transfer. The method is applied to blunted axisymmetric and three-dimensional ellipsoidal cones at angle of attack for the laminar flow of a perfect gas. The method is also applicable to turbulent and equilibrium-air conditions. The present technique predicts surface heating rates that compare favorably with experimental (ground-test and flight) data and numerical solutions of the Navier-Stokes (NS) and viscous shock-layer (VSL) equations. The new technique represents a significant improvement over current engineering aerothermal methods with only a modest increase in computational effort.

Riley, Christopher J.↗