Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamic 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.

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

At least 541 records · Page 30

A unified perspective on robot control - The energy Lyapunov function approach

A unified framework for the stability analysis of robot tracking control is presented. By using an energy-motivated Lyapunov function candidate, the closed-loop stability is shown for a large family of control laws sharing a common structure of proportional and derivative feedback and a model-based feedforward. The feedforward can be zero, partial or complete linearized dynamics, partial or complete nonlinear dynamics, or linearized or nonlinear dynamics with parameter adaptation. As result, the dichotomous approaches to the robot control problem based on the open-loop linearization and nonlinear Lyapunov analysis are both included in this treatment. Furthermore, quantitative estimates of the trade-offs between different schemes in terms of the tracking performance, steady state error, domain of convergence, realtime computation load and required a prior model information are derived.

Wen, John T.↗

Control System Modeling for the Thirty Meter Telescope Primary Mirror

The Thirty Meter Telescope primary mirror is composed of 492 segments that are controlled to high precision in the presence of wind and vibration disturbances, despite the interaction with structural dynamics. The higher bandwidth and larger number of segments compared with the Keck telescopes requires greater attention to modeling to ensure success. We focus here on the development and validation of a suite of quasi-static and dynamic modeling tools required to support the design process, including robustness verification, performance estimation, and requirements flowdown. Models are used to predict the dynamic response due to wind and vibration disturbances, estimate achievable bandwidth in the presence of control-structure-interaction (CSI) and uncertainty in the interaction matrix, and simulate and analyze control algorithms and strategies, e.g. for control of focus-mode, and sensor calibration. Representative results illustrate TMT performance scaling with parameters, but the emphasis is on the modeling framework itself.

uncetainty↗

The SMAP Level-4 ECO Product - Phase 1: Improving Vegetation Simulations Through Observation-Driven Parameter Estimation

Simulations of hydrological fields as well as water, energy and carbon fluxes from the land surface to the atmosphere are crucial for a wide range of applications, including agricultural advisories, forecasts of (short-term) atmospheric behavior and seasonal weather predictions including forecasts of extreme events, such as heatwaves or droughts. The NASA Soil Moisture Active Passive (SMAP) mission Level-4 (L4) Eco-Hydrology (ECO) project aims to improve modeled estimates of the terrestrial water, energy and carbon fluxes and states by developing a fully-coupled hydrology-vegetation data assimilation (DA) algorithm. The DA system is developed for the NASA Goddard Earth Observing System version 5 (GEOS-5) Catchment-CN land surface model, which combines land hydrology components of the GEOS-5 Catchment model with dynamic vegetation components of the Community Land Model version 4. Catchment-CN fully couples the terrestrial water, energy and carbon cycles, allowing feedbacks from the land hydrology to the biosphere and vice versa. For SMAP L4 ECO a calibration of the Catchment-CN vegetation parameterization against observations of the fraction of absorbed photosynthetically active radiation (FPAR) from the Moderate Resolution Imaging Spectroradiometer (MODIS) is implemented to improve the model's standalone skill. Next, the DA algorithm used to produce the SMAP L4 soil moisture product is adapted to Catchment-CN to assimilate SMAP brightness temperatures and inform the model's land hydrology component. The DA system is further extended to assimilate MODIS FPAR observations in order to constrain the model's dynamic vegetation component. In this presentation, we demonstrate that the Catchment-CN parameter calibration leads to more realistic vegetation simulations and reduces the root mean squared error between modeled and observed vegetation states across the model's various plant functional types. We also show that the assimilation of SMAP observations is able to improve the average correlation, bias and unbiased RMSE between the modeled surface and root zone soil moisture estimates, and ground observations from the SMAP core validation sites.

Kolassa, J.↗

ENSEMBLE KALMAN FILTERS FOR RELIABILITY ESTIMATION IN PERFUSION INFERENCE

