LARGE-EDDY SIMULATION OF TURBULENT FLOW CONSIDERING INFLOW WAKES IN A FRANCIS TURBINE BLADE PASSAGE

Turbulent flow in a 3-D blade passage of a Francis hydro turbine was simulated with the Large Eddy Simulation (LES) to investigate the spatial and temporal distributions of the turbulence when strongly distorted wakes in the inflow sweep over the passage. In a suitable consideration of the energy exchanging mechanism between the large and small scales in the complicated passage with a strong 3-D curvature,one-coefficient dynamic Sub-Grid-Scale (SGS) stress model was used in this article. The simulations show that the strong wakes in the inflow lead to a flow separation at the leading zone of the passage,and to form a primary vortex in the span-wise direction. The primary span-wise vortex evolves and splits into smaller vortex pairs due to the constraint of no-slip wall condition,which triggers losing stability of the flow in the passage. The computed pressures on the pressure and suction sides agree with the measured data for a working test turbine model.     

3-D hybrid LES-RANS model for simulation of open-channel T-diversion flows

The study of flow diversions in open channels plays an important practical role in the design and management of open-channel networks for irrigation or drainage.To accurately predict the mean flow and turbulence characteristics of open-channel dividing flows,a hybrid LES-RANS model,which combines the large eddy simulation (LES) model with the Reynolds-averaged Navier-Stokes (RANS) model,is proposed in the present study.The unsteady RANS model was used to simulate the upstream and downstream regions of a main channel,as well as the downstream region of a branch channel.The LES model was used to simulate the channel diversion region,where turbulent flow characteristics ate complicated.Isotropic velocity fluctuations were added at the inflow interface of the LES region to trigger the generation of resolved turbulence.A method based on the virtual body force is proposed to impose Reynolds-averaged velocity fields near the outlet of the LES region in order to take downstream flow effects computed by the RANS model into account and dissipate the excessive turbulent fluctuations.This hybrid approach saves computational effort and makes it easier to properly specify inlet and outlet boundary conditions.Comparison between computational results and experimental data indicates that this relatively new modeling approach can accurately predict open-channel T-diversion flows.     

NUMERICAL SIMULATION OF 3-D TURBULENT FLOWS OF PLUNGE POOL AND ENERGY DISSIPATION ANALYSIS

A 3-D numerical model for simulating the complicated turbulent flows was developed and the code was made. A numerical calculation of the plunge pool of Laxiwa project in China was carried on. Those 3-D distributions of velocity,pressure on bottom wall,turbulence kinetic energy and turbulence energy dissipation rate are revealed in detail and the detailed flow patterns of plunge pool were shown. By studying on the characteristics of turbulent diffusion and energy dissipation,the calculated results show that the major turbulence and energy dissipation are taken place near the axis of water jet. The calculated results also indicated that the calculated maximum impact pressures on the bottom wall of the plunge pool have a good agreement with those obtained by physical hydraulic model test.     

NUMERICAL SIMULATION OF THE ENERGY DISSIPATION CHARACTERISTICS IN STILLING BASIN OF MULTI-HORIZONTAL SUBMERGED JETS

3-D numerical simulation was carried out for the water flow in a stilling basin with multi-horizontal submerged jets by using two different turbulence models, namely, the VOF RNG k - ? and Mixture RNG k - ? turbulence models. The calculated water depth, velocity profile and pressure distribution are in good agreement with the data obtained in experiments. It indicates that the numerical simulation can effectively be used to study the water flow movement and the energy dissipation mechanism. The numerical simulation results show that the turbulent kinetic energy distribution obtained by using the Mixture turbulence model covers a region about 18% larger than that calculated by using the VOF turbulence model, and is in better agreement with the actual situation. Furthermore, the Mixture turbulence model is better than the VOF turbulence model in calculating the air entrainment.     

常曲率U型弯道水流紊动特性试验研究

采用三维声学多普勒流速仪,对水流三维流速进行了采集,对弯道水流的脉动进行了统计分析,研究了典型断面的紊动强度、紊动切应力、紊动能的分布。试验结果能够很好的揭示弯道水流的三维特性及分布规律。研究成果为进一步研究水流的紊动结构及泥沙的运动机理提供了参考资料。     

LARGE-EDDY AND DETACHED-EDDY SIMULATIONS OF THE SEPARATED FLOW AROUND A CIRCULAR CYLINDER

