Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic State 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 1,063 records · Page 59

Proton Diffusion in Liquid 1,2,3-Triazole Studied by Incoherent Quasi-Elastic Neutron Scattering

Improving the proton transport in polymer electrolytes impacts the performance of next-generation solid-state batteries. However, little is known about proton conductivity in nonaqueous systems due to the lack of an appropriate level of fundamental understanding. Here, we studied the proton transport in small molecules with dynamic hydrogen bonding, 1,2,3-triazole, as a model system of proton hopping in a nonaqueous environment using incoherent quasi-elastic neutron scattering. By using the jump-diffusion model, we identified the elementary jump-diffusion motion of protons at a much shorter length scale than those by nuclear magnetic resonance and impedance spectroscopy for the estimated long-range diffusion. In addition, a spatially restricted diffusive motion was observed, indicating that proton motion in 1,2,3-triazole is complex with various local correlated dynamics. In conclusion, these correlated dynamics will be important in elucidating the nature of the proton dynamics in nonaqueous systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

COFS 1: Beam dynamics and control technology overview

The Control of Flexible Structures (COFS) 1 Project provides the invaluable opportunity to test, validate, and measure the effectiveness of theories, structural concepts, control systems, and flight certification processes for future missions through a research program focusing on multiple issues in large flexible structures, dynamics, and controls. The COFS 1 Project consists of a series of ground and flight activities building progressively from modeling and dynamic characterization of large space systems to the more complex issues of flexible-body control. The program objectives are to: determine the degree to which theory and ground testing can predict flight performance of next-generation low-frequency structures; evaluate structural fidelity of representative next-generation large deployable precision structure; assess math modeling requirements for large lightweight complex systems on which ground test results are questionable; determine degree to which scale model analysis and tests can be correlated to full-scale performance; evaluate system identification and state estimation algorithms on complex lightweight structures in the space environment; evaluate and verify controls/structures modeling capability; evaluate control laws and control systems; and evaluate damping effects in micro-g environment.

Allen, John L.↗

Verification of stability and unraveling the electronic and physical properties of bulk and (001)-surfaces of newly synthesized Ti 2 ZnX (X = C, N) MAX phases

MAX phase family has been extended by the addition of late transition metals at the A-site with the expectation of diverse functional properties. Here, we present our systematic density functional investigation on the thermodynamic and phonon stabilities, elastic properties, including elastic constants, elastic moduli and elastic anisotropy of newly synthesized Ti 2 ZnX (X = C, N) phases in comparison with conventional Ti 2 AlX (X = C, N). Due to the smaller size of N as compared to C, the unit cell dimension is reduced when C atoms are replaced by N atoms at the X-site. Furthermore, the Ti 2 ZnC and Ti 2 ZnN are stable at the equilibrium volume of 110.84 Å 3 and 105.70 Å 3 . The thermodynamic, mechanical and dynamical stabilities are validated by estimating the formation energies, elastic constants and phonon dispersions, respectively. The elastic properties of Ti 2 ZnN are less anisotropic as compared to those of Ti 2 ZnC. To understand the thin-film characteristics in Ti 2 ZnX, the surface properties with (001)-terminated slabs are investigated. Both Ti 2 ZnX bulk and (001)-surfaces exhibit metal-like electronic structures. There is a strong covalent bonding between Ti-X and Ti-Zn atoms confirmed by the charge density map and Mulliken population analysis. Additional states are generated at the Fermi level (EF) due to the unusual d-p states hybridization between Ti and Zn atoms. The anisotropy in chemical bonding is confirmed by the cleavage energy difference between Ti-X and Ti-Zn. Here, Ti(X)-001 and Zn-001 terminations are stable surfaces; however, in terms of chemical potentials, Zn-001 termination is the most favourable in Ti 2 ZnX.

36 MATERIALS SCIENCE↗

Evaluation of obstacle modelling approaches for resource assessment and small wind turbine siting: case study in the northern Netherlands