We consider the solution of inverse problems in dynamic contrast–enhanced imaging by means of ensemble Kalman filters. Our quantity of interest is blood perfusion, i.e., blood flow rates in tissue. While existing approaches to compute blood perfusion parameters for given time series of radiological measurements mainly rely on deterministic, deconvolution–based methods, we aim at recovering probabilistic solution information for given noisy measurements. To this end, we model radiological image capturing as a sequential data assimilation process and solve it by an ensemble Kalman filter. Thereby, we recover deterministic results as an ensemble–based mean and are able to compute reliability information such as probabilities for the perfusion to be in a given range. Our target application is the inference of blood perfusion parameters in the human brain. A numerical study shows promising results for artificial measurements generated by a digital perfusion phantom.

Engineering↗

Further study of the atmospheric gravity waves over the eastern seaboard on 18 March 1969

Washington, D.C., microbarograph records for Mar. 18, 1969 reveal gravity-wave-associated pressure oscillations which appear to be directly related to upper tropospheric wave structure observed at the same time with a Wallops Island 10-cm wavelength radar. The consistency between the two sets of data provides observational support for the hypothesis that shear instability in the upper tropospheric flow is the mechanism responsible for the generation of such waves. The comparison presented suggests that microbarograph arrays might be useful adjuncts to future radar studies of upper tropospheric wave dynamics, supplying such wave parameters as phase velocity and wavelength in favorable cases. The radar data pertinent to this event show an apparent vertical wave phase variation, permitting a very approximate and somewhat uncertain estimate of the wave-associated vertical flux of horizontal momentum, which is found to be about 4 dyn per sq cm. Wave generation by shear instability in the upper tropospheric air flow and the resulting vertical momentum transport may be an important element of the global atmospheric momentum budget.

Hooke, W. H.↗

Control technology development

Static and dynamic control design approaches were developed for distributed parameter systems. A hardware flexible beam facility was constructed to demonstrate and verify the theoretical control concepts. Efforts were made in the area of model order estimation for control systems with uncertain or time varying parameters.

Rodriguez, G.↗

On-orbit Systems Identification of Flexible Spacecraft Technol., 1984, Pt. 2 p 465-481

Future spacecraft include configurations which are too flexible to be adequately tested prior to flight and which will require on-orbit systems identification to ensure safe operation of the flight control system. The structural dynamics model will evolve and its accuracy will improve in stages as ground tests of full-scale components and replica-scale models are performed. State Space Modeling and Conditional Maximum Likelihood Parameter Estimation methodology can provide the formal probability-based framework for the process of upgrading a model as additional test results are obtained. Although the number of unknown parameters can be reduced by the use of canonical forms for the stability matrix, the number of unknown model parameters quickly becomes unmanageable unless advantage is taken of the relationship of a much fewer number of global model parameters. Distributed parameter systems or partial differential equation models are one way to take advantage of such global parameters to reduce the number of unknown model parameters.

Taylor, L.↗

An engineering analysis of a closed cycle plant growth module

The SOLGEM model is a numerical engineering model which solves the flow and energy balance equations for the air flowing through a growing environment, assuming quasi-steady state conditions within the system. SOLGEM provides a dynamic simulation of the controlled environment system in that the temperature and flow conditions of the growing environment are estimated on an hourly basis in response to the weather data and the plant growth parameters. The flow energy balance considers the incident solar flux; incoming air temperature, humidity, and flow rate; heat exchange with the roof and floor; and heat and moisture exchange with the plants. A plant transpiration subroutine was developed based plant growth research facility, intended for the study of bioregenerative life support theories. The results of a performance analysis of the plant growth module are given. The estimated energy requirements of the module components and the total energy are given.

Stickford, G. H., Jr.↗

How Well Does the Latest Anthropomorphic Test Device Mimic Human Impact Responses?

