Evaluation of Different Turbulence Models Applied in Turbopump’s Hydraulic Turbine

2021 ◽  
Author(s):  
Daniel Ferreira Corrêa Barbosa ◽  
Daniel da Silva Tonon ◽  
Luiz Henrique Lindquist Whitacker ◽  
Jesuino Takachi Tomita ◽  
Cleverson Bringhenti

Abstract The aim of this work is an evaluation of different turbulence models applied in Computational Fluid Dynamics (CFD) techniques in the turbomachinery area, in this case, in an axial turbine stage used in turbopump (TP) application. The tip clearance region was considered in this study because it has a high influence in turbomachinery performance. In this region, due to its geometry and the relative movement between the rotor row and casing, there are losses associated with vortices and secondary flow making the flowfield even more turbulent and complex. Moreover, the flow that leaks in the tip region does not participate in the energy transfer between the fluid and rotor blades, degradating the machine efficiency and performance. In this work, the usual flat tip rotor blade geometry was considered. The modeling of turbulent flow based on Reynolds Averaged Navier-Stokes (RANS) equations predicts the variation of turbine operational characteristics that is sufficient for the present turbomachine and flow analysis. Therefore, the appropriate choice of the turbulence model for the study of a given flow is essential to obtain adequate results using numerical approximations. This comparison become important due to the fact that there is no general turbulence model for all engineering applications that has fluid and flow. The turbomachine considered in the present work, is the first stage of the hydraulic axial turbine used in the Low Pressure Oxidizer Turbopump (LPOTP) of the Space Shuttle Main Engine (SSME), considering the 3.0% tip clearance configuration relative to rotor blade height. The turbulence models evaluated in this work were the SST (Shear Stress Transport), the k-ε Standard and the k-ε RNG. The computational domain was discretized in several control volumes based on unstructured mesh. All the simulations were performed using the commercial software developed by ANSYS, CFX v15.0 (ANSYS). All numerical settings and how the boundary conditions were imposed at different surfaces are explained in the work. The boundary conditions settings follow the same rule used in the test facility and needs some attention during the simulations to vary the Blade-Jet-Speed ratio parameter adequately. The results from numerical simulations, were synthesized and compared with the experimental data published by National Aeronautics and Space Administration (NASA), in which the turbine efficiency and its jet velocity parameter are analyzed for each turbulence model result. The work fluid considered in this work was water, the same fluid used in the NASA test facility.

1995 ◽  
Author(s):  
S. Fan ◽  
B. Lakshminarayana

The objective of this paper is to investigate the three dimensional unsteady flow interactions in a turbomachine stage. A three-dimensional time accurate Euler code has been developed using an explicit four-stage Runge-Kutta scheme. Three-dimensional unsteady non-reflecting boundary conditions are formulated at the inlet and at the outlet of the computational domain to remove the spurious numerical reflections. The three-dimensional code is first validated for 2-D and 3-D cascades with harmonic vortical inlet distortions. The effectiveness of non reflecting boundary conditions is demonstrated. The unsteady Euler solver is then used to simulate the propagation of nozzle wake and secondary flow through rotor and the resulting unsteady pressure field in an axial turbine stage. The three dimensional and time dependent propagation of nozzle wakes in the rotor blade row and the effects of nozzle secondary flow on the rotor unsteady surface pressure and passage flow field are studied. It was found that the unsteady flow field in the rotor is highly three-dimensional and the nozzle secondary flow has significant contribution to the unsteady pressure on the blade surfaces. Even though the steady flow at the midspan is nearly two-dimensional, the unsteady flow is 3-D and the unsteady pressure distribution can not by predicted by a 2-D analysis.


Author(s):  
A. Namet-Allah ◽  
A. M. Birk

The current paper presents a cold flow simulation study of a low Mach number air-air ejector with a four ring entraining diffuser that is used in a variety of applications including infrared (IR) suppression of exhaust from helicopters and fixed wing aircraft. The main objectives of this investigation were to identify key issues that must be addressed in successful CFD modelling of such devices, and recognize opportunities to improve the performance of these devices. Two-dimensional CFD simulations were carried out using commercial software, Ansys14. Studies of mesh and domain size sensitivity were made to ensure the CFD results were independent of both factors. A turbulence model independence study using k-ε, k-ω and RSM turbulence models was performed to figure out the appropriate turbulence model that produced the best agreement with the experimental data for several of ejector performance criteria. The measured flow properties in the annulus were used as input boundary conditions for the CFD simulations. However, in the comprehensive turbulence model study, the measured flow parameters at the nozzle exit were also applied as inlet boundary conditions for the CFD simulations. The measured flow velocity at the nozzle exit, at one centerline section inside the mixing tube and at the diffuser exit and the system pressure recovery were compared with the CFD predictions. The ejector pumping ratios, back pressure coefficient and diffuser gap velocities were also compared. It was found that the RANS-based CFD predictions were sensitive to the changes in the ejector domain size, mesh refinement and inlet boundary condition locations. With the annulus inlet boundary conditions, the tested turbulence models under predicted the size of the core separation downstream of the system, back pressure, pumping ratio and pressure recovery in the mixing tube and diffuser. However, the ability of the RNG turbulence model to predict the ejector performance parameters was better than that of the other turbulence models at all inlet flow conditions. Nevertheless, applying the inlet boundary conditions at the nozzle exit enhanced the capability of the RANS-based turbulence model particularly in predicting the ejector pumping ratios, pressure recovery and the size of the core separation. Finally, the acceptable agreement between the experimental data and the CFD predictions provides a valid tool to continue improving these devices using CFD techniques.


