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 379 records · Page 21

Data-driven discovery of dynamics from time-resolved coherent scattering

Coherent X-ray scattering (CXS) techniques are capable of interrogating dynamics of nano- to mesoscale materials systems at time scales spanning several orders of magnitude. However, obtaining accurate theoretical descriptions of complex dynamics is often limited by one or more factors—the ability to visualize dynamics in real space, computational cost of high-fidelity simulations, and effectiveness of approximate or phenomenological models. In this work, we develop a data-driven framework to uncover mechanistic models of dynamics directly from time-resolved CXS measurements without solving the phase reconstruction problem for the entire time series of diffraction patterns. Our approach uses neural differential equations to parameterize unknown real-space dynamics and implements a computational scattering forward model to relate real-space predictions to reciprocal-space observations. This method is shown to recover the dynamics of several computational model systems under various simulated conditions of measurement resolution and noise. Moreover, the trained model enables estimation of long-term dynamics well beyond the maximum observation time, which can be used to inform and refine experimental parameters in practice. Finally, we demonstrate an experimental proof-of-concept by applying our framework to recover the probe trajectory from a ptychographic scan. Our proposed framework bridges the wide existing gap between approximate models and complex data.

36 MATERIALS SCIENCE↗

Membrane interactions of the globular domain and the hypervariable region of KRAS4b define its unique diffusion behavior

The RAS proteins are GTP-dependent switches that regulate signaling pathways and are frequently mutated in cancer. RAS proteins concentrate in the plasma membrane via lipid-tethers and hypervariable region side-chain interactions in distinct nano-domains. However, little is known about RAS membrane dynamics and the details of RAS activation of downstream signaling. Here, we characterize RAS in live human and mouse cells using single-molecule-tracking methods and estimate RAS mobility parameters. KRAS4b exhibits confined mobility with three diffusive states distinct from the other RAS isoforms (KRAS4a, NRAS, and HRAS); and although most of the amino acid differences between RAS isoforms lie within the hypervariable region, the additional confinement of KRAS4b is largely determined by the protein’s globular domain. To understand the altered mobility of an oncogenic KRAS4b, we used complementary experimental and molecular dynamics simulation approaches to reveal a detailed mechanism.

59 BASIC BIOLOGICAL SCIENCES↗

Analysis of Nonlinear Shrinkage for the Bound Metal Deposition Manufacturing using Multi-scale Approach

We consider problem of nonlinear shrinkage of the metal part during bound metal deposition manufacturing on the ground and in zero-G. To analyze this problem we developed multi-scale physics-based approach that spans atomistic dynamics at the scale of nanoseconds and the full part shrinkage at the time scale of hours. Using this approach we estimated the key parameters of the problem including grain boundary width, coefficient of surface diffusion, initial redistribution of particles during debinding stage, micro-structure evolution from round particles to densely packed grains and corresponding change of the total and chemical free energy, and sintering stress. The introduced method was used to predict shrinkage at the level of two particles, filament cross-section, sub-model, and the whole green, brown, and metal parts. To further improve accuracy and reliability of the shrinkage predictions we propose concept of intelligent additive manufacturing of metal powders in space that combines the strengths of both physics-based and data-driven methods of analysis of AM.

bound metal deposition↗

Dynamics of filaments during the edge-localized mode crash on NSTX

