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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Robust Airfoil Optimization to Achieve Consistent Drag Reduction Over a Mach Range

We prove mathematically that in order to avoid point-optimization at the sampled design points for multipoint airfoil optimization, the number of design points must be greater than the number of free-design variables. To overcome point-optimization at the sampled design points, a robust airfoil optimization method (called the profile optimization method) is developed and analyzed. This optimization method aims at a consistent drag reduction over a given Mach range and has three advantages: (a) it prevents severe degradation in the off-design performance by using a smart descent direction in each optimization iteration, (b) there is no random airfoil shape distortion for any iterate it generates, and (c) it allows a designer to make a trade-off between a truly optimized airfoil and the amount of computing time consumed. For illustration purposes, we use the profile optimization method to solve a lift-constrained drag minimization problem for 2-D airfoil in Euler flow with 20 free-design variables. A comparison with other airfoil optimization methods is also included.

Li, Wu↗
NASA NTRS1993

Interactive method for planning constrained, fuel-optimal orbital proximity operations

An interactive graphical method for planning fuel-efficient rendezvous trajectories in the multi-spacecraft environment of the space station is presented. The method allows the operator to compose a multi-burn transfer trajectory between arbitrary initial chaser and target trajectories. The available task time of the mission is limited and the maneuver is subject to various operational constraints, such as departure, arrival, plume impingement and spatial constraints. The maneuvers are described in terms of the relataive motion experienced in a Space-Station centered coordinate system. The optimization method is based on the primer vector and its extension to non-optimal trajectories. The visual feedback of trajectory shapes, operational constraints, and optimization functions, provided by user-transparaent and continuously active background computations, allows the operator to make fast, iterative design changes which rapidly converge to fuel-efficient solutions. The optimization functions are presented. A variety of simple design examples has been presented to demonstrate the usefulness of the method. In many cases the addition of a properly positioned intermediate waypoint resulted in fuel savings of up to 30%. Furthermore, due to the counter-intuitive character of the optimization functions, most fuel-optimal solutions could not have been found without the aid of the optimization tools. Operating the system was found to be very easy, and did not require any previous in-depth knowledge of orbital dynamics or trajectory. The planning tool is an example of operator assisted optimization of nonlinear cost-functions.

Abramovitz, Adrian↗
NASA NTRS2004

Above-Cloud Precipitable Water Retrievals using the MODIS 0.94 micron Band with Applications for Multi-Layer Cloud Detection

In order to perform satellite retrievals of cloud properties, it is important to account for the effect of the above-cloud atmosphere on the observations. The solar bands used in the operational MODIS Terra and Aqua cloud optical and microphysical algorithms (visible, NIR, and SWIR spectral windows) are primarily affected by water vapor, and to a lesser extent by well-mixed gases. For water vapor, the above-cloud column amount, or precipitable water, provides adequate information for an atmospheric correction; details of the vertical vapor distribution are not typically necessary for the level of correction required. Cloud-top pressure has a secondary effect due to pressure broadening influences. For well- mixed gases, cloud-top pressure is also required for estimates of above-cloud abundances. We present a method for obtaining above-cloud precipitable water over dark Ocean surfaces using the MODIS 0.94 pm vapor absorption band. The retrieval includes an iterative procedure for establishing cloud-top temperature and pressure, and is useful for both single layer water and ice clouds. Knowledge of cloud thermodynamic phase is fundamental in retrieving cloud optical and microphysical properties. However, in cases of optically thin cirrus overlapping lower water clouds, the concept of a single unique phase is ill- defined and depends, at least, on the spectral region of interest. We will present a method for multi-layer and multi-phase cloud detection which uses above-cloud precipitable water retrievals along with several existing MODIS operational cloud products (cloud-top pressure derived from a C02 slicing algorithm, IR and SWIR phase retrievals). Results are catagorized by whether the radiative signature in the MODIS solar bands is primarily that of a water cloud with ice cloud contamination, or visa-versa. Examples in polar and mid-latitude regions will be shown.

Platnick, S.↗
NASA NTRS1991

A combined finite element-boundary element formulation for solution of axially symmetric bodies

A new method is presented for the computation of electromagnetic scattering from axially symmetric bodies. To allow the simulation of inhomogeneous cross sections, the method combines the finite element and boundary element techniques. Interior to a fictitious surface enclosing the scattering body, the finite element method is used which results in a sparce submatrix, whereas along the enclosure the Stratton-Chu integral equation is enforced. By choosing the fictitious enclosure to be a right circular cylinder, most of the resulting boundary integrals are convolutional and may therefore be evaluated via the FFT with which the system is iteratively solved. In view of the sparce matrix associated with the interior fields, this reduces the storage requirement of the entire system to O(N) making the method attractive for large scale computations. The details of the corresponding formulation and its numerical implementation are described.

