Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear equations”

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 NTRS1989

Interaction between Tollmien-Schlichting waves and free-stream disturbances in boundary-layer flows

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small freestream disturbance and a small localized variation of the surface geometry has been carried out using finite difference methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T-S waves generated by the interaction between the freestream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the Tollmien-Schlichting waves.

Bodonyi, R. J.↗
NASA NTRS1989

The FORCE: A portable parallel programming language supporting computational structural mechanics

This project supports the conversion of codes in Computational Structural Mechanics (CSM) to a parallel form which will efficiently exploit the computational power available from multiprocessors. The work is a part of a comprehensive, FORTRAN-based system to form a basis for a parallel version of the NICE/SPAR combination which will form the CSM Testbed. The software is macro-based and rests on the force methodology developed by the principal investigator in connection with an early scientific multiprocessor. Machine independence is an important characteristic of the system so that retargeting it to the Flex/32, or any other multiprocessor on which NICE/SPAR might be imnplemented, is well supported. The principal investigator has experience in producing parallel software for both full and sparse systems of linear equations using the force macros. Other researchers have used the Force in finite element programs. It has been possible to rapidly develop software which performs at maximum efficiency on a multiprocessor. The inherent machine independence of the system also means that the parallelization will not be limited to a specific multiprocessor.

Jordan, Harry F.↗
NASA NTRS1989

Nonmigrating tides

Recently an interesting advancement in the study of nonmigrating tides has occurred. There have been two distinct approaches in works on this subject. One is based on mechanistic models as considering nonuniform global distribution of water vapor or heating only on land, solving a set of linear equations. It is found that insolation absorption of the nonuniformly distributed water vapor produces only weak nonmigrating tides in the lower thermosphere; the planetary boundary layer heating on land can explain the enhanced tides on land and those with short vertical wavelengths in the stratosphere. The other approach is novel and uses simulation on the general circulation model (GCM). This realistic model can reproduce tides globally and in many details. The enhancement of two nonmigrating modes as eastward traveling modes with a wave number 3 and westward traveling modes with a wave number 5 is in surprisingly good agreement with observation at sea level, at 700 mb and even at 300 mb.

Kato, S.↗
NASA NTRS1989

Parallel-vector solution of large-scale structural analysis problems on supercomputers

A direct linear equation solution method based on the Choleski factorization procedure is presented which exploits both parallel and vector features of supercomputers. The new equation solver is described, and its performance is evaluated by solving structural analysis problems on three high-performance computers. The method has been implemented using Force, a generic parallel FORTRAN language.

Storaasli, Olaf O.↗
NASA NTRS1989

Shock structure in non-circular jets

The shock-cell structure of supersonic jets with non-circular exit geometry is modeled using a linearized analysis. The model takes into account the finite thickness of the jet shear layer using realistic velocity and density profiles. The effects of the shear layer turbulence are included by incorporating eddy-viscosity terms. A finite-difference numerical method is used to solve the steady linearized equations of motion. A body-fitted coordinate system is used to describe the shear layer. The variation of the pressure fluctuation with downstream distance is given for circular jets and for an elliptic jet of aspect ratio 2.0. Comparisons with experimental data are made. Difficulties with the numerical technique are also discussed.

Morris, Philip J.↗
NASA NTRS1990

Output model-following control synthesis for an oblique-wing aircraft

Recent interest in oblique-wing aircraft has focused on the potential aerodynamic performance advantage of a variable-skew oblique wing over a conventional or symmetric sweep wing. Unfortunately, the resulting asymmetric configuration has significant aerodynamic and inertial cross-coupling between the aircraft longitudinal and lateral-directional axes. Presented here is a decoupling control law synthesis technique that integrates stability augmentation, decoupling, and the direct incorporation of desired handling qualities into an output feedback controller. The proposed design technique uses linear quadratic regulator concepts in the framework of explicit model following. The output feedback strategy used is a suboptimal projection from the state space to the output space. Dynamics are then introduced into the controller to improve steady-state performance and increase system robustness. Closed-loop performance is shown by application of the control laws to the linearized equations of motion and nonlinear simulation of an oblique-wing aircraft.

Pahle, Joseph W.↗
NASA NTRS1989

Thin-plate spline quadrature of geodetic integrals

