Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer algorithms”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Semigroup approximation and robust stabilization of distributed parameter systems

Theoretical results that enable rigorous statements of convergence and exponential stability of Galerkin approximations of LQR controls for infinite dimensional, or distributed parameter, systems have proliferated over the past ten years. In addition, extensive progress has been made over the same time period in the derivation of robust control design strategies for finite dimensional systems. However, the study of the convergence of robust finite dimensional controllers to robust controllers for infinite dimensional systems remains an active area of research. We consider a class of soft-constrained differential games evolving in a Hilbert space. Under certain conditions, a saddle point control can be given in feedback form in terms of a solution to a Riccati equation. By considering a related LQR problem, we can show a convergence result for finite dimensional approximations of this differential game. This yields a computational algorithm for the feedback gain that can be derived from similar strategies employed in infinite dimensional LQR control design problems. The approach described in this paper also inherits the additional properties of stability robustness common to game theoretic methods in finite dimensional analysis. These theoretical convergence and stability results are verified in several numerical experiments.

Kurdila, A. J.↗
NASA NTRS1994

Development of advanced Navier-Stokes solver

The objective of research was to develop and validate new computational algorithms for solving the steady and unsteady Euler and Navier-Stokes equations. The end-products are new three-dimensional Euler and Navier-Stokes codes that are faster, more reliable, more accurate, and easier to use. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible/incompressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. Convergence rates and the robustness of the codes are enhanced by the use of an implicit full approximation storage multigrid method.

Yoon, Seokkwan↗
NASA NTRS1994

Low computation vision-based navigation for a Martian rover

Construction and design details of the Mobot Vision System, a small, self-contained, mobile vision system, are presented. This system uses the view from the top of a small, roving, robotic vehicle to supply data that is processed in real-time to safely navigate the surface of Mars. A simple, low-computation algorithm for constructing a 3-D navigational map of the Martian environment to be used by the rover is discussed.

Gavin, Andrew S.↗
NASA NTRS1994

Effect of stratification and geometrical spreading on sonic boom rise time

The purpose of our investigation is to determine the effect of unsteadiness (not associated with turbulence) on rise time. The unsteadiness considered here is due to (1) geometrical spreading, (2) stratification, which includes variation in density, temperature, and relative humidity, and (3) N shaped waveform. A very general Burgers equation, which includes all these effects, is the propagation model for our study. The equation is solved by a new computational algorithm in which all the calculations are done in the time domain. The present paper is a progress report in which some of the factors contributing to unsteadiness are studied, namely geometrical spreading and variation in relative humidity. The work of Pierce and Kang, which motivated our study, is first reviewed. We proceed with a discussion of the Burgers equation model and the algorithm for solving the equation. Some comparison tests to establish the validity of the algorithm are presented. The algorithm is then used to determine the distance required for a steady-state shock, on encountering an abrupt change in relative humidity, to reach a new steady state based on the new humidity. It is found that the transition distance for plane shocks of amplitude 70 Pa is about 4 km when the change in relative humidity is 10 percent. Shocks of amplitude 140 Pa require less distance. The effect of spherical and cylindrical spreading is also considered. We demonstrate that a spreading shock wave never reaches steady state and that its rise time will be less than the equivalent steady state shock. Finally we show that an N wave has a slightly shorter rise time than a step shock of the same amplitude.

Cleveland, Robin O.↗
NASA NTRS1995

Earth horizon modeling and application to static Earth sensors on TRMM spacecraft

