Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iteration method”

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,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

An efficient iteration strategy for the solution of the Euler equations

A line Gauss-Seidel (LGS) relaxation algorithm in conjunction with a one-parameter family of upwind discretizations of the Euler equations in two-dimensions is described. The basic algorithm has the property that convergence to the steady-state is quadratic for fully supersonic flows and linear otherwise. This is in contrast to the block ADI methods (either central or upwind differenced) and the upwind biased relaxation schemes, all of which converge linearly, independent of the flow regime. Moreover, the algorithm presented here is easily enhanced to detect regions of subsonic flow embedded in supersonic flow. This allows marching by lines in the supersonic regions, converging each line quadratically, and iterating in the subsonic regions, thus yielding a very efficient iteration strategy. Numerical results are presented for two-dimensional supersonic and transonic flows containing both oblique and normal shock waves which confirm the efficiency of the iteration strategy.

Walters, R. W.↗
NASA NTRS1993

Extending substructure based iterative solvers to multiple load and repeated analyses

Direct solvers currently dominate commercial finite element structural software, but do not scale well in the fine granularity regime targeted by emerging parallel processors. Substructure based iterative solvers--often called also domain decomposition algorithms--lend themselves better to parallel processing, but must overcome several obstacles before earning their place in general purpose structural analysis programs. One such obstacle is the solution of systems with many or repeated right hand sides. Such systems arise, for example, in multiple load static analyses and in implicit linear dynamics computations. Direct solvers are well-suited for these problems because after the system matrix has been factored, the multiple or repeated solutions can be obtained through relatively inexpensive forward and backward substitutions. On the other hand, iterative solvers in general are ill-suited for these problems because they often must restart from scratch for every different right hand side. In this paper, we present a methodology for extending the range of applications of domain decomposition methods to problems with multiple or repeated right hand sides. Basically, we formulate the overall problem as a series of minimization problems over K-orthogonal and supplementary subspaces, and tailor the preconditioned conjugate gradient algorithm to solve them efficiently. The resulting solution method is scalable, whereas direct factorization schemes and forward and backward substitution algorithms are not. We illustrate the proposed methodology with the solution of static and dynamic structural problems, and highlight its potential to outperform forward and backward substitutions on parallel computers. As an example, we show that for a linear structural dynamics problem with 11640 degrees of freedom, every time-step beyond time-step 15 is solved in a single iteration and consumes 1.0 second on a 32 processor iPSC-860 system; for the same problem and the same parallel processor, a pair of forward/backward substitutions at each step consumes 15.0 seconds.

Farhat, Charbel↗
NASA NTRS2005

Real-Time Principal-Component Analysis

A recently written computer program implements dominant-element-based gradient descent and dynamic initial learning rate (DOGEDYN), which was described in Method of Real-Time Principal-Component Analysis (NPO-40034) NASA Tech Briefs, Vol. 29, No. 1 (January 2005), page 59. To recapitulate: DOGEDYN is a method of sequential principal-component analysis (PCA) suitable for such applications as data compression and extraction of features from sets of data. In DOGEDYN, input data are represented as a sequence of vectors acquired at sampling times. The learning algorithm in DOGEDYN involves sequential extraction of principal vectors by means of a gradient descent in which only the dominant element is used at each iteration. Each iteration includes updating of elements of a weight matrix by amounts proportional to a dynamic initial learning rate chosen to increase the rate of convergence by compensating for the energy lost through the previous extraction of principal components. In comparison with a prior method of gradient-descent-based sequential PCA, DOGEDYN involves less computation and offers a greater rate of learning convergence. The sequential DOGEDYN computations require less memory than would parallel computations for the same purpose. The DOGEDYN software can be executed on a personal computer.

Duong, Vu↗
NASA NTRS1994

A Parallel Incompressible Navier-Stokes Solver with Multigrid Iterations

We developed a parallel, numerically accurate and stable, and computationally efficient finate-difference incompressible Navier-Stokes (N-S) fluid flow solver. The solver runs on both sequential and massively parallel computers. The numerical method used here is a second-order projection method on a staggered grid. The code is highly modular and it can be used either as a stand-alone flow solver and or a template code which can be adapted or expanded to a specific application. Numerical results and parallel performances of our code on Intel Delta and Paragon are reported.

