Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “State Variable”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Application of FARM to an IES scenario within the FORCE ecosystem

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM supports the HERON software module in the evaluation of the optimal dispatch by evaluating feasible set-points for the different IES unit components. Set-points are required to satisfy limits on both production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). This problem is addressed by adopting a two-stage approach. First, the HERON power dispatcher determines set-points that meet the constraints on the former variables (e.g., power levels and power ramp rate limits). These constraints are called explicit constraints. Then, FARM adjusts these set-points to ensure the respect of the limits on the latter variables given the knowledge of the system physics 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). The original version of the FARM software module (FARM-Alpha) was released by Argonne National Laboratory in January 2021. In the latest version of the code released in July 2022 (FARM-Delta), the Reference Governor (RG) algorithm was upgraded to a Multi-Input Multi-Output version from its original Single-Input Single-Output form. The RG algorithm acts to enforce constraints. With this improvement an IES unit is now treated as a single dynamic system from the standpoint of control. The crosstalk among components in an IES unit is now fully considered thereby ensuring a true optimization is obtained for those units that have multiple set-points. In this report, the capabilities of FARM-Delta operating within the FORCE ecosystem are demonstrated for an IES test case. The specific configuration of IES unit for this case was selected by the IES team with consultation from the Advanced Reactor IES Expert Group. A full TEA analysis that invoked HERON, HYBRID, FARM, and RAVEN was performed and serves to demonstrate how the latest modification to FARM algorithms (i.e., state variable selection, state-space matrices derivation, set-point verification) can shape setpoints that might otherwise compromise the health of equipment through accelerated wear and tear. In this specific test case, it was demonstrated that these algorithms ensure a more efficient utilization of steam resources to be shared by two different subsystems, namely Balance of Plant (BOP) and High-Temperature Steam Electrolysis (HTSE). Finally, some code improvements that can further enhance the user-friendliness are suggested.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Quantum parity Hall effect in Bernal-stacked trilayer graphene

The quantum Hall effect has recently been generalized from transport of conserved charges to include transport of other approximately conserved-state variables, including spin and valley, via spin- or valley-polarized boundary states with different chiralities. In this paper, we report a class of quantum Hall effect in Bernal- or ABA-stacked trilayer graphene (TLG), the quantum parity Hall (QPH) effect, in which boundary channels are distinguished by even or odd parity under the system’s mirror reflection symmetry. At the charge neutrality point, the longitudinal conductance σ x x is first quantized to 4 e 2 / h at a small perpendicular magnetic field B ⊥ , establishing the presence of four edge channels. As B ⊥ increases, σ x x first decreases to 2 e 2 / h , indicating spin-polarized counterpropagating edge states, and then, to approximately zero. These behaviors arise from level crossings between even- and odd-parity bulk Landau levels driven by exchange interactions with the underlying Fermi sea, which favor an ordinary insulator ground state in the strong B ⊥ limit and a spin-polarized state at intermediate fields. The transitions between spin-polarized and -unpolarized states can be tuned by varying Zeeman energy. Our findings demonstrate a topological phase that is protected by a gate-controllable symmetry and sensitive to Coulomb interactions.

2D materials↗
DOE OSTI2015

Effect of point defects on the thermal conductivity of UO 2 : molecular dynamics simulations

The thermal conductivity of uranium dioxide (UO 2 ) fuel is an important materials property that affects fuel performance since it is a key parameter determining the temperature distribution in the fuel, thus governing, e.g., dimensional changes due to thermal expansion, fission gas release rates, etc. The thermal conductivity of UO 2 nuclear fuel is also affected by fission gas, fission products, defects, and microstructural features such as grain boundaries. Here, molecular dynamics (MD) simulations are carried out to determine quantitatively, the effect of irradiation induced point defects on the thermal conductivity of UO 2 , as a function of defect concentrations, for a range of temperatures, 300 – 1500 K. The results will be used to develop enhanced continuum thermal conductivity models for MARMOT and BISON by INL. These models express the thermal conductivity as a function of microstructure state-variables, thus enabling thermal conductivity models with closer connection to the physical state of the fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

The thermal conductivity of mixed fuel U x Pu 1-x O 2 : molecular dynamics simulations

