Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Equilibrium programming”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Monte Carlo Hauser-Feshbach computer code system to model nuclear reactions: YAHFC

A computer program framework, YAHFC, to model low-energy nuclear reactions is presented. The framework allows for reactions with incident particles ranging from protons/neutrons to alphas and is designed to address reactions that ultimately lead to the formation of compound nuclear systems that then decay statistically as outlined in concepts of Hauser and Feshbach. Additionally, instead of a reaction, it is also possible to model the decay of a nuclear system with an initial excitation and population. The code models nuclear decays with a Monte Carlo process that tracks the decay of each state. This allows for an exact representation of the spectra for all emitted particles in each of the final exit channels and the possibility of generating reaction data for simulation purposes. The program is interfaced with the optical model code system FRESCOX to calculate transmission coefficients as well as the effects of coupled channels and other direct excitations via the distorted wave Born approximation (DWBA). Modules are included to account for nuclear processes such as width corrections, pre-equilibrium emission, and fission. The program is controlled by a series of input commands and while a set of input parameters exists for each projectile and target, the input commands allow for complete control over each input parameter. Extensive data files are produced and a program is provided that converts YAHFC data files into nuclear data library entries in the generalized nuclear data structure (GNDS).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Development of a Performance Portable Non-Equilibrium Plasma Fluid Solver on Adaptive Grids

This presentation will describe the numerical techniques, programming paradigms, verification, and performance of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures. Our plasma fluid model solves the conservation equations for self-consistent electrostatic Poisson, electron and heavy species transport, and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive mesh management library, AMReX (Zhang et al., JOSS, 4 (37) 1370, 2019), and can be built and run on widely available vendor specific GPU architectures (NVIDIA/AMD/Intel). We utilize a non-subcycled second order semi-implicit time-stepping method where all adaptive mesh refinement (AMR) levels are advanced with the same time step. The composite multi-level multigrid solver from within AMReX is used for each of the governing equations that are cast into a Helmholtz equation form. We have also developed a python based chemical mechanism parser framework that uses a similar format as CANTERA (Goodwin et al., Zenodo, 2018) yaml files as input. Our custom parser reads the yaml file and provides C++ files with transport and production rate functions that can be executed on both host (CPU) and device (GPU). We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on low-pressure capacitive and high-pressure streamer discharges. Our initial performance studies indicate 10X speed-up using 20 NVIDIA GPUs versus 200 CPUs for an atmospheric streamer discharge problem solved on a 512 x 1024 x 512 grid.

graphics processing units↗
DOE OSTI2025

Convex relaxation for Fokker–Planck equation

We propose an approach to directly estimate the moments or marginals for a high-dimensional equilibrium distribution in statistical mechanics by solving the high-dimensional Fokker–Planck equation in terms of low-order cluster moments or marginals. With this approach, we bypass the exponential complexity of estimating the full high-dimensional distribution and directly solve the simplified partial differential equations for low-order moments/marginals. Moreover, the proposed moment/marginal relaxation is fully convex and can be solved via off-the-shelf solvers. We further propose a time-dependent version of the convex programs to study non-equilibrium dynamics. In a specific setting, we show the proposed method can recover a mean-field-type equilibrium density. Numerical results are provided to demonstrate the performance of the proposed algorithm for high-dimensional systems.

Chen, Yian↗
DOE OSTI2025

Adaptive PID Gain Scheduling Control for Hydropower Turbine Using Neural CDE and Stochastic Distribution Shaping

This paper introduces a gain-scheduling PID controller design strategy for hydroturbine frequency control mode. This scheme first uses real data to learn the nonlinear dynamics of the hydroturbine using neural controlled differential equations and then perturbs the obtained nonlinear system at different equilibrium points, based on which a static output feedback adaptive dynamic programming algorithm is then used to optimize the PID gains for each equilibrium point. Moreover, a continuous-time version of stochastic distribution control is proposed to further fine-tune the optimized PID gains. Finally, the controller is obtained by implementing linear interpolation between the optimized PID control gains. The simulation results show that the proposed gain-scheduling PID controller can control a larger range of operation points compared with the given fixed PID controller and the baseline method. Compared with the given fixed PID controller, the proposed gain-scheduling PID controller can regulate hydroturbine frequency against disturbances induced by power-load variation with over 50% less overshoot for some operation points.

