Search NASA⌕ Search

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.

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Ephemeral penalty functions for contact-impact dynamics

The use of penalty functions to treat a class of structural contact-impact problems is investigated, with emphasis on ones in which the impact phenomena are primarily nondestructive in nature and in which only the gross characterization of the response is required. The dynamic equations of motion are integrated by the difference method. The penalty is represented as an ephemeral fictitious nonlinear spring that is inserted on anticipation of contact. The magnitude and variation of the penalty force is determined through energy balancing considerations. The 'bell shape' of the penalty force function for positive gap was found to be satisfactory, as it depends on only two parameters that can be directly assigned the physical meaning of force and distance. The determination of force law parameters by energy balance worked well. The incorporation of restitution coefficients by the area balancing method yielded excellent results, and no substantial modifications are anticipated. Extensional penalty springs are obviously sufficient for the simple examples treated.

De La Fuente, Horacio M.↗
NASA NTRS2002

Very Large Scale Optimization

The purpose of this research under the NASA Small Business Innovative Research program was to develop algorithms and associated software to solve very large nonlinear, constrained optimization tasks. Key issues included efficiency, reliability, memory, and gradient calculation requirements. This report describes the general optimization problem, ten candidate methods, and detailed evaluations of four candidates. The algorithm chosen for final development is a modern recreation of a 1960s external penalty function method that uses very limited computer memory and computational time. Although of lower efficiency, the new method can solve problems orders of magnitude larger than current methods. The resulting BIGDOT software has been demonstrated on problems with 50,000 variables and about 50,000 active constraints. For unconstrained optimization, it has solved a problem in excess of 135,000 variables. The method includes a technique for solving discrete variable problems that finds a "good" design, although a theoretical optimum cannot be guaranteed. It is very scalable in that the number of function and gradient evaluations does not change significantly with increased problem size. Test cases are provided to demonstrate the efficiency and reliability of the methods and software.

Vanderplaats, Garrett↗
NASA NTRS2007

Simulating Space Capsule Water Landing with Explicit Finite Element Method

A study of using an explicit nonlinear dynamic finite element code for simulating the water landing of a space capsule was performed. The finite element model contains Lagrangian shell elements for the space capsule and Eulerian solid elements for the water and air. An Arbitrary Lagrangian Eulerian (ALE) solver and a penalty coupling method were used for predicting the fluid and structure interaction forces. The space capsule was first assumed to be rigid, so the numerical results could be correlated with closed form solutions. The water and air meshes were continuously refined until the solution was converged. The converged maximum deceleration predicted is bounded by the classical von Karman and Wagner solutions and is considered to be an adequate solution. The refined water and air meshes were then used in the models for simulating the water landing of a capsule model that has a flexible bottom. For small pitch angle cases, the maximum deceleration from the flexible capsule model was found to be significantly greater than the maximum deceleration obtained from the corresponding rigid model. For large pitch angle cases, the difference between the maximum deceleration of the flexible model and that of its corresponding rigid model is smaller. Test data of Apollo space capsules with a flexible heat shield qualitatively support the findings presented in this paper.

Wang, John T.↗
DOE OSTI2025

Consistent Second Moment Methods with Scalable Linear Solvers for Radiation Transport

Second moment methods (SMMs) are developed that are consistent with the discontinuous Galerkin spatial discretization of the discrete ordinates (or S\(_N\)) transport equations. The low-order (LO) diffusion system of equations is discretized with fully consistent P\(_1\), local discontinuous Galerkin (LDG), and interior penalty (IP) methods. A discrete residual approach is used to derive SMM correction terms that make each of the LO systems consistent with the high-order discretization. We show that the consistent methods are more accurate and have better solution quality than independently discretized LO systems, that they preserve the diffusion limit, and that the LDG and IP consistent SMMs can be scalably solved in parallel on a challenging, multimaterial benchmark problem.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1980

ACCESS 3. Approximation concepts code for efficient structural synthesis: User's guide

A user's guide is presented for ACCESS-3, a research oriented program which combines dual methods and a collection of approximation concepts to achieve excellent efficiency in structural synthesis. The finite element method is used for structural analysis and dual algorithms of mathematical programming are applied in the design optimization procedure. This program retains all of the ACCESS-2 capabilities and the data preparation formats are fully compatible. Four distinct optimizer options were added: interior point penalty function method (NEWSUMT); second order primal projection method (PRIMAL2); second order Newton-type dual method (DUAL2); and first order gradient projection-type dual method (DUAL1). A pure discrete and mixed continuous-discrete design variable capability, and zero order approximation of the stress constraints are also included.

