Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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,027 records · Page 57

A high Reynolds number numerical solution of the Navier-Stokes equations in stream function-vorticity form

Problems encoutered in investigations of high Reynolds number, incompressible flow are reviewed. A numerical solution computer program was modified to solve the stream function-vorticity form of the Navier-Stokes equations. Using a body fitted coordinate system with a U-shaped outer boundary, a simulation of incompressible flow at a Reynolds number of one million and a body angle of attack of zero was achieved.

Bearden, J. H.↗

Singular perturbation techniques for on-line optimal flight path control

This paper presents a partial evaluation on the use of singular perturbation methods for developing computer algorithms for on-line optimal control of aircraft. The evaluation is based on a study of the minimum time intercept problem using F-4 aerodynamic and propulsion data as a base line. The extensions over previous work on this subject are that aircraft turning dynamics (in addition to position and energy dynamics) are included in the analysis, the algorithm is developed for a moving end point and is adaptive to unpredictable target maneuvers, and short range maneuvers that do not have a cruise leg are included. Particular attention is given to identifying those quantities that can be precomputed and stored (as a function of aircraft total energy), thus greatly reducing the onboard computational load. Numerical results are given that illustrate the nature of the optimal intercept flight paths, and an estimate is given for the execution time and storage requirements of the control algorithm.

Calise, A. J.↗

A polymorphic reconfigurable emulator for parallel simulation

Microprocessor and arithmetic support chip technology was applied to the design of a reconfigurable emulator for real time flight simulation. The system developed consists of master control system to perform all man machine interactions and to configure the hardware to emulate a given aircraft, and numerous slave compute modules (SCM) which comprise the parallel computational units. It is shown that all parts of the state equations can be worked on simultaneously but that the algebraic equations cannot (unless they are slowly varying). Attempts to obtain algorithms that will allow parellel updates are reported. The word length and step size to be used in the SCM's is determined and the architecture of the hardware and software is described.

Parrish, E. A., Jr.↗

Numerical algorithms for finite element computations on arrays of microprocessors

The development of a multicolored successive over relaxation (SOR) program for the finite element machine is discussed. The multicolored SOR method uses a generalization of the classical Red/Black grid point ordering for the SOR method. These multicolored orderings have the advantage of allowing the SOR method to be implemented as a Jacobi method, which is ideal for arrays of processors, but still enjoy the greater rate of convergence of the SOR method. The program solves a general second order self adjoint elliptic problem on a square region with Dirichlet boundary conditions, discretized by quadratic elements on triangular regions. For this general problem and discretization, six colors are necessary for the multicolored method to operate efficiently. The specific problem that was solved using the six color program was Poisson's equation; for Poisson's equation, three colors are necessary but six may be used. In general, the number of colors needed is a function of the differential equation, the region and boundary conditions, and the particular finite element used for the discretization.

Ortega, J. M.↗

Vertical slender jets

The shape of a vertical slender jet of fluid falling steadily under the force of gravity is studied. The problem is formulated as a nonlinear free boundary-value problem for the potential. Surface tension effects are neglected. The use of perturbation expansions results in a system of equations that can be solved by an efficient numerical procedure. Computations were made for jets issuing from orifices in various shapes including an ellipse, a rectangle, and an equilateral triangle. Computational results are presented illustrating the propagation of discontinuities and the formation of thin sheets of fluid.

Geer, J. F.↗

Collapse and equilibrium of rotating, adiabatic clouds

A numerical hydrodynamics computer code analysis of the collapse and establishment of equilibrium of adiabatic gas clouds restricted to axial symmetry found that the clouds are originally uniform in density and rotation. The method can compare the dynamic collapse and approach to equilibrium with the data on incompressible uniformly rotating equilibrium clouds and on equilibrium structures of differentially rotating polytropes. It is concluded that the stellar formation theory indicates that the low alpha fragments produced at the termination of the dynamic isothermal collapse phase of interstellar clouds may undergo significant dynamic collapse in an adiabatic regime leading to transitory ring formation and additional fragmentation on a smaller scale.

Boss, A. P.↗

Research in nonlinear structural and solid mechanics

Recent and projected advances in applied mechanics, numerical analysis, computer hardware and engineering software, and their impact on modeling and solution techniques in nonlinear structural and solid mechanics are discussed. The fields covered are rapidly changing and are strongly impacted by current and projected advances in computer hardware. To foster effective development of the technology perceptions on computing systems and nonlinear analysis software systems are presented.

Mccomb, H. G., Jr.↗

The electromagnetic force field, fluid flow field and temperature profiles in levitated metal droplets

A mathematical representation was developed for the electromagnetic force field, the flow field, the temperature field (and for transport controlled kinetics), in a levitation melted metal droplet. The technique of mutual inductances was employed for the calculation of the electromagnetic force field, while the turbulent Navier - Stokes equations and the turbulent convective transport equations were used to represent the fluid flow field, the temperature field and the concentration field. The governing differential equations, written in spherical coordinates, were solved numerically. The computed results were in good agreement with measurements, regarding the lifting force, and the average temperature of the specimen and carburization rates, which were transport controlled.