Abstract. Growth in adoption of distributed wind turbines for energy generation is significantly impacted by challenges associated with siting and accurate estimation of the wind resource. Small turbines, at hub heights of 40 m or less, are greatly impacted by terrestrial obstacles such as built structures and vegetation that can cause complex wake effects. While some progress in high-fidelity complex fluid dynamics (CFD) models has increased the potential accuracy for modelling the impacts of obstacles on turbulent wind flow, these models are too computationally expensive for practical siting and resource assessment applications. To understand the efficacy of available models in situ, this study evaluates classic and commonly used methods alongside new state-of-the-art lower-order models derived from CFD simulations and machine learning approaches. This evaluation is conducted using a subset of an extensive original dataset of measurements from more than 300 operational wind turbines in the northern Netherlands. The results show that data-driven methods (e.g. machine learning and statistical modelling) are most effective at predicting production at real sites with an average error in annual energy production of 2.5 %. When sufficient data may not be available de novo to support these data-driven approaches, models derived from high-fidelity simulations show promise and reliably outperform classic methods. On average these models have 6.3 %–11.5 % error compared with 26 % for classic methods and 27 % baseline error for reanalysis data without obstacle correction. While more performant on average, these methods are also sensitive to the quality of obstacle descriptions and reanalysis inputs.

17 WIND ENERGY↗

Dynamic compression and volatile release of carbonates

Particle velocity profiles upon shock compression and isentropic releases were measured for polycrystalline calcite. The Solenhofen limestone release paths lie, close to the Hugoniot. Calcite 3 to 2 transition, upon release, was observed, but rarefaction shocks were not detected. The equation of state is used to predict the fraction of material devolatilized upon isentropic release as a function of shock pressure. The effect of ambient partial pressure of CO2 on the calculations is demonstrated and considered in models of atmospheric evolution by impact induced mineral devolatilization. The radiative characteristics of shocked calcite indicate that localization of thermal energy occurs under shock compression. Shock entropy calculations result in a minimum estimate of 90% devolatilization upon complete release from 10 GPa. Isentropic release paths from calculated continuum Hugoniot temperatures cross into the CaO (solid) + CO2 (vapor) field at improbably low pressures. It is found that release paths from measured shock temperatures cross into the melt plus vapor field at pressures greater than .5 GPa, which suggests that devolatilization is initiated at the shear banding sites.

Tyburczy, J. A.↗

Dynamic compression and volatile release of carbonates

Particle velocity profiles upon shock compression and isentropic releases were measured for polycrystalline calcite. The Solenhofen limestone release paths lie close to the Hugoniot. Calcite 3 to 2 transition, upon release was observed, but rarefaction shocks were not detected. The equation of state is used to predict the fraction of material devolatilized upon isentropic release as a function of shock pressure. The effect of ambient partial pressure of CO2 on the calculations is demonstrated and considered in models of atmospheric evolution by impact induced mineral devolatilization. The radiative characteristics of shocked calcite indicate that localization of thermal energy occurs under shock compression. Shock entropy calculations result in a minimum estimate of 90 percent devolatilization upon complete release from 10 GPa. Isentropic release paths from calculated continuum Hugoniot temperatures cross into the CaO (solid) + CO2 (vapor) field at improbably low pressures. It is found that release paths from measured shock temperatures cross into the melt plus vapor field at pressures greater than .5 GPa, which suggests that devolatilization is initiated at the shear banding sites.

Tyburczy, J. A.↗

A coupled channel study of HN2 unimolecular decay based on a global ab initio potential surface