Mixed oxides (MOX), in the context of nuclear fuels, are a mixture of the oxides of heavy actinide elements such as uranium, plutonium and thorium. The interest in the UO 2 -PuO 2 system arises from the fact that these oxides are used both in fast breeder reactors (FBRs) as well as in pressurized water reactors (PWRs). The thermal conductivity of UO 2 fuel is an important material property that affects fuel performance since it is the key parameter determining the temperature distribution in the fuel, thus governing, e.g., dimensional changes due to thermal expansion, fission gas release rates, etc. For this reason it is important to understand the thermal conductivity of MOX fuel and how it differs from UO 2 . Here, molecular dynamics (MD) simulations are carried out to determine quantitatively, the effect of mixing on the thermal conductivity of U x Pu 1-x O 2 , as a function of PuO 2 concentrations, for a range of temperatures, 300 – 1500 K. The results will be used to develop enhanced continuum thermal conductivity models for MARMOT and BISON by INL. These models express the thermal conductivity as a function of microstructure state-variables, thus enabling thermal conductivity models with closer connection to the physical state of the fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Statistical Analysis of the Limiting Dynamics of Two dimensional Boussinesq Turbulent Systems

In the current study, we investigate rotational and stratified turbulent systems within the framework of the Boussinesq equations. The effects of rotation and stratification give rise to a variety of dynamical regimes in atmospheric and oceanic turbulence. To develop a better understanding of the fundamental flow properties through a more straightforward mathematical framework, we commence our research with a simplified two-dimensional Boussinesq model that incorporates rotation and stratification effects. We explore two distinct limiting dynamics, considering pure rotation and pure stratification separately. We examine these intriguing dynamics through numerical investigations, applying the exact solution theory of Boussinesq equations and examining small-scale perturbations to the exact solution under rotation and stratification. We conduct an analysis of statistical quantities associated with the turbulent state variables. The results indicate that, under weak rotation at the equilibrium statistical state, the flow field exhibits vortex flow. However, strong rotation transforms the flow into vertical shear flow. We observe both downscale and upscale energy transfers with a decay rate of k −3 in the absence of external forcing and dissipation. Furthermore, we delve into the study of internal gravity wave mode interactions in the simplified two-dimensional Boussinesq system, which is a crucial aspect of geophysical turbulence due to gravity. Our findings reveal that the time series of the mode coefficients exhibit wave-like interactions between wave modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

An improved hindcast approach for evaluation and diagnosis of physical processes in global climate models

Here, we present an improved procedure of generating initial conditions (ICs) for climate model hindcast experiments with specified sea surface temperature and sea ice. The motivation is to minimize errors in the ICs and lead to a better evaluation of atmospheric parameterizations' performance in the hindcast mode. We apply state variables (horizontal velocities, temperature and specific humidity) from the operational analysis/reanalysis for the atmospheric initial states. Without a data assimilation system, we apply a two-step process to obtain other necessary variables to initialize both the atmospheric (e.g., aerosols and clouds) and land models (e.g., soil moisture). First, we nudge only the model horizontal velocities towards operational analysis/reanalysis values, given a 6-hour relaxation time scale, to obtain all necessary variables. Compared to the original strategy in which horizontal velocities, temperature and specific humidity are nudged, the revised approach produces a better representation of initial aerosols and cloud fields which are more consistent and closer to observations and model's preferred climatology. Second, we obtain land ICs from an offline land model simulation forced with observed precipitation, winds, and surface fluxes. This approach produces more realistic soil moisture in the land ICs. With this refined procedure, the simulated precipitation, clouds, radiation, and surface air temperature over land are improved in the Day 2 mean hindcasts. Following this procedure, we propose a “Core” integration suite which provides an easily repeatable test allowing model developers to rapidly assess the impacts of various parameterization changes on the fidelity of modelled cloud-associated processes relative to observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Degenerate resistive switching and ultrahigh density storage in resistive memory

