Next Article in Journal
An Environmental and Economic Analysis of Flocculation Technology Applied to a Corn-Based Ethanol Plant
Next Article in Special Issue
Thermally Optimum Spacing between Inner Plates in Natural Convection Flows in Cavities by Numerical Investigation
Previous Article in Journal
Simulation of Solid Oxide Fuel Cell Anode in Aspen HYSYS—A Study on the Effect of Reforming Activity on Distributed Performance Profiles, Carbon Formation, and Anode Oxidation Risk
Previous Article in Special Issue
Research on Rotordynamic Characteristics of Pump Annular Seals Based on a New Transient CFD Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cavitating Flow Suppression in the Draft Tube of a Cryogenic Turbine Expander through Runner Optimization

1
College of Mechanical and Transportation Engineering, China University of Petroleum (Beijing), Beijing 102249, China
2
Department of Energy and Power Engineering, State Key Laboratory of Hydro Science and Engineering, Tsinghua University, Beijing 100084, China
*
Author to whom correspondence should be addressed.
Processes 2020, 8(3), 270; https://doi.org/10.3390/pr8030270
Submission received: 19 January 2020 / Revised: 23 February 2020 / Accepted: 24 February 2020 / Published: 27 February 2020
(This article belongs to the Special Issue Advancement in Computational Fluid Mechanics and Optimization Methods)

Abstract

:
The application of a cryogenic liquefied natural gas expander can reduce the production of flash steam and improve the efficiency of natural gas liquefaction. Like traditional hydraulic machinery, cavitation will occur during the operation of a liquefied natural gas expander, in particular, there is a strong vortex flow in the draft tube, and the cavitation phenomenon is serious. In this paper, the energy loss coefficient of the draft tube is used to describe the cavitation flow in the draft tube, and the goal of reducing the cavitation in the draft tube is achieved through the optimization design of the runner. Different runner models within the range of design parameters were obtained using the Latin hypercube test, and the relationship between design parameters and objective functions is constructed by a second-order response surface model. Finally, the optimized runners were obtained using a genetic algorithm. The effects of blade loading distribution and blade lean angles on the cavitation in the draft tube were studied. According to the optimization results, the blade loading distribution and blade lean angles are recommended in the end.

1. Introduction

In the process of liquefaction, storage, and transportation of liquefied natural gas (LNG), it is necessary to reduce the pressure of high-pressure LNG. In the traditional process, a Joule–Thomson (J–T) valve is usually used to complete this. However, when reducing the pressure through the throttle valve, it is very easy to vaporize and flash, reducing the liquefaction efficiency and increasing the total energy consumption. In recent years, the cryogenic liquid expander has gradually replaced the J–T valve. The use of a cryogenic liquid expander can reduce the production of flash steam, and can recover additional high pressure to generate electricity and improve energy efficiency. The research in references [1,2] shows that the use of a cryogenic expander can increase the liquefaction rate by 5%.
In hydraulic machinery, cavitation occurs easily, which affects the operational stability of the whole unit. A cryogenic liquid expander is similar to conventional hydraulic machinery, cavitation may still occur [3]. The cavitation experiment of the cryogenic medium is difficult to operate, and the related research is less. Hord [4,5] carried out a series of detailed experiments with airfoils, blunt bodies and other experimental bodies of different sizes, and obtained the cavitation flow images of liquid nitrogen and liquid hydrogen, which has become an important index to verify the accuracy of numerical simulations.
Scholars have carried out a large number of numerical simulation studies on the cavitation flow simulation. Hsiao [6] proposed a multiscale two-phase flow model based on a coupled Eulerian/Lagrangian to capture the sheet cavitation formation and bubble cloud shedding on a hydrofoil. Du et al. [7] proposed a new cavitation model considering the evolution of bubble number density as an important factor. Örley [8] considered the compressibility of all phases in order to accurately capture the pressure wave dynamics of collapse events. Kähler [9] used the Boltzmann simulation method to calculate the cavitation flow of a liquid moving past a constraint. For the simulation of cryogenic fluid cavitation, Utturkar [10] used a modified cavitation model to study the steady-state cavitation characteristics of a cryogenic hydrofoil by calibrating the correlation model coefficients. Tailan et al. [11,12,13] introduced thermodynamic terms into the normal temperature cavitation model, which can better predict the temperature and pressure drop in the low-temperature cavitation zone. In the field of rotating machinery, research on cryogenic medium cavitation is mainly focused on the inducer in a cryogenic pump [14,15,16].
The flow in the draft tube of rotating machinery is also the focus of research. Liu [17] used the dynamic grid method to catch the pressure fluctuations in a pump–turbine draft tube. In the study by Zhang [18], the vortex identification methods were reviewed to reveal the complex vortex structures in hydroturbines. Arispe [19] used the computational fluid dynamics (CFD) method to obtain a draft tube geometry that improved the hydrodynamic performance. The pressure fluctuations generated by vortex ropes in the draft tube were studied when hydraulic turbines operated at off-design conditions [20].
Multiobjective evolutionary algorithms are widely used in runner optimization design systems [21,22,23]. The optimization strategy usually consists of a runner design method, design of experiment (DOE), CFD analysis, response surface methodology, and multiobjective genetic algorithm method. This strategy was successfully used in the design of a pump–turbine runner [24,25], after the optimization, the performance of the runner was obviously improved. Because of its simplicity, this optimization strategy could be used in the development of fluid machines.
Because the outflow from the runner contains a strong swirl component of velocity, the cavitation in the draft tube is serious. The flow in the draft tube is closely related to the flow condition of the runner, hence, optimizing the shape of the runner blades can achieve the purpose of minimizing the cavitation in the draft tube. This paper carries out the optimal design of the runner to restrain the cavitation in the draft tube and reduce the production of gas using the energy loss coefficient to describe the character of the draft tube.

