Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Dynamic process intensification: Fundamentals and implementation to ternary distillation

Distillation remains a key technology for separating liquid mixtures. Its versatility comes with the disadvantage of high energy consumption. We previously used empirical arguments to introduce dynamic process intensification (DPI) as a strategy for improving the energy efficiency of binary distillation. In this article, we focus on ternary distillation; we begin by providing a rigorous basis for DPI, then formulate the problem of identifying the operating states for DPI as a nonlinear optimization problem. Via an extensive case study considering a hydrocarbon mixture, we demonstrate that reboiler energy use can be reduced by more than 2.3% relative to an equivalent column operating at steady state. Here, we prove that this result is due to the fact that DPI aims to satisfy product flow and quality constraints on average in time, affording broader opportunities for optimizing column economics than steady-state operation, where constraints are met strictly at all times.

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

Coupled Microbial-Conversion and Computational-Fluid-Dynamics (CFD) Models for Butanediol Production in Micro-Aerated Reactors

Microbial conversion of substrates to macromolecules has been widely used in the synthesis of value-added products in pharmaceutical and biotechnology industries. These bioreactions are also being investigated in the production of low-value commodities such as biofuels [1]. Gas-liquid mass-transfer and transport-reaction coupling are important challenges when designing and scaling up these reactor systems. Experiments in wellmixed small-scale reactors have enabled characterization of microbial reactivity, while their coupling with macroscale transport remains relatively unexplored. In this work, we use a coupled metabolic-CFD model to study the action of a genetically engineered microbe Zymomonas mobilis [2] on sugars to produce 2,3-Butanediol (BDO). BDO is an important hydrocarbon intermediate that can be catalytically upgraded to several fuels and chemicals [3]. An important aspect to this particular microbial conversion is the need for micro-aerated environments as opposed to traditional aerobic fermentation. Slight variations in oxygen concentration can result in competing reaction pathways that disable BDO production. Hence, gas-liquid mass transfer and transport need to be optimized in large-scale reactors to maximize BDO production, for which CFD is a valuable tool. The aerobic-fermentation CFD model previously developed by the authors [5] for simulating bubble-column and airlift reactors at scale was used in this study. The Reynolds-averaged mass, momentum and energy transport equations for interpenetrating gas and liquid phase are solved in this model along with the transport and interphase mass transfer of oxygen. Our previous work used a phenomenological model for microbial oxygen uptake that neglected microbial growth and other reaction pathways. In this work, a detailed metabolic model enabled prediction of product formation and inhibition pathways. In order to manage computational cost, we used a subcycling technique [6] that takes advantage of the clear separation in transport (~ 200 sec) and reaction (~ 2-3 hours) timescales. The CFD model is first solved to steady state, after which the metabolic model is advanced at every cell in the computational domain using the local oxygen concentration. The CFD model is then run to achieve a new steady state that provides a new oxygen distribution for the metabolic model. This process, where reaction and fluid updates are interleaved together, is iterated until reactants are completely exhausted. This work will examine the performance of different reactor designs such as bubble column and airlift reactors at scale (250-500 m3). Oxygen mass-transfer coefficient and distribution are critically analyzed among reactors, and optimization studies pertaining to aeration is presented. Furthermore, it has been observed in experiments that high BDO production may be achieved by manipulating the aerobic environment over the course of reaction, such that oxygen concentration is high during the growth phase, and very low as sugar is depleted. This characteristic will be addressed by our simulations for which a timedependent scheduling strategy for aeration is presented that maximizes BDO production. [1] Humbird, D., Davis, R., and McMillan, J., Aeration costs in stirred-tank and bubble column bioreactors, Biochemical Engineering Journal, 127, 161—166, 2017 [2] Yang, S., Mohagheghi, A., Franden, M. A., Chou, Y.-C., Chen, X., Dowe, N., Himmel, M. E., and Zhang, M., Metabolic engineering of zymomonas mobilis for 2, 3-butanediol production from lignocellulosic biomass sugars. Biotechnology for biofuels, 9(1):189, 2016 [3] Kim, S. J., Sim, H. J., Kim, J. W., Lee, Y. G., Park, Y. C., and Seo, J. H., Enhanced production of 2,3-butanediol from xylose by combinatorial engineering of xylose metabolic pathway and cofactor regeneration in pyruvate decarboxylase-deficient Saccharomyces cerevisiae. Bioresource Technology, 245:1551–1557, 2017 [4] Weller, H., Tabor, G., Jasak, H. and Fureby, C., A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12, 6, 620--631, 1998

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Quantification of the importance of the modeling fidelity in the Best Estimate Plus Uncertainty analysis

