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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Identification and Online Updating of Dynamic Models for Demand Response of an Industrial Air Separation Unit

Demand-response operation of air separation units requires frequent changes in production rate(s), and scheduling calculations must explicitly consider process dynamics to ensure feasibility of the solutions. To this end, scale-bridging models (SBMs) approximate the scheduling-relevant dynamics of a process and its controller in a low-order representation. In contrast to previous works that have employed nonlinear SBMs, this paper proposes linear SBMs, developed using time-series analysis, to facilitate online scheduling computations. Using a year-long industrial dataset, we nd that compact linear SBMs are suitable approximations over typical scheduling horizons, but that their accuracies are unpredictable over time. We introduce a strategy for online updating of the SBMs, based on Kalman ltering schemes for online parameter estimation. The approach greatly improves the accuracy of SBM predictions and will enable the use of linear SBM-based demand-response scheduling in the future.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Greybox Thermal Parameter Identification of Electric Machine Stators

The parameters of electric machine thermal equivalent circuit networks are difficult to predict due to material and manufacturing uncertainties. In this paper, a Greybox system identification approach is used to identify parameters of electric machine stator lumped parameter thermal networks (LPTNs). LPTNs provide a low order, computationally efficient, dynamic model of temperatures at specific locations. Second and third order LPTN model structures are defined as state space equations with stator thermal parameters to be identified. To test the Greybox electric machine stator thermal system identification, five stator motorette prototypes were constructed with controlled variations in slot fill and slot liner thickness. The variation in the motorette thermal parameters and thermal time constants are detected using the Greybox identification. Special attention is given to the impact of sampling rate and Greybox data record length on parameter estimation accuracy.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Enhancing Gaussian Process Surrogates for Optimization and Posterior Approximation via Random Exploration

This paper proposes novel noise-free Bayesian optimization strategies that rely on a random exploration step to enhance the accuracy of Gaussian process surrogate models. The new algorithms retain the ease of implementation of the classical GP-UCB algorithm, but the additional random exploration step accelerates their convergence, nearly achieving the optimal convergence rate. Furthermore, to facilitate Bayesian inference with intractable likelihoods, we propose to utilize optimization iterates for maximum a posteriori estimation to build a Gaussian process surrogate model for the unnormalized log-posterior density. We provide bounds for the Hellinger distance between the true and the approximate posterior distributions in terms of the number of design points. We demonstrate the effectiveness of our Bayesian optimization algorithms in nonconvex benchmark objective functions, in a machine learning hyperparameter tuning problem, and in a black-box engineering design problem. The effectiveness of our posterior approximation approach is demonstrated in two Bayesian inference problems for parameters of dynamical systems.

Bayesian inference↗
NASA NTRS1978

Transonic static and dynamic stability characteristics of a finned projectile configuration

Static and dynamic stability tests were made of a finned projectile configuration with the aft-mounted fins arranged in a cruciform pattern. The tests were made at free stream Mach numbers of 0.7, 0.9, 1.1, and 1.2 in the Langley 8-foot transonic pressure tunnel. Some of the parameters measured during the tests were lift, drag, pitching moment, pitch damping, and roll damping. Configurations tested included the body with undeflected fins, the body with various fin deflections for control, and the body with fins removed. Theoretical estimates of the stability derivatives were made for the fins on configuration.

Boyden, R. P.↗
NASA NTRS1981

Inferences about the solar wind dynamics from observed distributions of electrons and ions

Recent measurements, mostly from Helios and Voyager, were reviewed. The observational material was compiled into a series of radial profiles for various macroscopic plasma parameters. Several examples for experimentally determined three-dimensional distribution functions of protons are shown. The physical properties of the solar wind electrons were surveyed. The still unresolved question of the driving mechanism of the solar wind is considered. The problem of the energy balance in pure, high velocity streams emanating from solar coronal holes was outlined. It is suggested that the suprathermal electrons may play an important role in driving the wind. Numerical estimates, referring to the lower corona, are offered in support of this hypothesis.

Olbert, S.↗
NASA NTRS1995

Interactive flutter analysis and parametric study for conceptual wing design

