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 217 records · Page 12

An automated water iodinating subsystem for manned space flight

Controlling microbial growth by injecting iodine (l2) into water supplies is a widely acceptable technique, but requires a specialized injection method for space flight. An electrochemical l2 injection method and l2 level monitor are discussed in this paper, which also describe iodination practices previously used in the manned space program and major l2 biocidal characteristics. The development and design of the injector and monitor are described, and results of subsequent experiments are presented. Also presented are expected vehicle penalties for utilizing the l2 injector in certain space missions, especially the Space Shuttle, and possible injector failure modes and their criticality.

Houck, O. K.↗

Optimal detumbling of a large manned spacecraft using an internal moving mass

This investigation deals with the use of a movable mass control system to stabilize a tumbling asymmetric spacecraft about the maximum inertia axis. A first-order gradient optimization technique is used to minimize angular velocity components along the intermediate and minimum inertia axes, thus, permitting a wide range of initial guesses for mass position history. Motion of the control mass is along a linear track fixed in the vehicle. The control variable is taken as mass acceleration with respect to body coordinates. Motion is limited to defined quantities and a penalty function is used to insure a given range of positions. Numerical solutions of the optimization equations verify that minimum time detumbling is achieved with the largest permissible movable mass, length of linear track, and positions of the mass on the two coordinates perpendicular to the linear motion. The optimal method permits detumbling in about one-fourth the time when compared to a force control law formulation available in the literature.

Kunciw, B. G.↗

Design optimization of space structures

The topology-shape-size optimization of space structures is investigated through Kikuchi's homogenization method. The method starts from a 'design domain block,' which is a region of space into which the structure is to materialize. This domain is initially filled with a finite element mesh, typically regular. Force and displacement boundary conditions corresponding to applied loads and supports are applied at specific points in the domain. An optimal structure is to be 'carved out' of the design under two conditions: (1) a cost function is to be minimized, and (2) equality or inequality constraints are to be satisfied. The 'carving' process is accomplished by letting microstructure holes develop and grow in elements during the optimization process. These holes have a rectangular shape in two dimensions and a cubical shape in three dimensions, and may also rotate with respect to the reference axes. The properties of the perforated element are obtained through an homogenization procedure. Once a hole reaches the volume of the element, that element effectively disappears. The project has two phases. In the first phase the method was implemented as the combination of two computer programs: a finite element module, and an optimization driver. In the second part, focus is on the application of this technique to planetary structures. The finite element part of the method was programmed for the two-dimensional case using four-node quadrilateral elements to cover the design domain. An element homogenization technique different from that of Kikuchi and coworkers was implemented. The optimization driver is based on an augmented Lagrangian optimizer, with the volume constraint treated as a Courant penalty function. The optimizer has to be especially tuned to this type of optimization because the number of design variables can reach into the thousands. The driver is presently under development.

Felippa, Carlos↗

Primer vector theory and applications

A method developed to compute two-body, optimal, N-impulse trajectories was presented. The necessary conditions established define the gradient structure of the primer vector and its derivative for any set of boundary conditions and any number of impulses. Inequality constraints, a conjugate gradient iterator technique, and the use of a penalty function were also discussed.

Jezewski, D. J.↗

Optimal placement of tuning masses for vibration reduction in helicopter rotor blades

Described are methods for reducing vibration in helicopter rotor blades by determining optimum sizes and locations of tuning masses through formal mathematical optimization techniques. An optimization procedure is developed which employs the tuning masses and corresponding locations as design variables which are systematically changed to achieve low values of shear without a large mass penalty. The finite-element structural analysis of the blade and the optimization formulation require development of discretized expressions for two performance parameters: modal shaping parameter and modal shear amplitude. Matrix expressions for both quantities and their sensitivity derivatives are developed. Three optimization strategies are developed and tested. The first is based on minimizing the modal shaping parameter which indirectly reduces the modal shear amplitudes corresponding to each harmonic of airload. The second strategy reduces these amplitudes directly, and the third strategy reduces the shear as a function of time during a revolution of the blade. The first strategy works well for reducing the shear for one mode responding to a single harmonic of the airload, but has been found in some cases to be ineffective for more than one mode. The second and third strategies give similar results and show excellent reduction of the shear with a low mass penalty.

Pritchard, Jocelyn I.↗

A LEAST ABSOLUTE SHRINKAGE AND SELECTION OPERATOR (LASSO) FOR NONLINEAR SYSTEM IDENTIFICATION