Navier-Stokes solver↗
NASA NTRS1977

A study of the optimization method used in the NAVY/NASA gas turbine engine computer code

Sources of numerical noise affecting the convergence properties of the Powell's Principal Axis Method of Optimization in the NAVY/NASA gas turbine engine computer code were investigated. The principal noise source discovered resulted from loose input tolerances used in terminating iterations performed in subroutine CALCFX to satisfy specified control functions. A minor source of noise was found to be introduced by an insufficient number of digits in stored coefficients used by subroutine THERM in polynomial expressions of thermodynamic properties. Tabular results of several computer runs are presented to show the effects on program performance of selective corrective actions taken to reduce noise.

Horsewood, J. L.↗
NASA NTRS1979

A consistent design procedure for supercritical airfoils in free air and a wind tunnel

A computational inverse procedure for transonic airfoils in which shapes are determined supporting prescribed pressure distributions is presented. The method uses the small disturbance equation and a consistent analysis-design differencing procedure at the airfoil surface. This avoids the intermediate analysis-design-analysis iterations. The effect of any openness at the trailing edge is taken onto account by adding an effective source term in the far field. The final results from a systematic expansion procedure which models the far field for solid, ideal slotted, and free jet tunnel walls are presented along with some design results for the associated boundary conditions and those for a free flight.

Shankar, V.↗
NASA NTRS1982

A search for spectroscopic binaries among the runaway O type stars

Numerous radial velocity measurements of medium dispersion were made for the 10 brighter stars given in Stone's list of very probable O type runaways. All plates were measured with the KPNO PDS microdensitometer, and a new iterative reductional analysis was used to derive plate velocities, which are estimated to be 1.6 times more accurate internally than those found by using the traditional method. Of thse stars, psi Per, alpha Cam, HD 188209, and 26 Cep are identified as probable velocity variables, while 9 Sge, lambda Cep, and HD 218915 are classed as possible variables. If the source of this variability is Keplerian rather than atmospheric, which cannot be established unequivocally from the observations of this paper, psi Per could be a spectroscopic binary with a black hole companion, and at least 1.2 solar mass. The detection of runaway binary systems from radial velocity measurements is discussed.

Stone, R. C.↗
NASA NTRS1986

Numerical and physical aspects of aerodynamic flows III; Proceedings of the Third Symposium, California State University, Long Beach, CA, January 21-24, 1985

The present symposium considers hydrodynamic stability and turbulent transition, transition calculations in three-dimensional flows, a quasi-simultaneous finite difference approach for strongly interacting flows, the significance of the thin layer Navier-Stokes approximation, unsteady airfoil boundary layers, predictions and experiments on airfoils at low Reynolds numbers, and a comparison of interactive boundary layer and thin layer Navier-Stokes procedures. Also discussed are a viscous-inviscid interaction method for computing unsteady transonic separation, massive separation and dynamic stall on a cusped trailing edge airfoil, the computation of turbulent separated flows over wings, an iterative scheme for three-dimensional transonic flows, and the computation of three-dimensional flows with shock wave/boundary layer interaction. Many of the papers in this symposium were abstracted previously (cf., A85-42951).

Cebeci, T.↗
NASA NTRS1988

Calculation of aerodynamic characteristics of airplane configurations at high angles of attack

Calculation of longitudinal and lateral directional aerodynamic characteristics of airplanes by the VORSTAB code is examined. The numerical predictions are based on the potential flow theory with corrections of high angle of attack phenomena; namely, vortex flow and boundary layer separation effects. To account for the vortex flow effect, vortex lift, vortex action point, augmented vortex lift and vortex breakdown effect through the method of suction analogy are included. The effect of boundary layer separation is obtained by matching the nonlinear section data with the three dimensional lift characteristics iteratively. Through correlation with results for nine fighter configurations, it is concluded that reasonably accurate prediction of longitudinal and static lateral directional aerodynamics can be obtained with the VORSTAB code up to an angle of attack at which wake interference and forebody vortex effect are not important. Possible reasons for discrepancy at higher angles of attack are discussed.

