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 505 records · Page 28

Impact of control signal phase noise on qubit fidelity

As qubit decoherence times are increased and readout technologies are improved, nonidealities in the drive signals, such as phase noise, are going to represent a growing limitation to the fidelity achievable at the end of complex control pulse sequencies. Here we study the impact on fidelity of phase noise affecting reference oscillators with the help of numerical simulations, which allow us to directly take into account the interaction between the phase fluctuations in the control signals and the evolution of the qubit state. Our method is based on the generation of phase noise realizations consistent with a given power spectral density, that are then applied to the pulse carrier in simulations, with Qiskit-Dynamics, of the qubit temporal evolution. By comparing the final state obtained at the end of a noisy pulse sequence with that in the ideal case and averaging over multiple noise realizations, we estimate the resulting degradation in fidelity, and exploiting an approximate analytical representation of a carrier affected by phase fluctuations, we discuss the contributions of the different spectral components of phase noise.

Barsotti, Agata [Pisa U.]↗

Impact of dynamics, entanglement and Markovian noise on the fidelity of few-qubit digital quantum simulation

Quantum algorithms have been proposed to accelerate the simulation of the chaotic dynamical systems that are ubiquitous in the physics of plasmas. Quantum computers without error correction might even use noise to their advantage to calculate the Lyapunov exponent by measuring the Loschmidt echo fidelity decay rate. For the first time, digital Hamiltonian simulations of the quantum sawtooth map, performed on the IBM-Q quantum hardware platform, show that the fidelity decay rate of a digital quantum simulation increases during the transition from dynamical localization to chaotic diffusion in the map. The observed error per CNOT gate increases by $1.5{\times }$ as the dynamics varies from localized to diffusive, while only changing the phases of virtual RZ gates and keeping the overall gate count constant. A gate-based Lindblad noise model that captures the effective change in relaxation and dephasing errors during gate operation qualitatively explains the effect of dynamics on fidelity as being due to the localization and entanglement of the states created. Specifically, highly delocalized states that are entangled with random phases show an increased sensitivity to dephasing and, on average, a similar sensitivity to relaxation as localized states. In contrast, delocalized unentangled states show an increased sensitivity to dephasing but a lower sensitivity to relaxation. This gate-based Lindblad model is shown to be a useful benchmarking tool by estimating the effective Lindblad coherence times during CNOT gates and finding a consistent $2\unicode{x2013}3{\times }$ shorter $T_2$ time than reported for idle qubits. Thus, the interplay of the dynamics of a simulation with the noise processes that are active can strongly influence the overall fidelity decay rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Land Surface in Earth System Reanalysis

Current global atmospheric reanalysis products such as ERA-5, MERRA-2, and JRA-55 provide estimates of land surface states and fluxes, including soil moisture, soil temperature, snow mass, latent and sensible heat fluxes, and runoff. These land surface estimates are based on land surface process models and, depending on the reanalysis product, on precipitation observations or the assimilation of land surface observations of soil moisture and snow conditions from satellites and in situ networks. The land models used in the current global atmospheric reanalysis systems are generally focused on water and energy balance processes, with prescribed static land cover and prescribed climatological vegetation conditions. That is, current global atmospheric reanalysis products do not reflect changes in land use or land cover, nor do they provide dynamic estimates of the land carbon cycle. These and other gaps, including the omission of anthropogenic processes such as irrigation in the land model, need to be addressed in the development of the next generation Earth system reanalysis systems. This presentation provides an overview of the state-of-the-art and future developments of land surface estimates in global reanalysis systems.

Rolf Reichle↗

Elastic deformation as a tool to investigate watershed storage connectivity

Storage-discharge relationships and dynamic changes in storage connectivity remain key unknowns in understanding and predicting watershed behavior. In this study, we use Global Positioning System measurements of load-induced Earth surface displacement as a proxy for total water storage change in four climatologically diverse mountain watersheds in the western United States. Comparing total water storage estimates with stream-connected storage derived from hydrograph analysis, we find that each of the investigated watersheds exhibits a characteristic seasonal pattern of connection and disconnection between total and stream-connected storage. We investigate how the degree and timing of watershed-scale connectivity is related to the timing of precipitation and seasonal changes in dominant hydrologic processes. Our results show that elastic deformation of the Earth due to water loading is a powerful new tool for elucidating dynamic storage connectivity and watershed discharge response across scales in space and time.

Environmental Sciences & Ecology↗

A CFD Approach to Modeling Spacecraft Fuel Slosh

