Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained optimization”

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 559 records · Page 31

Computational problems in autoregressive moving average (ARMA) models

The choice of the sampling interval and the selection of the order of the model in time series analysis are considered. Band limited (up to 15 Hz) random torque perturbations are applied to the human ankle joint. The applied torque input, the angular rotation output, and the electromyographic activity using surface electrodes from the extensor and flexor muscles of the ankle joint are recorded. Autoregressive moving average models are developed. A parameter constraining technique is applied to develop more reliable models. The asymptotic behavior of the system must be taken into account during parameter optimization to develop predictive models.

Agarwal, G. C.↗

Temporal Changes in Astronauts’ Muscle and Cardiorespiratory Physiology Before, During, and After Spaceflight

NASA’s planned space exploration missions will require astronauts to safely perform extravehicular activity (EVA) and to safely egress vehicles in a variety of landing scenarios. Prolonged exposure to spaceflight decreases cardiovascular and sensorimotor function, causes loss of bone mineral density, reduces muscle mass and strength, and ultimately diminishes tolerance for physical activity. Although exercise can help mitigate these spaceflight-induced physiological decrements, little is known regarding the time-course of changes in muscle and aerobic performance during spaceflight. Furthermore, these exercise countermeasures are not fully protective. For example, maximal aerobic capacity (VO2pk), lower body muscle cross-sectional area, and strength all decrease by about 10% to 15% after long-duration missions on the International Space Station (ISS). Future long-duration space missions beyond Low Earth Orbit will employ exploration vehicles with less robust exercise hardware and more constrained exercise capabilities (e.g., less operational volume, less active exercise time) than provided on the ISS. Thus, countermeasures will need to be optimized to protect crew health and performance on exploration-class missions lasting up to 3 years. This requires a more detailed understanding of the dynamic effects of spaceflight on human health and performance, the ability of exercise to protect against this deconditioning, and the interaction of exercise with interrelated factors like nutrition, sleep, and environmental conditions.

N.C. Strock↗

Temporal Changes in Astronauts’ Muscle and Cardiorespiratory Physiology Before, During, and After Spaceflight

NASA’s planned space exploration missions will require astronauts to safely perform extravehicular activity (EVA) and to safely egress vehicles in a variety of landing scenarios. Prolonged exposure to spaceflight decreases cardiovascular and sensorimotor function, causes loss of bone mineral density, reduces muscle mass and strength, and ultimately diminishes tolerance for physical activity. Although exercise can help mitigate these spaceflight-induced physiological decrements, little is known regarding the time-course of changes in muscle and aerobic performance during spaceflight. Furthermore, these exercise countermeasures are not fully protective. For example, maximal aerobic capacity (VO2pk), lower body muscle cross-sectional area, and strength all decrease by about 10% to 15% after long-duration missions on the International Space Station (ISS). Future long-duration space missions beyond Low Earth Orbit will employ exploration vehicles with less robust exercise hardware and more constrained exercise capabilities (e.g., less operational volume, less active exercise time) than provided on the ISS. Thus, countermeasures will need to be optimized to protect crew health and performance on exploration-class missions lasting up to 3 years. This requires a more detailed understanding of the dynamic effects of spaceflight on human health and performance, the ability of exercise to protect against this deconditioning, and the interaction of exercise with interrelated factors like nutrition, sleep, and environmental conditions.

N. Strock↗

State-constrained booster trajectory solutions via finite elements and shooting

This paper presents an extension of a FEM formulation based on variational principles. A general formulation for handling internal boundary conditions and discontinuities in the state equations is presented, and the general formulation is modified for optimal control problems subject to state-variable inequality constraints. Solutions which only touch the state constraint and solutions which have a boundary arc of finite length are considered. Suitable shape and test functions are chosen for a FEM discretization. All element quadrature (equivalent to one-point Gaussian quadrature over each element) may be done in closed form. The final form of the algebraic equations is then derived. A simple state-constrained problem is solved. Then, for a practical application of the use of the FEM formulation, a launch vehicle subject to a dynamic pressure constraint (a first-order state inequality constraint) is solved. The results presented for the launch-vehicle trajectory have some interesting features, including a touch-point solution.