In this work, the multi-physics coupling of CTF-PARCS for both steady-state and transient simulations is evaluated focusing on the quantification of the importance of the modeling fidelity when applied to the Best Estimate Plus Uncertainty safety analysis. For the steady-state, a Picard iteration approach is adopted in the coupled code while for the transient an explicit approach has been developed as a first step. The CTF-PARCS coupling is applied to the REA at HFP BOC scenario proposed in the exercise I-2c of LWR-UAM benchmark Phase III. Two different modeling choices are investigated for the nominal calculations. The first consists of including the 3D burnup information in CTF, and the second of modeling the spacers grids in both PARCS and CTF. It is observed that due to the degradation of the fuel thermal conductivity, the 3D burnup inclusion in CTF leads to significant higher centerline temperatures than the base case, with a difference up to 300 K. The spacer grids modeling shows the expected power and temperature decrease in the axial positions of the spacer grids. For the uncertainty quantification study, the multi-physics modeling including the burnup information in CTF is used. The results show that the inserted reactivity has a relative standard deviation of 13.6% mainly attributed to the nuclear data. The maximum linear power has a relative standard deviation of 16.2% due to the nuclear data as well, while the maximum centerline temperature has a relative standard deviation of 3.5% primarily caused by the fuel thermal conductivity. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Model-based and Model-free Designs for an Extended Continuous-time LQR with Exogenous Inputs

We present an extended linear quadratic regulator (LQR) design for continuous-time linear time-invariant (LTI) systems in the presence of exogenous inputs. We first propose a model-based solution with cost minimization guarantees for states and inputs using dynamic programming (DP). The control law consists of a combination of the optimal state feedback and an additional optimal term dependent on the exogenous inputs. The control gains for the two components are obtained by solving a set of matrix differential equations. We provide these solutions for both finite horizons and steady-state cases. In the second part of the paper, we formulate a reinforcement learning (RL) based algorithm which does not need any model information except the input matrix, and can compute an approximate steady-state LQR gain using measurements of the states, the control inputs, and the exogenous inputs. Both model-based and data-driven optimal control algorithms are tested with a numerical example under different exogenous inputs showcasing the effectiveness of the designs.

Mukherjee, Sayak↗
DOE OSTI2020

Cavity Swelling in Irradiated Materials

Cavity swelling, due to three-dimensional clustering of vacancies, is an undesirable isotropic volume expansion of materials under irradiation. This phenomenon occurs above stage-III recovery temperature where vacancies are mobile, between ~0.3–0.6 Tm where Tm is the absolute melting temperature. The primary mechanism for cavity swelling is “dislocation bias”, i.e., preferential absorption of interstitials by network dislocations, causing a vacancy supersaturation. Cavity swelling is highly temperature dependent and typically follows a bell-shaped curve with peak swelling occurring at an intermediate temperature. With increasing temperature, cavity size increases, and number density decreases logarithmically. With increasing dose, cavity swelling increases, but follows a hockey-stick type pattern, characterized by a low-swelling transient incubation period and then a steady-state regime. The incubation period, dominant in cavity nucleation, is dependent upon material’s chemistry and irradiation conditions such as dose rate/temperature. Cavity swelling increases linearly with dose in the steady-state regime, that is independent of chemistry and dose rate. This regime is dominated by cavity growth. Transmutation gases like He enhance cavity nucleation such that an early onset of cavity formation occurs. However, swelling magnitudes can be either increased or suppressed when He is present. Results in austenitic steels and Cu indicate that cavity swelling peaks at intermediate He/dpa ratios of ~10 appm He/dpa. The most effective method to mitigate cavity swelling is by increasing the fixed point-defect sink density, by promoting precipitation or dispersion strengthening.

Bhattacharya, Arunodaya↗
DOE OSTI2023

Air quality and comfort constrained energy efficient operation of multi-zone buildings

Maintaining indoor air quality (IAQ) through effective ventilation is essential for the well-being and productivity of building occupants. Control strategies aimed at improving the efficiency of heating, ventilation and air conditioning (HVAC) systems must jointly determine ventilation and heating and cooling processes. Here, in this paper, we study the problem of minimizing the energy consumption of the HVAC system in a multi-zone building, while meeting thermal comfort and IAQ requirements. We first perform a steady state analysis of the zonal carbon dioxide (CO 2 ) concentration and the temperature dynamics. The resulting expressions are convex in the zonal mass flow rates and zonal temperatures. Guided by the steady state solutions for meeting the thermal comfort constraints, we develop two control policies for improving the energy efficiency of building HVAC systems while jointly satisfying indoor temperature and IAQ constraints. We compare the performance of our proposed approaches with those of multiple baseline approaches which implement separate regimes for controlling zonal temperature and IAQ for a typical work-day in a multi-zone campus building. We have evaluated the performance of our proposed approaches under varying levels of flexibility in zonal temperatures. We have shown that zonal temperature flexibility can result in energy savings up to 32% (for the same control strategies) as compared to the case where no such flexibility is permitted. Our proposed approaches were seen to offer potential savings of nearly 29% compared to the baseline.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Design, fabrication, and testing of a water-cooled tungsten-foam high-heat flux plasma-facing module

