SEARCH · Engineering Papers

Results for “successive linear approximation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2006

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗

NASA NTRS1982

Analysis and design of steady transonic flow over airfoils by the method of parametric differentiation

The results reported here are based on applying the method of parametric differentiation (MPD) to transform the nonlinear differential equation governing small-disturbance transonic flow to a linear equation. Implicit approximate factorization and monotone methods were used to accelerate convergence of the linear problem by an order of magnitude over successive line over-relaxation. The relative merits of using MPD are discussed in comparison to conventional small-disturbance applications. Several MPD analyses are performed on an array of airfoils. A design procedure utilizing MPD is discussed and demonstrated for two nonlifting cases.

Halt, D. W.↗

NASA NTRS1989

Overview of Krylov subspace methods with applications to control problems

An overview of projection methods based on Krylov subspaces are given with emphasis on their application to solving matrix equations that arise in control problems. The main idea of Krylov subspace methods is to generate a basis of the Krylov subspace Span and seek an approximate solution the the original problem from this subspace. Thus, the original matrix problem of size N is approximated by one of dimension m typically much smaller than N. Krylov subspace methods have been very successful in solving linear systems and eigenvalue problems and are now just becoming popular for solving nonlinear equations. It is shown how they can be used to solve partial pole placement problems, Sylvester's equation, and Lyapunov's equation.

Saad, Youcef↗

NASA NTRS1985

Sensor/Actuator Selection for the Constrained Variance Control Problem

The problem of designing a linear controller for systems subject to inequality variance constraints is considered. A quadratic penalty function approach is used to yield a linear controller. Both the weights in the quadratic penalty function and the locations of sensors and actuators are selected by successive approximations to obtain an optimal design which satisfies the input/output variance constraints. The method is applied to NASA's 64 meter Hoop-Column Space Antenna for satellite communications. In addition the solution for the control law, the main feature of these results is the systematic determination of actuator design requirements which allow the given input/output performance constraints to be satisfied.

Delorenzo, M. L.↗

NASA NTRS1985

Nonlinear acoustic wave propagation in atmosphere

A model problem that simulates an atmospheric acoustic wave propagation situation that is nonlinear is considered. The model is derived from the basic Euler equations for the atmospheric flow and from the regular perturbations for the acoustic part. The nonlinear effects are studied by obtaining two successive linear problems in which the second one involves the solution of the first problem. Well posedness of these problems is discussed and approximations of the radiation boundary conditions that can be used in numerical simulations are presented.

Hariharan, S. I.↗

NASA NTRS1984

A rotor optimization using regression analysis

The design and development of helicopter rotors is subject to the many design variables and their interactions that effect rotor operation. Until recently, selection of rotor design variables to achieve specified rotor operational qualities has been a costly, time consuming, repetitive task. For the past several years, Kaman Aerospace Corporation has successfully applied multiple linear regression analysis, coupled with optimization and sensitivity procedures, in the analytical design of rotor systems. It is concluded that approximating equations can be developed rapidly for a multiplicity of objective and constraint functions and optimizations can be performed in a rapid and cost effective manner; the number and/or range of design variables can be increased by expanding the data base and developing approximating functions to reflect the expanded design space; the order of the approximating equations can be expanded easily to improve correlation between analyzer results and the approximating equations; gradients of the approximating equations can be calculated easily and these gradients are smooth functions reducing the risk of numerical problems in the optimization; the use of approximating functions allows the problem to be started easily and rapidly from various initial designs to enhance the probability of finding a global optimum; and the approximating equations are independent of the analysis or optimization codes used.

Giansante, N.↗

NASA NTRS2017

Understanding the Sun-Earth Libration Point Orbit Formation Flying Challenges For WFIRST and Starshade

In order to fly an occulter in formation with a telescope at the Sun-Earth L2 (SEL2) Libration Point, one must have a detailed understanding of the dy-namics that govern the restricted three body system. For initial purposes, a linear approximation is satisfactory, but operations will require a high-fidelity modeling tool along with strategic targeting methods in order to be successful. This paper focuses on the challenging dynamics of the transfer trajectories to achieve the relative positioning of two spacecraft to fly in formation at SEL2, in our case, the Wide-Field Infrared Survey Telescope (WFIRST) and a proposed Starshade. By modeling the formation transfers using a high fidelity tool, an accurate V approximation can be made to as-sist with the development of the subsystem design required for a WFIRST and Starshade formation flight mission.