Bless, Robert R.↗

An algorithm for the solution of dynamic linear programs

The algorithm's objective is to efficiently solve Dynamic Linear Programs (DLP) by taking advantage of their special staircase structure. This algorithm constitutes a stepping stone to an improved algorithm for solving Dynamic Quadratic Programs, which, in turn, would make the nonlinear programming method of Successive Quadratic Programs more practical for solving trajectory optimization problems. The ultimate goal is to being trajectory optimization solution speeds into the realm of real-time control. The algorithm exploits the staircase nature of the large constraint matrix of the equality-constrained DLPs encountered when solving inequality-constrained DLPs by an active set approach. A numerically-stable, staircase QL factorization of the staircase constraint matrix is carried out starting from its last rows and columns. The resulting recursion is like the time-varying Riccati equation from multi-stage LQR theory. The resulting factorization increases the efficiency of all of the typical LP solution operations over that of a dense matrix LP code. At the same time numerical stability is ensured. The algorithm also takes advantage of dynamic programming ideas about the cost-to-go by relaxing active pseudo constraints in a backwards sweeping process. This further decreases the cost per update of the LP rank-1 updating procedure, although it may result in more changes of the active set that if pseudo constraints were relaxed in a non-stagewise fashion. The usual stability of closed-loop Linear/Quadratic optimally-controlled systems, if it carries over to strictly linear cost functions, implies that the saving due to reduced factor update effort may outweigh the cost of an increased number of updates. An aerospace example is presented in which a ground-to-ground rocket's distance is maximized. This example demonstrates the applicability of this class of algorithms to aerospace guidance. It also sheds light on the efficacy of the proposed pseudo constraint relaxation scheme.

Psiaki, Mark L.↗

A parametric sensitivity study for single-stage-to-orbit hypersonic vehicles using trajectory optimization

The class of hypersonic vehicle configurations with single stage-to-orbit (SSTO) capability reflect highly integrated airframe and propulsion systems. These designs are also known to exhibit a large degree of interaction between the airframe and engine dynamics. Consequently, even simplified hypersonic models are characterized by tightly coupled nonlinear equations of motion. In addition, hypersonic SSTO vehicles present a major system design challenge; the vehicle's overall mission performance is a function of its subsystem efficiencies including structural, aerodynamic, propulsive, and operational. Further, all subsystem efficiencies are interrelated, hence, independent optimization of the subsystems is not likely to lead to an optimum design. Thus, it is desired to know the effect of various subsystem efficiencies on overall mission performance. For the purposes of this analysis, mission performance will be measured in terms of the payload weight inserted into orbit. In this report, a trajectory optimization problem is formulated for a generic hypersonic lifting body for a specified orbit-injection mission. A solution method is outlined, and results are detailed for the generic vehicle, referred to as the baseline model. After evaluating the performance of the baseline model, a sensitivity study is presented to determine the effect of various subsystem efficiencies on mission performance. This consists of performing a parametric analysis of the basic design parameters, generating a matrix of configurations, and determining the mission performance of each configuration. Also, the performance loss due to constraining the total head load experienced by the vehicle is evaluated. The key results from this analysis include the formulation of the sizing problem for this vehicle class using trajectory optimization, characteristics of the optimal trajectories, and the subsystem design sensitivities.

Lovell, T. Alan↗

A methodology based on reduced complexity algorithm for system applications using microprocessors

The paper considers a methodology on the analysis and design of a minimum mean-square error criterion linear system incorporating a tapped delay line (TDL) where all the full-precision multiplications in the TDL are constrained to be powers of two. A linear equalizer based on the dispersive and additive noise channel is presented. This microprocessor implementation with optimized power of two TDL coefficients achieves a system performance comparable to the optimum linear equalization with full-precision multiplications for an input data rate of 300 baud.

Yan, T. Y.↗

The Davidon-Fletcher-Powell penalty function method: A generalized iterative technique for solving parameter optimization problems