Here, we present a comprehensive study detailing the design, fabrication and testing a new modular “finger” component, featuring a hexagonal tungsten foam armor tile bonded to the top surface of the water-cooled finger unit. The mechanical, thermal, and particle transport properties of the tungsten foam armor were determined for a wide range of porosity and cell size, and then combined in an equivalent solid material with “effective” properties suitable for large-scale simulations. A series of 3-D multiphysics simulations for coupled thermal and mechanical behavior of the “finger” module are used to optimize the design before fabrication. The simulation results show that, for the expected experimental conditions, a safe range of heat flux and foam porosity is possible for both water-cooling and helium-cooling. Based on these multiphysics simulations, an optimum finger cooling module has been fabricated at Ultramet, Inc. A series of experiments were designed to study the thermomechanical performance of the “finger” unit under steady-state and transient thermal loading at surface heat flux up to 16 MW/m 2 in the High Energy Flux Test facility (HEFTY) at UCLA. Testing conditions included low-power, steady-state plasma exposure at 6 MW/m 2 , followed by severe high heat flux cycles at 16 MW/m 2 . Results of these tests are explained through comparisons with multiphysics simulations. To eliminate the effects of thermomechanical damage carried out on a module on its thermo-mechanical response on a subsequent test, three different yet identically designed and fabricated finger modules were tested. The results show that the current design is capable of cooling the walls of fusion devices at surface heat flux up to 6 MW/m 2 without failure, and that the failure heat flux limit is below 16 MW/m 2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Generic Primary-Control Model for Grid-Forming Inverters: Towards Interoperable Operation & Control

This paper outlines an architectural vision centered around the notion of interoperability to integrate grid-forming inverter-based resources in large-scale grids. With the underlying principle of interoperability guiding developments, we focus on modeling the characteristics of droop, virtual synchronous machine, and virtual oscillator controls. Emphasis is placed on these three controllers since they are leading grid-forming control candidates and are likely to be commonplace as primary-control schemes in future systems. We show that these controllers can each be considered as instantiations of a more generic model and that all these controllers exhibit similar droop-like relations between pertinent terminal variables in steady state. This commonality between controllers gives interoperability among them such that automatic synchronization, power sharing, and voltage regulation can be achieved. Simulation results validate the models and demonstrate how the steady-state droop characteristics of these control methods can be aligned with the aid of the developed modeling paradigm.

control↗
DOE OSTI2021

Deep Learning Approaches to Surrogates for Solving the Diffusion Equation for Mechanistic Real-World Simulations

In many mechanistic medical, biological, physical, and engineered spatiotemporal dynamic models the numerical solution of partial differential equations (PDEs), especially for diffusion, fluid flow and mechanical relaxation, can make simulations impractically slow. Biological models of tissues and organs often require the simultaneous calculation of the spatial variation of concentration of dozens of diffusing chemical species. One clinical example where rapid calculation of a diffusing field is of use is the estimation of oxygen gradients in the retina, based on imaging of the retinal vasculature, to guide surgical interventions in diabetic retinopathy. Furthermore, the ability to predict blood perfusion and oxygenation may one day guide clinical interventions in diverse settings, i.e., from stent placement in treating heart disease to BOLD fMRI interpretation in evaluating cognitive function (Xie et al., 2019; Lee et al., 2020). Since the quasi-steady-state solutions required for fast-diffusing chemical species like oxygen are particularly computationally costly, we consider the use of a neural network to provide an approximate solution to the steady-state diffusion equation. Machine learning surrogates, neural networks trained to provide approximate solutions to such complicated numerical problems, can often provide speed-ups of several orders of magnitude compared to direct calculation. Surrogates of PDEs could enable use of larger and more detailed models than are possible with direct calculation and can make including such simulations in real-time or near-real time workflows practical. Creating a surrogate requires running the direct calculation tens of thousands of times to generate training data and then training the neural network, both of which are computationally expensive. Often the practical applications of such models require thousands to millions of replica simulations, for example for parameter identification and uncertainty quantification, each of which gains speed from surrogate use and rapidly recovers the up-front costs of surrogate generation. We use a Convolutional Neural Network to approximate the stationary solution to the diffusion equation in the case of two equal-diameter, circular, constant-value sources located at random positions in a two-dimensional square domain with absorbing boundary conditions. Such a configuration caricatures the chemical concentration field of a fast-diffusing species like oxygen in a tissue with two parallel blood vessels in a cross section perpendicular to the two blood vessels. To improve convergence during training, we apply a training approach that uses roll-back to reject stochastic changes to the network that increase the loss function. The trained neural network approximation is about 1000 times faster than the direct calculation for individual replicas. Because different applications will have different criteria for acceptable approximation accuracy, we discuss a variety of loss functions and accuracy estimators that can help select the best network for a particular application. We briefly discuss some of the issues we encountered with overfitting, mismapping of the field values and the geometrical conditions that lead to large absolute and relative errors in the approximate solution.