Edge localized modes (ELMs) are routinely observed in H-mode plasma regimes of the National Spherical Torus Experiment (NSTX). Due to the explosive nature of the instability, only diagnostics with high temporal and spatial resolution could provide a detailed insight into the dynamics associated with the ELMs. Gas-puff imaging at NSTX provides 2D measurements of the magnetic field aligned fluctuations (e.g., ELM filaments) in the scrape-off layer and at the plasma edge with 2.5 μs temporal and 10 mm optical resolution. A novel analysis technique was developed to estimate the frame-by-frame velocities and the spatial parameters of the dominant structures associated with the ELMs. The analysis was applied to single ELM events to characterize the ELM crash dynamics and then extended to a database of 159 ELM events. Statistical analysis was performed in order to find the characterizing dynamics of the ELM crash. The results show that on average, an ELM crash consists of a filament with a circular cross section, which is propelled outward with a characterizing peak radial velocity of ~3.3 km/s. The radial velocity was found to be linearly dependent on the distance of the filament from the separatrix, which has never been seen before. The ELM filament is characterized by propagation in the ion-diamagnetic direction poloidally with a peak velocity of 11.4 km/s. The ELM crash lasts for approximately 100 μs until the radial propulsion settles back to the pre-ELM level. Furthermore, the experimental findings were compared with analytical theory. Two possible mechanisms were identified for explaining the observations: the curvature interchange model and the current–filament interaction model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Software for continuum modeling of controls-structures interactions

It is clear that computer software is needed to assist in the generation of the equations of motion for complex, flexible spacecraft. Daniel Poelaert of ESTEC has developed the software DISTEL with which he has modeled the structural dynamics for different satellites. He is interested in expanding the capabilities of DISTEL to include structural damping and control systems. Unfortunately, the software has not been released. The author has developed similar software, PDEMOD, which has been used to model the Spacecraft control Laboratory Experiment (SCOLE), the Solar Array Flight Experiment (SAFE), the Mini-MAST truss, and the LACE satellite. PDEMOD has been used also for optimal parameter estimation and integrated control-structures design. PDEMOD is also being extended to include structural damping and control systems which are imbedded into the same equations for the structural dynamics. This paper will address the formulation of the equations for the structural dynamics of spacecraft structures which are constructed of a 3-dimensional arrangement of rigid bodies and flexible beam elements. Control system dynamics are imbedded into the same equations so that model order reduction approximations are not necessary. The input data consists of the physical data of the elements and the topological information describing how the elements are connected. PDEMOD accomplishes the following: (1) automatically assembles the equations of motion for the entire structural model; (2) calculates the modal frequencies; (3) calculates the mode shapes; (4) generates perspective views of the mode shapes; and (5) forms selected transfer functions. The software PDEMOD continues to be developed to provide additional features to assist in analyzing and synthesizing control and structural systems for flexible spacecraft.

Taylor, Lawrence W., Jr.↗

Resonant Tidal Excitation of Internal Waves in the Earth's Fluid Core

It has long been speculated that there is a stably stratified layer below the core-mantle boundary, and two recent studies have improved the constraints on the parameters describing this stratification. Here we consider the dynamical implications of this layer using a simplified model. We first show that the stratification in this surface layer has sensitive control over the rate at which tidal energy is transferred to the core. We then show that when the stratification parameters from the recent studies are used in this model, a resonant configuration arrives whereby tidal forces perform elevated rates of work in exciting core flow. Specifically, the internal wave speed derived from the two independent studies (150 and 155 m/s) are in remarkable agreement with the speed (152 m/s) required for excitation of the primary normal mode of oscillation as calculated from full solutions of the Laplace Tidal Equations applied to a reduced-gravity idealized model representing the stratified layer. In evaluating this agreement it is noteworthy that the idealized model assumed may be regarded as the most reduced representation of the stratified dynamics of the layer, in that there are no non-essential dynamical terms in the governing equations assumed. While it is certainly possible that a more realistic treatment may require additional dynamical terms or coupling, it is also clear that this reduced representation includes no freedom for coercing the correlation described. This suggests that one must accept either (1) that tidal forces resonantly excite core flow and this is predicted by a simple model or (2) that either the independent estimates or the dynamical model does not accurately portray the core surface layer and there has simply been an unlikely coincidence between three estimates of a stratification parameter which would otherwise have a broad plausible range.

Geodynamo↗

Dynamics of filaments during the edge-localized mode crash on NSTX

