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

Efficient algorithms for use in probabilistic finite element analysis

This paper investigates the use of Fast Probability Integration (FPI) algorithms in a Finite Element environment. A method allowing the representation of correlated fields in terms of a vector of uncorrelated transformed variables, based on the spectral decomposition of the variance-covariance matrix is developed. The response of the deterministic model corresponding to selected perturbations of these uncorrelated variables is then obtained via a Newton-type iterative scheme. The results of the perturbed problems are used to construct a local representation of the model's behavior in the neighborhood of the deterministic state, which the FPI algorithm will use to estimate the reliability of the system. Although the proposed strategy has thus far only been applied to linear elastostatics, the extension of the method to a broader class of problems appears to be feasible.

Dias, J. B.↗
NASA NTRS1994

Multidisciplinary design optimization using genetic algorithms

Multidisciplinary design optimization (MDO) is an important step in the conceptual design and evaluation of launch vehicles since it can have a significant impact on performance and life cycle cost. The objective is to search the system design space to determine values of design variables that optimize the performance characteristic subject to system constraints. Gradient-based optimization routines have been used extensively for aerospace design optimization. However, one limitation of gradient based optimizers is their need for gradient information. Therefore, design problems which include discrete variables can not be studied. Such problems are common in launch vehicle design. For example, the number of engines and material choices must be integer values or assume only a few discrete values. In this study, genetic algorithms are investigated as an approach to MDO problems involving discrete variables and discontinuous domains. Optimization by genetic algorithms (GA) uses a search procedure which is fundamentally different from those gradient based methods. Genetic algorithms seek to find good solutions in an efficient and timely manner rather than finding the best solution. GA are designed to mimic evolutionary selection. A population of candidate designs is evaluated at each iteration, and each individual's probability of reproduction (existence in the next generation) depends on its fitness value (related to the value of the objective function). Progress toward the optimum is achieved by the crossover and mutation operations. GA is attractive since it uses only objective function values in the search process, so gradient calculations are avoided. Hence, GA are able to deal with discrete variables. Studies report success in the use of GA for aircraft design optimization studies, trajectory analysis, space structure design and control systems design. In these studies reliable convergence was achieved, but the number of function evaluations was large compared with efficient gradient methods. Applicaiton of GA is underway for a cost optimization study for a launch-vehicle fuel-tank and structural design of a wing. The strengths and limitations of GA for launch vehicle design optimization is studied.

Unal, Resit↗
NASA NTRS1976

A numerical solution of the flow field over a transonic airfoil including strong-shock-induced flow separation

An approximate numerical method has been developed for computing the entire transonic flow field over an airfoil, including the boundary layer, the shock-wave/boundary-layer interaction, and the shock-induced-separation. The inviscid flow, the turbulent boundary layer in the weak interaction zone, and the strong interaction elliptic region are solved separately. The inviscid flow is represented by the solution of the transonic small-disturbance equation, while the boundary layer is solved by the generalized Galerkin method, and the strong interaction region (including the wake) is computed from the time-averaged Navier-Stokes equations. A second-order interaction is computed iteratively using the inviscid pressure distribution as a boundary condition for the boundary layer and Navier-Stokes codes, and the viscous displacement thickness to modify the airfoil geometry for the inviscid computation. The computed pressure distributions and shock-wave profiles are compared with experimental data. Agreement is very good except that the pressure recovery on the trailing edge is greater than the experimental data.

Seginer, A.↗
NASA NTRS1976

Computer-aided methods for analysis and synthesis of supersonic cruise aircraft structures

Computer-aided methods are reviewed which are being developed by Langley Research Center in-house work and by related grants and contracts. Synthesis methods to size structural members to meet strength and stiffness (flutter) requirements are emphasized and described. Because of the strong interaction among the aerodynamic loads, structural stiffness, and member sizes of supersonic cruise aircraft structures, these methods are combined into systems of computer programs to perform design studies. The approaches used in organizing these systems to provide efficiency, flexibility of use in an iterative process, and ease of system modification are discussed.

