Numerical Investigation of Anisotropic Turbulence Effects Around a Film Cooled Turbine Vane

Author(s):  
Dieter E. Bohn ◽  
Christian Tu¨mmers

In this paper the 3D flow around a turbine vane with showerhead cooling is simulated with the anisotropic cubic-eddy-viscosity k-ε turbulence model of Craft et. al. The results are analyzed in detail and compared to calculations performed with the isotropic algebraic turbulence model by Baldwin & Lomax and the isotropic Low-Reynolds k-ε-model of Launder et. al. for the same test case. The computational domain consists of the coolant supply (plenum), the ejection holes and the main flow region around the vane. Periodic boundary conditions have been used in the radial direction. Thus, endwall effects have been excluded. The numerical investigations focus on the influence of the anisotropic effects in the flow field. The flow conditions are taken from experimental investigations conducted by other authors and the results have been documented as a test case for numerical calculations of ejection flow phenomena. The results show a similar distribution of the predicted total pressure loss in the kidney vortex and downstream the ejection holes for the cubic-eddy viscosity model in comparison to the other turbulence models. Furthermore it is shown, that the influence of the cooling jets on the main flow, predicted by the two-equation models of Launder et. al. and Craft. et al., seems to be slightly higher compared to the algebraic model of Baldwin & Lomax.

Author(s):  
Karsten Tawackolian ◽  
Martin Kriegel

AbstractThis study looks to find a suitable turbulence model for calculating pressure losses of ventilation components. In building ventilation, the most relevant Reynolds number range is between 3×104 and 6×105, depending on the duct dimensions and airflow rates. Pressure loss coefficients can increase considerably for some components at Reynolds numbers below 2×105. An initial survey of popular turbulence models was conducted for a selected test case of a bend with such a strong Reynolds number dependence. Most of the turbulence models failed in reproducing this dependence and predicted curve progressions that were too flat and only applicable for higher Reynolds numbers. Viscous effects near walls played an important role in the present simulations. In turbulence modelling, near-wall damping functions are used to account for this influence. A model that implements near-wall modelling is the lag elliptic blending k-ε model. This model gave reasonable predictions for pressure loss coefficients at lower Reynolds numbers. Another example is the low Reynolds number k-ε turbulence model of Wilcox (LRN). The modification uses damping functions and was initially developed for simulating profiles such as aircraft wings. It has not been widely used for internal flows such as air duct flows. Based on selected reference cases, the three closure coefficients of the LRN model were adapted in this work to simulate ventilation components. Improved predictions were obtained with new coefficients (LRNM model). This underlined that low Reynolds number effects are relevant in ventilation ductworks and give first insights for suitable turbulence models for this application. Both the lag elliptic blending model and the modified LRNM model predicted the pressure losses relatively well for the test case where the other tested models failed.


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.


2004 ◽  
Vol 126 (5) ◽  
pp. 844-850 ◽  
Author(s):  
Khaled S. Abdol-Hamid ◽  
S. Paul Pao ◽  
Steven J. Massey ◽  
Alaa Elmiligui

It is well known that the two-equation turbulence models under-predict mixing in the shear layer for high temperature jet flows. These turbulence models were developed and calibrated for room temperature, low Mach number, and plane mixing layer flows. In the present study, four existing modifications to the two-equation turbulence model are implemented in PAB3D and their effect is assessed for high temperature jet flows. In addition, a new temperature gradient correction to the eddy viscosity term is tested and calibrated. The new model was found to be in the best agreement with experimental data for subsonic and supersonic jet flows at both low and high temperatures.


Author(s):  
Chen Fu ◽  
C Patrick Bounds ◽  
Christian Selent ◽  
Mesbah Uddin

