scholarly journals Numerical study of the vortex breakdown and vortex reconnection in the flow path of high-pressure water turbine

2021 ◽  
Vol 2088 (1) ◽  
pp. 012040
Author(s):  
A V Sentyabov ◽  
D V Platonov ◽  
A V Minakov ◽  
A S Lobasov

Abstract The paper presents a study of the instability of the precessing vortex core in the model of the draft tube of a hydraulic turbine. The study was carried out using numerical modeling using various approaches: URANS, RSM, LES. The best agreement with the experimental data was shown by the RSM and LES methods with the modelling of the runner rotation by the sliding mesh method. In the regime under consideration, the precessing vortex rope is subject to instability, which leads to reconnection of its turns and the formation of an isolated vortex ring. Reconnection of the vortex core leads to aperiodic and intense pressure fluctuations recorded on the diffuser wall.

Author(s):  
Danielle Mason ◽  
Sean Clees ◽  
Mark Frederick ◽  
Jacqueline O’Connor

Abstract Many industrial combustion systems, especially power generation gas turbines, use fuel-lean combustion to reduce NOx emissions. However, these systems are highly susceptible to combustion instability, the coupling between combustor acoustics and heat release rate oscillations of the flame. It has been shown in previous work by the authors that a precessing vortex core (PVC) can suppress shear layer receptivity to external perturbations, reducing the potential for thermoacoustic coupling. The goal of this study is to understand the effect of combustor exit boundary condition on the flow structure of a swirling jet to increase fundamental understanding of how combustor design impacts PVC dynamics. The swirling jet is generated with a radial-entry, variable-angle swirler, and a quartz cylinder is fixed on the dump plane for confinement. Combustor exit constriction plates of different diameters are used to determine the impact of exit boundary condition on the flow field. Particle image velocimetry (PIV) is used to capture the velocity field inside the combustor. Spectral proper orthogonal decomposition, a frequency-resolved eigenvalue decomposition that can identify energetic structures in the flow, is implemented to identify the PVC at each condition in both energy and frequency space. We find that exit boundary diameter affects both the structure of the flow and the dynamics of the PVC. Higher levels of constriction (smaller diameters) force the downstream stagnation point of the vortex breakdown bubble upstream, resulting in greater divergence of the swirling jet. Further, as the exit diameter decreases, the PVC becomes less energetic and less spatially defined. Despite these changes in the base flow and PVC coherence, the PVC frequency is not altered by the exit boundary constriction. These trends will help inform our understanding of the impact of boundary conditions on both static and dynamic flame stability.


2016 ◽  
Vol 792 ◽  
pp. 620-657 ◽  
Author(s):  
O. Tammisola ◽  
M. P. Juniper

The large-scale coherent motions in a realistic swirl fuel-injector geometry are analysed by direct numerical simulations (DNS), proper orthogonal decomposition (POD), and linear global modes. The aim is to identify the origin of instability in this turbulent flow in a complex internal geometry. The flow field in the nonlinear simulation is highly turbulent, but with a distinguishable coherent structure: the precessing vortex core (a spiralling mode). The most energetic POD mode pair is identified as the precessing vortex core. By analysing the fast Fourier transform (FFT) of the time coefficients of the POD modes, we conclude that the first four POD modes contain the coherent fluctuations. The remaining POD modes (incoherent fluctuations) are used to form a turbulent viscosity field, using the Newtonian eddy model. The turbulence sets in from convective shear layer instabilities even before the nonlinear flow reaches the other end of the domain, indicating that equilibrium solutions of the Navier–Stokes are never observed. Linear global modes are computed around the mean flow from DNS, applying the turbulent viscosity extracted from POD modes. A slightly stable discrete $m=1$ eigenmode is found, well separated from the continuous spectrum, in very good agreement with the POD mode shape and frequency. The structural sensitivity of the precessing vortex core is located upstream of the central recirculation zone, identifying it as a spiral vortex breakdown instability in the nozzle. Furthermore, the structural sensitivity indicates that the dominant instability mechanism is the Kelvin–Helmholtz instability at the inflection point forming near vortex breakdown. Adjoint modes are strong in the shear layer along the whole extent of the nozzle, showing that the optimal initial condition for the global mode is localized in the shear layer. We analyse the qualitative influence of turbulent dissipation in the stability problem (eddy viscosity) on the eigenmodes by comparing them to eigenmodes computed without eddy viscosity. The results show that the eddy viscosity improves the complex frequency and shape of global modes around the fuel-injector mean flow, while a qualitative wavemaker position can be obtained with or without turbulent dissipation, in agreement with previous studies. This study shows how sensitivity analysis can identify which parts of the flow in a complex geometry need to be altered in order to change its hydrodynamic stability characteristics.


Author(s):  
Arnab Mukherjee ◽  
Nishanth Muthichur ◽  
Chaitali More ◽  
Saarthak Gupta ◽  
Santosh Hemchandra

Abstract The precessing vortex core (PVC) phenomenon in swirling jets is a helical instability in the flow driven by the coherent precession of the vortex breakdown bubble (VBB) around the flow axis, resulting in the helical rollup of the shear layer. This instabilitty is driven mainly by flow processes in the region upstream of the VBB. Centerbodies, commonly employed in combustor nozzles create a central wake recirculation zone (CWRZ) that can interfere with VBB precession and hence suppress the PVC. We study this phenomenon in a swirl nozzle with a centerbody whose end face is flush with the nozzle exit plane, using large eddy simulations (LES) and linear hydrodynamic stability analysis for flow Reynolds numbers Re = 48,767 and 82,751, based on nozzle exit diameter and bulk flow velocity. For one of the Re = 82,751 cases the centerbody end face diameter is halved resulting in the onset of coherent VBB precession. Linear stability analysis reveals a marginally unstable mode in this case. The same mode is found to be stable in the nominal cases. Structural sensitivity analysis for these two cases, shows that the VBB precession eigenmode is sensitive to changes in the time averaged flow in the VBB-CWRZ merger region. This suggests that the reduction in CWRZ length due to halving the centerbody end face diameter is the reason for the onset of VBB precession. These results suggest that in general, spatial separation between the CWRZ and VBB can result in the onset of VBB precession and the emergence of PVC oscillations in swirl flows.


2021 ◽  
Vol 143 (5) ◽  
Author(s):  
Arnab Mukherjee ◽  
Nishanth Muthichur ◽  
Chaitali More ◽  
Saarthak Gupta ◽  
Santosh Hemchandra

Abstract The precessing vortex core (PVC) phenomenon in swirling jets is a helical instability in the flow driven by the coherent precession of the vortex breakdown bubble (VBB) around the flow axis, resulting in the helical rollup of the shear layer. This instability is driven by flow processes in the region upstream of the VBB. Centerbodies, commonly employed in combustor nozzles, create a centerbody wake recirculation zone (CWRZ) that can interfere with VBB precession and hence suppress the PVC. We study this phenomenon in a swirl nozzle with a centerbody whose end face is flush with the nozzle exit plane, using large eddy simulations (LES) and linear hydrodynamic stability analysis for flow Reynolds numbers Re = 48,767 and 82,751, based on nozzle exit diameter and bulk flow velocity. For one of the Re = 82,751 cases, the centerbody end face diameter is halved, resulting in the onset of coherent VBB precession. Linear stability analysis reveals a marginally unstable mode in this case. The same mode is found to be stable in the nominal cases. Structural sensitivity analysis shows that the VBB precession eigenmode is sensitive to changes in the time-averaged flow in the VBB-CWRZ merger region. This suggests that the reduction in CWRZ length due to halving the centerbody end face diameter is the reason for the onset of VBB precession. These results suggest that in general, spatial separation between the CWRZ and VBB can result in the onset of VBB precession and the emergence of PVC oscillations in flows with swirl.


2001 ◽  
Vol 123 (3) ◽  
pp. 466-474 ◽  
Author(s):  
F. Shi ◽  
H. Tsukamoto

Two-dimensional and three-dimensional, unsteady state Reynolds-averaged Navier-Stokes (RANS) equations with standard k-ε turbulence models were solved within an entire stage of a diffuser pump to investigate pressure fluctuations due to the interaction between impeller and diffuser vanes. A complete solution of transient flows due to the interaction between components in the whole pump without approximating the blade count ratio of impeller to diffuser was obtained by employing an Arbitrary Sliding Mesh. The unsteady numerical results were compared with experimental data and values calculated by the singularity method. As a result of the present study, the Navier-Stokes code with the k-ε model is found to be capable of predicting pressure fluctuations in the diffuser. Furthermore, extensive numerical studies identified sources contributing to the pressure fluctuations in the diffuser, and helped to understand the mechanism of impeller-diffuser interactions in the diffuser pump.


