Search NASA⌕ Search

SEARCH · Search NASA

Results for “robust model predictive 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 109 records · Page 6

Aircraft Measurements of BrO, IO, Glyoxal, NO2, H2O, O2-O2 and Aerosol Extinction Profiles in the Tropics: Comparison with Aircraft-/Ship-Based in Situ and Lidar Measurements

Tropospheric chemistry of halogens and organic carbon over tropical oceans modifies ozone and atmospheric aerosols, yet atmospheric models remain largely untested for lack of vertically resolved measurements of bromine monoxide (BrO), iodine monoxide (IO) and small oxygenated hydrocarbons like glyoxal (CHOCHO) in the tropical troposphere. BrO, IO, glyoxal, nitrogen dioxide (NO2), water vapor (H2O) and O2-O2 collision complexes (O4/ were measured by the University of Colorado Airborne Multi-AXis Differential Optical Absorption Spectroscopy (CU AMAXDOAS) instrument, aerosol extinction by high spectral resolution lidar (HSRL), in situ aerosol size distributions by an ultra high sensitivity aerosol spectrometer (UHSAS) and in situ H2O by vertical-cavity surface-emitting laser (VCSEL) hygrometer. Data are presented from two research flights (RF12, RF17) aboard the National Science Foundation/ National Center for Atmospheric Research Gulfstream V aircraft over the tropical Eastern Pacific Ocean (tEPO) as part of the "Tropical Ocean tRoposphere Exchange of Reactive halogens and Oxygenated hydrocarbons" (TORERO) project (January/February 2012). We assess the accuracy of O4 slant column density (SCD) measurements in the presence and absence of aerosols. Our O4-inferred aerosol extinction profiles at 477 nm agree within 6% with HSRL in the boundary layer and closely resemble the renormalized profile shape of Mie calculations constrained by UHSAS at low (sub-Rayleigh) aerosol extinction in the free troposphere. CU AMAX-DOAS provides a flexible choice of geometry, which we exploit to minimize the SCD in the reference spectrum (SCDREF, maximize signal-to-noise ratio) and to test the robustness of BrO, IO and glyoxal differential SCDs. The RF12 case study was conducted in pristine marine and free tropospheric air. The RF17 case study was conducted above the NOAA RV Ka'imimoana (TORERO cruise, KA-12-01) and provides independent validation data from ship-based in situ cavity-enhanced DOAS and MAX-DOAS. Inside the marine boundary layer (MBL) no BrO was detected (smaller than 0.5 pptv), and 0.2-0.55 pptv IO and 32-36 pptv glyoxal were observed. The near-surface concentrations agree within 30% (IO) and 10% (glyoxal) between ship and air-craft. The BrO concentration strongly increased with altitude to 3.0 pptv at 14.5 km (RF12, 9.1 to 8.6 deg N; 101.2 to 97.4 deg W). At 14.5 km, 5-10 pptv NO2 agree with model predictions and demonstrate good control over separating tropospheric from stratospheric absorbers (NO2 and BrO). Our profile retrievals have 12-20 degrees of freedom (DoF) and up to 500m vertical resolution. The tropospheric BrO vertical column density (VCD) was 1.5 x 10(exp 13) molec cm(exp -2) (RF12) and at least 0.5 x 10(exp 13) molec cm(exp -2) (RF17, 0- 10 km, lower limit). Tropospheric IO VCDs correspond to 2.1 x 10(exp 12) molec cm(exp -2) (RF12) and 2.5 x 10(exp 12) molec cm(exp -2) (RF17) and glyoxal VCDs of 2.6 x 10(exp 14) molec cm(exp -2) (RF12) and 2.7 x 10(exp 14) molec cm(exp -2) (RF17). Surprisingly, essentially all BrO as well as the dominant IO and glyoxal VCD fraction was located above 2 km (IO: 58 plus or minus 5 %, 0.1-0.2 pptv; glyoxal: 52 plus or minus 5 %, 3-20 pptv). To our knowledge there are no previous vertically resolved measurements of BrO and glyoxal from aircraft in the tropical free troposphere. The atmospheric implications are briefly discussed. Future studies are necessary to better understand the sources and impacts of free tropospheric halogens and oxygenated hydrocarbons on tropospheric ozone, aerosols, mercury oxidation and the oxidation capacity of the atmosphere.

Volkamer, R.↗

Real-Time UAV Trajectory Prediction for Safety Monitoring in Low-Altitude Airspace

The rising number of small unmanned aerial vehicles (UAVs) expected in the next decade will enable a new series of commercial, service, and military operations in low altitude airspace as well as above densely populated areas. These operations may include on-demand delivery, medical transportation services, law enforcement operations, traffic surveillance and many more. Such unprecedented scenarios create the need for robust, efficient ways to monitor the UAV state in time to guarantee safety and mitigate contingencies throughout the operations. This work proposes a generalized monitoring and prediction methodology that utilizes realtime measurements of an autonomous UAV following a series of way-points. Two different methods, based on sinusoidal acceleration profiles and high-order splines, are utilized to generate the predicted path. The monitoring approach includes dynamic trajectory re-planning in the event of unexpected detour or hovering of the UAV during flight. It can be further extended to different vehicle types, to quantify uncertainty affecting the state variables, e.g., aerodynamic and other environmental effects, and can also be implemented to prognosticate safety-critical metrics which depend on the estimated flight path and required thrust. The proposed framework is implemented on a simplified, scalable UAV modeling and control system traversing 3D trajectories. Results presented include examples of real-time predictions of the UAV trajectories during flight and a critical analysis of the proposed scenarios under uncertainty constraints.

UAV trajectory prognosis↗

Modeling and Real-Time Process Monitoring of Organometallic Chemical Vapor Deposition of III-V Phosphides and Nitrides at Low and High Pressure

The purpose of this paper is to review modeling and real-time monitoring by robust methods of reflectance spectroscopy of organometallic chemical vapor deposition (OMCVD) processes in extreme regimes of pressure. The merits of p-polarized reflectance spectroscopy under the conditions of chemical beam epitaxy (CBE) and of internal transmission spectroscopy and principal angle spectroscopy at high pressure are assessed. In order to extend OMCVD to materials that exhibit large thermal decomposition pressure at their optimum growth temperature we have designed and built a differentially-pressure-controlled (DCP) OMCVD reactor for use at pressures greater than or equal to 6 atm. We also describe a compact hard-shell (CHS) reactor for extending the pressure range to 100 atm. At such very high pressure the decomposition of source vapors occurs in the vapor phase, and is coupled to flow dynamics and transport. Rate constants for homogeneous gas phase reactions can be predicted based on a combination of first principles and semi-empirical calculations. The pressure dependence of unimolecular rate constants is described by RRKM theory, but requires variational and anharmonicity corrections not included in presently available calculations with the exception of ammonia decomposition. Commercial codes that include chemical reactions and transport exist, but do not adequately cover at present the kinetics of heteroepitaxial crystal growth.

Bachmann, K. J.↗

Adaptive Inverse Control for Rotorcraft Vibration Reduction

This thesis extends the Least Mean Square (LMS) algorithm to solve the mult!ple-input, multiple-output problem of alleviating N/Rev (revolutions per minute by number of blades) helicopter fuselage vibration by means of adaptive inverse control. A frequency domain locally linear model is used to represent the transfer matrix relating the higher harmonic pitch control inputs to the harmonic vibration outputs to be controlled. By using the inverse matrix as the controller gain matrix, an adaptive inverse regulator is formed to alleviate the N/Rev vibration. The stability and rate of convergence properties of the extended LMS algorithm are discussed. It is shown that the stability ranges for the elements of the stability gain matrix are directly related to the eigenvalues of the vibration signal information matrix for the learning phase, but not for the control phase. The overall conclusion is that the LMS adaptive inverse control method can form a robust vibration control system, but will require some tuning of the input sensor gains, the stability gain matrix, and the amount of control relaxation to be used. The learning curve of the controller during the learning phase is shown to be quantitatively close to that predicted by averaging the learning curves of the normal modes. For higher order transfer matrices, a rough estimate of the inverse is needed to start the algorithm efficiently. The simulation results indicate that the factor which most influences LMS adaptive inverse control is the product of the control relaxation and the the stability gain matrix. A small stability gain matrix makes the controller less sensitive to relaxation selection, and permits faster and more stable vibration reduction, than by choosing the stability gain matrix large and the control relaxation term small. It is shown that the best selections of the stability gain matrix elements and the amount of control relaxation is basically a compromise between slow, stable convergence and fast convergence with increased possibility of unstable identification. In the simulation studies, the LMS adaptive inverse control algorithm is shown to be capable of adapting the inverse (controller) matrix to track changes in the flight conditions. The algorithm converges quickly for moderate disturbances, while taking longer for larger disturbances. Perfect knowledge of the inverse matrix is not required for good control of the N/Rev vibration. However it is shown that measurement noise will prevent the LMS adaptive inverse control technique from controlling the vibration, unless the signal averaging method presented is incorporated into the algorithm.

Jacklin, Stephen A.↗

Magnetization of the Lunar Crust

Magnetic fields measured by the satellite Lunar Prospector show large scale features resulting from remanently magnetized crust. Vector data synthesized at satellite altitude from a spherical harmonic model of the lunar crustal field, and the radial component of the magnetometer data, have been used to produce spatially continuous global magnetization models for the lunar crust. The magnetization is expressed in terms of localized basis functions, with a magnetization solution selected having the smallest root-mean square magnetization for a given fit to the data, controlled by a damping parameter. Suites of magnetization models for layers with thicknesses between 10 and 50 km are able to reproduce much of the input data, with global misfits of less than 0.5 nT (within the uncertainties of the data), and some surface field estimates. The magnetization distributions show robust magnitudes for a range of model thicknesses and damping parameters, however the magnetization direction is unconstrained. These global models suggest that magnetized sources of the lunar crust can be represented by a 30 km thick magnetized layer. Average magnetization values in magnetized regions are 30-40 mA/m, similar to the measured magnetizations of the Apollo samples and significantly weaker than crustal magnetizations for Mars and the Earth. These are the first global magnetization models for the Moon, providing lower bounds on the magnitude of lunar crustal magnetization in the absence of multiple sample returns, and can be used to predict the crustal contribution to the lunar magnetic field at a particular location.

Lunar↗

Dynamic Smagorinsky model on anisotropic grids

Large Eddy Simulation (LES) of complex-geometry flows often involves highly anisotropic meshes. To examine the performance of the dynamic Smagorinsky model in a controlled fashion on such grids, simulations of forced isotropic turbulence are performed using highly anisotropic discretizations. The resulting model coefficients are compared with a theoretical prediction (Scotti et al., 1993). Two extreme cases are considered: pancake-like grids, for which two directions are poorly resolved compared to the third, and pencil-like grids, where one direction is poorly resolved when compared to the other two. For pancake-like grids the dynamic model yields the results expected from the theory (increasing coefficient with increasing aspect ratio), whereas for pencil-like grids the dynamic model does not agree with the theoretical prediction (with detrimental effects only on smallest resolved scales). A possible explanation of the departure is attempted, and it is shown that the problem may be circumvented by using an isotropic test-filter at larger scales. Overall, all models considered give good large-scale results, confirming the general robustness of the dynamic and eddy-viscosity models. But in all cases, the predictions were poor for scales smaller than that of the worst resolved direction.

Scotti, A.↗

Assessment of Low Cycle Fatigue Behavior of Powder Metallurgy Alloy U720

The fatigue lives of modem powder metallurgy disk alloys are influenced by variabilities in alloy microstructure and mechanical properties. These properties can vary as functions of variables the different steps of materials/component processing: powder atomization, consolidation, extrusion, forging, heat treating, and machining. It is important to understand the relationship between the statistical variations in life and these variables, as well as the change in life distribution due to changes in fatigue loading conditions. The objective of this study was to investigate these relationships in a nickel-base disk superalloy, U720, produced using powder metallurgy processing. Multiple strain-controlled fatigue tests were performed at 538 C (1000 F) at limited sets of test conditions. Analyses were performed to: (1) assess variations of microstructure, mechanical properties, and LCF failure initiation sites as functions of disk processing and loading conditions; and (2) compare mean and minimum fatigue life predictions using different approaches for modeling the data from assorted test conditions. Significant variations in life were observed as functions of the disk processing variables evaluated. However, the lives of all specimens could still be combined and modeled together. The failure initiation sites for tests performed at a strain ratio R(sub epsilon) = epsilon(sub min)/epsilon(sub max) of 0 were different from those in tests at a strain ratio of -1. An approach could still be applied to account for the differences in mean and maximum stresses and strains. This allowed the data in tests of various conditions to be combined for more robust statistical estimates of mean and minimum lives.

Gabb, Tomothy P.↗

MSFC-TVC | TB-03: Derivation of Thrust Vector Control (TVC) Actuator-Force / Gimbal-Torque Transformation Matrix

Thrust vector control (TVC) systems for rocket engine propulsion traditionally use a simple linear relationship to convert between actuator forces and torques about the engine gimbal’s center-of-rotation (COR). As shown in Equation (1), the torque about the gimbal COR is proportional to the applied actuator force and the TVC moment arm (MA)—the perpendicular distance between the TVC actuator’s line-of-action (LOA) and the engine gimbal’s COR. While this fundamental relationship remains valid and accurate in a two-dimensional (2D), one-degree-of-freedom (1-DOF) context—particularly in its non-linear formulation as described in the ER63 Technical Bulletin TB-02 (Derivation of Thrust Vector Control (TVC) Engine-Gimbal / Actuator Moment-Arm Geometry)—it becomes limited when extended to three-dimensional (3D), two-degree-of-freedom (2-DOF) analyses. In 3D, out-of-plane angular displacements can arise, causing the gimbaled engine plane-of-motion to become non-coplanar with the respective actuator plane-of-motion. Such misalignments occur due to the rod-end (RE) and/or tail-stock (TE) mounting geometry of individual TVC actuators. These geometric complexities lead to inaccuracies in calculating engine torques and corresponding actuator forces if only the traditional TVC MA relationship is employed. For higher gimbal angular displacements (e.g., greater than 5 degrees), these inaccuracies become more pronounced, necessitating a more robust mathematical framework. This document presents a comprehensive 3D (2-DOF) derivation of a transformation matrix that accurately converts between gimbaled engine torques and TVC actuator forces. By incorporating the necessary geometric and rotational considerations, this new approach corrects the limitations of the traditional TVC MA method. Subsequent sections compare the newly formulated approach to the traditional method, demonstrating its enhanced accuracy and reliability for a broad range of gimbaled engine conditions.

Thrust Vector Control↗

Universal limiter for transient interpolation modeling of the advective transport equations: The ULTIMATE conservative difference scheme

A fresh approach is taken to the embarrassingly difficult problem of adequately modeling simple pure advection. An explicit conservative control-volume formation makes use of a universal limiter for transient interpolation modeling of the advective transport equations. This ULTIMATE conservative difference scheme is applied to unsteady, one-dimensional scalar pure advection at constant velocity, using three critical test profiles: an isolated sine-squared wave, a discontinuous step, and a semi-ellipse. The goal, of course, is to devise a single robust scheme which achieves sharp monotonic resolution of the step without corrupting the other profiles. The semi-ellipse is particularly challenging because of its combination of sudden and gradual changes in gradient. The ULTIMATE strategy can be applied to explicit conservation schemes of any order of accuracy. Second-order schemes are unsatisfactory, showing steepening and clipping typical of currently popular so-called high resolution shock-capturing of TVD schemes. The ULTIMATE third-order upwind scheme is highly satisfactory for most flows of practical importance. Higher order methods give predictably better step resolution, although even-order schemes generate a (monotonic) waviness in the difficult semi-ellipse simulation. Little is to be gained above ULTIMATE fifth-order upwinding which gives results close to the ultimate for which one might hope.

Leonard, B. P.↗

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the Summer Program - No. 7

The Seventh Summer Program of the Center for Turbulence Research took place in the four-week period, July 5 to July 31, 1998. This was the largest CTR Summer Program to date, involving thirty-six participants from the U. S. and nine other countries. Thirty-one Stanford and NASA-Ames staff members facilitated and contributed to most of the Summer projects. A new feature, and perhaps a preview of the future programs, was that many of the projects were executed on non-NASA computers. These included supercomputers located in Europe as well as those operated by the Departments of Defense and Energy in the United States. In addition, several simulation programs developed by the visiting participants at their home institutions were used. Another new feature was the prevalence of lap-top personal computers which were used by several participants to carry out some of the work that in the past were performed on desk-top workstations. We expect these trends to continue as computing power is enhanced and as more researchers (many of whom CTR alumni) use numerical simulations to study turbulent flows. CTR's main role continues to be in providing a forum for the study of turbulence for engineering analysis and in facilitating intellectual exchange among the leading researchers in the field. Once again the combustion group was the largest. Turbulent combustion has enjoyed remarkable progress in using simulations to address increasingly complex and practically more relevant questions. The combustion group's studies included such challenging topics as fuel evaporation, soot chemistry, and thermonuclear reactions. The latter study was one of three projects related to the Department of Energy's ASCI Program (www.llnl.gov/asci); the other two (rocket propulsion and fire safety) were carried out in the turbulence modeling group. The flow control and acoustics group demonstrated a successful application of the so-called evolution algorithms which actually led to a previously unknown forcing strategy for jets yielding increased spreading rate. A very efficient algorithm for flow in complex geometries with moving boundaries based on the immersed boundary forcing technique was tested with very encouraging results. Also a new strategy for the destruction of aircraft trailing vortices was introduced and tested. The Reynolds Averaged Modeling (RANS) group demonstrated that the elliptic relaxation concept for RANS calculations is also applicable to transonic flows with shocks; however, prediction of laminar/turbulent transition remains an important pacing item. A large fraction of the LES effort was devoted to the development and testing of a new algorithmic procedure (as opposed to phenomenological model) for subgrid scale modeling based on regularized de-filtering of the flow variables. This appears to be a very promising approach, and a significant effort is currently underway to assess its robustness in high Reynolds number flows and in conjunction with numerical methods for complex flows. As part of the Summer Program two review tutorials were given on Turbulent structures in hydrocarbon pool fires (Sheldon Tieszen), and Turbulent combustion modeling: from RANS to LES via DNS (Luc Vervisch); and two seminars entitled Assessment of turbulence models for engineering applications (Paul Durbin) and Subgrid-scale modeling for non-premixed, turbulent reacting flows (James Riley) were presented. A number of colleagues from universities, government agencies, and industry attended the final presentations of the participants on July 31 and participated in the discussions. There are twenty-six papers in this volume grouped in five areas. Each group is preceded with an overview by its coordinator.

Source record↗

Mathematical Modeling of Loop Heat Pipes with Multiple Capillary Pumps and Multiple Condensers: Stead State Stimulations - Part 1

Loop Heat Pipes (LHPs) have proven themselves as reliable and robust heat transport devices for spacecraft thermal control systems. So far, the LHPs in earth-orbit satellites perform very well as expected. Conventional LHPs usually consist of a single capillary pump for heat acquisition and a single condenser for heat rejection. Multiple pump/multiple condenser LHPs have shown to function very well in ground testing. Nevertheless, the test results of a dual pump/condenser LHP also revealed that the dual LHP behaved in a complicated manner due to the interaction between the pumps and condensers. Thus it is redundant to say that more research is needed before they are ready for 0-g deployment. One research area that perhaps compels immediate attention is the analytical modeling of LHPs, particularly the transient phenomena. Modeling a single pump/single condenser LHP is difficult enough. Only a handful of computer codes are available for both steady state and transient simulations of conventional LHPs. No previous effort was made to develop an analytical model (or even a complete theory) to predict the operational behavior of the multiple pump/multiple condenser LHP systems. The current research project offered a basic theory of the multiple pump/multiple condenser LHP operation. From it, a computer code was developed to predict the LHP saturation temperature in accordance with the system operating and environmental conditions.

Hoang, Triem T.↗

System Identification for Integrated Aircraft Development and Flight Testing [l'Identification Des Systemes Pour le Developpement Integre des Aeronefs et les Essais en Vol]

Over the last decades flight vehicles such as aircraft and helicopters entering service and requiring increased operational effectiveness have with few exceptions experienced prolonged flight test development to achieve full certification. In many cases the original requirements had later to be reduced to enable release to service. The impact on the customer, and manufacturer has been considerable leading to increased costs and or reduced operational capabilities. These costly experiences are largely a result of the flight vehicle not behaving as modelled and designed. The evaluation of flight test data can be used as a tool for validating windtunnel results and mathematical models describing the flight dynamical behaviour. In this sense the uncertainty of important aerodynamic stability and control parameters can be reduced and the confidence of aircraft mathematical models improved. An additional important factor comes from the implementation of active control systems offering the promise of significantly increased flight vehicle performance and operational capability. This approach extends the traditional trade-offs between aerodynamics, structures and propulsion systems to include full- time, full-authority fly-by-wire/light systems. It is imperative that the aerodynamic stability and control parameters of such integrated flight and propulsion control systems have to turn out inflight as predicted, since inherent stability margins will be lower and the flight control system must correct these deficiencies to provide flight critical redundancy and safety. With the methodology of system identification from flight tests it is possible to sense the control inputs and the flight vehicle reactions Such as accelerations, rates and attitudes. The mathematical model, e.g. the model structure and parameters, has to be determined from the relationship of the measured control inputs and the system's responses. The aim of this symposium was to review the present state of the art of flight vehicle system and parameter identification techniques, and to provide a critical appraisal of current methods developed and applied to flight test data in a number of NATO nations. Particular emphasis was placed on practical aspects and lessons learned in order to generate information useful to the flight test community in industry and government agencies. The technical papers share invaluable experience and emphasize the advances of flight vehicle system identification over the last years to the point where confidence and robustness level is now reasonably high. The symposium covered overviews of identification methodologies, flight test techniques, recent aircraft and helicopter application programs, and a session of short papers covering up-to-the-minute flight test results. A final discussion included prepared comments from experts and concluded with key issues learned in the application of system identification and future research needs. The essential benefits to NATO nations can be condensed as follows: More accurate mathematical models for high bandwidth flight control systems, Improved assessment and evaluation of flying qualities, High fidelity mathematical models for flight vehicle development and mission training simulators, and generally, Reduced flight test time and costs.

Advisory Group for Aerospace Research and Developm↗

Role of Subsurface Physics in the Assimilation of Surface Soil Moisture Observations

Root zone soil moisture controls the land-atmosphere exchange of water and energy and exhibits memory that may be useful for climate prediction at monthly scales. Assimilation of satellite-based surface soil moisture observations into a land surface model is an effective way to estimate large-scale root zone soil moisture. The propagation of surface information into deeper soil layers depends on the model-specific representation of subsurface physics that is used in the assimilation system. In a suite of experiments we assimilate synthetic surface soil moisture observations into four different models (Catchment, Mosaic, Noah and CLM) using the Ensemble Kalman Filter. We demonstrate that identical twin experiments significantly overestimate the information that can be obtained from the assimilation of surface soil moisture observations. The second key result indicates that the potential of surface soil moisture assimilation to improve root zone information is higher when the surface to root zone coupling is stronger. Our experiments also suggest that (faced with unknown true subsurface physics) overestimating surface to root zone coupling in the assimilation system provides more robust skill improvements in the root zone compared with underestimating the coupling. When CLM is excluded from the analysis, the skill improvements from using models with different vertical coupling strengths are comparable for different subsurface truths. Finally, the skill improvements through assimilation were found to be sensitive to the regional climate and soil types.

Reichle, R. H.↗

Ground Vibration Test of the Aerostructure Test Wing 2

The Aerostructures Test Wing (ATW) was developed to test unique concepts for flutter prediction and control synthesis. A follow-on to the successful ATW, denoted ATW2, was fabricated as a test bed to validate a variety of instrumentation in flight and to collect data for development of advanced signal processing algorithms for flutter prediction and aviation safety. As a means to estimate flutter speed, a ground vibration test (GVT) was performed. The results of a GVT are typically utilized to update structural dynamics finite element (FE) models used for flutter analysis. In this study, two GVT methodologies were explored to determine which nodes provide the best sensor locations: (i) effective independence and (ii) kinetic energy sorting algorithms. For measurement, ten and twenty sensors were used for three and 10 target test modes. A total of six accelerometer configurations measured frequencies and mode shapes. This included locations used in the original ATW GVT. Moreover, an optical measurement system was used to acquire data without mass effects added by conventional sensors. A considerable frequency shift was observed in comparing the data from the accelerometers to the optical data. The optical data provided robust data for use of the ATW2 finite element model update.

Herrera, Claudia↗

Potential for Integrating Entry Guidance into the Multi-Disciplinary Entry Vehicle Optimization Environment

The state-of-the-art in vehicle design decouples flight feasible trajectory generation from the optimization process of an entry spacecraft shape. The disadvantage to this decoupled process is seen when a particular aeroshell does not meet in-flight requirements when integrated into Guidance, Navigation, and Control simulations. It is postulated that the integration of a guidance algorithm into the design process will provide a real-time, rapid trajectory generation technique to enhance the robustness of vehicle design solutions. The potential benefit of this integration is a reduction in design cycles (possible cost savings) and increased accuracy in the aerothermal environment (possible mass savings). This work examines two aspects: 1) the performance of a reference tracking guidance algorithm for five different geometries with the same reference trajectory and 2) the potential of mass savings from improved aerothermal predictions. An Apollo Derived Guidance (ADG) algorithm is used in this study. The baseline geometry and five test case geometries were flown using the same baseline trajectory. The guided trajectory results are compared to separate trajectories determined in a vehicle optimization study conducted for NASA's Mars Entry, Descent, and Landing System Analysis. This study revealed several aspects regarding the potential gains and required developments for integrating a guidance algorithm into the vehicle optimization environment. First, the generation of flight feasible trajectories is only as good as the robustness of the guidance algorithm. The set of dispersed geometries modelled aerodynamic dispersions that ranged from +/-1% to +/-17% and a single extreme case was modelled where the aerodynamics were approximately 80% less than the baseline geometry. The ADG, as expected, was able to guide the vehicle into the aeroshell separation box at the target location for dispersions up to 17%, but failed for the 80% dispersion cases. Finally, the results revealed that including flight feasible trajectories for a set of dispersed geometries has the potential to save mass up to 430 kg.

