Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Load flow”

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 163 records · Page 9

Investigation of precooling unit design options in hydrogen refueling station for heavy-duty fuel-cell electric vehicles

Precooling gaseous hydrogen fuel to a cold temperature before refueling a heavy-duty (HD) hydrogen fuel cell electric vehicle (HFCEV) is essential to avoid overheating the vehicle tank, as well as achieving a high state of charge (SOC). Because a large volume of hydrogen is dispensed during each fill, the need for a shorter fill time amplifies the need to precool each load for refueling a HFCEV. Thus, the design and operation of a precooling unit (PCU), as well as the associated capital and operating costs, plays a pivotal role in any plans for heavy-duty hydrogen refueling stations. Here, in this paper, we present a thermodynamic and technoeconomic analysis of a PCU in a gaseous hydrogen refueling station (HRS) for HD HFCEVs. By employing Argonne National Laboratory's hydrogen station cost optimization and performance evaluation (H2SCOPE) model, the refueling of 50 kg of hydrogen on-board type IV tank at ambient temperatures of 15–45 °C and varying fill rates is simulated. The required degree of precooling temperature to obtain either 100% or the maximum possible SOC is obtained from the simulation. Additionally, the simulation results demonstrate that the average flow rate of hydrogen is approximately 40% lower than the maximum flow rate during a typical fill; which motivates further evaluation of the instantaneous hydrogen mass flow rate profile and suggests the scope of improving precooling unit design. Accordingly, a hybrid strategy of precooling hydrogen has been proposed to address the cooling load by sizing the refrigeration unit for the average flow rate of hydrogen, while supplementing the above average peak hydrogen flow cooling load through thermal buffering. The combined technique enables the downsizing of the original PCU capacity by 25–40% and demonstrates a potential cost reduction of the PCU by approximately 30%, which translates to an installed cost reduction of ∼$125,000 per dispenser.

25 ENERGY STORAGE↗

An Efficient Three-Dimensional CFD-Based Numerical Wave Tank for a Wave Energy Converter in Extreme Irregular Waves

A numerical wave tank approach for computational fluid dynamics (CFD) modelling of an extreme irregular seastate is presented. The technique couples a potential flow solution with a CFD solver for more efficient numerical predictions. This method has recently become attractive both for the research community and the industry working with offshore structures. The model is used to determine the response of a submerged pressure differential wave energy converter (WEC) in a fully nonlinear irregular wave condition using the high-fidelity CFD code, STAR-CCM+. Potential flow based numerical models are commonly used to predict motions and performance of wave energy converters. Wave kinematics can deviate from potential flow predictions for extreme wave conditions; the excitation loads on an absorber can also be increasingly influenced by viscous effects, not predicted by potential flow engineering level models. In these extreme conditions, a Reynolds-averaged Navier-Stokes CFD model can better predict motions and loads for a WEC. Long time series with varying random seed numbers can be used to identify singular extreme wave events from a stochastic irregular sea state. This approach simulates a more realistic wave series for a given sea state than a regular wave or a focused wave. However, it is computationally infeasible to run these long time series for three-dimensional (3D) CFD simulations. In this work, two-dimensional (2D) CFD simulations with a long domain allow the full development of an extreme nonlinear wave condition. The results are used to identify extreme events from a 50-year storm condition for the PacWave site off the coast of Oregon. A relatively short time window including this extreme event is then mapped to a 3D simulation using a user defined wave methodology. Convergence studies for domain length, wave forcing lengths, and time before the extreme event were conducted.

CFD↗

Real time computations of cryogenic He properties

The Fermilab PIP-II (proton improvement plan - II) project is being constructed at Fermilab to deliver $800\,MeV$ protons of $>1\,MW$ beam power to replace the present LINAC and provide protons to the remainder of the existing accelerator complex. The new LINAC consists of a warm front end, 23 superconducting RF cryomodules, and a beam transfer line to the existing complex. The cryomodules (CMs) are to be tested at Fermilab's CryoModule Test Facility (CMTF).An important measurement in cryogenic testing is the heat load of each CM. Traditionally, at Fermilab, these measurements were made collecting archived data offline and analyzing it. The new control system for PIP-II is being developed with the EPICS (Experimental Physics and Industrial Control System) framework, which allows us to compute the heat load in real time using the HePak library.We are exploring other $He$ properties, such as flow, where flow meters are not available, which can also be calculated in real time and fed back to the cryogenics engineers.This paper details the real time heat load calculation and $He$ flow software developed for CM testing at CMTF, as well as the first results from the prototype HB650 CM. Future plans for 2-phase $LHe$ flow will also be outlined.