2. Physical Model

The cryogenic expander studied in this paper is shown in Figure 1 and Figure 2, including the generator, shaft, runner, and guide vanes. The designed flow rate of the expander is 40 m3/h, the designed head of the runner is 98 m, the designed speed is 3000 rpm and the hydraulic efficiency is 74%. The parameters of the LNG cryogenic expander are listed in Table 1.

3. Numerical Method

3.1. Governing Equations

For cavitation flow, the homogeneous flow model is generally used. The homogeneous flow model defines the mixture of variables by averaging the parameters, and then the mixture conservation equation is used for the calculation. The governing equation is as follows:
( ρ m u j ) x j = 0
The momentum equation is:
( ρ m u i u j ) x i = p x i + x j [ μ e f f ( u i x j + u j x i 2 3 u k x k δ i j ) ]
The energy equation is:
( ρ m C p , m u j T ) x j = x j [ ( λ + λ t ) T x j ] m ˙ L
The fluid properties, ρ m , μ m , λ m , C p , m , are calculated from the local value of volume fractions of vapor and fluid phases:
ρ m = α v ρ v + ( 1 α v ) ρ l
μ m = α v μ v + ( 1 α v ) μ l
λ m = α v λ v + ( 1 α v ) λ l
C p , m = α v C p , v + ( 1 α v ) C p , l
where ρ is the density, u and T are the velocity and temperature, respectively, which consist of the two phases. p is the pressure and α is the volume fraction. μ and λ are the dynamic viscosity and thermal conductivity of the homogeneous mixture, respectively. C p is the specific heat and L is the latent heat due to the phase change. Subscripts v and l refer to vapor and liquid, respectively. m ˙ = m ˙ + m ˙ represents the mass transfer source term. m ˙ + and m ˙ refer to the mass transfer rate source term due to evaporation and condensation, respectively. The vapor volume fraction is solved using the following transport equation:
t ( ρ v α v ) + x i ( ρ v α v u i ) = m ˙ + + m ˙
Equations (1)–(8) are discretized by ANSYS-CFX software (SAS IP, Inc., Pittsburgh, PA, USA) to complete the simulation calculation of the expander.

3.2. Cavitation Model

The two-phase-flow mixture model and the Zwart–Gerber–Belamri mass transfer cavitation model [26] with a modified saturated pressure were executed after cavitation. The liquid–vapor evaporation and condensation rates for the present transport-based cavitation model are, respectively, shown as:
m ˙ + = F v a p 3 r n u c ( 1 α v ) ρ v R B 2 3 p v ( T ) p ρ l ,   if   p p v
m ˙ = F cond 3 α v ρ v R B 2 3 p p v ( T ) ρ l ,   if   p p v
where F v a p and F c o n d are the empirical coefficients of evaporation and condensation, respectively. In addition, r n u c is the nucleation site volume fraction and R B indicates the nucleation site radius. p v ( T ) stands for the liquid vaporization pressure as a function of the temperature. According to the researches in the simulation of cryogenic cavitation [27,28,29], F v a p and F c o n d have a significant impact on predictions of simulation, and the recommended values for those coefficients are F vap = 5 , F c o n d = 0.001 , r n u c = 5 × 10 4 and R B = 2 × 10 6 m . The effects of turbulence on the cavitation are considered by correcting the saturation pressure as:
p v = p s a t + p t / 2
p t = 0.39 ρ κ
where p s a t is the local saturation pressure dominated by the temperature, p t is the turbulent pressure and κ is local turbulence energy.
A series of cryogenic liquid cavitation experiments conducted by Hord in the 1970s has become the experimental basis to verify the accuracy of the cavitation model. Because of the difficulty in performing the cryogenic cavitation experiment, Hord’ results are usually used to verify the accuracy of the simulation. In this part of the study, the 290 C experimental results are used as the verification target, the grid model is shown in Figure 3. The temperature and vapor volume fraction distribution of simulation results are shown in Figure 4. Figure 5 compares the pressure and temperature between the experiment and simulation, and the results are in good agreement with the experimental data. In this way, the accuracy of the cavitation model in calculating cryogenic cavitation is verified. Then, the cavitation model is used in the expander numerical simulation.