One of the goals of the NASA Occupant Protection Group is to understand the human tolerance to dynamic loading. This knowledge has to come through indirect approaches such as existing human response databases, anthropometric test devices (ATD), animal testing, post‐mortem human subjects, and models. This study investigated the biofidelity of the National Highway Traffic Safety Administration's ATD named the THOR (test device for human occupant restraint). If THOR responds comparably to humans, then it could potentially be used as a human surrogate to help validate space vehicle requirements for occupant protection. The THOR responses to frontal and spinal impacts (ranging from 8 to 12 G with rise times of 40, 70, and 100 ms) were measured and compared to human volunteer responses (95 trials in frontal and 58 in spinal) previously collected by the U. S. Air Force on the same horizontal impact accelerator. The impact acceleration profiles tested are within the expected range of multi‐purpose crew vehicle (MPCV) landing dynamics. A correlation score was calculated for each THOR to human comparison using CORA (CORrelation and Analysis) software. A two‐parameter beta distribution model fit was obtained for each dependent variable using maximum likelihood estimation. For frontal impacts, the THOR head x‐acceleration peak response correlated with the human response at 8 and 10‐G 100 ms but not 10‐G 70 ms. The phase lagged the human response. Head z‐acceleration was not correlated. Chest x‐acceleration was in phase, had a higher peak response, and was well correlated with lighter subjects (Cora = 0.8 for 46 kg vs. Cora = 0.4 for 126 kg). Head x‐displacement had a leading phase. Several subjects responded with the same peak displacement but the mean of the group was lower. The shoulder x‐displacement was in phase but had higher peaks than the human response. For spinal impacts, the THOR head x‐acceleration was not well correlated. Head and chest z‐acceleration was in phase but had a higher peak response. Chest z‐acceleration was highly correlated with heavier subjects at lower G pulses (Cora = 0.86 for 125 kg at 8 G). The human response was variable in shoulder z‐displacement but the THOR was in phase and was comparable to the mean peak response. Head xand z‐displacement was in phase but had higher peaks. Seat pan forces were well correlated, were in phase, but had a larger peak response than most subjects. The THOR does not respond to frontal and spinal impacts exactly the same way that a human does. Some responses are well matched and others are not. Understanding the strengths and weaknesses of this ATD is an important first step in determining its usefulness in occupant protection at NASA

Newby, N.↗

How Well Does the Latest Anthropomorphic Test Device Mimic Human Impact Responses?

One of the goals of the NASA Occupant Protection Group is to understand the human tolerance to dynamic loading. This knowledge has to come through indirect approaches such as existing human response databases, anthropometric test devices (ATD), animal testing, post-­‐mortem human subjects, and models. This study investigated the biofidelity of the National Highway Traffic Safety Administration's ATD named the THOR (test device for human occupant restraint). If THOR responds comparably to humans, then it could potentially be used as a human surrogate to help validate space vehicle requirements for occupant protection. The THOR responses to frontal and spinal impacts (ranging from 8 to 12 G with rise times of 40, 70, and 100 ms) were measured and compared to human volunteer responses (95 trials in frontal and 58 in spinal) previously collected by the U. S. Air Force on the same horizontal impact accelerator. The impact acceleration profiles tested are within the expected range of multi-­‐purpose crew vehicle (MPCV) landing dynamics. A correlation score was calculated for each THOR to human comparison using CORA (CORrelation and Analysis) software. A two-­‐parameter beta distribution model fit was obtained for each dependent variable using maximum likelihood estimation. For frontal impacts, the THOR head x-­‐acceleration peak response correlated with the human response at 8 and 10-­‐G 100 ms but not 10-­‐G 70 ms. The phase lagged the human response. Head z-­‐acceleration was not correlated. Chest x-­‐acceleration was in phase, had a higher peak response, and was well correlated with lighter subjects (Cora = 0.8 for 46 kg vs. Cora = 0.4 for 126 kg). Head x-­‐displacement had a leading phase. Several subjects responded with the same peak displacement but the mean of the group was lower. The shoulder x-­‐displacement was in phase but had higher peaks than the human response. For spinal impacts, the THOR head x-­‐acceleration was not well correlated. Head and chest z-­‐acceleration was in phase but had a higher peak response. Chest z-­‐acceleration was highly correlated with heavier subjects at lower G pulses (Cora = 0.86 for 125 kg at 8 G). The human response was variable in shoulder z-­‐displacement but the THOR was in phase and was comparable to the mean peak response. Head x-­‐ and z-­‐displacement was in phase but had higher peaks. Seat pan forces were well correlated, were in phase, but had a larger peak response than most subjects. The THOR does not respond to frontal and spinal impacts exactly the same way that a human does. Some responses are well matched and others are not. Understanding the strengths and weaknesses of this ATD is an important first step in determining its usefulness in occupant protection at NASA