The characterization of a racecar’s aerodynamic behavior at various yaw and pitch configurations has always been an integral part of its on-track performance evaluation in terms of lap time predictions. Although computational fluid dynamics has emerged as the ubiquitous tool in motorsports industry, a clarity is still lacking about the prediction veracity dependence on the choice of turbulence models, which is central to the prediction variability and unreliability for the Reynolds Averaged Navier–Stokes simulations, which is by far the most widely used computational fluid dynamics methodology in this industry. Subsequently, this paper presents a comprehensive assessment of three commonly used eddy viscosity turbulence models, namely, the realizable [Formula: see text] (RKE), Abe–Kondoh–Nagano [Formula: see text], and shear stress transport [Formula: see text], in predicting the aerodynamic characteristics of a full-scale NASCAR Monster Energy Cup racecar under various yaw and pitch configurations, which was never been explored before. The simulations are conducted using the steady Reynolds Averaged Navier–Stokes approach with unstructured trimmer cells. The tested yaw and pitch configurations were chosen in consultation with the race teams such that they reflect true representations of the racecar orientations during cornering, braking, and accelerating scenarios. The study reiterated that the prediction discrepancies between the turbulence models are mainly due to the differences in the predictions of flow recirculation and separation, caused by the individual model’s effectiveness in capturing the evolution of adverse pressure gradient flows, and predicting the onset of separation and subsequent reattachment (if there be any). This paper showed that the prediction discrepancies are linked to the computation of the turbulent eddy viscosity in the separated flow region, and using flow-visualizations identified the areas on the car body which are critical to this analysis. In terms of racecar aerodynamic performance parameter predictions, it can be reasonably argued that, excluding the prediction of the %Front prediction, shear stress transport is the best choice between the three tested models for stock-car type racecar Reynolds Averaged Navier–Stokes computational fluid dynamics simulations as it is the only model that predicted directionally correct changes of all aerodynamic parameters as the racecar is either yawed from the 0° to 3° or pitched from a high splitter-ground clearance to a low one. Furthermore, the magnitude of the shear stress transport predicted delta force coefficients also agreed reasonably well with test results.


2008 ◽  
Vol 130 (9) ◽  
Author(s):  
G. Dufour ◽  
J.-B. Cazalbou ◽  
X. Carbonneau ◽  
P. Chassaing

Rotation and curvature (RC) effects on turbulence are expected to impact losses and flow structure in turbomachines. This paper examines two recent eddy-viscosity-model corrections devised to account for these effects: the Spalart and Shur (1997, “On the Sensitization of Turbulence Models to Rotation and Curvature,” Aerosp. Sci. Technol., 1(5), pp. 297–302) correction to the model of Spalart and Allmaras (1994, “A One-Equation Turbulence Model for Aerodynamic Flows,” Rech. Aerosp., 1, pp. 5–21) and the correction of Cazalbou et al. (2005, “Two-Equation Modeling of Turbulent Rotating Flows,” Phys. Fluids., 17, p. 055110) to the (k,ϵ) model. The method of verification and validation is applied to assess the impact of these corrections on the computation of a centrifugal-compressor test case. First, a review of RC effects on turbulence as they apply to centrifugal compressors is made. The two corrected models are then presented. Second, the Radiver open test case (Ziegler K. U., Gallus, H. E., and Niehuis R., 2003, “A Study on Impeller Diffuser Interaction Part 1: Influence on the Performance,” ASME J. Turbomach, 125, pp. 173–182) is used as a basis for the assessment of the two corrections. After a physical-consistency analysis, the Richardson extrapolation is applied to quantify the numerical errors involved in all the calculations. Finally, experimental data are used to perform validation for both global and local predictions. The consistency analysis shows that both corrections lead to significant changes in the turbulent field, in perfect agreement with the underlying theoretical considerations. The uncertainty analysis shows that the predictions of the global performances are more sensitive to grid refinement than they are to RC turbulence modeling. However, the opposite conclusion is drawn with regard to the prediction of some local flow properties: Improvements are obtained with the RC corrections, the best results being observed for the RC-corrected (k,ϵ) model.


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):  
Xin Deng ◽  
Harrison Gates ◽  
Brian Weaver ◽  
Houston Wood ◽  
Roger Fittro

