Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical optimal control”

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

Study of a Plug Nozzle for Supersonic Aircraft Concepts

Nozzles with external plugs are candidates for the exhaust system of supersonic aircraft concepts in an ongoing research effort at NASA. As part of that effort a fundamental study is being conducted at NASA Glenn Research Center (GRC). Experimental and limited numerical simulation results from this study, for a plug nozzle with a convergent nozzle (cowl),have been reported in recent publications[1,2]. Conventional method-of-characteristics (MoC) design for the plug requires that the flow is choked at the nozzle exit. This is why a sharp convergence of the cowl is invoked in such a design in most previous studies [3-5], leading to the earlier choice of the convergent cowl. Such a nozzle, however, is not suitable for supersonic flights because of large boat tail drag suffered by its outer surface. The boat tail drag consideration dictates that the cowl be more or less cylindrical in shape. While various other shapes have been studied in the past [6], essentially a cylindrical shape is adopted in the larger research program at NASA. The cylindrical shape has certain advantages. With the plug crown located somewhat inside the nozzle, the flow experiences convergence and divergence as with a C-D nozzle. By translating the plug, different throat-to-exit area ratio could be achieved to obtain fully expanded condition at different values of NPR. This offers a relatively easy active control strategy for the nozzle geometry throughout the flight regime from landing and takeoff (LTO) to cruise conditions. A plug nozzle model mimicking the geometry being considered in the NASA program, i.e., having a cylindrical exit, have been fabricated recently for continued fundamental studies. Some of the plugs used with the earlier convergent cowl [1] are readily adopted to this model. In addition, a few other plugs have been obtained following a numerical optimization study [7]. The various configurations are to be investigated for thrust, noise and flow fields, with complementing numerical simulations. Thrust data for a variety of configurations as well as limited noise data are in hand. Further explorations including schlieren visualization are to be conducted in the near future. The proposed paper is to summarize all these results. In the following, a few key results obtained so far are described.

jets↗

Optimal feedback control of turbulent channel flow

Feedback control equations were developed and tested for computing wall normal control velocities to control turbulent flow in a channel with the objective of reducing drag. The technique used is the minimization of a 'cost functional' which is constructed to represent some balance of the drag integrated over the wall and the net control effort. A distribution of wall velocities is found which minimizes this cost functional some time shortly in the future based on current observations of the flow near the wall. Preliminary direct numerical simulations of the scheme applied to turbulent channel flow indicates it provides approximately 17 percent drag reduction. The mechanism apparent when the scheme is applied to a simplified flow situation is also discussed.

Bewley, Thomas↗

A PCR Based Microbial Monitoring Alternative Method of Detection and Identification of Microbes Aboard ISS

Previous research has shown that microorganisms and potential human pathogens have been detected on the International Space Station (ISS) with additional introduction of new microflora occurring with every exchange of crew or addition of equipment and supplies. These microbes are readily transferred between crew and subsystems (i.e. ECLSS, environmental control and life support systems). As this can be detrimental to astronaut health and optimal performance of ISS systems, monitoring of systems such as ECLSS to include identification of microbial contaminants could prevent adverse effects on human health and life support systems. Current monitoring on ISS is laborious and utilizes culture based methods followed by sample return to Earth for complete analysis. Future, long-distance spaceflight missions will require real-time monitoring capabilities that enable efficient and rapid assessments of the microbial environment allowing for expedited decisions and more targeted response to cope with anomalies. Polymerase chain reaction (PCR), a molecular microbial monitoring method was chosen and numerous PCR instruments investigated for their potential to perform in microgravity conditions. Using ISS as a test bed for PCR verification in microgravity will enable NASA to assess whether molecular based microbiological sensors may be components of reliable, closed-loop life support and habitation systems in spacecraft, enhancing infrastructure capabilities through increased efficiency, reliability, and time savings by enabling sample analysis on orbit. NASA selected the Water Monitoring Suite as one of the rapid spaceflight hardware demonstration activities utilizing a streamlined process to minimize the time required to fly experimental flight hardware. The RAZOR EX (BioFire Defense, Salt Lake City, UT) system was part of the water monitoring suite and is a commercial off-the-shelf (COTS) real-time PCR instrument designed for field work. The RAZOR EX was originally designed for Department of Defense (DoD) under a small business innovative research (SBIR) grant and is ruggedized, compact and provides a rapid, sample to answer in less than an hour. PCR assays using a fluorescent probe were optimized and spiked with known concentrations of DNA (Pseudomonas aeruginosa) ranging from 0.002 to 20 ng. PCR reagents were lyophilized and configured in customized pouches and tested for flight readiness. Three types of water were used to rehydrate the reagents and demonstrate the fidelity of the PCR reaction in microgravity. Molecular grade deionized water served as a control while filtered and unfiltered ISS potable water served to test for chemical or biological inhibitors. All three types were compared to parallel ground test results. Nine tests were run on ISS (3 of each water type) and the critical threshold cycle (Ct) was compared to parallel ground tests completed at Kennedy Space Center, FL and Johnson Space Center, TX. All concentrations of Pseudomonas aeruginosa DNA were detected. A comparison of the Ct produced in real time PCR indicated similarity between flight and ground samples. There appeared to be no significant difference between flight or ground PCR reactions or between any of the three water types. This testing demonstrated the ability to perform molecular testing during spaceflight operations with similar sensitivity. It will allow for future ground development of molecular protocols and minimize the need for spaceflight testing. Future testing will include development of additional targets including environmental and health related organisms.

