Results for “distributed state estimation” — Page 34

DOE OSTI2022

Femtosecond X-ray Spectroscopy Directly Quantifies Transient Excited-State Mixed Valency

Quantifying charge delocalization associated with short-lived photoexcited states of molecular complexes in solution remains experimentally challenging, requiring local element specific femtosecond experimental probes of time-evolving electron transfer. In this study, we quantify the evolving valence hole charge distribution in the photoexcited charge transfer state of a prototypical mixed valence bimetallic iron–ruthenium complex, [(CN) 5 Fe II CNRu III (NH 3 ) 5 ]-, in water by combining femtosecond X-ray spectroscopy measurements with time-dependent density functional theory calculations of the excited-state dynamics. We estimate the valence hole charge that accumulated at the Fe atom to be 0.6 ± 0.2, resulting from excited-state metal-to-metal charge transfer, on an ~ 60 fs time scale. Our combined experimental and computational approach provides a spectroscopic ruler for quantifying excited-state valency in solvated complexes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Real Time Applications Using Linear State Estimation Technology (RTA/LSE)

Electric Power Group, LLC (EPG) was awarded DOE OE0000849 to design, develop, and demonstrate three real-time applications for monitoring power system stability. The applications use phasor measurements and processing by a linear state estimator (LSE) to improve system visibility and accuracy. These applications leverage the large and expanding phasor measurement systems being deployed by utilities throughout the world, and specifically in the United States. Phasor measurements provide a much more detailed view of the power systems than traditional SCADA systems, so can provide better visibility of events and an early warning of developing problems. These three applications are real-time contingency analysis (RTCA), area angle monitoring (AAM), and a voltage stability index (VSI). The project was fully completed as proposed. All three applications were successfully implemented and demonstrated. There were challenges the caused project delays, however. While the RTCA is based on well-known concepts and readily available tools, its deployment using phasor measurements meant it needed to be adapted to sometimes sparse measurements and limited coverage areas. EPG met these challenges with extensive modeling and testing. The AAM is based on unproven research, so EPG had to deal with a lot of unknowns and develop additional methods for practical implementation. These activities took more time than anticipated, so required extra time for development. In addition, there was difficulty getting adequate measurement locations for successful deployment and utility installation delays. Consequently, the project required a time extension to complete. These three real-time applications, RTCA, AAM, and VSI, use only phasor measurement input. Since they do not require the EMS or other measurement systems, they can serve as a backup and security check for the EMS. Phasor measurements are validated and extended by an LSE. The LSE will always produce an output even when some inputs are impaired, so the RT applications will operate even when the EMS state estimator does not converge and the applications it supports cannot function. Consequently, these RT applications can provide valuable additional reliability for control center operation.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Daily Forecasting of Regional Epidemics of Coronavirus Disease with Bayesian Uncertainty Quantification, United States

To increase situational awareness and support evidence-based policymaking, we formulated a mathematical model for coronavirus disease transmission within a regional population. This compartmental model accounts for quarantine, self-isolation, social distancing, a nonexponentially distributed incubation period, asymptomatic persons, and mild and severe forms of symptomatic disease. We used Bayesian inference to calibrate region-specific models for consistency with daily reports of confirmed cases in the 15 most populous metropolitan statistical areas in the United States. We also quantified uncertainty in parameter estimates and forecasts. This online learning approach enables early identification of new trends despite considerable variability in case reporting.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2003

Formation of Supermassive Black Holes in Galactic Bulges: A Rotating Collapse Model Consistent with the M(sub BH-sigma) Relation

Motivated by the observed correlation between black hole masses M(sub BH) and the velocity dispersion sigma of host galaxies, we develop a theoretical model of black hole formation in galactic bulges (this paper generalizes an earlier ApJ Letter). The model assumes an initial state specified by a uniform rotation rate OMEGA and a density distribution of the form rho = a(sup 2)(sub eff)per2piGR(sup 2)(so that a(sub eff)is an effective transport speed). The black hole mass is determined when the centrifugal radius of the collapse flow exceeds the capture radius of the central black hole (for Schwarzschild geometry). This model reproduces the observed correlation between the estimated black hole masses and the velocity dispersions of galactic bulges, i.e., M(sub BH) approximately equal to 10(sup 8) solar mass(sigma per 200 kilometers per second)(sup 4) where sigma = the square root of 2a(sub eff). To obtain this normalization, the rotation rate OMEGA approximately equal to 2 x 10(exp -15) rad per second. The model also defines a bulge mass scale M(sub B). If we identify the scale M(sub B) with the bulge mass, the model determines the ratio mu(sub B) of black hole mass to the host mass: mu(sub B) approximately equal to 0.0024(sigma per 200 kilometer per second), again in reasonable agreement with observed values. In this scenario, supermassive black holes form quickly (in approximately 10(exp 5) yr) and are born rapidly rotating (with a per M approximately 0.9). This paper also shown how these results depend on the assumed initial conditions; the most important quantity is the initial distribution of specific angular momentum in the precollapse state.