Edge localized modes (ELMs) are routinely observed in H-mode plasma regimes of the National Spherical Torus Experiment (NSTX). Due to the explosive nature of the instability, only diagnostics with high temporal and spatial resolution could provide a detailed insight into the dynamics associated with the ELMs. Gas-puff imaging (GPI) at NSTX provides 2D measurements of the magnetic field aligned fluctuations (e.g. ELM filaments) in the scrape-off layer and the at the plasma edge with 2.5 us temporal and 10 mm optical resolution.A novel analysis technique was developed to estimate the frame-by-frame velocities and the spatial parameters of the dominant structures associated with the ELMs. The analysis was applied to single ELM events to characterize the ELM crash dynamics, and then extended to a database of 169 ELM events.Statistical analysis was performed in order to find the characterizing dynamics of the ELM crash. The results show that on average an ELM crash consists of a filament with a circular cross-section which is propelled outwards with a characterizing peak radial velocity of ~3.3 km/s. The radial velocity was found to be linearly dependent on the distance of the filament from the separatrix, which has never been seen before. The ELM filament is characterized by propagation in the ion-diamagnetic direction poloidally with a peak velocity of 11.4 km/s. The ELM crash lasts for approximately 100us until the radial propulsion settles back to the pre-ELM level. The experimental findings were compared with analytical theory. Two possible mechanisms were identified for explaining the observations: the curvature interchange model and the current-filament interaction model.

correlation↗

Estimation of the longitudinal and lateral-directional aerodynamic parameters from flight data for the NASA F/A-18 HARV

This progress report presents the results of an investigation focused on parameter identification for the NASA F/A-18 HARV. This aircraft was used in the high alpha research program at the NASA Dryden Flight Research Center. In this study the longitudinal and lateral-directional stability derivatives are estimated from flight data using the Maximum Likelihood method coupled with a Newton-Raphson minimization technique. The objective is to estimate an aerodynamic model describing the aircraft dynamics over a range of angle of attack from 5 deg to 60 deg. The mathematical model is built using the traditional static and dynamic derivative buildup. Flight data used in this analysis were from a variety of maneuvers. The longitudinal maneuvers included large amplitude multiple doublets, optimal inputs, frequency sweeps, and pilot pitch stick inputs. The lateral-directional maneuvers consisted of large amplitude multiple doublets, optimal inputs and pilot stick and rudder inputs. The parameter estimation code pEst, developed at NASA Dryden, was used in this investigation. Results of the estimation process from alpha = 5 deg to alpha = 60 deg are presented and discussed.

Napolitano, Marcello R.↗

Dealing with Measurement Uncertainties as Nuisance Parameters in Bayesian Model Calibration

In the presence of model discrepancy, the calibration of physics-based models for physical parameter inference is a challenging problem. Lack of identifiability between calibration parameters and model discrepancy requires additional identifiability constraints to be placed on the model discrepancy to obtain unique physical parameter estimates. If these assumptions are violated, the inference for the calibration parameters can be systematically biased. In many applications, such as in dynamic material property experiments, many of the calibration inputs refer to measurement uncertainties. Here, we develop a metric for identifying overfitting of these measurement uncertainties, propose a prior capable of reducing this overfitting, and show how this leads to a diagnostic tool for validation of physical parameter inference. The approach is demonstrated for a benchmark example and applied for a material property application to perform inference on the equation of state parameters of tantalum.

97 MATHEMATICS AND COMPUTING↗

Improved semiclassical model for real-time evaporation of matrix black holes

In this paper, we study real-time classical matrix mechanics of a simplified 2 × 2 matrix model inspired by the black hole evaporation problem. This is a step towards making a quantitative model of real-time evaporation of a black hole, which is realized as a bound state of D0-branes in string theory. The model we study is the reduction of Yang–Mills in 2 + 1 dimension to 0 + 1 dimensions, which has been corrected with an additional potential that can be interpreted as a zero-point energy for fermions. Our goal is to understand the lifetime of such a classical bound state object in the classical regime. To do so, we pay particular attention to when D-particles separate to check that the “off-diagonal modes” of the matrices become adiabatic and use that information to improve on existing models of evaporation. It turns out that the naive expectation value of the lifetime with the fermionic correction is infinite. This is a logarithmic divergence that arises from very large excursions in the separation between the branes near the threshold for classical evaporation. The adiabatic behavior lets us get some analytic control of the dynamics in this regime to get this estimate. This divergence is cutoff in the quantum theory due to quantization of the adiabatic parameter, resulting in a long lifetime of the bound state, with a parametric dependence of order log(1/ℏ).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Isotherm Modeling and Techno-Economic Analysis of Contactor Technologies for New Tetraamine-Appended MOF for NGCC Applications

Recently, a family of tetraamine-functionalized metal-organic frameworks (MOF) has been reported as promising sorbent materials for capturing CO2 from flue gas conditions relevant to natural gas combined cycle (NGCC) applications. The main advantages of these materials are their two-step cooperative CO2 adsorption, which gives rise to unusual two step-shaped CO2 adsorption profiles and their high thermal stability. This work presents the modelling of the two-transition isotherm of the tetraamine-appended MOF, N,N'-bis(3-aminopropyl)-1,4-diaminobutane (3-4-3)-appended Mg2(dobpdc), and the techno-economic analysis (TEA) of carbon capture processes utilizing this sorbent. Due to the unusual isotherm shapes of the experimental CO2 adsorption data for tetraamine-appended Mg2(dobpdc) and the strong nonlinearity of CO2 loading with respect to temperature and pressure, we tested two different models which use logistic functions for representing the different isotherm behaviors in the different pressure ranges. The first model uses the quadratic isotherm model in the low-pressure region, the Langmuir isotherm model in the middle pressure range, and the dual site Langmuir isotherm model in the high-pressure range. To model the transition between regions we used the arctangent functions independent of temperature and the thermal effect was accounted by using the Clausius-Clapeyron relation]. The second model is an extension of the weighted dual-site Langmuir isotherm model presented by Hughes et al.. In this extended model, the dual-site Langmuir isotherm is employed in the three transition regions, using temperature-dependent logistic functions to activate or deactivate the isotherm model in the low, middle, and high-pressure ranges. Both models fit the experimental data quite well with root mean squared errors (RMSE’s) of 0.41 and 0.17 for model 1 and model 2, respectively. Since model 2 resulted in a lower RMSE, it was leveraged for the development of the gas/solid contactor models used by the TEA. Specifically, two different contactor models, an axial-flow fixed bed and moving bed contactor, were developed as part of this work. These models are dynamic, pressure-driven, and consist of mass, energy, and momentum conservation equations. A kinetic model was also developed by performing parameter estimation using experimental fixed bed breakthrough data. These models are then used to simulate CO2 capture processes from the flue gas generated from a ~600 gross MW NGCC power plant. A cost model was developed which considers the capital cost of the reactors and the significant operating costs such as steam and electricity. Using NETL’s Framework for the Optimization and Quantification of Uncertainty of Uncertainty and Surrogates tool (FOQUS), which has the capability of linking models built using numerous modelling platforms with derivative-free optimization solvers, a techno-economic optimization of the carbon capture processes was performed which minimizes the cost of capture.

Caballero, Daison↗

Parameter identification and modeling of longitudinal aerodynamics

Using a comprehensive flight test database and a parameter identification software program produced at NASA Ames Research Center, a math model of the longitudinal aerodynamics of the Harrier aircraft was formulated. The identification program employed the equation error method using multiple linear regression to estimate the nonlinear parameters. The formulated math model structure adhered closely to aerodynamic and stability/control theory, particularly with regard to compressibility and dynamic manoeuvring. Validation was accomplished by using a three degree-of-freedom nonlinear flight simulator with pilot inputs from flight test data. The simulation models agreed quite well with the measured states. It is important to note that the flight test data used for the validation of the model was not used in the model identification.

Aksteter, J. W.↗