An interactive computer program was developed for wing flutter analysis in the conceptual design stage. The objective was to estimate the flutter instability boundary of a flexible cantilever wing, when well defined structural and aerodynamic data are not available, and then study the effect of change in Mach number, dynamic pressure, torsional frequency, sweep, mass ratio, aspect ratio, taper ratio, center of gravity, and pitch inertia, to guide the development of the concept. The software was developed on MathCad (trademark) platform for Macintosh, with integrated documentation, graphics, database and symbolic mathematics. The analysis method was based on nondimensional parametric plots of two primary flutter parameters, namely Regier number and Flutter number, with normalization factors based on torsional stiffness, sweep, mass ratio, aspect ratio, center of gravity location and pitch inertia radius of gyration. The plots were compiled in a Vaught Corporation report from a vast database of past experiments and wind tunnel tests. The computer program was utilized for flutter analysis of the outer wing of a Blended Wing Body concept, proposed by McDonnell Douglas Corporation. Using a set of assumed data, preliminary flutter boundary and flutter dynamic pressure variation with altitude, Mach number and torsional stiffness were determined.

Mukhopadhyay, Vivek↗
DOE OSTI2026

An Ab Initio Molecular Dynamics Study of Key Thermodynamic Input Parameters for Computer Simulation of U-6Nb Solidification

The key to metallic fuel development is the fabrication of uranium metal and alloys into fuel forms. U-Nb alloys are one of the best candidates for a metallic fuel alloy with high-temperature strength sufficient to support the core, acceptable nuclear properties, good fabricability, and compatibility with usable coolant media. Melt processing has been a key component of the metallic fuel cycle, and process models require thermophysical parameters at elevated temperatures, particularly above the melting temperatures, regarding which experimental data are scarce, for accurate simulations and process development. By means of ab initio density-functional theory (DFT) quantum molecular dynamics (QMD), we have calculated the main thermophysical parameters—the density, thermal expansion coefficient, specific heat, thermal conductivity, melting temperature, latent heat of fusion, and viscosity—used in the modeling of the U-6 wt.% Nb alloy casting. The melting temperature of the U-6 wt.% Nb alloy at ambient pressure is obtained by means of QMD simulations using the Z-method. The ambient volume change and latent heat of melting of U-6 wt.% Nb are also derived from QMD simulations in conjunction with analytical fitting for the energy and pressure. The thermal conductivity for the solid U-Nb alloy is calculated from the semi-classical Boltzmann transport equation combined with an estimate of the electron relaxation time obtained from DFT simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

A New Computational Method for Membrane Compressibility: Bilayer Mechanical Thickness Revisited

Because lipid bilayers can bend and stretch in ways similar to thin elastic sheets, physical models of bilayer deformation have utilized mechanical constants such as the moduli for bending rigidity (κ C ) and area compressibility (K A ). However, the use of these models to quantify the energetics of membrane deformation associated with protein-membrane interactions, and the membrane response to stress is often hampered by the shortage of experimental data suitable for the estimation of the mechanical constants of various lipid mixtures. Although computational tools such as molecular dynamics simulations can provide alternative means to estimate K A values, current approaches suffer significant technical limitations. Here, in this study, we present a novel, to our knowledge, computational framework that allows for a direct estimation of K A values for individual bilayer leaflets. The theory is based on the concept of elasticity and derives K A from real-space analysis of local thickness fluctuations sampled in molecular dynamics simulations. We explore and validate the model on a large set of single and multicomponent bilayers of different lipid compositions and sizes, simulated at different temperatures. The calculated bilayer compressibility moduli agree with values estimated previously from experiments and those obtained from a standard computational method based on a series of constrained tension simulations. We further validate our framework in a comparison with an existing polymer brush model and confirm the polymer brush model’s predicted linear relationship with proportionality coefficient of 24, using elastic parameters calculated from the simulation trajectories. The robustness of the results that emerge from the method allows us to revisit the origins of the bilayer mechanical (compressible) thickness and in particular its dependence on acyl-chain unsaturation and the presence of cholesterol.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Enhancement of detachment control with simplified real-time modelling on the KSTAR tokamak

