Search NASA⌕ Search

SEARCH · Search NASA

Results for “quadratic programming”

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.

232 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Modal methods in optimal control synthesis

Efficient algorithms for solving linear smoother-follower problems with quadratic criteria are presented. For time-invariant systems, the algorithm consists of one backward integration of a linear vector equation and one forward integration of another linear vector equation. Furthermore, the backward and forward Riccati matrices can be expressed in terms of the eigenvalues and eigenvectors of the Euler-Lagrange equations. Hence, the gains of the forward and backward Kalman-Bucy filters and of the optimal state-feedback regulator can be determined without integration of matrix Riccati equations. A computer program has been developed, based on this method of determining the gains, to synthesize the optimal time-invariant compensator in the presence of random disturbance inputs and random measurement errors. The program also computes the rms state and control variables of the optimal closed-loop system.

Bryson, A. E., Jr.↗
NASA NTRS1986

Mission analysis flow sequencing optimization

This investigation is an extension of a project dealing with the problem of optimal use of ground resources for future space missions. This problem was formulated as a linear programming problem using an indirect approach. Instead of minimizing the inventory level of needed ground resources, the overlapping periods during which the same types of resources are used by various flights are minimized. The model was built upon the assumption that during the time interval under consideration, the costs of various needed resources remain constant. Under other assumptions concerning costs of resources, the objective function, in general, assumes a non-linear form. In this study, one case where the form of the objective function turns out to be quadratic is considered. Also, disadvantages and limitations of the approach used are briefly discussed.

Scott, M.↗
NASA NTRS2001

Application of multivariate statistics to vestibular testing: discriminating between Meniere's disease and migraine associated dizziness

Meniere's disease (MD) and migraine associated dizziness (MAD) are two disorders that can have similar symptomatologies, but differ vastly in treatment. Vestibular testing is sometimes used to help differentiate between these disorders, but the inefficiency of a human interpreter analyzing a multitude of variables independently decreases its utility. Our hypothesis was that we could objectively discriminate between patients with MD and those with MAD using select variables from the vestibular test battery. Sinusoidal harmonic acceleration test variables were reduced to three vestibulo-ocular reflex physiologic parameters: gain, time constant, and asymmetry. A combination of these parameters plus a measurement of reduced vestibular response from caloric testing allowed us to achieve a joint classification rate of 91%, independent quadratic classification algorithm. Data from posturography were not useful for this type of differentiation. Overall, our classification function can be used as an unbiased assistant to discriminate between MD and MAD and gave us insight into the pathophysiologic differences between the two disorders.

NASA Discipline Neuroscience↗
NASA NTRS1980

An improved panel method for the solution of three-dimensional leading-edge vortex flows. Volume 1: Theory document

An improved panel method for the solution of three dimensional flow and wing and wing-body combinations with leading edge vortex separation is presented. The method employs a three dimensional inviscid flow model in which the configuration, the rolled-up vortex sheets, and the wake are represented by quadratic doublet distributions. The strength of the singularity distribution as well as shape and position of the vortex spirals are computed in an iterative fashion starting with an assumed initial sheet geometry. The method calculates forces and moments as well as detail surface pressure distributions. Improvements include the implementation of improved panel numerics for the purpose of elimination the highly nonlinear effects of ring vortices around double panel edges, and the development of a least squares procedure for damping vortex sheet geometry update instabilities. A complete description of the method is included. A variety of cases generated by the computer program implementing the method are presented which verify the mathematical assumptions of the method and which compare computed results with experimental data to verify the underlying physical assumptions made by the method.

Johnson, F. T.↗
NASA NTRS1987

SHARP simulation of discontinuities in highly convective steady flow

For steady multidimesional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗
NASA NTRS1993

The Middeck Active Control Experiment (MACE)

The Middeck Active Control Experiment (MACE) is a NASA In-Step and Control Structure Interaction (CSI) Office funded Shuttle middeck experiment. The objective is to investigate the extent to which closed-loop behavior of flexible spacecraft in zero-gravity (0-g) can be predicted. This prediction becomes particularly difficult when dynamic behavior during ground testing exhibits extensive suspension and direct gravity coupling. On-orbit system identification and control reconfiguration is investigated to improve performance which would otherwise be limited due to errors in prediction. The program is presently in its preliminary design phase with launch expected in the summer of 1994. The MACE test article consists of three attitude control torque wheels, a two axis gimballing payload, inertial sensors and a flexible support structure. With the acquisition of a second payload, this will represent a multiple payload platform with significant structural flexibility. This paper presents on-going work in the areas of modelling and control of the MACE test article in the zero and one-gravity environments. Finite element models, which include suspension and gravity effects, and measurement models, derived from experimental data, are used as the basis for Linear Quadratic Gaussian controller designs. Finite element based controllers are analytically used to study the differences in closed-loop performance as the test article transitions between the 0-g and 1-g environments. Measurement based controllers are experimentally applied to the MACE test article in the 1-g environment and achieve over an order of magnitude improvement in payload pointing accuracy when disturbed by a broadband torque disturbance. The various aspects of the flight portion of the experiment are also discussed.