Fleury, C.↗
NASA NTRS1994

CFD assessment of the carbon monoxide and nitric oxide formation from RD-170 hot-fire testing at MSFC

Computational Fluid Dynamics (CFD) technology has been used to assess the exhaust plume pollutant environment of the RD-170 engine hot-firing on the F1 Test Stand at Marshall Space Flight Center. Researchers know that rocket engine hot-firing has the potential for forming thermal nitric oxides (NOx), as well as producing carbon monoxide (CO) when hydrocarbon fuels are used. Because of the complicated physics involved, however, little attempt has been made to predict the pollutant emissions from ground-based engine testing, except for simplified methods which can grossly underpredict and/or overpredict the pollutant formations in a test environment. The objective of this work, therefore, has been to develop a technology using CFD to describe the underlying pollutant emission physics from ground-based rocket engine testing. This resultant technology is based on a three-dimensional, viscous flow, pressure-based CFD formulation, where wet CO and thermal NOx finite-rate chemistry mechanisms are solved with a Penalty function method.

Wang, Ten-See↗
NASA NTRS1995

CFD assessment of the pollutant environment from RD-170 propulsion system testing

Computational Fluid Dynamics (CFD) technology has been used to assess the exhaust plume pollutant environment of the RD-170 engine hot-firing on the F1 Test Stand at Marshall Space Flight Center. Researchers know that rocket engine hot-firing has the potential for forming thermal nitric oxides (NO(x)), as well as producing carbon monoxide (CO) when hydrocarbon fuels are used. Because of the complicated physics involved, however, little attempt has been made to predict the pollutant emissions from ground-based engine testing, except for simplified methods which can grossly underpredict and/or overpredict the pollutant formations in a test environment. The objective of this work, therefore, has been to develop a technology using CFD to describe the underlying pollutant emission physics from ground-based rocket engine testing. This resultant technology is based on a three-dimensional (3D), viscous flow, pressure-based CFD formulation, where wet CO and thermal NO finite-rate chemistry mechanisms are solved with a Penalty Function method. A nominal hot-firing of a RD-170 engine on the F1 stand has been computed. Pertinent test stand flow physics such as the multiple-nozzle clustered engine plume interaction, air aspiration from base and aspirator, plume mixing with entrained air that resulted in contaminant dilution and afterburning, counter-afterburning due to flame bucket water-quenching, plume impingement on the flame bucket, and restricted multiple-plume expansion and turning have been captured. The predicted total emission rates compared reasonably well with those of the existing hydrocarbon engine hot-firing test data.

Wang, Ten-See↗
NASA NTRS1999

Structural Optimization for Reliability Using Nonlinear Goal Programming

This report details the development of a reliability based multi-objective design tool for solving structural optimization problems. Based on two different optimization techniques, namely sequential unconstrained minimization and nonlinear goal programming, the developed design method has the capability to take into account the effects of variability on the proposed design through a user specified reliability design criterion. In its sequential unconstrained minimization mode, the developed design tool uses a composite objective function, in conjunction with weight ordered design objectives, in order to take into account conflicting and multiple design criteria. Multiple design criteria of interest including structural weight, load induced stress and deflection, and mechanical reliability. The nonlinear goal programming mode, on the other hand, provides for a design method that eliminates the difficulty of having to define an objective function and constraints, while at the same time has the capability of handling rank ordered design objectives or goals. For simulation purposes the design of a pressure vessel cover plate was undertaken as a test bed for the newly developed design tool. The formulation of this structural optimization problem into sequential unconstrained minimization and goal programming form is presented. The resulting optimization problem was solved using: (i) the linear extended interior penalty function method algorithm; and (ii) Powell's conjugate directions method. Both single and multi-objective numerical test cases are included demonstrating the design tool's capabilities as it applies to this design problem.

El-Sayed, Mohamed E.↗
NASA NTRS1992

Global/local analysis of laminated composite plates using variable kinematic finite elements

A finite element modeling methodology is developed for the hierarchical, global/local analysis of laminated composite plates. The method incorporates a new variable kinematics, displacement-based, finite element that is developed using a multiple assumed displacement field approach. The variable kinematic elements provide a great degree of flexibility in defining the transverse (through thickness) variation of the assumed displacement field. The resulting finite element model permits different subregions of the computational domain to be described by different mathematical models. Enforcing displacement continuity along subregion boundaries requires only the specification of certain homogeneous essential boundary conditions, thus avoiding the inconvenience of multi-point constraints, penalty function methods, or special transition elements.