We show that in tantalum oxide resistive memories, activation power provides a multi-level variable for information storage that can be set and read separately from the resistance. These two state variables (resistance and activation power) can be precisely controlled in two steps: (1) the possible activation power states are selected by partially reducing resistance, then (2) a subsequent partial increase in resistance specifies the resistance state and the final activation power state. We show that these states can be precisely written and read electrically, making this approach potentially amenable for ultra-high density memories. We provide a theoretical explanation for information storage and retrieval from activation power and experimentally demonstrate information storage in a third dimension related to the change in activation power with resistance.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Entanglement structures in quantum field theories. II. Distortions of vacuum correlations through the lens of local observers

When observing a quantum field via detectors with access to only the mixed states of spatially separated, local regions—a ubiquitous experimental design—the capacity to access the full extent of distributed entanglement can be limited, shrouded by classical correlations. By performing projective measurements of the field external to two detection patches and classically communicating the results, underlying pure states may be identified for which entanglement quantification is clear. In the Gaussian continuous-variable states of the free scalar field vacuum, this protocol uncovers a disparity between the spacelike entanglement established within the field and that which is locally detectable. Further, this discrepancy is found to grow exponentially with the separation between observation regions. The protocol developed herein offers insight and practical guidance for clarifying the unavoidable distortion of quantum field correlations when viewed from the vantage of a pair of local observers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Enhanced Tensor Completion Based Approaches for State Estimation in Distribution Systems

Grid state estimation is essential for effective control and management of distribution systems. While weighted least squares has been the conventional method for state estimation, sparsity-aware methods have become popular due to their superior performance with limited data. Matrix completion and compressed sensing-based state estimation approaches exploit the underlying smoothness in the state variables. However, classic matrix completion methods do not take into account the temporal correlation of system states. Compressed sensing methods, on the other hand, require an appropriate choice of sparsifying basis that may not be easy to identify. This paper proposes a blocktensor completion based framework which uses an alternative approach to estimate voltage phasor, power injections and branch currents. This approach utilizes the temporal correlation of the system states in a tensor trace-norm minimization formulation with power flow equations as constraints. Herein, feature scaling is introduced in the problem formulation to benefit from the improved sensitivity of the tensor trace norm to the matrix columns in the scaled unfoldings of the tensor. Weighted tensor norm is utilized to exploit the structures of the different unfoldings of the state measurement tensor to improve the voltage estimation. The estimation accuracy is further improved by alternatively estimating the tensor columns and increasing the available data at each stage in the tensor completion process. The proposed methods are evaluated on the IEEE-33, 37 test systems and a 100- node test system. The proposed methods are shown to provide significant performance gains relative to the classic matrix and tensor completion based approaches.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Schrödinger cat states of a nuclear spin qudit in silicon

High-dimensional quantum systems are a valuable resource for quantum information processing. They can be used to encode error-correctable logical qubits, which has been demonstrated using continuous-variable states in microwave cavities or the motional modes of trapped ions. For example, high-dimensional systems can be used to realize ‘Schrödinger cat’ states, which are superpositions of widely displaced coherent states that can be used to illustrate quantum effects at large scales. Recent proposals have suggested encoding qubits in high-spin atomic nuclei, which are finite-dimensional systems that can host hardware-efficient versions of continuous-variable codes. Here, in this study, we demonstrate the creation and manipulation of Schrödinger cat states using the spin-7/2 nucleus of an antimony atom embedded in a silicon nanoelectronic device. We use a multi-frequency control scheme to produce spin rotations that preserve the symmetry of the qudit, and we constitute logical Pauli operations for qubits encoded in the Schrödinger cat states. Our work demonstrates the ability to prepare and control non-classical resource states, which is a prerequisite for applications in quantum information processing and quantum error correction, using our scalable, manufacturable semiconductor platform.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effective permeability of an immiscible fluid in porous media determined from its geometric state

Based on the phenomenological extension of Darcy's law, two-fluid flow is dependent on a relative permeability function of saturation only that is process or path dependent with an underlying dependency on pore structure and wettability. For applications (fuel cells to underground CO 2 storage), it is imperative to determine the effective phase permeability relationships where the traditional approach is based on the inverse modeling of time-consuming experiments. The underlying reason is that the fundamental upscaling step from pore to Darcy scale, which links the pore structure of the porous medium to the continuum hydraulic conductivities, is not solved. Herein, we develop an artificial neural network (ANN) that relies on fundamental geometrical relationships to determine the mechanical energy dissipation during creeping immiscible two-fluid flow. Finally, the developed ANN is based on a prescribed set of state variables based on physical insights that predict the effective permeability of 4500 unseen pore-scale geometrical states with R 2 = 0.98.

42 ENGINEERING↗
DOE OSTI2025

Karhunen–Loève deep learning method for surrogate modeling and approximate Bayesian parameter estimation

We evaluate the performance of the Karhunen-Loève Deep Neural Network (KL-DNN) framework for surrogate modeling and approximate Bayesian parameter estimation in partial differential equation models. In the surrogate model, the Karhunen-Loève (KL) expansions are used for the dimensionality reduction of the number of unknown parameters and variables, and a deep neural network is employed to relate the reduced space of parameters to that of the state variables. The KL-DNN surrogate model is used to formulate a maximum-a-posteriori-like least-squares problem, which is randomized to draw samples of the posterior distribution of the parameters. We test the proposed framework for a hypothetical unconfined aquifer via comparison with the forward MODFLOW and inverse PEST++ iterative ensemble smoother (IES) solutions as well as the state-of-the-art Fourier neural operator (FNO) and deep operator networks (DeepONets) operator learning surrogate models. Our results show that the KL-DNN surrogate model outperforms FNO and DeepONet for forward predictions. For solving inverse problems, the randomized algorithm provides the same or more accurate Bayesian predictions of the parameters than IES as evidenced by the higher log-predictive probability of both the estimated parameter field and the forecast hydraulic head. The posterior mean obtained from the randomized algorithm is closer to the reference parameter field than that obtained with FNO as the maximum a posteriori estimate.

Approximate Bayesian inference↗
DOE OSTI2022

Millimeter Light Curves of Sagittarius A* Observed during the 2017 Event Horizon Telescope Campaign

The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5–11 in the 1.3 mm wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very long baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare and exhibit strongly enhanced variability. All of the light curves are consistent with a red-noise process, with a power spectral density (PSD) slope measured to be between −2 and −3 on timescales between 1 minute and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 hr. The spectral energy distribution is flat at 220 GHz, and there are no time lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*’s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Dependence of calculated postshock thermodynamic variables on vibrational equilibrium and input uncertainty

The purpose of this article is to explore the dependence of calculated postshock thermodynamic properties in shock tube experiments upon the vibrational state of the test gas and upon the uncertainties inherent to calculation inputs. This paper first offers a comparison between state variables calculated according to a Rankine–Hugoniot–equation-based algorithm, known as FROSH, and those derived from shock tube experiments on vibrationally nonequilibrated gases. It is shown that incorrect vibrational relaxation assumptions could lead to errors in temperature as large as 8% for 25% oxygen/argon mixtures at 3500 K. Following this demonstration, this article employs the algorithm to show the importance of correct vibrational equilibration assumptions, noting, for instance, that errors in temperature of up to about 2% at 3500 K may be generated for 10% nitrogen/argon mixtures if vibrational relaxation is not treated properly. Lastly, this article presents an extensive uncertainty analysis, showing that postshock temperatures can be calculated with root-of-sum-of-square errors of better than ±1% given sufficiently accurate experimentally measured input parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Interacting mesons as degrees of freedom in a chiral model

Here, we study the equation of state of hot and dense hadronic matter using an extended chiral mean field (CMF) model framework where the addition is the inclusion of interactions of thermally excited mesons. This is implemented by calculating the in-medium masses of pseudoscalar and vector mesons, obtained through the explicit chiral symmetry-breaking and vector-interaction terms in the Lagrangian, respectively, prior to applying the mean-field approximation. As a result, the in-medium meson contributions generate a feedback term to the CMF’s equations of motion, which then modifies the equation of state. With this improvement, we quantify the effect on the equation of state of strongly interacting matter through comparisons with state-of-the-art lattice QCD results and other hadronic models like the hadron resonance gas model. We find that the results of the updated hadronic CMF model with an improved meson description (mCMF) provide a better agreement with lattice-QCD data for thermodynamic state variables across a wide range of temperatures and baryon chemical potentials.

phenomenology↗
DOE OSTI2023

Variable gap thermal conductivity apparatus and method

An apparatus and a method for determining the thermal conductivity of a fluid specimen are provided. The apparatus and the method include determining thermal conductivity using a quasi-steady state variable gap axial flow technique. The fluid specimen is heated on one side by a heat source with a known power output and cooled on the other side. After reaching steady state, a resulting temperature drop through the fluid specimen exists. This temperature drop, the known fluid specimen thickness (or gap distance), and the known power output are used to calculate the thermal resistance of the fluid specimen. The thermal conductivity of the fluid specimen is then determined using a curve fit of thermal resistance with respect to gap distance.

Gallagher, Ryan C.↗
DOE OSTI2022

Autoencoder Neural Network for chemically reacting systems

Incorporating detailed chemical kinetic models is critical for accurate simulations of reacting flows. However, detailed models involve a large number of thermochemical (TC) state variables. Solving the governing equations to evolve these TC variables becomes impractical for real-world applications. In this work, we propose an autoencoder (AE) neural network (NN)-based reduced model to accelerate such simulations. The AE NN is first trained to find a low-dimensional latent representation of the TC states. Then, the evolving state of a chemical system can be tracked by solving the equations of the latent variables instead of the original TC equations. We demonstrate the reduced model in a syngas CO/H 2 combustion system, using training data collected from canonical perfectly stirred reactors (PSRs). It is found that the AE model can reduce the dimension of the combustion system from 12 to 2 while maintaining low reconstruction error and excellent elemental mass conservation for the test dataset. In the a posteriori test, the combustion states obtained from solving the two latent equations are compared to those from solving the 12 equations of the full model. The AE reduced method is found to be able to capture the diverse combustion states on the top two branches of the S-curve well including the extinction turning point, but with higher prediction errors for states near the ignition turning point.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Data Driven Temperature Control of Heat Exchangers

The following poster presents a couple of approaches for temperature control at heat exchanger’s exit. One of the approaches is based on using n-step ahead neural network while the other approach utilizes a combination of control + learning module for temperature control. Designing robust controllers for typical heat exchangers is very challenging primarily because of complexities associated with the synthesis of dynamics of industrial heat exchangers. Even if a dynamics model is synthesized, the accuracy of it would be called into question as there can be lot of unmodeled system dynamics and random noise in input parameters. In many situations, the degree of uncertainty in the model of system being controlled limits the utility of optimal control design. Controller can be manually tuned on the field but it is very difficult to determine manual adjustments that result in overall improvement. These practical challenges call for learning based approaches for temperature control that would result in optimal heat exchanger performance. We utilize the existing mathematical model for air-water system synthesized by Underwood and Crawford for testing our algorithms. Underwood and Crawford developed a model of heating coil by fitting a second-order non-linear equations to measurements of air and water temperature and flow rates obtained from the actual coil. The uncertainties in the model include air and water inlet temperatures and air-flow rate. These uncertainties are modeled as random variables that are sampled from a uniform distribution. The control variable is the water flow rate which is controlled for regulating the air-outlet temperature. Neural networks are frequently used in various literatures for learning time-series. We propose a modified version of neural networks known as n-step ahead neural network that take the state variables and output at time t as input to successfully predict output at time t+n. The variable n is a hyper-parameter which can be optimized based on performance constraints. The n-step ahead neural network was then trained on training data generated from experimental model to high training accuracies. This trained model performed appreciably in its objective to regulate outlet temperature of a heat exchanger. The second approach uses a combination of control and learning module for temperature regulation. The control module used is a proportional controller while the learning module used is steady state neural network. The steady state neural network is used to learn the controller output as a function of random variables and temperature set-point. The proportional controller is then used to compensate for transient behavior of linearized system around various operating points. The algorithm works effectively for temperature regulation as neural network captures the non-linear input-output behavior of closed loop system leaving proportional controller to compensate for steady state errors of a linearized plant. The proposed algorithm’s performances were evaluated against the performance of a PI controller. Further, the robustness of the algorithms were also studied in great detail under small perturbations around tuned parameters of the proposed algorithms.

Sundar, Sreenath↗