Identification of parametric nonlinear models involves estimating unknown parameters and detecting its underlying structure. Structure computation is concerned with selecting a subset of parameters to give a parsimonious description of the system which may afford greater insight into the functionality of the system or a simpler controller design. In this study, a least absolute shrinkage and selection operator (LASSO) technique is investigated for computing efficient model descriptions of nonlinear systems. The LASSO minimises the residual sum of squares by the addition of a 1 penalty term on the parameter vector of the traditional 2 minimisation problem. Its use for structure detection is a natural extension of this constrained minimisation approach to pseudolinear regression problems which produces some model parameters that are exactly zero and, therefore, yields a parsimonious system description. The performance of this LASSO structure detection method was evaluated by using it to estimate the structure of a nonlinear polynomial model. Applicability of the method to more complex systems such as those encountered in aerospace applications was shown by identifying a parsimonious system description of the F/A-18 Active Aeroelastic Wing using flight test data.

Kukreja, Sunil L.↗

A variational method for finite element stress recovery: Applications in one-dimension

It is well-known that stresses (and strains) calculated by a displacement-based finite element analysis are generally not as accurate as the displacements. In addition, the calculated stress field is typically discontinuous at element interfaces. Because the stresses are typically of more interest than the displacements, several procedures have been proposed to obtain a smooth stress field, given the finite element stresses, and to improve the accuracy. Hinton and Irons introduced global least squares smoothing of discrete data defined on a plane using a finite element formulation. Tessler and co-workers recently developed a conceptually similar formulation for smoothing of two-dimensional data based on a discrete least square approximation with a penalty constraint. The penalty constraint results in a stress field which is C(exp 1)-continuous, a result not previously obtained. The approach requires additional, 'smoothing' finite element analysis and for their two-dimensional application, they used a conforming C(exp 0)-continuous triangular finite element based on a conforming plate element. This paper presents the results of a detailed investigation into the application of Tessler's smoothing procedure to the smoothing of finite element stresses from one-dimensional problems. Although the one-dimensional formulation has some practical applicability, such as in truss, beam, axisymmetric mechanics, and one-dimensional heat conduction, the primary motivation for developing the one-dimensional smoothing case is to explore the characteristics of the general smoothing strategy. In particular, it is used to describe the behavior of the method and to explore the suitability of criteria proposed for the smoothing analysis. Prior to presenting numerical results, the variational formulation of the smoothing strategy is presented and a criterion for the smoothing analysis is described.

Riggs, H. Ronald↗

Interior noise control prediction study for high-speed propeller-driven aircraft

An analytical model was developed to predict the noise levels inside propeller-driven aircraft during cruise at M = 0.8. The model was applied to three study aircraft with fuselages of different size (wide body, narrow body and small diameter) in order to determine the noise reductions required to achieve the goal of an A-weighted sound level which does not exceed 80 dB. The model was then used to determine noise control methods which could achieve the required noise reductions. Two classes of noise control treatments were investigated: add-on treatments which can be added to existing structures, and advanced concepts which would require changes to the fuselage primary structure. Only one treatment, a double wall with limp panel, provided the required noise reductions. Weight penalties associated with the treatment were estimated for the three study aircraft.

Rennison, D. C.↗

Implementation of natural frequency analysis and optimality criterion design

A description is presented of an effective optimality criterion computer design approach for member size selection to improve frequency characteristics for moderately large structure models. It is shown that the implementation of the simultaneous iteration method within a natural frequency structural design optimization provides a method which is more efficient in isolating the lowest natural frequency modes than the frequently applied Stodola method. Additional computational advantages are derived by using previously converged eigenvectors at the start of the iterations during the second and the following design cycles. Vectors with random components can be used at the first design cycle, which, in relation to the entire computer time for the design program, results in only a moderate computational penalty.

Levy, R.↗

Error and Complexity Analysis for a Collocation-Grid-Projection Plus Precorrected-FFT Algorithm for Solving Potential Integral Equations with LaPlace or Helmholtz Kernels

In this paper we derive error bounds for a collocation-grid-projection scheme tuned for use in multilevel methods for solving boundary-element discretizations of potential integral equations. The grid-projection scheme is then combined with a precorrected FFT style multilevel method for solving potential integral equations with 1/r and e(sup ikr)/r kernels. A complexity analysis of this combined method is given to show that for homogeneous problems, the method is order n natural log n nearly independent of the kernel. In addition, it is shown analytically and experimentally that for an inhomogeneity generated by a very finely discretized surface, the combined method slows to order n(sup 4/3). Finally, examples are given to show that the collocation-based grid-projection plus precorrected-FFT scheme is competitive with fast-multipole algorithms when considering realistic problems and 1/r kernels, but can be used over a range of spatial frequencies with only a small performance penalty.

Phillips, J. R.↗

A Poisson equation formulation for pressure calculations in penalty finite element models for viscous incompressible flows

