Computational Modeling of Three-Dimensional Endwall Flow Through a Turbine Rotor Cascade With Strong Secondary Flows

Author(s):  
Y.-H. Ho ◽  
B. Lakshminarayana

A steady, three-dimensional Navier-Stokes solver which utilizes a pressure-based technique for incompressible flows is used to simulate the three-dimensional flow field in a turbine cascade. A new feature of the numerical scheme is the implementation of a second-order plus fourth-order artificial dissipation formulation, which provides a precise control of the numerical dissipation. A low-Reynolds-number form of a two-equation turbulence model is used to account for the turbulence effects. Comparison between the numerical predictions and the experimental data indicates that the numerical model is able to capture most of the complex flow phenomena in the endwall region of a turbine cascade, except the high gradient region in the secondary vortex core. The effects of inlet turbulence intensity and turbulence length scale on secondary vortices, total pressure loss, and turbulence kinetic energy inside the passage are presented and interpreted. It is found that higher turbulence intensity energizes the vortical motions and tends to move the passage vortex away from the endwall. With a larger turbulence length scale the secondary flow inside the passage is reduced. However, the total pressure loss increases due to higher turbulence kinetic energy production.

1996 ◽  
Vol 118 (2) ◽  
pp. 250-261 ◽  
Author(s):  
Y.-H. Ho ◽  
B. Lakshminarayana

A steady, three-dimensional Navier–Stokes solver that utilizes a pressure-based technique for incompressible flows is used to simulate the three-dimensional flow field in a turbine cascade. A new feature of the numerical scheme is the implementation of a second-order plus fourth-order artificial dissipation formulation, which provides a precise control of the numerical dissipation. A low-Reynolds-number form of a two-equation turbulence model is used to account for the turbulence effects. Comparison between the numerical predictions and the experimental data indicates that the numerical model is able to capture most of the complex flow phenomena in the endwall region of a turbine cascade, except the high gradient region in the secondary vortex core. The effects of inlet turbulence intensity and turbulence length scale on secondary vortices, total pressure loss, and turbulence kinetic energy inside the passage are presented and interpreted. It is found that higher turbulence intensity energizes the vortical motions and tends to move the passage vortex away from the endwall. With a larger turbulence length scale, the secondary flow inside the passage is reduced. However, the total pressure loss increases due to higher turbulence kinetic energy production.


Author(s):  
A. Chernobrovkin ◽  
B. Lakshminarayana

A viscous flow solver based on the Runge-Kutta scheme has been modified for the numerical investigation of the aerothermal field due to the leading edge film cooling at a compound angle. An existing code has been modified to incorporate multi-block capabilities. Good agreement with the measured data has been achieved. Results of the numerical investigation have been used to analyze the vortex structure associated with the coolant jet-freestream interaction to understand the contribution of different vortices on the cooling effectiveness and aerothermal losses. Two counter-rotating vortices generated by the interaction between the mainflow and the coolant jet have been found to have a major influence in decreasing the cooling efficiency through strong entrainment of the hot fluid. Numerical simulation was carried out to investigate the influence of the inlet Mach number, inlet turbulence intensity, and length scale on the aerothermal field due to the leading edge film cooling. Variation of the inlet Mach number leads to a minor modification of the cooling effectiveness, and this is predominantly caused by the modified pressure gradient. Increased turbulence intensity has profound effect on the cooling near the leading edge. Adiabatic effectiveness downstream of the second row of coolant holes is less sensitive to a change in turbulence intensity. Results of the numerical simulation indicate that the turbulence length scale has a significant effect on the accuracy of the numerical prediction of film cooling. Not only the inlet turbulence intensity but also the turbulence length scale should be accurately set to achieve a reliable numerical prediction of the heat and mass transfer due to film cooling.


2018 ◽  
Vol 3 (2) ◽  
pp. 533-543 ◽  
Author(s):  
Mark Kelly