Data from Earth sensor assemblies (ESA's) often are used in the attitude determination (AD) for both spinning and Earth-pointing spacecraft. The ESA's on previous such spacecraft for which the ground-based AD operation was performed by the Flight Dynamics Division (FDD) used the Earth scanning method. AD on such spacecraft requires a model of the shape of the Earth disk as seen from the spacecraft. AD accuracy requirements often are too severe to permit Earth oblateness to be ignored when modeling disk shape. Section 2 of this paper reexamines and extends the methods for Earth disk shape modeling employed in AD work at FDD for the past decade. A new formulation, based on a more convenient Earth flatness parameter, is introduced, and the geometric concepts are examined in detail. It is shown that the Earth disk can be approximated as an ellipse in AD computations. Algorithms for introducing Earth oblateness into the AD process for spacecraft carrying scanning ESA's have been developed at FDD and implemented into the support systems. The Tropical Rainfall Measurement Mission (TRMM) will be the first spacecraft with AD operation performed at FDD that uses a different type of ESA - namely, a static one - containing four fixed detectors D(sub i) (i = 1 to 4). Section 3 of this paper considers the effect of Earth oblateness on AD accuracy for TRMM. This effect ideally will not induce AD errors on TRMM when data from all four D(sub i) are present. When data from only two or three D(sub i) are available, however, a spherical Earth approximation can introduce errors of 0.05 to 0.30 deg on TRMM. These oblateness-induced errors are eliminated by a new algorithm that uses the results of Section 2 to model the Earth disk as an ellipse.

Keat, J.↗
NASA NTRS1995

High-order essentially non-oscillatory methods for computational aeroacoustics

The desire to obtain acoustic information from the numerical solution of a nonlinear system of equations is a demanding proposition for a computational algorithm. High-order accuracy is required for the propagation of high-frequency, low-amplitude waves. In addition, it is desirable to highly resolve discontinuities that can develop in the solutions of the Euler or Navier-Stokes equations. The class of essentially non-oscillatory (ENO) shock-capturing schemes has been designed to have both of these properties. The dual capacity of ENO schemes for high-order accuracy and non-oscillatory shock-capturing is achieved through the use of adaptive stenciling, which makes these schemes highly nonlinear. These schemes are briefly described and referenced herein. A fourth-order algorithm is then applied to the solution of an acoustic wave in a quasi-one-dimensional converging-diverging nozzle.

Casper, Jay↗
NASA NTRS1993

Computational mechanics analysis tools for parallel-vector supercomputers

Computational algorithms for structural analysis on parallel-vector supercomputers are reviewed. These parallel algorithms, developed by the authors, are for the assembly of structural equations, 'out-of-core' strategies for linear equation solution, massively distributed-memory equation solution, unsymmetric equation solution, general eigen-solution, geometrically nonlinear finite element analysis, design sensitivity analysis for structural dynamics, optimization algorithm and domain decomposition. The source code for many of these algorithms is available from NASA Langley.

Storaasli, O. O.↗
NASA NTRS1995

Program Implements Variable-Sampling Procedures

MIL-STD-414 Variable Sampling Procedures (M414) computer program developed to automate calculations and acceptance/rejection procedures of MIL-STD-414, "Sampling Procedures and Tables for Inspection by Variables for Percent Defective." M414 automates entire calculation-and-decision process by use of computational algorithms determining threshold acceptability values for lots. Menu-driven and user-friendly. Reduces burden of manual operations, promoting variable-sampling practice in industry in lieu of "go/no-go" inspection. Written in BASIC.

Huang, Zhaofeng↗
NASA NTRS1995

Device Would Monitor Body Parameters Continuously

Proposed miniature electronic circuit continuously measures temperature of human subject. Once mounted on subject's skin with medical adhesive tape, electronic thermometer remains in thermal equilibrium with subject's body; thereafter, no need to wait until thermometer reaches body temperature before taking reading. Design provides for switches used to set alarm alerting medical attendants if subject's temperature exceeds critical level. For use on very young child, electronic thermometer sewed into shirt or other suitable garment; device held in contact with skin, and child could not swallow it. Replacement of sensor and computing algorithm changes temperature monitor to cardiorespiratory monitor.

Cook, Joseph S., Jr.↗
NASA NTRS1996

Computational considerations for the simulation of shock-induced sound

The numerical study of aeroacoustic problems places stringent demands on the choice of a computational algorithm, because it requires the ability to propagate disturbances of small amplitude and short wavelength. The demands are particularly high when shock waves are involved, because the chosen algorithm must also resolve discontinuities in the solution. The extent to which a high-order-accurate shock-capturing method can be relied upon for aeroacoustics applications that involve the interaction of shocks with other waves has not been previously quantified. Such a study is initiated in this work. A fourth-order-accurate essentially nonoscillatory (ENO) method is used to investigate the solutions of inviscid, compressible flows with shocks in a quasi-one-dimensional nozzle flow. The design order of accuracy is achieved in the smooth regions of a steady-state test case. However, in an unsteady test case, only first-order results are obtained downstream of a sound-shock interaction. The difficulty in obtaining a globally high-order-accurate solution in such a case with a shock-capturing method is demonstrated through the study of a simplified, linear model problem. Some of the difficult issues and ramifications for aeroacoustics simulations of flows with shocks that are raised by these results are discussed.

Casper, Jay↗
NASA NTRS1996

Solwnd: A 3D Compressible MHD Code for Solar Wind Studies. Version 1.0: Cartesian Coordinates

Solwnd 1.0 is a three-dimensional compressible MHD code written in Fortran for studying the solar wind. Time-dependent boundary conditions are available. The computational algorithm is based on Flux Corrected Transport and the code is based on the existing code of Zalesak and Spicer. The flow considered is that of shear flow with incoming flow that perturbs this base flow. Several test cases corresponding to pressure balanced magnetic structures with velocity shear flow and various inflows including Alfven waves are presented. Version 1.0 of solwnd considers a rectangular Cartesian geometry. Future versions of solwnd will consider a spherical geometry. Some discussions of this issue is presented.

Deane, Anil E.↗
NASA NTRS1995

Effects of Meteoroid Erosion in Planetary Rings

This grant supported continuing studies of the effects of ballistic transport on the evolution of Saturn's rings. Ballistic transport, as used in this context, refers to the net transport of mass and angular momentum caused by the exchange of meteoroid impact ejecta between neighboring ring regions (Ip 1983, 1984, Morfill et al. 1983, Lissauer 1984, Durisen 1984a,b). The characteristic time scale associated with this process is the gross erosion time t(sub g) the time it would take a ring region to be completed eroded if all impact ejecta were lost. This time scale is estimated to be about 10(exp 5) to 10(exp 6) years for a ring region with normal optical depth tau approximately 1. Earlier work by myself and collaborators developed the physical theory and simulation techniques to model this process (Durisen et al. 1989, Cuzzi and Durisen 1990). Detailed simulations, supported in part by this grant, have demonstrated that ballistic transport can produce observed structures in Saturn's rings, especially at and near the inner edges of the A and B Rings (Durisen et al. 1992, 1996). The structures of interest in the real rings are illustrated in Figures 1 and 2. Most of these structures were previously unexplained. The computational results plus analytic treatments place useful constraints on fundamental ring properties, including the indication of a relatively young ring age less than or equal 10(exp 8) years (see reviews by Nicholson and Dones 1991, Esposito 1993, Cuzzi 1995, and Porco 1995). This grant also supported development of the faster computational algorithms necessary to permit longer evolutions. Resulting simplications in the ballistic transport equations permitted an analytic linear stability analysis (Durisen 1995), which has provided considerable insight into ballistic transport processes and applications to Saturn's rings. All these accomplishments are described in more detail below.

Durisen, Richard H.↗
NASA NTRS1996

Integrated Control Using the SOFFT Control Structure

The need for integrated/constrained control systems has become clearer as advanced aircraft introduced new coupled subsystems such as new propulsion subsystems with thrust vectoring and new aerodynamic designs. In this study, we develop an integrated control design methodology which accomodates constraints among subsystem variables while using the Stochastic Optimal Feedforward/Feedback Control Technique (SOFFT) thus maintaining all the advantages of the SOFFT approach. The Integrated SOFFT Control methodology uses a centralized feedforward control and a constrained feedback control law. The control thus takes advantage of the known coupling among the subsystems while maintaining the identity of subsystems for validation purposes and the simplicity of the feedback law to understand the system response in complicated nonlinear scenarios. The Variable-Gain Output Feedback Control methodology (including constant gain output feedback) is extended to accommodate equality constraints. A gain computation algorithm is developed. The designer can set the cross-gains between two variables or subsystems to zero or another value and optimize the remaining gains subject to the constraint. An integrated control law is designed for a modified F-15 SMTD aircraft model with coupled airframe and propulsion subsystems using the Integrated SOFFT Control methodology to produce a set of desired flying qualities.

Halyo, Nesim↗
NASA NTRS1996

Searching for Faint Companions to Nearby Stars with the Hubble Space Telescope

A search for faint companions (FC's) to selected stars within 5 pc of the Sun using the Hubble Space Telescope's Planetary Camera (PC) has been initiated. To assess the PC's ability to detect FCs, we have constructed both model and laboratory-simulated images and compared them to actual PC images. We find that the PC's point-spread function (PSF) is 3-4 times brighter over the angular range 2-5 sec than the PSF expected for a perfect optical system. Azimuthal variations of the PC's PSF are 10-20 times larger than expected for a perfect PSF. These variations suggest that light is scattered nonuniformly from the surface of the detector. Because the anomalies in the PC's PSF cannot be precisely simulated, subtracting a reference PSF from the PC image is problematic. We have developed a computer algorithm that identifies local brightness anomalies within the PSF as potential FCs. We find that this search algorithm will successfully locate FCs anywhere within the circumstellar field provided that the average pixel signal from the FC is at least 10 sigma above the local background. This detection limit suggests that a comprehensive search for extrasolar Jovian planets with the PC is impractical. However, the PC is useful for detecting other types of substellar objects. With a stellar signal of 10(exp 9) e(-), for example, we may detect brown dwarfs as faint as M(sub I) = 16.7 separated by 1 sec from alpha Cen A.

Schroeder, Daniel J.↗
NASA NTRS1996

Symmetries in Turbulent Boundary Layer Flows

The objective is the development of a new theory which enables the algorithmic computation of all self-similar mean velocity profiles. The theory is based on Liegroup analysis and unifies a large set of self-similar solutions for the mean velocity of stationary parallel turbulent shear flows. The results include the logarithmic law of the wall, an algebraic law, the viscous sublayer, the linear region in the middle of a Couette flow and in the middle of a rotating channel flow, and a new exponential mean velocity profile not previously reported. Experimental results taken in the outer parts of a high Reynolds number flat-plate boundary layer, strongly support the exponential profile. From experimental as well as from DNS data of a turbulent channel flow the algebraic scaling law could be confirmed in both the center region and in the near wall region. In the case of the logarithmic law of the wall, the scaling with the wall distance arises as a result of the analysis and has not been assumed in the derivation. The crucial part of the derivation of all the different mean velocity profiles is to consider the invariance of the equation for the velocity fluctuations at the same time as the invariance of the equation for the velocity product equations. The latter is the dyad product of the velocity fluctuations with the equation for the velocity fluctuations. It has been proven that all the invariant solutions are also consistent with similarity of all velocity moment equations up to any arbitrary order.

Oberlack, M.↗
NASA NTRS1997

Using High-Order Methods on Lower-Order Geometries

The desire to obtain acoustic information from the numerical solution of a nonlinear system of equations is a demanding proposition for a computational algorithm. High-order accuracy is required for the propagation of high-frequency, low-amplitude waves. The accuracy of an algorithm can be compromised by low-order errors that naturally occur in the solution of a particular problem. Such errors arise from two sources: the presence of discontinuities in the flow field or because the geometry on which the problem is defined is not everywhere smooth to the order of the scheme. The performance of high-order accurate essentially non-oscillatory (ENO) schemes on piecewise smooth solutions is well documented. Herein, the performance of these methods on smooth solutions defined on piecewise smooth geometries is investigated. The propagation of sound in a quasi-one-dimensional nozzle is considered as a test case. Some of the issues involved in the extension to two spatial dimensions are discussed.

Casper, Jay↗
NASA NTRS1997

Computational Considerations for the Simulation of Discontinuous Flows

The numerical study of aeroacoustic problems places stringent demands on the choice of a computational algorithm, because it requires the ability to propagate disturbances of small amplitude and short wavelength. The demands are particularly high when shock waves are involved, because the chosen algorithm must also resolve discontinuities in the solution. In a previous work the capabilities and deficiencies of shock-capturing methods for aeroacoustic problems were demonstrated using a high-order essentially nonoscillatory (ENO) numerical method. It was shown that first-order results are obtained when simulating time-dependent flows with discontinuities. The present study reaffirms this conclusion by comparing the ENO results with those obtained using a conventional linear scheme. A sixth-order-accurate compact implicit finite difference scheme is used to investigate various discontinuous flows. The design order of accuracy is achieved in the smooth regions of a steady-state, quasi-one-dimensional Euler test case, as well as in the time-dependent Burgers' equation. However, in the unsteady Euler sound-shock interaction, first-order results are obtained downstream of the shock. A comparison is made between the linear and nonlinear results, noting the advantages of each method. A discontinuous linear model problem is then used to identify the cause of the first-order results. Here, the nature of the solution error is quantified as being predominantly a numerical phase shift, and a post-processing procedure is demonstrated which increases the solution accuracy downstream of the discontinuity to second-order.

Carpenter Mark H.↗
NASA NTRS1997

A Coupled/Uncoupled Computational Scheme for Deformation and Fatigue Damage Analysis of Unidirectional Metal-Matrix Composites

A fatigue damage computational algorithm utilizing a multiaxial, isothermal, continuum-based fatigue damage model for unidirectional metal-matrix composites has been implemented into the commercial finite element code MARC using MARC user subroutines. Damage is introduced into the finite element solution through the concept of effective stress that fully couples the fatigue damage calculations with the finite element deformation solution. Two applications using the fatigue damage algorithm are presented. First, an axisymmetric stress analysis of a circumferentially reinforced ring, wherein both the matrix cladding and the composite core were assumed to behave elastic-perfectly plastic. Second, a micromechanics analysis of a fiber/matrix unit cell using both the finite element method and the generalized method of cells (GMC). Results are presented in the form of S-N curves and damage distribution plots.

Wilt, Thomas E.↗