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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Contribution of the National Aeronautics and Space Administration Langley Research Center

As part of a special international effort, three nozzles were designed and tested on single nacelle models in wind tunnels of several nations belonging to the North Atlantic Treaty Organization. All three of these nozzles were investigated in the Langley 16-foot transonic wind tunnel at the National Aeronautics and Space Administration's Langley Research Center. Langley Research Center also contributed theoretical calculations of the jet plume boundary and afterbody pressures. The calculations were obtained using an iterative solution which combined the inviscid Douglas Neumann method for the external flow with the method of characteristics for the flow in the jet plume. For the investigation, the nozzles were mounted on a single nacelle model 15.24 centimeters in diameter and 162.56 centimeters long. Tests were made at free stream Mach number from 0.4 to 1.2, and at Reynolds numbers per meter from 7.38 million to 13.78 million depending on the Mach number. Four types of data were recorded: afterbody pressure data, afterbody force data, model boundary layer data, and tunnel wall pressure data. The ratio of jet total pressure to free stream static pressure ranged up to 8.5. A description of the wind tunnel, model, and test procedure is included.

Compton, W. B., III↗
NASA NTRS1976

Electrostatic analysis of charge-coupled structures

An analysis is presented which is based upon a numerical solution of Gauss's law for the multiple dielectric semiconductor and electrode structure. The formulation resulted in a large set of equations, usually nonlinear, which were solved by an iterative scheme based upon the Gauss-Seidel method employing a relaxation parameter. Proper consideration of the nonlinearity of the space-charge in the semiconductor allowed convergent solutions which can give reasonable approximations to the channel potential profiles though not as good for estimating the space charge itself. The results of this work, which are three computer programs, are listed. None of the programs give the electric field explicitly as output data. However, when this information is desired, the programs can be expanded to include field estimations based on a polynominal fit to the potential field, or they may serve as subprograms for main programs dealing with general aspects of charge transfer.

Gassaway, J. D.↗
NASA NTRS1976

Numerical computation of transonic flows with shock waves

Some recent developments in numerical methods for calculating solutions to the transonic potential flow equation are reviewed, including (1) the construction of the stable coordinate independent difference schemes; (2) the use of conservation form to insure proper shock jump conditions; (3) analysis of the relaxation method by the time dependent analogy; (4) accelerated iterative schemes.

Jameson, A.↗
NASA NTRS1978

Dynamic programming and direct interaction for the optimum design of skeletal towers

A computer technique is proposed for automatically designing tower structures. Dynamic programming was used to find the optimum geometric configuration of the structural members, while the member sizes were proportioned by direct iteration. Tower structures are particularly suited to this method of automatic design since the rapidity of the analysis and design depends primarily upon substructuring. Substructuring of towers was comparatively simple because interaction between adjacent substructures is simulated with reasonable accuracy. Typical examples are presented to illustrate the method.

Howell, G. C.↗
NASA NTRS1979

The design of supercritical wings by the use of three-dimensional transonic theory

A procedure was developed for the design of transonic wings by the iterative use of three dimensional, inviscid, transonic analysis methods. The procedure was based on simple principles of supersonic flow and provided the designer with a set of guidelines for the systematic alteration of wing profile shapes to achieve some desired pressure distribution. The method was generally applicable to wing design at conditions involving a large region of supercriterical flow. To illustrate the method, it was applied to the design of a wing for a supercritical maneuvering fighter that operates at high lift and transonic Mach number. The wing profiles were altered to produce a large region of supercritical flow which was terminated by a weak shock wave. The spanwise variation of drag of this wing and some principles for selecting the streamwise pressure distribution are also discussed.

Mann, M. J.↗
NASA NTRS1982

Profiles of electrified drops and bubbles

Axisymmetric equilibrium shapes of conducting drops and bubbles, (1) pendant or sessile on one face of a circular parallel-plate capacitor or (2) free and surface-charged, are found by solving simultaneously the free boundary problem consisting of the augmented Young-Laplace equation for surface shape and the Laplace equation for electrostatic field, given the surface potential. The problem is nonlinear and the method is a finite element algorithm employing Newton iteration, a modified frontal solver, and triangular as well as quadrilateral tessellations of the domain exterior to the drop in order to facilitate refined analysis of sharply curved drop tips seen in experiments. The stability limit predicted by this computer-aided theoretical analysis agrees well with experiments.

Basaran, O. A.↗
NASA NTRS1983

A time-dependent approach for calculating steady inverse boundary-layer flows with separation