Modeling the Launch Abort Vehicle's Subsonic Aerodynamics from Free Flight Testing

An investigation into the aerodynamics of the Launch Abort Vehicle for NASA's Constellation Crew Launch Vehicle in the subsonic, incompressible flow regime was conducted in the NASA Langley 20-ft Vertical Spin Tunnel. Time histories of center of mass position and Euler Angles are captured using photogrammetry. Time histories of the wind tunnel's airspeed and dynamic pressure are recorded as well. The primary objective of the investigation is to determine models for the aerodynamic yaw and pitch moments that provide insight into the static and dynamic stability of the vehicle. System IDentification Programs for AirCraft (SIDPAC) is used to determine the aerodynamic model structure and estimate model parameters. Aerodynamic models for the aerodynamic body Y and Z force coefficients, and the pitching and yawing moment coefficients were identified.

Hartman, Christopher L.↗

An Interactive MATLAB Program for Fitting Transfer Functions to Frequency Responses

A computer program called FRFit (Frequency Response Fitting) is described for fitting single-input single-output transfer function models to empirical frequency response data. The program is interactive in that the user specifies ``elementary factors'' (gain, delay, pure differentiators and integrators, and first- and second-order zeros and poles) by entering numerical values or moving sliders in a graphical user interface. A nonlinear optimization can then be performed to obtain maximum likelihood estimates of transfer function parameters and uncertainties to provide feedback on the modeling and refine estimates. Several examples are discussed, including the identification of aircraft pitch dynamics from simulation data and data reported in the literature, approximating Theodorsen's function of unsteady aerodynamics, and obtaining a reduced-order model of a computational fluid dynamics code describing the unsteady aerodynamics around an aeroelastic wing. The program has some usefulness as a teaching aid, and can be applied to model structure determination, reduced-order modeling, preliminary analysis, and simple system identification problems. The program was written in MATLAB and is planned for public release through the NASA Software Catalog.

System identification↗

Reduced-order model to estimate safe flammable refrigerant charge limits for a one-room air conditioner

The Heating, Ventilation, Air-Conditioning and Refrigeration (HVAC&R) industry is currently transitioning away from high global warming potential (GWP) hydrofluorocarbon (HFC) refrigerants. However, many lower-GWP alternatives are flammable. Current safety codes restrict use of all flammable refrigerants; revision of the codes is needed to facilitate wider use of flammable alternatives. Hence, information is needed to enable reliable estimates of safe charge limits for flammable refrigerant options for different HVAC&R applications. A reduced-order model (ROM) was developed and validated to help estimate safe charge limits based on computational fluid dynamics (CFD) simulations of refrigerant releases in a single room for a range of parameters. Validation of the ROM vs. the CFD results showed good general agreement for mean refrigerant concentration. The normalized mean bias error (NMBE) was 0.67% with a maximum error of 15%.

Edwards, Dean↗

Active Aeroelastic Wing Aerodynamic Model Development and Validation for a Modified F/A-18A

A new aerodynamic model has been developed and validated for a modified F/A-18A used for the Active Aeroelastic Wing (AAW) research program. The goal of the program was to demonstrate the advantages of using the inherent flexibility of an aircraft to enhance its performance. The research aircraft was an F/A-18A with wings modified to reduce stiffness and a new control system to increase control authority. There have been two flight phases. Data gathered from the first flight phase were used to create the new aerodynamic model. A maximum-likelihood output-error parameter estimation technique was used to obtain stability and control derivatives. The derivatives were incorporated into the National Aeronautics and Space Administration F-18 simulation, validated, and used to develop new AAW control laws. The second phase of flights was used to evaluate the handling qualities of the AAW aircraft and the control law design process, and to further test the accuracy of the new model. The flight test envelope covered Mach numbers between 0.85 and 1.30 and dynamic pressures from 600 to 1250 pound-force per square foot. The results presented in this report demonstrate that a thorough parameter identification analysis can be used to improve upon models that were developed using other means. This report describes the parameter estimation technique used, details the validation techniques, discusses differences between previously existing F/A-18 models, and presents results from the second phase of research flights.

