Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLOW RATE”

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 73 records · Page 4

Hanford Waste Treatment Plant Effluent Management Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant Effluent Management Facility (EMF) stack monitor location was qualified using a combination of scale model stacks to mitigate the risk of identifying that the sampling location does not meet the qualification criteria on the full-scale stack. The LV-S1 scale model stack was used as a baseline, augmented by the LB-S1 and LV-S2 scale model stacks to address the Direct Feed Low Activity Waste Effluent Management Facility Vessel Vent Process (DVP) injection into the main Active Confinement Ventilation (ACV) system duct. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling locations were geometrically similar to the actual stack, and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. The LV-S1, LB-S1, and LV-S2 scale model stack tests have met the criteria of the ANSI/HPS N13.1-1999 standard to demonstrate the stack sampling locations are well mixed. Verification tests of the EMF stack were performed at normal operating conditions. The minimum 1/6 DV value and the maximum 6 DV value from the scale model testing determine the range of stack flow rates for which the full-scale stack may be operated while remaining in compliance with the stack verification criterion. A practical range for the full-scale stack qualification uses the average DV through 6 DV from the scale model tests to compute the corresponding flow rates. Table S1 lists the operating flow rate along with the average and maximum qualified stack flow rate based on the LV-S1 scale model DV values. The operating flow is below the maximum qualified stack flow, which means that the scale model test results are acceptable for stack qualification. The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be ≤20°. Second, the velocity uniformity at the full-scale stack must be ≤20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack test at the EMF. Flow angle results were <5°; all flow angle results were within the ≤20° criterion. The velocity uniformity results for each test condition ranged between 2.2% COV and 4.3% COV, all of which were within the range of the target % COV values from the scale model tests on the LV-S1, LB-S1, and LV S2 scale models. Based on these stack verification test results, the EMF filtered exhaust stack sampling location meets the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all planned fan operating configurations. This includes each combination of ACV fans with DVP exhausters. Further changes to the system configuration or operating conditions that are outside the qualified flow rates described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

In-situ Observations of Directed Energy Deposition Additive Manufacturing Using High-Speed X-ray Imaging

In laser-based directed energy deposition (DED) additive manufacturing, interactions among the laser beam, particle flow, and melt pool influence the properties of the solidified final part. Two separate DED systems, one with high powder flow rates to represent industrial-scale DED processing and the other with low powder flow rates for individual particle tracking, were synchronized with the high-speed imaging setup at the Advanced Photon Source in Argonne National Laboratory. In-situ x-ray imaging of the DED process using both systems highlighted the influence of powder flow rates. Increased powder flow rates resulted in less laser absorption into the melt pool, leading to a transition from a keyhole mode to a melt pool without a keyhole but with surface fluctuations due to powder flow. Increased velocities of particles during powder flow resulted in a decrease in particle melting times and a greater propensity for porosity formation. Altogether, better understanding of the interactions that occur during various scales of the DED process will enable flexibility, control, and new materials development in DED-based additive manufacturing.

36 MATERIALS SCIENCE↗

Quantitatively Monitoring Bubble-Flow at a Seep Site Offshore Oregon: Field Trials and Methodological Advances for Parallel Optical and Hydroacoustical Measurements

Two lander-based devices, the Bubble-Box and GasQuant-II, were used to investigate the spatial and temporal variability and total gas flow rates of a seep area offshore Oregon, United States. The Bubble-Box is a stereo camera–equipped lander that records bubbles inside a rising corridor with 80 Hz, allowing for automated image analyses of bubble size distributions and rising speeds. GasQuant is a hydroacoustic lander using a horizontally oriented multibeam swath to record the backscatter intensity of bubble streams passing the swath plain. The experimental set up at the Astoria Canyon site at a water depth of about 500 m aimed at calibrating the hydroacoustic GasQuant data with the visual Bubble-Box data for a spatial and temporal flow rate quantification of the site. For about 90 h in total, both systems were deployed simultaneously and pressure and temperature data were recorded using a CTD as well. Detailed image analyses show a Gaussian-like bubble size distribution of bubbles with a radius of 0.6–6 mm (mean 2.5 mm, std. dev. 0.25 mm); this is very similar to other measurements reported in the literature. Rising speeds ranged from 15 to 37 cm/s between 1- and 5-mm bubble sizes and are thus, in parts, slightly faster than reported elsewhere. Bubble sizes and calculated flow rates are rather constant over time at the two monitored bubble streams. Flow rates of these individual bubble streams are in the range of 544–1,278 mm 3 /s. One Bubble-Box data set was used to calibrate the acoustic backscatter response of the GasQuant data, enabling us to calculate a flow rate of the ensonified seep area (~1,700 m 2 ) that ranged from 4.98 to 8.33 L/min (5.38 × 10 6 to 9.01 × 10 6 CH 4 mol/year). Such flow rates are common for seep areas of similar size, and as such, this location is classified as a normally active seep area. For deriving these acoustically based flow rates, the detailed data pre-processing considered echogram gridding methods of the swath data and bubble responses at the respective water depth. The described method uses the inverse gas flow quantification approach and gives an in-depth example of the benefits of using acoustic and optical methods in tandem.