Hanlet, Pierrick [Fermilab]↗

Leading Edge Microspoilers for Load Dumping in Wind Turbine Rotors

As rotor designs have grown in size and flexibility, new design challenges have arisen that could benefit from novel control strategies. Leading edge microspoilers (LEMS) have been shown to effectively reduce forces on airfoils by inducing stall at lower angles of attack. This paper investigates how LEMS can be used to reduce blade loading on large, flexible, downwind rotors during design-limiting design load cases. During shutdown in turbulent flows near cutout speed, LEMS could reduce the tip deflection spikes toward the tower by over 50% and keep flapwise root bending moments below normal operating conditions, thus eliminating these as design-driving load cases. In power producing load cases, deploying LEMS for a short period of time during a load spike can be effective to reduce toward-tower tip deflections without significantly impacting power produced, assuming these load spikes can be identified by a controller. Leveraging this technology could allow for lower specific power machines to be designed.

17 WIND ENERGY↗

Ethanol-fueled metal supported solid oxide fuel cells with a high entropy alloy internal reforming catalyst

High-performance metal supported solid oxide fuel cells (MS-SOFC) with an integrated high entropy alloy (HEA) internal reforming catalyst (IRC) are demonstrated for transportation applications using ethanol and methanol as fuels. Addition of the HEA IRC dramatically improves cell performance and stability when using ethanol/water blend fuel. Absence of carbon deposition predicted by thermodynamic calculations is confirmed by Raman spectroscopy analysis of posttest anodes. Optimal catalyst processing (deposition technique, loading, firing temperature) and cell operation conditions (flow rates, temperature, fuel compositions) are explored. Infiltrated HEA reforming catalyst provides a highly porous structure and low catalyst loading (6 mg cm –2 ). The designed structure and catalysts achieve small mass transport resistances in the fuel electrode (26.2 s m –1 ) and oxygen electrode (41.6 s m –1 ). The best ethanol concentration (60:40 v% ethanol: water) provides 0.83 W cm –1 at 700 °C, without carbon deposition. The ethanol-fueled MS-SOFC is operated for 500 h, including five thermal cycles. As a result, cell evolution is similar to that reported previously for hydrogen fuel; nickel aggregation and chromia deposition were the major observed changes, and carbon formation can be avoided even after long-term operation.

30 DIRECT ENERGY CONVERSION↗

One Dimensional Penstock Flow Models for Hydropower Digital Twin

This report details the inelastic, unsteady, one-dimensional flow equations that are used to model flow through a penstock in a hydropower facility. The flow model is coupled to a mechanistic turbine model to compute the rotation rate of the turbine, which is dependent on the flow rate, gravitational head, and various parameters of the machine, such as inlet vane angle. Results are given for an example hydropower facility in which the turbine is subjected to a specified time-varying load, and the results on the flow rate are shown. Future applications for the flow model are also briefly reviewed.

13 HYDRO ENERGY↗

Studying Wind Loading on CSP Collectors to Improve Performance and Reliability

