Modeling the Unsteady Cavitating Flow in a Cross-Flow Water Turbine

2010 ◽  
Vol 132 (7) ◽  
Author(s):  
E. Sansone ◽  
C. Pellone ◽  
T. Maitre

The noncavitating and cavitating flows over a cross-flow water turbine are simulated by using an unsteady Navier–Stokes formulation. For the cavitating flow case, a homogeneous mixture with a varying density is considered and one additional transport equation is explicitly solved in time for the liquid volume fraction. The instantaneous rate of vapor production and absorption appearing as a source term is governed by a hydrodynamic model based on a simplified bubble dynamic equation. The spatial discretization is achieved by a 2D multiblock technique consisting of fixed and rotating blocks, which were especially adapted for Darrieus geometry. Several test cases corresponding to experiments performed on fixed and rotating blades are selected to compare the numerical results with experimental data. Finally, a calculation of a monobladed cavitating cross-flow turbine is presented. The effect of cavitation on the dynamic stall phenomenon and on the turbine performance is analyzed. In particular, it is shown that cavitation earlier reveals the stall phenomenon on the blades and magnifies the size of the shedding vortex structures in the turbine.

2005 ◽  
Vol 40 ◽  
pp. 179-184 ◽  
Author(s):  
C. Bock ◽  
H. Eicken

AbstractThe microstructural evolution of brine inclusions in granular and columnar sea ice has been investigated through magnetic resonance imaging (MRI) for temperatures between –28 and –3˚C. Thin-section and salinity measurements were completed on core samples obtained from winter sea ice near Barrow, Alaska, USA. Subsamples of granular (2–5cm depth in core) and columnar sea ice (20–23 cm depth) were investigated with morphological spin-echo and diffusion-weighted imaging in a Bruker 4.7T MRI system operating at field gradients of 200 mTm–1 at temperatures of approximately –28, –15, –6 and –3˚C. Average linear pore dimensions range from 0.2 to 1 mm and increase with bulk liquid volume fraction as temperatures rise from –15 to –3˚C. Granular ice pores are significantly larger than columnar ice pores and exhibit a higher degree of connectivity. No evidence is found of strongly non-linear increases in pore connectivity based on the MRI data. This might be explained by shortcomings in resolution, sensitivity and lack of truly three-dimensional data, differences between laboratory and field conditions or the absence of a percolation transition. Pore connectivity increases between –6 and –3˚C. Pore-number densities average at 1.4±1.2mm–2. The pore-number density distribution as a function of cross-sectional area conforms with power-law and lognormal distributions previously identified, although significant variations occur as a function of ice type and temperature. At low temperatures (< –26˚C), pore sizes were estimated from 1H self-diffusivity measurements, with self-diffusivity lower by up to an order of magnitude than in the free liquid. Analysis of diffusional length scales suggests characteristic pore dimensions of <1 μm at < –26˚C.


Author(s):  
Luis San Andrés ◽  
Xueliang Lu

Wet gas compression systems and multiphase pumps are enabling technologies for the deep sea oil and gas industry. This extreme environment determines both machine types have to handle mixtures with a gas in liquid volume fraction (GVF) varying over a wide range (0 to 1). The gas (or liquid) content affects the system pumping (or compression) efficiency and reliability, and places a penalty in leakage and rotordynamic performance in secondary flow components, namely seals. In 2015, tests were conducted with a short length smooth surface annular seal (L/D = 0.36, radial clearance = 0.127 mm) operating with an oil in air mixture whose liquid volume fraction (LVF) varied to 4%. The test results with a stationary journal show the dramatic effect of a few droplets of liquid on the production of large damping coefficients. This paper presents further measurements and predictions of leakage, drag power, and rotordynamic force coefficients conducted with the same test seal and a rotating journal. The seal is supplied with a mixture (air in ISO VG 10 oil), varying from a pure liquid to an inlet GVF = 0.9 (mostly gas), a typical range in multiphase pumps. For operation with a supply pressure (Ps) up to 3.5 bar (a), discharge pressure (Pa) = 1 bar (a), and various shaft speed (Ω) to 3.5 krpm (ΩR = 23.3 m/s), the flow is laminar with either a pure oil or a mixture. As the inlet GVF increases to 0.9 the mass flow rate and drag power decrease monotonically by 25% and 85% when compared to the pure liquid case, respectively. For operation with Ps = 2.5 bar (a) and Ω to 3.5 krpm, dynamic load tests with frequency 0 < ω < 110 Hz are conducted to procure rotordynamic force coefficients. A direct stiffness (K), an added mass (M) and a viscous damping coefficient (C) represent well the seal lubricated with a pure oil. For tests with a mixture (GVFmax = 0.9), the seal dynamic complex stiffness Re(H) increases with whirl frequency (ω); that is, Re(H) differs from (K-ω2M). Both the seal cross coupled stiffnesses (KXY and −KYX) and direct damping coefficients (CXX and CYY) decrease by approximately 75% as the inlet GVF increases to 0.9. The finding reveals that the frequency at which the effective damping coefficient (CXXeff = CXX-KXY/ω) changes from negative to positive (i.e., a crossover frequency) drops from 50% of the rotor speed (ω = 1/2 Ω) for a seal with pure oil to a lesser magnitude for operation with a mixture. Predictions for leakage and drag power based on a homogeneous bulk flow model match well the test data for operation with inlet GVF up to 0.9. Predicted force coefficients correlate well with the test data for mixtures with GVF up to 0.6. For a mixture with a larger GVF, the model under predicts the direct damping coefficients by as much as 40%. The tests also reveal the appearance of a self-excited seal motion with a low frequency; its amplitude and broad band frequency (centered at around ∼12 Hz) persist and increase as the gas content in the mixture increase. The test results show that an accurate quantification of wet seals dynamic force response is necessary for the design of robust subsea flow assurance systems.