Giles, G. L.↗
NASA NTRS1980

Supersonic flow development in slotted wind tunnels

The development of test section slot shapes for achieving smooth supersonic Mach number distribution without overexpansion or waviness has, in the past, been largely an experimentally iterative or 'cutand try' procedure for each wind tunnel. To overcome the obvious disadvantages of time and expense involved in such an experimental approach, a simple analytical method has been developed to predict the supersonic flow development in a two-dimensional slotted tunnel given only the variation of open area ratio with downstream distance and the Mach number corresponding to the plenum static pressure. The well known method of characteristics is used with the constraint that it be compatible with the quadratic cross-flow pressure drop boundary condition at the slotted wall. The predicted results from this method agree remarkably well with the experimental calibration data available for some of the existing facilities. The flow mechanism responsible for causing overexpansion in the centerline Mach number distribution with some slot shapes has been brought to light.

Ramaswamy, M. A.↗
NASA NTRS1986

Separated flows near the nose of a body of revolution

The solution of the Navier-Stokes equations for the problem of cross-flow separataion about a deforming cylinder was achieved by iteration. It was shown that the separation starts at the rear stagnation point and the point of primary separation moves upstram along the cylinder surface. A general method of linear stability analysis for nonparallel external flows was constructed, which consists of representing the eigenfunctions with complete orthogonal sets and forms characteristic equations with the Galerkin method. The method was applied to the Kovasznay flow which is an exact solution of the Navier-Stokes equation. The results show that when the critical parameter is exceeded, there are only a few isolated unstable eigen-frequencies. Another exact solution is shown to be absolutely and monotonically stable with respect to infinitesimal disturbances of all frequencies. The flow is also globally, asymptotically, and monotonically stable in the mean with respect o three-dimensional disturbances. This result forms the sound foundation of rigorous stability analysis for nonparallel flows, and provides an invaluable test ground for future studies of nonparallel flows in which the basic states do not posses exact solutions. The application of this method to the study of the formation of spiral vorticies near the nose of a rotating body of revolution is underway. The same method will be applied to the stability analysis of reversed flow over a plate with suction.

Lin, S. P.↗
NASA NTRS1989

Multi-grid for structures analysis

In structural analysis the amount of computational time necessary for a solution is proportional to the number of degrees of freedom times the bandwidth squared. In implicit time analysis, this must be done at each discrete point in time. If, in addition, the problem is nonlinear, then this solution must be iterated at each point in time. If the bandwidth is large, the size of the problem that can be analyzed is severely limited. The multi-grid method is a possible algorithm which can make this solution much more computationally efficient. This method has been used for years in computational fluid mechanics. It works on the fact that relaxation is very efficient on the high frequency components of the solution (nearest neighbor interactions) and not very good on low frequency components of the solution (far interactions). The multi-grid method is then to relax the solution on a particular model until the residual stops changing. This indicates that the solution contains the higher frequency components. A coarse model is then generated for the lower frequency components to the solution. The model is then relaxed for the lower frequency components of the solution. These lower frequency components are then interpolated to the fine model. In computational fluid mechanics the equations are usually expressed as finite differences.

Kascak, Albert F.↗
NASA NTRS1991

Calculating satellite umbra/penumbra entry and exit positions and times

The problem of calculating earth satellite entry/exit positions and times through the earth's umbra and penumbra for elliptical satellite orbits is solved in closed form; i.e., without iteration, by reducing the problem to that of finding the roots of a quartic polynomial in a cartesian coordinate. Typical results from an algorithm constructed from the method are shown.

Mullins, Larry D.↗
NASA NTRS1955

Theoretical Prediction of Pressure Distributions on Nonlifting Airfoils at High Subsonic Speeds

