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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Diffusion coefficients predicting facet-dependent crystallization in amorphous silicon nitride

Amorphous silicon nitride is a common material in microelectronics devices, which acts as an insulating barrier. Extended annealing times at elevated temperature can initiate crystallization of α-Si 3 ⁢ N 4 , which does not possess the same barrier properties. Molecular dynamics can resolve the fundamental mechanism for α-Si 3 ⁢N 4 crystallization and the influence of local environments. Here, we compare two interatomic potentials and conclude that these models predict structural features (e.g., angular distributions and densities) which span the range of experimental measurements. We confirmed these models reproduce experimental estimates of activation energy and leveraged these models to identify crystallization drivers. We conclude that near-T g , facet-dependent silicon nitride crystal growth rates can be predicted directly by either bulk or interfacial diffusion properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

CFDverify

Estimating the discretization error of computational fluid dynamics (CFD) or other scientific codes as part of solution verification is often a non-trivial part of the analysis process. Methods can be complicated and may include assumptions/qualifications that need to be checked during analysis. CFD analysts, therefore, are likely to make errors when trying to conduct this necessary analysis on their own in not knowing about the best method for their problem, not correctly implementing a method, or in not have diagnostic tools to determine if the method was correctly applied.

Weinmeister, Justin [Oak Ridge National Laboratory↗
DOE OSTI2017

Imputation for multisource data with comparison and assessment techniques

Missing data are prevalent issue in analyses involving data collection. The problem of missing data is exacerbated for multisource analysis, where data from multiple sensors are combined to arrive at a single conclusion. In this scenario, it is more likely to occur and can lead to discarding a large amount of data collected; however, the information from observed sensors can be leveraged to estimate those values not observed. We propose two methods for imputation of multisource data, both of which take advantage of potential correlation between data from different sensors, through ridge regression and a state-space model. These methods, as well as the common median imputation, are applied to data collected from a variety of sensors monitoring an experimental facility. Performance of imputation methods is compared with the mean absolute deviation; however, rather than using this metric to solely rank themethods,we also propose an approach to identify significant differences. Imputation techniqueswill also be assessed by their ability to produce appropriate confidence intervals, through coverage and length, around the imputed values. Finally, performance of imputed datasets is compared with a marginalized dataset through a weighted k-means clustering. In general, we found that imputation through a dynamic linearmodel tended to be the most accurate and to produce the most precise confidence intervals, and that imputing the missing values and down weighting them with respect to observed values in the analysis led to the most accurate performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Global Data-Driven Determination of Baryon Transition Form Factors

Hadronic resonances emerge from strong interactions encoding the dynamics of quarks and gluons. The structure of these resonances can be probed by virtual photons parametrized in transition form factors. Here, in this study, twelve N* and Δ transition form factors at the pole are extracted from data with the center-of-mass energy from πN threshold to 1.8 GeV, and the photon virtuality 0 ≤ Q 2 /GeV 2 ≤ 8. For the first time, these results are determined from a simultaneous analysis of more than one state, i.e., ~10 5 π⁢N, η⁢N, and K⁢Λ electroproduction data. In addition, about 5×10 4 data in the hadronic sector as well as photoproduction serve as boundary conditions. For the Δ⁡(1232) and N⁡(1440) states our results are in qualitative agreement with previous studies, while the transition form factors at the poles of some higher excited states are estimated for the first time. Realistic uncertainties are determined by further exploring the parameter space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Dynamics of Qubits in the Field of Unipolar Pulses: Magnus Propagator, Generalized Area Theorem, and Motion on Groups

The problem of speeding up quantum computations through controlling the states of qubits by short large-amplitude unipolar pulses is discussed. A new method is proposed to describe the dynamics of qubits that is based on the Magnus representation for the evolution operator (propagator) of qubits, which allows one to find a solution to the quantum control problem outside the framework of perturbation theory. The evolution of the states of a multiqubit system is represented as a motion on a special unitary group SU(N), similar to the rotation of the Bloch vector on the group SU(2) in the case of a single qubit. As an example, unitary transformations of single-qubit states under different excitation modes are considered. A general formula is found that establishes a relation between the shape of a pulse and the result of its action, similar to the well-known “area theorem” for Rabi pulses. The effect of unipolar pulses on a two-qubit system is analyzed in detail. Two situations are considered: first, a symmetric configuration, when the parameters of qubits do not differ and are acted upon by the same magnetic field from a fluxon, and second, when the configuration of levels arises in which, for one of the qubits, a transition to the nearest level becomes possible under perturbation. The general case of two interacting qubits—a four-level system—is considered analytically. A direct numerical simulation of the dynamics of a multiqubit system is performed to estimate the accuracy of approximation, and fidelity is used as a criterion of accuracy of operations. It is shown that the expressions obtained for the propagator allow one to formulate conditions for the pulse parameters that are necessary to perform logical operations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Need A Boost? A Comparison of Traditional Commuting Models with the XGBoost Model for Predicting Commuting Flows

Commuting models estimate the number of commuting trips from home to work locations in a given area. Since their infancy, they have been increasingly used in a variety of fields to reduce traffic and pollution, drive infrastructure choices, and solve a variety of other problems. Traditional commuting models, such as gravity and radiation models, typically have a strict structural form and limited number of input variables, which may limit their ability to predict commuting flows as well as machine learning models that might better capture the complex dynamics of the commuting process. To determine whether machine learning models might add value to the field of commuter flow prediction, we compare and discuss the performance of two standard traditional models with the XGBoost machine learning algorithm for predicting home to work commuter flows from a well-known United States commuting dataset. We find that the XGBoost model outperforms the traditional models on three commonly used metrics, indicating that machine learning models may add value to the field of commuter flow prediction.

Morton, April↗
DOE OSTI2022

Development and Performance Evaluation of of Active Insulation Systems using Solid-State Thermal Switches

Traditional building envelopes have passive insulation systems that cannot respond to dynamic changes in the environment. An Active Insulation System (AIS) consists of Active Insulation Materials (AIMs), which dynamically vary the thermal conductivity of the insulation system. Several researchers have evaluated the impact of AIS on building thermal and energy performance by using simulation tools. Up to 70% savings in annual heating and cooling energy and significant reductions in peak demand have been predicted for some climates with wall systems employing AIS. Sensitivity studies have provided requirements for the AIS properties and control schemes that can lead to those savings. However, materials and assembly development have not yet achieved a cost-effective product that achieves the required performance. In this study, we present the process to develop an AIS that we will install in a test hut for their performance evaluation. Minimum performance criteria of the AIS system are developed based on Rmin/Rmax ratio, required time and efficiency to switch states and cost estimate. The next step carried out during this study is creating the concept to meet the requirements, predicting the performance by simulations, developing the experimental setup for bench-scale testing, and finally, constructing a full-scale wall assembly and monitoring the performance when exposed to natural weather conditions. The selected approach uses off-the-shelf products to create a materials system that can switch R-value between R~=1 ft2·°F·h/BTU (0.18 m2·K/W) and R~=7 ft2·°F·h/BTU (1.23 m2·K/W and have a switching time of less than one minute between R-high and R-low.

Iffa, Emishaw↗
DOE OSTI2026

CO 2 Sorption in Moisture Swing Anion Exchange Resins for Direct Air Capture: Experimental Isotherm Determination and Modeling

Moisture swing (MS) sorption is a promising direct air capture (DAC) technology to achieve negative CO 2 emissions and counter global warming. Here, in this study, MS CO 2 sorption in a model MS sorbent, IRA900, was investigated. IRA900 is a macroporous commercial strong-base anion exchange resin (AER) with quaternary ammonium functional groups. A rigorous CO 2 desorption process was developed to desorb all CO 2 from the sample and obtain CO 2 sorption isotherms as a function of relative humidity, CO 2 partial pressure, and temperature. CO 2 sorption increased with decreasing relative humidity, decreasing temperature, and increasing CO 2 partial pressure. Remarkably, the total CO 2 uptake from the sorption isotherms matched the ion exchange capacity (IEC) of the material, suggesting a stoichiometry of one CO 2 molecule reacting per active site. An isotherm model for CO 2 sorption starting with empty sorption sites (i.e., with the AER in the OH – form) described the experimental data well, supporting the hypothesis that the bicarbonate-loaded AER fully unloaded to the OH – state following rigorous CO 2 desorption. CO 2 sorption kinetics were studied, and carbon diffusion coefficients were estimated as a function of %RH and CO 2 partial pressure. Overall, this study provides a roadmap for the systematic evaluation of MS CO 2 DAC materials.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

The Energy Exascale Earth System Model Version 3: 1. Overview of the Atmospheric Component

This paper describes the atmospheric component of the US Department of Energy's Energy Exascale Earth System Model (E3SM) version 3. Significant updates have been made to the atmospheric physics compared to earlier versions. Specifically, interactive gas chemistry has been implemented, along with improved representations of aerosols and dust emissions. A new stratiform cloud microphysics scheme more physically treats ice processes and aerosol‐cloud interactions. The deep convection parameterization has been largely improved with sophisticated microphysics for convective clouds, making model convection sensitive to large‐scale dynamics, and incorporating the dynamical and physical effects of organized mesoscale convection. Improvements in aerosol wet removal processes and parameter re‐tuning of key aerosol and cloud processes have improved model aerosol radiative forcing. The model's vertical resolution has increased from 72 to 80 layers with the extra eight layers added in the lower stratosphere to better simulate the Quasi‐Biennial Oscillation. These improvements have enhanced E3SM's capability to couple aerosol, chemistry, and biogeochemistry and reduced some long‐standing biases in simulating tropical variability. Compared to its predecessors, the model shows a much stronger signal for the Madden‐Julian Oscillation, Kelvin waves, mixed Rossby‐gravity waves, and eastward inertia‐gravity waves. Aerosol radiative forcing has been considerably reduced and is now better aligned with community best estimates, leading to significantly improved skill in simulating historical temperature records. Its simulated mean‐state climate is largely comparable to E3SMv2, but with some notable degradation in shortwave cloud radiative effect, precipitable water, and surface wind stress, which will be addressed in future updates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Tri-Level Linear Programming Model for Automatic Load Shedding Using Spectral Clustering

Traditional load shedding schemes can be inadequate in grids with high renewable penetration, leading to unstable events and unnecessary grid islanding. Although for both manual and automatic operating modes load shedding areas have been predefined by grid operators, they have remained fixed, and may be sub-optimal due to dynamic operating conditions. In this work, a distributed tri-level linear programming model for automatic load shedding to avoid system islanding is presented. Preventing islanding is preferred because it reduces the need for additional load shedding besides the disconnection of transmission lines between islands. This is crucial as maintaining the local generation-demand balance is necessary to preserve frequency stability. Furthermore, uneven distribution of generation resources among islands can lead to increased load shedding, causing economic and reliability challenges. This issue is further compounded in modern power systems heavily dependent on non-dispatchable resources like wind and solar. The upper-level model uses complex power flow measurements to determine the system areas to shed load depending on actual operating conditions using a spectral clustering approach. The mid-level model estimates the area system state, while the lower-level model determines the locations and load values to be shed. The solution is practical and promising for real-world applications.

Baquedano-Aguilar, Mario D.↗
DOE OSTI2020

Inferring the Dynamics of the State Evolution During Quantum Annealing

To solve an optimization problem using a commercial quantum annealer, one has to represent the problem of interest as an Ising or a quadratic unconstrained binary optimization (QUBO) problem and submit its coefficients to the annealer, which then returns a user-specified number of low-energy solutions. It would be useful to know what happens in the quantum processor during the anneal process so that one could design better algorithms or suggest improvements to the hardware. However, existing quantum annealers are not able to directly extract such information from the processor. Hence, in this work we propose to use advanced features of D-Wave 2000Q to indirectly infer information about the dynamics of the state evolution during the anneal process. Specifically, D-Wave 2000Q allows the user to customize the anneal schedule, that is, the schedule with which the anneal fraction is changed from the start to the end of the anneal. Furthermore, using this feature, we design a set of modified anneal schedules whose outputs can be used to generate information about the states of the system at user-defined time points during a standard anneal. With this process, called "slicing", we obtain approximate distributions of lowest-energy anneal solutions as the anneal time evolves. We use our technique to obtain a variety of insights into the annealer, such as the state evolution during annealing, when individual bits in an evolving solution flip during the anneal process and when they stabilize, and we introduce a technique to estimate the freeze-out point of both the system as well as of individual qubits.

42 ENGINEERING↗
DOE OSTI2024

FARM User Guidance and Instructions

The FARM (Feasible Actuator Range Modifier) module is a component of the RAVEN-based FORCE framework for the analysis of Integrated Energy Systems (IES). FARM aids HERON in the solution of the power dispatch problem by evaluating feasible set-point signals to be issued to the control systems of the different IES unit components. Set-points need to satisfy limits on both production variables (i.e., the variables to be optimized such as the electrical power, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). To enforce all these limits, a two-stage approach is adopted. First, the power dispatcher algorithm in HERON module estimates set-points that meet the constraints on the production variables, e.g., power levels and power ramp rate limits. These constraints are called explicit constraints. Then, if necessary, FARM adjusts these set-points to ensure the respect of the limits on the process variables of interest, given the knowledge of the system dynamics acquired through machine learning algorithms. These constraints are called implicit constraints. From this standpoint, FARM constitutes a bridge between the HERON power dispatcher that adopts a simplified description of the IES unit (low-resolution physics) and the HYBRID high-fidelity models (high-resolution physics). In this report, the general workflow and the software requirements of FARM module are summarized, and the detailed instructions for installing FARM software, running built-in example cases, deriving Linear Parameter-Varying (LPV) state-space models, and using FARM for user-defined power dispatch problems are provided.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

The Foundational Industrial Energy Dataset (FIED): Open-Source Data on Industrial Facilities

The state of data on industrial energy use has co-evolved over several decades with the demands of industrial energy analysis. The most recent development - analysis in support of decarbonizing the industrial sector - has changed the characteristics of industrial data that are useful for analysts and model developers. Although data and its collection processes may be cast from a conventional viewpoint as objective and free from the influence of social dynamics, this provides an incomplete picture of not only the processes by which information is generated, but also the limitations and opportunities of data to be useful for analysis. The foundational industry energy data set (FIED) is a result of the confluence of trends in open data and the demand for higher resolution industrial energy analysis. The general approach to compiling the FIED involves accessing, filtering, and formatting data published by federal organizations on the Internet for public use. Unlike most industrial energy datasets, which are published by the U.S. Energy Information Administration (EIA), the FIED relies on core datasets from the U.S. Environmental Protection Agency (EPA). The FIED addresses several of the areas of growing disconnect between the demands of industrial energy analysis and the state of industrial energy data by providing unit-level characterization - including estimates of energy use, greenhouse gas emissions, and design capacities - for facilities that are identified by latitude and longitude. This enables local-level analysis of existing combustion equipment, as well as regional comparisons with traditional industrial energy data estimates. The report summarizes the general logic behind compiling the FIED. The FIED itself and its Python code are available from OpenEI and GitHub, respectively.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid–solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Quantum magnetism in the frustrated square lattice oxyhalides YbBi 2 IO 4 and YbBi 2 ClO 4

Square-lattice systems offer a direct route for realizing two-dimensional (2D) quantum magnetism with frustration induced by competing interactions. In this work, the square lattice materials YbBi 2 IO 4 and YbBi 2 ClO 4 were investigated using a combination of magnetization and specific-heat measurements on polycrystalline samples. Specific-heat measurements provide evidence for long-range magnetic order below $T$ N = 0.21 K (0.25 K) for YbBi 2 IO 4 (YbBi 2 ClO 4 ). On the other hand, a rather broad maximum is found in the temperature-dependent magnetic susceptibility, located at $T$ max = 0.33 K (0.38 K) in YbBi 2 IO 4 (YbBi 2 ClO 4 ), consistent with the quasi-2D magnetism expected for the large separation between the magnetic layers. Estimation of the magnetic entropy supports the expected Kramers' doublet ground state for Yb 3+ and the observed paramagnetic behavior is consistent with a well-isolated doublet. Roughly two-thirds of the entropy is consumed above $T$ N , due to a combination of the quasi-2D behavior and magnetic frustration. The impact of frustration is examined from the viewpoint of a simplified $J$ 1 -$J$ 2 square lattice model, which is frustrated for antiferromagnetic interactions. Specifically, a high-temperature series expansion analysis of the temperature-dependent specific-heat and magnetization data yields $J$ 2 /$J$ 1 = 0.30 (=0.23) for YbBi 2 IO 4 (YbBi 2 ClO 4 ). In conclusion, this simplified analysis suggests strong frustration that should promote significant quantum fluctuations in these compounds, and thus motivates future work on the static and dynamic magnetic properties of these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Palladium at high pressure and high temperature: A combined experimental and theoretical study

Palladium is one of the most important technological materials, yet its phase diagram remains poorly understood. At ambient conditions, its solid phase is face-centered cubic (fcc). However, another solid phase of Pd, body-centered cubic (bcc), was very recently predicted in two independent theoretical studies to occur at high pressures and temperatures. In this work, we report an experimental study on the room-temperature equation of state (EOS) of Pd to a pressure of 80 GPa, as well as a theoretical study on the phase diagram of Pd including both fcc-Pd and bcc-Pd. Our theoretical approach consists in ab initio quantum molecular dynamics (QMD) simulations based on the Z methodology which combines both direct Z method for the simulation of melting curves and inverse Z method for the calculation of solid–solid phase transition boundaries. We obtain the melting curves of both fcc-Pd and bcc-Pd and an equation for the fcc–bcc solid– solid phase transition boundary as well as the thermal EOS of Pd which is in agreement with experimental data and QMD simulations. We uncover the presence of another solid phase of Pd on its phase diagram, namely, random hexagonal close-packed (rhcp), and estimate the location of the rhcp-bcc solid–solid phase transition boundary and the rhcp–fcc–bcc triple point. We also discuss the topological similarity of the phase diagrams of palladium and silver, the neighbor of Pd in the periodic table. We argue that Pd is a reliable standard for shock-compression studies and present the analytic model of its principal Hugoniot in a wide pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The effect of initial conditions on mixing transition of the Richtmyer–Meshkov instability

We investigate the late-time Richtmyer–Meshkov instability (RMI) growth of sinuous perturbations on an air/sulphur hexafluoride interface (Atwood number, A ~ 0.67 ) subjected to a Mach 1.2 planar shock wave at Los Alamos National Laboratory's vertical shock tube facility. Interface perturbations are established using a novel membraneless technique where cross-flowing air and SF 6 separated by an oscillating splitter plate create a perturbed density interface. The interface formed has multi-modal features and residual small perturbations, however, a dominant mode is still noticeable. The late-time perturbation growths scale with $ka_0$ initial conditions (where $k$ is the wavenumber and $a_0$ is the initial amplitude of the dominant mode) as measured at the pre-shock interface. Past nonlinear models based on potential-flow theory, heuristic/interpolation approaches, Padé approximants and numerical simulations are evaluated against present experimental results. Accounting for an explicit $ka_0$ dependence in Sadot et al.'s model, we propose an empirical rational function that captures the asymptotic behaviour of perturbation growth for a broad range of initial conditions ( 0.30 $≤ ka_0 ≤$ 0.86 ). The onset of mixing transition and its initial condition dependence are investigated with respect to the minimum state criterion ($Re =$ 1.6 $\times$ 10 5 ) for unsteady flows by Zhou. Earlier mixing transitions for higher ka0 initial conditions are noted from local and global Reynolds number estimates which are corroborated by the existence of an inertial sub-range and formation of mixing regions indicating the physical significance of the minimum state criterion in RMI flows. The transition is accompanied by the increasing teapot-like appearance of joint probability density functions of $p–q$ (invariants of the reduced velocity gradient tensor), establishing the technique as a useful tool for turbulence detection in two-dimensional diagnostics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