Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equation of State”

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,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Power System Frequency Dynamics Modeling, State Estimation, and Control using Neural Ordinary Differential Equations (NODEs) and Soft Actor-Critic (SAC) Machine Learning Approaches

With the global energy transition of the electric power system, grid control, supervision, and protection is becoming more challenging. With the increasing integration of renewable energy sources (RES), the system dynamics are changing, causing traditional power system dynamic modeling with swing equation-based modeling approaches to fail. Additionally, the converter-dominated power grid is decreasing the system inertia, making the power system more fragile to the frequency swings. This paper first investigates and compares the application of a model-based Kalman filter state estimation approach with (i) a model-free machine learning approach --- neural ordinary differential equations (NODEs) --- and (ii) a data-driven system identification (SysId) approach to model and infer critical state values of the power system frequency dynamics. Then a model predictive control (MPC) framework is compared to a model-free Soft Actor-Critic (SAC) reinforcement learning (RL) control algorithm in providing efficient fast frequency response (FFR) to the power system frequency dynamics. The approaches are compared in terms of their performance goals as well as their per-timestep computational efficiency. Furthermore, the comparative study for state estimation shows that for the model-free requirement, both NODEs and SysId can provide accurate state estimates; however, with increasing model complexity, NODEs can be a better choice for model identification. Similarly, the results from the FFR comparative study show that the SAC RL-based FFR, once trained, outperforms MPC with better control signals and faster computation time, making the SAC RL-based FFR better option for providing FFR to the power system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The development of natural gas and hydrogen pipeline capital cost estimating equations

Natural gas pipeline cost data collected by the Oil and Gas Journal (O&GJ) for interstate pipelines constructed from 1980 through 2017 were used to develop capital cost estimating equations that are a function of pipeline diameter, length, and U.S. region. Equations were developed for material, labor, miscellaneous, and right-of-way costs, the four cost components in the O&GJ data, for six different regions of the United States (U.S.). Each equation is a function of pipeline diameter and length. Adjustment mechanisms were then developed for converting the natural gas pipeline equations into equations for estimating the costs of hydrogen pipelines. These adjustments were based in part on an analysis completed by the National Institute for Standards and Technology (NIST). Here the results of this work were used to update cost models in the Hydrogen Delivery Scenario Analysis Model (HDSAM), developed by Argonne National Laboratory for the U.S. Department of Energy’s Hydrogen Program. Our analysis shows a wide range of pipeline cost across different U.S. regions, especially with respect to labor and right-of-way costs. The developed cost formulas for hydrogen pipelines are both important and timely as hydrogen is being considered as a zero-carbon energy carrier with the potential to decarbonize all energy sectors, and the cost of hydrogen transportation is essential for techno-economic analysis of its potential use in these sectors.

03 NATURAL GAS↗
DOE OSTI2024

Parallel expansion of a fuel pellet plasmoid

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

Physics↗
DOE OSTI2021

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Demonstration of x-ray fluorescence spectroscopy as a sensitive temperature diagnostic for high-energy-density physics experiments

We present the use of x-ray fluorescence spectroscopy (XFS) to a sensitive temperature diagnostic in shocked foams at temperatures of 30–75 eV. Cobalt-doped foams were shock compressed using a planar drive at the OMEGA laser facility and photo-pumped with a Zn He⁢𝛼 x-ray source. Analysis of the resulting cobalt 𝐾⁢𝛽 x-ray fluorescence spectra using collisional radiative codes allows the temperature to be determined in the shocked foams. Furthermore, this method provides a sensitive and robust technique to determine temperatures in high-energy-density physics experiments in the tens of electronvolts temperature range. In these experiments, we find that radiation hydrodynamic simulations predict a lower temperature in the shocked foams compared to analysis of the XFS data using collisional radiative models. Although additional experiments with an independent temperature diagnostic to absolutely calibrate XFS spectra for these conditions will be required to resolve this discrepancy, these results demonstrate the excellent temperature sensitivity of XFS spectra for high-energy-density physics experiments.

Atomic spectra↗
DOE OSTI2021

Quantum dynamics of a nucleon in the Fermi accelerator