The Fletcher-Powell version of the Davidon variable metric unconstrained minimization technique is described. Equations that have been used successfully with the Davidon-Fletcher-Powell penalty function technique for solving constrained minimization problems and the advantages and disadvantages of using them are discussed. The experience gained in the behavior of the method while iterating is also related.

Johnson, I. L., Jr.↗

A real-time approximate optimal guidance law for flight in a plane

A real-time guidance scheme is presented for the problem of maximizing the payload into orbit subject to the equations of motion of a rocket over a nonrotating spherical earth. The flight is constrained to a path in the equatorial plane while reaching an orbital altitude at orbital injection speeds. The dynamics of the problem can be separated into primary and perturbation effects by a small parameter, epsilon, which is the ratio of the atmospheric scale height to the radius of the earth. The Hamilton-Jacobi-Bellman or dynamic programming equation is expanded in an asymptotic series where the zeroth-order term (epsilon = 0) can be obtained in closed form. The neglected perturbation terms are included in the higher-order terms of the expansion, which are determined from the solution of first-order linear partial differential equations requiring only integrations which are quadratures. The quadratures can be performed rapidly with emerging computer capability, so that real-time approximate optimization can be used to construct the launch guidance law. The application of this technique to flight in three-dimensions is made apparent from the solution presented.

Feeley, Timothy S.↗

Ankle and hip postural strategies defined by joint torques

Previous studies have identified two discrete strategies for the control of posture in the sagittal plane based on EMG activations, body kinematics, and ground reaction forces. The ankle strategy was characterized by body sway resembling a single-segment-inverted pendulum and was elicited on flat support surfaces. In contrast, the hip strategy was characterized by body sway resembling a double-segment inverted pendulum divided at the hip and was elicited on short or compliant support surfaces. However, biomechanical optimization models have suggested that hip strategy should be observed in response to fast translations on a flat surface also, provided the feet are constrained to remain in contact with the floor and the knee is constrained to remain straight. The purpose of this study was to examine the experimental evidence for hip strategy in postural responses to backward translations of a flat support surface and to determine whether analyses of joint torques would provide evidence for two separate postural strategies. Normal subjects standing on a flat support surface were translated backward with a range of velocities from fast (55 cm/s) to slow (5 cm/s). EMG activations and joint kinematics showed pattern changes consistent with previous experimental descriptions of mixed hip and ankle strategy with increasing platform velocity. Joint torque analyses revealed the addition of a hip flexor torque to the ankle plantarflexor torque during fast translations. This finding indicates the addition of hip strategy to ankle strategy to produce a continuum of postural responses. Hip torque without accompanying ankle torque (pure hip strategy) was not observed. Although postural control strategies have previously been defined by how the body moves, we conclude that joint torques, which indicate how body movements are produced, are useful in defining postural control strategies. These results also illustrate how the biomechanics of the body can transform discrete control patterns into a continuum of postural corrections.

Non-NASA Center↗

Limitations of Separate Cloud and Rain Categories in Parameterizing Collision-Coalescence for Bulk Microphysics Schemes

Warm rain collision-coalescence has been persistently difficult to parameterize in bulk microphysics schemes. We use a flexible bulk microphysics scheme with bin scheme process parameterizations, called AMP, to investigate reasons for the difficulty. AMP is configured in a variety of ways to mimic bulk schemes and is compared to simulations with the bin scheme upon which AMP is built. We find that an important limitation in traditional bulk schemes is the use of separate cloud and rain categories. When the drop size distribution is instead represented by a continuous distribution, the simulation of cloud-to-rain conversion is substantially improved. We also find large sensitivity to the threshold size to distinguish cloud and rain in traditional schemes; substantial improvement is found by decreasing the threshold from 40 to 25 μm. Neither the use of an assumed functional form for the size distribution nor the choice of predicted distribution moments has a large impact on the ability of AMP to simulate rain production. When predicting four total moments of the liquid drop size distribution, either with a traditional two-category, two-moment scheme with a reduced size threshold, or a four-moment single-category scheme, errors in the evolution of mass and the cloud size distribution are similar, but the single-category scheme has a substantially better representation of the rain size distribution. Optimal moment combinations for the single-category approach are investigated and appear to be linked more to the information content they provide for constraining the size distributions than to their correlation with collision-coalescence rates.