Cumming, Stephen B.↗

Active Aeroelastic Wing Aerodynamic Model Development and Validation for a Modified F/A-18A Airplane

A new aerodynamic model has been developed and validated for a modified F/A-18A airplane used for the Active Aeroelastic Wing (AAW) research program. The goal of the program was to demonstrate the advantages of using the inherent flexibility of an aircraft to enhance its performance. The research airplane was an F/A-18A with wings modified to reduce stiffness and a new control system to increase control authority. There have been two flight phases. Data gathered from the first flight phase were used to create the new aerodynamic model. A maximum-likelihood output-error parameter estimation technique was used to obtain stability and control derivatives. The derivatives were incorporated into the National Aeronautics and Space Administration F-18 simulation, validated, and used to develop new AAW control laws. The second phase of flights was used to evaluate the handling qualities of the AAW airplane and the control law design process, and to further test the accuracy of the new model. The flight test envelope covered Mach numbers between 0.85 and 1.30 and dynamic pressures from 600 to 1250 pound-force per square foot. The results presented in this report demonstrate that a thorough parameter identification analysis can be used to improve upon models that were developed using other means. This report describes the parameter estimation technique used, details the validation techniques, discusses differences between previously existing F/A-18 models, and presents results from the second phase of research flights.

Cumming, Stephen B.↗

Impacts of Light Use Efficiency and fPAR Parameterization on Gross Primary Production Modeling

This study examines the impact of parameterization of two variables, light use efficiency (LUE) and the fraction of absorbed photosynthetically active radiation (fPAR or fAPAR), on gross primary production(GPP) modeling. Carbon sequestration by terrestrial plants is a key factor to a comprehensive under-standing of the carbon budget at global scale. In this context, accurate measurements and estimates of GPP will allow us to achieve improved carbon monitoring and to quantitatively assess impacts from cli-mate changes and human activities. Spaceborne remote sensing observations can provide a variety of land surface parameterizations for modeling photosynthetic activities at various spatial and temporal scales. This study utilizes a simple GPP model based on LUE concept and different land surface parameterizations to evaluate the model and monitor GPP. Two maize-soybean rotation fields in Nebraska, USA and the Bartlett Experimental Forest in New Hampshire, USA were selected for study. Tower-based eddy-covariance carbon exchange and PAR measurements were collected from the FLUXNET Synthesis Dataset. For the model parameterization, we utilized different values of LUE and the fPAR derived from various algorithms. We adapted the approach and parameters from the MODIS MOD17 Biome Properties Look-Up Table (BPLUT) to derive LUE. We also used a site-specific analytic approach with tower-based Net Ecosystem Exchange (NEE) and PAR to estimate maximum potential LUE (LUEmax) to derive LUE. For the fPAR parameter, the MODIS MOD15A2 fPAR product was used. We also utilized fAPAR chl, a parameter accounting for the fAPAR linked to the chlorophyll-containing canopy fraction. fAPAR chl was obtained by inversion of a radiative transfer model, which used the MODIS-based reflectances in bands 1-7 produced by Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm. fAPAR chl exhibited seasonal dynamics more similar with the flux tower based GPP than MOD15A2 fPAR, especially in the spring and fall at the agricultural sites. When using the MODIS MOD17-based parameters to estimate LUE, fAPAR chl generated better agreements with GPP (r2= 0.79-0.91) than MOD15A2 fPAR (r2= 0.57-0.84).However, underestimations of GPP were also observed, especially for the crop fields. When applying the site-specific LUE max value to estimate in situ LUE, the magnitude of estimated GPP was closer to in situ GPP; this method produced a slight overestimation for the MOD15A2 fPAR at the Bartlett forest. This study highlights the importance of accurate land surface parameterizations to achieve reliable carbon monitoring capabilities from remote sensing information.

grossprimary↗