Sun-Earth-Moon↗

NASA NTRS1986

Nonlinear acoustic wave propagation in atmosphere

In this paper a model problem is considered that simulates an atmospheric acoustic wave propagation situation that is nonlinear. The model is derived from the basic Euler equations for the atmospheric flow and from the regular perturbations for the acoustic part. The nonlinear effects are studied by obtaining two successive linear problems in which the second one involves the solution of the first problem. Well-posedness of these problems is discussed and approximations of the radiation boundary conditions that can be used in numerical simulations are presented.

Hariharan, S. I.↗

NASA NTRS1987

Nonlinear acoustic wave propagation in atmosphere TLSP - Final report

Hariharan, S. I.↗

NASA NTRS1987

Model-size reduction for the analysis of symmetric structures with asymmetric boundary conditions

A simple computational procedure is presented for reducing the size of the analysis model for a symmetric structure with asymmetric boundary conditions to that of the corresponding structure with symmetric boundary conditions. The procedure is based on approximating the asymmetric response of the structure by a linear combination of symmetric and antisymmetric global approximation vectors (or modes). The key elements of the procedure are (1) restructuring the governing finite-element equations to delineate the contributions to the symmetric and antisymmetric components of the asymmetric response, (2) successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite-element method is first used to generate a few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The effectiveness of the computational procedure is demonstrated by means of numerical examples of linear static problems of shells, and its potential for solving nonlinear problems is discussed.

Noor, Ahmed K.↗

NASA NTRS1977

Method and apparatus for measuring web material wound on a reel

The method and apparatus for measuring the number of layers of a web material of known thickness wound on a storage or take-up reel is presented. The method and apparatus are based on the principle that, at a relatively large radius, the loci of layers of a thin web wound on the reel approximate a family of concentric circles having radii respectively successively increasing by a length equal to the web thickness. Tachometer pulses are generated in response to linear movement of the web and reset pulses are generated in response to rotation of the reel. A digital circuit, responsive to the tachometer and reset pulses, generates data indicative of the layer number of any layer of the web and of position of the web within the layer without requiring numerical interpolation.

Muller, R. M.↗

NASA NTRS1946

Numerical procedures for the calculation of the stresses in monocoques III : calculation of the bending moments in fuselage frames

This report deals with the calculation of the bending moments in and the distortions of fuselage rings upon which known concentrated and distributed loads are acting. In the procedure suggested, the ring is divided into a number of beams each having a constant radius of curvature. The forces and moments caused in the end sections of the beams by individual unit displacements of the end sections are listed in a table designated as the operations table in conformity with Southwell's nomenclature. The operations table and the external loads are equivalent to a set of linear equations. For their solution the following three procedures are presented: 1) Southwell's method of systematic relaxations. This is a step-by-step approximation procedure guided by the physical interpretation of the changes in the values of the unknown. 2) The growing unit procedure in which the individual beams are combined successively into beams of increasing length until finally the entire ring becomes a single beam. In each step of the procedure a set of not more than three simultaneous linear equations is solved. 3) Solution of the entire set of simultaneous equations by the methods of the matrix calculus. In order to demonstrate the manner in which the calculations may be carried out, the following numerical examples are worked out: 1) Curved beam with both its end sections rigidly fixed. The load is a concentrated force. 2) Egg-shape ring with symmetric concentrated loads. 3) Circular ring with antisymmetric concentrated loads and shear flow (torsion of the fuselage). 4) Same with V-braces incorporated in the ring. 5) Egg-shape ring with antisymmetric concentrated loads and shear flow (torsion of the fuselage). 6) Same with V-braces incorporated in the ring. The results of these calculations are checked, whenever possible, by calculations carried out according to known methods of analysis. The agreement is found to be good. The amount of work necessary for the solution of ring problems by the methods described in the present report is practically independent of the degree of redundancy of the structure. For this reason the methods are recommended for use particularly in problems of rings having one or more internal bracing elements.

LOADS - FUSELAGES↗

NASA NTRS2003

Genetic Algorithm Tuned Fuzzy Logic for Gliding Return Trajectories

The problem of designing and flying a trajectory for successful recovery of a reusable launch vehicle is tackled using fuzzy logic control with genetic algorithm optimization. The plant is approximated by a simplified three degree of freedom non-linear model. A baseline trajectory design and guidance algorithm consisting of several Mamdani type fuzzy controllers is tuned using a simple genetic algorithm. Preliminary results show that the performance of the overall system is shown to improve with genetic algorithm tuning.