spacecraft guidance↗

Greedy Sampling and Incremental Surrogate Model-Based Tailoring of Aeroservoelastic Model Database for Flexible Aircraft

This paper presents a data analysis and modeling framework to tailor and develop linear parameter-varying (LPV) aeroservoelastic (ASE) model database for flexible aircrafts in broad 2D flight parameter space. The Kriging surrogate model is constructed using ASE models at a fraction of grid points within the original model database, and then the ASE model at any flight condition can be obtained simply through surrogate model interpolation. The greedy sampling algorithm is developed to select the next sample point that carries the worst relative error between the surrogate model prediction and the benchmark model in the frequency domain among all input-output channels. The process is iterated to incrementally improve surrogate model accuracy till a pre-determined tolerance or iteration budget is met. The methodology is applied to the ASE model database of a flexible aircraft currently being tested at NASA/AFRC for flutter suppression and gust load alleviation. Our studies indicate that the proposed method can reduce the number of models in the original database by 67%. Even so the ASE models obtained through Kriging interpolation match the model in the original database constructed directly from the physics-based tool with the worst relative error far below 1%. The interpolated ASE model exhibits continuously-varying gains along a set of prescribed flight conditions. More importantly, the selected grid points are distributed non-uniformly in the parameter space, a) capturing the distinctly different dynamic behavior and its dependence on flight parameters, and b) reiterating the need and utility for adaptive space sampling techniques for ASE model database compaction. The present framework is directly extendible to high-dimensional flight parameter space, and can be used to guide the ASE model development, model order reduction, robust control synthesis and novel vehicle design of flexible aircraft.