3.3. Numerical Simulation

The whole calculation basin of the LNG cryogenic expander was meshed by the ANSYS ICEM software (SAS IP, Inc., Pittsburgh, PA, USA). The meshing is shown in Figure 6, and Figure 7 presents the grid independence. It is shown that the head changes a little when the total grid number is higher than 2.57 × 106. The final calculation grid model is guide vane 1.37 × 106, runner 1.20 × 106, draft tube 0.50 × 106.
The commercial software ANSYS CFX was used to simulate the internal flow characteristics of the expander flow field. The turbulence model is selected as the SST k–ω model. The second-order upwind scheme was employed for the convective terms, and the central difference scheme was used for the diffusive terms in the governing equations in the numerical simulation. The frozen rotor model is used for the rotor–stator interfaces (guide vane–runner and runner–draft tube interfaces). The second wind scheme was used for the convective energy terms. The convergence criterion was set to 1 × 10 6 . In this study, the velocity inlet and the pressure outlet conditions were set for the calculation. In addition, a smooth no-slip wall condition was imposed for the solid surfaces. In the cavitation simulation of liquid hydrogen, the thermodynamic properties of both liquid and vapor are updated from NIST database [30]. The head H and efficiency η are defined in Equations (13) and (14), respectively
H = P I n f l o w P O u t f l o w ρ g
η = M ω ρ g Q H
where Q is discharge, and ω is the angular velocity, which can be calculated using the rotational speed. The torque acting on the runner M was calculated through ANSYS CFX simulation.
Because LNG is very dangerous, in our experiment, using liquid nitrogen as the medium, the flow rate of the expander is constant, and the expander performance under different operating conditions is obtained by changing the speed of the expander. The enclosed layout of the expander model test rig is presented in Figure 8. The density of the liquid nitrogen is 808 kg/m3 and the temperature is 77 K. The numerical simulation results of the LNG expander are compared with the experimental results in Figure 9. From the comparative analysis, we can see that in the range of working conditions, the numerical simulation results are in good agreement with the test results, the error between the simulation result and the test result at the design point is 1.5%, and the errors at low-speed and high-speed conditions are both less than 10%. Therefore, the numerical simulation results are reliable.

4. Optimization Design System

4.1. Runner Design Method

The 3D design software TURBOdesign 5.2 (Advanced Design Technology, London, UK) [31,32] was used to describe the LNG cryogenic expander runner parametrically. For incompressible potential flow, the pressure distribution can be expressed as follows:
p + p = 2 π B ρ W b l ( r V ¯ θ ) m ,  
where p + p is the pressure difference across the blade, W b l is the relative meridional velocity on the blade surface, r V ¯ θ is swirl velocity and B is the number of blades, ( r V ¯ θ ) / m is also called the blade loading distribution.
When using TURBOdesign 5.2 to design runner blades, blade loadings were the most important parameters in determining the blade shape [31,32,33]. Blade loading distributions were given along the hub and shroud streamlines. The runner blade loadings of the hub and shroud were determined by using linear interpolation. The blade loading distribution is shown in Figure 10. Along each streamline, a three-segment distribution was adopted, including the loading at the leading edge of the runner blade DRVT, the linear slope SLOPE and connection point locations NC and ND. In the optimization process, extensive trial designs were made to check whether runners with reasonable shapes could be obtained. After these trial designs, three variables were fixed at N C h = 0.4 , N D h = 0.8 , and N D s = 0.8 , where subscripts h and s refer to hub and shroud, respectively, and the variation ranges of the other runner blade variables were determined as shown in Table 2.

4.2. Optimization Strategy