Burchett, Bradley T.↗

NASA NTRS2003

Capacitive Sensors for Measuring Masses of Cryogenic Fluids

An effort is under way to develop capacitive sensors for measuring the masses of cryogenic fluids in tanks. These sensors are intended to function in both microgravitational and normal gravitational settings, and should not be confused with level sensors, including capacitive ones. A sensor of this type is conceptually simple in the sense that (1) it includes only one capacitor and (2) if properly designed, its single capacitance reading should be readily convertible to a close approximation of the mass of the cryogenic fluid in the tank. Consider a pair of electrically insulated electrodes used as a simple capacitive sensor. In general, the capacitance is proportional to the permittivity of the dielectric medium (in this case, a cryogenic fluid) between the electrodes. The success of design and operation of a sensor of the present type depends on the accuracy of the assumption that to a close approximation, the permittivity of the cryogenic fluid varies linearly with the density of the fluid. Data on liquid nitrogen, liquid oxygen, and liquid hydrogen, reported by the National Institute of Standards and Technology, indicate that the permittivities and densities of these fluids are, indeed, linearly related to within a few tenths of a percent over the pressure and temperature regions of interest. Hence, ignoring geometric effects for the moment, the capacitance between two electrodes immersed in the fluid should vary linearly with the density, and, hence, with the mass of the fluid. Of course, it is necessary to take account of the tank geometry. Because most cryogenic tanks do not have uniform cross sections, the readings of level sensors, including capacitive ones, are not linearly correlated with the masses of fluids in the tanks. In a sensor of the present type, the capacitor electrodes are shaped so that at a given height, the capacitance per unit height is approximately proportional to the cross-sectional area of the tank in the horizontal plane at that height (see figure).

Nurge, Mark↗

NASA NTRS1993

Neural network-based control using Lyapunov functions

We have successfully demonstrated how the problem of stabilization of plants can be reduced to a problem of approximation of functions. Neural networks have been shown to have approximating and interpolating properties. This approach is good for linear and nonlinear plants. Software has been generated to demonstrate this approach.

Luxemburg, Leon A.↗

NASA NTRS1994

Nonperturbative methods in HZE propagation

An analytical solution to the perturbative multiple collision series of a fragmenting HZE ion beam has limited usefulness since the first collision term has several hundred contributions, the second collision term has tens of thousands of contributions, and each successive collision term progresses to unwieldy computational proportions. Our previous work has revealed the multiple collision terms in the straight-ahead approximation to be simple products of a spatially dependent factor times a linear energy-dependent factor of limited domain and unit normalization. The properties of these forms allow the development of the nonperturbative summation of the series to all orders assuming energy-independent nuclear cross sections as matrix products of a scaled Green's function described herein. This nonperturbative Green's function with multiple scattering correction factors compares well with experiments using 670 MeV/u neon-20 ion beams in thick water targets.

Non-NASA Center↗

NASA NTRS1995

Modeling of collision induced absorption spectra of CO2-CO2 pairs for planetary atmosphere of Venus

The objective of the proposal was to model the rototranslational and rotovibrational collision induced absorption spectral bands of importance for the radiative transfer analysis of the atmosphere of Venus. Our main task has involved CO2 pairs. The approach is not straightforward: whereas computational techniques to compute CIA spectra of small linear molecules exist, and were successfully applied to molecules like H2 or N2, they fail when applied to large molecules like CO2. For small molecules one can safely assume that the interaction potential is isotropic. The same approximation does not work for CO2, and when employed, it gives an incorrect band shape and only 50 percent of the CIA intensity.

Borysow, Aleksandra↗

NASA NTRS1991

Implicit solvers for unstructured meshes

Implicit methods for unstructured mesh computations are developed and tested. The approximate system which arises from the Newton-linearization of the nonlinear evolution operator is solved by using the preconditioned generalized minimum residual technique. These different preconditioners are investigated: the incomplete LU factorization (ILU), block diagonal factorization, and the symmetric successive over-relaxation (SSOR). The preconditioners have been optimized to have good vectorization properties. The various methods are compared over a wide range of problems. Ordering of the unknowns, which affects the convergence of these sparse matrix iterative methods, is also investigated. Results are presented for inviscid and turbulent viscous calculations on single and multielement airfoil configurations using globally and adaptively generated meshes.

Venkatakrishnan, V.↗