Collins, Jeffrey D.↗
NASA NTRS1985

Recent Developments in Nonlinear Normal Mode Initialization

The importance of a balanced initial condition upon GLAS GCM forecasts and assimilation cycles was assessed. An effort to combine previous work on normal mode initialization at GLA is underway to develop an initialization process for the production version of the GLAS 4th order GCM. The major aspects of this work fall into two parts: vectorization of the linear projector code and the insertion of the mode projector and Machenhauer iteration algorithm into the full GLAS GCM. Memory and paging constraints place restrictions on the number of horizontal modes stored for initialization purposes, and on the manner in which they are stored. Only the first five vertical structures of the gravity modes are used. Differing phase and normalization conventions provided many elusive coding errors. A Machenhauer nonlinear normal mode initialization technique is used. This method entails the insertion of a modified version of the mode projector into the full GCM, and the modification of the GCM to allow for iterative calls to the projector.

Bloom, S. C.↗
NASA NTRS2004

Bayesian Vision for Shape Recovery

We present a new Bayesian vision technique that aims at recovering a shape from two or more noisy observations taken under similar lighting conditions. The shape is parametrized by a piecewise linear height field, textured by a piecewise linear irradiance field, and we assume Gaussian Markovian priors for both shape vertices and irradiance variables. The observation process. also known as rendering, is modeled by a non-affine projection (e.g. perspective projection) followed by a convolution with a piecewise linear point spread function. and contamination by additive Gaussian noise. We assume that the observation parameters are calibrated beforehand. The major novelty of the proposed method consists of marginalizing out the irradiances considered as nuisance parameters, which is achieved by Laplace approximations. This reduces the inference to minimizing an energy that only depends on the shape vertices, and therefore allows an efficient Iterated Conditional Mode (ICM) optimization scheme to be implemented. A Gaussian approximation of the posterior shape density is computed, thus providing estimates both the geometry and its uncertainty. We illustrate the effectiveness of the new method by shape reconstruction results in a 2D case. A 3D version is currently under development and aims at recovering a surface from multiple images, reconstructing the topography by marginalizing out both albedo and shading.

Jalobeanu, Andre↗
NASA NTRS1971

Digital simulation error curves for a spring-mass-damper system

Plotting digital simulation errors for a spring-mass-damper system and using these error curves to select type of integration, feedback update method, and number of samples per cycle at resonance reduces excessive number of samples per cycle and unnecessary iterations.

Knox, L. A.↗
NASA NTRS1972

Computation of optimal output-feedback compensators for linear time-invariant systems

The control of linear time-invariant systems with respect to a quadratic performance criterion was considered, subject to the constraint that the control vector be a constant linear transformation of the output vector. The optimal feedback matrix, f*, was selected to optimize the expected performance, given the covariance of the initial state. It is first shown that the expected performance criterion can be expressed as the ratio of two multinomials in the element of f. This expression provides the basis for a feasible method of determining f* in the case of single-input single-output systems. A number of iterative algorithms are then proposed for the calculation of f* for multiple input-output systems. For two of these, monotone convergence is proved, but they involve the solution of nonlinear matrix equations at each iteration. Another is proposed involving the solution of Lyapunov equations at each iteration, and the gradual increase of the magnitude of a penalty function. Experience with this algorithm will be needed to determine whether or not it does, indeed, possess desirable convergence properties, and whether it can be used to determine the globally optimal f*.

Platzman, L. K.↗
NASA NTRS1977

Identification of state variables and dynamic inflow from rotor model dynamic tests

The paper describes methods for extracting unknown state variables and parameters from dynamic rotor model tests given transient cyclic pitch stirring inputs, blade root flap-bending measurements, and the form of the dynamic rotor equations, including a rotor dynamic inflow description, when none of the physical parameters are known. A simplified version of the maximum likelihood method seems best suited for this purpose. The measurement equation error covariance matrix is assumed constant during each iteration, but updated for the subsequent iteration. A detailed analysis of the suitability of the derived techniques for studying various rotor dynamic inflow effects is provided.

Banerjee, D.↗
NASA NTRS1977

A Navier-Stokes solution of the three-dimensional viscous compressible flow in a centrifugal compressor impeller