The response surface model (RSM) model is mainly used to construct the approximate relationship between optimization variables and optimization objectives to provide a basis for further optimization. The design of experiment (DOE) is an important mathematical method for systematically optimizing the sampling generation. In this study, the runner design parameters are taken as optimization variables, and the Latin hypercube test design method [34] is adopted to obtain different runner shapes in the space of the design parameters. The second-order response surface model is used to establish the approximate correspondence between independent variables and corresponding dependent variables, the expression is as follows:
y = β 0 + j = 1 n β j x j + j = 1 n β j j x j 2 + i j β i j x i x j ,  
where y is the optimization objective, n is the number of runner design parameters, and x j , β j are the runner parameters and constants, respectively.
The improved non-dominated sorting genetic algorithm (NSGA-II) was used to optimize the runner parameters in the design space. This algorithm uses a fast non-dominated sorting method to reduce the complexity of calculation, expands the optimization space by introducing elite strategy, ensures that the improved species will not be mistakenly abandoned in the process of evolution, and improves the optimization accuracy. The parameter settings for NSGA-II are shown in Table 3.
For a description of the flow loss of the draft tube, the energy loss coefficient (ζ) and the average pressure recovery coefficient ( C p m ) in the draft tube are commonly used to characterize the flow in the draft tube, the expressions are as follows:
ζ = 1 A i n i n P t d A 1 A o u t o u t P t d A 0.5 ρ ( Q A i n ) 2 ,
C p m = 1 A o u t o u t P d A 1 A i n i n P d A 0.5 ρ ( Q A i n ) 2 ,  
where A i n ,   A o u t are the inlet and outlet areas of the draft tube, respectively, and P t is total pressure, P is static pressure, P t = P + 0.5 ρ ( u 2 + v 2 + w 2 ) .
In Chen’s [35] multiobjects optimization strategy, the loss coefficient is chosen as an objective to optimize the structure of a Francis turbine draft tube. In the study by Galva’n [36], the loss coefficient of the draft tube is more sensitive to the pressure distribution of the draft tube inlet, which can better describe the velocity distribution of the draft tube inlet and the flow loss caused by the change of the draft tube inlet velocity. Therefore, the loss coefficient ζ is chosen as one of the objective functions of multiobjective optimization. At the same time, to ensure the performance of the runner, the head of the single runner in the optimization process is also set as the objective function.
In the process of optimization, the single channel of the runner is used for the simulation calculation, the boundary conditions are set as shown in Figure 11, with liquefied nitrogen used in the experiment as the medium. Compared with the full-channel model, the grid of the single channel is greatly reduced, which can save calculation time.
The runner test model is obtained by changing the runner design parameters, and the head and energy loss coefficient in the draft tube under different runner design parameters are obtained by CFD numerical simulation. The explicit relationship between the runner design parameters and the objective function is obtained using the second-order RSM, and then the optimal solutions are obtained using the NSGA-II optimization method. According to the optimization results, the Pareto curve of the optimization target is obtained, and then the optimal design parameters of the runner are chosen. iSIGHT platform is used to complete the modeling of RSM and the subsequent optimization process.

5. Optimization Result

Figure 12 presents the black plots of the multiobjective optimization results, and the blue plots are the results on the Pareto front. Models A, B, and C (the red plots) were selected on the Pareto curve for further analysis. A comparison between the predicted values of the genetic algorithm and the numerical simulation results of CFD is shown in Table 4. There are some errors in the head and energy loss coefficients of the three models, but the results of different runners have the same trend, which verifies that the genetic algorithm can provide a reference for the design of runners.

5.1. Effect of Blade Loading

Table 5 shows the runner blade loading of the optimized models chosen on the Pareto curve. As shown in Figure 13, the blade loading distributions, which show great similarity, are aft-loading and fore-loading on the hub and shroud for runners A, B, and C. Figure 14 shows the runners’ shapes.
The pressure and gas volume fraction distribution in the draft tube of the three runner models and the original model are shown in Figure 15 and Figure 16. When the low-pressure area in the draft tube is large, the cavitation is serious, and the gas distribution area is large. Through the calculation of the numerical simulation results, when the cavitation is serious, the energy loss coefficient in the draft tube is larger, when the energy loss coefficient is small, the distribution of the low-pressure area in the draft tube is smaller, and the cavitation can be better suppressed. It is proven that the energy loss coefficient can reflect the severity of cavitation in the draft tube.
According to the calculation results, the cavitation of the original model is the most serious, the energy loss coefficient of the draft tube is the largest, the inhibition effect of model A is the best, and the energy loss coefficient is the smallest.
According to the optimization results and CFD analysis, it can be found that runners with aft-loading and fore-loading on the hub and shroud have better performance considering the effects on the runner head and draft tube cavitation.