Abstract Detachment control based on ion saturation current I sat measurements from Langmuir probes (LPs) is implemented in the KSTAR tokamak and shown to be capable of following dynamic and constant target trajectories with good accuracy, in H-mode, by moderating the flow rate of nitrogen or deuterium. I sat controllers normalize I sat in order to form attachment fraction ( A frac ) as their control parameter. The KSTAR implementation of A frac control differs from previous work in that it continuously calculates a model for attached I sat and uses that as the denominator in A frac , whereas prior implementations either record peak I sat at rollover as they pass it or take estimated I s a t , r o l l o v e r as a manual input prior to the shot. The KSTAR controller therefore does not need to keep track of rollover status and keep separate targets for pre- and post-rollover states, and it can automatically adapt to changes in scenario at any time. It is also less vulnerable to noise as it will not lock in an outlier as a rollover point.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Predictions of Cell-to-Cell Propagation and Vent Gas Production in the Thermal Runaway of Lithium-Ion Battery Stacks

This work presents the thermal runaway propagation model LIM1TR (Lithium-ion Modeling with 1-D Thermal Runaway) as an efficient tool to predict different cell-to-cell thermal runaway propagation scenarios. Here, we explored the vent gas volume production and reaction duration highlighting the relationship between these parameters and thermal runaway propagation due to convection by the vented gases. Two metrics based on gas production rate and heating rate are utilized as good indicators of the start and end of thermal runaway. LIM1TR results are compared with and validated by experiments from the literature for single-cell and multicell array experiments of 5 Ah and 10 Ah cells. By accounting for intraparticle diffusion of reacting species in the electrodes, we were able to capture the general dynamics of thermal runaway propagation and estimate acceptable reaction durations compared with the experimental values. Simulation results further demonstrated that varying heating modes lead to distinct reaction durations, consistent with experimental observations. Vent gas volume predictions indicate the need to consider both full and partial oxidation of the electrolyte. The outcomes of this work are building blocks for further investigations of module-to-module propagation by vented gases through convective heat transfer.

25 ENERGY STORAGE↗
DOE OSTI2020

Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in-situ and remote sensing datasets

The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1991

Optimal interaction of indenter with inhomogeneous plate

Consideration is given to a new class of problems dealing with an optimal design of an inhomogeneous plate during dynamic penetration of the rigid indenter. The quality criterion of the process is defined by the specific mass of the target, which absorbs the given kinetic mass of the indenter. Parameters of control are expressed in terms of mechanical characteristics, i.e., distribution of density and the related hardness across the plate thickness. The maximum principle of Pontryagin is used to search for the piecewise continuous control function. With consideration of impact conditions and characteristics for a given class of material, an optimal target structure criterion was estimated for engineering applications.

Aptukov, Valery N.↗
NASA NTRS1992

The evolution equation for the flame surface density in turbulent premixed combustion

One central ingredient in flamelet models for turbulent premixed combustion is the flame surface density. This quantity conveys most of the effects of the turbulence on the rate of energy release and is obtained via a modeled transport equation, called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact, but unclosed, formulation for the turbulent Sigma-equation. In this exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the flame surface density and the turbulent flame stretch are not available from experiments and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the obvious, complementary approach to get basic information on these fundamental quantities. Three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, the effects of the Lewis number on the rate of production of flame surface area are described in great detail and meaningful comparisons with flamelet models can be performed. The analysis reveals in particular the tendency of the models to overpredict flame surface dissipation as well as their inability to reproduce variations due to thermo-diffusive phenomena. Thanks to the detailed information produced by a DNS-based analysis, this type of comparison not only underscores the shortcomings of current models but also suggests ways to improve them.

Trouve, A.↗
NASA NTRS2002

An Investigation of the Influence of Urban Areas on Rainfall Using the TRMM Satellite and a Cloud-Mesoscale Model