Adams, Fred C.↗

DOE OSTI2023

HVAC control fine-grained occupancy pattern estimation

A thermal state of a plurality of zones of the building is updated according to a building thermal model and information received from temperature sensors of the building. Predicted occupant counts for an upcoming plurality of time slots for each of the plurality of zones are updated using actual occupancy counts for each of the plurality of zones. A misprediction type distribution for the upcoming plurality of time slots for each of the plurality of zones is updated, the misprediction type distribution indicating misprediction for true negatives, false positives, false negatives, and true positives. A total misprediction cost expectation is updated according to the predicted occupant counts and the misprediction type distribution. HVAC power for each of the plurality of zones is determined to optimize occupant thermal comfort weighted according to the predicted occupant counts while minimizing the total misprediction cost expectation. HVAC operation is controlled per the HVAC power.

Lin, Shan↗

NASA NTRS1980

Small impact craters in the lunar regolith - Their morphologies, relative ages, and rates of formation

The size frequency distributions, relative ages, and absolute age determination of lunar impact craters ranging in diameter from 1 to 100 m are discussed in relation to present-day and past impact rates. The size frequency distributions of craters produced by impacts on the lunar surface and of craters remaining after degradation with time and later impacts are examined, together with the steady state frequency curves for various crater morphologies. Relations between the relative ages of craters, their morphologies and their diameters are discussed, and the calibration of this kind of relationship with absolute crater ages by the location of individual craters with known absolute ages, diameters and frequency distribution statistics is considered. Estimates of present and past cratering rates based on crater size-frequency distributions are presented which suggest that the rate of impacting objects producing small lunar craters could have been constant over the last 100 million years.

Moore, H. J.↗

NASA NTRS1992

Accurate Inventories Of Irrigated Land

System for taking land-use inventories overcomes two problems in estimating extent of irrigated land: only small portion of large state surveyed in given year, and aerial photographs made on 1 day out of year do not provide adequate picture of areas growing more than one crop per year. Developed for state of California as guide to controlling, protecting, conserving, and distributing water within state. Adapted to any large area in which large amounts of irrigation water needed for agriculture. Combination of satellite images, aerial photography, and ground surveys yields data for computer analysis. Analyst also consults agricultural statistics, current farm reports, weather reports, and maps. These information sources aid in interpreting patterns, colors, textures, and shapes on Landsat-images.

Wall, S.↗

NASA NTRS1998

Acquisition and Early Losses of Rare Gases from the Deep Earth