Electricity generation through Concentrating Solar-thermal Power (CSP) adds the advantage of thermal energy storage and heat production. CSP collector costs constitute nearly one-third of total plant costs. Managing wind loading, especially dynamic wind loading caused by the turbulent wind flow, is a significant design challenge. Traditional collector designs rely on wind tunnel experiments and numerical simulations, which do not fully capture dynamic effects. Here, we provide an overview of NREL's activities on wind loading on CSP collectors, aiming to improve their reliability and cost-efficiency. We present insights from field measurements at operational parabolic trough and power-tower CSP plants. Over two years, we measured turbulent wind fields and resulting structural loads in a parabolic trough plant, revealing wind, turbulence, and collector field interactions. Similar measurements are ongoing at the Crescent Dunes heliostat field. These campaigns show how atmospheric turbulence, wind direction, and collector orientation affect dynamic wind loading. Upwind collector structures can generate turbulent structures, causing fluctuating loads at downstream collectors, impacting fatigue lifetime and optical efficiency. Preliminary results indicate that interior collectors experience higher turning moments than those at the edges, suggesting higher drive wear. At Crescent Dunes, we use instrumented heliostats at the field's edge and interior to study the translation of turbulent wind to dynamic structural loads. The campaign is still ongoing, and we will present first findings. Another focus of our work is tying the experimental wind loading findings to optical performance of the collectors, using simulation tools. For example, we derive torsional errors of the parabolic troughs - the angular offset from the sun position along a collector row. This error is impacted by wind-induced mirror deformations. Future work will focus on continuing the measurement efforts and combining them with computational fluid dynamic simulations. By providing detailed field measurements and validated models, our efforts aim to enhance understanding of wind loading on CSP collectors, improving their structural integrity and optical performance.

CSP↗

Development of Real-Time System Identification to Detect Abnormal Operations in a Gas Turbine Cycle

Here, we present a novel online system identification methodology for monitoring the performance of power systems. This methodology was demonstrated in a gas turbine recuperated power plant designed for a hybrid configuration. A 120-kW Garrett microturbine modified to test dynamic control strategies for hybrid power systems designed at the National Energy Technology Laboratory (NETL) was used to implement and validate this online system identification methodology. The main component of this methodology consists of an empirical transfer function model implemented in parallel to the turbine speed operation and the fuel control valve, which can monitor the process response of the gas turbine system while it is operating. During fully closed-loop operations or automated control, the output of the controller, fuel valve position, and the turbine speed measurements were fed for a given period of time to a recursive algorithm that determined the transfer function parameters during the nominal condition. After the new parameters were calculated, they were fed into the transfer function model for online prediction. The turbine speed measurement was compared against the transfer function prediction, and a control logic was implemented to capture when the system operated at nominal or abnormal conditions. To validate the ability to detect abnormal conditions during dynamic operations, drifting in the performance of the gas turbine system was evaluated. A leak in the turbomachinery working fluid was emulated by bleeding 10% of the airflow from the compressor discharge to the atmosphere, and electrical load steps were performed before and after the leak. This tool could detect the leak 7 s after it had occurred, which accounted for a fuel flow increase of approximately 15.8% to maintain the same load and constant turbine speed operations.

algorithms↗

Voltage Optimization

Voltage optimization refers to a volt-var optimization technique which was originally designed to minimize energy consumption and improve end-use efficiency on the distribution system. This reduces source voltage at a substation, which lowers generation demand, reduces system losses, and improves system stability. This has been the conventional implementation of voltage optimization: focusing on regulation of voltages throughout the distribution system via coordinated adjustments of load tap changers, line-voltage regulators, switched shunt compensation, and the targeted placement of new shunt compensation. These methods have received increased focus following a growth in distributed-energy resources (DERs) interconnecting at the distribution level, given their developing capability to participate in voltage regulation. This is paired with the intermittency challenges they pose, load imbalances, and voltage fluctuations. Voltage optimization has long been an area of extensive research within the distribution system, with various strategies and philosophies demonstrating a range of effects. The conventional focus on the distribution system owes to the typical operation of transmission systems, which assumes firm local resources and flows strictly from generation resources to loads. The strategies employed in distribution have also been optimized for this typically radial configuration and operation of the distribution system.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

Measurement of Wind Loading on Heliostats at the Crescent Dunes Power Plant: An Overview

The cost of solar collectors constitutes almost one third of the total cost of a CSP plant. One of the ongoing challenges in the design of these collectors is wind loading on mirrors, support structures, and drives. A particular challenge is dynamic wind loading, caused by the turbulent wind flow. To date, the design of solar collector structures has relied on wind tunnel experiments and numerical simulations that do not entirely capture the dynamic effects observed at scale. The CSP industry has shown increased interest in validating the idealized assumptions with measurements obtained in operational settings to improve wind load assumptions and increase reliability and cost-efficiency of the collector design. Further, performance models need realistic assumptions about wind loading and its impact on optical performance. In a parabolic trough field campaign, NREL successfully collected a wealth of long-term, high-resolution wind and loads data [2], that are publicly available and that can be used for the above-mentioned purposes. Heliostats are impacted differently by wind than parabolic troughs, due to their different shape, size, and field layout. To study the impact of wind and turbulence on heliostats, we initiated another field campaign in an operational power- tower plant, Crescent Dunes, in Nevada, USA. In this work, we present an overview of the measurements, first results, and potential implications of wind driven loads on Heliostats.