60 APPLIED LIFE SCIENCES↗
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 OSTI2020

Photoinduced phase separation in the lead halides is a polaronic effect

We present a perspective on recent observations of the photoinduced phase separation of halides in multi-component lead-halide perovskites. The spontaneous phase separation of an initial homogeneous solid solution under steady-state illumination conditions is found experimentally to be reversible, stochastic, weakly dependent on morphology, yet strongly dependent on composition and thermodynamic state. Regions enriched in a specific halide species that form upon phase separation are self-limiting in size, pinned to specific compositions, and grow in number in proportion to the steady-state carrier concentration until saturation. These empirical observations of robustness rule out explanations based on specific defect structures and point to the local modulation of an existing miscibility phase transition in the presence of excess charge carriers. A model for rationalizing existing observations based on the coupling between composition, strain, and charge density fluctuations through the formation of polarons is reviewed.

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

Relative permeability as a stationary process: Energy fluctuations in immiscible displacement

Relative permeability is commonly used to model immiscible fluid flow through porous materials. In this work, we derive the relative permeability relationship from conservation of energy, assuming that the system to be non-ergodic at large length scales and relying on averaging in both space and time to homogenize the behavior. Explicit criteria are obtained to define stationary conditions: (1) there can be no net change for extensive measures of the system state over the time averaging interval; (2) the net energy inputs into the system are zero, meaning that the net rate of work done on the system must balance with the heat removed; and (3) there is no net work performed due to the contribution of internal energy fluctuations. Results are then evaluated based on direct numerical simulation. Dynamic connectivity is observed during steady-state flow, which is quantitatively assessed based the Euler characteristic. We show that even during steady-state flow at low capillary number ([Formula: see text]), typical flow processes will explore multiple connectivity states. The residence time for each connectivity state is captured based on the time-and-space average. The distribution for energy fluctuations is shown to be multi-modal and non-Gaussian when terms are considered independently. However, we demonstrate that their sum is zero. Given an appropriate choice of the thermodynamic driving force, we show that the conventional relative permeability relationship is sufficient to model the energy dissipation in systems with complex pore-scale dynamics that routinely alter the structure of fluid connected pathways.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2021

A Novel Multiphysics Multiscale Multiporosity Shale Gas Transport Model for Geomechanics/Flow Coupling in Steady and Transient States

Summary A novel multiphysics multiscale multiporosity shale gas transport (M3ST) model was developed to investigate shale gas transport in both transient and steady states. The microscale model component contains a kerogen domain and an inorganic matrix domain, and each domain has its own geomechanical and gas transport properties. Permeabilities of various shale cores were measured in the laboratory using a pulse decay permeameter (PDP) with different pore pressure and confining stress combinations. The PDP-measured apparent permeability as a function of pore pressure under two effective stresses was fitted using the microscale M3ST model component based on nonlinear least squares fitting (NLSF), and the fitted model parameters were able to provide accurate model predictions for another effective stress. The parameters and petrophysical properties determined in the steady state were then used in the transient-state, continuum-scale M3ST model component, which performed history matching of the evolutions of the upstream and downstream gas pressures. In addition, a double-exponential empirical model was developed as a powerful alternative to the M3ST model to fit laboratory-measured apparent permeability under various effective stresses and pore pressures. The developed M3ST model and the research findings in this study provided critical insights into the role of the multiphysics mechanisms, including geomechanics, fluid dynamics and transport, and the Klinkenberg effect on shale gas transport across different spatial scales in both steady and transient states.

Engineering↗
DOE OSTI2020

Enabling Seamless Integration of EV Charging Infrastructure with Weak Electric Grids