Theoretical pressure distributions on nonlifting circular-arc airfoils in two-dimensional flows with high subsonic free-stream velocity are found by determining approximate solutions, through an iteration process, of an integral equation for transonic flow proposed by Oswatitsch. The integral equation stems directly from the small-disturbance theory for transonic flow. This method of analysis possesses the advantage of remaining in the physical, rather than the hodograph, variable and can be applied in airfoils having curved surfaces. After discussion of the derivation of the integral equation and qualitative aspects of the solution, results of calculations carried out for circular-arc airfoils in flows with free-stream Mach numbers up to unity are described. These results indicate most of the principal phenomena observed in experimental studies.

Spreiter, John R↗
NASA NTRS1994

Recovering the real density field of galaxies from redshift space

This paper proposes a number of methods for recovering the true, real-space large-scale density distributions from redshift data sets. On the largest scales, the effects of peculiar velocities is to compress the galaxy clustering along the line of sight and consequently the density gradient in the radial direction is enhanced in comparison with the density gradient perpendicular to the line of sight. We calculate the relation between the density gradients in real and redshift space in the linear regime and show that it is a simple function of the density parameter, Omega. By considering the anisotropy of density gradients in redshift space it is possible, in principle, to obtain information about the value of Omega. Using N-body simulations we examine different methods for correcting the effects of large-scale velocities to obtain the true density distributions. We derive the peculiar velocity field from the redshift data and iterate to place each galaxy at a distance consistent with its measured redshift and its derived velocity. Linear theory, as well as second-order solutions for the peculiar velocity field are tested. The density distributions in low-density regions predicted by both linear theory and second-order theory are similar, and in good agreement with the true density field in the simulation. The use of the second-order corrections improves the prediction of the density field in high-density regions.

Gramann, Mirt↗
NASA NTRS1997

Aerodynamic Shape Sensitivity Analysis and Design Optimization of Complex Configurations Using Unstructured Grids

A three-dimensional unstructured grid approach to aerodynamic shape sensitivity analysis and design optimization has been developed and is extended to model geometrically complex configurations. The advantage of unstructured grids (when compared with a structured-grid approach) is their inherent ability to discretize irregularly shaped domains with greater efficiency and less effort. Hence, this approach is ideally suited for geometrically complex configurations of practical interest. In this work the nonlinear Euler equations are solved using an upwind, cell-centered, finite-volume scheme. The discrete, linearized systems which result from this scheme are solved iteratively by a preconditioned conjugate-gradient-like algorithm known as GMRES for the two-dimensional geometry and a Gauss-Seidel algorithm for the three-dimensional; similar procedures are used to solve the accompanying linear aerodynamic sensitivity equations in incremental iterative form. As shown, this particular form of the sensitivity equation makes large-scale gradient-based aerodynamic optimization possible by taking advantage of memory efficient methods to construct exact Jacobian matrix-vector products. Simple parameterization techniques are utilized for demonstrative purposes. Once the surface has been deformed, the unstructured grid is adapted by considering the mesh as a system of interconnected springs. Grid sensitivities are obtained by differentiating the surface parameterization and the grid adaptation algorithms with ADIFOR (which is an advanced automatic-differentiation software tool). To demonstrate the ability of this procedure to analyze and design complex configurations of practical interest, the sensitivity analysis and shape optimization has been performed for a two-dimensional high-lift multielement airfoil and for a three-dimensional Boeing 747-200 aircraft.

Taylor, Arthur C., III↗
NASA NTRS2004

Procedure for Determining One-Dimensional Flow Distributions in Arbitrarily Connected Passages Without the Influence of Pumping

A calculation procedure is presented which allows the one-dimensional determination of flow distributions in arbitrarily connected (branching) flow passages having multiple inlets and exits. The procedure uses an adaptation of the finite element technique, iteratively coupled with an accurate one-dimensional flow solver. The procedure eliminates the usual restrictions inherent with finite element flow calculations. Unlike existing one-dimensional methods, which require simplifications to the flow equations (uncoupling the momentum and energy equations), to allow for arbitrary branching and multiple inlets and exits, the only limitation of the described methodology is that, at present, it can only accommodate non-rotating configurations (no pumping effects). The calculation procedure is robust, and will always converge for physically possible flow. The procedure is described, and its use is illustrated by an example.