concentrated solar power↗

Numerical Investigation of Two-Stage Vapor Compression System with Simultaneous Vapor and Liquid Injection

Two-stage vapor compression systems can be advantageous over single-stage systems by providing improved system performance, lower discharge temperature and reduced throttling losses. These systems employ various intermediate configurations, such as liquid injection and vapor injection. This paper presents a configuration for a R1234ze(E) two-stage air-conditioning system using a turbo compressor with two injection ports, one for vapor injection and the second for liquid injection. The liquid injection is used to cool the motor and electronics. A component-based representation and solution approach was used to simulate the two-stage compression system with simultaneous vapor and liquid injection at steady state. The turbo compressor was represented using a customized performance map. The condenser and the evaporator were modeled using finite-volume approach. A parametric study was conducted to assess the impact of the following three variables on the system performance: vapor injection ratio, condenser air flow rate, and discharge pipe pressure drop. The simulation results show that as the vapor injection ratio increased, the system performance undergoes a tradeoff between an enhanced subcooling effect (and thus enhanced unit refrigeration capacity) and a decreased suction mass flow rate. Maximum COP occurs when the vapor injection ratio was 0.1. The results also show that as the condenser air flow rate increased, both the capacity and power consumption (including fan power) increased monotonically, and COP increased first and then decreased. At 75% load, the COP improvement at the optimum flow rate was marginally less than 0.5%. Lastly, higher discharge pipe pressure drop increased the discharge pressure. It showed very small effect on the overall system performance at the condition selected for the current study.

component-based steady-state simulation↗

Oregon Direct Use Geothermal Well Data

This excel file provides direct use geothermal data for various wells located in Oregon. Data includes well name, location, heating application, temperature and flow of fluid, capacity, energy use, and load factor. The data is from the State of Oregon Department of Geology and Mineral Industries (linked below), where additional data downloads and metadata may be found.

15 GEOTHERMAL ENERGY↗

Off-design operation and performance of pumped thermal energy storage

In this article, we describe off-design models and control strategies for a Pumped Thermal Energy Storage (PTES) system that uses liquid thermal energy storage: specifically molten salt for hot storage and methanol for cold storage. Off-design conditions arise when load-following, or due to variations in storage tank temperatures or ambient temperatures. We propose a control strategy that uses inventory control to manage the mass flow rate in the thermodynamic cycles, which facilitates load following. We also propose a control strategy for the storage fluid mass flow rates, which are varied to ensure the molten salt is maintained at its design temperature. This maximizes efficiency and minimizes problems with salt freezing or degradation. The cold storage fluid mass flow rate is varied so that the cold tanks have the same state-of-charge as the hot tanks. This leads to variations in cold fluid temperature, but these variations are shown to be acceptably small (e.g. 7.5% increase), and this control method is shown to be simpler and more efficient than an alternative strategy where tanks become unbalanced. The ambient temperature and storage tank temperatures are moved ±50 °C from the design values and the impact on power, duration, and tank temperatures is quantified. Results demonstrate that the proposed control strategy is stable and self-correcting - that is, storage temperatures converge on stable values after two-to-three charge-discharge cycles. When inputs return to design values, the system returns to its design point after two charge-discharge cycles. We also demonstrate that inventory control enables delivery of the target power output even when off-design conditions exist that would normally reduce the power output.

25 ENERGY STORAGE↗

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗

High Solids Performance Testing in a Scaled TSCR System