An unsteady inverse boundary-layer method is developed which can be used to calculate steady flows with separation. Two versions of Keller's box method with the Mechul function formulation developed by Cebeci (1976) are employed, depending on the complexity of the flow. The regular box is employed in regions of positive streamwise velocity component u, whereas the zigzag box is employed in regions where u becomes negative (t greater than 0). The regular box with the FLARE approximation is employed when t = 0 and u becomes negative in some region across the layer. Results of calculations show that the use of a time-dependent inverse boundary-layer method in which time is used as an iteration parameter provides a good approach in improving the accuracy of the solutions obtained from the FLARE approximation.

Cebeci, T.↗
NASA NTRS1985

Nonlinear transient responses of initially stressed composite plates

The nonlinear transient response of initially stressed composite plates is investigated using the finite element method. A nine-node isoparametric quadrilateral element is developed to model laminated plates under initial deformation and initial stress according to the Mindlin plate theory and von Karman large deflection assumptions. In the time integration, the Newmark constant acceleration method in conjunction with an efficient and accurate iteration scheme is used. Numerical results for deflections and bending moments for isotropic and laminated plates are obtained.

Chen, J. K.↗
NASA NTRS1987

Vortex breakdown simulation

In this paper, steady, axisymmetric inviscid, and viscous (laminar) swirling flows representing vortex breakdown phenomena are simulated using a stream function-vorticity-circulation formulation and two numerical methods. The first is based on an inverse iteration, where a norm of the solution is prescribed and the swirling parameter is calculated as a part of the output. The second is based on direct Newton iterations, where the linearized equations, for all the unknowns, are solved simultaneously by an efficient banded Gaussian elimination procedure. Several numerical solutions for inviscid and viscous flows are demonstrated, followed by a discussion of the results. Some improvements on previous work have been achieved: first order upwind differences are replaced by second order schemes, line relaxation procedure (with linear convergence rate) is replaced by Newton's iterations (which converge quadratically), and Reynolds numbers are extended from 200 up to 1000.

Hafez, M.↗
NASA NTRS1989

Inclusion of filter modulation in synthetic-discriminant-function construction

A technique in which the filter modulation is included in the synthesis of a synthetic-discriminant-function (SDF) matched spatial filter is presented. In the filter synthesis, a system of simultaneous nonlinear equations is solved with an iteration procedure. A computer simulation of the new method using thresholded images of the Space Shuttle over a range of aspect angles was performed for phase-only filters (POFs) and binary-phase-only filters (BPOFs). The filters constructed are capable of obtaining the specified peak-correlation response to within 1 percent with a high signal-to-clutter-ratio for the one-class problem, the two-class problem, and the multilevel problem. In contrast, conventional projection SDF POFs and BPOFs are unable to produce the desired peak-correlation response.

Jared, David A.↗
NASA NTRS1990

Monitoring of contact stresses in advanced propulsion systems

This paper discusses a new and innovative numerical procedure for evaluating and monitoring contact stresses in the mating surfaces of propulsion systems. The procedure is intended to extend the range of applicability and the accuracy of contact stress analyses. Originally proposed by Burton Paul, et. al., (1981) the procedure is based upon fundamental solutions from the theory of elasticity. It is an iterative numerical procedure. It is believed that the method has distinct advantages over the classical Hertz method, the finite-element method and over existing approaches with the boundary element method. Several examples are discussed. Favorable comparisons are obtained with classical results. Applications with mating element of propulsion systems are discussed.

Somprakit, Paisan↗
NASA NTRS1992

A new scheme for multidimensional line transfer. II - ETLA method in one dimension with application to iron K-alpha lines

The 1D version of the radiation hydrodynamics computer code ALTAIR is discribed. It is an implementation of the equivalent two-level atom method, and also incorporates a variety of other iterative techniques to achieve a fast, accurate solution of the problem of coupled atomic kinetics and radiative transfer. All the techniques are described in sufficient detail to enable their use by others. The application of the program is illustrated by a set of calculations of the iron X-ray emission spectrum from intercloud hot gas that may exist in the broad-line region of active galactic nuclei. These calculations demonstrate that such spectra may be a rich load to mine for information about the hot gas in AGNs, and also that the concept of 'the iron K-alpha line' inadequately depicts the spectrum.

Castor, John I.↗
NASA NTRS1992

Neural network error correction for solving coupled ordinary differential equations