Christina Khodadad↗

DYCAST: A finite element program for the crash analysis of structures

DYCAST is a nonlinear structural dynamic finite element computer code developed for crash simulation. The element library contains stringers, beams, membrane skin triangles, plate bending triangles and spring elements. Changing stiffnesses in the structure are accounted for by plasticity and very large deflections. Material nonlinearities are accommodated by one of three options: elastic-perfectly plastic, elastic-linear hardening plastic, or elastic-nonlinear hardening plastic of the Ramberg-Osgood type. Geometric nonlinearities are handled in an updated Lagrangian formulation by reforming the structure into its deformed shape after small time increments while accumulating deformations, strains, and forces. The nonlinearities due to combined loadings are maintained, and stiffness variation due to structural failures are computed. Numerical time integrators available are fixed-step central difference, modified Adams, Newmark-beta, and Wilson-theta. The last three have a variable time step capability, which is controlled internally by a solution convergence error measure. Other features include: multiple time-load history tables to subject the structure to time dependent loading; gravity loading; initial pitch, roll, yaw, and translation of the structural model with respect to the global system; a bandwidth optimizer as a pre-processor; and deformed plots and graphics as post-processors.

Pifko, A. B.↗

Drag reduction strategies

previously a description was given of an active control scheme using wall transpiration that leads to a 15% reduction in surface skin friction beneath a turbulent boundary layer, according to direct numerical simulation. In this research brief further details of that scheme and its variants are given together with some suggestions as to how sensor/actuator arrays could be configured to reduce surface drag. The research which is summarized here was performed during the first half of 1994. This research is motivated by the need to understand better how the dynamics of near-wall turbulent flow can be modified so that skin friction is reduced. The reduction of turbulent skin friction is highly desirable in many engineering applications. Experiments and direct numerical simulations have led to an increased understanding of the cycle of turbulence production and transport in the boundary layer and raised awareness of the possibility of disrupting the process with a subsequent reduction in turbulent skin friction. The implementation of active feedback control in a computational setting is a viable approach for the investigation of the modifications to the flow physics that can be achieved. Bewley et al. and Hill describe how ideas from optimal control theory are employed to give 'sub-optimal' drag reduction schemes. The objectives of the work reported here is to investigate in greater detail the assumptions implicit within such schemes and their limitations. It is also our objective to describe how an array of sensors and actuators could be arranged and interconnected to form a 'smart' surface which has low skin friction.

Hill, D. Christopher↗

Single-Mode WGM Resonators Fabricated by Diamond Turning