Tseng, J. B.↗
NASA NTRS1987

Fast simulation techniques for switching converters

Techniques for simulating a switching converter are examined. The state equations for the equivalent circuits, which represent the switching converter, are presented and explained. The uses of the Newton-Raphson iteration, low ripple approximation, half-cycle symmetry, and discrete time equations to compute the interval durations are described. An example is presented in which these methods are illustrated by applying them to a parallel-loaded resonant inverter with three equivalent circuits for its continuous mode of operation.

King, Roger J.↗
NASA NTRS1989

Influence of impeller shroud forces on turbopump rotor dynamics

The shrouded-impeller leakage path forces calculated by Childs (1987) have been analyzed to answer two questions. First, because of certain characteristics of the results of Childs, the forces could not be modeled with traditional approaches. Therefore, an approach has been devised to include the forces in conventional rotordynamic analyses. The forces were approximated by traditional stiffness, damping and inertia coefficients with the addition of whirl-frequency-dependent direct and cross-coupled stiffness terms. The forces were found to be well-modeled with this approach. Finally, the effect these forces had on a simple rotor-bearing system was analyzed, and, therefore, they, in addition to seal forces, were applied to a Jeffcott rotor. The traditional methods of dynamic system analysis were modified to incorporate the impeller forces and yielded results for the eigenproblem, frequency response, critical speed, transient response and an iterative technique for finding the frequency of free vibration as well as system stability. All results lead to the conclusion that the forces have little influence on natural frequency but can have appreciable effects on system stability. Specifically, at higher values of fluid swirl at the leakage path entrance, relative stability is reduced. The only unexpected response characteristics that occurred are attributed to the nonlinearity of the model.

Williams, Jim P.↗
NASA NTRS1991

Concurrent optimization of airframe and engine design parameters

An integrated system for the multidisciplinary analysis and optimization of airframe and propulsion design parameters is being developed. This system is known as IPAS, the Integrated Propulsion/Airframe Analysis System. The traditional method of analysis is one in which the propulsion system analysis is loosely coupled to the overall mission performance analysis. This results in a time consuming iterative process. First, the engine is designed and analyzed. Then, the results from this analysis are used in a mission analysis to determine the overall aircraft performance. The results from the mission analysis are used as a guide as the engine is redesigned and the entire process repeated. In IPAS, the propulsion system, airframe, and mission are closely coupled. The propulsion system analysis code is directly integrated into the mission analysis code. This allows the propulsion design parameters to be optimized along with the airframe and mission design parameters, significantly reducing the time required to obtain an optimized solution.

Lavelle, Thomas M.↗
NASA NTRS1993

Nonlinear rescaling of control values simplifies fuzzy control

Traditional control theory is well-developed mainly for linear control situations. In non-linear cases there is no general method of generating a good control, so we have to rely on the ability of the experts (operators) to control them. If we want to automate their control, we must acquire their knowledge and translate it into a precise control strategy. The experts' knowledge is usually represented in non-numeric terms, namely, in terms of uncertain statements of the type 'if the obstacle is straight ahead, the distance to it is small, and the velocity of the car is medium, press the brakes hard'. Fuzzy control is a methodology that translates such statements into precise formulas for control. The necessary first step of this strategy consists of assigning membership functions to all the terms that the expert uses in his rules (in our sample phrase these words are 'small', 'medium', and 'hard'). The appropriate choice of a membership function can drastically improve the quality of a fuzzy control. In the simplest cases, we can take the functions whose domains have equally spaced endpoints. Because of that, many software packages for fuzzy control are based on this choice of membership functions. This choice is not very efficient in more complicated cases. Therefore, methods have been developed that use neural networks or generic algorithms to 'tune' membership functions. But this tuning takes lots of time (for example, several thousands iterations are typical for neural networks). In some cases there are evident physical reasons why equally space domains do not work: e.g., if the control variable u is always positive (i.e., if we control temperature in a reactor), then negative values (that are generated by equal spacing) simply make no sense. In this case it sounds reasonable to choose another scale u' = f(u) to represent u, so that equal spacing will work fine for u'. In the present paper we formulate the problem of finding the best rescaling function, solve this problem, and show (on a real-life example) that after an optimal rescaling, the un-tuned fuzzy control can be as good as the best state-of-art traditional non-linear controls.