Widespread adoption of Electric Vehicles (EVs) for light, medium, and heavy-duty applications has gained significant interest. Based on ownership and user preferences, light duty and local delivery medium-duty EVs are typically connected to grid distribution networks whereas larger medium duty and heavy-duty EVs are typically connected to distribution or even sub-transmission networks. Challenges with at-scale EVs and charging infrastructure supporting them include interoperability, distribution network upgrades, surge in demand charges, power quality issues, voltage stability, etc. Based on real-world data, steady-state and dynamic assessments of light-duty EV charging with a focus on DC Fast Charging (DCFC) for weak distribution grids is presented. Three main areas of contributions in this paper include - steady state and dynamic power quality measurements and assessments using real world data, transient assessments of light duty EV integration with weak distribution networks (IEEE 13 node feeder system) under grid-connected and microgrid modes, and formulation of an intelligent charging decision algorithm to enable an impact-minimal charging.

actionable intelligence↗
DOE OSTI2024

Data-Driven Digital Twin for Reliability Assessment of DC/DC Buck Converter

In commercial applications, the operation of DC/DC converters significantly impacts overall system performance and long-term reliability. This study introduces a data-driven digital twin (DT) approach for estimating critical degradation parameters of DC/DC BUCK converter under steady-state condition. Initially, a circuit-level MATLAB/Simulink digital model (DM C ) is refined against a hardware prototype’s switching model dataset using offline particle swarm optimization. The optimized digital model’s steady-state response is then verified with its average model response while varying the duty and load. Subsequently, degradation profiles are imposed on the inductor, capacitor, MOSFET in the DMC. A large dataset is generated from this model, allowing training, validation, and testing of machine learning (ML) models for component health regression tasks. The proposed method employs random forest ML models, achieving impressive regression results with a squared R value as high as 0.99978 and a root mean square error of 4.2× 10 –6 . The method is further validated on a medium power level DC/DC BUCK prototype with varying load conditions, and includes the analysis of MOSFET’s on-resistance under degradation conditions. This data-driven DT method shows promise for identifying parasitic degradation and ohmic loss parameters, enhancing converter reliability assessments in a non-invasive, generalized, and computationally efficient manner.

14 SOLAR ENERGY↗
DOE OSTI2021

HyPerPy (Hydrogen Extraction and Parabolic Trough Plant Performance Models using Python) [SWR-21-53]

The HyPerPy package consists of three scripts: 1) Model for Hydrogen Tracking in Parabolic Trough Power Plants, 2) Model for Receiver Performance, and 3) Model for Hydrogen Extraction Process. The power plant model (1) tracks hydrogen generation and transport within the circulating heat transfer fluid (HTF) of the power plant. This script is a transient, initial value simulation, in which the hydrogen concentration in the circulating HTF is 0 moles per cubic meter everywhere at time 0 seconds. Hydrogen concentration is calculated at discrete locations within the circulating HTF with 4-second time resolution. During each time step, the change in hydrogen concentration due to hydrogen generation and permeation is calculated for each location according to the local HTF temperature, vessel or piping properties, hydrogen concentration and partial pressure. This model predicts hourly hydrogen concentrations for typical operating days in the spring, summer, fall, and winter seasons. Data is used to create hourly mappings of hydrogen concentrations for a typical operating year. The receiver model (2) uses a 1-hour time step to estimate getter loading and annulus pressure. For each time step, the model uses HTF and ambient temperature data to estimate absorber tube, bellows, and getter temperatures for the time step. In addition, the model uses hourly HTF hydrogen concentrations that are generated by the power plant model (1), and annulus hydrogen pressure from the previous time step. With these data, the model calculates the moles of hydrogen permeating across the absorber tube and bellows during the time step. The net change in moles hydrogen is added or subtracted from the getter loading for the previous time step, and the hydrogen pressure is re-calculated based on getter loading and temperature. The model uses this algorithm to simulate hydrogen permeation and loading 24 hours per day, 365 days per year using seasonal temperature data. The model repeats these calculations for 25 years to create a mapping of receiver getter loading and annulus hydrogen pressure for four seasons of each year. The model for hydrogen extraction (3) estimates hydrogen extraction rates for a specific separation module configuration. The rate depends primarily on membrane area, vacuum pump performance, headspace gas flowrate to the membrane, and headspace gas hydrogen partial pressure. This model has two versions. The steady-state version predicts hydrogen extraction rates when the module is operating in separation mode. The dynamic version predicts hydrogen transfer through the membrane when the module is operating in sensor mode. The steady-state version is used with the plant model (1) to predict hydrogen partial pressures in the power plant when the extraction process is operating.

Glatzmaier, Gregory↗