A two-dimensional time-dependent computer code was utilized to calculate the three-dimensional steady flow within the impeller blading. The numerical method is an explicit time marching scheme in two spatial dimensions. Initially, an inviscid solution is generated on the hub blade-to-blade surface by the method of Katsanis and McNally (1973). Starting with the known inviscid solution, the viscous effects are calculated through iteration. The approach makes it possible to take into account principal impeller fluid-mechanical effects. It is pointed out that the second iterate provides a complete solution to the three-dimensional, compressible, Navier-Stokes equations for flow in a centrifugal impeller. The problems investigated are related to the study of a radial impeller and a backswept impeller.

Harp, J. L., Jr.↗
NASA NTRS1979

Application of the AMI C sub l sub max prediction method to a number of airfoils

A method for calculating the flow about airfoils up to and beyond the stall is described. It is an iterative procedure between potential flow and boundary layer solutions. The separated region is modeled in the potential flow analysis using free vortex sheets which require an inner iteration to establish their shapes. The free vortex sheet length is an important parameter in the potential flow calculation. Results so far indicate a possible correlation between wake length and airfoil thickness/chord ratio. Calculated and experimental results are compared for a series of airfoils.

Dvorak, F. A.↗
NASA NTRS1983

Multi-grid calculation of three-dimensional transonic potential flows

A multi-grid algorithm has been developed to speed the iterative convergence of calculations for the transonic potential flow past swept wings and wing-fuselage combinations. The method is based upon a fully-conservative, finite-volume approximation to the steady potential equation which is second-order accurate everywhere in the flow field except near shock waves. The multi-grid scheme is incorporated within the framework of an alternating successive-line-overrelaxation (SLOR) solver of the difference equations. Computed results confirm the second-order accuracy of the scheme, and demonstrate the effectiveness of the multi-grid procedure.

Caughey, D. A.↗
NASA NTRS1988

Spectral element multigrid. Part 2: Theoretical justification

A multigrid algorithm is analyzed which is used for solving iteratively the algebraic system resulting from tha approximation of a second order problem by spectral or spectral element methods. The analysis, performed here in the one dimensional case, justifies the good smoothing properties of the Jacobi preconditioner that was presented in Part 1 of this paper.

Maday, Yvon↗
NASA NTRS1988

Spectral element miltigrid. II - Theoretical justification

A multigrid algorithm is analyzed which is used for solving iteratively the algebraic system resulting from the approximation of a second order problem by spectral or spectral element methods. The analysis, performed here in the one-dimensional case, justifies the good smoothing properties of the Jacobi preconditioner that was presented in Part 1 of this paper.

Maday, Yvon↗
NASA NTRS1990

Predicting Unsteady Aeroelastic Behavior

New method for predicting subsonic flutter, static deflections, and aeroelastic divergence developed. Unsteady aerodynamic loads determined by unsteady-vortex-lattice method. Accounts for aspect ratio and angle of attack. Equations for motion of wing and flow field solved iteratively and simultaneously. Used to predict transient responses to initial disturbances, and to predict steady-state static and oscillatory responses. Potential application for research in such unsteady structural/flow interactions as those in windmills, turbines, and compressors.

Strganac, Thomas W.↗
NASA NTRS1987

Program VSAERO theory document: A computer program for calculating nonlinear aerodynamic characteristics of arbitrary configurations

The VSAERO low order panel method formulation is described for the calculation of subsonic aerodynamic characteristics of general configurations. The method is based on piecewise constant doublet and source singularities. Two forms of the internal Dirichlet boundary condition are discussed and the source distribution is determined by the external Neumann boundary condition. A number of basic test cases are examined. Calculations are compared with higher order solutions for a number of cases. It is demonstrated that for comparable density of control points where the boundary conditions are satisfied, the low order method gives comparable accuracy to the higher order solutions. It is also shown that problems associated with some earlier low order panel methods, e.g., leakage in internal flows and junctions and also poor trailing edge solutions, do not appear for the present method. Further, the application of the Kutta conditions is extremely simple; no extra equation or trailing edge velocity point is required. The method has very low computing costs and this has made it practical for application to nonlinear problems requiring iterative solutions for wake shape and surface boundary layer effects.

Maskew, Brian↗
NASA NTRS1993

Fluctuations at the blue edge of saturated wind lines in IUE spectra of O-type stars

We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-Lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles and then derived from formal solution integration using this source function. The more appropriate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10 percent or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.

Owocki, Stanley P.↗