5.2. Effect of The Blade Lean

Based on the optimization results, more runners with the same blade loading distribution and different blade lean were designed to be investigated. Table 6 shows the blade loading distribution, with the blade lean ranging from −10° to 10°, as shown in Figure 17. Figure 18 shows the gas volume fraction in the draft tube, and the energy loss coefficients are also listed in Table 6.
The results show that the runners with large positive or negative blade lean angles have large energy loss coefficients and the gas volume fraction distribution in the draft tube is worsening. When the blade lean angles are between 0° and 5°, the cavitation in the draft tube is improved, which can also be inferred from the optimization results on the Pareto curve. Therefore, large blade lean angles are not recommended to be used for the expander runner considering the cavitation in the draft tube.

6. Conclusions

A strategy for cavitation suppression in an LNG cryogenic expander draft tube by a multiobjective optimization method was employed in the present study. The process combines the runner design method, DOE, RSM, CFD analysis and genetic algorithm.
During the optimization, runners with different design parameters were obtained using the Latin hypercube test. The second-order RSM was used to describe the relationship between runner design parameters and objective functions, and then the final optimization models were obtained using the NSGA-II genetic algorithm.
The accuracy of the simulation of cryogenic liquid cavitation is verified by comparing with the experimental results of Hord, and the cavitation model is applied to the simulation of the cryogenic expander to obtain the cavitation flow in the draft tube.
The energy loss coefficient of the draft tube was used to evaluate the flow in the draft tube. Through the numerical simulation calculation, it was proven that the energy loss coefficient of the draft tube can be used to evaluate the cavitation severity of the draft tube. The more serious the cavitation is, the greater the energy loss coefficient is, as shown in Figure 16.
During the optimization process, the single runner head and the energy loss coefficient of the draft tube were taken as the optimization objectives. After optimization, the head is obviously increased, and the energy loss coefficient is reduced. The low-pressure area in the draft tube is obviously reduced and the cavitation is restrained to a certain extent. The runners on the Pareto curve have similar blade loadings. Considering the expander head and cavitation in the draft tube, it is recommended to design the expander runner with aft-loading and fore-loading on the hub and shroud.
The effects of blade lean angles for the cavitation in the draft tube were studied. It was found that large positive and negative angles are not recommended, and the runners with blade lean angles ranging from 0° to 5° have better performance for cavitation in the draft tube.

Author Contributions

