SEARCH · Engineering Papers

Results for “model parameter estimation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS2022

Exploration of an Adaptive Routine for Battery Modeling

The purpose of this document is to explore the use of adaptive routines in battery modeling. The adaptive routines consist of real-time state estimators combined with battery parameter model components that are adjusted in real-time as battery data becomes available. Several aspects are explored. It is shown that model parameter identification is possible for simple battery models using available input/output data measurements. The online system identification used is recursive least squares. Model identification may be combined with a state observer such as the extended Kalman filter or the unscented Kalman filter to form an adaptive model combined with state estimation. However, such a combination is found to be problematic due to uncertainty, observability and stability issues. This paper is organized as follows. Section 1 introduces adaptive routines and possible roles they play in battery modeling. In Section 2 real-time parameter identification is described with results based on battery data. Section 3 reviews various state estimators and results using a simple battery model. In Section 4 parameter identification and state estimation are combined to form an adaptive routine. Finally, in Section 5 conclusions are drawn and future work is suggested.

Adaptive↗

NASA NTRS1983

Current status of system identification methodology

Large space structures, system identification, model formulation, experimental design, model order and structure determination, parameter estimation, reduced order modeling, and closed loops are considered.

Larimore, W. E.↗

NASA NTRS2013

Modeling Non-homologous End Joining

Non-homologous end joining (NHEJ) is the dominant DNA double strand break (DSB) repair pathway and involves several NHEJ proteins such as Ku, DNA-PKcs, XRCC4, Ligase IV and so on. Once DSBs are generated, Ku is first recruited to the DNA end, followed by other NHEJ proteins for DNA end processing and ligation. Because of the direct ligation of break ends without the need for a homologous template, NHEJ turns out to be an error-prone but efficient repair pathway. Some mechanisms have been proposed of how the efficiency of NHEJ repair is affected. The type of DNA damage is an important factor of NHEJ repair. For instance, the length of DNA fragment may determine the recruitment efficiency of NHEJ protein such as Ku [1], or the complexity of the DNA breaks [2] is accounted for the choice of NHEJ proteins and subpathway of NHEJ repair. On the other hand, the chromatin structure also plays a role of the accessibility of NHEJ protein to the DNA damage site. In this talk, some mathematical models of NHEJ, that consist of series of biochemical reactions complying with the laws of chemical reaction (e.g. mass action, etc.), will be introduced. By mathematical and numerical analysis and parameter estimation, the models are able to capture the qualitative biological features and show good agreement with experimental data. As conclusions, from the viewpoint of modeling, how the NHEJ proteins are recruited will be first discussed for connection between the classical sequential model [4] and recently proposed two-phase model [5]. Then how the NHEJ repair pathway is affected, by the length of DNA fragment [6], the complexity of DNA damage [7] and the chromatin structure [8], will be addressed

Li, Yongfeng↗

NASA NTRS2008

Identifying Bearing Rotodynamic Coefficients Using an Extended Kalman Filter

An Extended Kalman Filter is developed to estimate the linearized direct and indirect stiffness and damping force coefficients for bearings in rotor dynamic applications from noisy measurements of the shaft displacement in response to imbalance and impact excitation. The bearing properties are modeled as stochastic random variables using a Gauss-Markov model. Noise terms are introduced into the system model to account for all of the estimation error, including modeling errors and uncertainties and the propagation of measurement errors into the parameter estimates. The system model contains two user-defined parameters that can be tuned to improve the filter's performance; these parameters correspond to the covariance of the system and measurement noise variables. The filter is also strongly influenced by the initial values of the states and the error covariance matrix. The filter is demonstrated using numerically simulated data for a rotor bearing system with two identical bearings, which reduces the number of unknown linear dynamic coefficients to eight. The filter estimates for the direct damping coefficients and all four stiffness coefficients correlated well with actual values, whereas the estimates for the cross-coupled damping coefficients were the least accurate.

Miller, Brad A.↗

NASA NTRS2008

Identifying Bearing Rotordynamic Coefficients using an Extended Kalman Filter

An Extended Kalman Filter is developed to estimate the linearized direct and indirect stiffness and damping force coefficients for bearings in rotor-dynamic applications from noisy measurements of the shaft displacement in response to imbalance and impact excitation. The bearing properties are modeled as stochastic random variables using a Gauss-Markov model. Noise terms are introduced into the system model to account for all of the estimation error, including modeling errors and uncertainties and the propagation of measurement errors into the parameter estimates. The system model contains two user-defined parameters that can be tuned to improve the filter s performance; these parameters correspond to the covariance of the system and measurement noise variables. The filter is also strongly influenced by the initial values of the states and the error covariance matrix. The filter is demonstrated using numerically simulated data for a rotor-bearing system with two identical bearings, which reduces the number of unknown linear dynamic coefficients to eight. The filter estimates for the direct damping coefficients and all four stiffness coefficients correlated well with actual values, whereas the estimates for the cross-coupled damping coefficients were the least accurate.