Miller, David W.↗
NASA NTRS2000

Primate translational vestibuloocular reflexes. II. Version and vergence responses to fore-aft motion

To maintain binocular fixation on near targets during fore-aft translational disturbances, largely disjunctive eye movements are elicited the amplitude and direction of which should be tuned to the horizontal and vertical eccentricities of the target. The eye movements generated during this task have been investigated here as trained rhesus monkeys fixated isovergence targets at different horizontal and vertical eccentricities during 10 Hz fore-aft oscillations. The elicited eye movements complied with the geometric requirements for binocular fixation, although not ideally. First, the corresponding vergence angle for which the movement of each eye would be compensatory was consistently less than that dictated by the actual fixation parameters. Second, the eye position with zero sensitivity to translation was not straight ahead, as geometrically required, but rather exhibited a systematic dependence on viewing distance and vergence angle. Third, responses were asymmetric, with gains being larger for abducting and downward compared with adducting and upward gaze directions, respectively. As frequency was varied between 4 and 12 Hz, responses exhibited high-pass filter properties with significant differences between abduction and adduction responses. As a result of these differences, vergence sensitivity increased as a function of frequency with a steeper slope than that of version. Despite largely undercompensatory version responses, vergence sensitivity was closer to ideal. Moreover, the observed dependence of vergence sensitivity on vergence angle, which was varied between 2.5 and 10 MA, was largely linear rather than quadratic (as geometrically predicted). We conclude that the spatial tuning of eye velocity sensitivity as a function of gaze and viewing distance follows the general geometric dependencies required for the maintenance of foveal visual acuity. However, systematic deviations from ideal behavior exist that might reflect asymmetric processing of abduction/adduction responses perhaps because of different functional dependencies of version and vergence eye movement components during translation.

NASA Program Biomedical Research and Countermeasur↗
NASA NTRS1996

An optimized index of human cardiovascular adaptation to simulated weightlessness

Prolonged exposure to weightlessness is known to produce a variety of cardiovascular changes, some of which may influence the astronaut's performance during a mission. In order to find a reliable indicator of cardiovascular adaptation to weightlessness, we analyzed data from nine male subjects after a 24-hour period of normal activity and after a period of simulated weightlessness produced by two hours in a launch position followed by 20 hours of 6 degrees head-down tilt plus pharmacologically induced diuresis (furosemide). Heart rate, arterial pressure, thoracic fluid index, and radial flow were analyzed. Autoregressive spectral estimation and decomposition were used to obtain the spectral components of each variable from the subjects in the supine position during pre- and post-simulated weightlessness. We found a significant decrease in heart rate power and an increase in thoracic fluid index power in the high frequency region (0.2-0.45 Hz) and significant increases in radial flow and arterial pressure powers in the low frequency region (<0.2 Hz) in response to simulated weightlessness. However, due to the variability among subjects, any single variable appeared limited as a dependable index of cardiovascular adaptation to weightlessness. The backward elimination algorithm was then used to select the best discriminatory features from these spectral components. Fisher's linear discriminant and Bayes' quadratic discriminant were used to combine the selected features to obtain an optimal index of adaptation to simulated weightlessness. Results showed that both techniques provided improved discriminant performance over any single variable and thus have the potential for use as an index to track adaptation and prescribe countermeasures to the effects of weightlessness.

NASA Discipline Number 14-10↗
NASA NTRS2001

Integration of Libration Point Orbit Dynamics into a Universal 3-D Autonomous Formation Flying Algorithm

The autonomous formation flying control algorithm developed by the Goddard Space Flight Center (GSFC) for the New Millennium Program (NMP) Earth Observing-1 (EO-1) mission is investigated for applicability to libration point orbit formations. In the EO-1 formation-flying algorithm, control is accomplished via linearization about a reference transfer orbit with a state transition matrix (STM) computed from state inputs. The effect of libration point orbit dynamics on this algorithm architecture is explored via computation of STMs using the flight proven code, a monodromy matrix developed from a N-body model of a libration orbit, and a standard STM developed from the gravitational and coriolis effects as measured at the libration point. A comparison of formation flying Delta-Vs calculated from these methods is made to a standard linear quadratic regulator (LQR) method. The universal 3-D approach is optimal in the sense that it can be accommodated as an open-loop or closed-loop control using only state information.

Folta, David↗
NASA NTRS2017

QSPIN: A High Level Java API for Quantum Computing Experimentation