Meitner, Peter L.↗
NASA NTRS1998

A Numerical Optimization Approach for Tuning Fuzzy Logic Controllers

This paper develops a method to tune fuzzy controllers using numerical optimization. The main attribute of this approach is that it allows fuzzy logic controllers to be tuned to achieve global performance requirements. Furthermore, this approach allows design constraints to be implemented during the tuning process. The method tunes the controller by parameterizing the membership functions for error, change-in-error and control output. The resulting parameters form a design vector which is iteratively changed to minimize an objective function. The minimal objective function results in an optimal performance of the system. A spacecraft mounted science instrument line-of-sight pointing control is used to demonstrate results.

Woodard, Stanley E.↗
NASA NTRS2009

Three Dimensional Sector Design with Optimal Number of Sectors

The concept of dynamic sector design suggests a strategic approach to ease air traffic congestion, which is predicted to become a serious problem in the national airspace system by 2025. Considerable research has been conducted to address the sectorization problem. In previous work, an approach that combines the Voronoi diagrams, Genetic Algorithms (GA), and the iterative deepening algorithm was proposed. However, as originally formulated, the number of sectors used was predefined and only two-dimensional partitions were allowed, which constrained the method's ability to achieve good designs. The current work extends the earlier Voronoi-based method by treating the number of sectors as an additional decision variable, allowing 3D partitions, and developing more comprehensive costs.

Xue, Min↗
NASA NTRS1975

Simultaneous calculation of aircraft design loads and structural member sizes

A design process which accounts for the interaction between aerodynamic loads and changes in member sizes during sizing of aircraft structures is described. A simultaneous iteration procedure is used wherein both design loads and member sizes are updated during each cycle yielding converged, compatible loads and member sizes. A description is also given of a system of programs which incorporates this process using lifting surface theory to calculate aerodynamic pressure distributions, using a finite-element method for structural analysis, and using a fully stressed design technique to size structural members. This system is tailored to perform the entire process with computational efficiency in a single computer run so that it can be used effectively during preliminary design. Selected results, considering maneuver, taxi, and fatigue design conditions, are presented to illustrate convergence characteristics of this iterative procedure.

Giles, G. L.↗
NASA NTRS1996

Iterative LQG Controller Design Through Closed-Loop Identification

This paper presents an iterative Linear Quadratic Gaussian (LQG) controller design approach for a linear stochastic system with an uncertain open-loop model and unknown noise statistics. This approach consists of closed-loop identification and controller redesign cycles. In each cycle, the closed-loop identification method is used to identify an open-loop model and a steady-state Kalman filter gain from closed-loop input/output test data obtained by using a feedback LQG controller designed from the previous cycle. Then the identified open-loop model is used to redesign the state feedback. The state feedback and the identified Kalman filter gain are used to form an updated LQC controller for the next cycle. This iterative process continues until the updated controller converges. The proposed controller design is demonstrated by numerical simulations and experiments on a highly unstable large-gap magnetic suspension system.

Hsiao, Min-Hung↗
NASA NTRS1985

Full two-dimensional transient solutions of electrothermal aircraft blade deicing

Two finite difference methods are presented for the analysis of transient, two-dimensional responses of an electrothermal de-icer pad of an aircraft wing or blade with attached variable ice layer thickness. Both models employ a Crank-Nicholson iterative scheme, and use an enthalpy formulation to handle the phase change in the ice layer. The first technique makes use of a 'staircase' approach, fitting the irregular ice boundary with square computational cells. The second technique uses a body fitted coordinate transform, and maps the exact shape of the irregular boundary into a rectangular body, with uniformally square computational cells. The numerical solution takes place in the transformed plane. Initial results accounting for variable ice layer thickness are presented. Details of planned de-icing tests at NASA-Lewis, which will provide empirical verification for the above two methods, are also presented.

Masiulaniec, K. C.↗