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 55 records · Page 3

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↗

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

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

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

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

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

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

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↗

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↗

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

A strongly implicit procedure for the compressible Navier-Stokes equations

The strongly implicit procedure of Stone is used to obtain numerical solutions of the compressible Navier-Stokes equations in conservative form. In contrast to the spatially split Douglas-Gunn type methods, the method is shown to be numerically stable for the three-dimensional wave equation. The method is applied to a variety of external and internal two-dimensional flow problems involving shock wave boundary-layer interaction for both laminar and turbulent flows. The results are in good agreement with other methods and/or experiments. The storage penalty associated with the method is discussed and a simple, yet effective, means of minimizing the problem is presented.

Walters, R. W.↗

Using Statistical Penalties In The Tsai-Wu Failure Criterion

Improved methods of applying statistical penalties when using Tsai-Wu failure criterion lead to more accurate predictions of failures of composite-material structural components under stress, and provide better safety factors for designing such components. Intended to ensure proper use of statistical penalties with respect to failure hypersurface.

Richardson, D. E.↗

Robust Strategy for Rocket Engine Health Monitoring

Monitoring the health of rocket engine systems is essentially a two-phase process. The acquisition phase involves sensing physical conditions at selected locations, converting physical inputs to electrical signals, conditioning the signals as appropriate to establish scale or filter interference, and recording results in a form that is easy to interpret. The inference phase involves analysis of results from the acquisition phase, comparison of analysis results to established health measures, and assessment of health indications. A variety of analytical tools may be employed in the inference phase of health monitoring. These tools can be separated into three broad categories: statistical, rule based, and model based. Statistical methods can provide excellent comparative measures of engine operating health. They require well-characterized data from an ensemble of "typical" engines, or "golden" data from a specific test assumed to define the operating norm in order to establish reliable comparative measures. Statistical methods are generally suitable for real-time health monitoring because they do not deal with the physical complexities of engine operation. The utility of statistical methods in rocket engine health monitoring is hindered by practical limits on the quantity and quality of available data. This is due to the difficulty and high cost of data acquisition, the limited number of available test engines, and the problem of simulating flight conditions in ground test facilities. In addition, statistical methods incur a penalty for disregarding flow complexity and are therefore limited in their ability to define performance shift causality. Rule based methods infer the health state of the engine system based on comparison of individual measurements or combinations of measurements with defined health norms or rules. This does not mean that rule based methods are necessarily simple. Although binary yes-no health assessment can sometimes be established by relatively simple rules, the causality assignment needed for refined health monitoring often requires an exceptionally complex rule base involving complicated logical maps. Structuring the rule system to be clear and unambiguous can be difficult, and the expert input required to maintain a large logic network and associated rule base can be prohibitive.

Santi, L. Michael↗

Design and Test of an Improved Crashworthiness Small Composite Airframe

The purpose of this small business innovative research (SBIR) program was to evaluate the feasibility of developing small composite airplanes with improved crashworthiness. A combination of analysis and half scale component tests were used to develop an energy absorbing airframe. Four full scale crash tests were conducted at the NASA Impact Dynamics Research Facility, two on a hard surface and two onto soft soil, replicating earlier NASA tests of production general aviation airplanes. Several seat designs and restraint systems including both an air bag and load limiting shoulder harnesses were tested. Tests showed that occupant loads were within survivable limits with the improved structural design and the proper combination of seats and restraint systems. There was no loss of cabin volume during the events. The analysis method developed provided design guidance but time did not allow extending the analysis to soft soil impact. This project demonstrated that survivability improvements are possible with modest weight penalties. The design methods can be readily applied by airplane designers using the examples in this report.

Terry, James E.↗

Light transport and general aviation aircraft icing research requirements

A short term and a long term icing research and technology program plan was drafted for NASA LeRC based on 33 separate research items. The specific items listed resulted from a comprehensive literature search, organized and assisted by a computer management file and an industry/Government agency survey. Assessment of the current facilities and icing technology was accomplished by presenting summaries of ice sensitive components and protection methods; and assessments of penalty evaluation, the experimental data base, ice accretion prediction methods, research facilities, new protection methods, ice protection requirements, and icing instrumentation. The intent of the research plan was to determine what icing research NASA LeRC must do or sponsor to ultimately provide for increased utilization and safety of light transport and general aviation aircraft.

Breeze, R. K.↗

Simulation of Hypervelocity Impact on Aluminum-Nextel-Kevlar Orbital Debris Shields

An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.

Fahrenthold, Eric P.↗

Bounded Error Schemes for the Wave Equation on Complex Domains

This paper considers the application of the method of boundary penalty terms ("SAT") to the numerical solution of the wave equation on complex shapes with Dirichlet boundary conditions. A theory is developed, in a semi-discrete setting, that allows the use of a Cartesian grid on complex geometries, yet maintains the order of accuracy with only a linear temporal error-bound. A numerical example, involving the solution of Maxwell's equations inside a 2-D circular wave-guide demonstrates the efficacy of this method in comparison to others (e.g. the staggered Yee scheme) - we achieve a decrease of two orders of magnitude in the level of the L2-error.

Abarbanel, Saul↗