Oil-lubricated bearings are widely used in high speed rotating machines such as those found in the aerospace and automotive industries. However, environmental issues and risk-averse operations are resulting in the removal of oil and the replacement of all sealed oil bearings with reliable water-lubricated bearings. Due to the different fluid properties between oil and water, the low viscosity of water increases Reynolds numbers drastically and therefore makes water-lubricated bearings prone to turbulence effects. This requires finer meshes when compared to oil-lubricated bearings as the low-viscosity fluid produces a very thin lubricant film. Analyzing water-lubricated bearings can also produce convergence and accuracy issues in traditional oil-based analysis codes. Thermal deformation largely affects oil-lubricated bearings, while having limited effects on water lubrication; mechanical deformation largely affects water lubrication, while its effects are typically lower than thermal deformation with oil. One common turbulence model used in these analysis tools is the eddy-viscosity model. Eddy-viscosity depends on the wall shear stress, therefore effective wall shear stress modeling is necessary in determining an appropriate turbulence model. Improving the accuracy and efficiency of modeling approaches for eddy-viscosity in turbulence models is of great importance. Therefore, the purpose of this study is to perform mesh refinement for water-lubricated bearings based on methodologies of eddy-viscosity modeling to improve their accuracy. According to Szeri [1], εm/v for the Boussinesq hypothesis is given by Reichardt’s formula. Fitting the velocity profile with experiments having a y+ in the range of 0–1,000 results in Ng-optimized Reichardt’s constants k = 0.4 and δ+ = 10.7. He clearly states that for y+ > 1000 theoretical predictions and experiments have a greater variance. Armentrout and others [2] developed an equation for δ+ as a function of the pivot Reynolds number, which they validated with CFD simulations. The definition of y+ can be used to approximate the first layer thickness calculated for a uniform mesh. Together with Armentrout’s equation, the number of required elements across the film thickness can be obtained. For typical turbulence models, the y+ must be within a certain range to be accurate. On the condition that the y+ is fixed to that of a standard oil bearing for which an oil bearing code was validated, the number of elements across the film thickness and coefficients used in the eddy-viscosity equation can be adjusted to allow for convergence with other fluids other than that which the traditional oil bearing code was designed for. In this study, the number of required elements across the film for improved prediction quality was calculated based on the proposed eddy-viscosity model mesh correction from the known literature. A comparison between water lubrication using the parameter correction and oil lubrication was also made. The results of this study could aid in improving future designs and models of water-lubricated bearings.


Author(s):  
Jian Wang ◽  
Yong Wang ◽  
Houlin Liu ◽  
Haoqin Huang ◽  
Linglin Jiang

Purpose – The purpose of this paper is to study the unsteady caivitating flows in centrifugal pump, especially for improving the turbulence model to obtain highly resolution results-capable of predicting the cavitation inception, shedding off and collapse procedures. Design/methodology/approach – Both numerical simulations and experimental visualizations were performed in the present paper. An improved RCD turbulence models was proposed by considering three corrected methods: the rotating corrected method, the compressible corrected method and the turbulent viscosity corrected method. Unsteady RANS computations were conducted to compare with the experiments. Findings – The comparison of pump cavitation performance showed that the RCD turbulence model obtained better performance both in non-cavitation and cavitation conditions. The visualization of the cavitation evolution was recorded to validate the unsteady simulations. Good agreement was noticed between calculations and visualizations. It is indicated the RCD model can successfully capture the bubbles detachment and collapse at the rear of the cavity region, since it effectively reduces the eddy viscosity in the multiphase region of liquid and vapor. Furthermore, the eddy viscosity, the instantaneous pressure and density distribution were investigated. The effectiveness of the compressibility was found. Meanwhile, the influence of the rotating corrected method on prediction was explored. It is found that the RCD model solved more unsteady flow characteristics. Originality/value – The current work presented a turbulence model which was much more suitable for predicting the cavitating flow in centrifugal pump.


2021 ◽  
Vol 5 ◽  
pp. 28-38
Author(s):  
Wenqiang Zhang ◽  
Mehdi Vahdati

Experimental studies have shown that tip injection upstream of the rotor can extend its operational range when subjected to circumferential inlet distortion. Typically, injectors are placed uniformly around the annulus. However, such arrangement consumes a large amount of high-pressure air and decreases the overall efficiency of the compression system. The aim of this paper is to minimise the amount of the injected air by determining the most effective circumferential location for the injector. In this study, NASA stage 35 was used as the test case. The experiment was conducted with a circumferential total pressure distortion of 120 degrees. In the first part of this paper, numerical simulations were compared against the experimental data and good match was obtained. In the second part, tip injection at three different positions were tested: the clean flow region (Position 1), the distorted region (Position 2) and the border between the clean and distorted regions (Position 3). It was found that a mild injection (0.66% of the main flow) at Position 2 and Position 3 can extend the stall margin by 1.8% and 2.7%, respectively. No obvious improvement was observed for the injection at Position 1. With a larger injection of 1.5% of main flow at Position 3, the stall margin improved further with no efficiency loss.


Sign in / Sign up

Export Citation Format

Share Document