Energy dissipation and resonant coupling from sloshing fuel in spacecraft fuel tanks is a problem that occurs in the design of many spacecraft. In the case of a spin stabilized spacecraft, this energy dissipation can cause a growth in the spacecrafts' nutation (wobble) that may lead to disastrous consequences for the mission. Even in non-spinning spacecraft, coupling between the spacecraft or upper stage flight control system and an unanticipated slosh resonance can result in catastrophe. By using a Computational Fluid Dynamics (CFD) solver such as Fluent, a model for this fuel slosh can be created. The accuracy of the model must be tested by comparing its results to an experimental test case. Such a model will allow for the variation of many different parameters such as fluid viscosity and gravitational field, yielding a deeper understanding of spacecraft slosh dynamics. In order to gain a better understanding of the dynamics behind sloshing fluids, the Launch Services Program (LSP) at the NASA Kennedy Space Center (KSC) is interested in finding ways to better model this behavior. Thanks to past research, a state-of-the-art fuel slosh research facility was designed and fabricated at Embry Riddle Aeronautical University (ERAU). This test facility has produced interesting results and a fairly reliable parameter estimation process to predict the necessary values that accurately characterize a mechanical pendulum analog model. The current study at ERAU uses a different approach to model the free surface sloshing of liquid in a spherical tank using Computational Fluid Dynamics (CFD) methods. Using a software package called Fluent, a model was created to simulate the sloshing motion of the propellant. This finite volume program uses a technique called the Volume of Fluid (VOF) method to model the interaction between two fluids [4]. For the case of free surface slosh, the two fluids are the propellant and air. As the fuel sloshes around in the tank, it naturally displaces the air. Using the conservation of mass, momentum, and energy equations, as well as the VOF equations, one can predict the behavior of the sloshing fluid and calculate the forces, pressure gradients, and velocity field for the entire liquid as a function of time.

Marsell, Brandon↗

IM3/HyperFACETS Thermodynamic Global Warming (TGW) Simulation Datasets

Publication For a thorough description of the methods, see the peer-reviewed paper: Jones, A.D., Rastogi, D., Vahmani, P. et al. Continental United States climate projections based on thermodynamic modification of historical weather. Sci Data 10, 664 (2023). https://doi.org/10.1038/s41597-023-02485-5 Overview The IM3 / HyperFACETS climate simulations provide 40-year historical (1980-2019) as well as four 80-year future simulations (2020-2099) over the U.S. The future simulations are split into near (2020-2059) and far future (2060-2099) segments. The future scenarios span a range of plausible changes in future climate (both Global Circulation Model (GCM) and Representative Concentration Pathways/Shared Socioeconomic Pathway (RCP/SSP) dimensions). The simulations provide climate variables with high spatiotemporal resolution (25 hourly variables and 207 3-hourly variables at 12 km2). The datasets are generated using dynamical downscaling with the WRF (Weather Research and Forecasting) model (version 4.2.1) and therefore preserve physical consistency across variables. WRF is a state-of-the-art, fully compressible, non-hydrostatic, mesoscale numerical weather prediction model. WRF is coupled with an urban canopy model (UCM), which resolves urban surfaces. The future scenarios were developed using a thermodynamic global warming approach where past events are replayed under a range of future warming conditions. These scenarios therefore provide a perspective on potential increases in extreme event intensity, geographic scope, and duration, with previously non-extreme conditions potentially crossing new thresholds to be considered extreme by today's standards. This approach is not intended to estimate future changes in extreme event frequency that might result from changes in large-scale atmospheric dynamics. This dataset has NOT been bias corrected. A bias corrected version of selected variables is under development and will be released here when available. Scenarios Files Data for each scenario is provided in weekly NetCDF files. 25 variables are available at hourly resolution, and 207 variables are available at three-hourly resolution. Spatial resolution is 12km and spans the conterminous United States (CONUS), including some areas of Canada and Mexico, resulting in a grid of 424 by 299 cells. The spatial projection is a Lambert Conformal Conic with the following proj-string: "+proj=lcc +lat_0=40.0000076293945 +lon_0=-97 +lat_1=30 +lat_2=45 +x_0=0 +y_0=0 +R=6370000 +units=m +no_defs". The available scenarios and simulation periods are listed below: historical | 1980 - 2019 rcp45cooler | 2020 - 2059 rcp45cooler | 2060 - 2099 rcp45hotter | 2020 - 2059 rcp45hotter | 2060 - 2099 rcp85cooler | 2020 - 2059 rcp85cooler | 2060 - 2099 rcp85hotter | 2020 - 2059 rcp85hotter | 2060 - 2099 * The first year (1979, 2019, and 2059) of data within each scenario represents a model warmup period and should not be used. These are located in the `spinup_files` directory. Historical year 2020 is considered an extra year of data beyond the simulation period and can be found in the `additional_files` directory. For information on specific variables and a more in-depth discussion of methodology, please refer to the data landing page at https://tgw-data.msdlive.org. Delta Warming Files The global and CONUS warming deltas for each scenario are provided in degrees Celsius annually and monthly. Restart Files Yearly restart files are provided for each scenario which can be used to restart the WRF model at a particular point in time. Spinup Files The first year of data within each simulation period represents a model warmup period and should not be used. The files are provided here for the sake of reproducibility. Additional Files Additional years of data are provided as an extension of the historic simulation.