Miller, Bard A.↗

NASA NTRS2008

Identifying Bearing Rotordynamic Coefficients using an Extended Kalman Filter

Miller, Brad A.↗

NASA NTRS1994

Adaptive model reduction for continuous systems via recursive rational interpolation

A method for adaptive identification of reduced-order models for continuous stable SISO and MIMO plants is presented. The method recursively finds a model whose transfer function (matrix) matches that of the plant on a set of frequencies chosen by the designer. The algorithm utilizes the Moving Discrete Fourier Transform (MDFT) to continuously monitor the frequency-domain profile of the system input and output signals. The MDFT is an efficient method of monitoring discrete points in the frequency domain of an evolving function of time. The model parameters are estimated from MDFT data using standard recursive parameter estimation techniques. The algorithm has been shown in simulations to be quite robust to additive noise in the inputs and outputs. A significant advantage of the method is that it enables a type of on-line model validation. This is accomplished by simultaneously identifying a number of models and comparing each with the plant in the frequency domain. Simulations of the method applied to an 8th-order SISO plant and a 10-state 2-input 2-output plant are presented. An example of on-line model validation applied to the SISO plant is also presented.

Lilly, John H.↗

NASA NTRS2000

A statistical model of the human core-temperature circadian rhythm

We formulate a statistical model of the human core-temperature circadian rhythm in which the circadian signal is modeled as a van der Pol oscillator, the thermoregulatory response is represented as a first-order autoregressive process, and the evoked effect of activity is modeled with a function specific for each circadian protocol. The new model directly links differential equation-based simulation models and harmonic regression analysis methods and permits statistical analysis of both static and dynamical properties of the circadian pacemaker from experimental data. We estimate the model parameters by using numerically efficient maximum likelihood algorithms and analyze human core-temperature data from forced desynchrony, free-run, and constant-routine protocols. By representing explicitly the dynamical effects of ambient light input to the human circadian pacemaker, the new model can estimate with high precision the correct intrinsic period of this oscillator ( approximately 24 h) from both free-run and forced desynchrony studies. Although the van der Pol model approximates well the dynamical features of the circadian pacemaker, the optimal dynamical model of the human biological clock may have a harmonic structure different from that of the van der Pol oscillator.

NASA Discipline Regulatory Physiology↗

NASA NTRS1996

Modeling of Aircraft Unsteady Aerodynamic Characteristics/Part 3 - Parameters Estimated from Flight Data: Parameters Estimated from Flight Data - Part 3

A nonlinear least squares algorithm for aircraft parameter estimation from flight data was developed. The postulated model for the analysis represented longitudinal, short period motion of an aircraft. The corresponding aerodynamic model equations included indicial functions (unsteady terms) and conventional stability and control derivatives. The indicial functions were modeled as simple exponential functions. The estimation procedure was applied in five examples. Four of the examples used simulated and flight data from small amplitude maneuvers to the F-18 HARV and X-31A aircraft. In the fifth example a rapid, large amplitude maneuver of the X-31 drop model was analyzed. From data analysis of small amplitude maneuvers ft was found that the model with conventional stability and control derivatives was adequate. Also, parameter estimation from a rapid, large amplitude maneuver did not reveal any noticeable presence of unsteady aerodynamics.

Klein, Vladislav↗

NASA NTRS1997

Multiwavelength Modeling of Nova Atmospheres

We have analyzed the early optically thick ultraviolet spectra of Nova OS And 1986 using a grid of spherically symmetric, non-LTE, line-blanketed, expanding model atmospheres and synthetic spectra with the following set of parameters: 5,000 less than or equal to T(sub model) less than or equal to 60,000K, solar abundances, (rho)(alpha) r(sup -3), v(sub max) = 2000 km/s, L = 6 x 10(exp 4) solar luminosity, and a statistical or microturbulent velocity of 50 km/s. We used the synthetic spectra to estimate the model parameters corresponding to the observed IUE spectra. The fits to the observations were then iteratively improved by changing the parameters of the model atmospheres, in particular T(sub model) and the abundances, to arrive at the best fits to the optically thick pseudo-continuum and the features found in the IUE spectra. The IUE spectra show two different optically thick subphases. The earliest spectra, taken a few days after maximum optical light, show a pseudo-continuum created by overlapping absorption lines. The later observations, taken approximately 3 weeks after maximum light, show the simultaneous presence of allowed, semi-forbidden, and forbidden lines in the observed spectra. Analysis of these phases indicate that OS And 86 had solar metallicities except for Mg which showed evidence of being underabundant by as much as a factor of 10. We determine a distance of 5.1 kpc to OS And 86 and derive a peak bolometric luminosity of approximately 5 x 10(exp 4) solar luminosity. The computed nova parameters provide insights into the physics of the early outburst and explain the spectra seen by IUE. Lastly, we find evidence in the later observations for large non-LTE effects of Fe II which, when included, lead to much better agreement with the observations.