Methodology, N.H. and B.Z.; Software, B.Z.; resources, Z.L.; writing—original draft preparation, N.H.; writing—review and editing, B.Z.; All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (grant numbers: 51736008, 51679122); the National Key Research and Development Program of China (grant number: SQ2018YFB060154-02); the High Technology ship Research of Ministry of Industry and Information Technology of China (Grant No. [2014]503).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kimmel, H.E. Cryogenic LNG Expanders Reduce Natural Gas Liquefaction Costs; Natural Gas Logistics, Handling and Contracts: Bali, Indonesia, 2011; pp. 18–22. [Google Scholar]
  2. Chen, J.; Hua, Y.; Su, Q.; Wan, X.; Li, Z. Development and industrial tests of the first LNG hydraulic turbine system in China. J. Nat. Gas Ind. B 2016, 3, 283–290. [Google Scholar] [CrossRef]
  3. Song, P.; Sun, J.; Wang, K. Swirling and cavitating flow suppression in a cryogenic liquid turbine expander through geometric optimization. J. Power Energy 2015, 229, 628–646. [Google Scholar] [CrossRef]
  4. Hord, J. Cavitation in Liquid Cryogens II-Hydrofoils; National Aeronautics and Space Administration: Washington, DC, USA, 1973; CR-2156. [Google Scholar]
  5. Hord, J. Cavitation in Liquid Cryogens III-Ogives; National Aeronautics and Space Administration: Washington, DC, USA, 1973; CR-2242. [Google Scholar]
  6. Hsiao, C.T.; Ma, J.; Chahine, G.L. Multiscale tow-phase flow modeling of sheet and cloud cavitation. Int. J. Multiph. Flow 2017, 90, 102–117. [Google Scholar] [CrossRef] [Green Version]
  7. Du, T.Z.; Wang, Y.W.; Liao, L.J.; Huang, C.G. A numerical model for the evolution of internal structure of cavitation cloud. J. Phys. Fluids 2016, 28, 3–20. [Google Scholar] [CrossRef]
  8. Örley, F.; Trummler, T.; Hickel, S.; Mihatsch, M.S.; Schmidt, S.J.; Adams, N.A. Large-eddy simulation of cavitating nozzle flow and primary jet break-up. J. Phys. Fluids 2015, 27, 086101. [Google Scholar] [CrossRef] [Green Version]
  9. Kähler, G.; Bonelli, F.; Gonnella, G.; Lamura, A. Cavitation inception of a van der Waals fluid at a sack-wall obstacle. J. Phys. Fluids 2015, 27, 123307. [Google Scholar] [CrossRef] [Green Version]
  10. Utturkar, Y.; Wu, J.; Wang, G.; Shyy, W. Recent progress in modeling of cryogenic cavitation for liquid rocket propulsion. J. Prog. Aerosp. Sci. 2005, 41, 558–608. [Google Scholar] [CrossRef]
  11. Chen, T.R.; Huang, B.; Wang, G.Y. Numerical study of cavitating flows in a wide range of water temperatures with special emphasis on two typical cavitation dynamics. Int. J. Heat Mass Transr. 2016, 101, 886–900. [Google Scholar] [CrossRef]
  12. Chen, T.; Huang, B.; Wang, G.; Zhang, H.; Wang, Y. Numerical investigation of thermo-sensitive cavitating flows in a wide range of free-stream temperatures and velocities in fluoroketone. Int. J. Heat Mass Transf. 2017, 112, 125–136. [Google Scholar] [CrossRef]
  13. Chen, T.; Chen, H.; Liang, W.; Huang, B.; Xiang, L. Experimental investigation of liquid nitrogen cavitating flows in converging-diverging nozzle with special emphasis on thermal transition. Int. J. Heat Mass Transf. 2019, 132, 618–630. [Google Scholar] [CrossRef]
  14. Cao, X.L.; Zhang, X.B.; Qiu, L.M.; Gan, Z.H. Validation of full cavitation model in cryogenic fluids. J. Chin. Sci. Bull. 2009, 54, 1633–1640. [Google Scholar] [CrossRef] [Green Version]
  15. Huang, B.; Wu, Q.; Wang, G. Numerical investigation of cavitating flow in liquid hydrogen. Int. J. Hydrogen Energy 2014, 39, 1698–1709. [Google Scholar] [CrossRef]
  16. Fei, T.; Wen, L.J. Influence of thermodynamics effect on inducer rotating cavitation under low temperature condition. J. Rocket Propuls. 2013, 39, 29–34. [Google Scholar]
  17. Liu, Q.Z.; Su, W.T.; Li, X.B.; Zhang, Y.N. Dynamic characteristics of load rejection process in a reversible pump-turbine. J. Renew. Energy 2020, 146, 1922–1931. [Google Scholar] [CrossRef]
  18. Zhang, Y.; Liu, K.; Xian, H.; Du, X. A review of methods for vortex identification in hydroturbines. J. Renew. Sustain. Energy Rev. 2018, 81, 1269–1285. [Google Scholar] [CrossRef]
  19. Arispe, T.M.; De Oliveira, W.; Ramirez, R.G. Francis turbine draft tube parameterization and analysis of performance characteristics using CFD techniques. J. Renew. Energy 2018, 127, 114–124. [Google Scholar] [CrossRef]
  20. Yu, A.; Zou, Z.; Zhou, D.; Zheng, Y.; Luo, X. Investigation of the correlation mechanism between cavitation rope behavior and pressure fluctuations in a hydraulic turbine. J. Renew. Energy 2020, 147, 1199–1208. [Google Scholar] [CrossRef]
  21. Zhu, B.; Wang, X.; Tan, L.; Zhou, D.; Zhao, Y.; Cao, S. Optimization design of a reversible pump–turbine runner with high efficiency and stability. J. Renew. Energy 2015, 81, 366–376. [Google Scholar] [CrossRef]
  22. Chirkov, D.V.; Ankudinova, A.S.; Kryukov, A.E.; Cherny, S.G.; Skorospelov, V.A. Multi-objective shape optimization of a hydraulic turbine runner using efficiency, strength and weight criteria. J. Struct. Multidiscip. Optim. 2018, 58, 627–640. [Google Scholar] [CrossRef]
  23. Lu, Y.M.; Wang, X.F.; Wang, W.; Zhou, F.M. Application of the Modified Inverse Design Method in the Optimization of the Runner Blade of a Mixed-Flow Pump. Chin. J. Mech. Eng. 2018, 31, 105. [Google Scholar] [CrossRef] [Green Version]
  24. Liu, L.; Zhu, B.; Bai, L.; Liu, X.; Zhao, Y. Parametric Design of an Ultrahigh-Head Pump-Turbine Runner Based on Multiobjective Optimization. Energies 2017, 10, 1169. [Google Scholar] [CrossRef] [Green Version]
  25. Ma, Z.; Zhu, B.; Rao, C.; Shangguan, Y. Comprehensive Hydraulic Improvement and Parametric Analysis of a Francis Turbine Runner. Energies 2019, 12, 307. [Google Scholar] [CrossRef] [Green Version]
  26. Zwart, P.J.; Gerber, A.G.; Belamri, T. A two-phase flow model for predicting cavitation dynamics. In Proceedings of the Fifth International Conference on Multiphase Flow, Yokohama, Japan, 30 May–4 June 2004. [Google Scholar]
  27. Goel, T.; Thakur, S.; Haftka, R.T.; Shyy, W.; Zhao, J. Surrogate model-based strategy for cryogenic cavitation model validation and sensitivity evaluation. Int. J. Numer. Methods Eng. 2008, 58, 969–1007. [Google Scholar] [CrossRef] [Green Version]
  28. Senocak, I.; Shyy, W. Interfacial dynamics-based modeling of turbulent cavitating flows Part-1: Model development and steady-state computations. Int. J. Numer. Methods Fluids 2004, 44, 975–995. [Google Scholar] [CrossRef]
  29. Sun, T.Z.; Ma, X.F.; Wei, Y.J.; Wang, C. Computational modeling of cavitating flows in liquid nitrogen by an extended transport-based cavitation model. Sci. China Technol. Sci. 2016, 59, 337–346. [Google Scholar] [CrossRef]
  30. Lemmon, E.W.; McLinden, M.O.; Huber, M.L. Reference Fluid Thermodynamic and Transport Properties, Standard Database 23 version 7.0; NIST: Gaithersburg, MD, USA, 2002. [Google Scholar]
  31. Zangeneh, M. A compressible three-dimensional design method for radial and mixed flow turbomachinery blades. Int. J. Numer. Methods Fluids 1991, 13, 599–624. [Google Scholar] [CrossRef]
  32. Zangeneh, M.; Goto, A.; Harada, H. On the design criteria for suppression of secondary flows in centrifugal and mixed flow impellers. J. Turbomach. 1998, 120, 723–735. [Google Scholar] [CrossRef]
  33. Bonaiuti, D.; Zangeneh, M. On the coupling of inverse design and optimization techniques for the multiobjective, multipoint design of turbomachinery blades. J. Turbomach. 2009, 131, 21014–21029. [Google Scholar] [CrossRef]
  34. Myers, R.H.; Montgomery, D.C.; Anderson-Cook, C.M. Response Surface Methodology: Process and Product Optimization Using Designed Experiments; John Wiley & Sons: Hoboken, NJ, USA, 2009. [Google Scholar]
  35. Chen, Z.; Baek, S.H.; Cho, H.; Choi, Y.D. Optimal design of J-groove shape on the suppression of unsteady flow in the Francis turbine draft tube. J. Mech. Sci. Technol. 2019, 33, 2211–2218. [Google Scholar] [CrossRef]
  36. Galván, S.; Rubio, C.; Pacheco, J.; Mendoza, C.; Toledo, M. Optimization methodology assessment for the inlet velocity profile of a hydraulic turbine draft tube: Part I—Computer optimization techniques. J. Global Optim. 2013, 55, 53–72. [Google Scholar] [CrossRef]