El-Kaddah, N.↗

Experimental boundary integral equation applications in speckle interferometry

A complete data analysis system utilizing laser speckle interferometry and Experimental Integral Equation Method (EBIE) is described. A host computer provides the optical data analysis and serves as an input device to a PDP11/VAX computer for numerical analysis. The basic theory of pointwise filtering and digital correlation is described as experimental data input to the BIE method.

Fraley, J. E.↗

Determination of rotor wake induced empennage airloads

A computer program has been developed that predicts the unsteady aerodynamic forces that impinge on the empennage surface due to its interactions with the main rotor wake. The program was utilized to determine the vibration airloads acting on the Black Hawk horizontal stabilizer and on the CH-53A stabilizer under high speed forward flight conditions. The results demonstrate that it is possible to compute the high frequency empennage vibration airloads efficiently by utilizing suitable analysis techniques; a CH-53A case required only 1.3 minutes of computer time. Numerical problems associated with the unsteady Kutta condition have been minimized. Even though the analysis is of transient nature, the results for harmonic cases show a very fast convergence to periodicity. In general, good correlation was shown between the predicted vibratory airloads and the measured stabilizer airloads for a CH-53A helicopter operating at high speed flight conditions.

Gangwani, S. T.↗

Closed loop control performance sensitivity to parameter variations

A very efficient technique for computing the closed loop performance sensitivities to parameter variations of a dynamic system with a reduced order controller has been developed. The eigensystem of the closed loop system is computed once. With this information, the closed loop filter and state rms responses, and the first and second derivatives of these rms values with respect to given parameters are computed. Detailed numerical examples using the JPL flexible beam and a 55 meter offset fed, wrap-rib antenna are included.

Schaechter, D. B.↗

Study of tethered satellite active attitude control

Existing software was adapted for the study of tethered subsatellite rotational dynamics, an analytic solution for a stable configuration of a tethered subsatellite was developed, the analytic and numerical integrator (computer) solutions for this "test case' was compared in a two mass tether model program (DUMBEL), the existing multiple mass tether model (SKYHOOK) was modified to include subsatellite rotational dynamics, the analytic "test case,' was verified, and the use of the SKYHOOK rotational dynamics capability with a computer run showing the effect of a single off axis thruster on the behavior of the subsatellite was demonstrated. Subroutines for specific attitude control systems are developed and applied to the study of the behavior of the tethered subsatellite under realistic on orbit conditions. The effect of all tether "inputs,' including pendular oscillations, air drag, and electrodynamic interactions, on the dynamic behavior of the tether are included.

Colombo, G.↗

CO 1-0 band isotopic lines as intensity standards

Fundamental band vibration-rotational line intensities are reported for the (C-12)(O-16), (C-13)(O-16), and (C-12)(O-18) molecules. The intensities were computed using numerically obtained matrix elements of a recently, determined electric dipole moment function. New intensity measurements for some of these lines obtained via high-resolution Fourier spectra are also reported. The overall consistency (approximately 2.2 percent) between calculations, which are based on recent experimental determinations of CO 1-0 band intensities, not including our own, those of the authors, and independent experimental results of the authors suggests that the lines of the fundamental band of CO can be used as absolute intensity standards.

Chackerian, C., Jr.↗

Dynamics analysis of electrodynamic satellite tethers. Equations of motion and numerical solution algorithms for the tether

The equations of motion are developed for a perfectly flexible, inelastic tether with a satellite at its extremity. The tether is attached to a space vehicle in orbit. The tether is allowed to possess electrical conductivity. A numerical solution algorithm to provide the motion of the tether and satellite system is presented. The resulting differential equations can be solved by various existing standard numerical integration computer programs. The resulting differential equations allow the introduction of approximations that can lead to analytical, approximate general solutions. The differential equations allow more dynamical insight of the motion.

Nacozy, P. E.↗

A new minimax algorithm

The representation min t s.t. F(I)(x). - t less than or equal to 0 for all i is examined. An active set strategy is designed of functions: active, semi-active, and non-active. This technique will help in preventing zigzagging which often occurs when an active set strategy is used. Some of the inequality constraints are handled with slack variables. Also a trust region strategy is used in which at each iteration there is a sphere around the current point in which the local approximation of the function is trusted. The algorithm is implemented into a successful computer program. Numerical results are provided.

Vardi, A.↗

An outer ion coma model in two dimension

A simple cometary outer ion coma model is constructed by adopting published results from numerical MHD computations. The ion composition and density distributions so derived are used to highlight some of the basic features of cometary plasma flows and the need for cometary ion observations with absolute brightness calibration.

Ip, W.-H.↗

Numerical algorithms for finite element computations on concurrent processors

The work of several graduate students which relate to the NASA grant is briefly summarized. One student has worked on a detailed analysis of the so-called ijk forms of Gaussian elemination and Cholesky factorization on concurrent processors. Another student has worked on the vectorization of the incomplete Cholesky conjugate method on the CYBER 205. Two more students implemented various versions of Gaussian elimination and Cholesky factorization on the FLEX/32.

Ortega, J. M.↗