Search NASA⌕ Search

SEARCH · Search NASA

Results for “Penalty 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Preliminary design of composite wing-box structures for global damage tolerance

A procedure is presented that incorporates the influence of potential global damage conditions into the design process for minimum-mass wing-box structures. The procedure is based on mathematical-programming optimization techniques. Material-strength, minimum-gage, and panel-buckling constraints are introduced by penalty functions, and Newton's method with approximate second derivatives of the penalty terms is used as the search algorithm to obtain minimum-mass designs. A potential global damage condition is represented by a structural model with the damaged components removed. Example minimum-mass designs are obtained that simultaneously satisfy the constraints of the damaged and undamaged configurations of both graphite-epoxy and aluminum wing-box structural models. These examples are designed with and without the influence of potential damage conditions, and results indicate that for equal mass cases the residual strength of a damaged structure is higher when the influence of potential damage is properly included in the design from the outset. Results of these examples also identify the minimum structural mass increase required to increase residual strength levels.

Starnes, J. H., Jr.↗
NASA NTRS1982

Longitudinal aerodynamic characteristics of an elliptical body with a horizontal tail at Mach numbers from 2.3 to 4.63

Longitudinal aerodynamic characteristics of a configuration consisting of an elliptical body with an in plane horizontal tail were investigated. The tests were conducted at Mach numbers of 2.3, 2.96, 4.0, and 4.63. In some cases, the configuration with negative tail deflections yielded higher values of maximum lift drag ratio than did the configuration with an undeflected tail. This was due to body upwash acting on the tail and producing an additional lift increment with essentially no drag penalty. Linear theory methods used to estimate some of the longitudinal aerodynamic characteristics of the model yielded results which compared well with experimental data for all Mach numbers in this investigation and for both small angles of attack and larger angles of attack where nonlinear (vortex) flow phenomena were present.

Shrout, B. L.↗
NASA NTRS1984

Manual of phosphoric acid fuel cell power plant optimization model and computer program

An optimized cost and performance model for a phosphoric acid fuel cell power plant system was derived and developed into a modular FORTRAN computer code. Cost, energy, mass, and electrochemical analyses were combined to develop a mathematical model for optimizing the steam to methane ratio in the reformer, hydrogen utilization in the PAFC plates per stack. The nonlinear programming code, COMPUTE, was used to solve this model, in which the method of mixed penalty function combined with Hooke and Jeeves pattern search was chosen to evaluate this specific optimization problem.

Lu, C. Y.↗
NASA NTRS1978

Preliminary design of composite wings for buckling, strength and displacement constraints

An unstiffened panel buckling constraint for balanced, symmetric laminated composites is included on the global design level in a mathematical programming structural optimization procedure for designing wing structures. Constraints are introduced by penalty functions, and Newton's method based on approximate second derivatives of the penalty terms is used as the search algorithm to obtain minimum-mass designs. Constraint approximations used during the optimization process contribute to the computational efficiency of the procedure. A criterion is developed that identifies the appropriate conservative form of the constraint approximations that are used with the optimization procedure. Minimum-mass design results are obtained for a multispar high-aspect-ratio wing subjected to material strength, minimum-gage, displacement, panel buckling and twist constraints. The material systems considered for the examples are all graphite-epoxy, graphite-epoxy with boron-epoxy spar caps, and all aluminum. The composite material designs are shown to have an advantage over the aluminum designs since they can often satisfy additional constraints with only small mass increases.

Starnes, J. H., Jr.↗
NASA NTRS1977

Optimization of flexible wing structures subject to strength and induced drag constraints

An optimization procedure for designing wing structures subject to stress, strain, and drag constraints is presented. The optimization method utilizes an extended penalty function formulation for converting the constrained problem into a series of unconstrained ones. Newton's method is used to solve the unconstrained problems. An iterative analysis procedure is used to obtain the displacements of the wing structure including the effects of load redistribution due to the flexibility of the structure. The induced drag is calculated from the lift distribution. Approximate expressions for the constraints used during major portions of the optimization process enhance the efficiency of the procedure. A typical fighter wing is used to demonstrate the procedure. Aluminum and composite material designs are obtained. The tradeoff between weight savings and drag reduction is investigated.

Haftka, R. T.↗
NASA NTRS1984

Orbital transfer vehicle propulsion issues

The development of a reusable and space-based orbital transfer vehicle (OTV) necessitates an integral approach toward structural and propulsion subsystems design. A single engine installation necessitates moving the engine further aft and/or relocation of the engine gimbal point to accommodate vehicle control requirements. Penalties associated with gimbal point relocation without increasing stage length or modifying typical advanced engine concepts, as well as a method for minimizing such penalties, are presented for a single engine toroidal tank OTV configuration. Alternative integrated vehicle structure/engine concepts are also presented for multi-engine configurations. Features of these potential concepts are given which indicate the need for substantial additional study of feedline gimbal alternatives before firmly establishing advanced engine design. The issue of vehicle/engine integration is addressed in three areas: interfaces (physical and functional), installation requirements, and reliability apportionment (i.e., number of engines required to assure mission completion).