Vanlangingham, H.↗
NASA NTRS1992

Concurrent optimization of airframe and engine design parameters

An integrated system for the multidisciplinary analysis and optimization of airframe and propulsion design parameters is being developed. This system is known as IPAS, the Integrated Propulsion/Airframe Analysis System. The traditional method of analysis is one in which the propulsion system analysis is loosely coupled to the overall mission performance analysis. This results in a time consuming iterative process. First, the engine is designed and analyzed. Then, the results from this analysis are used in a mission analysis to determine the overall aircraft performance. The results from the mission analysis are used as a guide as the engine redesigned and the entire process repeated. In IPAS, the propulsion system, airframe, and mission are closely coupled. The propulsion system analysis code is directly integrated into the mission analysis code. This allows the propulsion design parameters to be optimized along with the airframe and mission design parameters, significantly reducing the time required to obtain an optimized solution.

Lavelle, Thomas M.↗
NASA NTRS1993

Influence of impeller shroud forces on turbopump rotor dynamics

The shrouded-impeller leakage path forces calculated by Childs have been analyzed to answer two questions. First, because of certain characteristics or the results of Childs, the forces could not be modeled with traditional approaches. Therefore, an approach has been devised to include the forces in conventional rotordynamic analyses. The forces were found to be well-modeled with this approach. Finally, the effect these forces had on a simple rotor-bearing system was analyzed, and, therefore, they, in addition to seal forces, were applied to a Jeffcott rotor. The traditional methods of dynamic system analysis were modified to incorporate the impeller forces and yielded results for the eigenproblem, frequency response, critical speed, transient response, and an iterative technique for finding the frequency of free vibration as well as system stability. All results lead to the conclusion that the forces have little influence on natural frequency but can have appreciable effects on system stability. Specifically, at higher values of fluid swirl at the leakage path entrance, relative stability is reduced. The only unexpected response characteristics that occurred are attributed to the nonlinearity of the model.

Williams, J. P.↗
NASA NTRS1987

Procedures for Computing Transonic Flows for Control of Adaptive Wind Tunnels

Numerical methods are developed for control of three dimensional adaptive test sections. The physical properties of the design problem occurring in the external field computation are analyzed, and a design procedure suited for solution of the problem is worked out. To do this, the desired wall shape is determined by stepwise modification of an initial contour. The necessary changes in geometry are determined with the aid of a panel procedure, or, with incident flow near the sonic range, with a transonic small perturbation (TSP) procedure. The designed wall shape, together with the wall deflections set during the tunnel run, are the input to a newly derived one-step formula which immediately yields the adapted wall contour. This is particularly important since the classical iterative adaptation scheme is shown to converge poorly for 3D flows. Experimental results obtained in the adaptive test section with eight flexible walls are presented to demonstrate the potential of the procedure. Finally, a method is described to minimize wall interference in 3D flows by adapting only the top and bottom wind tunnel walls.

Rebstock, Rainer↗
NASA NTRS1986

A numerical method for three-dimensional incompressible flows using nonorthogonal body-fitted coordinate systems

In the present paper, a numerical method for solving the equations of motion of three-dimensional incompressible flows in nonorthogonal body-fitted coordinate (BEC) systems has been developed and evaluated. The equations of motion are transformed to a generalized curvilinear coordinate system from which the transformed equations are discretized using finite difference approximations in the transformed domain. The hybrid scheme and a central differencing plus artificial dissipation scheme are used to approximate the convection terms in the governing equations. Effects of these two schemes on the accuracy of numerical predictions are studied. Solutions of the finite difference equations are obtained iteratively by using a pressure-velocity correction algorithm, SIMPLE-C. Numerical examples of two- and three-dimensional, laminar and turbulent flow problems are employed to evaluate the accuracy and efficiency of the present numerical method.

Chen, Y. S.↗