Hauschildt, P. H.↗

NASA NTRS2007

Astrophysics to z approx. 10 with Gravitational Waves

The most useful characterization of a gravitational wave detector's performance is the accuracy with which astrophysical parameters of potential gravitational wave sources can be estimated. One of the most important source types for the Laser Interferometer Space Antenna (LISA) is inspiraling binaries of black holes. LISA can measure mass and spin to better than 1% for a wide range of masses, even out to high redshifts. The most difficult parameter to estimate accurately is almost always luminosity distance. Nonetheless, LISA can measure luminosity distance of intermediate-mass black hole binary systems (total mass approx.10(exp 4) solar mass) out to z approx.10 with distance accuracies approaching 25% in many cases. With this performance, LISA will be able to follow the merger history of black holes from the earliest mergers of proto-galaxies to the present. LISA's performance as a function of mass from 1 to 10(exp 7) solar mass and of redshift out to z approx. 30 will be described. The re-formulation of LISA's science requirements based on an instrument sensitivity model and parameter estimation will be described.

Stebbins, Robin↗

NASA NTRS1977

Value of Landsat in urban water resources planning

The reported investigation had the objective to evaluate the utility of satellite multispectral remote sensing in urban water resources planning. The results are presented of a study which was conducted to determine the economic impact of Landsat data. The use of Landsat data to estimate hydrologic model parameters employed in urban water resources planning is discussed. A decision regarding an employment of the Landsat data has to consider the tradeoff between data accuracy and cost. Bayesian decision theory is used in this connection. It is concluded that computer-aided interpretation of Landsat data is a highly cost-effective method of estimating the percentage of impervious area.

Jackson, T. J.↗

NASA NTRS2008

Estimation of Nutation Time Constant Model Parameters for On-Axis Spinning Spacecraft

Calculating an accurate nutation time constant for a spinning spacecraft is an important step for ensuring mission success. Spacecraft nutation is caused by energy dissipation about the spin axis. Propellant slosh in the spacecraft fuel tanks is the primary source for this dissipation and can be simulated using a forced motion spin table. Mechanical analogs, such as pendulums and rotors, are typically used to simulate propellant slosh. A strong desire exists for an automated method to determine these analog parameters. The method presented accomplishes this task by using a MATLAB Simulink/SimMechanics based simulation that utilizes the Parameter Estimation Tool.

Schlee, Keith↗

NASA NTRS1987

Probabilistic risk analysis of flying the space shuttle with and without fuel turbine discharge temperature redline protection

An exact mathematical model is presented of the Space Shuttle Main Engine computer voting logic in the presence of High Pressure Fuel Turbine (HPFT) overtemperature events and fuel turbine temperature sensor failures. The model provides estimates of the probability of erroneous engine shutdown and the probability of not detecting a HPFT overtemperature event. Because it is believed that the likelihood of sensor failures and overtemperature events in the HPFT greatly overshadows those in the High Pressure Oxygen Turbine (HPOT), this modeling effort focused on the HPFT. However, because the redline protection logic is the same for both turbines, estimation of the model parameters using relevant HPOT data would provide estimates of the risk of erroneous engine shutdown and an undetected overtemperature event in the HPOT. Because of the complexity of the model, it was necessary to program the solution which thus makes it feasible to accommodate a changing data base. This is considered to be of great interest because of the subjective nature in determining the relevancy of certain failures and the fact that the data base is constantly changing as a result of the frequent engine tests.

Howell, Leonard↗

NASA NTRS1997

Adaptation of a Fast Optimal Interpolation Algorithm to the Mapping of Oceangraphic Data

A fast, recently developed, multiscale optimal interpolation algorithm has been adapted to the mapping of hydrographic and other oceanographic data. This algorithm produces solution and error estimates which are consistent with those obtained from exact least squares methods, but at a small fraction of the computational cost. Problems whose solution would be completely impractical using exact least squares, that is, problems with tens or hundreds of thousands of measurements and estimation grid points, can easily be solved on a small workstation using the multiscale algorithm. In contrast to methods previously proposed for solving large least squares problems, our approach provides estimation error statistics while permitting long-range correlations, using all measurements, and permitting arbitrary measurement locations. The multiscale algorithm itself, published elsewhere, is not the focus of this paper. However, the algorithm requires statistical models having a very particular multiscale structure; it is the development of a class of multiscale statistical models, appropriate for oceanographic mapping problems, with which we concern ourselves in this paper. The approach is illustrated by mapping temperature in the northeastern Pacific. The number of hydrographic stations is kept deliberately small to show that multiscale and exact least squares results are comparable. A portion of the data were not used in the analysis; these data serve to test the multiscale estimates. A major advantage of the present approach is the ability to repeat the estimation procedure a large number of times for sensitivity studies, parameter estimation, and model testing. We have made available by anonymous Ftp a set of MATLAB-callable routines which implement the multiscale algorithm and the statistical models developed in this paper.