Thin-plate spline functions (known for their flexibility and fidelity in representing experimental data) are especially well-suited for the numerical integration of geodetic integrals in the area where the integration is most sensitive to the data, i.e., in the immediate vicinity of the evaluation point. Spline quadrature rules are derived for the contribution of a circular innermost zone to Stoke's formula, to the formulae of Vening Meinesz, and to the recursively evaluated operator L(n) in the analytical continuation solution of Molodensky's problem. These rules are exact for interpolating thin-plate splines. In cases where the integration data are distributed irregularly, a system of linear equations needs to be solved for the quadrature coefficients. Formulae are given for the terms appearing in these equations. In case the data are regularly distributed, the coefficients may be determined once-and-for-all. Examples are given of some fixed-point rules. With such rules successive evaluation, within a circular disk, of the terms in Molodensky's series becomes relatively easy. The spline quadrature technique presented complements other techniques such as ring integration for intermediate integration zones.

Vangysen, Herman↗
NASA NTRS1989

Modal identities for multibody elastic spacecraft: An aid to selecting modes for simulation

The question: Which set of modes furnishes a higher fidelity math model of dynamics of a multibody, deformable spacecraft (hinges-free or hinges-locked vehicle modes) is answered. Two sets of general, discretized, linear equations of motion of a spacecraft with an arbitrary number of deformable appendages, each articulated directly to the core body, are obtained using the above two families of modes. By a comparison of these equations, ten sets of modal identities are constructed which involve modal momenta coefficients and frequencies associated with both classes of modes. The sums of infinite series that appear in the identities are obtained in terms of mass, and first and second moments of inertia of the appendages, core body, and vehicle by using certain basic identities concerning appendage modes. Applying the above identities to a four-body spacecraft, the hinges-locked vehicle modes are found to yield a higher fidelity model than hinges-free modes, because the latter modes have nonconverging modal coefficients; a characteristic proved and illustrated.

Hablani, Hari B.↗
NASA NTRS1988

A Strassen-Newton algorithm for high-speed parallelizable matrix inversion

Techniques are described for computing matrix inverses by algorithms that are highly suited to massively parallel computation. The techniques are based on an algorithm suggested by Strassen (1969). Variations of this scheme use matrix Newton iterations and other methods to improve the numerical stability while at the same time preserving a very high level of parallelism. One-processor Cray-2 implementations of these schemes range from one that is up to 55 percent faster than a conventional library routine to one that is slower than a library routine but achieves excellent numerical stability. The problem of computing the solution to a single set of linear equations is discussed, and it is shown that this problem can also be solved efficiently using these techniques.

Bailey, David H.↗
NASA NTRS1989

A numerical study of the interaction between unsteady free-stream disturbances and localized variations in surface geometry

A numerical study of the generation of Tollmien-Schlichting (T-S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite-difference and spectral methods. The nonlinear steady flow is of the viscous-inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier-Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T-S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T-S waves.

Bodonyi, R. J.↗
NASA NTRS1990

Effect of nose bluntness on boundary layer stability and transition

The effect of nose bluntness on boundary layer instability is studied theoretically for a Mach 8 flow past a 7 degree semivertex cone. The basic flow is computed by solving the parabolized Navier-Stokes equations. Linear stability analysis of the basic flow reveals that, with small amount of bluntness, the critical Reynolds number for the onset of instability increases by an order of magnitude compared to the sharp cone value. The computed second mode frequencies are also in reasonable agreement with the experimental results. The results are used to explain the effect of unit Reynolds number on transition present in the quiet aeroballistic range data.

Malik, M. R.↗
NASA NTRS1990

Theoretical study of turbulent channel flow - Bulk properties, pressure fluctuations, and propagation of electromagnetic waves

In this paper, we apply two theoretical turbulence models, DIA and the recent GISS model, to study properties of a turbulent channel flow. Both models provide a turbulent kinetic energy spectral function E(k) as the solution of a non-linear equation; the two models employ the same source function but different closures. The source function is characterized by a rate n sub s (k) which is derived from the complex eigenvalues of the Orr-Sommerfeld (OS) equation in which the basic flow is taken to be of a Poiseuille type. The O-S equation is solved for a variety of Reynolds numbers corresponding to available experimental data. A physical argument is presented whereby the central line velocity characterizing the basic flow, U0 sup L, is not to be identified with the U0 appearing in the experimental Reynolds number. The theoretical results are compared with two types of experimental data: (1) turbulence bulk properties, and (2) properties that depend strongly on the structure of the turbulence spectrum at low wave numbers. The only existing analytical expression for Pi (k) cannot be used in the present case because it applies to the case of a flat plate, not a finite channel.

Canuto, V. M.↗
NASA NTRS1990

A parallel-vector algorithm for rapid structural analysis on high-performance computers

A fast, accurate Choleski method for the solution of symmetric systems of linear equations is presented. This direct method is based on a variable-band storage scheme and takes advantage of column heights to reduce the number of operations in the Choleski factorization. The method employs parallel computation in the outermost DO-loop and vector computation via the 'loop unrolling' technique in the innermost DO-loop. The method avoids computations with zeros outside the column heights, and as an option, zeros inside the band. The close relationship between Choleski and Gauss elimination methods is examined. The minor changes required to convert the Choleski code to a Gauss code to solve non-positive-definite symmetric systems of equations are identified. The results for two large-scale structural analyses performed on supercomputers, demonstrate the accuracy and speed of the method.

Storaasli, Olaf O.↗
NASA NTRS1990

The pdf approach to turbulent flow

This paper provides a detailed discussion of the theory and application of probability density function (pdf) methods, which provide a complete statistical description of turbulent flow fields at a single point or a finite number of points. The basic laws governing the flow of Newtonian fluids are set up in the Eulerian and the Lagrangian frame, and the exact and linear equations for the characteristic functionals in those frames are discussed. Pdf equations in both frames are derived as Fourier transforms of the equations of the characteristic functions. Possible formulations for the nonclosed terms in the pdf equation are discussed, their properties are assessed, and closure modes for the molecular-transport and the fluctuating pressure-gradient terms are reviewed. The application of pdf methods to turbulent combustion flows, supersonic flows, and the interaction of turbulence with shock waves is discussed.

Kollmann, W.↗
NASA NTRS1989

Parallel solution of pentadiagonal systems using generalized odd-even elimination

A method for the solution of pentadiagonal systems of linear equations is presented. The method is a generalization of ordinary odd-even elimination used for tridiagonal systems. Using n processors, an n x n pentadiagonal system can be solved using the new method (generalized odd-even elimination) in time proportional to log(2) n.

Levit, Creon↗
NASA NTRS1990

Optimal orbital rendezvous using high and low thrust

Optimal control theory is used to examine a specific class of spacecraft trajectory problems where high- and low-thrust propulsion systems are utilized. These problems assume a spacecraft is in an established orbit about a planet. It is desired to execute an intercept of a pre-determined position in space in a specified amount of time using an optimal high-thrust program. The spacecraft then returns to the original orbit station using the low-thrust propulsion system in an optimal fashion. A minimum fuel solution is sought using the linearized equations of motion, known as the CW equations, which simplify the necessary computations. Solutions are obtained for problems with a fixed final time. However, for the time-open case, the optimal solution is for the final time to be infinite. With a weighted function of the final time in the performance index, a limited range of optimal single impulse solutions for the time -open case can also be found.

Prussing, John E.↗
NASA NTRS1990

Predicted torque equilibrium attitude utilization for Space Station attitude control

An approximate knowledge of the torque equilibrium attitude (TEA) is shown to improve the performance of a control moment gyroscope (CMG) momentum management/attitude control law for Space Station Freedom. The linearized equations of motion are used in conjunction with a state transformation to obtain a control law which uses full state feedback and the predicted TEA to minimize both attitude excursions and CMG peak and secular momentum. The TEA can be computationally determined either by observing the steady state attitude of a 'controlled' spacecraft using arbitrary initial attitude, or by simulating a fixed attitude spacecraft flying in desired orbit subject to realistic environmental disturbance models.

Kumar, Renjith R.↗
NASA NTRS1990

Automated terminal guidance for a Shuttle rendezvous to Space Station Freedom

An automated terminal guidance for a Shuttle rendezvous based on the Clohessy-Wiltshire (CW) linear equations of relative motion is described. The algorithm guides the Shuttle from the last rendezvous midcourse maneuver (MC4) through docking with Space Station Freedom (SSF). The uniqueness of this algorithm is that it makes it possible to use the CW equations to fly a line-of-sight (LOS) Vbar or Rbar in the final-approach docking phase. The algorithm is made of two parts, in and out-of-plane, and can also be used for station keeping during final approaches. Simulation results of the guidance integrated with the Shuttle's flight control system in the Systems Engineering Simulator (SES) at NASA Johnson are discussed. Plans to add a laser radar docking sensor to the SES and integrate it with the Shuttle's rendezvous navigation are examined.

Olszewski, Oscar W.↗