Jones, Andrew D.↗

Exploring an ‘ideal hill': how lithology and transport mechanisms affect the possibility of a steady state during weathering and erosion

We introduce a model of chemical reaction within hills to explore how evolving porosity (and by inference, permeability) affects flow pathways and weathering. The model consists of hydrologic and reactive-transport equations that describe alteration of ferrous minerals and feldspar. These reactions were chosen because previous work emphasized that oxygen- and acid-driven weathering affects porosity differently in felsic and mafic rocks. A parameter controlling the order of the fronts is presented. In the absence of erosion, the two reaction fronts move at different velocities and the relative depths depend on geochemical conditions and starting composition. In turn, the fronts and associated changes in porosity drastically affect both the vertical and lateral velocities of water flow. For these instances, estimates of weathering advance rates based on simple models that posit unidirectional constant-velocity advection do not apply. In the model hills, weathering advance rates diminish with time as the Darcy velocities decrease with depth. The system can thus attain a dynamical steady state at any erosion rate where the regolith thickness is constant in time and velocities of both fronts become equal to one another and to the erosion rate. The slower the advection velocities in a system, the faster it attains a steady state. For example, a massive rock with relatively fast-dissolving minerals such as diabase -- where solute transport across the reaction front mainly occurs by diffusion -- can reach a steady state more quickly than granitoid rocks in which advection contributes to solute transport. The attainment of a steady state is controlled by coupling between weathering and hydrologic processes that force water to flow horizontally above reaction fronts where permeability changes rapidly with depth and acts as a partial barrier to fluid flow.

54 ENVIRONMENTAL SCIENCES↗

Error Bounds for Dynamical Spectral Estimation

Dynamical spectral estimation is a well-established numerical approach for estimating eigenvalues and eigenfunctions of the Markov transition operator from trajectory data. Although the approach has been widely applied in biomolecular simulations, its error properties remain poorly understood. Here we analyze the error of a dynamical spectral estimation method called “the variational approach to conformational dynamics" (VAC). We bound the approximation error and estimation error for VAC estimates. Our analysis establishes VAC's convergence properties and suggests new strategies for tuning VAC to improve accuracy.

97 MATHEMATICS AND COMPUTING↗

Estimating time in quantum chaotic systems and black holes

We characterize new universal features of the dynamics of chaotic quantum many-body systems, by considering a hypothetical task of "time estimation". Most macroscopic observables in a chaotic system equilibrate to nearly constant late-time values. Intuitively, it should become increasingly difficult to estimate the precise value of time by making measurements on the state. We use a quantity called the Fisher information from quantum metrology to quantify the minimum uncertainty in estimating time. Due to unitarity, the uncertainty in the time estimate does not grow with time if we have access to optimal measurements on the full system. Restricting the measurements to act on a small subsystem or to have low computational complexity leads to results expected from equilibration, where the time uncertainty becomes large at late times. With optimal measurements on a subsystem larger than half of the system, we regain the ability to estimate the time very precisely, even at late times. Hawking's calculation for the reduced density matrix of the black hole radiation in semiclassical gravity contradicts our general predictions for unitary quantum chaotic systems. Hawking's state always has a large uncertainty for attempts to estimate the time using the radiation, whereas our general results imply that the uncertainty should become small after the Page time. This gives a new version of the black hole information loss paradox in terms of the time estimation task. By restricting to simple measurements on the radiation, the time uncertainty becomes large. This indicates from a new perspective that the observations of computationally bounded agents are consistent with the semiclassical effective description of gravity.

Black holes↗

Improving Regional Air Quality Forecasting Through Chemical Data Assimilation and Dynamic Emissions Adjustment