A diamond turning process has made possible a significant advance in the art of whispering-gallery-mode (WGM) optical resonators. By use of this process, it is possible to fashion crystalline materials into WGM resonators that have ultrahigh resonance quality factors (high Q values), are compact (ranging in size from millimeters down to tens of microns), and support single electromagnetic modes. This development combines and extends the developments reported in "Few- Mode Whispering-Gallery-Mode Resonators" (NPO-41256), NASA Tech Briefs, Vol. 30, No. 1 (January 2006), page 16a and "Fabrication of Submillimeter Axisymmetric Optical Components" (NPO-42056), NASA Tech Briefs, Vol. 31, No. 5 (May 2007), page 10a. To recapitulate from the first cited prior article: A WGM resonator of this special type consists of a rod, made of a suitable transparent material, from which protrudes a thin circumferential belt of the same material. The belt is integral with the rest of the rod and acts as a circumferential waveguide. If the depth and width of the belt are made appropriately small, then the belt acts as though it were the core of a single-mode optical fiber: the belt and the rod material adjacent to it support a single, circumferentially propagating mode or family of modes. To recapitulate from the second cited prior article: A major step in the fabrication of a WGM resonator of this special type is diamond turning or computer numerically controlled machining of a rod of a suitable transparent crystalline material on an ultrahigh-precision lathe. During the rotation of a spindle in which the rod is mounted, a diamond tool is used to cut the rod. A computer program is used to control stepping motors that move the diamond tool, thereby controlling the shape cut by the tool. Because the shape can be controlled via software, it is possible to choose a shape designed to optimize a resonator spectrum, including, if desired, to limit the resonator to supporting a single mode. After diamond turning, a resonator can be polished to increase its Q. By virtue of its largely automated, computer-controlled nature, the process is suitable for mass production of nominally identical single-mode WGM resonators. In a demonstration of the capabilities afforded by this development, a number of WGM resonators of various designs were fabricated side by side on the surface of a single CaF2 rod (see figure).

Grudinin, Ivan↗

An Analysis of the Optimal Control Modification Method Applied to Flutter Suppression

Unlike basic Model Reference Adaptive Control (MRAC)l, Optimal Control Modification (OCM) has been shown to be a promising MRAC modification with robustness and analytical properties not present in other adaptive control methods. This paper presents an analysis of the OCM method, and how the asymptotic property of OCM is useful for analyzing and tuning the controller. We begin with a Lyapunov stability proof of an OCM controller having two adaptive gain terms, then the less conservative and easily analyzed OCM asymptotic property is presented. Two numerical examples are used to show how this property can accurately predict steady state stability and quantitative robustness in the presence of time delay, and relative to linear plant perturbations, and nominal Loop Transfer Recovery (LTR) tuning. The asymptotic property of the OCM controller is then used as an aid in tuning the controller applied to a large scale aeroservoelastic longitudinal aircraft model for flutter suppression. Control with OCM adaptive augmentation is shown to improve performance over that of the nominal non-adaptive controller when significant disparities exist between the controller/observer model and the true plant model.

aeroservoelasti↗

Application of numerical methods to extend capabilities for optimal rocket guidance: report on reentry guidance of shuttle orbiter

Analytical models are presented for optimal trajectories and reentry guidance of the space shuttle orbiter. Major emphasis is placed on the development of a "footprint', which refers to a set of reachable destination positions attainable by the shuttle at a specified terminal altitude. An unconstrained reentry footprint was calculated for a shuttle vehicle which enters the earth's atmosphere at 93 km initial altitude after a deboost from a near earth orbit. The method of computation is briefly described, and graphs are presented which illustrate the footprint and the variation of state and control variables along it. The effects of constraints and of variations in initial state upon the footprint are discussed.

Source record↗

Optimal discrete-time dynamic output-feedback design - A w-domain approach

An alternative method for optimal digital control design is described in this paper. The method is based on the usage of the w-transform and has many attractive design features. One of these is its immediate connection with frequency loop-shaping techniques that are now popular and effective for multivariable control synthesis in continuous-time domain. Furthermore, any design algorithms originally developed for continuous-time systems can now be immediately extended to the discrete-time domain. The main results presented in this paper are the exact problem formulation and solution of an optimal discrete-time dynamic output-feedback design in the w-domain involving a quadratic performance index to random disturbances. In addition, necessary conditions for optimality are obtained for the numerical solution of the optimal output-feedback compensator design. A numerical example is presented illustrating its application to the design of a low-order dynamic compensator in a stability augmentation system of a commercial transport.

Ha, Cheolkeun↗

The isolation limits of stochastic vibration

The vibration isolation problem is formulated as a 1D kinematic problem. The geometry of the stochastic wall trajectories arising from the stroke constraint is defined in terms of their significant extrema. An optimal control solution for the minimum acceleration return path determines a lower bound on platform mean square acceleration. This bound is expressed in terms of the probability density function on the significant maxima and the conditional fourth moment of the first passage time inverse. The first of these is found analytically while the second is found using a Monte Carlo simulation. The rms acceleration lower bound as a function of available space is then determined through numerical quadrature.

Knopse, C. R.↗

An analysis and comparison of several trajectory optimization methods