The calculation of pressures when the penalty-function approximation is used in finite-element solutions of laminar incompressible flows is addressed. A Poisson equation for the pressure is formulated that involves third derivatives of the velocity field. The second derivatives appearing in the weak formulation of the Poisson equation are calculated from the C0 velocity approximation using a least-squares method. The present scheme is shown to be efficient, free of spurious oscillations, and accurate. Examples of applications are given and compared with results obtained using mixed formulations.

Sohn, J. L.↗

Theoretical Performance Characteristics of Sharp Lip Inlets at Subsonic Speeds

A method is presented for the estimation of the subsonic-flight-speed characteristics of sharp-lip inlets applicable to supersonic aircraft. The analysis, based on a simple momentum balance consideration, permits the computation of inlet pressure recovery - mass-flow relations and additive-drag coefficients for forward velocities from zero to the speed of sound. The penalties for operation of a sharp-lip inlet at velocity ratios other than 1.0 may be severe; at lower velocity ratios an additive drag is incurred that is not cancelled by lip suction, while at higher velocity ratios, unavoidable losses in inlet total pressure will result. In particular, at the take-off condition, the total pressure and the mass flow for a choked inlet are only 79 percent of the values ideally attainable with a rounded lip. Experimental data obtained at zero speed with a sharp-lip supersonic inlet model were in substantial agreement with the theoretical results.

Fradenburgh, Evan A↗

GCR Transport in the Brain: Assessment of Self-Shielding, Columnar Damage, and Nuclear Reactions on Cell Inactivation Rates

Radiation shield design is driven by the need to limit radiation risks while optimizing risk reduction with launch mass/expense penalties. Both limitation and optimization objectives require the development of accurate and complete means for evaluating the effectiveness of various shield materials and body-self shielding. For galactic cosmic rays (GCR), biophysical response models indicate that track structure effects lead to substantially different assessments of shielding effectiveness relative to assessments based on LET-dependent quality factors. Methods for assessing risk to the central nervous system (CNS) from heavy ions are poorly understood at this time. High-energy and charge (HZE) ion can produce tissue events resulting in damage to clusters of cells in a columnar fashion, especially for stopping heavy ions. Grahn (1973) and Todd (1986) have discussed a microlesion concept or model of stochastic tissue events in analyzing damage from HZE's. Some tissues, including the CNS, maybe sensitive to microlesion's or stochastic tissue events in a manner not illuminated by either conventional dosimetry or fluence-based risk factors. HZE ions may also produce important lateral damage to adjacent cells. Fluences of high-energy proton and alpha particles in the GCR are many times higher than HZE ions. Behind spacecraft and body self-shielding the ratio of protons, alpha particles, and neutrons to HZE ions increases several-fold from free-space values. Models of GCR damage behind shielding have placed large concern on the role of target fragments produced from tissue atoms. The self-shielding of the brain reduces the number of heavy ions reaching the interior regions by a large amount and the remaining light particle environment (protons, neutrons, deuterons. and alpha particles) may be the greatest concern. Tracks of high-energy proton produce nuclear reactions in tissue, which can deposit doses of more than 1 Gv within 5 - 10 cell layers. Information on rates of cell killing from GCR, including patterns of cell killing from single particle tracks. can provide useful information on expected differences between proton and HZE tracks and clinical experiences with photon irradiation. To model effects on cells in the brain, it is important that transport models accurately describe changes in the GCR due to interactions in the cranium and proximate tissues. We describe calculations of the attenuated GCR particle fluxes at three dose-points in the brain and associated patterns of cell killing using biophysical models. The effects of the brain self-shielding and bone-tissue interface of the skull in modulating the GCR environment are considered. For each brain dose-point, the mass distribution in the surrounding 4(pi) solid angle is characterized using the CAM model to trace 512 rays. The CAM model describes the self-shielding by converting the tissue distribution to mass-equivalent aluminum, and nominal values of spacecraft shielding is considered. Particle transport is performed with the proton, neutron, and heavy-ion transport code HZETRN with the nuclear fragmentation model QMSFRG. The distribution of cells killed along the path of individual GCR ions is modeled using in vitro cell inactivation data for cells with varying sensitivity. Monte Carlo simulations of arrays of inactivated cells are considered for protons and heavy ions and used to describe the absolute number of cell killing events of various magnitude in the brain from the GCR. Included are simulations of positions of inactivated cells from stopping heavy ions and nuclear stars produced by high-energy ions most importantly, protons and neutrons.

Shavers, M. R.↗

Theoretical performance characteristics of sharp-lip inlets at subsonic speeds