54 ENVIRONMENTAL SCIENCES↗

Testing of a 40-kWth Counterflow Particle-Supercritical Carbon Dioxide Narrow-Channel, Fluidized Bed Heat Exchanger

Particle-based primary heat exchangers (HXs) must deliver sCO2 fluid temperatures above 700°C to couple particle-based concentrating solar receivers and thermal energy storage (TES) sub-systems with efficient sCO2 power cycles. Particle-sCO2 HX designs have struggled to meet DOE cost targets (≤ $150/kWth) due to the amount of expensive nickel alloys necessary for manufacturing full-scale, particle-sCO2 HXs. Our team has demonstrated that mild bubbling fluidization of falling particles in a counterflow narrow-channel fluidized bed can reduce required HX surface area and thus, costs by increasing particle-wall heat transfer coefficients hT,w > 800 W m-2 K-1. This paper reports on the fabrication and testing of a stainless steel, particle-sCO2 HX with 12 fluidized-bed channels approximately 10.5 mm deep spaced between diffusion-bonded, micro-channel sCO2 plates. The HX with a core length of ≈0.56 m is fed with CARBOBEAD HSP particles through a short, fluidized freeboard zone just above the core. Testing to date in the National Solar Thermal Test Facility (NSTTF) at Sandia National Laboratories has shown that parallel bed fluidization maintains uniform particle inventory across the instrumented channels. Heat transfer thermal duty between the particle and sCO2 flows exceeds 30 kWth with sCO2 inlet temperatures of 200ºC and particle inlet temperatures up to 440ºC and mass flow rates of 0.2 kg s-1 fluidized by counterflowing gas flow rates of 0.005 kg s-1. Tests at higher particle and sCO2 inlet temperatures (600ºC and 400ºC respectively) are targeted to achieve > 40 kWth with model-predicted overall heat transfer coefficients U > 400 W m-2 K-1.

14 SOLAR ENERGY↗

Experimental Evaluation of Refrigerant Leak Characteristics for Different HVAC&R Equipment Types–Phase 2

