V. Nagabhushana Rao

Numerical Investigation of Contrasting Flow Physics in Different Zones of a High-Lift Low-Pressure Turbine Blade

Using a range of high-fidelity large eddy simulations (LES), the contrasting flow physics on the suction surface, pressure surface, and endwalls of a low-pressure turbine (LPT) blade (T106A) was studied. The current paper attempts to provide an improved understanding of the flow physics over these three zones under the influence of different inflow boundary conditions. These include: (a) the effect of wakes at low and high turbulence intensity on the flow at midspan and (b) the impact of the state of the incoming boundary layer on endwall flow features. On the suction surface, the pressure fluctuations on the aft portion significantly reduced at high freestream turbulence (FST). The instantaneous flow features revealed that this reduction at high FST (HF) is due to the dominance of “streak-based” transition over the “Kelvin–Helmholtz” (KH) based transition. Also, the transition mechanisms observed over the turbine blade were largely similar to those on a flat plate subjected to pressure gradients. On pressure surface, elongated vortices were observed at low FST (LF). The possibility of the coexistence of both the Gortler instability and the severe straining of the wakes in the formation of these elongated vortices was suggested. While this was true for the cases under low turbulence levels, the elongated vortices vanished at higher levels of background turbulence. At endwalls, the effect of the state of the incoming boundary layer on flow features has been demonstrated. The loss cores corresponding to the passage vortex and trailing shed vortex were moved farther from the endwall with a turbulent boundary layer (TBL) when compared to an incoming laminar boundary layer (LBL). Multiple horse-shoe vortices, which constantly moved toward the leading edge due to a low-frequency unstable mechanism, were captured.

Internal instability of thin liquid sheets

Linear stability analysis of an inviscid liquid sheet with different velocity profiles across its thickness is reported. The velocity profiles for which there is a progressive increase or decrease in velocities between the two interfaces are demonstrated to be inherently unstable even in the absence of the destabilizing aerodynamic shear at the liquid-gas interfaces. Compared to a flat velocity profile, a linear or a parabolic profile, symmetric at the center line of the sheet reduced both the maximum growth rate and the wavelength range over which the waves grow. The convective acceleration from the velocity gradient is found to stabilize longer waves while the growth of shorter waves is hampered by the combined effect of the surface tension and a decrease in the interface velocity between gas and liquid media. The wave forms are dominantly sinuous for symmetric velocity profiles; however, with larger velocity gradients the dilatational modes are observed. The inherent instability of liquid sheets with a...

Simulation studies on multi-mode heat transfer from an open cavity with a flush-mounted discrete heat source

Prominent results of a simulation study on conjugate convection with surface radiation from an open cavity with a traversable flush mounted discrete heat source in the left wall are presented in this paper. The open cavity is considered to be of fixed height but with varying spacing between the legs. The position of the heat source is varied along the left leg of the cavity. The governing equations for temperature distribution along the cavity are obtained by making energy balance between heat generated, conducted, convected and radiated. Radiation terms are tackled using radiosity-irradiation formulation, while the view factors, therein, are evaluated using the crossed-string method of Hottel. The resulting non-linear partial differential equations are converted into algebraic form using finite difference formulation and are subsequently solved by Gauss–Seidel iterative technique. An optimum grid system comprising 111 grids along the legs of the cavity, with 30 grids in the heat source and 31 grids across the cavity has been used. The effects of various parameters, such as surface emissivity, convection heat transfer coefficient, aspect ratio and thermal conductivity on the important results, including local temperature distribution along the cavity, peak temperature in the left and right legs of the cavity and relative contributions of convection and radiation to heat dissipation in the cavity, are studied in great detail.

Investigation of Wake Induced Transition in Low-Pressure Turbines Using Large Eddy Simulation

Modern ‘high-lift’ blade designs incorporated into the low pressure turbine (LPT) of aero-engines typically exhibit a separation bubble on the suction surface of the airfoil. The size of the bubble and the loss it generates is governed by the transition process in the separated shear layer. However, the wakes shed by the upstream blade rows, the turbulent fluctuations in the free-stream and the roughness over the blade complicates the transition process.The current paper numerically investigates the transition of a separated shear layer over a flat plate with an elliptic leading edge using large eddy simulations (LES). The upper wall of the test section is inviscid and specifically contoured to impose a streamwise pressure distribution over the flat plate to simulate the suction surface of a LPT blade. The influences of free-stream turbulence (FST), periodic wake passing and streamwise pressure distribution (blade loading) are considered.The simulations were carried out at a Reynolds number of 83,000 based on the length of the flat plate (S0 = 0.5m) and the velocity at the nominal trailing edge (UTE ∼ 2.55 m/s). A high turbulence intensity of 4% and a dimensionless wake passing frequency (fr = fwakeS0/UTE, where fwake is the dimensional wake frequency) of 0.84 is chosen for the study. Two different distributions representative of a ‘high-lift’ and an ‘ultra-high-lift’ turbine blade are examined. An in-house, high order, flow solver is used for the Large Eddy Simulations (LES). The Variational Multi-scale approach is used to account for the sub-grid scale stresses.Results obtained from the current LES compare favorably with the extensive experimental data previously obtained for the test cases considered. The LES results are then used to further explore the flow physics involved in the transition process, in particular the role of Klebanoff streaks and their influence on performance. The additional effect of surface roughness of the blade has also been studied for one of the blade loadings. The benefit that roughness can offer for highly loaded turbine blades is demonstrated.Copyright

Differential Equation Specification of Integral Turbulence Length Scales

A Hamilton–Jacobi differential equation is used to naturally and smoothly (via Dirichlet boundary conditions) set turbulence length scales in separated flow regions based on traditional expected length scales. Such zones occur for example in rim-seals. The approach is investigated using two test cases, flow over a cylinder at a Reynolds number of 140,000 and flow over a rectangular cavity at a Reynolds number of 50,000. The Nee–Kovasznay turbulence model is investigated using this approach. Predicted drag coefficients for the cylinder test-case show significant (15%) improvement over standard steady RANS and are comparable with URANS results. The mean flow-field also shows a significant improvement over URANS. The error in re-attachment length is improved by 180% compared with the steady RANS k-ω model. The wake velocity profile at a location downstream shows improvement and the URANS profile is inaccurate in comparison. For the cavity case, the HJ–NK approach is generally comparable with the other RANS models for measured velocity profiles. Predicted drag coefficients are compared with large eddy simulation. The new approach shows a 20–30% improvement in predicted drag coefficients compared with standard one and two equation RANS models. The shape of the recirculation region within the cavity is also much improved.

Role of axisymmetric axial velocity gradients in nonlinear disintegration of liquid sheets

The nonlinear growth of sinuous and dilatational disturbances on thin inviscid liquid sheets with a velocity profile across the thickness is determined. Variations of the velocity profile in the liquid sheet, with maximum velocity at the centre linearly decreasing to a minimum at the interfaces, show an optimum value of the velocity gradient, below which disintegration of the liquid sheet is accelerated. The ligaments are formed in local dilatational and global dilatational modes and are asymmetric, unlike the symmetric ligaments formed with a liquid sheet issuing at constant velocity. The presence of velocity gradients is seen to yield smaller ligaments when breakup occurs in the local dilatational mode.