SEARCH · Search NASA

Results for “penalty methods”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

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↗

DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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.↗

NASA NTRS1986

New methodology for shaft design based on life expectancy

The design of power transmission shafting for reliability has not historically received a great deal of attention. However, weight sensitive aerospace and vehicle applications and those where the penalties of shaft failure are great, require greater confidence in shaft design than earlier methods provided. This report summarizes a fatigue strength-based, design method for sizing shafts under variable amplitude loading histories for limited or nonlimited service life. Moreover, applications factors such as press-fitted collars, shaft size, residual stresses from shot peening or plating, corrosive environments can be readily accommodated into the framework of the analysis. Examples are given which illustrate the use of the method, pointing out the large life penalties due to occasional cyclic overloads.

Loewenthal, S. H.↗

NASA NTRS1979

A cubic extended interior penalty function for structural optimization

This paper describes an optimization procedure for the minimum weight design of complex structures. The procedure is based on a new cubic extended interior penalty function (CEIPF) used with the sequence of unconstrained minimization technique (SUMT) and Newton's method. The Hessian matrix of the penalty function is approximated using only constraints and their derivatives. The CEIPF is designed to minimize the error in the approximation of the Hessian matrix, and as a result the number of structural analyses required is small and independent of the number of design variables. Three example problems are reported. The number of structural analyses is reduced by as much as 50 per cent below previously reported results.

Prasad, B.↗

NASA NTRS1988

Stabilization of computational procedures for constrained dynamical systems

A new stabilization method of treating constraints in multibody dynamical systems is presented. By tailoring a penalty form of the constraint equations, the method achieves stabilization without artificial damping and yields a companion matrix differential equation for the constraint forces; hence, the constraint forces are obtained by integrating the companion differential equation for the constraint forces in time. A principal feature of the method is that the errors committed in each constraint condition decay with its corresponding characteristic time scale associated with its constraint force. Numerical experiments indicate that the method yields a marked improvement over existing techniques.

Park, K. C.↗

NASA NTRS1978

A new formulation for the epsilon method applied to the minimum-time-to-climb problem

Balakrishnan's epsilon technique is used to compute minimum-time profiles for the F-104 airplane. This technique differs from the classical gradient method in that a quadratic penalty on the error in satisfying the equation of motion is included in the cost function to be minimized as a means of eliminating the requirement of satisfying the equations of motion. Although the number of unknown independent functions is increased to include the state variables, the evaluation of the gradient of the cost function is simplified, resulting in considerable computational savings, thereby making it appear feasible to use the epsilon method for real-time application.

Taylor, L. W., Jr.↗

NASA NTRS1979

A multi-disciplinary approach to structural design for stochastic loads

The development of aircraft gust design techniques is reviewed with emphasis on those techniques applicable to large scale multidisciplinary closed-loop optimization methods. An overall view is presented of the optimization methodology utilized. An analysis is presented wherein the bivariate normal distribution function for a particular loading condition is considered. By use of maximization techniques, functions are derived for stress components along equal probability contours which allow for directionality and auto-correlation effects to be considered. These functions are then applied to the problem being considered and the effects on the structural sizing and weight penalties evaluated against more conventional methods.

Gross, D. W.↗

NASA NTRS1984

Penalty function finite element analysis of steady viscous incompressible flow in rotating coordinates

Finite element methods for incompressible viscous flow in turbomachines have not been presented in the literature previously. This paper develops a penalty function primitive variable method including Coriolis and centrifugal force terms for steady flow in a rotating coordinate system. Simplex elements are used with the result of solution times comparable to equivalent finite different solutions. Example cases considered are Couette flow, Poiseuile flow, flow over a step and flow in a rotating channel. Both laminar and turbulent flows are discussed. The accuracy of computed solutions compares well with theoretical solutions and experimenal measurements.

Rosen, M. C.↗

DOE OSTI2024

Dynamic flux surrogate-based partitioned methods for interface problems

Loosely coupled partitioned methods for multiphysics problems treat each subproblem as a separate entity and advance them independently in time. In so doing these methods enable code reuse, increase concurrency and provide a convenient framework for plug-and-play multiphysics simulations. However, mathematically loosely coupled schemes are equivalent to a single step of an iterative solution method, which can compromise their accuracy and stability. We present a new data-driven partitioned method for coupled parametric PDEs that can improve upon the accuracy of traditional loosely coupled methods without incurring a performance penalty. To that end, we replace conventional field transfers across the interface by a surrogate for the dynamics of the interface flux exchanged between the subdomains. To develop this surrogate we apply dynamic mode decomposition to a non-standard staggered-in-time state, comprising the interface flux and small solution patches near the interface. The new approach shifts the main computational burden to an offline training phase, whereas application of the surrogate in the online phase amounts to a single matrix–vector multiplication. In conclusion, we provide stability analysis of the surrogate-based partitioned scheme and include numerical results that demonstrate its potential.

Dynamic mode decomposition (DMD)↗

NASA NTRS1972

A discrete variable approximation to minimum weight panel designs subject to a supersonic flutter-speed constraint.

A numerical technique is presented that determines the minimum-weight thickness distribution for a one-dimensional, simply supported panel subjected to a parallel supersonic flow on one side. An aeroelastic eigenvalue, which characterizes flutter speed, is held constant. The governing differential equations are approximated by sets of difference equations adjoined to the weight function via a penalty function. A conjugate gradient method then solves the resulting sequence of unconstrained minimization problems. Numerical results are obtained without imposing a minimum thickness constraint; the effect of constant inplane tensile and compressive stresses is also investigated.

Pierson, B. L.↗

NASA NTRS1978

Programs for analysis and resizing of complex structures

The paper describes the PARS (Programs for Analysis and Resizing of Structures) system. PARS is a user oriented system of programs for the minimum weight design of structures modeled by finite elements and subject to stress, displacement, flutter and thermal constraints. The system is built around SPAR - an efficient and modular general purpose finite element program, and consists of a series of processors that communicate through the use of a data base. An efficient optimizer based on the Sequence of Unconstrained Minimization Technique (SUMT) with an extended interior penalty function and Newton's method is used. Several problems are presented for demonstration of the system capabilities.

Haftka, R. T.↗

NASA NTRS1983

Thrust performance of a variable-geometry, divergent exhaust nozzle on a turbojet engine at altitude

A variable geometry, low aspect ratio, nonaxisymmetric, two dimensional, convergent-divergent exhaust nozzle was tested at simulated altitude on a turbojet engine to obtain baseline axial, dry thrust performance over wide ranges of operating nozzle pressure ratios, throat areas, and internal expansion area ratios. The thrust data showed good agreement with theory and scale model test results after the data were corrected for seal leakage and coolant losses. Wall static pressure profile data were also obtained and compared with one dimensional theory and scale model data. The pressure data indicate greater three dimensional flow effects in the full scale tests than with models. The leakage and coolant penalties were substantial, and the method to determine them is included.

Straight, D. M.↗