The unimolecular decay lifetimes of several vibrational states of HN2 are determined on the basis of an accurate coupled channel dynamics study using a global analytical potential surface. The surface reproduces the ab initio points with an rms error of 0.08 kcal/mol for energies below 20 kcal/mol. Modifications to the potential that describe the effect of improving the basis set in the ab initio calculations are provided. Converged coupled channel calculations are performed for the ground rotational state of HN2 to determine the lifetimes of the lowest ten vibrational states. Only the ground vibrational state (000) and first excited bend (001) are found to have lifetimes longer than 1 ps. The lifetimes of these states are estimated at 3 x 10 to the -9th and 2 x 10 to the -10th s, respectively. Variation of these results with quality of the ab initio calculations is not more than a factor of 5.

Koizumi, Hiroyasu↗

United States Forest Disturbance Trends Observed Using Landsat Time Series

Disturbance events strongly affect the composition, structure, and function of forest ecosystems; however, existing U.S. land management inventories were not designed to monitor disturbance. To begin addressing this gap, the North American Forest Dynamics (NAFD) project has examined a geographic sample of 50 Landsat satellite image time series to assess trends in forest disturbance across the conterminous United States for 1985-2005. The geographic sample design used a probability-based scheme to encompass major forest types and maximize geographic dispersion. For each sample location disturbance was identified in the Landsat series using the Vegetation Change Tracker (VCT) algorithm. The NAFD analysis indicates that, on average, 2.77 Mha/yr of forests were disturbed annually, representing 1.09%/yr of US forestland. These satellite-based national disturbance rates estimates tend to be lower than those derived from land management inventories, reflecting both methodological and definitional differences. In particular the VCT approach used with a biennial time step has limited sensitivity to low-intensity disturbances. Unlike prior satellite studies, our biennial forest disturbance rates vary by nearly a factor of two between high and low years. High western US disturbance rates were associated with active fire years and insect activity, while variability in the east is more strongly related to harvest rates in managed forests. We note that generating a geographic sample based on representing forest type and variability may be problematic since the spatial pattern of disturbance does not necessarily correlate with forest type. We also find that the prevalence of diffuse, non-stand clearing disturbance in US forests makes the application of a biennial geographic sample problematic. Future satellite-based studies of disturbance at regional and national scales should focus on wall-to-wall analyses with annual time step for improved accuracy.

Masek, Jeffrey G.↗

Impact of Different Thermal Gradients on the Dynamics of Cylindrical Lithium-ion Cells Subject to Accelerated Aging and on Module Performance

This study investigates the impacts of applying different thermal gradient patterns to cylindrical lithium-ion cells in a module on cell dynamics (temperatures, current flows, state of charge), module performance (evolution of resistance, capacity, and energy versus cycle number), and module lifetime. The thermal gradients were generated using cooling plates (CPs) with three different flow-field designs, namely, straight, perpendicular, and U-turn. The study uses computational fluid dynamics (CFD), the pseudo-two-dimensional (P2D) battery model, capacity loss and increased impedance due to the growth of a solid-electrolyte-interphase, and the electric current distribution from module terminals to cells that depends on the series-parallel electrical connections among the cells. The impact of the thermal gradient (resulting from the CP designs) on the variability in resistance, current, state of charge, and voltage among the cells was analyzed and linked to differences in the module's performance. Applying a thermal gradient to parallel-connected strings of series-connected cells led to variation in the current through each parallel string and an imbalance in the voltage of series-connected cells. Module performance is poorer when the thermal gradient causes a voltage imbalance than when it causes a current imbalance. Module performance becomes the worst when both current variation and voltage imbalance happen together. For instance, the module's lifetime (estimated as reaching 80% of its initial capacity) varied by 5% to 17.5%, depending on the magnitude and pattern of the imposed thermal gradient. As the relative orientation between thermal gradients and cells' electrical connectivity influences the module's performance, appropriate consideration should be given to the choice of the CP, especially if large thermal gradients are allowed.

Battery thermal management↗

Combined Satellite and Surface-Based Estimation of the Intracloud/ Cloud-to-Ground Lightning Ratio Over the Continental United States