QSPIN is a high level Java language API for experimentation in QC models used in the calculation of Ising spin glass ground states and related quadratic unconstrained binary optimization (QUBO) problems. The Java API is intended to facilitate research in advanced QC algorithms such as hybrid quantum-classical solvers, automatic selection of constraint and optimization parameters, and techniques for the correction and mitigation of model and solution errors. QSPIN includes high level solver objects tailored to the D-Wave quantum annealing architecture that implement hybrid quantum-classical algorithms [Booth et al.] for solving large problems on small quantum devices, elimination of variables via roof duality, and classical computing optimization methods such as GPU accelerated simulated annealing and tabu search for comparison. A test suite of documented NP-complete applications ranging from graph coloring, covering, and partitioning to integer programming and scheduling are provided to demonstrate current capabilities.

Quantu↗
NASA NTRS1986

The design of a turboshaft speed governor using modern control techniques

The objectives of this program were: to verify the model of off schedule compressor variable geometry in the T700 turboshaft engine nonlinear model; to evaluate the use of the pseudo-random binary noise (PRBN) technique for obtaining engine frequency response data; and to design a high performance power turbine speed governor using modern control methods. Reduction of T700 engine test data generated at NASA-Lewis indicated that the off schedule variable geometry effects were accurate as modeled. Analysis also showed that the PRBN technique combined with the maximum likelihood model identification method produced a Bode frequency response that was as accurate as the response obtained from standard sinewave testing methods. The frequency response verified the accuracy of linear models consisting of engine partial derivatives and used for design. A power turbine governor was designed using the Linear Quadratic Regulator (LQR) method of full state feedback control. A Kalman filter observer was used to estimate helicopter main rotor blade velocity. Compared to the baseline T700 power turbine speed governor, the LQR governor reduced droop up to 25 percent for a 490 shaft horsepower transient in 0.1 sec simulating a wind gust, and up to 85 percent for a 700 shaft horsepower transient in 0.5 sec simulating a large collective pitch angle transient.

Delosreyes, G.↗
NASA NTRS1996

ARTSN: An Automated Real-Time Spacecraft Navigation System

As part of the Deep Space Network (DSN) advanced technology program an effort is underway to design a filter to automate the deep space navigation process.The automated real-time spacecraft navigation (ARTSN) filter task is based on a prototype consisting of a FORTRAN77 package operating on an HP-9000/700 workstation running HP-UX 9.05. This will be converted to C, and maintained as the operational version. The processing tasks required are: (1) read a measurement, (2) integrate the spacecraft state to the current measurement time, (3) compute the observable based on the integrated state, and (4) incorporate the measurement information into the state using an extended Kalman filter. This filter processes radiometric data collected by the DSN. The dynamic (force) models currently include point mass gravitational terms for all planets, the Sun and Moon, solar radiation pressure, finite maneuvers, and attitude maintenance activity modeled quadratically. In addition, observable errors due to troposphere are included. Further data types, force and observable models will be ncluded to enhance the accuracy of the models and the capability of the package. The heart of the ARSTSN is a currently available continuous-discrete extended Kalman filter. Simulated data used to test the implementation at various stages of development and the results from processing actual mission data are presented.

Burkhart, P. Daniel↗
NASA NTRS1988

Velocity-pressure integrated versus penalty finite element methods for high Reynolds number flows

Velocity-pressure integrated and consistent penalty finite element computations of high Reynolds number, laminar flows are presented. In both of the methods, the pressure has been interpolated using linear shape functions for a triangular element. The triangular element is contained inside the bi-quadratic isoparametric element. It has been reported previously that the pressure interpolation method, when used in the velocity-pressure integrated method, yielded accurate computational results for high Reynolds number flows. It is shown that use of the same pressure interpolation method in the consistent penalty finite element method yielded accurate velocity and pressure fields which were comparable to those obtained using the velocity-pressure integrated method. Accuracy of the two finite element methods has been demonstrated by comparing the computational results with available experimental data and/or fine-grid finite difference computational results. Advantages and disadvantages of the two methods are discussed on the basis of accuracy and convergence nature. Example problems considered include a lid-driven cavity flow for Reynolds number of 10,000, a laminar backward-facing step flow, a laminar flow through a nest of cylinders, and a channel flow with an internal blockage. A finite element computer program (NSFLOW/P) for the 2-D, incompressible Navier-Stokes equations is also presented.

Kim, Sang-Wook↗
NASA NTRS2003

Biological Effectiveness of Accelerated Particles for the Induction of Chromosome Damage Measured in Metaphase and Interphase Human Lymphocytes