The sensitivities of the convergence characteristics of the methods to initially assumed parameters and trial solution, convergence times, computer logic, and storage requirements are discussed. Numerical comparison of the convergence characteristics is made by considering a minimum time, low thrust, Earth-Mars transfer trajectory. A modified quasi-linearization method reduces convergence time by approximately 70% when compared with the generalized Newton-Raphson method and allows the terminal boundary to be specified by a general function of the problem variables. A uniquely specified and easily determined, time dependent weighting matrix for the gradient techniques accelerates the shaping of the optimal control program and improves the convergence characteristics during the terminal iterations. Convergence envelopes, indicating how sensitive the convergence characteristics are to initially assumed parameters, are plotted for the perturbation and quasi-linearization methods. Several iteration schemes are proposed which increase the size of the convergence envelopes and decrease the sensitivity of the method to initially assumed parameters.

Lewallen, J. M.↗

Considerations of open-loop, closed-loop, and adaptive multicyclic control systems

Four different types of self-tuning regulators were studied for multicyclic control of helicopter vibration. A numerical simulation of the helicopter is made, using a multivariable frequency-domain model, in terms of transfer function with six input control harmonics and six output harmonics. The model characteristics vary with flight speed. An off-line identification of model characteristics is made, using the least-squared-error method and using a succession of input and output measurements. The on-line identification of model characteristics is made using the Kalman filter solution. The optimal controls are calculated from the minimization of quadratic performance function based on response and multicyclic inputs. The performance of various regulators or controllers is judged from the stability, transient response, convergence time, and amplitude of the steady state.

Chopra, I.↗

Optimal orbit transfer suitable for large flexible structures

The problem of continuous low-thrust planar orbit transfer of large flexible structures is formulated as an optimal control problem with terminal state constraints. The dynamics of the spacecraft motion are treated as a point-mass central force field problem; the thrust-acceleration magnitude is treated as an additional state variable; and the rate of change of thrust-acceleration is treated as a control variable. To ensure smooth transfer, essential for flexible structures, an additional quadratic term is appended to the time cost functional. This term penalizes any abrupt change in acceleration. Numerical results are presented for the special case of a planar transfer.

Chatterjee, Alok K.↗

HTMT-class Latency Tolerant Parallel Architecture for Petaflops Scale Computation

Computational Aero Sciences and other numeric intensive computation disciplines demand computing throughputs substantially greater than the Teraflops scale systems only now becoming available. The related fields of fluids, structures, thermal, combustion, and dynamic controls are among the interdisciplinary areas that in combination with sufficient resolution and advanced adaptive techniques may force performance requirements towards Petaflops. This will be especially true for compute intensive models such as Navier-Stokes are or when such system models are only part of a larger design optimization computation involving many design points. Yet recent experience with conventional MPP configurations comprising commodity processing and memory components has shown that larger scale frequently results in higher programming difficulty and lower system efficiency. While important advances in system software and algorithms techniques have had some impact on efficiency and programmability for certain classes of problems, in general it is unlikely that software alone will resolve the challenges to higher scalability. As in the past, future generations of high-end computers may require a combination of hardware architecture and system software advances to enable efficient operation at a Petaflops level. The NASA led HTMT project has engaged the talents of a broad interdisciplinary team to develop a new strategy in high-end system architecture to deliver petaflops scale computing in the 2004/5 timeframe. The Hybrid-Technology, MultiThreaded parallel computer architecture incorporates several advanced technologies in combination with an innovative dynamic adaptive scheduling mechanism to provide unprecedented performance and efficiency within practical constraints of cost, complexity, and power consumption. The emerging superconductor Rapid Single Flux Quantum electronics can operate at 100 GHz (the record is 770 GHz) and one percent of the power required by convention semiconductor logic. Wave Division Multiplexing optical communications can approach a peak per fiber bandwidth of 1 Tbps and the new Data Vortex network topology employing this technology can connect tens of thousands of ports providing a bi-section bandwidth on the order of a Petabyte per second with latencies well below 100 nanoseconds, even under heavy loads. Processor-in-Memory (PIM) technology combines logic and memory on the same chip exposing the internal bandwidth of the memory row buffers at low latency. And holographic storage photorefractive storage technologies provide high-density memory with access a thousand times faster than conventional disk technologies. Together these technologies enable a new class of shared memory system architecture with a peak performance in the range of a Petaflops but size and power requirements comparable to today's largest Teraflops scale systems. To achieve high-sustained performance, HTMT combines an advanced multithreading processor architecture with a memory-driven coarse-grained latency management strategy called "percolation", yielding high efficiency while reducing the much of the parallel programming burden. This paper will present the basic system architecture characteristics made possible through this series of advanced technologies and then give a detailed description of the new percolation approach to runtime latency management.

Sterling, Thomas↗

