Iterative solution of nonlinear optimal control problems.
Mathematical programming methods for solving nonlinear state-constrained discrete optimal control problems
SEARCH · Search NASA
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.
Mathematical programming methods for solving nonlinear state-constrained discrete optimal control problems
Hill-Brown differential equations for satellite coordinates transformed for integration with analytical programming language
A program for the measurement and analysis of the depolarization and differential attenuation that occur when millimeter wave signals propagate through rain is described. Initial data are taken along a 1.43 km path at 17.65 GHz and a supporting theoretical model is developed to relate the propagation effects to rainfall rate and wind velocity. A block diagram of the overall experiment is included. It consists of: (1) an RF system (millimeter wave transmitter and receiver), (2) transmitting and receiving antennas, (3) a weather system with rain gauges, wind sensors, and drop counters, and (4) a digital control, processing, and data storage system.
To understand the flame flash-back phenomena more extensively, an analytical model was formed and a numerical program was written and tested to solve the set of differential equations describing the model. Results show that under a given set of conditions flame propagates in the boundary layer on a flat plate when the free stream is at or below 1.8 m/s.
Barrier surface construction in the joint space of the differential turning game for thrust-vectored vs. conventional aircraft is discussed. Differential-turn studies are based on modifications of existing computer programs including an energy-turn program, and one which generates hodograph data. Optimal turning flight in energy approximation is discussed for the conventional aircraft configurations. It is concluded that any advantages realized from thrust-vectoring are minor, unless hover is possible, where advantages would be major at low energies, and affect tactics at high energies as well.
Low bit rates suitable for low-cost satellite communications. LOOP computer program written to simulate automatic frequency control (AFC) subsystem of differential minimum-shift keying (DMSK) receiver with bit rate of 2,400 baud. AFC simulated by LOOP is first-order loop configuration with first-order resistance-and-capacitance filter. Written in MS SuperSoft FORTRAN 77.
Report describes some path-following techniques for solution of nonlinear equations and compares with other methods. Use of multipurpose techniques applicable at more than one stage of path-following computation results in system relatively simple to understand, program, and use. Comparison of techniques with method of parametric differentiation (MPD) reveals definite advantages for path-following methods. Emphasis in investigation on multiuse techniques being applied at more than one stage of path-following computation. Incorporation of multipurpose techniques results in concise computer code relatively simple to use.
Automatic Differentiation (AD) is a tool that systematically implements the chain rule of differentiation to obtain the derivatives of functions calculated by computer programs. AD is assessed as a tool for engineering design. The forward and reverse modes of AD, their computing requirements, as well as approaches to implementing AD are discussed. The application of two different tools to two medium-size structural analysis problems to generate sensitivity information typically necessary in an optimization or design situation is also discussed. The observation is made that AD is to be preferred to finite differencing in most cases, as long as sufficient computer storage is available; in some instances, AD may be the alternative to consider in lieu of analytical sensitivity analysis.
Automatic Differentiation (AD) is a tool that systematically implements the chain rule of differentiation to obtain the derivatives of functions calculated by computer programs. In this paper, it is assessed as a tool for engineering design. The paper discusses the forward and reverse modes of AD, their computing requirements, and approaches to implementing AD. It continues with application to two different tools to two medium-size structural analysis problems to generate sensitivity information typically necessary in an optimization or design situation. The paper concludes with the observation that AD is to be preferred to finite differencing in most cases, as long as sufficient computer storage is available.
Genetic and evolutionary algorithms have been applied to solve numerous problems in engineering design where they have been used primarily as optimization procedures. These methods have an advantage over conventional gradient-based search procedures became they are capable of finding global optima of multi-modal functions and searching design spaces with disjoint feasible regions. They are also robust in the presence of noisy data. Another desirable feature of these methods is that they can efficiently use distributed and parallel computing resources since multiple function evaluations (flow simulations in aerodynamics design) can be performed simultaneously and independently on ultiple processors. For these reasons genetic and evolutionary algorithms are being used more frequently in design optimization. Examples include airfoil and wing design and compressor and turbine airfoil design. They are also finding increasing use in multiple-objective and multidisciplinary optimization. This lecture will focus on an evolutionary method that is a relatively new member to the general class of evolutionary methods called differential evolution (DE). This method is easy to use and program and it requires relatively few user-specified constants. These constants are easily determined for a wide class of problems. Fine-tuning the constants will off course yield the solution to the optimization problem at hand more rapidly. DE can be efficiently implemented on parallel computers and can be used for continuous, discrete and mixed discrete/continuous optimization problems. It does not require the objective function to be continuous and is noise tolerant. DE and applications to single and multiple-objective optimization will be included in the presentation and lecture notes. A method for aerodynamic design optimization that is based on neural networks will also be included as a part of this lecture. The method offers advantages over traditional optimization methods. It is more flexible than other methods in dealing with design in the context of both steady and unsteady flows, partial and complete data sets, combined experimental and numerical data, inclusion of various constraints and rules of thumb, and other issues that characterize the aerodynamic design process. Neural networks provide a natural framework within which a succession of numerical solutions of increasing fidelity, incorporating more realistic flow physics, can be represented and utilized for optimization. Neural networks also offer an excellent framework for multiple-objective and multi-disciplinary design optimization. Simulation tools from various disciplines can be integrated within this framework and rapid trade-off studies involving one or many disciplines can be performed. The prospect of combining neural network based optimization methods and evolutionary algorithms to obtain a hybrid method with the best properties of both methods will be included in this presentation. Achieving solution diversity and accurate convergence to the exact Pareto front in multiple objective optimization usually requires a significant computational effort with evolutionary algorithms. In this lecture we will also explore the possibility of using neural networks to obtain estimates of the Pareto optimal front using non-dominated solutions generated by DE as training data. Neural network estimators have the potential advantage of reducing the number of function evaluations required to obtain solution accuracy and diversity, thus reducing cost to design.
A pulsed, 2-m coherent Differential Absorption Lidar (DIAL) / Integrated Path Differential Absorption (IPDA) transceiver, developed under the Laser Risk Reduction Program (LRRP) at NASA, is integrated into a fully functional lidar instrument. This instrument will measure atmospheric CO2 profiles (by DIAL) initially from a ground platform, and then be prepared for aircraft installation to measure the atmospheric CO2 column densities in the atmospheric boundary layer (ABL) and lower troposphere. The airborne prototype CO2 lidar can measure atmospheric CO2 column density in a range bin of 1km with better than 1.5% precision at horizontal resolution of less than 50km. It can provide the image of the pooling of CO2 in lowlying areas and performs nighttime mass balance measurements at landscape scale. This sensor is unique in its capability to study the vertical ABL-free troposphere exchange of CO2 directly. It will allow the investigators to pursue subsequent in science-driven deployments, and provides a unique tool for Active Sensing of CO2 Emissions over Night, Days, and Seasons (ASCENDS) validation that was strongly advocated in the recent ASCENDS Workshop.
A unique method for solving the design problem of determining the areas required to meet an imposed pressure differential limit is presented. The associated analytical technique for determining the transient pressure differentials in a multicompartment series-parallel connected system is also presented. In an effort to verify the results obtained by using the analytical methods presented, a test was conducted using a seven compartment system. It was found in the investigations that a computer program based on the analytical model described can determine the pressure, temperature, pressure differentials, and mass flow rates into and out of each compartment at subsonic or sonic velocities.
A jointly planned U.S. Air Force-NASA program has been established to investigate the spacecraft charging phenomenon that has caused electronic anomalies in satellites in geosynchronous orbits. The objectives of this program are to provide design criteria, techniques, and test methods to insure control of absolute and differential charging of spacecraft surfaces. These objectives will be updated continuously over the next four years as data become available from the combined contractual and in-house programs. The geosynchronous altitude environment will be defined, ground and flight tests will be conducted, and materials and charge control techniques will be developed as required. The ultimate output of the program will be a spacecraft charging design criteria and test specification document. The program will be coordinated by a spacecraft charging program review group which has both Air Force and NASA representation.
Abstract Power system dynamics are generally modeled by high dimensional non‐linear differential‐algebraic equations (DAEs) given a large number of components forming the network. These DAEs' complexity can grow exponentially due to the increasing penetration of distributed energy resources, whereas their computation time becomes sensitive due to the increasing interconnection of the power grid with other energy systems. This paper demonstrates the use of quantum computing algorithms to solve DAEs for power system dynamic analysis. We leverage a symbolic programming framework to equivalently convert the power system's DAEs into ordinary differential equations (ODEs) using index reduction methods and then encode their data into qubits using amplitude encoding. The system non‐linearity is captured by Hamiltonian simulation with truncated Taylor expansion so that state variables can be updated by a quantum linear equation solver. Our results show that quantum computing can solve the power system's DAEs accurately with a computational complexity polynomial in the logarithm of the system dimension. We also illustrate the use of recent advanced tools in scientific machine learning for implementing complex computing concepts, that is, Taylor expansion, DAEs/ODEs transformation, and quantum computing solver with abstract representation for power engineering applications.
The CMS Inner Tracker, made of silicon pixel modules, will be entirely replaced prior to the start of the High Luminosity LHC period. One of the crucial components of the new Inner Tracker system is the readout chip, being developed by the RD53 Collaboration, and in particular its analogue front-end, which receives the signal from the sensor and digitizes it. Three different analogue front-ends (Synchronous, Linear, and Differential) were designed and implemented in the RD53A demonstrator chip. A dedicated evaluation program was carried out to select the most suitable design to build a radiation tolerant pixel detector able to sustain high particle rates with high efficiency and a small fraction of spurious pixel hits. The test results showed that all three analogue front-ends presented strong points, but also limitations. The Differential front-end demonstrated very low noise, but the threshold tuning became problematic after irradiation. Moreover, a saturation in the preamplifier feedback loop affected the return of the signal to baseline and thus increased the dead time. The Synchronous front-end showed very good timing performance, but also higher noise. For the Linear front-end all of the parameters were within specification, although this design had the largest time walk. This limitation was addressed and mitigated in an improved design. The analysis of the advantages and disadvantages of the three front-ends in the context of the CMS Inner Tracker operation requirements led to the selection of the improved design Linear front-end for integration in the final CMS readout chip.
This talk covers wide-ranging information on the Nuclear Criticality Safety Program (NCSP) Data Program. The Lecture includes slides on the history and mission of the NCSP along with charts and graphs on budgets from 2005 onwards. This lecture additionally mentions NCSP critical assemblies, differential experiments, and data measurements and evaluation work. The talk concludes with a gaze into NCSP benefits and successes.
The fuel optimal control problem arising in coplanar orbital transfer employing aeroassisted technology is addressed. The mission involves the transfer from high energy orbit (HEO) to low energy orbit (LEO) without plane change. The basic approach here is to employ a combination of propulsive maneuvers in space and aerodynamic maneuvers in the atmosphere. The basic sequence of events for the coplanar aeroassisted HEO to LEO orbit transfer consists of three phases. In the first phase, the transfer begins with a deorbit impulse at HEO which injects the vehicle into a elliptic transfer orbit with perigee inside the atmosphere. In the second phase, the vehicle is optimally controlled by lift and drag modulation to satisfy heating constraints and to exit the atmosphere with the desired flight path angle and velocity so that the apogee of the exit orbit is the altitude of the desired LEO. Finally, the second impulse is required to circularize the orbit at LEO. The performance index is maximum final mass. Simulation results show that the coplanar aerocapture is quite different from the case where orbital plane changes are made inside the atmosphere. In the latter case, the vehicle has to penetrate deeper into the atmosphere to perform the desired orbital plane change. For the coplanar case, the vehicle needs only to penetrate the atmosphere deep enough to reduce the exit velocity so the vehicle can be captured at the desired LEO. The peak heating rates are lower and the entry corridor is wider. From the thermal protection point of view, the coplanar transfer may be desirable. Parametric studies also show the maximum peak heating rates and the entry corridor width are functions of maximum lift coefficient. The problem is solved using a direct optimization technique which uses piecewise polynomial representation for the states and controls and collocation to represent the differential equations. This converts the optimal control problem into a nonlinear programming problem which is solved numerically by using a modified version of NPSOL. Solutions were obtained for the described problem for cases with and without heating constraints. The method appears to be more robust than other optimization methods. In addition, the method can handle complex dynamical constraints.
The derivation of the differential equations of motion of a 3 Degrees of Freedom body joined to a 3 Degrees of Freedom body by an elastic tether. The tether is represented by a spring and dashpot in parallel. A computer program which integrates the equations of motion is also described. Although the derivation of the equations of motions are for a general system, the computer program is written for defining loads in large boosters recovered by parachutes.