Robbins, D. H., Jr.↗
DOE OSTI2025

A time-parallel multiple-shooting method for large-scale quantum optimal control

Quantum optimal control plays a crucial role in quantum computing by providing the interface between compiler and hardware. Solving the optimal control problem is particularly challenging for multi-qubit gates, due to the exponential growth in computational complexity with the system's dimensionality and the deterioration of optimization convergence. To ameliorate the computational complexity of time-integration, this paper introduces a multiple-shooting approach in which the time domain is divided into multiple windows and the intermediate states at window boundaries are treated as additional optimization variables. Further, this enables parallel computation of state evolution across time-windows, significantly accelerating objective function and gradient evaluations. Since the initial state matrix in each window is only guaranteed to be unitary upon convergence of the optimization algorithm, the conventional gate trace infidelity is replaced by a generalized infidelity that is convex for non-unitary state matrices. Continuity of the state across window boundaries is enforced by equality constraints. A quadratic penalty optimization method is used to solve the constrained optimal control problem, and an efficient adjoint technique is employed to calculate the gradients in each iteration. We demonstrate the effectiveness of the proposed method through numerical experiments on quantum Fourier transform gates in systems with 2, 3, and 4 qubits, noting a speedup of 80x for evaluating the gradient in the 4-qubit case, highlighting the method's potential for optimizing control pulses in multi-qubit quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

Optimal rocket thrust profile shaping using third degree spline function interpolation

Optimal solid-rocket thrust profiles for the parallel-burn, solid-rocket-assisted space shuttle are investigated. Solid-rocket thrust profiles are simulated by using third-degree spline functions, with the values of the thrust ordinates defined as parameters. The profiles are optimized parametrically, using the Davidon-Fletcher-Powell penalty function method, by minimizing propellant weight subject to state and control inequality constraints and to terminal boundary conditions. This study shows that optimizing a control variable parametrically by using third-degree spline function interpolation allows the control to be shaped so that inequality constraints are strictly adhered to and all corners are eliminated. The absence of corners, which is realistic in nature, makes this method attractive from the viewpoint of solid rocket grain design.

Johnson, I. L.↗
NASA NTRS1983

Coolant passage heat transfer with rotation. A progress report on the computational aspects

Turbine airfoils are subjected to increasingly higher heat loads which escalate the cooling requirements in order to satisfy life goals for the component materials. If turbine efficiency is to be maintained, however, cooling requirements should be as low as possible. To keep the quantity of cooling air bounded, a more efficient internal cooling scheme must be developed. One approach is to employ airfoils with multipass cooling passages that contain devices to augment internal heat transfer while limiting pressure drop. Design experience with multipass cooling passage airfoils has shown that a surplus of cooling air must be provided as a margin of safety. This increased cooling air leads to a performance penalty. Reliable methods for predicting the internal thermal and aerodynamic performance of multipass cooling passage airfoils would reduce or eliminate the need for the safety margin of surplus cooling air. The objective of the program is to develop and verify improved analytical methods that will form the basis for design technology which will result in efficient turbine components with improved durability without sacrificing performance. The objective will be met by: (1) establishing a comprehensive experimental data base that can form the basis of an empirical design system; (2) developing computational fluid dynamic techniques; and (3) analyzing the information in the data base with both phenomenological modeling and mathematical modeling to derive a suitable design and analysis procedure.

Aceto, L. D.↗
NASA NTRS1990

Large scale nonlinear numerical optimal control for finite element models of flexible structures

This paper discusses the development of large scale numerical optimal control algorithms for nonlinear systems and their application to finite element models of structures. This work is based on our expansion of the optimal control algorithm (DDP) in the following steps: improvement of convergence for initial policies in non-convex regions, development of a numerically accurate penalty function method approach for constrained DDP problems, and parallel processing on supercomputers. The expanded constrained DDP algorithm was applied to the control of a four-bay, two dimensional truss with 12 soft members, which generates geometric nonlinearities. Using an explicit finite element model to describe the structural system requires 32 state variables and 10,000 time steps. Our numerical results indicate that for constrained or unconstrained structural problems with nonlinear dynamics, the results obtained by our expanded constrained DDP are significantly better than those obtained using linear-quadratic feedback control.

Shoemaker, Christine A.↗
NASA NTRS2000

Design of the Active Elevon Rotor for Low Vibration