Highlights: • In nuclear physics, the Fermi accelerator is a coherent open quantum system. • There is not a true equilibrium for a nucleon‘s mean-energy in the Fermi accelerator. • Occupation of higher energy eigen-states remains even after the mean-energy plateaus. • This work paves the way for work on quantum tunnelling in the Fermi accelerator. The quantum dynamics of a particle in a one-dimensional box with an oscillating wall (the Fermi accelerator) is investigated. The model is applied to the motion of a single nucleon in the mean-field potential of a heavy atomic nucleus whose surface vibrates. By directly solving the time-dependent Schrödinger equation, both the state of the particle and its mean-energy are studied. The effects of the frequency of the wall oscillation on the nucleon’s energy are addressed. Its energy oscillates in phase with the moving wall for all frequencies, showing no chaotic behaviour. There is a large initial peak of the nucleon’s energy as the particle adjusts to the sudden change in the size of the box and a varying relaxation time as it plateaus towards lower energy and a partial equilibrium. Small oscillations in energy continue, since there cannot be a true equilibrium while the wall is moving. The quantum coherence between the different parts of the nucleon’s wave-function in real space is very much preserved. This research lays the foundation for future investigations into quantum tunnelling in the Fermi accelerator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Investigation of the ignition processes of a multi-injection flame in a Diesel engine environment using the flamelet model

In this paper, the first flamelet analysis is conducted of a highly resolved DNS of a multi-injection flame with both auto-ignition and ignition induced by flame-flame interaction. A novel method is proposed to identify the different combustion modes of ignition processes using generalized flamelet equations. The state-of-the-art DNS database generated by Rieth et al. (US National Combustion Meeting, 2019) for a multi-injection flame in a Diesel engine environment is investigated. Three-dimensional flamelets are extracted from the DNS at different time instants with a focus on auto-ignition and interaction-ignition processes. The influences of mixture field interactions and the scalar dissipation rate on the ignition process are investigated by varying the species composition boundary conditions of the transient flamelet equations. Budget analyses of the generalized flamelet equations show that the transport along the mixture fraction iso-surface is insignificant during the auto-ignition process, but becomes important when interaction-ignition occurs, which is further confirmed through a flamelet regime classification method.

42 ENGINEERING↗
DOE OSTI2022

State elimination for mixed‐integer optimal control of partial differential equations by semigroup theory

Abstract Mixed‐integer optimal control problems governed by partial differential equations (MIPDECOs) are powerful modeling tools but also challenging in terms of theory and computation. We propose a highly efficient state elimination approach for MIPDECOs that are governed by partial differential equations that have the structure of an abstract ordinary differential equation in function space. This allows us to avoid repeated calculations of the states for all time steps, and our approach is applied only once before starting the optimization. The presentation of theoretical results is complemented by numerical experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Exotic phases in finite-density ℤ 3 theories

Lattice ℤ 3 theories with complex actions share many key features with finite- density QCD including a sign problem and $\mathcal{CK}$ symmetry. Complex ℤ 3 spin and gauge models exhibit a generalized Kramers-Wannier duality mapping them onto chiral ℤ 3 spin and gauge models, which are simulatable with standard lattice methods in large regions of parameter space. The Migdal-Kadanoff real-space renormalization group (RG) preserves this duality, and we use it to compute the approximate phase diagram of both spin and gauge ℤ 3 models in dimensions one through four. Chiral ℤ 3 spin models are known to exhibit a Devil’s Flower phase structure, with inhomogeneous phases that can be thought of as ℤ 3 analogues of chiral spirals. Out of the large class of models we study, we find that only chiral spin models and their duals have a Devil’s Flower structure with an infinite set of inhomogeneous phases, a result we attribute to Elitzur’s theorem. We also find that different forms of the Migdal-Kadanoff RG produce different numbers of phases, a violation of the expectation for universal behavior from a real-space RG. We discuss extensions of our work to ℤ N models, SU(N) models and nonzero temperature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Diastereomers and Low-Temperature Oxidation

Diastereomers have historically been ignored when building kinetic mechanisms for combustion. Low-temperature oxidation kinetics, which continues to gain interest in both combustion and atmospheric communities, may be affected by the inclusion of diastereomers in radical chain-branching pathways. In this work, key intermediates and transition states lacking stereochemical specification in an existing diethyl ether low-temperature oxidation mechanism were replaced with their diastereomeric counterparts. Rate coefficients for reactions involving diastereomers were computed with ab initio transition state theory master equation calculations. The presence of diastereomers increased rate coefficients by factors of 1.2-1.6 across various temperatures and pressures. Ignition delay simulations incorporating these revised rate coefficients indicate that the diastereomers enhanced the overall reactivity of the mechanism by almost 15% and increased the peak ketohydroperoxide concentration by 30% in the negative temperature coefficient region at combustion-relevant pressures. Furthermore, these results provide an illustrative indication of the important role of stereomeric effects in oxidation kinetics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Parallelized POD-based suboptimal economic model predictive control of a state-constrained Boussinesq approximation