Author(s):  
N S Tachos ◽  
A E Filios ◽  
D P Margaris

The analysis of the near and far flow fields of an experimental National Renewable Energy Laboratory (NREL) rotor, which has been used as the reference rotor for the Viscous and Aeroelastic Effects on Wind Turbine Blades (VISCEL) research program of the European Union, is described. The horizontal axis wind turbine (HAWT) flow is obtained by solving the steady-state Reynolds-averaged Navier—Stokes (RANS) equations, which are combined with one of four turbulence models (Spalart—Allmaras, k—∊, k—∊ renormalization group, and k—ω shear stress transport (SST)) aiming at validation of these models through a comparison of the predictions and the free field experimental measurements for the selected rotor. The computational domain is composed of 4.2×106 cells merged in a structured way, taking care of refinement of the grid near the rotor blade in order to enclose the boundary layer approach. The constant wind condition 7.2 m/s, which is the velocity of the selected experimental data, is considered in all calculations, and only the turbulence model is altered. It is confirmed that it is possible to analyse a HAWT rotor flow field with the RANS equations and that there is good agreement with experimental results, especially when they are combined with the k—ω SST turbulence model.


2021 ◽  
pp. 24-30
Author(s):  
S. V. Kartashev ◽  
◽  
Yu. V. Kozhukhov ◽  

The paper considers the issue of improving the quality of the numerical experiment in the calculation of viscous gas in the flowing part of a low-flow centrifugal compressor stage. The choice of turbulence model in creating a calculation model for calculations by methods of computational fluid dynamics is substantiated. As object of research is chosen low-flow stage with conditional flow coefficient Ф=0,008 and relative width at impeller outlet b2 /D2 =0,0133. The issue of qualitative modeling of friction losses in low-flow stages is of fundamental importance and is directly related to the choice of turbulence model. It is shown that the choice of low-Reynolds turbulence models in the case of unloaded and discontinuous low-flow stages can be made from the main common models (SpalartAllmaras, SST, k-ω) based on the economy of calculations, speed of convergence, solution stability and adequacy of the obtained results. For models with wall functions, the quality of the mesh model and the observance of the dimensionless distance to the wall y+ throughout the calculation domain are particularly important. For highReynolds turbulence models, at values of y+=25...50 on all friction surfaces of the computational domain in the optimal mode of operation, the grid independence of the solution for the entire gas-dynamic characteristic is ensured. It is unacceptable for y+ to fall into the transition region of 4...15 between the viscous sublayer and the region of the logarithmic velocity profile


Author(s):  
D. Scott Holloway ◽  
James H. Leylek

This paper documents the computational investigation of the unsteady rollup and breakdown of a turbulent separated shear layer. This complex phenomenon plays a key role in many applications, such as separated flow at the leading edge of an airfoil at off-design conditions; flow through the tip clearance of a rotor in a gas turbine; flow over the front of an automobile or aircraft carrier; and flow through turbulated passages that are used to cool turbine blades. Computationally, this problem poses a significant challenge in the use of traditional RANS-based turbulence models for the prediction of unsteady flows. To demonstrate this point, a series of 2-D and 3-D unsteady simulations have been performed using a variety of well-known turbulence models, including the “realizable” k-ε model, a differential Reynolds stress model, and a new model developed by the present authors that contains physics that account for the effects of local unsteadiness on turbulence. All simulations are fully converged and grid independent in the unsteady framework. A proven computational methodology is used that takes care of several important aspects, including high-quality meshes (2.5 million finite volumes for 3-D simulations) and a discretization scheme that will minimize the effects of numerical diffusion. To isolate the shear layer breakdown phenomenon, the well-studied flow over a blunt leading edge (Reynolds number based on plate half-thickness of 26,000) is used for validation. Surprisingly, none of the traditional eddy-viscosity or Reynolds stress models are able to predict an unsteady behavior even with modifications in the near-wall treatment, repeated adaption of the mesh, or by adding small random perturbations to the flow field. The newly developed unsteady-based turbulence model is shown to predict some important features of the shear layer rollup and breakdown.


1983 ◽  
Vol 105 (4) ◽  
pp. 439-446 ◽  
Author(s):  
D. R. Boyle ◽  
M. W. Golay