Author(s):  
Tan Dung Tran ◽  
Bernd Nennemann ◽  
Thi Cong Vu ◽  
François Guibault

The objective of this paper is to evaluate the applicability of different cavitation models and determine appropriate numerical parameters for cavitating flows around a hydrofoil. The simulations are performed for a NACA 66 foil at 6 degrees angle of attack, Reynolds number of 750 000 and for a cavitation number of 1.49 corresponding to the partial sheet cavitating regime. The incompressible, multiphase Reynolds-averaged Navier-Stokes (RANS) equations are solved by the CFD solver CFX with Kubota and Merkle cavitation models. As part of the work, the Merkle model is implemented into CFX by User Fortran code because this model has shown good cavitation prediction capability according to the literature. The effects of the k-ε and SST turbulence models on the cavitating flow dynamics are compared. Also, an investigation on structured and hybrid meshes with different mesh sizes and concentrations is carried out in order to better understand the mesh influence for this cavitation simulation. The local compressibility effect is considered by correcting the turbulent eddy viscosity inside the mixture vapor/liquid zones. The numerical results are validated by experiments conducted in a cavitation tunnel at the French Naval Academy.


Author(s):  
Feng Hong ◽  
Jianping Yuan ◽  
Banglun Zhou ◽  
Zhong Li

Compared to non-cavitating flow, cavitating flow is much complex owing to the numerical difficulties caused by cavity generation and collapse. In the present work, cavitating flow around a two-dimensional Clark-Y hydrofoil is studied numerically with particular emphasis on understanding the cavitation structures and the shedding dynamics. A cavitation model, coupled with the mixture multi-phase approach, and the modified shear stress transport k-ω turbulence model has been developed and implemented in this study to calculate the pressure, velocity, and vapor volume fraction of the hydrofoil. The cavitation model has been implemented in ANSYS FLUENT platform. The hydrofoil has a fixed angle of attack of α = 8° with a Reynolds number of Re = 7.5 × 105. Simulations have been carried out for various cavitation numbers ranging from non-cavitating flows to the cloud cavitation regime. In particular, we compared the lift and drag coefficients, the cavitation dynamics, and the time-averaged velocity with available experimental data. The comparisons between the numerical and experimental results show that the present numerical method is capable to predict the formation, breakup, shedding, and collapse of the sheet/cloud cavity. The periodical formation, shedding, and collapse of sheet/cloud cavity lead to substantial increase in turbulent velocity fluctuations in the cavitation regimes around the hydrofoil and in the wake flow.


2019 ◽  
Vol 245 ◽  
pp. 106859
Author(s):  
Joaquín López ◽  
Julio Hernández ◽  
Pablo Gómez ◽  
Claudio Zanzi ◽  
Rosendo Zamora

2020 ◽  
Vol 61 (2) ◽  
Author(s):  
Lukas Weiss ◽  
Michael Wensing ◽  
Joonsik Hwang ◽  
Lyle M. Pickett ◽  
Scott A. Skeen

2019 ◽  
Vol 21 (2) ◽  
pp. 248-262 ◽  
Author(s):  
David Markt ◽  
Ashish Pathak ◽  
Mehdi Raessi ◽  
Seong-Young Lee ◽  
Roberto Torelli

This article uniquely characterizes the secondary droplets formed during the impingement of a train of ethanol drops, using three-dimensional direct numerical simulations performed under conditions studied experimentally by Yarin and Weiss. Our numerical results have been previously validated against experimental data demonstrating the ability to accurately capture the splashing dynamics. In this work, the predictive ability of the model is leveraged to gain further insight into secondary droplet formation. We present a robust post-processing algorithm, which scrutinizes the liquid volume fraction field in the volume-of-fluid method and quantifies the number, volume and velocity of secondary droplets. The high-resolution computational simulations enable secondary droplet characterization within close proximity of the impingement point at small length and time scales, which is extremely challenging to achieve experimentally. By studying the temporal evolution of secondary droplet formation, direct connections are made between liquid structures seen in the simulation and the instantaneous distribution of secondary droplets, leading to detailed insight into the instability-driven breakup process of lamellae. Time-averaged secondary droplet characteristics are also studied to describe the global distribution of secondary droplets. Such analysis is vital to understanding fuel drop impingement in direct injection engines, facilitating the development of highly accurate spray–wall interaction models for use in Lagrangian solvers.


Sign in / Sign up

Export Citation Format

Share Document