Four years of observations from the NASA Optical Transient Detector (OTD) and Global Atmospherics National Lightning Detection Network (NLDN) are combined to determine the geographic distribution of the intracloud/cloud to ground lightning ratio over the continental United States. The mean ratio over this region is 2.64-2.94, with a standard deviation of 1.1-1.3 and anomalies as low as 1.0 or less over the Rocky and Appalachian mountains and as high as 8-10 in the central-upper midwest. There is some indication that the ratio covaries with ground elevation, although the relationship is nonunique. Little evidence is found to support a latitudinal covariance, despite significant variation in the climatological mean tropopause pressure over the latitudes considered. The dynamic range of local variability is comparable to the range of values cited by previous studies for latitudinal variation from the deep tropics to midlatitudes. Local high anomalies of this ratio in the midwest are coincident with anomalies in the climatological percentage of positive CG occurrence, as well as in the occurrence of large positive CGs characteristic of MCS convective and trailing stratiform regions. This suggests that storm type, morphology and level of organization may dominate over environmental cofactors in the local determination of this ratio.

Boccippio, D. J.↗

Detecting outbreaks using a spatial latent field

In this paper, we present a method for estimating the infection-rate of a disease as a spatial-temporal field. Our data comprises time-series case-counts of symptomatic patients in various areal units of a region. We extend an epidemiological model, originally designed for a single areal unit, to accommodate multiple units. The field estimation is framed within a Bayesian context, utilizing a parameterized Gaussian random field as a spatial prior. We apply an adaptive Markov chain Monte Carlo method to sample the posterior distribution of the model parameters condition on COVID-19 case-count data from three adjacent counties in New Mexico, USA. Our results suggest that the correlation between epidemiological dynamics in neighboring regions helps regularize estimations in areas with high variance (i.e., poor quality) data. Using the calibrated epidemic model, we forecast the infection-rate over each areal unit and develop a simple anomaly detector to signal new epidemic waves. Our findings show that anomaly detector based on estimated infection-rates outperforms a conventional algorithm that relies solely on case-counts.

