Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “State Variable”

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 199 records · Page 11

The compressed state Kalman filter for nonlinear state estimation: Application to large-scale reservoir monitoring

Reservoir monitoring aims to provide snapshots of reservoir conditions and their uncertainties to assist operation management and risk analysis. These snapshots may contain millions of state variables, e.g., pressures and saturations, which can be estimated by assimilating data in real time using the Kalman filter (KF). However, the KF has a computational cost that scales quadratically with the number of unknowns, m, due to the cost of computing and storing the covariance and Jacobian matrices, along with their prod-ucts. The compressed state Kalman filter (CSKF) adapts the KF for solving large-scale monitoring problems. The CSKF uses N preselected orthogonal bases to compute an accurate rank-N approximation of the covari-ance that is close to the optimal spectral approximation given by SVD. The CSKF has a computational cost that scales linearly in m and uses an efficient matrix-free approach that propagates uncertainties using N + 1 forward model evaluations, where N << m. Here in this paper we present a generalized CSKF algorithm for nonlin-ear state estimation problems such as CO 2 monitoring. For simultaneous estimation of multiple types of state variables, the algorithm allows selecting bases that represent the variability of each state type. Through synthetic numerical experiments of CO 2 monitoring, we show that the CSKF can reproduce the Kal-man gain accurately even for large compression ratios (m/N). For a given computational cost, the CSKF uses a robust and flexible compression scheme that gives more reliable uncertainty estimates than the ensemble Kalman filter, which may display loss of ensemble variability leading to suboptimal uncertainty estimates.

97 MATHEMATICS AND COMPUTING↗

Application of Consider Covariance to the Extended Kalman Filter

The extended Kalman filter (EKF) is the basis for many applications of filtering theory to real-time problems where estimates of the state of a dynamical system are to be computed based upon some set of observations. The form of the EKF may vary somewhat from one application to another, but the fundamental principles are typically unchanged among these various applications. As is the case in many filtering applications, models of the dynamical system (differential equations describing the state variables) and models of the relationship between the observations and the state variables are created. These models typically employ a set of constants whose values are established my means of theory or experimental procedure. Since the estimates of the state are formed assuming that the models are perfect, any modeling errors will affect the accuracy of the computed estimates. Note that the modeling errors may be errors of commission (errors in terms included in the model) or omission (errors in terms excluded from the model). Consequently, it becomes imperative when evaluating the performance of real-time filters to evaluate the effect of modeling errors on the estimates of the state.

Lundberg, John B.↗

Analytic error analysis of cross section interpolation methods in nodal diffusion codes - II: Numerical results

This paper is the second part of a two-part paper that documents the numerical results for the partial derivatives model presented in part I. In this paper, we derive the error bounds for the analytical point-wise error expression and verify our bounds with numerical experiments. The point-wise error expressions make available, and bound, the sources that contribute to the total error of the interpolated cross section in terms of the Lagrange interpolation errors and the model form error. MPACT is used to generate two-group homogenized cross sections for Westinghouse's AP1000 Region 4 lattice to evaluate the accuracy of the bounds. Error bounds calculated over a grid are compared to numerical data for uni-variate and multi-variate interpolation. The point-wise error bounds of a typical case matrix - two branches in each state variable - are displayed for bi-variate interpolation in the state variables: moderator density, fuel temperature, and boron concentration. The error bounds are shown to be highly accurate compared to numerical results, and in accordance with the underlying physics. We then discuss and show how the sources of error contribute to the total error, and consider the improvement of each error source. Finally, we mention future work such as propagating our cross section error bounds through a reactivity calculation. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Synthesis of hover autopilots for rotary-wing VTOL aircraft

The practical situation is considered where imperfect information on only a few rotor and fuselage state variables is available. Filters are designed to estimate all the state variables from noisy measurements of fuselage pitch/roll angles and from noisy measurements of both fuselage and rotor pitch/roll angles. The mean square response of the vehicle to a very gusty, random wind is computed using various filter/controllers and is found to be quite satisfactory although, of course, not so good as when one has perfect information (idealized case). The second part of the report considers precision hover over a point on the ground. A vehicle model without rotor dynamics is used and feedback signals in position and integral of position error are added. The mean square response of the vehicle to a very gusty, random wind is computed, assuming perfect information feedback, and is found to be excellent. The integral error feedback gives zero position error for a steady wind, and smaller position error for a random wind.