Turbulent flow measurements have been performed in a two-dimensional flow cell which is a 1/15-scale model of the Fast Flux Test Facility nuclear reactor outlet plenum. In a steady water flow, maps of the mean velocity field, turbulence kinetic energy, and Reynolds stress have been obtained using a laser doppler anemometer. The measurements are compared to numerical simulations using both the K–ε and K–σ two-equation turbulence models. A relationship between K–σ and K–ε turbulence models is derived, and the two models are found to be nearly equivalent. The steady-state mean velocity data are predicted well through-out most of the test cell. Calculated spatial distributions of the scalar turbulence quantities are qualitatively similar for both models; however, the predicted distributions do not match the data over major portions of the flow area. The K–σ model provides better estimates of the turbulence quantity magnitudes. The predicted results are highly sensitive to small changes in the turbulence model constants and depend heavily on the levels of inlet turbulence. However, important differences between prediction and measurement cannot be significantly reduced by simple changes to the model’s constants.


2016 ◽  
Vol 2016 ◽  
pp. 1-9
Author(s):  
Masami Suzuki

In designing a wind turbine, the validation of the mathematical model’s result is normally carried out by comparison with wind tunnel experiment data. However, the Reynolds number of the wind tunnel experiment is low, and the flow does not match fully developed turbulence on the leading edge of a wind turbine blade. Therefore, the transition area from laminar to turbulent flow becomes wide under these conditions, and the separation point is difficult to predict using turbulence models. The prediction precision decreases dramatically when working with tip speed ratios less than the maximum power point. This study carries out a steadiness calculation with turbulence model and an unsteadiness calculation with laminar model for a three-blade horizontal axis wind turbine. The validation of the calculations is performed by comparing with experimental results. The power coefficients calculated without turbulence models are in agreement with the experimental data for a tip speed ratio greater than 5.


Author(s):  
Bob Mischo ◽  
Thomas Behr ◽  
Reza S. Abhari

In axial turbine the tip clearance flow occurring in rotor blade rows is responsible for about one third of the aerodynamic losses in the blade row and in many cases is the limiting factor for the blade lifetime. The tip leakage vortex forms when the leaking fluid crosses the gap between the rotor blade tip and the casing from pressure to suction side and rolls up into a vortex on the blade suction side. The flow through the tip gap is both of high velocity and high temperature, with the heat transfer to the blade from the hot fluid being very high in the blade tip area. In order to avoid blade tip burnout and a failure of the machine, blade tip cooling is commonly used. This paper presents the physical study and an improved design of a recessed blade tip for a highly loaded axial turbine rotor blade with application in high pressure axial turbines in aero engine or power generation. With use of three-dimensional Computational Fluid Dynamics (CFD), the flow field near the tip of the blade for different shapes of the recess cavities is investigated. Through better understanding and control of cavity vortical structures, an improved design is presented and the differences to the generic flat tip blade are highlighted. It is observed that by an appropriate profiling of the recess shape, the total tip heat transfer Nusselt Number was significantly reduced, being 15% lower than the flat tip and 7% lower than the baseline recess shape. Experimental results also showed an overall improvement of 0.2% in the one-and-1/2-stage turbine total efficiency with the improved recess design compared to the flat tip case. The CFD analysis conducted on single rotor row configurations predicted a 0.38% total efficiency increase for the rotor equipped with the new recess design compared to the flat tip rotor.


Author(s):  
D. Hasen ◽  
S. Elangovan ◽  
M. Sundararaj ◽  
K.M. Parammasivam

In this study, the effects of different turbulence models on the decay characteristics of round jets were studied. The turbulence models considered for the current study is SST, k-ε, k-ω, RNG kε. For the entire turbulence model mesh density and boundary conditions were mentioned same. By comparing the simulated results with the experiments interesting results were obtained. SST predicts the flow better than the other models in this flow regime.


2008 ◽  
Vol 130 (2) ◽  
Author(s):  
Bob Mischo ◽  
Thomas Behr ◽  
Reza S. Abhari

In axial turbine, the tip clearance flow occurring in rotor blade rows is responsible for about one-third of the aerodynamic losses in the blade row and in many cases is the limiting factor for the blade lifetime. The tip leakage vortex forms when the leaking fluid crosses the gap between the rotor blade tip and the casing from pressure to suction side and rolls up into a vortex on the blade suction side. The flow through the tip gap is both of high velocity and of high temperature, with the heat transfer to the blade from the hot fluid being very high in the blade tip area. In order to avoid blade tip burnout and a failure of the machine, blade tip cooling is commonly used. This paper presents the physical study and an improved design of a recessed blade tip for a highly loaded axial turbine rotor blade with application in high pressure axial turbines in aero engine or power generation. With use of three-dimensional computational fluid dynamics (CFD), the flow field near the tip of the blade for different shapes of the recess cavities is investigated. Through better understanding and control of cavity vortical structures, an improved design is presented and its differences from the generic flat tip blade are highlighted. It is observed that by an appropriate profiling of the recess shape, the total tip heat transfer Nusselt number was significantly reduced, being 15% lower than the flat tip and 7% lower than the base line recess shape. Experimental results also showed an overall improvement of 0.2% in the one-and-a-half-stage turbine total efficiency with the improved recess design compared to the flat tip case. The CFD analysis conducted on single rotor row configurations predicted a 0.38% total efficiency increase for the rotor equipped with the new recess design compared to the flat tip rotor.


Sign in / Sign up

Export Citation Format

Share Document