Abstract. In wind energy, the effect of turbulence upon turbines is typically simulated using wind “input” time series based on turbulence spectra. The velocity components' spectra are characterized by the amplitude of turbulent fluctuations, as well as the length scale corresponding to the dominant eddies. Following the IEC standard, turbine load calculations commonly involve use of the Mann spectral-tensor model to generate time series of the turbulent three-dimensional velocity field. In practice, this spectral-tensor model is employed by adjusting its three parameters: the dominant turbulence length scale LMM (peak length scale of an undistorted isotropic velocity spectrum), the rate of dissipation of turbulent kinetic energy ε, and the turbulent eddy-lifetime (anisotropy) parameter Γ. Deviation from “ideal” neutral sheared turbulence – i.e., for non-zero heat flux and/or heights above the surface layer – is, in effect, captured by setting these parameters according to observations. Previously, site-specific {LMM,ε,Γ} values were obtainable through fits to measured three-dimensional velocity component spectra recorded with sample rates resolving the inertial range of turbulence (≳1 Hz); however, this is not feasible in most industrial wind energy projects, which lack multi-dimensional sonic anemometers and employ loggers that record measurements averaged over intervals of minutes. Here a form is derived for the shear dependence implied by the eddy-lifetime prescription within the Mann spectral-tensor model, which leads to derivation of useful forms of the turbulence length scale. Subsequently it is shown how LMM can be calculated from commonly measured site-specific atmospheric parameters, namely mean wind shear (dU∕dz) and standard deviation of streamwise fluctuations (σu). The derived LMM can be obtained from standard (10 min average) cup anemometer measurements, in contrast with an earlier form based on friction velocity. The new form is tested across several different conditions and sites, and it is found to be more robust and accurate than estimates relying on friction velocity observations. Assumptions behind the derivations are also tested, giving new insight into rapid-distortion theory and eddy-lifetime modeling – and application – within the atmospheric boundary layer. The work herein further shows that distributions of turbulence length scale, obtained using the new form with typical measurements, compare well with distributions P(LMM) obtained by fitting to spectra from research-grade sonic anemometer measurements for the various flow regimes and sites analyzed. The new form is thus motivated by and amenable to site-specific probabilistic loads characterization.


Author(s):  
Changfu You ◽  
Haiying Qi ◽  
Xuchang Xu

Effect of turbulence on drag force in gas-particle two-phase flow had been investigated using numerical simulation. In order to select an accurate turbulence model, some promising models, such as standard k-ε model, RNG k-ε model and Realizable k-ε model, had been examined through calculating the flow over a backward-facing step. RNG k-ε model performing better than others had been used to simulate the turbulence flow over a spherical particle. In computation, the turbulence intensity was ranged from 10% to 80%, and the turbulence length scale from 10−5m to 4m. Results show that the turbulence length scale had a small effect on the drag force, except at small length scale. Comparing with the drag on a particle in laminar flow, the turbulence intensity enhances it comparatively, especially at small particle Reynolds numbers, which differs from the previous publications.


Author(s):  
F. Mumic ◽  
B. Sunden

In the present work, a numerical study has been performed to simulate the effect of free-stream turbulence, length scale and variations in rotational speed of the rotor on heat transfer and fluid flow for a transonic high-pressure turbine stage with tip clearance. The stator and rotor rows interact via a mixing plane, which allows the stage to be computed in a steady manner. The focus is on turbine aerodynamics and heat transfer behavior at the mid-span location, and at the rotor tip and casing region. The results of the fully 3D CFD simulations are compared with experimental results available for the so-called MT1 turbine stage. The predicted heat transfer and static pressure distributions show reasonable agreement with the experimental data. In general, the local Nusselt number increases, at the same turbulence length scale, as the turbulence intensity increases, and the location of the suction side boundary layer transition moves upstream towards the blade leading edge. Comparison of the different length scales at the same turbulence intensity shows that the stagnation heat transfer was significantly increased as the length scale increased. However, the length scale evidenced no significant effects on blade tip or rotor casing heat transfer. Also, the results presented in this paper show that the rotational speed in addition to the turbulence intensity and length scale has an important contribution to the turbine blade aerodynamics and heat transfer.


Author(s):  
Zhigang Li ◽  
Bo Bai ◽  
Luxuan Liu ◽  
Jun Li ◽  
Shuo Mao ◽  
...  