Hall, W. E.↗

Network synthesis

A discussion, with numerous examples, on the application of state variable methods to network analysis and synthesis is reported. The state variable point of view is useful in the design of control circuits for regulators because, unlike frequency domain methods, it is applicable to linear and nonlinear problems. The reported are intended as an introduction to this theory.

Brockett, R. W.↗

A continuum damage model of fatigue-induced damage in laminated composites

A model is presented which predicts the stress-strain behavior of continuous fiber reinforced laminated composites in the presence of microstructural damage. The model is based on the concept of continuum damage mechanics and uses internal state variables to characterize the various damage modes. The associated internal state variable growth laws are mathematical models of the loading history induced development of microstructural damage. The model is demonstrated by using it to predict the response of damaged AS-4/3502 graphite/epoxy laminate panels.

Harris, Charles E.↗

Survey of Methods to Predict Controller Workload for Real-Time Monitoring of Airspace Safety

The Real Time Safety Monitoring (RTSM) approach allows for assessment and prediction of the safety margin in the National Airspace System (NAS) to help preempt incidents and accidents, rather than having to reactively mitigate them. In RTSM, the NAS is modeled using state variables, and safety metrics are defined in terms of these state variables. The safety metrics have been classified as weather-related, airspace-related, and human-related. Many of the formulated human-related safety metrics need an estimate of the controller workload for their computation. However, this computation is not trivial. Hence, in this report, we perform a literature survey to identify the different factors that enable the computation of controller workload and categorize these factors. Next, we describe studies undertaken to determine a minimal set of factors that pro- vide a correct assessment of controller workload. Lastly, we survey approaches for evaluating how well the selected factors correlate with the controllers' subjective assessment of their workload. Based on this survey, we present factors beneficial to computing and predicting controller workload in real time, and discuss the status of data sources necessary for these computations.

National Airspace System↗

Science in Exploration: From the Moon to Mars and Back Home to Earth

NASA is embarking on a grand journey of exploration that naturally integrates the past successes of the Apollo missions to the Moon, as well as robotic science missions to Mars, to Planet Earth, and to the broader Universe. The US Vision for Space Exporation (VSE) boldly lays out a plan for human and robotic reconnaissance of the accessible Universe, starting with the surface of the Moon, and later embracing the surface of Mars. Sustained human and robotic access to the Moon and Mars will enable a new era of scientific investigation of our planetary neighbors, tied to driving scientific questions that pertain to the evolution and destiny of our home planet, but which also can be related to the search habitable worlds across the nearby Universe. The Apollo missions provide a vital legacy for what can be learned from the Moon, and NASA is now poised to recapture the lunar frontier starting with the flight of the Lunar Reconnaissance Orbiter (LRO) in late 2008. LRO will provide a new scientific context from which joint human and robotic exploration will ensue, guided by objectives some of which are focused on the grandest scientific challenges imaginable : Where did we come from? Are we alone? and Where are we going? The Moon will serve as an essential stepping stone for sustained human access and exploration of deep space and as a training ground while robotic missions with ever increasing complexity probe the wonders of Mars. As we speak, an armada of spacecraft are actively investigating the red planet both from orbit (NASA's Mars Reconnaissance Orbiter and Mars Odyssey Orbiter, plus ESA's Mars Express) and from the surface (NASA's twin Mars Exploration Rovers, and in 2008 NASA's Phoenix polar lander). The dramatically changing views of Mars as a potentially habitable world, with its own flavor of global climate change and unique climate records, provides a new vantage point from which to observe and question the workings of our own planet Earth. By 2010 NASA will have its first mobile analytical laboratory operating on the surface of Mars (Mars Science Laboratory) in search of potentially subtle expressions of past life or at least of life-hospitable environments. Meanwhile back here on Planet Earth, NASA will be continuing to implement an increasingly comprehensive program of robotic missions that address major issues associated with global climate variability, and the "state variables" that affect the quality of human life on our home planet. Ultimately, the fmits of NASA's emergent program of Exploration (VSE) will provide never-beforepossible opportunities for scientific leadership and advancement, culminating in a new state of awareness from which to better plan for the sustainability of life on Earth and for extending Earth life to the Moon and eventually to Mars. As NASA nears its 50th anniversary, the unimaginable and unexpected wealth of strategic knowledge its missions have generated about Earth, the Universe, and our local Solar System boggles the mind and serves as a legacy of knowledge for Educators to inspire future generations.

Garvin, James B.↗

A viscoplastic theory with thermodynamic considerations

A thermodynamic foundation using the concept of internal state variables is given for a general theory of viscoplasticity for initially isotropic materials. Three fundamental, internal, state variables are admitted; they are: a tensorial back stress for kinematic effects, and scalar drag and yield strengths for isotropic effects. All three are considered to evolve phenomenologically according to competitive processes between strain hardening, deformation-induced dynamic recovery, and thermally induced static recovery. Within this phenomenological framework, a thermodynamically admissible set of evolution equations is proposed. The theory allows each of the three internal variables to be composed as a sum of independently evolving constituents. The evolution of internal state can also include terms that vary linearly with the external variable rates, whose presence affects the energy dissipation properties of a material.

Freed, A. D.↗

A viscoplastic theory with thermodynamic considerations

A thermodynamic foundation using the concept of internal state variables is given for a general theory of viscoplasticity for initially isotropic materials. Three, fundamental, internal, state variables are admitted; they are: a tensorial back stress for kinematic effects, and scalar drag and yield strengths for isotropic effects. All three are considered to evolve phenomenologically according to competitive processes between strain hardening, deformation induced dynamic recovery, and thermally induced static recovery. Within this phenomenological framework, a thermodynamically admissible set of evolution equations is proposed. The theory allows each of the three internal variables to be composed as a sum of independently evolving constituents. The evolution of internal state can also include terms that vary linearly with the external variable rates, whose presence affects the energy dissipation properties of a material.

Freed, A. D.↗

Optimal management of electric hotel loads in mild hybrid heavy duty truck

Here, the problem of engine idling for heavy-duty trucks has been under study for decades with Auxiliary Power Units (APUs) and Truck Stop Electrification (TSE) as the most compelling solutions. With the electrification of trucks approaching feasibility in terms of cost-effective technology, hybridization offers another “degree of freedom” to tackle the problem. Here, this work aims at exploiting a battery pack of a 48 V mild-hybrid heavy-duty truck to store sufficient onboard energy for powering the auxiliary loads during the hoteling. This problem is not trivial, as the battery packs typically cannot recover the entire energy required through regeneration alone; hence an optimal energy management strategy needs to be employed to charge the battery through the engine during drive operation. This strategy optimizes powertrain performances among the four modes: (i) Engine of Coasting (EOC), (ii) Regeneration by braking, (iii) Regeneration by engine, and (iv) engine idling. This paper presents the development of a Dynamic Programming (DP) framework that employs a multi-objective cost function to minimize the fuel consumption and maximize the regeneration using the above-mentioned four modes. A typical heavy-duty truck drive cycle is used to represent the drive phase, with mandatory hoteling stops as per regulations. A comprehensive powertrain model is developed using validated components’ model. The DP employs two state variables: battery State-of-Charge (SOC) and engine mode, and three control inputs: (i) the engine ON–OFF state, (ii) clutch engagement state, and (iii) power request at the Electric Machine (EM) for calculating optimal SOC trajectory. The framework also tackles rapid engine ON–OFF scenarios to avoid the challenges associated with DP and the compromises in fuel cost with those approaches. Finally, the effectiveness of the proposed framework is tested for potential fuel savings on two different battery packs by performing the full cycle simulations. The results show 6.47% of fuel consumption reduction as compared to traditional APU-based heavy-duty truck.

33 ADVANCED PROPULSION SYSTEMS↗

A modified model parametrization algorithm for solving a special type of heat and mass transfer systems

A new method for solving nonlinear heat and mass transfer design tasks was considered. Systems using the Number of Transfer Units (NTU) method are a special type of mathematical model of heat and mass exchangers. It was observed, that the NTU models in a form of differential-algebraic equations (DAEs) cannot be directly solved with higher values of NTU. The requirements for consistent initial conditions, as well as numerical limitations of DAEs solvers, result, that the solution to the considered design problems that cannot be obtained by a classical direct shooting procedure. To overcome the presented difficulties, the αDAE model optimization algorithm was adjusted for solving NTU-based models. The new approach consists of 3 main steps: 1) task discretization by a multiple-shooting approach, 2) design an appropriate function $f_{NTU}$(α) to effectively influence the variability of the state variables described by dynamical relations, 3) the iterative numerical optimization algorithm for the new parametrized system. Moreover, computations can be performed by a chosen numerical optimization approach, which can be communicated with an available outer procedure for solving differential-algebraic equations. The presented algorithm was implemented and applied to solve the design task with the NTU model of a counter-flow exchanger. Here, the new approach was used to modify the system dynamics to influence the difficulty of the considered problem. Finally, the presented method enabled failure-free numerical computations for the higher values of the NTU parameter.

97 MATHEMATICS AND COMPUTING↗

Propagating Parameter Uncertainty in Power System Nonlinear Dynamic Simulations Using a Koopman Operator-Based Surrogate Model

In this work, we propose a Koopman operator-based surrogate model for propagating parameter uncertainties in power system nonlinear dynamic simulations. First, we augment a priori known state-space model by reformulating parameters deemed uncertain as pseudo-state variables. Then, we apply the Koopman operator theory to the resulting state-space model and obtain a linear dynamical system model. This transformation allows us to analyze the evolution of the system dynamics through its Koopman eigenfunctions, eigenvalues, and modes. Of particular importance for this letter, the obtained linear dynamical system is a surrogate that enables the evaluation of parameter uncertainties by simply perturbing the initial conditions of the Koopman eigenfunctions associated with the pseudo-state variables. Simulations carried out on the New England test system reveal the excellent performance of the proposed method in terms of accuracy and computational efficiency.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Systems theoretic model of stratospheric pollution on ozone

Discrete state space analysis is applied to the development of a model for predicting the reduction in the amount of stratospheric ozone due to stratospheric pollution. The modeling approach is to segment the stratosphere into a number of altitudinal bands and to associate state variables with the state of ozone reduction in each band due to NO(x) injections. The ozone reductions in each band are then combined into an estimate of total global ozone reduction. The model proposed predicts changes in stratospheric ozone based on altitude of pollutant injection and volume of pollutant injected at that altitude. The present approach is not a detailed analysis of stratospheric constituents and their interactions, but rather an analysis procedure for modeling the results of other researchers in the field.

Wilkins, R. D.↗

Consequences of Different Air-Sea Feedbacks on Ocean Using MITgcm and MERRA-2 Forcing: Implications for Coupled Data Assimilation Systems

Ocean surface flux estimates from atmospheric and oceanic reanalyses contain errors that compensate for inaccuracies in the respective atmosphere and ocean models used to generate these reanalyses. A conundrum for climate studies is the discrepancy between surface fluxes that minimize model-data differences for an atmosphere-only model vs surface fluxes that minimize model-data differences for an ocean model. As a first step towards a consistent coupled ocean-atmosphere data-assimilation (DA) system, we compare surface net heat flux from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to net heat flux from a state-of-the-art ocean state estimate, the Estimating the Circulation and Climate of the Ocean Version 4 (ECCO-v4). The possible impacts of the MERRA-2 and ECCO-v4 air-sea net heat flux difference in a coupled DA system were assessed using a set of experiments designed to imitate different “flavors” of a coupled DA system in an ocean-only setup. This was done by forcing the ECCO-v4 underlying ocean model - the Massachusetts Institute of Technology general circulation model (MITgcm) - with different sets of MERRA-2 fields and utilizing different forcing methods. By doing so we were able to turn off different air-sea feedbacks which, in a coupled DA setup, are partially muted by the constraining observations. The set of experiments, therefore, represents a range of active feedbacks in different “flavors” of coupled data-assimilation systems. For the period 1992–2011, MERRA-2 net heat flux has a global mean difference of -4.9 Wm(exp -2) relative to ECCO-v4. When MERRA-2 surface fields are used to force MITgcm, imbalances in the energy and the hydrological cycles of MERRA-2, which are directly related to the fact that MERRA-2 was created without an interactive ocean, propagate to the ocean. The experiment in which MITgcm is forced with MERRA-2 fluxes (MERRA-2-flux experiment) results in a 2.5°C global mean Sea Surface Temperature (SST) cooling, a 1m reduction in global mean sea level, and other drastic changes in the large scale ocean circulation relative to those resulting when the MITgcm is forced with the optimized ECCO-v4 net heat flux (the ECCO-v4 experiment itself). When MITgcm is forced with MERRA-2 state variables (MERRA-2-state experiment), the SST is somewhat restored to the observed SST, but the errors are shifted to the water cycle, resulting in a global mean sea level increase of 2.7 m. To further explore the pros and cons of these two approaches, we introduce a new intermediate forcing method in which the ocean is forced with turbulent fluxes but has a long wave feedback. This method, unlike MERRA-2 state, preserves the MERRA-2 water and salinity cycles, and it reduces the SST error compared to the MERRA-2-flux experiment, but the SST is not as good as that in the MERRA-2-state experiment. Our results have implications for ocean-model forcing recipes and clearly reveal the undesirable consequences of limiting the feedbacks in either these types of experiments or in coupled DA.