Menemenlis, Dimitris↗

NASA NTRS2012

Global Burned Area and Biomass Burning Emissions from Small Fires

In several biomes, including croplands, wooded savannas, and tropical forests, many small fires occur each year that are well below the detection limit of the current generation of global burned area products derived from moderate resolution surface reflectance imagery. Although these fires often generate thermal anomalies that can be detected by satellites, their contributions to burned area and carbon fluxes have not been systematically quantified across different regions and continents. Here we developed a preliminary method for combining 1-km thermal anomalies (active fires) and 500 m burned area observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) to estimate the influence of these fires. In our approach, we calculated the number of active fires inside and outside of 500 m burn scars derived from reflectance data. We estimated small fire burned area by computing the difference normalized burn ratio (dNBR) for these two sets of active fires and then combining these observations with other information. In a final step, we used the Global Fire Emissions Database version 3 (GFED3) biogeochemical model to estimate the impact of these fires on biomass burning emissions. We found that the spatial distribution of active fires and 500 m burned areas were in close agreement in ecosystems that experience large fires, including savannas across southern Africa and Australia and boreal forests in North America and Eurasia. In other areas, however, we observed many active fires outside of burned area perimeters. Fire radiative power was lower for this class of active fires. Small fires substantially increased burned area in several continental-scale regions, including Equatorial Asia (157%), Central America (143%), and Southeast Asia (90%) during 2001-2010. Globally, accounting for small fires increased total burned area by approximately by 35%, from 345 Mha/yr to 464 Mha/yr. A formal quantification of uncertainties was not possible, but sensitivity analyses of key model parameters caused estimates of global burned area increases from small fires to vary between 24% and 54%. Biomass burning carbon emissions increased by 35% at a global scale when small fires were included in GFED3, from 1.9 Pg C/yr to 2.5 Pg C/yr. The contribution of tropical forest fires to year-to-year variability in carbon fluxes increased because small fires amplified emissions from Central America, South America and Southeast Asia-regions where drought stress and burned area varied considerably from year to year in response to El Nino-Southern Oscillation and other climate modes.

Randerson, J. T.↗

NASA NTRS1985

On-orbit systems identification of flexible spacecraft

The process of on-orbit systems identification of flexible spacecraft is examined in terms of the difficulties that are expected because of the potentially unmanageable number of unknown model parameters due to the very high order system models involved. A Jordon block canonical form and global model parameters are used to reduce the number of unknown parameters to manageable numbers. A Bayesian approach is discussed which enables the merging of theoretical models with ground or on-orbit test results by using Unconditional Maximum Likelihood Parameter Estimation. A Modified Newton-Raphson technique is proposed for determining the model parameter estimates and an expected fit error criterion is recommended for the determination of the model structure and order reduction.

Taylor, L. W., Jr.↗

NASA NTRS1997

Thermal and Hydrologic Signatures of Soil Controls on Evaporation: A Combined Energy and Water Balance Approach with Implications for Remote Sensing of Evaporation

The overall goal of this research is to examine the feasibility of applying a newly developed diagnostic model of soil water evaporation to large land areas using remotely sensed input parameters. The model estimates the rate of soil evaporation during periods when it is limited by the net transport resulting from competing effects of capillary rise and drainage. The critical soil hydraulic properties are implicitly estimated via the intensity and duration of the first stage (energy limited) evaporation, removing a major obstacle in the remote estimation of evaporation over large areas. This duration, or "time to drying" (t(sub d)), is revealed through three signatures detectable in time series of remote sensing variables. The first is a break in soil albedo that occurs as a small vapor transmission zone develops near the surface. The second is a break in either surface to air temperature differences or in the diurnal surface temperature range, both of which indicate increased sensible heat flux (and/or storage) required to balance the decrease in latent heat flux. The third is a break in the temporal pattern of near surface soil moisture. Soil moisture tends to decrease rapidly during stage 1 drying (as water is removed from storage), and then become more or less constant during soil limited, or "stage 2" drying (as water is merely transmitted from deeper soil storage). The research tasks address: (1) improvements in model structure, including extensions to transpiration and aggregation over spatially variable soil and topographic landscape attributes; and (2) applications of the model using remotely sensed input parameters.

Salvucci, Guido D.↗