Bergeron, R. P.↗
NASA NTRS1997

Nonlinear Rescaling and Proximal-Like Methods in Convex Optimization

The nonlinear rescaling principle (NRP) consists of transforming the objective function and/or the constraints of a given constrained optimization problem into another problem which is equivalent to the original one in the sense that their optimal set of solutions coincides. A nonlinear transformation parameterized by a positive scalar parameter and based on a smooth scaling function is used to transform the constraints. The methods based on NRP consist of sequential unconstrained minimization of the classical Lagrangian for the equivalent problem, followed by an explicit formula updating the Lagrange multipliers. We first show that the NRP leads naturally to proximal methods with an entropy-like kernel, which is defined by the conjugate of the scaling function, and establish that the two methods are dually equivalent for convex constrained minimization problems. We then study the convergence properties of the nonlinear rescaling algorithm and the corresponding entropy-like proximal methods for convex constrained optimization problems. Special cases of the nonlinear resealing algorithm are presented. In particular a new class of exponential penalty-modified barrier functions methods is introduced.

Polyak, Roman↗
NASA NTRS1989

Convergence results for pseudospectral approximations of hyperbolic systems by a penalty type boundary treatment

A new method of imposing boundary conditions in the pseudospectral approximation of hyperbolic systems of equations is proposed. It is suggested to collocate the equations, not only at the inner grid points, but also at the boundary points and use the boundary conditions as penalty terms. In the pseudo-spectral Legrendre method with the new boundary treatment, a stability analysis for the case of a constant coefficient hyperbolic system is presented and error estimates are derived.

Funaro, Daniele↗
NASA NTRS1989

The art and science of flow control

The ability to actively or passively manipulate a flow field to effect a desired change is of immense technological importance. In this article, methods of control to achieve transition delay, separation postponement, lift enhancement, drag reduction, turbulence augmentation, or noise suppression are considered. Emphasis is placed on external boundary-layer flows although applicability of some of the methods reviewed for internal flows will be mentioned. Attempts will be made to present a unified view of the different methods of control to achieve a variety of end results. Performance penalties associated with a particular method such as cost, complexity, or trade-off will be elaborated.

Gad-El-hak, Mohamed↗
NASA NTRS1999

Low-Storage, Explicit Runge-Kutta Schemes for the Compressible Navier-Stokes Equations

The derivation of storage explicit Runge-Kutta (ERK) schemes has been performed in the context of integrating the compressible Navier-Stokes equations via direct numerical simulation. Optimization of ERK methods is done across the broad range of properties, such as stability and accuracy efficiency, linear and nonlinear stability, error control reliability, step change stability, and dissipation/dispersion accuracy, subject to varying degrees of memory economization. Following van der Houwen and Wray, 16 ERK pairs are presented using from two to five registers of memory per equation, per grid point and having accuracies from third- to fifth-order. Methods have been assessed using the differential equation testing code DETEST, and with the 1D wave equation. Two of the methods have been applied to the DNS of a compressible jet as well as methane-air and hydrogen-air flames. Derived 3(2) and 4(3) pairs are competitive with existing full-storage methods. Although a substantial efficiency penalty accompanies use of two- and three-register, fifth-order methods, the best contemporary full-storage methods can be pearl), matched while still saving two to three registers of memory.

Kennedy, Chistopher A.↗
NASA NTRS1972

Synthesis of stiffened conical shells.

The development of a method to effect the automated minimum weight design of ring and stringer stiffened shells is presented. Membrane theory is used for the shell prebuckling analysis. The buckling analysis is based upon an arbitrary shell of revolution computer program. The structural analysis includes both buckling and yielding modes of failure. The synthesis involves the coupling of an exterior penalty function with a method for the unconstrained minimization of a function comprised of a sum of squares. Results of the application of the method to the design of the Viking Aeroshell cone are presented. The least weight Viking Aeroshell appears to be an all magnesium shell with ring stiffeners of hollow circular cross section. Because the method incorporates a general shell of revolution buckling analysis, it can be readily modified and applied to the design of any axisymmetrically loaded uniformly stiffened shell of revolution for which a membrane prebuckling solution exists.

Thornton, W. A.↗
NASA NTRS1986

Impact of mismodeled idle engine performance on calculation and tracking of optimal 4-D descent trajectories

Advanced flight management systems are being developed which are capable of calculating optimal 3-D and 4-D flight trajectories for arbitrary fuel and time costs. These systems require mathematical models of airplane performance in order to compute the optimal profiles. Mismodeled idle engine characteristics can result in descent trajectories requiring excessive throttle and/or speedbrake activity in order to achieve the desired end conditions. This paper evaluates the cost and fuel penalties, trajectory variations, and flight control requirements associated with typical idle engine modeling errors for a twin-jet transport airplane. Variations in idle power setting, thrust, fuel flow, and surge bleed operation were evaluated for a cruise/descent flight segment. The results of this analysis provide insight into the penalties associated with uncertainties in idle engine performance and suggest methods of modeling which minimize these penalties.