A recent paper by Shepherd and Pierce (in press at Journal of Applied Meteorology) used rainfall data from the Precipitation Radar on NASA's Tropical Rainfall Measuring Mission's (TRMM) satellite to identify warm season rainfall anomalies downwind of major urban areas. Data (PR) were employed to identify warm season rainfall (1998-2000) patterns around Atlanta, Montgomery, Nashville, San Antonio, Waco, and Dallas. Results reveal an average increase of approx. 28% in monthly rainfall rates within 30-60 kilometers downwind of the metropolis with a modest increase of 5.6% over the metropolis. Portions of the downwind area exhibit increases as high as 51%. The percentage changes are relative to an upwind control area. It was also found that maximum rainfall rates in the downwind impact area exceeded the mean value in the upwind control area by 48%-116%. The maximum value was generally found at an average distance of 39 km from the edge of the urban center or 64 km from the center of the city. Results are consistent with METROMEX studies of St. Louis almost two decades ago and with more recent studies near Atlanta. A convective-mesoscale model with extensive land-surface processes is currently being employed to (a) determine if an urban heat island (UHI) thermal perturbation can induce a dynamic response to affect rainfall processes and (b) quantify the impact of the following three factors on the evolution of rainfall: (1) urban surface roughness, (2) magnitude of the UHI temperature anomaly, and (3) physical size of the UHI temperature anomaly. The sensitivity experiments are achieved by inserting a slab of land with urban properties (e.g. roughness length, albedo, thermal character) within a rural surface environment and varying the appropriate lower boundary condition parameters. The study will discuss the feasibility of utilizing satellite-based rainfall estimates for examining rainfall modification by urban areas on global scales and over longer time periods. The talk also introduces very preliminary results from the modeling component of the study. Such research has implications for weather forecasting, urban planning, water resource management, and understanding human impact on the environment and climate.

Shepherd, J. Marshall↗
DOE OSTI2020

Influence of Vegetation Canopy Architecture on Dynamics of Spatiotemporal Light Distribution and its Effect on Land-Biosphere-Atmosphere Fluxes of Energy, Carbon, and Water

The dynamic spatiotemporal distribution of light and micro climatic variables (wind speed, atmospheric CO 2 concentration, humidity, pressure, and air temperature) within the plant canopy as determined by the canopy architecture play a critical role in land-biosphere-atmosphere (LBA) photosynthetic carbon uptake, transpiratory water loss and energy partitioning. Vegetation models that include effects of sun and shade flecks, leaf temperature and energy balance partitioning between sensible and latent heat fluxes, and reduced turbulent mixing in dense vegetation canopies can accurately predict canopy performance, particularly its shaded leaves, which can constitute 80% of leaf area and contribute up to 50% of carbon update and 65% of water loss. However, current vegetation models with their bulk characterization of shaded leaves poorly represent the complex three dimensional (3D) canopy architecture are insufficient to capture the spatiotemporal dynamics of shaded leaves. Here, we propose to employ an x-ray computed tomography (XCT) imaging technique in combination with an automated, high throughput, image processing pipeline (3DArch, to be developed) to obtain an explicit 3D triangulated surface mesh of a plant canopy (Fig 1). 3DArch would provide a much needed technological jump in 3D LBA modeling where computer generated artificial 3D models suffer from lack of high resolution detailed canopy architecture parameters. The 3D canopy surface mesh output from 3DArch is coupled with an explicit 3D forward ray tracing model (fastTracer) to obtain a high resolution representation of the dynamic spatiotemporal light environment within the canopy over the course of a day (Fig 1). The results from the light attenuation model are then coupled with a process based multi-scale bio-physical canopy and leaf model of photosynthesis (PlantGro) to get fluxes of carbon, water and energy at the leaf and canopy scale, in addition to estimating whole plant carbon uptake, and water loss. Our proposed model framework would be the first of its kind to provide a breakthrough capability in the field of process-based plant modeling that bridges the key technological gap in shaded leaf modeling through explicit representation of 3D canopy architecture. Computer generated artificial 3D model canopy simulations on grasses show that the canopy light distribution does not follow an exponential behavior as assumed in most 1D vegetation models, but rather exhibits a heterogeneous behavior, with a small but significant portion of interior canopy leaves experiencing peak light levels due to sun flecks while the rest of the canopy receives substantially lower light levels. These results indicate a bias in traditional models that overestimates shaded leaf contribution towards whole canopy photosynthesis by up to 15% and underestimate transpiration by up to 12% with an associated increase in sensible heat fluxes since photosynthesis, transpiration and energy partitioning are tightly coupled processes. However, these test results need to be validated using the 3D architecture obtained from a real plant imaged in 3D using 3DArch, rather than an artificial model plant. Our proposed model framework would allow us to identify inefficiencies in canopy light distribution and prescribe optimal canopy architectures (leaf optical properties, leaf area, leaf angles, row spacing, planting density, etc.) to traditional and genetic engineering based crop breeders, for maximizing target functions such as light use efficiency, water use efficiency, drought tolerance, crop yield, crop resilience etc., with important applications in biofuel crops such as sorghum, and switchgrass. Our model can also be applied for designing optimal canopy architecture that increase surface albedo as a mitigation strategy for global warming and projections of crop yield for future food and water security under climate change.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1982