Helicopter fuselages vibrate more than desired, and traditional solutions have limited effectiveness and can impose an appreciable weight penalty. Alternative methods of combating high vibration, including Higher Harmonic Control (HHC) via harmonic swashplate motion and Individual Blade Control (IBC) via active pitch links, have been studied for several decades. HHC via an on-blade control surface was tested in 1977 on a full scale rotor using a secondary active swashplate and a mechanical control system. Recent smart material advances have prompted new research into the use of on-blade control concepts. Recent analytical studies have indicated that the use of on-blade control surfaces produces vibration reduction comparable to swashplate-based HHC but for less power. Furthermore, smart materials (such as piezoceramics) have been shown to provide sufficient control authority for preliminary rotor experiments. These experiments were initially performed at small scale for reduced tip speeds. More recent experiments have been conducted at or near full tip speeds, and a full-scale active rotor is under development by Boeing with Eurocopter et al. pursuing a similarly advanced full-scale implementation. The US Army Aeroflightdynamics Directorate has undertaken a new research program called the Active Elevon Rotor (AER) Focus Demo. This program includes the design, fabrication, and wind. tunnel testing of a four-bladed, 12.96 ft diameter rotor with one or two on-blade elevons per blade. The rotor, which will be Mach scaled, will use 2-5/rev elevon motion for closed-loop control and :will be tested in late 2001. The primary goal of the AER Focus Demo is the reduction of vibratory hub loads by 80% and the reduction of vibratory blade structural loads. A secondary goal is the reduction of rotor power. The third priority is the measurement and possible reduction of Blade Vortex Interaction (BVI) noise. The present study is focused on elevon effectiveness, that is, the elevon's ability to reduce all six components of the nonrotating 4/rev hub loads. Some design parameters have been kept fixed in this study, while others have been varied to determine their influence - on elevon effectiveness. The fixed parameters include all blade structural properties except for torsion stiffness; the varied parameters include torsion stiffness, elevon aerodynamic location, and the number and individual authority of elevon aerodynamic surfaces. This paper describes the preliminary design process being used for the AER, and describes and quantifies the emerging active rotor characteristics.

Fulton, Mark V.↗
NASA NTRSDate not supplied

Low-cost Quantification of Fluid Flow Parameter Sensitivity using Reduced-order Modeling

Uncertainty quantification of computational fluid dynamics (CFD) simulations is a complicated procedure which still relies in many cases on engineering judgment and factors of safety. This is in part because the computational cost of measuring the simulation's sensitivity to all meaningful parameters (e.g., body surface roughness) and hyperparameters (e.g., subiteration convergence criterion) is intractable for even a single simulation. Reduced-order modeling dramatically lowers this computational cost of simulating fluid flows, but usually only where similar data is already available. In this work, fluid reduced-order models are utilized to quantify a flow's sensitivity to certain physical parameters for the purposes of improved uncertainty quantification. Characteristic observability and sensitivity are both explored. The resulting sensitivity quantifications enable more informed CFD frameworks and more rigorous uncertainty bounds on the resulting data.

uncertainty quantification↗
NASA NTRSDate not supplied

Low-cost Quantification of Fluid Flow Parameter Sensitivity using Reduced-order Modeling

Uncertainty quantification of computational fluid dynamics (CFD) simulations is a complicated procedure which still relies in many cases on engineering judgment and factors of safety. This is in part because the computational cost of measuring the simulation's sensitivity to all meaningful parameters (e.g., body surface roughness) and hyperparameters (e.g., subiteration convergence criterion) is intractable for even a single simulation. Reduced-order modeling dramatically lowers this computational cost of simulating fluid flows, but usually only where similar data is already available. In this work, fluid reduced-order models are utilized to quantify a flow's sensitivity to certain physical parameters for the purposes of improved uncertainty quantification. Characteristic observability and sensitivity are both explored. The resulting sensitivity quantifications enable more informed CFD frameworks and more rigorous uncertainty bounds on the resulting data.

uncertainty quantification↗
NASA NTRS1996

Implementation of Hybrid V-Cycle Multilevel Methods for Mixed Finite Element Systems with Penalty

The goal of this paper is the implementation of hybrid V-cycle hierarchical multilevel methods for the indefinite discrete systems which arise when a mixed finite element approximation is used to solve elliptic boundary value problems. By introducing a penalty parameter, the perturbed indefinite system can be reduced to a symmetric positive definite system containing the small penalty parameter for the velocity unknown alone. We stabilize the hierarchical spatial decomposition approach proposed by Cai, Goldstein, and Pasciak for the reduced system. We demonstrate that the relative condition number of the preconditioner is bounded uniformly with respect to the penalty parameter, the number of levels and possible jumps of the coefficients as long as they occur only across the edges of the coarsest elements.

Lai, Chen-Yao G.↗