Motivated by an energy efficient building application, we want to optimize a quadratic cost functional subject to the Boussinesq approximation of the Navier-Stokes equations and to bilateral state and control constraints. Since the computation of such an optimal solution is numerically costly, we design an efficient strategy to compute a sub-optimal (but applicationally acceptable) solution with significantly reduced computational effort. We employ an economic Model Predictive Control (MPC) strategy to obtain a feedback control. The MPC sub-problems are based on a linear-quadratic optimal control problem subjected to mixed control and state constraints and a convection-diffusion equation, reduced with proper orthogonal decomposition. Finally, to solve each sub-problem, we apply a primal-dual active set strategy. The method can be fully parallelized, which enables the solution of large problems with real-world parameters.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Multireference Equation-of-Motion Driven Similarity Renormalization Group: Theoretical Foundations and Applications to Ionized States

We present a formulation and implementation of an equation-of-motion (EOM) extension of the multireference driven similarity renormalization group (MR-DSRG) formalism for ionization potentials (IP-EOM-DSRG). The IP-EOM-DSRG formalism results in a Hermitian generalized eigenvalue problem, delivering accurate ionization potentials for strongly correlated systems. The EOM step scales as O(N 5 ) with the basis set size N, allowing for efficient calculation of spectroscopic properties, such as transition energies and intensities. The IP-EOM-DSRG formalism is combined with three truncation schemes of the parent MR-DSRG theory: an iterative nonperturbative method with up to two-body excitations [MR-LDSRG(2)] and second- and third-order perturbative approximations [DSRG-MRPT2/3]. We benchmark these variants by computing (1) the vertical valence ionization potentials of a series of small molecules at both equilibrium and stretched geometries; (2) the spectroscopic constants of several low-lying electronic states of the OH, CN, N 2 + , and CO + radicals; and (3) the binding curves of low-lying electronic states of the CN radical. A comparison with experimental data and theoretical results shows that all three IP-EOM-DSRG methods accurately reproduce the vertical ionization potentials and spectroscopic constants of these systems. Notably, the DSRG-MRPT3 and MR-LDSRG(2) versions outperform several state-of-the-art multireference methods of comparable or higher cost.

Hamiltonians↗
DOE OSTI2021

Equation-of-Motion Coupled-Cluster Cumulant Green’s Function for Excited States and X-ray Spectra

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states and spectra. Conventional methods using the Dyson equation or the cumulant expansion are typically based on the GW self-energy approximation. In order to extend this approximation in molecular systems, a non-perturbative real-time coupled-cluster cumulant Green’s function approach has been introduced, where the cumulant is obtained as the solution to a set of coupled first order, non-linear differential equations. This approach naturally includes non-linear corrections to conventional cumulant Green’s function techniques where the cumulant is linear in the GW self-energy. The method yields the spectral function associated with the core Green’s function, which is directly related to the x-ray photoemission spectra (XPS) of molecular systems. The approach also yields very good results for binding energies and satellite excitations. The x-ray absorption spectrum (XAS) is then obtained as a convolution of the core spectral function and an effective one-body absorption spectrum. Here this approach is extended to include the full coupled-cluster-singles (CCS) core Green’s function by including the complete form of the non-linear contributions to the cumulant as well as all single, double, and triple cluster excitations in the CC amplitude equations. This approach naturally builds in orthogonality and shake-up effects analogous to those in the Mahan-Noizeres-de Dominicis edge singularity corrections that enhance the XAS near the edge. Themethod is illustrated for the XPS and XAS of NH 3 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

The repulsive Casimir-type forces in binary mixture of weakly interacting Bose–Einstein condensates gases