A neural network is presented to learn errors generated by a numerical algorithm for solving coupled nonlinear differential equations. The method is based on using a neural network to correctly learn the error generated by, for example, Runge-Kutta on a model molecular dynamics (MD) problem. The neural network programs used in this study were developed by NASA. Comparisons are made for training the neural network using backpropagation and a new method which was found to converge with fewer iterations. The neural net programs, the MD model and the calculations are discussed.

Shelton, R. O.↗
NASA NTRS2024

Iterated Sigma-Point Kalman Filtering for Trajectory Reconstruction

Two second-order Kalman filtering methods are added to the New Statistical Trajec-tory Estimation Program (NewSTEP). The two filters are of the class of Sigma-Point Kalman Filters, namely the Unscented Kalman Filter and the Divided Difference Filter. Both filters are implemented with an iterative measurement update step. Results from simulation test case and a flight data test case are provided to show the filter performance. The results indicate that the new second-order filters pro-duce trajectory estimates with reduced uncertainties when compared to the existing first-order filtering method.

Sophia Vedvik↗
NASA NTRS1972

Finite element analysis of elasto-plastic soils. Report no. 4: Finite element analysis of elasto-plastic frictional materials for application to lunar earth sciences

The behavior of finite element models employing different constitutive relations to describe the stress-strain behavior of soils is investigated. Three models, which assume small strain theory is applicable, include a nondilatant, a dilatant and a strain hardening constitutive relation. Two models are formulated using large strain theory and include a hyperbolic and a Tresca elastic perfectly plastic constitutive relation. These finite element models are used to analyze retaining walls and footings. Methods of improving the finite element solutions are investigated. For nonlinear problems better solutions can be obtained by using smaller load increment sizes and more iterations per load increment than by increasing the number of elements. Suitable methods of treating tension stresses and stresses which exceed the yield criteria are discussed.

Marr, W. A., Jr.↗
NASA NTRS1980

Studies of leading-edge thrust phenomena

A study of practical limitations on achievement of theoretical leading-edge thrust has been made and an empirical method for estimation of attainable thrust has been developed. The method is based on a theoretical analysis of a set of two-dimensional airfoils to define thrust dependence on airfoil geometric characteristics and arbitrarily defined limiting pressures, an examination of two-dimensional airfoil experimental data to provide an estimate of limiting pressure dependence on local Mach number and Reynolds number, and employment of simple sweep theory to adapt the method to three-dimensional wings. Because the method takes into account the spanwise variation of airfoil section characteristics, an opportunity is afforded for design by iteration to maximize the attainable thrust and the attendant performance benefits. The applicability of the method was demonstrated by comparisons of theoretical and experimental aerodynamic characteristics for a series of wing-body configurations. Generally, good predictions of the attainable thrust and its influence on lift and drag characteristics were obtained over a range of Mach numbers from 0.24 to 2.0.

Carlson, H. W.↗
NASA NTRS1990

A finite element conjugate gradient FFT method for scattering

An extension of a two dimensional formulation is presented for a three dimensional body of revolution. With the introduction of a Fourier expansion of the vector electric and magnetic fields, a coupled two dimensional system is generated and solved via the finite element method. An exact boundary condition is employed to terminate the mesh and the fast fourier transformation (FFT) is used to evaluate the boundary integrals for low O(n) memory demand when an iterative solution algorithm is used. By virtue of the finite element method, the algorithm is applicable to structures of arbitrary material composition. Several improvements to the two dimensional algorithm are also described. These include: (1) modifications for terminating the mesh at circular boundaries without distorting the convolutionality of the boundary integrals; (2) the development of nonproprietary mesh generation routines for two dimensional applications; (3) the development of preprocessors for interfacing SDRC IDEAS with the main algorithm; and (4) the development of post-processing algorithms based on the public domain package GRAFIC to generate two and three dimensional gray level and color field maps.

Collins, Jeffery D.↗
NASA NTRS1991

A numerical algorithm for optimal feedback gains in high dimensional linear quadratic regulator problems

A hybrid method for computing the feedback gains in linear quadratic regulator problem is proposed. The method, which combines use of a Chandrasekhar type system with an iteration of the Newton-Kleinman form with variable acceleration parameter Smith schemes, is formulated to efficiently compute directly the feedback gains rather than solutions of an associated Riccati equation. The hybrid method is particularly appropriate when used with large dimensional systems such as those arising in approximating infinite-dimensional (distributed parameter) control systems (e.g., those governed by delay-differential and partial differential equations). Computational advantages of the proposed algorithm over the standard eigenvector (Potter, Laub-Schur) based techniques are discussed, and numerical evidence of the efficacy of these ideas is presented.

Banks, H. T.↗