2005 ◽  
Vol 127 (2) ◽  
pp. 363-371 ◽  
Author(s):  
Kitano Majidi

Computational fluid dynamics (CFD) analysis has been used to solve the unsteady three-dimensional viscous flow in the entire impeller and volute casing of a centrifugal pump. The results of the calculations are used to predict the impeller/volute interaction and to obtain the unsteady pressure distribution in the impeller and volute casing. The calculated unsteady pressure distribution is used to determine the unsteady blade loading. The calculations at the design point and at two off-design points are carried out with a multiple frame of reference and a sliding mesh technique is applied to consider the impeller/volute interaction. The results obtained show that the flow in the impeller and volute casing is periodically unsteady and confirm the circumferential distortion of the pressure distribution at the impeller outlet and in the volute casing. Due to the interaction between impeller blades and the tongue of the volute casing the flow is characterized by pressure fluctuations, which are strong at the impeller outlet and in the vicinity of the tongue. These pressure fluctuations are died away in the casing as the advancement angle increases. These reduced pressure fluctuations are spread to the discharge nozzle; the pressure fluctuations are also reflected to the impeller inlet and they affect the mass flow rate through the blade passages.


2002 ◽  
Vol 471 ◽  
pp. 51-70 ◽  
Author(s):  
M. PÉREZ-SABORID ◽  
M. A. HERRADA ◽  
A. GÓMEZ-BAREA ◽  
A. BARRERO

We present a numerical study of the downstream evolution (mechanical and thermal) of vortex-jet cores whose velocity and temperature fields far from the axis match a family of inviscid and non-conducting vortices. The far-velocity field is rotational, except for a particular case which corresponds to the well-known Long's vortex. The evolution of the vortex core depends on both the conditions at a certain upstream station, characterized by the dimensionless value of the velocity at the axis, and a dimensionless swirling parameter L defined as the ratio of the values of the azimuthal and axial velocities outside the vortex core. This numerical study, based on the quasi-cylindrical approximation (QC) of the Navier–Stokes equations, determines the conditions under which the vortex evolution proceeds smoothly, eventually reaching an asymptotic self-similar behaviour as described in the literature (Fernández-Feria, Fernández de la Mora & Barrero 1995; Herrada, Pérez-Saborid & Barrero 1999), or breaks in a non-slender solution (vortex breakdown). In particular, the critical value L = Lb(a) beyond which vortex breakdown occurs downstream is a function of a dimensionless parameter a characterizing the axial momentum of the vortex jet at an initial upstream station. It is found numerically that for very large values of a this vortex breakdown criterion tends to an asymptote which is precisely the value L = L* predicted by the self-similar analysis, and beyond which a self-similar structure of the vortex core does not exist. In addition, the computation of the total temperature field provides useful information on the physical mechanisms responsible for the thermal separation phenomenon observed in Ranque–Hilsch tubes and other swirling jet devices. In particular, the mechanical work of viscous forces which gives rise to an intense loss of kinetic energy during the initial stages of the evolution has been identified as the physical mechanism responsible for thermal separation.


Author(s):  
Elena Schneider ◽  
Amsini Sadiki ◽  
Alexander Maltsev ◽  
Johannes Janicka

Swirl flows play an important role in modern combustion systems such as gas turbines, aero propulsion systems etc. Next to desirable effects such as enhanced mixing such flows often exhibit aerodynamical instabilities called precessing vortex core. The configuration under study here represents a model Gas Turbine(GT) combustion chamber and features the main properties of real gas turbine combustors: a confined swirled flow with multiple recirculation zones and reattachment points, resulting in reacting case in a partially premixed methane/air aerodynamically stabilised flame. This flame exibits also precessing vortex core (PVC). The present study especially concentrates on an evaluation of the performance of different URANS-based model-combinations in predicting this confined swirling reacting flow exhibiting such aerodynamic instabilities. For this purpose an extended Bray-Moss-Libby model and a G-equation based approach, both coupled to the mixture fraction transport equation to account for partially premixed effects, are applied. Their prediction potential in capturing partially premixed combustion properties is appraised by comparison with LDV, Raman and PLIF measurements. It turns out that the influence of the combustion model on simulation results of the flame front stabilisation or mean flow field is not obvious. Nevertheless it could be mentioned that the computation time with G-equation was approximately three times longer than with BML model due to the reinitialization needed in steady case calculations and 2 times longer in case of unsteady calculations.


Sign in / Sign up

Export Citation Format

Share Document