13 HYDRO ENERGY↗
DOE OSTI2022

Q-POP-Thermo: A general-purpose thermodynamics solver for ferroelectric materials

We report that Q-POP-Thermo is a program designed to compute thermodynamic monodomain equilibrium states and their properties for ferroelectric single crystals and thin films based on the Landau-Ginzburg-Devonshire (LGD) Theory. Utilizing symbolic manipulation with the SymPy Library, the governing equations along with appropriate boundary conditions are solved for speedy minimization of the free energy of a crystal. Utilizing the popular Differential Evolution algorithm, with appropriate hybridization, multiple phase diagrams, such as the pressure-temperature phase diagram for bulk single crystals and the common strain-temperature phase diagram for monodomain thin-film systems can be readily generated. Furthermore, a variety of material properties of stable ferroelectric phases, including dielectric, piezoelectric, and electrocaloric properties, can simultaneously be calculated. Validation studies are presented for both thin-film and single crystal systems to test the effectiveness and capability of the open-source program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Compiling Adiabatic Quantum Programs

We develop a non-cooperative game-theoretic model for the problem of graph minor-embedding to show that optimal compiling of adiabatic quantum programs in the sense of Nash equilibrium is possible.

Khan, Faisal Shah↗
DOE OSTI2019

Cryogenic Carbon Capture Development (Final Report)

This project extends and leverages the work completed under funding from ARPA-E and other sources by exploring aspects of the Cryogenic Carbon Capture™ (CCC) process and issues discovered during those projects, but that were not within the scope or budget of previous projects or that became apparent during field tests. Our successful completion of the tasks outlined in this project has increased the reliability, efficiency, and scalability of the CCC process and we are better prepared for pilot-scale testing. The first phase of this project improved key areas of the process through iterative design and experiment, culminating with recommendations for improvements to the existing skid-scale system developed under non-federal funding. The second phase of this project implemented the improved unit operations into the Cryogenic Carbon Capture External Cooling Loop (CCC ECL™) skid and field tests the improved system at a coal-fired power plant. The specific areas of work in this project are: Budget Period 1: 1) Selection of the most energy-efficient drying system for ambient-pressure flue gas. Sustainable Energy Solutions, LLC (SES) has successfully demonstrated a dryer that removes bulk water using a low-temperature brine to remove about 99% of the water followed by a mol sieve 3A desiccant dryer to remove the remaining 1% water; 2) Selection and demonstration of efficient solution to fouling due to dissolved CO 2 . SES has tested several commercial heat exchangers that purport to be able to handle fouling issues similar to those experienced in the contact liquid chiller. Replicated tests using the CCC ECL™ skid and dimpled brazed plate heat exchangers demonstrate that quick throttling of refrigerant combined with a burst of clean gas mitigate the fouling caused by CO 2 solubility in contact liquid with only minor increases in parasitic load. SES models indicate that parasitic load would increase by approximately 0.0125 GJ/tonne, an order of magnitude less than what was determined to be the maximum allowable additional parasitic load caused by this system; 3) Selection and demonstration of efficient and scalable solid–liquid separation system. The screw press technology met or exceeded the 70% solids loading benchmark during several test runs exceeding 24 hours on the CCC ECL™ skid. The average solids loading during the entire 50-hour run was 67%. The energy demand for this system is less than 0.1 GJ e /tonne of CO 2 processed; 4) Optimize design and demonstration of patented desublimating heat exchangers. We designed and built a hybrid spray tower that was tested over many hundreds of hours and a 50-hour continuous test run was completed capturing over 0.5 tonne CO 2 /day at >90% capture; (5) Plan for safe dilution of oxygen-depleted, cleaned flue-gas stream. The AERSCREEN simulator developed by the EPA modeled the dispersion of oxygen-depleted flue gas from the CCC process at a commercial plant. The results indicate that this plan effectively mitigates the flue gas hazards; (6) Simulation of and experimental work related to multi-pollutant capture. SES programmed an upgrade to its solid/vapor/liquid equilibrium (SVLE) algorithm. This enables simultaneous modeling of multiple pollutants (i.e., CO 2 , SO x , NO x , Hg, HCl). Budget Period 2; (7) Modify the CCC ECL™ skid for field testing and complete shakedown testing. The plans for modifying the CCC ECL™ skid were finalized and put into place. Shakedown testing of the updated CCC ECL™ skid occurred at SES using closed-loop mode (i.e., using a synthetic flue gas consisting of only N 2 and CO 2 ) and at the Hunter Power Plant using flue gas from the coal-fired boiler. Shakedown testing took longer than expected due to delays in receiving equipment and unexpected equipment performance in the field; (8) Operate the CCC ECL™ skid in the field. Longer testing of the CCC ECL™ skid occurred at the Hunter Power Plant. Multiple tests exceeding 35 hours were completed, with a total of over 600 hours of CO 2 capture completed using flue gas at the Hunter Power Plant with an average CO 2 capture rate of over 91% during those hours. Extensive notes and documentation were acquired informing improved methods for controls, startup, operation, and shut down; (9) Determine figures of merit and quantify process improvement. SES defined figures of merit for unit operation performance and overall CCC ECL™ process performance. All figures of merit were achieved or surpassed, except for 500 continuous hours of demonstration at the Hunter Power Plant; (10) Complete an Environmental Health and Safety Risk assessment by a third party. Tri-State completed an Environmental Health and Safety Risk assessment; (11) Update techno-economic analysis with third-party validation. The techno-economic analysis has been completed. It confirms the potential large reduction in energy and cost of the CCC process compared with leading alternatives and that the estimates and estimation techniques are reliable.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2018