autoconversion↗

Radar-Based Bayesian Estimation of Ice Crystal Growth Parameters within a Microphysical Model

The potential for polarimetric Doppler radar measurements to improve predictions of ice microphysical processes within an idealized model–observational framework is examined. In an effort to more rigorously constrain ice growth processes (e.g., vapor deposition) with observations of natural clouds, a novel framework is developed to compare simulated and observed radar measurements, coupling a bulk adaptive-habit model of vapor growth to a polarimetric radar forward model. Bayesian inference on key microphysical model parameters is then used, via a Markov chain Monte Carlo sampler, to estimate the probability distribution of the model parameters. The statistical formalism of this method allows for robust estimates of the optimal parameter values, along with (non-Gaussian) estimates of their uncertainty. To demonstrate this framework, observations from Department of Energy radars in the Arctic during a case of pristine ice precipitation are used to constrain vapor deposition parameters in the adaptive habit model. The resulting parameter probability distributions provide physically plausible changes in ice particle density and aspect ratio during growth. A lack of direct constraint on the number concentration produces a range of possible mean particle sizes, with the mean size inversely correlated to number concentration. Consistency is found between the estimated inherent growth ratio and independent laboratory measurements, increasing confidence in the parameter PDFs and demonstrating the effectiveness of the radar measurements in constraining the parameters. The combined Doppler and polarimetric observations produce the highest-confidence estimates of the parameter PDFs, with the Doppler measurements providing a stronger constraint for this case.

Robert S. Schrom↗

Structural optimization via a design space hierarchy

Mathematical programming techniques provide a general approach to automated structural design. An iterative method is proposed in which design is treated as a hierarchy of subproblems, one being locally constrained and the other being locally unconstrained. It is assumed that the design space is locally convex in the case of good initial designs and that the objective and constraint functions are continuous, with continuous first derivatives. A general design algorithm is outlined for finding a move direction which will decrease the value of the objective function while maintaining a feasible design. The case of one-dimensional search in a two-variable design space is discussed. Possible applications are discussed. A major feature of the proposed algorithm is its application to problems which are inherently ill-conditioned, such as design of structures for optimum geometry.

Vanderplaats, G. N.↗

LMI-Based Multiobjective Optimization and Control of Flexible Aircraft Using VCCTEF

This paper considers the control of coupled aeroelastic aircraft model with Variable Camber Continuous Trailing Edge Flap (VCCTEF) system. The relative motion between two adjacent flaps is constrained and this actuation constraint problem is converted into an output covariance constraint problem, and therefore can be formulated using linear matrix inequalities (LMIs). A set of LMI conditions is derived for the design of an observer-based dynamic output feedback controller for VCCTEF configured aeroelastic aircraft model. The proposed controller is then applied to the NASA Generic Transport Model (GTM) for simulation, and the results demonstrate the efficacy of the proposed approach.

LMI control↗

Dipole and quadrupole synthesis of electric potential fields

A general technique for expanding an unknown potential field in terms of a linear summation of weighted dipole or quadrupole fields is described. Computational methods were developed for the iterative addition of dipole fields. Various solution potentials were compared inside the boundary with a more precise calculation of the potential to derive optimal schemes for locating the singularities of the dipole fields. Then, the problem of determining solutions to Laplace's equation on an unbounded domain as constrained by pertinent electron trajectory data was considered.

Tilley, D. G.↗

Optimal estimation of undersea acoustic transponder locations

Using principles from multilateration and optimal estimation theories an approach is derived for estimating the relative positions of three or more submerged and anchored acoustic transponders. The procedure is not constrained to processing range data collected at special points or on special trajectories. While the data normally collected over transponders and between transponder pairs can be processed, simultaneous ranges from anywhere on the surface to three or more transponders can also be processed. Simulated examples involving four stations in different geometries with different range collection schemes demonstrate the effectiveness of the procedure.

Carta, D. G.↗