The separated turbulent flow around a circular cylinder is investigated using Large-Eddy Simulation (LES), Detached-Eddy Simulation (DES, or hybrid RANS/LES methods), and Unsteady Reynolds-Averaged Navier-Stokes (URANS). The purpose of this study is to examine some typical simulation approaches for the prediction of complex separated turbulent flow and to clarify the capability of applying these approaches to a typical case of the separated turbulent flow around a circular cylinder. Several turbulence models, i.e. dynamic Sub-grid Scale (SGS) model in LES, the DES-based Spalart-Allmaras (S-A) and k ? ω Shear-Stress-Transport (SST) models in DES, and the S-A and SST models in URANS, are used in the calculations. Some typical results, e.g., the mean pressure and drag coefficients, velocity profiles, Strouhal number, and Reynolds stresses, are obtained and compared with previous computational and experimental data. Based on our extensive calculations, we assess the capability and performance of these simulation approaches coupled with the relevant turbulence models to predict the separated turbulent flow.     

2-D eddy resolving simulations of flow past a circular array of cylindrical plant stems

In the present study, 2-D large eddy simulations(LES) are conducted for flow past a porous circular array with a solid volume fraction(SVF) of 8.8%, 15.4% and 21.5%. Such simulations are relevant to understanding flow in natural streams and channels containing patches of emerged vegetation. In the simulations discussed in the paper, the porous cylinder of diameter D contains a variable number of identical solid circular cylinders(rigid plant stems) of diameter d= 0.048 D. Most of the simulations are conducted at a Reynolds number of 2 100 based on the diameter D and the velocity of the steady uniform incoming flow. Though in all cases wake billows are shed in the regions where the separated shear layers(SSLs) forming on the sides of the porous cylinder interact, the effect of these wake billows on the mean drag is different. While in the high SVF case(21.5%), the total drag force oscillates quasi-regularly in time, similar to the canonical case of a large solid cylinder, in the cases with a lower SVF the shedding of the wake billows takes place sufficiently far from the cylinder such that the unsteady component of the total drag force is negligible. The mean amplitude of the oscillations of the drag force on the individual cylinders is the largest in a streamwise band centered around the center of the porous cylinder, where the wake to wake interactions are the strongest. In all cases the maximum drag force on the individual cylinders is the largest for the cylinders directly exposed to the flow, but this force is always smaller than the one induced on a small isolated cylinder and the average magnitude of the force on the cylinders directly exposed to the flow decreases monotonically with the increase in the SVF. Predictions of the global drag coefficients, Strouhal numbers associated with the wake vortex shedding and individual forces on the cylinders in the array from the present LES are in very good agreement with those of 2-D direct numerical simulations conducted on finer meshes, which suggests LES is a better option to numerically investigate flow in channels containing canopy patches, given that LES is computationally much less expensive than DNS at high Reynolds number. To prove this point, the paper also discusses results of 2-D LES conducted at a much higher Reynolds number, where the near-wake flow is strongly turbulent. For the higher Reynolds number cases, where the influence of the turbulence model is important, the effect of the sub-grid scale model and the predictive capabilities of the unsteady Reynolds averaged Navier-Stokes(RANS) approach to predict flow past porous cylinders are discussed.     

NUMERICAL SIMULATION OF 3-D TURBULENT FLOWS IN A CUT-OFF VALVE AND ANALYSIS OF FLOW CHARACTERISTICS IN PIPE

1 .　INTRODUCTIONThevalvehasbeenplayinganimportantroleinthe pipeengineeringandisalsooneofmaincomponentsgeneratingtheheadloss .Becauseofthecomplexgeometry ,thedesignwasusuallyeval uatedempiricallyinthepast,whichcausesmuchenergyloss .Toreducetheenergyloss ,itisessen tialforadesignertoknowthedetailsofflowfieldinthevalvepipe .Recently ,withtheremarkableprogressintheComputationalFluidDynamics(CFD)andcomputerhardware ,itisexpectedthatCFDwillbecomeapowerfultoolinthedesignofvalve .Intheactual…     

NUMERICAL PREDICTION OF SILT ABRASIVE EROSION IN HYDRAULIC TURBINE

1.INTRODUCTIONTherearemanyhydraulicturbinesoperatedinsiltladenriversintheworld.Thepresenceofsolidparticlesintheflowoftencauseabrasivedegradationofcertainpartsofthehydraulicturbinewhichareespeciallyexposedtotheflow.Theerosivewearduetohydroabrasiondependsprimarilyonthesurfacematerialandonthematerialoftheabrasiveparticlesinwater.Thefluiddynamicsoftheparticlemotionalsoplayaveryimportantroleinthistribologicalsystem.Therelativevelocityofparticlesandthefrequencyofimpactarekeyfactorsindetermining…     