Comparison of gas dynamic model with steady solar wind flow around Venus

A gas dynamic model for solar wind flow around Venus has been compared with Pioneer Venus orbiter plasma analyzer measurements from times when the solar wind flow seemed steadiest. The comparisons were made near the terminator. When the observed and model bow shock locations are matched, the model agrees fairly well with the observed parameters: the components of the flow velocity and magnetic field, and the proton number density and isotropic temperature. However, the Mach numbers required to fit the observed bow shock locations are less than 2/3 those estimated from the measured parameters of the free-stream flow, because the bow shock locations (near the terminator) are farther from Venus than expected. Sometimes the measured flow speeds appear to be retarded near the ionopause.

Mihalov, J. D.↗
NASA NTRS1985

Solar wind variations and geomagnetic storms - A study of individual storms based on high time resolution ISEE 3 data

Two independent methods are employed to determine the relationship between the parameter epsilon and total energy dissipation rate of the magnetosphere U sub T by selecting disturbed periods from the same data d set used by Baker et al. (1983). Specifically, four storms are examined in detail, since the accuracy of estimating U sub T is significantly improved during disturbed periods. The first method assumes that U sub T = M sub A exp.2- alpha(epsilon) where M sub A is the Alfven Mach number and alpha varies with time. The second method considers a linear, time-invariant dynamic system with epsilon as input and U sub T as output. This means that U sub T = W(asterisk)epsilon where asterisk is the convolution and W is a transfer function characteristic of the system. It is found that alpha values fluctuate mainly between 0 and -0.25. The transfer function analysis indicates that W often resembles a delta-function or a narrow rectangular impulse. Both results give the same implication (namely that U sub T is approximately equal to epsilon) and thus are consistent with the view that the magnetosphere is primarily a directly driven system during disturbed periods.

Akasofu, S.-I.↗
DOE OSTI2022

Feasibility of Vascular Parameter Estimation for Assessing Hypertensive Pregnancy Disorders

Abstract Hypertensive pregnancy disorders (HPDs), such as pre-eclampsia, are leading sources of both maternal and fetal morbidity in pregnancy. Noninvasive imaging, such as ultrasound (US) and magnetic resonance imaging (MRI), is an important tool for predicting and monitoring these high risk pregnancies. While imaging can measure hemodynamic parameters, such as uterine artery pulsatility and resistivity indices (PI and RI), the interpretation of such metrics for disease assessment relies on ad hoc standards, which provide limited insight to the physical mechanisms underlying the emergence of hypertensive pregnancy disorders. To provide meaningful interpretation of measured hemodynamic data in patients, advances in computational fluid dynamics can be brought to bear. In this work, we develop a patient-specific computational framework that combines Bayesian inference with a reduced-order fluid dynamics model to infer parameters, such as vascular resistance, compliance, and vessel cross-sectional area, known to be related to the development of hypertension. The proposed framework enables the prediction of hemodynamic quantities of interest, such as pressure and velocity, directly from sparse and noisy MRI measurements. We illustrate the effectiveness of this approach in two systemic arterial network geometries: an aorta with branching carotid artery and a maternal pelvic arterial network. For both cases, the model can reconstruct the provided measurements and infer parameters of interest. In the case of the maternal pelvic arteries, the model can make a distinction between the pregnancies destined to develop hypertension and those that remain normotensive, expressed through the value range of the predicted absolute pressure.

Biophysics↗