Highlights: • The equations of the state of a binary mixture weakly interacting Bose gases are established based on the Lagrangian of the system. • The formulae of the Casimir-type energy and the resulting Casimir-type force are produced from the quantum part of the grand canonical energy. • The results for the Casimir-type effect in the one-field case can be re-obtained by using the results of the two-field case. Within the framework of the quantum field theory (QFT) in the one-loop approximation, we investigate the Casimir-type effect in a mixture of two-component weakly interacting Bose-Einstein condensates (BECs) gases confined between two parallel plates at zero temperature. The Casimir-type forces in the miscible state of BECs for different boundary conditions (BCs) are discovered and compared to each other. The repulsive nature of the Casimir-type force in BECs is defined for the first time. Our calculations are finished in both grand canonical ensemble (GCE) and canonical ensemble (CE).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Uncertainty quantification for Joule heating processes in fibrous pore-resolved media

Joule heating (JH) is an energy-efficient and sustainable technique for heating materials. Its application for industrial heating, particularly, has been gaining attention due to its potential for increasing the yield of various chemical products. The process involves the use of heating elements (materials that are highly conductive electrically and thermally) to heat up other materials or substances. These conductors, however, can exhbit varying degrees of uncertainty due to non-linearities in their temperature-dependent properties, which could result in variable material behavior. In this work, we carry out uncertainty quantification (UQ) at the pore scale to describe the uncertainty of such materials. In so doing, we applied the non-intrusive polynomial chaos expansion (PCE) technique to quantify the uncertainty within the system. The steady state Joule heating equation was solved numerically at the pore scale mimicking conditions within a heating chamber for propane dehydrogenation, and various electro-thermal profiles were obtained. We also examined the effect of the number of sampling points (20 – 100) and order of the PCE coefficients (2 – 5) on the accuracy of the temperature evaluations. The results were then benchmarked with the standard Monte Carlo (MC) method. The average temperature of the 4th-order global PCE showed good agreement with the MC results (which were positively skewed). Orders greater than 4 gave an underestimation of the temperatures while predictions for the peak temperature improved as the number of sampling points increased.

Fagbemi, Samuel [ORNL] (ORCID:0000000236995025)↗
DOE OSTI2023

Design of a prototypical natural circulation water-based reactor cavity cooling system (RCCS) for a pebble-bed generic FHR

A prototypical natural circulation water-based reactor cavity cooling system (RCCS) for the UC Berkeley Mark-1 pebble-bed gFHR is designed based on one-dimensional thermal hydraulics modeling and optimization implemented in an in-house MATLAB code. The model employs a lumped core consisting of fuel pebble and graphite pebble regions, while the graphite reflector, vessel, and RCCS are represented as separate but energy-coupled regions. The model is derived based on steady state energy balance equations accounting for conductive heat transfer from the lumped core to the vessel and subsequent radiative heat transfer to the RCCS and convective heat transfer in the water. Further, mass flow rate is calculated based on momentum and thermal energy balance in the RCCS. Starting with ANL’s water NSTF as baseline, effects of RCCS design parameters such as size of the plate between pipes, pipe diameter, source-to-sink distance, surface emissivity, and pitch from the core are examined in an integral effects framework which accounts for subsequent variation of the mass flow rate, number of pipes, convection coefficient, mean region and surface temperatures, radiative view factors, and water outlet temperature. The importance of including a prototypical reactor design in RCCS design calculations is emphasized as the calculations show that temperature drops from the core to the RCCS fluid to transfer the heat tend to be the limitation not the heat removal capacity of the fluid. A closed-loop design is obtained with capability to safely remove up to 0.72% of the nominal reactor power at maximum estimated peak conditions and 0.40% at shutdown. The design is based on physics calculations and does not account for economic optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Including the vacuum energy in stellarator coil design

Being three-dimensional, stellarators have the advantage that plasma currents are not essential for creating rotational-transform; however, the external current-carrying coils in stellarators can have strong geometrical shaping, which can complicate the construction. Reducing the inter-coil electromagnetic forces acting on strongly shaped 3D coils and the stress on the support structure while preserving the favorable properties of the magnetic field is a design challenge. In this work, we recognize that the inter-coil forces are the gradient of the vacuum magnetic energy. We introduce an objective functional built on the usual quadratic flux on a prescribed target surface together with a weighed penalty on the vacuum energy. The Euler–Lagrange equation for stationary states is derived, and numerical illustrations are computed using a modern stellarator optimization framework. A study of the effect of the energy functional on the inter-coil forces is conducted and the energy is shown to be a promising quantity in producing coils with low forces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