Poor air quality (AQ) is one of the most important human-health and environmental problems facing the United States (US). In addition to the detrimental impacts on human- and environmental-health, poor AQ has an economic cost of ~5% of the US gross domestic product (~$790 billion). AQ managers use AQ analyses and modeling to better understand, anticipate, and avoid poor AQ events. Our research focuses on improving AQ analysis/forecast skill, predictability, and emission estimates through improved and more efficient: (i) modeling and data assimilation strategies; (ii) dynamic emissions adjustment strategies; and (iii) use of satellite remote-sensing Earth observations (e.g., MOPITT, IASI, MODIS, OMI, TROPOMI, TEMPO, etc.). This seminar will review: (i) regional chemical weather forecasting/data assimilation with dynamic emissions adjustment with WRF-Chem/DART; (ii) strategies for efficiently assimilating satellite retrieval profiles with ‘compact phase space retrievals’ (CPSRs); (iii) results from joint assimilation of multiple satellite retrievals at medium (12 km × 12 km) and high (4 km × 4 km) spatial resolutions; and (iv) results from observing system simulation experiments (OSSEs) to investigate whether we can recover COVID-period anthropogenic emissions by assimilating synthetic TEMPO NO2 tropospheric column retrievals with dynamic emissions adjustment. Biographical Sketch: Dr. Mizzi is a Senior Research Fellow working and Dr. Johnson at the NASA Ames Research Center. He holds BA and MS degrees in Environmental Science from the University of Virginia, MS and PhD degrees in Applied Mathematics from the University of Colorado at Boulder (CUB), and a JD degree (with an emphasis in Environmental Law) from the University of Colorado School of Law. He worked at the National Center for Atmospheric Research for nearly 25 years on global atmospheric modeling, dynamic and physical initialization, regional hybrid data assimilation, and most recently on regional chemical data assimilation. He also worked as an environmental attorney and consultant for nearly 15 years. He is an expert in numerical modeling and is recognized internationally as a leading expert in regional, chemical data assimilation with dynamic emissions adjustment. Dr. Mizzi became affiliated with NASA Ames in March 2020 to work on improving AQ analysis/forecast skill, predictability, and ‘top-down’ emissions adjustment though the assimilation of Earth observations. An emphasis of his current work is developing methods for assimilating synthetic TEMPO retrievals to quantify the expected benefits of TEMPO relative to existing AQ observations.

Arthur P. Mizzi↗

Selection of Slow and Fast Variables in the Three-Dimensional Flight Dynamics

A transformation from the altitude and velocity state variables of three-dimensional flight mechanics to a new set of more desirable variables is found. The new variables provide a greater time-scale separation, decrease system coupling, and give better estimates of the fast-variable values along the reduced solution. One of the new variables is the often-used specific energy, whereas the other variable changes along a given trajectory, depending on the nature of the local reduced solution. Numerical examples are included.

Ardema, M.↗

Earth Observation System Flight Dynamics System Covariance Realism

This presentation applies a covariance realism technique to the National Aeronautics and Space Administration (NASA) Earth Observation System (EOS) Aqua and Aura spacecraft based on inferential statistics. The technique consists of three parts: collection calculation of definitive state estimates through orbit determination, calculation of covariance realism test statistics at each covariance propagation point, and proper assessment of those test statistics.

A-Train↗

Model-free control for frequency response support in microgrids utilizing wind turbines

Islanded microgrids integrated with renewable energy resources bring additional challenges to the primary frequency control due to their power-electronic interfaced generations. Developing an efficient control strategy for inertia emulation and frequency response support of microgrids is becoming an important issue. Wind turbine dynamics, however, will degrade the emulated inertial and primary responses with most traditional control strategies. On the other hand, compensation methods often require the knowledge of system models at the design stage and state estimation for feedback once deployed online. To overcome these challenges, this paper explores a model-free control (MFC) strategy to emulate the desired inertia and support the frequency response of a diesel-wind microgrid system. Additionally, the proposed MFC strategy utilizes a model reference concept to obtain guaranteed inertial response and improve the frequency response of microgrids. This work considers different control strategies with different operating points. The proposed controller is implemented and verified using the modified 33-bus full nonlinear model of the three-phase diesel-wind system in Simulink. Simulation results demonstrate the improved frequency response by precisely emulating the desired inertia using the MFC strategy.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Path Properties of Atmospheric Transitions: Illustration with a Low-Order Sudden Stratospheric Warming Model