Newby, Nate↗

Sketch-to-Solution: A Case Study in RCS Aerodynamic Interaction

Thanks to recent advances in the fields of anisotropic grid adaptation, error estimation, and geometry modeling, a sketch-to-solution work flow is now possible for viscous computational fluid dynamic (CFD) simulations. With this workflow, a CFD application engineer provides geometry, boundary conditions, and flow parameters; and the sketch-to-solution process yields a CFD simulation through automatic, error-based, grid adaptation. To explore the benefits of this nascent capability, a conventional manual grid generation work flow is compared to this new automatic grid generation work flow for a given engineering question: What are the aerodynamic interactions caused by the reaction control system (RCS) on an entry vehicle? This case study indicates that while the automatic grid generation sketch-to- solution process is not yet mature, it is preferred over a manual grid generation work flow because it greatly reduces manual labor, eliminates many opportunities for human error, and provides grid sensitivity information.

Bill Kleb↗

Results of an Automated GPS Tracking System in Support of Topex/Poseidon and GPSMet

A fully automated near real-time GPS tracking system has been developed around JPL's GIPSY/OASIS II software. The system produces <25 cm (3D rms) GPS orbits and one-half nanosecond (15 cm) clock estimates. The process starts automatically when a favorable global distribution of ground data from the IGS network (International GPS Service for Geodynamics) becomes available. Ionospherically corrected phase and pseudorange data are optimally combined to remove satellite and ground receiver clock errors, including selective availability. After the GPS orbits are determined within the data arc, they are then propagated with empirically determined dynamic force models. Real-time <2 meter (3D rms) GPS orbits are always available. As a by-product of this process, other calibration estimates such as station clocks, troposphere estimates, and Earth orientation parameters are also produced. For daily arc fits, the process requires 7-8 hours of CPU time on an HP9000/735 workstation.

GPS↗

Determination of Celestial Body Principal Axes via Gravity Field Estimation

The determination of physical parameters of celestial bodies provides critical navigational and scientific information. Determining the mass, rotation state, and density distribution is an important task of the navigation team on an operational mission. One representation of the internal density distribution, spherical harmonics, may be leveraged to compute the orientation of the principal axis frame of the celestial body, which in turn informs the dynamics of the body’s motion. Using data from NASA’s Dawn mission to the massive asteroids (1) Ceres and (4) Vesta, a method is proposed that refines the knowledge of the principal axis frame by utilizing estimated spherical harmonic coefficients. Applying an iterative rotation scheme to the solved-for body-fixed frame leads to convergence on a frame that is very close to the actual dynamical principal axis frame.

Takahashi, Yu↗

Control of CuPt Ordering in GaInP Grown by Hydride Vapor Phase Epitaxy

We investigate the degree of ordering present in GaInP layers grown by dynamic hydride vapor phase epitaxy (D-HVPE) as a function of various growth conditions. We assess order parameter from both transmission electron and x-ray diffraction and compare these trends to the order parameter estimated by the photoluminescence-determined band gap in the GaInP layers. The degree of ordering increases with deposition temperature and the ratio of gas-phase group V to group III precursors in the reactor. We also investigate the effect of growth rate, demonstrating the capability of HVPE to suppress order through rapid growth rates. Such ordering dependences are consistent with mechanisms established for ordering observed in GaInP grown by organometallic vapor phase epitaxy (OMVPE), but with weaker ordering than GaInP grown by that method.