Hedging Strategies: Electricity Investment Decisions under Policy Uncertainty

Given uncertainty in long-term carbon reduction goals, how much non-carbon generation should be developed in the near-term? This research investigates the optimal balance between the risk of overinvesting in non-carbon sources that are ultimately not needed and the risk of underinvesting in non-carbon sources and subsequently needing to reduce carbon emissions dramatically. We employ a novel framework that incorporates a computable general equilibrium (CGE) model of the U.S. into a two-stage stochastic approximate dynamic program (ADP) focused on decisions in the electric power sector. We solve the model using an ADP algorithm that is computationally tractable while exploring the decisions and sampling the uncertain carbon limits from continuous distributions. The results of the model demonstrate that an optimal hedge is in the direction of more non-carbon investment in the near-term, in the range of 20–30% of new generation. We also demonstrate that the optimal share of non-carbon generation is increasing in the variance of the uncertainty about the long-term carbon targets, and that with greater uncertainty in the future policy regime, a balanced portfolio of non-carbon, natural gas, and coal generation is desirable.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2026

Supervisory Control and Data Acquisition for Electrochemical Separation Experimentation

The Python-based program is a laboratory automation tool designed to control and monitor electrochemical systems. The tool was developed for capacitive deionization (CDI) experiments, but it can be used for any system that requires controlled voltage or current segments and multi-parameter monitoring. The program integrates hardware components to run user-defined experimental parameters, providing operational control of a programmable power supply, peristaltic pump, and data acquisition devices. Currently, the program is structured with a workflow that includes an initialization (or pre-run) phase, a main loop, and a post-experiment stabilization (or post-run) phase. The initialization phase prepares and stabilizes the cell, ensuring that the electrodes and solution reach a baseline state before the experiment begins. The main loop consists of multiple voltage segments that repeat, controlling the experiment while recording key parameters such as time, voltage, current, pH, and conductivity. Finally, the post-experiment stabilization phase allows the system to stabilize after the experiment, returning the cell and solution to equilibrium conditions before ending the sequence. The program is designed with four variations, each tailored to different experimental needs. All variations include both the initialization and post-experiment stabilization stages, which run for a set amount of time, voltage, current, and flow rate before and after the main experiment block. The main loop runs for a set number of cycles, as defined by the user input, and each cycle is composed of 2 or 4 segments. The 4 program variations are described as follows: Program 1: The main program includes 2 segments. Each segment is defined to have a set duration, flow rate, voltage, and current. This program measures conductivity, flow rate, voltage, and current. Program 2: The main program expands Program 1 to include 4 segments. Each segment has a specified duration, flow rate, voltage, and current. Like Program 1, it measures conductivity, flow rate, voltage, and current. Program 3: The main program consists of 2 segments, each defined by time, flow rate, voltage, and current. In addition to conductivity, flow rate, voltage, and current, Program 3 collects pH and temperature data through a 4-channel data acquisition device. Program 4: This program independently controls two channels of a multi-channel power supply simultaneously. While conductivity can only be measured for one cell at a time, the dual-channel control makes it possible to operate two cells simultaneously under different voltage/current conditions. The main program includes 2 segments.For each program, all measurements are automatically logged and integrated into a single Excel output file. Data are displayed in numerical format and plotted, both in real time, to track system performance. A key feature of the program is its ability to synchronize all outputs so that every measurement shares a single timestamp, ensuring accurate alignment of voltage, current, pH, conductivity, and pH data.By combining hardware control, real-time monitoring, and unified data collection, this program significantly reduces manual workload and minimizes errors, making it a reliable platform for researchers, engineers, and laboratory technicians conducting CDI experiments, among other electrochemical tests.