Abstract In gas turbine engines, the first-stage vanes usually suffer harsh incoming flow conditions from the combustor with high pressure, high temperature and high turbulence. The combustor-generated high freestream turbulence and strong secondary flows in a gas turbine vane passage have been reported to augment the endwall thermal load significantly. This paper presents a detailed numerical study on the effects of high freestream turbulence intensity, turbulence length scale, and exit Reynolds number on the endwall secondary flow pattern and heat transfer distribution of a transonic linear turbine vane passage at realistic engine Mach numbers, with a flat endwall no cooling. Numerical simulations were conducted at a range of different operation conditions: six freestream turbulence intensities (Tu = 1%, 5%, 10%, 13%, 16% and 20%), six turbulence length scales (normalized by the vane pitch of Λ/P = 0.01, 0.04, 0.07, 0.12, 0.24, 0.36), and three exit isentropic Mach number (Maex = 0.6, 0.85 and 1.02 corresponding exit Reynolds number Reex = 1.1 × 106, 1.7 × 106 and 2.2 × 106, respectively, based on the vane chord). Detailed comparisons were presented for endwall heat transfer coefficient distribution, endwall secondary flow field at different operation conditions, while paying special attention to the link between endwall thermal load patterns and the secondary flow structures. Results show that the freestream turbulence intensity and length scale have a significant influence on the endwall secondary flow field, but the influence of the exit Reynolds number is very weak. The Nusselt number patterns for the higher turbulence intensities (Tu = 16%, 20%) appear to be less affected by the endwall secondary flows than the lower turbulence cases. The thermal load distribution in the arc region around the vane leading edge and the banded region along the vane pressure side are influenced most strongly by the freestream turbulence intensity. In general, the higher freestream turbulence intensities make the vane endwall thermal load more uniform. The Nusselt number distribution is only weakly affected by the turbulence length scale when Λ/P is larger than 0.04. The heat transfer level appears to have a significant uniform augmentation over the whole endwall region with the increasing Maex. The endwall thermal load distribution is classified into four typical regions, and the effects of freestream turbulence, exit Reynolds number in each region were discussed in detail.


2018 ◽  
Author(s):  
Mark Kelly

Abstract. In wind energy, the the effect of turbulence upon turbines is typically simulated using wind ‘input’ time-series based on turbulence spectra. The velocity components’ spectra are characterized by the amplitude of turbulent fluctuations, as well as the length scale corresponding to the dominant eddies. Following the IEC standard, turbine loads calculations commonly involve use of the Mann spectral-tensor model to generate timeseries of the turbulent three-dimensional velocity field. In practice, this spectral-tensor model is employed by adjusting its three parameters: the dominant turbulence length scale LMM (peak length scale of an undistorted isotropic velocity spectrum), the rate of dissipation of turbulent kinetic energy ε, and the turbulent eddy-lifetime (anisotropy) parameter Γ. Deviation from ‘ideal’ neutral sheared turbulence – i.e. for non-zero heat flux and/or heights above the surface layer – is, in effect, captured by setting these parameters according to observations. Previously, site-specific {LMM, ε, Γ} were obtainable through fits to measured three-dimensional velocity component spectra recorded with sample rates resolving the inertial range of turbulence (>~ 1 Hz); however, this is not feasible in most industrial wind energy projects, which lack multi-dimensional sonic anemometers and employ loggers that record measurements averaged over intervals of minutes. Here a form is derived for the shear dependence implied by the eddy-lifetime prescription within the Mann spectral-tensor model, which leads to derivation of useful forms of the turbulence length scale. Subsequently it is shown how LMM can be calculated from commonly-measured site-specific atmospheric parameters, namely mean wind shear (dU/dz) and standard deviation of streamwise fluctuations (σu). The derived LMM can be obtained from standard (10-minute average) cup anemometer measurements, in contrast with forms based on friction velocity. The new form is tested across several different conditions and sites, and is found to be more robust and accurate than estimates relying on friction velocity observations. Assumptions behind the derivations are also tested, giving new insight into rapid-distortion theory and eddy-lifetime modelling – and application – within the atmospheric boundary layer. The work herein further shows that distributions of turbulence length scale, obtained using the new form with typical measurements, compare well with distributions P(LMM) obtained by fitting to spectra from research-grade sonic anemometer measurements for the various flow regimes and sites analyzed. The new form is thus motivated by and amenable to site-specific probabilistic loads characterization.


Sign in / Sign up

Export Citation Format

Share Document