growth↗

Assimilation of a knowledge base and physical models to reduce errors in passive-microwave classifications of sea ice

An expert system framework has been developed to classify sea ice types using satellite passive microwave data, an operational classification algorithm, spatial and temporal information, ice types estimated from a dynamic-thermodynamic model, output from a neural network that detects the onset of melt, and knowledge about season and region. The rule base imposes boundary conditions upon the ice classification, modifies parameters in the ice algorithm, determines a `confidence' measure for the classified data, and under certain conditions, replaces the algorithm output with model output. Results demonstrate the potential power of such a system for minimizing overall error in the classification and for providing non-expert data users with a means of assessing the usefulness of the classification results for their applications.

Maslanik, J. A.↗

Computations of Axisymmetric Flows in Hypersonic Shock Tubes

A time-accurate two-dimensional fluid code is used to compute test times in shock tubes operated at supersonic speeds. Unlike previous studies, this investigation resolves the finer temporal details of the shock-tube flow by making use of modern supercomputers and state-of-the-art computational fluid dynamic solution techniques. The code, besides solving the time-dependent fluid equations, also accounts for the finite rate chemistry in the hypersonic environment. The flowfield solutions are used to estimate relevant shock-tube parameters for laminar flow, such as test times, and to predict density and velocity profiles. Boundary-layer parameters such as bar-delta(sub u), bar-delta(sup *), and bar-tau(sub w), and test time parameters such as bar-tau and particle time of flight t(sub f), are computed and compared with those evaluated by using Mirels' correlations. This article then discusses in detail the effects of flow nonuniformities on particle time-of-flight behind the normal shock and, consequently, on the interpretation of shock-tube data. This article concludes that for accurate interpretation of shock-tube data, a detailed analysis of flowfield parameters, using a computer code such as used in this study, must be performed.

Sharma, Surendra P.↗

Computational Simulations of the 10-MW TP3 Arc-Jet Facility Flow

This paper reports computational simulations and analysis in support of calibration tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using stagnation calorimeters and two different blunted wedge models with calibration plates at a wide range of conditions in the NASA Ames 10-megawatt TP3 (Test Position 3) facility. Data were obtained using four different conical nozzles with the same test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot and null-point heat flux probes were also performed at several arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested, and comparisons with the experimental measurements. These computational simulations provide estimates of the arc-jet test environment parameters that are not measured but are needed to evaluate the performance of thermal protection system materials, along with further valuable insights into the arc-jet testing environment. Simulation results are used to estimate centerline total enthalpy, surface shear, boundary layer thickness, and boundary layer edge Mach number and to verify that specific test requirements from the Orion program are met.

CFD↗

Derivative-based SINDy (DSINDy): Addressing the challenge of discovering governing equations from noisy data

Recent advances in the field of data-driven dynamics allow for the discovery of ODE systems using state measurements. One approach, known as Sparse Identification of Nonlinear Dynamics (SINDy), assumes the dynamics are sparse within a predetermined basis in the states and finds the expansion coefficients through linear regression with sparsity constraints. This approach requires an accurate estimation of the state time derivatives, which is not necessarily possible in the high-noise regime without additional constraints. We present an approach called Derivative-based SINDy (DSINDy) that combines two novel methods to improve ODE recovery at high-noise levels. First, we denoise the state variables by applying a projection operator that leverages the assumed basis for the system dynamics. Second, we use a second order cone program (SOCP) to find the derivative and governing equations simultaneously. We derive theoretical results for the projection-based denoising step, which allow us to estimate the values of hyperparameters used in the SOCP formulation. This underlying theory helps limit the number of required user-specified parameters. Finally, we present results demonstrating that our approach leads to improved system recovery for the Van der Pol oscillator, the Duffing oscillator, the Rössler attractor, and the Lorenz 96 model.

97 MATHEMATICS AND COMPUTING↗