Valentino, Lauren [Argonne National Laboratory (AN↗
DOE OSTI2023

An energy-optimization method to study gel-swelling in confinement

Using our open-source programming environment Morpho, an energy-optimization method is developed to compute equilibrium shapes of an arbitrarily confined swollen hydrogel. The internal strains and contact forces are measured, probing its mechanics.

Joshi, Chaitanya↗
DOE OSTI2023

Extending the thermodynamic form factor bootstrap program: multiple particle-hole excitations, crossing symmetry, and reparameterization invariance

In this study, we further the thermodynamic bootstrap program which involves a set of recently developed ideas used to determine thermodynamic form factors of local operators in integrable quantum field theories. These form factors are essential building blocks for dynamic correlation functions at finite temperatures or non-equilibrium stationary states. In this work we extend this program in three ways. Firstly, we demonstrate that the conjectured annihilation pole axiom is valid in the low energy particle-hole excitations. Secondly, we introduce a crossing relation, which establishes a connection between form factors with different excitation content. Typically, the crossing relation is a consequence of Lorentz invariance, but due to the finite energy density of the considered states, Lorentz invariance is broken. Nonetheless a crossing relation involving excitations with both particles and holes can established using the finite volume representation of the thermodynamic form factors. Finally, we demonstrate that the thermodynamic form factors satisfy a reparameterization invariance, an invariance which encompasses crossing. Reparameterization invariance exploits the fact that the details of the representation of the thermodynamic state are unimportant. In the course of developing these results, we demonstrate the internal consistency of the thermodynamic form factor bootstrap program in a number of ways. Finally, we provide explicit computations of form factors of conserved charges and densities with crossed excitations and show our results can be used to infer information about thermodynamic form factors in the Lieb-Liniger model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Solving Stochastic Dynamic Programming Problems: A Mixed Complementarity Approach

We detail a mixed complementarity problem (MCP) formulation of continuous state dynamic programming problems (DP-MCP). We write the solution to projection methods in value function iteration (VFI) as a joint set of optimality conditions that characterize maximization of the Bellman equation; and approximation of the value function. The MCP approach replaces the iterative component of projection based VFI with a one-shot solution to a square system of complementary conditions. We provide three numerical examples to illustrate our approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Optimal Equilibrium Selection of Price-Maker Agents in Performance-Based Regulation Market

This paper analyzes the oligopolistic equilibrium of multiple price-maker agents in performance-based regulation (PBR) markets. In these markets, there are price-maker agents representing some frequency regulation (FR) providers and a number of independent price-taker FR providers. An equilibrium problem with equilibrium constraints (EPEC) model is employed in this paper to study the equilibria of a PBR market in the presence of price-maker agents and pricetaker FR providers. Due to the incorporation of the FR providers' dynamics, the proposed model is reformulated as a mixed-integer linear programming (MILP) problem over innovative mathematical techniques. An optimal equilibrium point is also selected for the market, where none of the agents is unique deviator and the dynamic performance of power system is improved simultaneously. The effectiveness of proposed optimal equilibrium point is evaluated in the numerical results section by comparing the outputs with conventional optimal dispatches of the FR providers.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Off-of-equilibrium effects on Kurtosis Along Strangeness-Neutral Trajectories

The Beam Energy Scan program at the Relativistic Heavy Ion Collider (RHIC) is searching for the QCD critical point. The main signal for the critical point is the kurtosis of the distribution of proton yields obtained on an event by event basis where one expects a peak at the critical point. However, its exact behavior is still an open question due to out-of-equilibrium effects and uncertainty in the equation of state. Here we use a simplistic hydrodynamic model that enforces strangeness-neutrality, selecting trajectories that pass close to the critical point. We vary the initial conditions to estimate the effect of out-of-equilibrium hydrodynamics on the kurtosis signal.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Spin-selective photocatalysis and quantum transport using ab initio density-matrix dynamics

The overarching goal of this research program is to quantitatively predict out-of-equilibrium quantum dynamics and transport with electrons and spin degrees of freedom coupled to phonons, photons, defects and liquid environments, specifically targeting spin-selective chemistry. Spin polarized carriers, which can be generated at room temperature by circularly-polarized light, provide unique opportunities for creating new reaction pathways and increasing selectivity towards specific pathways. If realized, this can substantially increase the level of reaction control possible in the field of heterogeneous catalysis. Spin, a quantum mechanical entity without a classical counterpart, requires quantum dynamics. Capturing both coherent and incoherent dynamics and accurate correlation of spins, electrons and phonon environments are essential to quantitative description of generation, relaxation and transport, as well as chemical reaction of spin-polarized carriers in molecules and materials promising for spin-photocatalysis.

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

Mix and match—a versatile equilibrium approach for hybrid colloidal synthesis

Colloidal synthesis is a powerful bottom-up approach for programmed self-assembly which holds promise for both research and industry. While diverse, each synthetic process is typically restricted to a specific chemistry. Many applications however require composite materials, whereas a chemical equilibrium can typically only match one material but not the other. Here, a scalable general approach is presented, alleviating the dependency on a specific chemical reaction, by resorting to a mechanical equilibrium; an isopycnic density-gradient-step is tailored to form clusters with prescribed composition. Valence control is demonstrated, making dimers, trimers, and tetramers with purity as high as 96%. The measured kinetics shows a scaleable throughput. The density gradient step plays a dual role of both filtering out undesired products and concentrating the target structures. The “Mix-and-Match” approach is general, and applies to a broad range of colloidal matter: diverse material compositions (plastics, glasses, and emulsions); a range of colloidal interactions (van der Waals, Coulomb, and DNA hybridization); and a spectrum of sizes (nanoscale to multiple micrometers). Lastly, the strength of the method is displayed by producing a monodisperse suspension from a highly polydisperse emulsion. The ability to combine colloids into architectures of hybrid materials has applications in pharmaceuticals, cosmetics, and photonics.

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

Thermodynamic Model of COREX Melter Gasifier Using FactSage™ and Macro Facility

In the current study, the thermochemical model of the COREX melter gasifier has been developed using the FactSage™ thermodynamic software, and macro programming facility. The melter gasifier was conceptually divided into various equilibrium zones dealing with the reactions happening in the reactor. The coal char and the coke being used are divided among them optimally. The model has been validated with respect to the plant data for the hot-metal, slag, and reducing gas composition. The effects of the variation in the degree of metallization, on the final hot metal, slag, and reducing gas, have been predicted using this model. The variations in the process parameters for the different industrial coal grades used in the COREX plant have been determined. The effects of the volatile, ash, and the moisture contents of the coal on the hot-metal and the reducing gas compositions have been reported.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Power-traffic coordinated operation for bi-peak shaving and bi-ramp smoothing – A hierarchical data-driven approach

The power distribution system and urban transportation system are two networked system bare their own operation constraints, such peak load in power systems and traffic congestion in transportation system. With the increasing number of electrical vehicles and charging/discharging stations, two systems are become tightly coupled. However, to optimize the two systems target using electrical vehicles as decision control variables cannot be easily solved using a uniformed optimization frame work. Thus we propose a hierarchical optimization approach to address this problem, which consists of a higher and a lower level. In the higher level, the power distribution system and urban transportation system are treated together to minimize the social cost. Meanwhile, the electrical vehicles and the charging/discharging stations are treated as customers to minimize their own expenditures. Then, an equilibrium is designed to determine the optimal charging/discharging price. In the lower level, the models of power distribution system and urban transportation system are developed to provide a detailed analysis. Specifically, in power distribution system, the three-phase unbalanced optimal power flow problem is relaxed with the semidefinite relaxation programming, and solved with alternating direction method of multiplier. A dynamic user equilibrium problem is formulated for the urban transportation system. For electrical vehicles, the state of charge is considered to optimize the charging/discharging schedule and reduce the impacts of power distribution systems. We conducted the simulation and numerical analysis using the IEEE 8500-bus distribution system and the Sioux Falls system with about 10,000 cars. The results demonstrate the feasibility and effectiveness of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