LARGE EDDY SIMULATION OF FREE SURFACE TURBULENT CHANNEL FLOW WITH HEAT TRANSFER

1. INTRODUCTIONFreesurface turbulent flow exists widely in engineering. Due to difficulties in measurement and simulation, very little has been understood about the structures and transport mechanism of turbulence near the freesurface. Komori and Ueda…     

NUMERICAL MODELING OF COMPOUND CHANNEL FLOWS

A numerical model capable of predicting flow characteristics in a compound channel was established with the 3-D steady continuity and momentum equations along with the transport equations for turbulence kinetic energy and dissipation rate. Closure was achieved with the aid of algebraic relations for turbulent shear stresses. The above equations were discretized with implicit difference approach and solved with a step method along the flow direction. The computational results showing the lateral distribution of vertical average velocities and the latio of total flow in the compound channel agree well with the available experimental data.     

三维圆柱绕流数值模拟湍流方法的选择

研究湍流模拟方法对三维圆柱绕流数值模拟精度的影响,分别采用雷诺平均法（RANS）中的κ-ω模型、SST模型以及大涡模拟法（LES）对亚临界区内雷诺数Re=3 900时的三维圆柱绕流进行数值计算,分析了圆柱体表面的受力、圆柱后流场时均速度特性与瞬时涡量分布情况。结果表明,当流体流过圆柱表面时,圆柱表面出现的与流速方向相反的负压力差区域使流体从圆柱表面分离,引起了不稳定的周期性交替脱落的湍流涡泄,从而在圆柱表面产生周期性波动的升力,同时在圆柱后近流场区域形成回流区。研究还发现,LES法对圆柱体的受力以及流场时均速度特性的模拟效果要优于κ-ω模型与SST模型;相较于前人利用浸入边界法得到的模拟结果,LES法的模拟精度也有了较大提升;同样,通过对瞬时流场涡量等值线图的分析,并与已有的模拟结果进行对比,发现LES法不但可以从整体上表现出漩涡的周期性脱落,而且对流场中不同位置的、复杂的小尺度湍流涡泄也描述得非常细致,得到的自由分离剪切层长度与湍流涡泄的卷曲度更符合湍流涡泄的特征。所以,在亚临界区,LES法对湍流的模拟效果相对较好。     

水垫塘三维紊流数值模拟及消能分析

开发了能模拟复杂紊流的三维数学模型，编制了大型计算程序，针对拉西瓦水垫塘进行了数值模拟，揭示了水垫塘内速度，压力，紊动能及其耗散率的三维分布规律，获得了详细的水垫塘流态特性。研究了水流紊动扩主消能特性，反映出紊动和消能大量地发生在射流轴线附近，计算所得的底板最大冲击动压与实验值进行了比较，二者符合良好。     

CFD SIMULATION OF 3D FLOW IN LARGE-BORE AXIAL-FLOW PUMP WITH HALF-ELBOW SUCTION SUMP

1. INTRODUCTION Axial-flow pumps are usually designed for pumping large volumes of water against low pressure head. Suction sump is one of key elements in pump stations. For pumps with larger than 2 m, closed type suction sump plans are used to provide pr…     

Vortex shedding in the flow around two side-by-side circular cylinders of different diameters

In this paper,the 3-D turbulent flow around two side-by-side circular cylinders of different diameters,at sub-critical Reynolds number(Re=3 900),is numerically simulated by the large eddy simulation(LES).The spacing ratios(T/D)between the two cylinders are considered in four cases(T/D=1.2,1.5,1.8 and 2.7)to study the vortex shedding and turbulent properties in the flow field.The main results are focused on the drag and lift coefficients,the vortex shedding frequency,the coherent structure,and the scale properties.It is shown that when T/D is equal to 1.2,the vortex shedding of the main cylinder is strongly suppressed by the small cylinder,the drag and lift coefficients of the main cylinder are smaller than those in other three cases.While T/D is equal to 1.5,the vortex shedding of the main cylinder can be improved,the drag and lift coefficients of the main cylinder are larger than those in other three cases.The empirical mode decomposition(EMD)method is applied to decompose the velocity signals traced by the LES.It is shown that there is a linear relationship between the mean period and the mode in the semi-log coordinates.The vortex shedding period of the main cylinder is consistent with the period of the restructured coherent structures quantitatively.     

Numerical simulation of turbulent flow around a forced moving circular cylinder on cut cells