Ehud Strobach↗

Differential Flatness and Cooperative Tracking in the Lorenz System

In this paper the control of the Lorenz system for both stabilization and tracking problems is studied via feedback linearization and differential flatness. By using the Rayleigh number as the control, only variable physically tunable, a barrier in the controllability of the system is incidentally imposed. This is reflected in the appearance of a singularity in the state transformation. Composite controllers that overcome this difficulty are designed and evaluated. The transition through the manifold defined by such a singularity is achieved by inducing a chaotic response within a boundary layer that contains it. Outside this region, a conventional feedback nonlinear control is applied. In this fashion, the authority of the control is enlarged to the whole. state space and the need for high control efforts is mitigated. In addition, the differential parametrization of the problem is used to track nonlinear functions of one state variable (single tracking) as well as several state variables (cooperative tracking). Control tasks that lead to integrable and non-integrable differential equations for the nominal flat output in steady-state are considered. In particular, a novel numerical strategy to deal with the non-integrable case is proposed. Numerical results validate very well the control design.

Crespo, Luis G.↗

Quantification of Variation in Combustion Instability Amplitude in a Multi-Nozzle Can Combustor

In this work, we quantify the level of variation in unstable combustion oscillation amplitudes and identify the source of this variation in a multi-nozzle can combustor. At conditions where pollutant emissions are reduced, lean-premixed combustors can undergo thermoacoustic instability when pressure and heat release rate oscillations couple. A commonly used method for suppressing instability is fuel staging, which is a method where fuel is unevenly distributed between nozzles in a multi-nozzle combustor. Our work follows others who have characterized the effect of fuel staging on combustion instability and the mechanisms by which it works during both steady-state and transient operation. One of the outcomes from our previous work was that certain instability operating points display a high level of pressure oscillation amplitude variation. Instead of oscillating at a constant limit-cycle amplitude, the pressure oscillation amplitude varies significantly in time. In this work, we use the concept of permutation entropy to quantify the level of variation in the pressure fluctuation amplitude. We correlate the level of variation with a number of state variables over a range of operating conditions, 291 test cases in all. These state variables include mixture equivalence ratio; transient timescale, amplitude, and direction; hardware temperatures; gas temperatures; and thermoacoustic damping and growth rates. Significant pressure oscillation amplitude variation occurs when the thermoacoustic damping rate is low. The damping and growth rates can be low for a number of reasons, but they are highly correlated with the metal temperature of the centerbody, where the flames are anchored; lower temperatures result in lower damping rates during stable operation and lower growth rates during unstable operation. These results show the importance of the thermal boundary condition on the time-dependent behavior of the thermoacoustic instability.

Westfall, Seth↗

Decoupled and linear quadratic regulator control of a large, flexible space antenna with an observer in the control loop

An analysis is performed to compare decoupled and linear quadratic regulator (LQR) procedures for the control of a large, flexible space antenna. Control objectives involve: (1) commanding changes in the rigid-body modes, (2) nulling initial disturbances in the rigid-body modes, or (3) nulling initial disturbances in the first three flexible modes. Control is achieved with two three-axis control-moment gyros located on the antenna column. Results are presented to illustrate various effects on control requirements for the two procedures. These effects include errors in the initial estimates of state variables, variations in the type, number, and location of sensors, and deletions of state-variable estimates for certain flexible modes after control activation. The advantages of incorporating a time lag in the control feedback are also illustrated. In addition, the effects of inoperative-control situations are analyzed with regard to control requirements and resultant modal responses. Comparisons are included which show the effects of perfect state feedback with no residual modes (ideal case). Time-history responses are presented to illustrate the various effects on the control procedures.

Hamer, H. A.↗