Figure 1. Schematic of a liquefied natural gas (LNG) cryogenic expander.
Figure 1. Schematic of a liquefied natural gas (LNG) cryogenic expander.
Processes 08 00270 g001
Figure 2. Computations domain of an LNG cryogenic expander.
Figure 2. Computations domain of an LNG cryogenic expander.
Processes 08 00270 g002
Figure 3. Hydrofoil computational domain grid.
Figure 3. Hydrofoil computational domain grid.
Processes 08 00270 g003
Figure 4. The computed temperature and vapor volume fraction contours.
Figure 4. The computed temperature and vapor volume fraction contours.
Processes 08 00270 g004
Figure 5. Comparison of coefficients between the simulation and experiment.
Figure 5. Comparison of coefficients between the simulation and experiment.
Processes 08 00270 g005
Figure 6. Grid profile of the flow passage of the LNG cryogenic expander.
Figure 6. Grid profile of the flow passage of the LNG cryogenic expander.
Processes 08 00270 g006aProcesses 08 00270 g006b
Figure 7. The grid independence.
Figure 7. The grid independence.
Processes 08 00270 g007
Figure 8. Layout of the expander model test rig.
Figure 8. Layout of the expander model test rig.
Processes 08 00270 g008
Figure 9. Performance comparison between the numerical and experimental results.
Figure 9. Performance comparison between the numerical and experimental results.
Processes 08 00270 g009
Figure 10. Blade loading distributions.
Figure 10. Blade loading distributions.
Processes 08 00270 g010
Figure 11. Single-channel flow passage of runner.
Figure 11. Single-channel flow passage of runner.
Processes 08 00270 g011
Figure 12. Optimization results.
Figure 12. Optimization results.
Processes 08 00270 g012
Figure 13. Blade loading distributions of the optimized runners.
Figure 13. Blade loading distributions of the optimized runners.
Processes 08 00270 g013
Figure 14. Optimized Model A, B, and C blade shapes.
Figure 14. Optimized Model A, B, and C blade shapes.
Processes 08 00270 g014
Figure 15. Pressure distribution in the draft tube of the optimized and original runners.
Figure 15. Pressure distribution in the draft tube of the optimized and original runners.
Processes 08 00270 g015
Figure 16. Gas volume fraction distribution in the draft tube of the optimized and original runners.
Figure 16. Gas volume fraction distribution in the draft tube of the optimized and original runners.
Processes 08 00270 g016
Figure 17. Runners with different blade lean angles.
Figure 17. Runners with different blade lean angles.
Processes 08 00270 g017
Figure 18. Gas volume fraction distribution in the draft tube.
Figure 18. Gas volume fraction distribution in the draft tube.
Processes 08 00270 g018
Table 1. Parameters of the LNG cryogenic expander.
Table 1. Parameters of the LNG cryogenic expander.
ParametersValue
Runner inlet diameter D 2 / m m 216
Runner outlet diameter D 1 / m m 96
Number of runner blades Z b 7
Number of guide vanes Z s 9
Table 2. Variation range of runner blade parameters.
Table 2. Variation range of runner blade parameters.
Optimized InputsParametersValue
Blade loading D R V T h −0.2 to 0.2
D R V T s −2.0 to 2.0
S L O P E h −2.0 to 2.0
S L O P E s −0.2 to 0.2
N C s 0.2 to 0.4
Blade lean angle θ −10.0° to 10.0°
Table 3. Parameter settings for NSGA-II.
Table 3. Parameter settings for NSGA-II.
ParametersValue
Population size100
Number of generations100
Crossover probability0.9
Crossover distribution index10
Mutation distribution index20
Initialization modeRandom
Table 4. Comparison of the selected runners and the initial runner.
Table 4. Comparison of the selected runners and the initial runner.
Head (m)Loss Coefficient
RSMCFDRSMCFD
Original19.8540.00
Model A21.6521.889.009.55
Model B23.1023.6015.0014.50
Model C23.7024.1120.0020.45
Table 5. Optimized parameters of the chosen models.
Table 5. Optimized parameters of the chosen models.
D R V T h D R V T s S L O P E h S L O P E s N C s θ
Model A−0.0680−0.2000.884−2.000.3000.00
Model B−0.0840−0.1901.52−2.000.3262.20
Model C0.00−0.2001.70−1.900.2803.76
Table 6. Blade loading distribution of the designed models.
Table 6. Blade loading distribution of the designed models.
D R V T h D R V T s S L O P h S L O P s N C s θ Loss Coefficient
Model B10.00−0.2001.70−2.000.300−10.0012.98
Model B20.00−0.2001.70−2.000.300−5.0011.44
Model B30.00−0.2001.70−2.000.3000.009.14
Model B40.00−0.2001.70−2.000.3005.008.66
Model B50.00−0.2001.70−2.000.30010.0021.4

Share and Cite

MDPI and ACS Style

Huang, N.; Li, Z.; Zhu, B. Cavitating Flow Suppression in the Draft Tube of a Cryogenic Turbine Expander through Runner Optimization. Processes 2020, 8, 270. https://doi.org/10.3390/pr8030270

AMA Style

Huang N, Li Z, Zhu B. Cavitating Flow Suppression in the Draft Tube of a Cryogenic Turbine Expander through Runner Optimization. Processes. 2020; 8(3):270. https://doi.org/10.3390/pr8030270

Chicago/Turabian Style

Huang, Ning, Zhenlin Li, and Baoshan Zhu. 2020. "Cavitating Flow Suppression in the Draft Tube of a Cryogenic Turbine Expander through Runner Optimization" Processes 8, no. 3: 270. https://doi.org/10.3390/pr8030270

APA Style

Huang, N., Li, Z., & Zhu, B. (2020). Cavitating Flow Suppression in the Draft Tube of a Cryogenic Turbine Expander through Runner Optimization. Processes, 8(3), 270. https://doi.org/10.3390/pr8030270

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop