Some problems and some recent results in stochastic control.
Stochastic optimal control applications in dynamic programming for linear and nonlinear systems
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.
Stochastic optimal control applications in dynamic programming for linear and nonlinear systems
The following areas related to landing trajectory optimization research were discussed: (1) programming and modifying the steepest descent optimization procedure, (2) successfully iterating toward the optimum for a four-mile trajectory, (3) beginning optimization runs for a twenty-mile trajectory, and (4) adapt wind tunnel data for computer usage. Other related areas were discussed in detail in the two previous annual reports.
Framework for analysis and optimization problems (FRANOPP) is a software aid for the study and solution of design (optimization) problems which provides the driving program and plotting capability for a user generated programming system. In addition to FRANOPP, the programming system also contains the optimization code CONMIN, and two user supplied codes, one for analysis and one for output. With FRANOPP the user is provided with five options for studying a design problem. Three of the options utilize the plot capability and present an indepth study of the design problem. The study can be focused on a history of the optimization process or on the interaction of variables within the design problem.
An interactive computer program (AESOP) which solves quadratic optimal control and is discussed. The program can also be used to perform system analysis calculations such as transient and frequency responses, controllability, observability, etc., in support of the control and filter design computations.
Variable Metric Algorithm for Constrained Optimization (VMACO) is nonlinear computer program developed to calculate least value of function of n variables subject to general constraints, both equality and inequality. First set of constraints equality and remaining constraints inequalities. Program utilizes iterative method in seeking optimal solution. Written in ANSI Standard FORTRAN 77.
The following subject areas are covered: (1) the science program (star formation and origins of planetary systems; structure and evolution of the interstellar medium; stellar population; the galactic and extragalactic distance scale; nature of galaxy nuclei, AGNs, and QSOs; formation and evolution of galaxies at high redshifts; and cosmology); (2) implementation of the science program; (3) the observatory-class missions (HST; LST - the 6m successor to HST; and next-generation 16m telescope); (4) moderate and small missions (Delta-class Explorers; imaging astrometric interferometer; small Explorers; optics development and demonstrations; and supporting ground-based capabilities); (5) prerequisites - the current science program (Lyman-FUSE; HTS optimization; the near-term science program; data analysis, modeling, and theory funding; and archives); (6) technologies for the next century; and (7) lunar-based telescopes and instruments.
Iron deficiency is the leading nutritional deficiency in the U.S. and the rest of the world, with its highest prevalences in the developing world. Iron fortification of food has been proposed as a strategy to reduce the high prevalence of iron deficiency. Poor consumer acceptance, unacceptable taste, and discoloration of the iron-fortified foods have been frequently listed as causes of unsuccessful iron fortification programs. An excellent prospect for improving consumer acceptance of iron-fortified foods is the incorporation of a thorough, organized, and unified approach to sensory evaluation practices into iron fortification programs for product optimization. The information gained from systematic sensory evaluation allows for the manipulation of the sensory attributes, and thus improvement of the sensory properties of the fortified food. However, iron fortification programs have not systematically measured the effect of fortification on the sensory quality of the food. Because sensory evaluation is an important criterion in successful iron fortification, an integrated approach is necessary. Therefore, nutritionists and sensory scientists should work closely with each other to select the most suitable sensory tests and methods. The objectives of this article are to: (1) critically review and discuss some traditional and contemporary approaches and applications of sensory evaluation practices in iron fortification programs, and (2) demonstrate the importance of incorporating a multidisciplinary, systematic sensory evaluation approach in iron fortification programs.
A tool has been developed to address the growing Design for Manufacturing needs for composite structures, specifically those manufactured with automated fiber placement. This manufacturing approach presents unique challenges associated with puckers and wrinkling, tow overlaps and gaps, fiber deviation, and laminate strength. Achieving a satisfactory laminate design usually requires finding compromises between those four areas. The developed tool, dubbed the Central Optimizer, assists with the process of balancing competing design metrics in automated fiber placement. This tool was developed under the NASA Advanced Composites Consortium with input from industry partners. Under this program, the Central Optimizer has been exercised on three different structures to verify its functionality and validate ability to improve the design process for automated fiber placement structures.
Model-based design of experiments (MBDoE) is a powerful framework for selecting and calibrating science-based mathematical models from data. Here, this work extends popular MBDoE workflows by proposing a convex mixed integer (non)linear programming (MINLP) to optimize the selection of measurements. The solver MindtPy is modified to support calculating the D-optimality objective and its gradient via an external package, scipy, using the grey-box module in Pyomo. The new approach is demonstrated in two case studies: estimating highly correlated kinetics from a batch reactor and estimating transport parameters in a large-scale rotary packed bed for CO 2 capture. Both case studies show how examining the Pareto optimal trade-offs between information content measured by A- and D-optimality versus measurement budget offers practical guidance for selecting measurements for scientific experiments.
Computer program based on pair of algorithms selects approximately optimal locations of antennas and approximately optimal number of elements in each antenna of multiple-antenna communication system. Obscuration in field of view at given antenna location taken into account in choice of number of antenna elements mounted there. Directional coverages of combinations of up to four antenna elements computed in search for combination to cover clear portion of field of view. Developed to aid design of antenna system of conceptual space station. Applied to system aboard ship or aircraft, on building in city, or in any location where transmission and reception blocked in some directions from each potential antenna-mounting point.
A quadratic extended interior penalty function formulation especially well suited for second-order unconstrained optimization procedures is presented. Analytical derivatives of constraints and an approximate analysis technique are used. Minimum-mass design results are presented which indicate that the combination of these procedures can help make mathematical programming a useful optimization tool for large-order structural design problems with a large number of design variables and multiple constraints. Examples include statically loaded high- and low-aspect-ratio wings simultaneously subjected to stress, displacement, minimum gage and, in some cases, maximum strain constraints.
This report provides a description and performance evaluation of the optimization techniques for high performance computing (HPC) implementation of the open source Computational Materials mini-application Aladyn (https://github.com/nasa/aladyn). Aladyn is a basic molecular dynamics code written in FORTRAN 2003, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. While achieving orders of magnitude faster computational performance than DFT, the ANN-based approach was still very computationally demanding compared to the conventional approach of using empirically fitted energy functions. After its initial development, Aladyn was evaluated and optimized by experts at the NASA Advanced Supercomputing (NAS) division to exploit modern supercomputer architectures. The code has been optimized for execution on multicore central processing units (CPUs), including Intel® Skylake microarchitecture, and on graphic accelerators, such as Nvidia® V100 graphic processing units (GPUs), using Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) programming interfaces. The optimization achieved a speedup of 4.7 times the baseline version on CPU performance and an additional 2.4 times on CPU+GPU performance. Atomistic computer simulations are a fundamental tool in materials research to model material properties form physics-based first principles. Atomic interaction, governed by Quantum Mechanics (QM) require sophisticated and highly computationally demanding mathematical models to calculate [1]. Classical methods use approximate functional forms, empirically fitted through a set of variable parameters to emulate atomic energies as direct functions of atomic coordinates [2]. While empirical potentials are computationally much simpler, allowing simulations of large-scale systems of up to a trillion (1012) atoms [3], they are substantially less accurate compared to quantum calculations and applicable only to very specific atomic configurations or predefined crystallographic phases. A recently suggested approach is to use heuristic machine learning methods [4], such as those based on Adaptive Neural Networks (ANNs) to predict atomic energies, after being trained on a sufficiently large database of QM-calculated structures [5,6]. This approach reduces significantly the computational complexity, allowing for simulations of orders of magnitude larger systems compared to QM-based methods without compromising accuracy. Still, compared to classical methods using empirical energy functions, ANN methods remain two- to three orders of magnitude more computationally demanding. Hence, the computational cost of simulations, together with the need for extensive training of ANNs, still makes the practical implementation of ANN-based methods quite challenging. The purpose of the Aladyn mini-application software [7], available as open source at https://github.com/nasa/aladyn, is to be a testbed for exploring possible optimization strategies to develop highly scalable parallel algorithms for ANN-based atomistic simulations. Aladyn is aimed at utilizing the architecture of the high-end modern highperformance computing (HPC) hardware based on multicore central processing units (CPUs) equipped with graphic processing unit (GPU) accelerators. Specifically, the goal is to optimize the performance on a single HPC compute node, before implementing scaling to multi-node parallelization using message passing interface (MPI). At the same time, the open source code of Aladyn can serve as a training model for students and professors in academia.
Report treats optimization and effectiveness separately. Report illustrates example of dynamic programming solution to system optimization. Computer algorithm has been developed to solve effectiveness problem and is included in report.
The Magnetospheric Multiscale Mission (MMS) is a NASA mission intended to make fundamental advancements in our understanding of the Earth's Magnetosphere. There are three processes that MMS will study including magnetic reconnection, charged particle acceleration, and turbulence. There are four phases in the nominal mission and this work addresses some of the outstanding issues in phase I. The nominal phase I orbit is 1.2 x 12 R(sub e) highly elliptic orbit with four spacecraft nominally forming a regular tetrahedron. In this paper we investigate the relative dynamics of the four MMS spacecraft about an assumed reference orbit. There are several tetrahedron dimensions required in Phase I of the mission and in this work we design optimal tetrahedrons for the 10 km baseline. The performance metric used in the optimization process is directly related to the science return, and is based on an extension of previous work performed by Glassmeier. The optimizer we use is a commercially available Sequential Quadratic Programming (SQP) routine. Multiple optimal solutions are found, and we characterize how the performance of the formation varies between different regions of the reference orbit.
The proposed U.S. High Speed Civil Transport (HSCT) will revolutionize commercial air travel by providing economical supersonic passenger service to destinations worldwide. Unlike the high-bypass turbofan engines that propel today's subsonic airliners, HSCT engines will have much higher jet exhaust speeds. Jet noise, caused by the turbulent mixing of high-speed exhaust with the surrounding air, poses a significant challenge for HSCT engine designers. To resolve this challenge, engineers have designed advanced mixer rejector nozzles that reduce HSCT jet noise to airport noise certification levels by entraining and mixing large quantities of ambient air with the engines' jet streams. Although this works well during the first several minutes of flight, far away from the airport, as the HSCT gains speed and climbs, poor ejector inlet recovery and ejector ram drag contribute to poor thrust, making it advantageous to turn off the ejector. Doing so prematurely, however, can cause unacceptable noise levels to propagate to the ground, even when the aircraft is many miles from the airport. This situation lends itself ideally to optimization, where the aircraft trajectory, throttle setting, and ejector setting can be varied (subject to practical aircraft constraints) to minimize the noise propagated to the ground. A method was developed at the NASA Lewis Research Center that employs a variation of the classic energy state approximation: a trajectory analysis technique historically used to minimize climb time or fuel burned in many aircraft problems. To minimize the noise on the ground at any given throttle setting, high aircraft altitudes are desirable; but the HSCT may either climb quickly to high altitudes using a high, noisy throttle setting or climb more slowly at a lower, quieter throttle setting. An optimizer has been programmed into NASA's existing aircraft and noise analysis codes to balance these options by dynamically choosing the best altitude-velocity path and throttle setting history. The noise level standard, or metric, used in the optimizer should be one that accurately reflects the subjective annoyance levels of ground-based observers under the flight path. A variety of noise metrics are available, many of which are practical for airport-vicinity noise certification. Unlike airport noise, however, the HSCT's climb noise will be characterized by relatively low noise levels, long durations, and low-frequency spectra. The noise metrics used in these calculations are based on the recommendations of researchers at the NASA Langley Research Center, who have correlated the flyover noise annoyance levels of actual laboratory subjects with a variety of measurements. Analysis of data from this optimizer has shown that significant reductions in noise may be obtained with trajectory optimization. And since throttling operations are performed in the subsonic portion of the climb path (where thrust is plentiful), only small penalties in HSCT range or fuel performance occur.
This document reports on the findings of an inventory of the educational and workforce development resources in a four-county area on the I-91 corridor in Connecticut. The overarching goal of this exercise is to help the local workforce leaders in Bridgeport, Connecticut to optimize their training programs to support the county's clean energy transition. The inventory generally finds that Fairfield County, where Bridgeport is located, lacks adequate education and training programs compared to other surrounding counties with similar populations. More than half of Fairfield County's programs occur through high school or career technical education programs, suggesting that there are minimal opportunities for workers who are not currently high school students. Much of the existing focus within Fairfield County is on general construction or plumbing and electrical skills. There is room to expand by offering more programs focusing on energy efficiency and renewable energy skills. The document outlines some potential next steps and questions for consideration for the local workforce leaders in Bridgeport. Bridgeport is a city located within Fairfield County in southwestern Connecticut. It is the largest city in Connecticut with an estimated population of 150,000. Through the Department of Energy-funded Communities LEAP program, NREL conducted an analysis of existing education and workforce development (EWD) programs located in Fairfield County that align with (or could support) the worker pipeline for occupations related to building energy efficiency (EE), renewable energy (RE), and clean energy manufacturing. The analysis also looked at EWD offerings in other counties in Connecticut as a point of comparison.
Research related to reliable aircraft design is summarized. Topics discussed include systems reliability optimization, failure detection algorithms, analysis of nonlinear filters, design of compensators incorporating time delays, digital compensator design, estimation for systems with echoes, low-order compensator design, descent-phase controller for 4-D navigation, infinite dimensional mathematical programming problems and optimal control problems with constraints, robust compensator design, numerical methods for the Lyapunov equations, and perturbation methods in linear filtering and control.
A systems model reflecting the design concepts of Satellite Power Stations (SPS) was developed. The model is of sufficient scope to include the interrelationships of the following major design parameters: the transportation to and between orbits; assembly of the SPS; and maintenance of the SPS. The systems model is composed of a set of equations that are nonlinear with respect to the system parameters and decision variables. The model determines a figure of merit from which alternative concepts concerning transportation, assembly, and maintenance of satellite power stations are studied. A hybrid optimization model was developed to optimize the system's decision variables. The optimization model consists of a random search procedure and the optimal-steepest descent method. A FORTRAN computer program was developed to enable the user to optimize nonlinear functions using the model. Specifically, the computer program was used to optimize Satellite Power Station system components.