The Tank Side Cesium Removal (TSCR) project is a technology demonstration that will pretreat Hanford tank waste supernatant in support of the Direct Feed Low-Activity Waste (DFLAW) mission. The TSCR system employs two key separation technologies: dead-end filtration (DEF) and ion exchange (IX) using crystalline silicotitanate (CST) media. DEF will be used to remove undissolved solids from tank waste to protect the functionality of the IX columns and the IX system will remove Cs-137 from tank waste. The separation technologies (DEF and IX) used in TSCR are technically mature and have also been successfully deployed at the Savannah River site in a similar facility known as the Tank Closure Cesium Removal (TCCR) system. While testing with simulants and real waste has been successfully performed under conditions expected during the initial operation of TSCR, test data is absent for assessing off normal high solids loading that may be in the TSCR waste feed. Normal TSCR treatment operations are expected to handle wastes with solids content on the order of 200 ppm, and off normal solids loading could be much larger than the nominal level. The testing program described in this report was conducted to understand the consequence of operating the TSCR system at elevated solids loadings up to the high-solids limit of 15,000 ppm [i.e., 1.5-wt%] identified in the TSCR design basis. Although the system is not required to make throughput above the nominal solids loading, the testing was intended to provide important information related to potential off normal operations. At off normal levels near the high-solids limit, there are potential implications for TSCR performance in the areas of throughput, DEF pressure drop, filter backflush frequency, and IX column pressure drop. In addition, intrusion of solids into the IX column was postulated to impact the Cs-137 loading behavior by promoting channeling or flow maldistribution in the column; since the magnitude of the postulated effect was unknown, assessing it was also of interest. The testing was performed using representative waste simulants and a prototypic, integrated TSCR system designed and assembled specifically to conduct the high solids performance assessment. Overall, the scaled TSCR testing demonstrated that full-scale unit operations can succeed in fulfilling their processing objectives in the presence of solids up to 3,000 ppm, but there are potential performance challenges to filtration operations at solids loadings as low as ~500 ppm. The severity of the challenge is likely to be dependent on the type and size distribution of solids, of which the current testing only examined a single type and size distribution. To provide some flexibility for future full-scale operations, the results of the testing suggest two possible risk reduction strategies that can be implemented without any changes in TSCR design or configuration. One option would be to enact an administrative limit on the solids loading to protect TSCR from feeds that are likely to require a high DEF swap frequency. Another option is to permit operation of the DEFs at differential pressures greater than 2 psid before swapping filters. The selection of a higher differential pressure target is not anticipated to adversely impact DEF backflushing efficacy and would reduce both swap frequency and the amount of waste sent to AP-108.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

MHD stability and disruptions in the SPARC tokamak

SPARC is being designed to operate with a normalized beta of $β_N$ = 1.0, a normalized density of $n_G$ = 0.37 and a safety factor of $q_{95}$ ≈ 3.4, providing a comfortable margin to their respective disruption limits. Further, a low beta poloidal $β_p$ = 0.19 at the safety factor $\textit{q}$ = 2 surface reduces the drive for neoclassical tearing modes, which together with a frozen-in classically stable current profile might allow access to a robustly tearing-free operating space. Although the inherent stability is expected to reduce the frequency of disruptions, the disruption loading is comparable to and in some cases higher than that of ITER. The machine is being designed to withstand the predicted unmitigated axisymmetric halo current forces up to 50 MN and similarly large loads from eddy currents forced to flow poloidally in the vacuum vessel. Runaway electron (RE) simulations using GO+CODE show high flattop-to-RE current conversions in the absence of seed losses, although NIMROD modelling predicts losses of ~80 %; self-consistent modelling is ongoing. A passive RE mitigation coil designed to drive stochastic RE losses is being considered and COMSOL modelling predicts peak normalized fields at the plasma of order 10 –2 that rises linearly with a change in the plasma current. Massive material injection is planned to reduce the disruption loading. A data-driven approach to predict an oncoming disruption and trigger mitigation is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evaluation of Hopper Flow Prediction based on Cubical Triaxial Tester Data and FEM

A cubical triaxial tester (CTT) laboratory device was used to collect key bulk flow characteristics of Crumbles® 6mm rotary sheared woody biomass and Crumbles® 2mm rotary sheard corn stover. The flow properties and a CAD model of a hopper and conveyor system to predict the consistency of flow of material from the hopper. Prediction included the likeliness of jamming, bridging, or ratholing. Each material was then loaded into the hopper and the actual bulk flow was measured. Measured flow closely matched the predicted flow.

09 BIOMASS FUELS↗