Solid Fuel Ignition and Extinction (SoFIE) Project on ISS

The Solid Fuel Ignition and Extinction (SoFIE) project studies ignition and flammability of solid spacecraft materials (fuels) in practical geometries and realistic atmospheric conditions. It is an experiment insert designed for use within the existing Combustion Integrated Rack (CIR). The CIR chamber provides a level of containment and permits testing at variable oxygen concentrations and pressures representative of current and planned NASA Space Exploration Atmospheres. The applications of SoFIE include: (1) Determining safer selection of cabin materials and validating NASA materials flammability selection using 1-g test protocols for low-gravity fires, (2) Improving understanding of early fire growth behavior, (3) Validating material flammability numerical models, (4) Determining optimal suppression techniques for burning materials by diluents, flow reduction, and venting, (5) Obtaining burning behavior of actual engineering materials planned for spacecraft, (6) Developing corresponding models of microgravity flame spread, flammability, and extinction, and use the results to improve normal gravity combustion models for terrestrial applications. The hardware permits a wide range of solid-material combustion and fire suppression studies. It supports multiple investigations using common infrastructure including sample holders, flow control, test sections, external radiant heaters, igniters, and diagnostics. SoFIE has been developed to meet the requirements of five unique investigations. It is currently being built and slated to begin operations on the ISS in July 2021. Given the general capabilities of the hardware insert, it is intended to be used as a facility for future researchers who can propose to NASA for related solid combustion studies.

Ferkul, Paul↗

The dynamics and optimal control of spinning spacecraft and movable telescoping appendages, part B: Effect of gravity-gradient torques on the dynamics of a spinning spacecraft with telescoping appendages

The effects of gravity gradient torques during boom deployment maneuvers of a spinning spacecraft are examined. Configurations where the booms extended only along the hub principal axes and where one or two booms are offset from the principal axes were considered. For the special case of symmetric deployment (principal axes booms) the stability boundaries are determined, and a stability chart is used to study the system behavior. Possible cases of instability during this type of maneuver are identified. In the second configuration an expression for gravity torque about the hub center of mass was developed. The nonlinear equations of motion are solved numerically, and the substantial influence of the gravity torque during asymmetric deployment maneuvers is indicated.

Bainum, P. M.↗

Simulation of a Controlled Airfoil with Jets

Numerical simulations of a two-dimensional airfoil, controlled by an applied moment in pitch and an airfoil controlled by jets, were investigated. These simulations couple the Reynolds-averaged Navier-Stokes equations and Euler's equations of rigid body motion, with an active control system. Controllers for both systems were designed to track altitude commands and were evaluated by simulating a closed-loop altitude step response using the coupled system. The airfoil controlled by a pitching moment used an optimal state feedback controller. A closed-loop simulation, of the airfoil with an applied moment, showed that the trajectories compared very well with quasi-steady aerodynamic theory, providing a measure of validation. The airfoil with jets used a controller designed by robust control methods. A linear plant model for this system was identified using open-loop data generated by the nonlinear coupled system. A closed-loop simulation of the airfoil with jets, showed good tracking of an altitude command. This simulation also showed oscillations in the control input as a result of dynamics not accounted for in the control design. This research work demonstrates how computational fluid dynamics, coupled with rigid body dynamics, and a control law can be used to prototype control systems in problematic nonlinear flight regimes.

Allan, Brian G.↗

Constraining Numerical Geodynamo Modeling with Surface Observations

Numerical dynamo solutions have traditionally been generated entirely by a set of self-consistent differential equations that govern the spatial-temporal variation of the magnetic field, velocity field and other fields related to dynamo processes. In particular, those solutions are obtained with parameters very different from those appropriate for the Earth s core. Geophysical application of the numerical results therefore depends on correct understanding of the differences (errors) between the model outputs and the true states (truth) in the outer core. Part of the truth can be observed at the surface in the form of poloidal magnetic field. To understand these differences, or errors, we generate new initial model state (analysis) by assimilating sequentially the model outputs with the surface geomagnetic observations using an optimal interpolation scheme. The time evolution of the core state is then controlled by our MoSST core dynamics model. The final outputs (forecasts) are then compared with the surface observations as a means to test the success of the assimilation. We use the surface geomagnetic data back to year 1900 for our studies, with 5-year forecast and 20-year analysis periods. We intend to use the result; to understand time variation of the errors with the assimilation sequences, and the impact of the assimilation on other unobservable quantities, such as the toroidal field and the fluid velocity in the core.

Kuang, Weijia↗