Chromosome aberrations were investigated in human lymphocytes after in vitro exposure to 1H-, 3He-, 12C-, 40Ar-, 28Si-, 56Fe-, or 197Au-ion beams, with LET ranging from approximately 0.4-1393 keV/microm in the dose range of 0.075-3 Gy. Dose-response curves for chromosome exchanges, measured at the first mitosis postirradiation using fluorescence in situ hybridization (FISH) with whole-chromosome probes, were fitted with linear or linear-quadratic functions. The relative biological effectiveness (RBE) was estimated from the initial slope of the dose-response curve for chromosomal damage with respect to low- or high-dose-rate gamma rays. Estimates of RBEmax values for mitotic spreads, which ranged from near 0.7 to 11.1 for total exchanges, increased with LET, reaching a maximum at about 150 keV/microm, and decreased with further increase in LET. RBEs for complex aberrations are undefined due to the lack of an initial slope for gamma rays. Additionally, the effect of mitotic delay on RBE values was investigated by measuring chromosome aberrations in interphase after chemically induced premature chromosome condensation (PCC), and values were up to threefold higher than for metaphase analysis.

NASA Discipline Radiation Health↗
DOE OSTI2026

Classical combinatorial optimization scaling for random Ising models on 2D heavy-hex graphs

Motivated by near term quantum computing hardware limitations, combinatorial optimization problems that can be addressed by current quantum algorithms and noisy hardware with little or no overhead are used to probe capabilities of quantum algorithms such as the quantum approximate optimization algorithm. In this study, a specific class of near term quantum computing hardware defined combinatorial optimization problems, Ising models on heavy-hex graphs both with and without geometrically local cubic terms, are examined for their classical computational hardness via empirical computation time scaling quantification. Specifically the time-to-solution (TTS) metric using the classical heuristic simulated annealing is measured for finding optimal variable assignments (ground states), as well as the time required for the optimization software Gurobi to find an optimal variable assignment. Because of the sparsity of these Ising models, the classical algorithms are able to find optimal solutions efficiently even for large instances (i.e. 100 000 spin variables). The Ising models both with and without geometrically local cubic terms exhibit average-case linear-time or weakly quadratic scaling when solved exactly using Gurobi, and the Ising models with no cubic terms show evidence of exponential-time TTS scaling when sampled using simulated annealing. These findings point to the necessity of developing and testing more complex, namely more densely connected, optimization problems in order for quantum computing to ever have a practical advantage over classical computing. Our results are another illustration that different classical algorithms can indeed have exponentially different running times, thus making the identification of the best practical classical technique important in any quantum computing vs. classical computing comparison.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2016

Closed Loop Control Compact Exercise Device for Use on MPCV

Long duration space travel to Mars or to an asteroid will expose astronauts to extended periods of reduced gravity. To combat spaceflight physiological deconditioning, astronauts will use resistive and aerobic exercise regimens for the duration of the space flight to minimize the loss of bone density, muscle mass and aerobic capacity that occurs during exposure to a reduced gravity environment. Unlike the International Space Station (ISS), the mass and volume available for an exercise device in the next generation of spacecraft is limited. Therefore, compact exercise device prototypes are being developed for human in the loop evaluations. The NASA Human Research Program (HRP) is managing Advanced Exercise Concepts (AEC) requirements development and candidate technology maturation for all exploration mission profiles from Multi-Purpose Crew Vehicle (MPCV) exploration missions (e.g., EM-2, up to 21 day) to Mars Transit (up to 1000 day) missions. Numerous technologies have been considered and evaluated against HRP-approved functional requirements and include flywheel, pneumatic and closed-loop microprocessor-controlled motor driven power plants. Motor driven technologies offer excellent torque density and load accuracy characteristics as well as the ability to create custom mechanical impedance (the dynamic relationship between force and velocity) and custom load versus position exercise algorithms. Further, closed-loop motor-driven technologies offer the ability to monitor exercise dose parameters and adapt to the needs of the crewmember for real time optimization of exercise prescriptions. A simple proportional-integral-derivative (PID) controller is demonstrated in a prototype motor driven exercise device with comparison to resistive static and dynamic load set points and aerobic work rate targets. The resistive load term in the algorithm includes a constant force component (Fcmg) as well as inertial component (Fima) and a discussion of system tuning is presented in terms of addressing key functional requirements and human interfaces. The device aerobic modality is modelled as a rowing exercise using ground data sets obtained from Concept 2 rowers as well as competitive rowing1. A discussion of software and electronic implementations are presented which demonstrate unique approaches to meeting the constrained mass, volume and power requirements of the MPCV. . In addition to utilizing traditional PID control, controllers utilizing state feedback with gains solved using a Linear Quadratic Regulator will be developed. Controllability and observability will be utilized to investigate the need for state measurement in the design. As the control system directly interacts with human test subjects, robust methods such as H-infinity are also being investigated.1. Kleshnev V. Biomechanics. In: Rowing, Handbook of Sports Medicine and Science. ed. by Secher N., Voliantis S. IOC Medical Commission, Blackwell Pub. pp. 22-34, 2007

microgravity↗