Many rare weather events, including hurricanes, droughts, and floods, dramatically impact human life. To accurately forecast these events and characterize their climatology requires specialized mathematical techniques to fully leverage the limited data that are available. Here we describe transition path theory (TPT), a framework originally developed for molecular simulation, and argue that it is a useful paradigm for developing mechanistic understanding of rare climate events. TPT provides a method to calculate statistical properties of the paths into the event. As an initial demonstration of the utility of TPT, we analyze a low-order model of sudden stratospheric warming (SSW), a dramatic disturbance to the polar vortex that can induce extreme cold spells at the surface in the midlatitudes. SSW events pose a major challenge for seasonal weather prediction because of their rapid, complex onset and development. Climate models struggle to capture the long-term statistics of SSW, owing to their diversity and intermittent nature. We use a stochastically forced Holton–Mass-type model with two stable states, corresponding to radiative equilibrium and a vacillating SSW-like regime. In this stochastic bistable setting, from certain probabilistic forecasts TPT facilitates estimation of dominant transition pathways and return times of transitions. These “dynamical statistics” are obtained by solving partial differential equations in the model’s phase space. With future application to more complex models, TPT and its constituent quantities promise to improve the predictability of extreme weather events through both generation and principled evaluation of forecasts.

54 ENVIRONMENTAL SCIENCES↗

Short- and Intermediate-Range Structure and Dynamics of Fe-Ni-C Liquid Under Compression

Properties of liquid Fe alloys under high-pressure conditions are crucial for understanding the composition, thermal state, and dynamics of Earth’s core. Experiments on such liquids, however, are often performed under pressures far below those of the outer core, necessitating long extrapolations of experimental results to core conditions. Such estimates can be complicated by light elements possibly forming pressure-dependent molecular clusters that can significantly affect the physical properties of liquids as core conditions are approached. First-principles molecular dynamics simulations were employed to compute the properties of an Fe-Ni-C liquid with a composition of Fe 3.7 Ni 0.37 C at 1673 K and pressures from 0 to 67 GPa to benchmark computational methods on pressure effects on the structure and properties of the liquid relative to low pressure experimental results. The short-range structure is manifested by the coordination number (CN) of Fe/Ni-Fe/Ni being around 12, indicative of a nearly close-packed structure in the pressure range, and the CN of C-Fe/Ni gradual increasing from 6.5 to 8.5, indicative of an approximately octahedral to cubic transition as pressure increases. The Fe/Ni-Fe/Ni bond distance, however, is found to be 10 times more compressible than the C-Fe/Ni distance. The intermediate-range structure of Fe/Ni-Fe/Ni and C-Fe/Ni subsystems, described by a partial configurationally decomposed distribution function, undergoes substantial changes, characterized by a significant increase of the number of polyhedra that share 3 atoms with each other. Such dense configurations are related to an increased bulk modulus, decreased diffusion coefficient, decreased activation volume for diffusion, and increased shear viscosity. Reproducing the experimental observations at low pressures provides important support for modeling the liquid under the conditions relevant to the outer core.

58 GEOSCIENCES↗

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.↗

Situation awareness and dynamic ensemble forecasting of abnormal behavior in cyber-physical system

A plurality of monitoring nodes may each generate a time-series of current monitoring node values representing current operation of a cyber-physical system. A feature-based forecasting framework may receive the time-series of and generate a set of current feature vectors using feature discovery techniques. The feature behavior for each monitoring node may be characterized in the form of decision boundaries that separate normal and abnormal space based on operating data of the system. A set of ensemble state-space models may be constructed to represent feature evolution in the time-domain, wherein the forecasted outputs from the set of ensemble state-space models comprise anticipated time evolution of features. The framework may then obtain an overall features forecast through dynamic ensemble averaging and compare the overall features forecast to a threshold to generate an estimate associated with at least one feature vector crossing an associated decision boundary.

Abbaszadeh, Masoud↗

Two self-test methods applied to an inertial system problem

The paper describes two self-test procedures applied to the problem of estimating the biases in accelerometers and gyroscopes on an inertial platform. The first technique is the weighted sum-squared residual (WSSR) test, with which accelerator bias jumps are easily isolated, but gyro bias jumps are difficult to isolate. The WSSR method does not take full advantage of the knowledge of system dynamics. The other technique is a multiple hypothesis method developed by Buxbaum and Haddad (1969). It has the advantage of directly providing jump isolation information, but suffers from computational problems. It might be possible to use the WSSR to detect state jumps and then switch to the BH system for jump isolation and estimate compensation.

Willsky, A. S.↗