Direct observations show that the deep Earth contains rare gases of solar composition distinct from those in the atmosphere. We examine the implications of mantle rare gas characteristics on acquisition of rare gases from the solar nebula and subsequent losses due to a large impact. Deep mantle rare gas concentrations and isotopic compositions can be obtained from a model of transport and distribution of mantle rare gases. This model assumes the lower mantle closed early, while the upper mantle is open to subduction from the atmosphere and mass transfer from the lower mantle. Constraints are derived that can be incorporated into models for terrestrial volatile acquisition: (1) Calculated lower-mantle Xe-isotopic ratios indicate that the fraction of radiogenic Xe produced by I-129 and Pu-244 during the first about 10(exp 8) yr was lost, a conclusion also drawn for atmospheric Xe. Thus, either the Earth was made from materials that had lost >99% of rare gases about (0.7-2) x 10(exp 8) yr after the solar system formed, or gases were then lost from the fully formed Earth. (2) Concentrations of 3He and 20Ne in the lower mantle were established after these losses. (3) Neon-isotopic data indicates that mantle Ne has solar composition. The model allows for solar Ar/Ne and Xe/Ne in the lower mantle if a dominant fraction of upper mantle Ar and Xe are subduction-derived. If Earth formed in the presence of the solar nebula, it could have been melted by accretional energy and the blanketing effect of a massive, nebula-derived atmosphere. Gases from this atmosphere would have been sequestered within the molten Earth by dissolution at the surface and downward mixing. It was found that too much Ne would be dissolved in the Earth unless the atmosphere began to escape when the Earth was only partially assembled. Here we consider conditions required to initially dissolve sufficient rare gases to account for the present lower mantle concentrations after subsequent losses at 10(exp 8) yr. It is assumed that equilibration of the atmosphere with a thoroughly molten mantle was rapid, so that initial abundances of gases retained in any mantle layer reflected surface conditions when the layer solidified. For subsequent gas loss of 99.5% and typical solubility coefficients, a total pressure of 100 atm was required for an atmosphere of solar composition. Calculations of the pressure at the base of a primordial atmosphere indicate that this value might be exceeded by an order of magnitude or more for an atmosphere supported by accretional energy. Surface temperatures of about 4000 K would have been produced, probably high enough to melt the deep mantle. Initial distributions of retained rare gases would then be determined by the history of surface pressure and temperature during mantle cooling and solidification, i.e., the coupled cooling of Earth and atmosphere. The Earth's thermal state was determined by its surface temperature and the efficiency of convection in the molten mantle, estimated to be sufficient to maintain an adiabatic gradient. Because the melting curve is steeper than the adiabat, solidification of the mantle proceeded outward from the interior. Incorporation of atmospheric gases in the mantle therefore occurred over a range in surface temperature of a few thousand degrees Kelvin. The thermal state of the atmosphere was controlled by total luminosity of the Earth (energy) released by accreting planetesimals and the cooling Earth), nebular temperature and pressure, and atmospheric opacity. The energy released by accretion declined with time as did nebular pressure. Analytical solutions for an idealized (constant opacity radiative atmosphere show that declining energy sources under constant nebular conditions result in slowly diminishing surface temperature but dramatically increasing surface pressure. For such an atmosphere with declining nebular pressure but constant total luminosity, surface pressure decreases gradually with decreasing temperaure. A decline in accretion luminosity might be compensated by energy released as the mantle cools for about 10(exp 5) year, after which luminosity must decline. The total complement of dissolved rare gases will depend on the particular evolutionary path determined by the declining accretional luminosity, the Earth thermal history, removal of the nebula, and opacity variations of the atmosphere. Models for these coupled evolutionary histories for Earth's acquisition of nebular-derived noble gases are in progress. The later losses required at about 10(exp 8) yr (depleting the interior concentrations of the sequestered solar gases by a factor of > 100) were presumably related to the major impact in which the Moon formed.

Porcelli, D.↗

NASA NTRSDate not supplied

Highlights of NASA’s Orbital Debris Program Office In Situ and Laboratory Measurements

NASA’s Orbital Debris Program Office (ODPO) maintains various returned spacecraft materials, capabilities, and facilities used for in situ and laboratory measurements that directly support orbital debris environmental models. In situ measurements include the analysis of exposed and returned hardware surfaces. These surfaces serve as passive sensors for the small-sized micrometeoroid and orbital debris (MMOD) flux below the sensitivity of ground-based radar and optical sensors. Various instruments and techniques are used to determine the size and depth of selected impact features, and – if feasible – the composition of the projectile material. Analysis of the impactor residues enables the differentiation of MM and OD for debris below 1 mm to support modeling the OD environment. In addition, projectiles identified as OD can be further differentiated in low-, medium-, and high- density impactors based on chemical analyses. In addition to in situ measurements, the ODPO has also worked in collaboration with the U.S. Space Force Space Systems Command (formerly the U.S. Air Force Space and Missile Systems Center), the Aerospace Corporation, and the University of Florida on a laboratory-based hypervelocity impact test, DebriSat, conducted at the Air Force Arnold Engineering Development Complex in 2014. The resulting data from this impact test series are being analyzed to assess the fragments’ sizes/masses, materials/densities, shapes, and other parameters of interest. The DebriSat project provides the data needed to update NASA’s breakup models and size estimation models using the simulated orbital breakup of a modern, low Earth orbit spacecraft. Ultimately, over 200,000 fragments from this impact test will be stored at NASA Johnson Space Center (JSC) and further analyzed by the ODPO. This project will also use machine learning techniques to infer physical parameters of fragments embedded in the soft-catch foam used in the impact experiment. Applied to X-ray imagery of the foam panels, these techniques promise to minimize human-in-the-loop processes for fragment extraction and physical characterization. A brief overview of this project and data collected will be presented. Lastly, the ODPO provides various capabilities hosted at NASA JSC for optical inspections and measurements using a variety of techniques and scientific instrumentation to support both in situ and laboratory measurements. The ODPO’s Optical Measurement Center (OMC) is an advanced facility for photometric and spectroscopic laboratory measurements of targets, including fragments from the DebriSat project. The OMC simulates telescopic observations by using space-like illumination conditions and source-target-sensor orientation techniques. Additionally, the OMC is uniquely equipped to acquire pseudo-bidirectional reflectance distribution data for broadband photometric measurements, thus removing aspect angle dependencies that can affect target size estimates using the optical size estimation model. Narrow-band surface material characterization using spectroscopic instrumentation gives insight into how the albedo parameter – also important in the optical size estimation model – may vary depending on the state of the material. The OMC also performs simulations of photometric measurements using optical ray-tracing software to model the OMC optical throughput. In addition to the OMC, the ODPO houses a start-of-the-art Fragment Analysis Facility that uses multiple microscopic inspection instruments to support in situ measurements and material characterization. An overview of both facilities will be highlighted in this paper.

Orbital Debris↗

NASA NTRSDate not supplied

Highlights of NASA’s Orbital Debris Program Office In Situ and Laboratory Measurements

NASA’s Orbital Debris Program Office (ODPO) maintains various returned spacecraft materials, capabilities, and facilities used for in situ and laboratory measurements that directly support orbital debris environmental models. In situ measurements include the analysis of exposed and returned hardware surfaces. These surfaces serve as passive sensors for the small-sized micrometeoroid and orbital debris (MMOD) flux below the sensitivity of ground-based radar and optical sensors. Various instruments and techniques are used to determine the size and depth of selected impact features, and – if feasible – the composition of the projectile material. Analysis of the impactor residues enables the differentiation of MM and OD for debris below 1 mm to support modeling the OD environment. In addition, projectiles identified as OD can be further differentiated in low-, medium-, and high- density impactors based on chemical analyses. In addition to in situ measurements, the ODPO has also worked in collaboration with the U.S. Space Force Space Systems Command (formerly the U.S. Air Force Space and Missile Systems Center), the Aerospace Corporation, and the University of Florida on a laboratory-based hypervelocity impact test, DebriSat, conducted at the Air Force Arnold Engineering Development Complex in 2014. The resulting data from this impact test series are being analyzed to assess the fragments’ sizes/masses, materials/densities, shapes, and other parameters of interest. The DebriSat project provides the data needed to update NASA’s breakup models and size estimation models using the simulated orbital breakup of a modern, low Earth orbit spacecraft. Ultimately, over 200,000 fragments from this impact test will be stored at NASA Johnson Space Center (JSC) and further analyzed by the ODPO. This project will also use machine learning techniques to infer physical parameters of fragments embedded in the soft-catch foam used in the impact experiment. Applied to X-ray imagery of the foam panels, these techniques promise to minimize human-in-the-loop processes for fragment extraction and physical characterization. A brief overview of this project and data collected will be presented. Lastly, the ODPO provides various capabilities hosted at NASA JSC for optical inspections and measurements using a variety of techniques and scientific instrumentation to support both in situ and laboratory measurements. The ODPO’s Optical Measurement Center (OMC) is an advanced facility for photometric and spectroscopic laboratory measurements of targets, including fragments from the DebriSat project. The OMC simulates telescopic observations by using space-like illumination conditions and source-target-sensor orientation techniques. Additionally, the OMC is uniquely equipped to acquire pseudo-bidirectional reflectance distribution data for broadband photometric measurements, thus removing aspect angle dependencies that can affect target size estimates using the optical size estimation model. Narrow-band surface material characterization using spectroscopic instrumentation gives insight into how the albedo parameter – also important in the optical size estimation model – may vary depending on the state of the material. The OMC also performs simulations of photometric measurements using optical ray-tracing software to model the OMC optical throughput. In addition to the OMC, the ODPO houses a start-of-the-art Fragment Analysis Facility that uses multiple microscopic inspection instruments to support in situ measurements and material characterization. An overview of both facilities will be highlighted in this paper.

Orbital Debris↗

NASA NTRSDate not supplied

Bond Line Thickness Estimation in Composite Structures Using Multiple Inspection Techniques

Imaging and other nondestructive evaluation techniques are commonly used for material characterization and defect recognition in safety critical aerospace applications, with data fusion providing the framework for uncertainty quantification in these contexts. Most commonly, forward physics-based modeling predicts the response conditioned on material properties and defect assumptions, and probabilistic methods are used to infer the hidden state of the subject of the inspection from a combination of prior information, likelihoods, and inspection data. In this paper Bayesian methods are used to estimate bond thickness in lap joints comprised of aluminum adherends using a combination of infrared thermography and ultrasound. The concept of the conflation of probability distributions is applied to combine the posterior distributions derived from thermography and ultrasound and the quality of the fused estimates are compared against the individual estimates against synthetic data that was created to mimic the inspection of a lap joint comprised of aluminum adherends.

thermal nondestructive evaluation↗

DOE OSTI2023

Large-Eddy Simulation of Flow Over Boeing Gaussian Bump Using Multiagent Reinforcement Learning Wall Model: Preprint

We develop a wall model for large-eddy simulation (LES) that takes into account various pressure-gradient effects using multi-agent reinforcement learning. The model is trained using low-Reynolds-number flow over periodic hills with agents distributed on the wall at various computational grid points. It utilizes a wall eddy-viscosity formulation as the boundary condition to apply the modeled wall shear stress. Each agent receives states based on local instantaneous flow quantities at an off-wall location, computes a reward based on the estimated wall-shear stress, and provides an action to update the wall eddy viscosity at each time step. The trained wall model is validated in wall-modeled LES of flow over periodic hills at higher Reynolds numbers, and the results show the effectiveness of the model on flow with pressure gradients. The analysis of the trained model indicates that the model is capable of distinguishing between the various pressure gradient regimes present in the flow. To further assess the robustness of the developed wall model, simulations of flow over the Boeing Gaussian bump are conducted at a Reynolds number of 2 x 10^6, based on the free-stream velocity and the bump width. The results of mean skin friction and pressure on the bump surface, as well as the velocity statistics of the flow field, are compared to those obtained from equilibrium wall model (EQWM) simulations and published experimental data sets. The developed wall model is found to successfully capture the acceleration and deceleration of the turbulent boundary layer on the bump surface, providing better predictions of skin friction near the bump peak and exhibiting comparable performance to the EQWM with respect to the wall pressure and velocity field. We also conclude that the subgrid-scale model is crucial to the accurate prediction of the flow field, in particular the prediction of separation.

boundary layer↗

DOE OSTI2025

In-Situ Synchrotron X-Ray Diffraction of Ultrasonic Microstructural Refinement During Solidification in a Commercial Al–Si–Mg Alloy

This study reports the first use of in-situ synchrotron X-ray diffraction (SXRD) to study the effects of ultrasonic melt processing (USMP) on phase and grain size evolution during solidification in a commercial Al–Si–Mg casting alloy. USMP is a technique that, when applied to aluminum as it solidifies, can be used to refine the local microstructure of large-scale castings. Analysis of the in-situ SXRD data to estimate the average grain size of primary α-Al grains during USMP demonstrates that USMP slows the growth rate of α-Al grains and reduces grain size by 36 pct. Furthermore, there is also evidence that USMP causes the primary α-Al grains to move relative to the X-ray beam; such motion increases the probability of primary α-Al grains colliding and fragmenting. This movement becomes constrained at the onset of the Al–Si binary eutectic, suggesting that USMP ceases to effectively refine the microstructure once the Al–Si binary eutectic begins to form. Complementary laboratory-scale X-ray diffraction (XRD) data were used to correlate the lattice parameters of the α-Al and Si (D-A4) phases with temperature to estimate cooling rate during solidification. Thus, this study can guide the design of novel castings with spatially distributed fine-grained regions produced using local ultrasonic processing.

Aluminum Alloys↗

DOE OSTI2026

Microscopic theory of angular momentum distributions across the full range of fission fragments

Modern nuclear theory provides qualitative insights into the fundamental mechanisms of nuclear fission and is increasingly capable of making reliable quantitative predictions. Most quantities of interest pertain to the primary fission fragments, whose subsequent decay is typically modeled using statistical reaction models. Consequently, a key objective of fission theory is to inform these models by predicting the initial conditions of the primary fragments. Here, in this work, we employ a framework that combines joint angular momentum and particle number projection with time-dependent configuration mixing to calculate the angular momentum distributions of primary fragments. Focusing on the benchmark cases of neutron-induced fission of 235 U and 239 Pu , we predict—for the first time—microscopic angular momentum distributions for all fragments observed in experiments. Our results reveal a pronounced sawtooth pattern in the average angular momentum as a function of fragment mass, consistent with recent measurements. Additionally, we observe substantial variations in angular momentum distributions along isobaric chains, indicating that commonly used empirical formulas lack sufficient accuracy. We also quantify a strong correlation between the angular momentum and the deformation of the fragments at scission, and a weak correlation in the magnitude of the angular momentum between fragment partners. The generated data will enable estimation of the impact of microscopic distributions on fission spectra, paving the way toward fission modeling based on microscopic inputs.

Physics - Nuclear physics and radiation physics↗

NASA NTRS2012

Combined Use of Airborne Lidar and DBInSAR Data to Estimate LAI in Temperate Mixed Forests

The objective of this study was to determine whether leaf area index (LAI) in temperate mixed forests is best estimated using multiple-return airborne laser scanning (lidar) data or dual-band, single-pass interferometric synthetic aperture radar data (from GeoSAR) alone, or both in combination. In situ measurements of LAI were made using the LiCor LAI-2000 Plant Canopy Analyzer on 61 plots (21 hardwood, 36 pine, 4 mixed pine hardwood; stand age ranging from 12-164 years; mean height ranging from 0.4 to 41.2 m) in the Appomattox-Buckingham State Forest, Virginia, USA. Lidar distributional metrics were calculated for all returns and for ten one meter deep crown density slices (a new metric), five above and five below the mode of the vegetation returns for each plot. GeoSAR metrics were calculated from the X-band backscatter coefficients (four looks) as well as both X- and P-band interferometric heights and magnitudes for each plot. Lidar metrics alone explained 69% of the variability in LAI, while GeoSAR metrics alone explained 52%. However, combining the lidar and GeoSAR metrics increased the R2 to 0.77 with a CV-RMSE of 0.42. This study indicates the clear potential for X-band backscatter and interferometric height (both now available from spaceborne sensors), when combined with small-footprint lidar data, to improve LAI estimation in temperate mixed forests.

leaf area index↗

NASA NTRS1994

Finite element approximation of an optimal control problem for the von Karman equations

This paper is concerned with optimal control problems for the von Karman equations with distributed controls. We first show that optimal solutions exist. We then show that Lagrange multipliers may be used to enforce the constraints and derive an optimality system from which optimal states and controls may be deduced. Finally we define finite element approximations of solutions for the optimality system and derive error estimates for the approximations.

Hou, L. Steven↗

NASA NTRS2015

Architecture, Voltage, and Components for a Turboelectric Distributed Propulsion Electric Grid (AVC-TeDP)

The purpose of this effort was to advance the selection, characterization, and modeling of a propulsion electric grid for a Turboelectric Distributed Propulsion (TeDP) system for transport aircraft. The TeDP aircraft would constitute a miniature electric grid with 50 MW or more of total power, two or more generators, redundant transmission lines, and multiple electric motors driving propulsion fans. The study proposed power system architectures, investigated electromechanical and solid state circuit breakers, estimated the impact of the system voltage on system mass, and recommended DC bus voltage range. The study assumed an all cryogenic power system. Detailed assumptions within the study include hybrid circuit breakers, a two cryogen system, and supercritical cyrogens. A dynamic model was developed to investigate control and parameter selection.

power↗

NASA NTRS2018

Fuzzy Modeling and Parallel Distributed Compensation for Aircraft Flight Control from Simulated Flight Data

A method is described that combines fuzzy system identification techniques with Parallel Distributed Compensation (PDC) to develop nonlinear control methods for aircraft using minimal a priori knowledge, as part of NASA’s Learn-to-Fly initiative. A fuzzy model was generated with simulated flight data, and consisted of a weighted average of multiple linear time invariant state-space cells having parameters estimated using the equation-error approach and a least-squares estimator. A compensator was designed for each subsystem using Linear Matrix Inequalities (LMI) to guarantee closed-loop stability and performance requirements. This approach is demonstrated using simulated flight data to automatically develop a fuzzy model and design control laws for a simplified longitudinal approximation of the F-16 nonlinear flight dynamics simulation. Results include a comparison of flight data with the estimated fuzzy models and simulations that illustrate the feasibility and utility of the combined fuzzy modeling and control approach.

Weinstein, Rose↗