A method is presented for the estimation of the subsonic-flight-speed characteristics of sharp-lip inlets applicable to supersonic aircraft. The analysis, based on a simple momentum balance consideration, permits the computation of inlet-pressure-recovery mass-flow relations and additive-drag coefficients for forward velocities from zero to the speed of sound. The penalties for operation of a sharp-lip inlet at velocity ratios other than 1.0 may be severe; at lower velocity ratios an additive drag is incurred that is not cancelled by lip suction, while at higher velocity ratios, unavoidable losses in inlet total pressure will result. In particular, at the take-off condition, the total pressure and the mass flow for a choked inlet are only 79 percent of the values ideally attainable with a rounded lip. The test specimens were polished and unnotched. The manufacturer of the material, the Aluminum Company of America, has made axial-load tests on 24S-T4 and 75S-T6 rod material. The test techniques used at the three laboratories are described in detail; the test results are presented and are compared with each other and with results obtained on unpolished sheet by the National Bureau of Standards. Experimental data obtained at zero speed with a sharp-lip supersonic inlet model were in substantial agreement with the theoretical results. (author)

Fradenburgh, Evan A↗

Supersonic cruise research aircraft structural studies: Methods and results

NASA Langley Research Center SCAR in-house structural studies are reviewed. In methods development, advances include a new system of integrated computer programs called ISSYS, progress in determining aerodynamic loads and aerodynamically induced structural loads (including those due to gusts), flutter optimization for composite and metal airframe configurations using refined and simplified mathematical models, and synthesis of active controls. Results given address several aspects of various SCR configurations. These results include flutter penalties on composite wing, flutter suppression using active controls, roll control effectiveness, wing tip ground clearance, tail size effect on flutter, engine weight and mass distribution influence on flutter, and strength and flutter optimization of new configurations. The ISSYS system of integrated programs performed well in all the applications illustrated by the results, the diversity of which attests to ISSYS' versatility.

Sobieszczanski-Sobieski, J.↗

Transform methods for precision continuum and control models of flexible space structures

An open loop optimal control algorithm is developed for general flexible structures, based on Laplace transform methods. A distributed parameter model of the structure is first presented, followed by a derivation of the optimal control algorithm. The control inputs are expressed in terms of their Fourier series expansions, so that a numerical solution can be easily obtained. The algorithm deals directly with the transcendental transfer functions from control inputs to outputs of interest, and structural deformation penalties, as well as penalties on control effort, are included in the formulation. The algorithm is applied to several structures of increasing complexity to show its generality.

Lupi, Victor D.↗

Performance and Flowfield Measurements on a 10-inch Ducted Rotor VTOL UAV

A ducted fan VTOL UAV with a 10-inch diameter rotor was tested in the US Army 7-by 10-Foot Wind Tunnel. The test conditions covered a range of angle of attack from 0 to 110 degrees to the freestream. The tunnel velocity was varied from 0 (simulating a hover condition) to 128 ft/sec in propeller mode. A six-component internal balance measured the aerodynamic loads for a range of model configurations. including the isolated rotor, the isolated duct, and the full configuration of the duct and rotor. For some conditions, hotwire velocity surveys were conducted along the inner and outer surface of the duct and across the downstream wake. In addition, fluorescent oil flow visualization allowed the flow separation patterns inside and outside of the duct to be mapped for a few test conditions. Two different duct shapes were tested to determine the performance effects of leading edge radius. For each duct, a range of rotor tip gap from 1%R to 4.5%R was tested to determine the performance penalty in hover and axial flight. Measured results are presented in terms of hover performance, hover performance in a crosswind, and high angle of attack performance in propeller mode. In each case, the effects of both tip gap and duct leading edge radius are illustrated using measurements. Some of the hover performance issues were also studied using a simple analytical method, and the results agreed with the measurements.

Martin, Preston↗

The cost of applying current helicopter external noise reduction methods while maintaining realistic vehicle performance

Analytical methods were developed and/or adopted for calculating helicopter component noise, and these methods were incorporated into a unified total vehicle noise calculation model. Analytical methods were also developed for calculating the effects of noise reduction methodology on helicopter design, performance, and cost. These methods were used to calculate changes in noise, design, performance, and cost due to the incorporation of engine and main rotor noise reduction methods. All noise reduction techniques were evaluated in the context of an established mission performance criterion which included consideration of hovering ceiling, forward flight range/speed/payload, and rotor stall margin. The results indicate that small, but meaningful, reductions in helicopter noise can be obtained by treating the turbine engine exhaust duct. Furthermore, these reductions do not result in excessive life cycle cost penalties. Currently available main rotor noise reduction methodology, however, is shown to be inadequate and excessively costly.

Bowes, M. A.↗