numerical analysi↗

The Langley Stability and Transition Analysis Code (LASTRAC) : LST, Linear and Nonlinear PSE for 2-D, Axisymmetric, and Infinite Swept Wing Boundary Layers

During the past two decades, our understanding of laminar-turbulent transition flow physics has advanced significantly owing to, in a large part, the NASA program support such as the National Aerospace Plane (NASP), High-speed Civil Transport (HSCT), and Advanced Subsonic Technology (AST). Experimental, theoretical, as well as computational efforts on various issues such as receptivity and linear and nonlinear evolution of instability waves take part in broadening our knowledge base for this intricate flow phenomenon. Despite all these advances, transition prediction remains a nontrivial task for engineers due to the lack of a widely available, robust, and efficient prediction tool. The design and development of the LASTRAC code is aimed at providing one such engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. LASTRAC was written from scratch based on the state-of-the-art numerical methods for stability analysis and modem software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory (LST) or linear parabolized stability equations (LPSE) method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. Coupled with the built-in receptivity model that is currently under development, the nonlinear PSE method offers a synergistic approach to predict transition onset for a given disturbance environment based on first principles. This paper describes the governing equations, numerical methods, code development, and case studies for the current release of LASTRAC. Practical applications of LASTRAC are demonstrated for linear stability calculations, N-factor transition correlation, non-linear breakdown simulations, and controls of stationary crossflow instability in supersonic swept wing boundary layers.

Chang, Chau-Lyan↗

The dynamics and control of the orbiting spacecraft control laboratory experiment (SCOLE) during station keeping

A mathematical model is developed to predict the dynamics of the proposed orbiting Spacecraft Control Laboratory Experiment during the station keeping phase. The Shuttle as well as the reflector are assumed to be rigid, the mast is flexible and is assumed to undergo elastic displacements very small as compared with its length. The equations of motion are derived using a Newton-Euler formulation. The model includes the effects of gravity, flexibility, and orbital dynamics. The control is assumed to be provided to the system through the Shuttle's three torquers, and through six actuators located by pairs at two points on the mast and at the mass center of the reflector. At each of the locations, an actuator acts parallel to the roll axis while the other one acts parallel to the pitch axis. It is seen that, in the presence of gravity-gradient torques in the system dynamics, the system assumes a new equilibrium position about which the equations must be linearized, primarily due to the offset in the mast attachment point to the reflector. The linear regulator theory is used to derive control laws for the linear model of the SCOLE including the first four flexible modes. Numerical results confirm the robustness of this control strategy for station keeping with maximum control efforts significantly below saturation levels.

Bainum, Peter M.↗