The US Department of Energy (DOE) Building Technology Office (BTO) and the Air-conditioning, Heating & Refrigeration Technology Institute (AHRTI) collaborated to sponsor an experimental study of refrigerant leak characteristics at Oak Ridge National Laboratory (ORNL); AHRTI project 9012. The project has been conducted in two phases with objectives to conduct refrigerant leak tests on several heating, ventilating, air-conditioning and refrigeration (HVAC&R) systems under operating conditions representative of actual applications to document the pressure decay rate and the mass flow rate of leaked refrigerant as a function of time. Phase 1 of the project ran from late 2017 to late 2018 with refrigerant leaks imposed on five systems and the collection of data to meet the project objectives. Systems chosen covered three air-conditioning (AC) applications (packaged terminal AC or PTAC, 3-ton residential split system AC, and 5-ton packaged rooftop AC) and two refrigeration applications (split system unit cooler, and standalone display case). AHRTI members donated the test systems used for the project. The systems were tested at two different leak rates (catastrophic high rate simulating full line break and a lower flow rate), two leak locations (high-pressure side and low-pressure side), and two operating conditions (compressor ON or OFF). Phase 1 results were summarized in a project report by Baxter, et al. (2019)1. The leak rates imposed in Phase 1 led to much more rapid charge releases than the more typical much slower rates seen in practice that would take days, weeks, or months to empty a system of its entire charge. For Phase 2, the project sponsors asked that the “higher” leak rate target a total charge release time of about 4 min. This was intended to represent a “reasonable” worst case release scenario to match the refrigerant leak rate assumption used in the International Electrotechnical Commission (IEC) standard IEC 60335-2-40 (2016 version [IEC 2016]. A charge release time of about 20–40 min (leak release rate 5 to 10 times lower than higher rate) was targeted for the “lower” leak rate. This report provides a summary of the Phase 2 effort. Two additional HVAC systems were selected for testing in this phase. First was a multi-split air-conditioning (MS-AC) system with one outdoor section and two indoor units and a nameplate refrigerant charge of 1.75 kg (3.86 lb.) of R-410A. The MS-AC was tested in cooling only or AC operation. Second was a split system heat pump coupled to two different representative supply duct systems (HP-duct) and a nameplate charge of 3.83 kg (8.44 lb.) of R-410A. This system was tested in both space heating and AC modes. As in Phase 1 of the project, all tests were conducted using R-410A.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Reduced Order Model to Predict Dispersion of Flammable Refrigerant into a Space

As the HVAC&R industry mobilizes to deploy more low-GWP refrigerants, relevant standards are being continually reviewed and updated. Those include the general safety standards ISO 5149 and ASHRAE 15, and the equipment standards IEC and UL. The standards systematically set the allowable maximum amount of refrigerant that should be used in different equipment types and different applications. To do so, they rely on predictions of how a leaked refrigerant mass will disperse into a space. Dispersion characteristics, such as total flammable volume and its residence time, determine the risk associated with the presence of the flammable refrigerant. The standards have included provisions for the use of flammable refrigerants for approximately two decades. They relied on limited analytical analyses and test cases in their development. Dispersion of a refrigerant into a space is complex. Computational fluid dynamics (CFD) are the most accurate in predicting a given problem. However, CFD is computationally expensive and requires specialized expertise and resources and is not suitable for use by standards development working group as prediction tool. This paper presents the development of a reduced order model (ROM) that predicts the key dispersion characteristics relevant to the dispersion of a leaked refrigerant into a space for any combination of input variables. The inputs are the refrigerant release height, the total released refrigerant mass and its release flow rate, the refrigerant molecular weight, the ventilation flow rate, the floor area and height of the space, recirculation air flow rate, and the tightness of the space. The outputs are histograms of volume fraction of the room in prescribed concentration bins and the total mass of the refrigerant in each bin normalized by the total refrigerant charge at 13 prescribed simulation time stamps between 1 and 900 seconds. The ROM is constructed from a set of CFD simulations with carefully chosen combinations of input parameters. The selection if done using a multidimensional sparse grid which is a generalization of the classical tensor approach but offers additional flexibility and thus can be more carefully tuned towards a specific model. The tuning is done to improve the accuracy, measured in the difference between the output values of the ROM and the CFD model, while minimizing the computational cost, measured in number of CFD simulations which is orders of magnitude more expensive than the processing the training data.

Edwards, Dean↗

Fluid Flow Observations of The Spray Near-Field Using High-Speed X-Ray Imaging

Sprays are important to may industrial processes, but the near-field region where spray formation takes place is difficult to visualize and characterize through visible light imaging because this region is typically optically dense. X-ray imaging is one technique that can penetrate this optically dense region and provide unique observations of the spray formation process. This study uses high-speed white beam X-ray imaging conducted at the Advanced Photon Source to capture spray formation dynamics from a canonical airblast atomizer. In this work, a range of momentum flux ratios are visualized by varying the gas flow rate while the liquid flow rate is held constant. The use of X-rays show overlapping flow structures including bags, ligaments, mushrooms, and webs. Air bubbles are captured within the liquid regions, and a very dynamic crown forms at the liquid needle exit at higher momentum flux ratios. Finally, at the highest momentum flux ratios, air bubbles and liquid droplets are observed migrating upstream in the center of the crown, and this motion is attributed to a toroidal recirculation cell at the nozzle exit.

47 OTHER INSTRUMENTATION↗

Experimental Investigation of Fuel Staging Effect on Modal Dynamics of Thermoacoustic Azimuthal Instabilities in a Multi-Nozzle Can Combustor

Abstract This paper analyzes the dynamics of unstable azimuthal thermoacoustic modes in a lean premixed combustor. Azimuthal modes can be decomposed into two counter rotating waves where they can either compete and potentially suppress one of them (spinning) or coexist (standing), depending on the operating conditions. This paper describes experimental results of the dynamical behaviors of these two waves. The experimental data were taken at different mass flow rates as well as different azimuthal fuel staging in a multi-nozzle can combustor. It is shown that at a low flow rate with uniform fuel distribution, the two waves have similar amplitudes, giving rise to a standing wave. However, the two amplitudes are slowly oscillating out of phase to each other, and the phase difference between the two waves also shows oscillatory behavior. For an intermediate flow rate, the dynamics show intermittency between standing and spinning waves, indicating that the system is bistable. In addition, the phase difference dramatically shifts when the mode switches between standing and spinning waves. For a high flow rate, the system stabilizes at a spinning wave most of the time. These experimental observations demonstrate that not only the amplitudes of two waves but also the phase difference plays an important role in the dynamics of azimuthal mode. For non-uniform azimuthal fuel staging, the modal dynamics exhibit only an oscillatory standing wave behavior regardless of the mass flow rate. Compared to the uniform fuel staging, however, the pressure magnitude is considerably reduced, which provides a potential strategy to mitigate and/or suppress the instabilities.

Kim, J.↗

Effect of Initial Water Saturation on The Performance of Fracturing Fluids With and Without Polyallylamine under Simulated EGS Conditions

Objectives/Scope: StimuFrac (US Patents 9,873,828 B2 and 9,447,315 B2), a CO2-reactive polymer aqueous solution [polyallylamine (PAA) 1wt% in water] combined with CO2, can be used as a less water-intensive fracturing fluid for enhanced geothermal systems (EGS). Our previous results show that in hot dry rock (HDR), PAA/CO2 fracturing fluids outperformed other fluids such as water, CO2, and CO2/water in generating large fractures with less fluid consumed. The objective of this work is to study the effect of initial water saturation of rock on the performance of StimuFrac fluid in ½ foot cubic rock samples and under representative EGS pressure/temperature conditions using cyclic and constant flow rate injection strategies. The fracturing results are compared with results using different fracturing fluids in terms of controlling fracture propagation rates, fracture hydraulic conductivity, breakdown pressures and volumes of fluids required. Methods/Procedures/Process: In all tests, water was initially injected into the rock to increase the water saturation before the fracturing processes to simulate actual geothermal reservoir conditions. For the cyclic injection, one complete cycle consisted of (1) a PAA slug (or water slug) injection followed by (2) CO2 injection to initiate the fracture. In the second step of the first cycle, the pressure of CO2 is increased until a maximum pressure is reached (fracture is initiated at this moment), and then the injection of CO2 is allowed for another 30 seconds to propagate the fracture. Then, the two-step cycle of PAA followed by CO2 injection (up to 2-4 mL/min) was continued. For the constant flow rate injection strategies, the initial water saturation was increased by injecting water at 1000 psi and 200°C for three days. After that, an initial slug of water, CO2, or PAA was injected and then fracturing was initiated and propagated by injecting CO2 at a constant flow rate. Applications/Significance/Novelty: The results of this study suggest that water saturation, especially near the wellbore region, will significantly affect the fracturing fluid transmission into the rock porous media by changing the relative permeability of CO2 or water, thus affecting the fracture initiation and propagation. In this study, fracturing with cyclic injection or constant flow rate injection methods were performed using three different kinds of fluids systems. These fluids are water, CO2, or CO2 with PAA. Splitting the rock samples in half after fracturing reveals that the fracture propagation is significantly limited under these high water saturated conditions compared to dry initial conditions: The fractures propagate less than 1/3 length of the distance from the wellbore to rock surface, and in some cases no fracture is generated. This may be caused by the fact that leak-off is dominating the fracturing process and the injected fluid flow rate is not high enough to overcome the leak-off even under high flow rate injection conditions. Additionally, CO2 could be leaking off into the wellbore annulus and this may be making it more difficult to generate pressure gradients away from the near-wellbore region.

Jian, Guoqing↗

A Design Method for Low-Pressure Venturi Nozzles

The purpose of this work is to provide empirical design models for low-pressure, subsonic Venturi nozzles. Experimentally validated simulations were used to determine the effect of nozzle geometry and operating conditions on the suction ratio (ratio of suction mass flow rate to motive mass flow rate) of low-pressure, subsonic Venturi nozzles, over a wide range of geometries and operating conditions, through a parametric study. The results of the parametric study were used to develop seven empirical models, each with a different range of applicability or calculating a different indicator of nozzle performance (i.e., suction ratio, momentum ratio, or dynamic pressure ratio), of the Venturi nozzles using a constrained multi-variable global optimization method. Of the seven empirical models, the best models were found to be those for low- (less than one) and high-suction ratios (greater than one), with mean absolute percentage errors of 5% and 18%, respectively. These empirical models provide a design tool for subsonic, low-pressure Venturi nozzles that is more than an order of magnitude more accurate than a governing equation approach or conventional flow head calculations. These newly-developed empirical models can be applied for initial nozzle design when precise suction ratios are required.

O’Hern, Hannah↗

Chemical looping air separation with Sr 0.8 Ca0.2Fe 0.9 Co 0.1 O 3-δ perovskite sorbent: Packed bed modeling, verification, and optimization

Chemical looping air separation (CLAS) represents a promising approach for efficient O 2 production from the air. This present study aims at optimizing the absorber/desorber operations and the separation process with extensive experimental validation. Specifically, a one-dimensional packed bed model was developed to investigate the CLAS operation with a Sr 0.8 Ca 0.2 Fe 0.9 Co 0.1 O 3-δ perovskite sorbent. The redox thermodynamics of perovskite sorbent was measured by TGA and then incorporated into a linear driving force model to describe the O 2 absorption and desorption rates. Both 4-step and 5-step air separation cycle configurations, with various cyclic structures, were performed in a subpilot-scale packed bed. The model predicted O2 purity and productivity were consistent with experimental results, supporting its accuracy and applicability. Parametric analysis and multi-objective optimization were further carried out to assess the performance of CLAS. Both O 2 purity and recovery increased monotonically with the cycle time, airflow rate, steam flow rate, and absorption pressure. Meanwhile, optimal O 2 productivity and power consumption can only be achieved by specific combinations of these parameters. The optimized results showed that CLAS can be highly competitive when compared to conventional pressure swing adsorption (PSA) or cryogenic distillation. The 5-step cycle configuration achieved a minimum power consumption of 118 kW·h for producing 1 ton O 2 with ≥ 95% purity. The maximum O 2 productivity reached 0.0932 g O2 /(g sorbent ·h) with 390 kW·h/ton O 2 of energy consumption (95% pure). The optimization results also indicate that CLAS can potentially be more efficient than cryogenic distillation even when the required O 2 purity is above 99%.

42 ENGINEERING↗

Margin to onset of nucleate boiling and flow instability studies for preliminary MITR design-demonstration element thermal-hydraulics

This study covers the onset of flow instability (OFI) preliminary results obtained from leveraging correlations, in addition to the preliminary thermal hydraulics results such as pressure, flow velocity, temperature, and oxide layer over the design demonstration experiment (DDE) for the Massachusetts Institute of Technology Reactor (MITR). Current computational fluid dynamics (CFD) models in fluid structure interaction (FSI) have added the capability of assessing margins to onset of nucleate boiling (ONB). This study initiates the capability to model the margin to OFI and ONB presented for the MITR. Such study is supportive of the United States High Performance Research Reactor (USHPRR) program. Previous studies provided preliminary thermal-hydraulic and mechanical analyses of the hydrodynamic effects in the MITR DDE under conservative approximations for plate power distribution. In conclusion, this study focuses on providing insights into the OFI future research direction optimizing the transport of thermal energy, mass-flow rates, flow-channel geometries, and boundary conditions.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A dimensionless number for high-throughput design of multi-principal element alloys in directed energy deposition

The present study demonstrates the development of a dimensionless number to predict the build height in the additive manufacturing technique of directed energy deposition (DED). The build height can also be used to estimate the dendrite arm spacing and, thus, the cooling rate in the fabrication of samples. A baseline sample, 316L stainless steel, was used to fit the build height to the dimensionless number. A range of process parameters, including laser power, laser feed rate, powder flow rate, layer thickness, and hatch spacing, were varied. Based upon dendrite arm spacing, the estimated cooling rate varied between 102 and 104 K/s. Using the fitted relationship for the stainless steel, high-throughput (HT) processing of multi-principal element alloys (MPEAs) was performed. For this study, HT is the ability to fabricate a batch of 25 bulk samples (∼1 cm3) with different compositions within a 5-h period with ±10 at. % accuracy. A range of compositions using in situ alloying of elemental powders in the Fe–Ni–Cr–Mo system were made. The MPEAs' build height followed the same relationship to the dimensionless number as the 316L alloy. The dimensionless number predicts both macro and meso-scale features in HT processing, thus offering a design tool for choosing process parameters in DED additive manufacturing. Also, the ability to control or increase cooling rates can enhance the ability to promote metastability as well as control meso-scale chemical distributions of alloy samples.

Islam, Zahabul (ORCID:0000000169613209)↗

Advanced Multi-Dimensional Capacitance Sensors Based Subsea Multiphase Mass Flow Meter to Measure and Monitor Offshore Enhanced Oil

This innovation is based on advanced multi-dimensional extensions of ECVT sensors that involve ECVT, Displacement Current Phase Tomography (DCPT), Maxwell-Wagner-Sillars polarization effect (MWS), and velocimetry which exploit the variation of electric properties between the oil, water, and gas coming out of the well. Capacitance sensors are embedded on the inside of a pipe spool and placed in-line with the subsea oil line piping on the extraction end of the well. The difference in dielectric, dielectric loss, surface polarization, and velocity of each phase are used in the multi-dimensional algorithm to measure the volume fraction, distribution, velocity, mass flow rate, and flow regime of the mixture moving through the pipe.

02 PETROLEUM↗

Steel Creek, Pen Branch, and D-Area Watershed Stream Gauging Stations

A network of stream gauging systems were installed in the Steel Creek, Pen Branch, D-Area Discharge Canal, and the D006 Stream in support of the groundwater modeling efforts for the P-Area Groundwater Operable Unit (OU); Chemical, Metals, and Pesticides (CMP) Pits OU; and the D-Area Watershed, respectively. Each location is monitored by a MACE Floseries3 FloPro data logger and a MACE doppler ultrasonic area/velocity sensor. Each stream gauging system is powered by an internal 12-volt battery supplied by a solar panel with a trickle charger. Information collected by each data logger is logged internally and telecommunicated via a cellular network to an online server for real time analysis and monitoring. The MACE doppler ultrasonic area/velocity sensor can measure stream depth and velocities to give output values of flow rates, total flow, net flow, and volumes. The water depth is measured by a ceramic pressure transducer located on the top of the sensor. The velocity is measured by a continuous wave doppler sensor to give an average velocity across the whole stream profile. This report discusses the equipment and methods used to install continuous stream gauging stations and provides a summary of data collected through FY2025.

54 ENVIRONMENTAL SCIENCES↗

Machine learning modeling and model predictive control of a closed-circuit reverse osmosis system

Closed-circuit reverse osmosis (CCRO) offers a flexible and energy-efficient alternative to conventional reverse osmosis by operating in a semi-batch mode that recycles brine, enabling higher recovery rates and reduced specific energy consumption (SEC). However, developing accurate, system-level dynamic models for CCRO remains challenging due to its nonlinear, multi-phase operation and sensitivity to variable feed water conditions. Traditional modeling approaches, such as NARMAX (nonlinear autoregressive moving average with exogenous inputs), often struggle to generalize across varying inlet feed concentrations, necessitating frequent parameter re-estimation and limiting their utility for real-time control applications. To address these limitations, we developed a long short-term memory (LSTM) neural network model trained on an extensive experimental data set from a CCRO pilot plant. The model accepts three inputs, feed flow rate, recirculation flow rate, and initial feed conductivity, and predicts three key outputs: reject conductivity, feed pump power draw, and recirculation pump power draw. We validated the LSTM model against experimental data, demonstrating its ability to distinguish between different feed conductivities and adapt to variable flow rates. Subsequently, we incorporated the LSTM model within a nonlinear model predictive control (MPC) scheme and conducted closed-loop simulations to optimize the integrated SEC (iSEC). In conclusion, the results project up to a 6% reduction in iSEC by using MPC to optimize performance over the entire experiment duration, without requiring any random excitation for data collection or parameter re-estimation.

Desalination↗

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↗