Safta, Cosmin [Sandia National Laboratories (SNL-C↗

Revisiting the pseudo-supercritical path method: An improved formulation for the alchemical calculation of solid–liquid coexistence

Alchemical free energy calculations via molecular dynamics have been applied to obtain thermodynamic properties related to solid–liquid equilibrium conditions, such as melting points. In recent years, the pseudo-supercritical path (PSCP) method has proved to be an important approach to melting point prediction due to its flexibility and applicability. In the present work, we propose improvements to the PSCP alchemical cycle to make it more compact and efficient through a concerted evaluation of different potential energies. The multistate Bennett acceptance ratio (MBAR) estimator was applied at all stages of the new cycle to provide greater accuracy and uniformity, which is essential concerning uncertainty calculations. In particular, for the multistate expansion stage from solid to liquid, we employed the MBAR estimator with a reduced energy function that allows affine transformations of coordinates. Free energy and mean derivative profiles were calculated at different cycle stages for argon, triazole, propenal, and the ionic liquid 1-ethyl-3-methyl-imidazolium hexafluorophosphate. Comparisons showed a better performance of the proposed method than the original PSCP cycle for systems with higher complexity, especially the ionic liquid. Furthermore, a detailed study of the expansion stage revealed that remapping the centers of mass of the molecules or ions is preferable to remapping the coordinates of each atom, yielding better overlap between adjacent states and improving the accuracy of the methodology.

25 ENERGY STORAGE↗

Digital adaptive model following control

Development of an implementable digital adaptive control system requires a control law that can be efficiently adjusted during system operation and an identification algorithm that is capable of supplying parameter estimates at an accuracy and rate specified by the controller characteristics. To this effect a control law has been derived based upon the minimization of a weighted combination of the instantaneous squared error between the states of a linear plant and model at the next stage and the present control energy. The result is a control signal that is a linear combination of plant states, model states, and any model inputs. Explicit formula are available for computing the pertinent gains. To illustrate the utility of such a control law, its application to typical linearized lateral aircraft dynamics is currently being considered.

Alag, G.↗

A methodology for designing robust multivariable nonlinear feedback control systems

This paper presents an outline of a methodology for the design of nonlinear dynamic compensators for nonlinear multivariable systems with guarantees of closed-loop stability, robustness, and performance. The method is an extension of the Linear-Quadratic-Gaussian with Loop-Transfer-Recovery (LQG/LTR) methodology for linear systems, thus hinging upon the idea of constructing an approximate nonlinear inverse operator for the plant. A major feature of the method is a unification of both the state-space and Input-Output formulations. The major results of the paper are: (1) Properties of the solution of certain optimal control problems; (2) New results on nonlinear estimation, including a new nonlinear observer with certain guaranteed properties; (3) The Loop Operator Recovery (LOR) procedure, and (4) Inclusion of command-following and disturbance rejection with a goal of a closed-loop stable, robust system.

Grunberg, D. B.↗

Adjoints and Low-rank Covariance Representation

Quantitative measures of the uncertainty of Earth System estimates can be as important as the estimates themselves. Second moments of estimation errors are described by the covariance matrix, whose direct calculation is impractical when the number of degrees of freedom of the system state is large. Ensemble and reduced-state approaches to prediction and data assimilation replace full estimation error covariance matrices by low-rank approximations. The appropriateness of such approximations depends on the spectrum of the full error covariance matrix, whose calculation is also often impractical. Here we examine the situation where the error covariance is a linear transformation of a forcing error covariance. We use operator norms and adjoints to relate the appropriateness of low-rank representations to the conditioning of this transformation. The analysis is used to investigate low-rank representations of the steady-state response to random forcing of an idealized discrete-time dynamical system.

Tippett, Michael K.↗

HLPW-4/GMGW-3: Mesh Adaptation for RANS Technology Focus Group Workshop Summary

The AIAA workshop series helps guide the Computational Fluid Dynamics (CFD) community through publicly available test cases and impartial evaluation to establish the state-of-the-art in applied CFD. A summary of the Mesh Adaptation Technical Focus Group submissions for the joint Fourth CFD High Lift Prediction and Third Geometry & Mesh Generation Workshop is provided. The Common Research Model High-Lift variant (CRM-HL) is the target of this investigation. A 2D high-lift airfoil solution verification case was extracted from the 3D CRM-HL. Low variation is shown between 2D solutions with mesh adaptation to control estimated solution interpolation and output errors. As expected, variation for these mesh adaptation techniques is larger for the 3D CRM-HL than the 2D airfoil. More importantly, the 3D mesh adapted results have a lower variation than expert-crafted meshes of equal or larger size. Suspected multiple solutions to Reynolds-averaged Navier-Stokes with the Spalart-Allmaras turbulence model appear to resolve into a common solution with mesh resolution.

HLPW-4↗

Quantum utility-scale error mitigation for quantum quench dynamics in Heisenberg spin chains

Here, we implement a quantum error mitigation method termed self-mitigation, which is comparable to zero-noise extrapolation, at large scales to achieve quantum utility on near-term, noisy quantum computers. We investigate the effectiveness of several quantum error mitigation strategies, including self-mitigation, by simulating quantum quench dynamics for Heisenberg spin chains with system sizes up to 104 qubits using IBM quantum processors. In particular, we discuss the limitations of zero-noise extrapolation and the advantages offered by self-mitigation at large scales. The self-mitigation method demonstrates stable accuracy with large systems of 104 qubits comprising more than 3,000 CNOT gates. Also, we combine the discussed quantum error mitigation methods with practical entanglement entropy measuring methods, and it shows a good agreement with the theoretical estimation. Our study illustrates the usefulness of near-term noisy quantum hardware in examining the quantum quench dynamics of many-body systems at large scales and lays the groundwork for surpassing classical simulations with quantum methods prior to the development of fault-tolerant quantum computers.

97 MATHEMATICS AND COMPUTING↗