Fixed and forced moving circular cylinders in turbulent flows are studied by using the Large Eddy Simulation （LES） and two-equation based Detached Eddy Simulation （DES） turbulence models. The Cartesian cut cell approach is adopted to track the body surface across a stationary background grid covering the whole computational domain. A cell-centered finite volume method of second-order accuracy in both time and space is developed to solve the flow field in fluid cells, which is also modified accordingly in cut cells and merged cells. In order to compare different turbulence models, the current flow past a fixed circular cylinder at a moderate Reynolds number, Re = 3 900, is tested first. The model is also applied to the simulation of a forced oscillating circular cylinder in the turbulent flow, and the influences of different oscillation amplitudes, frequencies and free stream velocities are discussed. The numerical results indicate that the present numerical model based on the Cartesian cut cell approach is capable of solving the turbulent flow around a body undergoing motions, which is a foundation for the possible future study on wake induced oscillation and vortex induced vibration.     

环形波纹钢管涵洞式鱼道水力特性数值模拟

环形波纹钢管涵洞式鱼道内波纹的粗糙度可以使边界附近产生足够低的流速,鱼类可以由此游向上游。对3种流量0.05,0.07,0.09 m3/s,3种埋深0D,0.1D,0.2D,坡度0.4%工况下,环形波纹钢管涵洞式鱼道内水面线、流速场和紊流场分布分别进行了数值模拟研究。结果表明环形波纹钢管涵洞中,高流速区（无量纲流速大于0.9）位于涵洞过水断面中心区域,高紊流区（紊流强度大于0.2）位于水面中心附近。在涵洞底部边壁附近较大区域内,由于流速较小,紊流强度也较低（紊流强度小于0.1）,鱼类可以由此顺利完成上溯。相比无埋深式涵洞,嵌入式涵洞内平均流速较小,紊流强度变化较为平缓,更有利于鱼类的洄游。     

Turbulence characteristics in free-surface flow over two-dimensional dunes

The turbulent structure of open-channel flows over two-dimensional dunes is investigated numerically using large-eddy simulation (LES), in order to improve our understanding of the interaction between the dune-generated turbulence and the free surface dynamics. The filtered Navier–Stokes equations in the LES model have been discretised using the finite volume method, with a dynamic sub-grid model being employed for the unresolved scales of turbulence. The partial cell treatment has been implemented in a Cartesian grid form to deal with the dune topography. Both the volume of fluid method and rigid lid approach have been employed in the numerical framework to investigate the effects of the free surface treatment on the flow characteristics. The numerical model predicted mean flow velocities, turbulence intensities and Reynolds stresses have been compared with experimental measurements published in the literature, with a detailed analysis being undertaken to assess the accuracy of the model results and the effects of the free surface treatment on the velocity and turbulence predictions. The instantaneous flow structure has been investigated, with emphasis being focused on the free surface dynamics and coherent structures.     

NUMERICAL PREDICTION OF VORTEX FLOW IN HYDRAULIC TURBINE DRAFT TUBE FOR LES

The three-dimensional unsteady turbulent flow is studied numerically in the whole flow passage of hydraulic turbine, and vortex flow in the draft tube is predicted accurately in this paper. The numerical prediction is based on the Navier-Stokes equations and Large-Eddy Simulation （LES） model. The SIMPLE algorithm with the body fitted coordinate and tetrahedroid grid system is applied for the solution of the discretization governing equations.     

APPLICATION OF LBM-SGS MODEL TO FLOWS IN A PUMPING-STATION FOREBAY

A 2-D Lattice Boltzmann Method (LBM) coupled with a Sub-Grid Stress (SGS) model is proposed and validated by flows around a non-submerged spur dike in a channel. And then the LBM-SGS model is further applied to flows in a pumping-station forebay. Shallow water equations are numerically solved by the LBM and the turbulence can be taken into account and modeled efficiently by the Large Eddy Simulation (LES) model. The bounce-back scheme of the non-equilibrium part of the distribution function is used at the inlet boundary, the normal gradient of the distribution function is set as zero at the outlet boundary and the bounce-back scheme is applied to the solid wall to ensure non-slip boundary conditions. Firstly, the model successfully predicts the flow characteristics around a spur dike, such as circulating flow, velocity and water depth distributions. The results are verified by the experimental data and compared to the results obtained by conventional Smagoringsky Model (SM) of LES. Finally, the LBM-SGS model is used to further predict the flow characteristics in a forebay, such as secondary flow and water level. The comparisons show that the model scheme has the capacity to simulate complex flows in shallow water with reasonable accuracy and reliability.