Williams, D. H.↗
NASA NTRS1977

Comparison of some optimal control methods for the design of turbine blades

This paper attempts a comparative study of some numerical methods for the optimal control design of turbine blades whose vibration characteristics are approximated by Timoshenko beam idealizations with shear and incorporating simple boundary conditions. The blade was synthesized using the following methods: (1) conjugate gradient minimization of the system Hamiltonian in function space incorporating penalty function transformations, (2) projection operator methods in a function space which includes the frequencies of vibration and the control function, (3) epsilon-technique penalty function transformation resulting in a highly nonlinear programming problem, (4) finite difference discretization of the state equations again resulting in a nonlinear program, (5) second variation methods with complex state differential equations to include damping effects resulting in systems of inhomogeneous matrix Riccatti equations some of which are stiff, (6) quasi-linear methods based on iterative linearization of the state and adjoint equation. The paper includes a discussion of some substantial computational difficulties encountered in the implementation of these techniques together with a resume of work presently in progress using a differential dynamic programming approach.

Desilva, B. M. E.↗
NASA NTRS1990

A balloon borne 19-GHz radiometer

The design and performance of a 19-GHz one-horn maser radiometer is discussed. The antenna is a lensed, conical horn with 3-deg angular resolution. The radiometer is stabilized by chopping against an internal cold load. A novel chopping scheme allows continual calibration of maser gain fluatuations with only a modest penalty in noise. The method of calibrating this instrument is described. The instrument has been used to map the sky with an equivalent blackbody temperature sensitivity of 1 mK per 3 x 3 deg resolution element.

Boughn, S. P.↗
NASA NTRS1991

Airfoil optimization with efficient gradient calculations

The viscous airfoil design analysis code XFOIL was extended to allow optimization using conformal mapping coefficients as design variables. The optimization technique used was the Steepest Descent method applied to a Penalty Function. The gradients of the aerodynamic variables with respect to the design variables were cheaply calculated as by-products of XFOIL's integral boundary layer Newton solver. The speed of the optimization process further increased by updating the Newton system boundary layer variables after each optimization step using the available gradient information. Two examples are presented.

Sorensen, Thomas↗
NASA NTRS2004

Penalty-Based Interface Technology for Prediction of Delamination Growth in Laminated Structures

An effective interface element technology has been developed for connecting and simulating crack growth between independently modeled finite element subdomains (e.g., composite plies). This method has been developed using penalty constraints and allows coupling of finite element models whose nodes do not necessarily coincide along their common interface. Additionally, the present formulation leads to a computational approach that is very efficient and completely compatible with existing commercial software. The present interface element has been implemented in the commercial finite element code ABAQUS as a User Element Subroutine (UEL), making it easy to test the approach for a wide range of problems. The interface element technology has been formulated to simulate delamination growth in composite laminates. Thanks to its special features, the interface element approach makes it possible to release portions of the interface surface whose length is smaller than that of the finite elements. In addition, the penalty parameter can vary within the interface element, allowing the damage model to be applied to a desired fraction of the interface between the two meshes. Results for double cantilever beam DCB, end-loaded split (ELS) and fixed-ratio mixed mode (FRMM) specimens are presented. These results are compared to measured data to assess the ability of the present damage model to simulate crack growth.

Averill, Ronald C.↗
NASA NTRS1976

Evaluation of a cost-effective loads approach

A shock spectra/impedance method for loads prediction is used to estimate member loads for the Viking Orbiter, a 7800-lb interplanetary spacecraft that has been designed using transient loads analysis techniques. The transient loads analysis approach leads to a lightweight structure but requires complex and costly analyses. To reduce complexity and cost a shock spectra/impedance method is currently being used to design the Mariner Jupiter Saturn spacecraft. This method has the advantage of using low-cost in-house loads analysis techniques and typically results in more conservative structural loads. The method is evaluated by comparing the increase in Viking member loads to the loads obtained by the transient loads analysis approach. An estimate of the weight penalty incurred by using this method is presented. The paper also compares the calculated flight loads from the transient loads analyses and the shock spectra/impedance method to measured flight data.

Garba, J. A.↗
NASA NTRS1976

Evaluation of a cost-effective loads approach

A shock spectra/impedance method for loads predictions is used to estimate member loads for the Viking Orbiter, a 7800-lb interplanetary spacecraft that has been designed using transient loads analysis techniques. The transient loads analysis approach leads to a lightweight structure but requires complex and costly analyses. To reduce complexity and cost, a shock spectra/impedance method is currently being used to design the Mariner Jupiter Saturn spacecraft. This method has the advantage of using low-cost in-house loads analysis techniques and typically results in more conservative structural loads. The method is evaluated by comparing the increase in Viking member loads to the loads obtained by the transient loads analysis approach. An estimate of the weight penalty incurred by using this method is presented. The paper also compares the calculated flight loads from the transient loads analyses and the shock spectra/impedance method to measured flight data.

Garba, J. A.↗