基于WMLES方法的复式断面明渠三维紊流数值模拟

针对基于雷诺时均Navier-Stokes方程的数值模型在精确模拟复式断面明渠三维水流结构方面存在困难和传统大涡模拟(LES)方法计算成本相对较高的问题,采用壁面建模的大涡模拟(WMLES)方法建立三维数值模型,对水深比(滩地水深与主槽水深之比)为0.5的复式断面明渠水流结构进行模拟计算。数值模拟计算结果与已有试验测量结果和LES模拟结果的对比表明,WMLES方法能够准确地模拟时均速度、床面切应力、紊动强度、紊动能、雷诺应力分布、二次漩涡结构,该方法在复式断面明渠水流三维模拟中是可靠的;在保证计算精度的前提下,与传统LES方法相比,采用WMLES方法能够显著降低计算成本。     

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.     

用SST-DES和SST-URANS方法数值模拟亚临界雷诺数下三维圆柱绕流问题

基于剪切应力运输(shear stress transport,SST)的分离涡模拟(detached-eddy simulation,DES)在近壁面的流动区域采用SST模型,而在其他区域采用亚格子模型求解流场。本研究中分别应用SST-DES和SST-URANS方法对亚临界雷诺数Re=3900下的三维圆柱绕流进行了数值模拟,比较分析了两种数值方法下圆柱表面的时均压力分布、圆柱后方的漩涡脱落特征、圆柱下游流场的时均速度剖面分布等,并将数值模拟结果同前人物理试验比较,结果表明SST-DES方法在亚临界雷诺数下的三维圆柱绕流数值模拟问题上比SST-RANS方法具有更好的效果。     

An hybrid RANS/LES model for simulation of complex turbulent flow

A non-linear eddy viscosity model(NLEVM) and a scalable hybrid Reynolds averaged Navier-Stokes/ large eddy simulation(RANS/LES) strategy are developed to improve the capability of the eddy viscosity model(EVM) to simulate complex flows featuring separations and unsteady motions. To study the performance of the NLEVM, numerical simulations around S809 airfoil are carried out and the results show that the NLEVM performs much better when a large separation occurs. Calculated results of the flow around a triangular cylinder show that the NLEVM can improve the precision of the flow fields to some extents, but the error is still considerable, and the small turbulence structures can not be clearly captured as the EVM. Whereas the scalable hybrid RANS/LES model based on the NLEVM is fairly effective on resolving the turbulent structures and can give more satisfactory predictions of the flow fields.     

PIV and PTV measurements in hydro-sciences with focus on turbulent open-channel flows

PIV is one of the most popular measurement techniques in hydraulic engineering as well as in fluid sciences. It has been applied to study various turbulent phenomena in laboratory experiments related to natural rivers, e.g., bursting phenomena near the bed, mixing layers observed at confluences, wake turbulence around dikes and piers, and so on. In these studies, PIV plays important roles in revealing the space-time structure of velocity fluctuations and coherent vortices. This review article focuses particularly on the applications of PIV to turbulent open-channel flows, which have been conducted for the past decade in Hydraulics Laboratory of Kyoto University. In Section 2, we introduce our experimental setup and PIV/PTV algorithm. In Section 3, we apply the PIV measurements to reveal turbulence characteristics and coherent structures in open-channel flows as well as in vegetated canopy flows. For complex flow situations, various applications of PIV to compound open-channel flows and wind-induced water waves are considered to reveal coherent vortices. In Section 4, we discuss some advanced PIV measurements in open-channel flows. The free-surface-elevation fluctuations and velocity components were measured simultaneously with two sets of cameras to examine phase-averaged parameters of turbulence. A multi-layer scanning PIV was developed to reveal 3D turbulence structure in compound open-channel flows. Our discriminator PIV/PTV was applied successfully to sediment-laden open-channel flows and revealed the fluid/particle interaction and the relationship between coherent structures and sediment concentration. Finally, we conducted simultaneous measurements of velocity and dye concentration with a combination of PIV and LIF in vegetated open-channel flow, which enables us to examine turbulent scalar flux of a passive contaminant.     

A hybrid RANS-LES model for combining flows in open-channel T-junctions

In the present study, a hybrid RANS-LES (Reynolds Averaged Navier-Stokes equation — Large Eddy Simulation) model is developed for the simulation of open-channel T-junction flows. The hybrid approach can save computational effort comparing to the LES approach. The comparison between the computational results and the detailed experimental data shows that this relatively new modeling approach is more accurate than the RANS approach in the prediction of the combining flows in T-junctions.

     

子格（SGS）模型在内流湍流中的大涡模拟

本文将子格模型应用于大涡模拟，对直方管，弯方管和弯曲槽道内充分展湍流，进行了数值模拟，并将所得的结果与实验数据作比较，计算表明，应用于子格模型的大涡模拟可以用来模拟比较复杂的内流，研究湍流运动的结构，检验湍流模型。     

On the numerical simulations of vortical cavitating flows around various hydrofoils

This paper reviews the numerical models of various cavitating flows around hydrofoils. Numerical models relating to cavitation flows, including mass transfer models and turbulence models, are summarized at first. Then numerical results and analysis of flow characteristics for the cavitating flows around twisted hydrofoils, truncated hydrofoil and tip leakage are discussed respectively. For mean flow fields, Reynolds averaged Navier-Stokes(RANS) simulation associated with a kind of nonlinear turbulence model is found to be an economic and robust numerical approach for different kinds of cavitating flows including cloud cavitation, tip cavitation and tip leakage cavitation. To predict the fluctuations of pressure and velocity, large eddy simulation(LES) and detached eddy simulation(DES) are two effective approaches. Finally, a few open questions are proposed for future research.     

A Partially-Averaged Navier-Stokes model for hill and curved duct flow

Turbulent flows past hill and curved ducts exist in many engineering applications. Simulations of the turbulent flow are carried out based on a newly developed technique, the Partially-Averaged Navier-Stokes (PANS) model, including separation, recirculation, reattachment, turbulent vortex mechanism. The focus is on how to accurately predict typical separating, reattaching and secondary motion at a reasonable computational expense. The effect of the parameter, the unresolved-to-total ratio of kinetic energy (f_k), is examined with a given unresolved-to-total ratio of dissipation (f_?) for the hill flow with a much coarser grid system than required by the LES. An optimal value of f_k can be obtained to predict the separation and reattachment locations and for more accurate simulation of the resolved turbulence. In addition, the turbulent secondary motions are captured by a smaller f_k as compared with the RANS method with the same grid.     

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.     

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 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…     

On the turbulence characteristics in developed and developing rough narrow open-channel flow

An analytical power-law for velocity distribution is developed and validated using experimental results for rough, narrow open-channel flows. In addition, an empirical power-law for the mean velocity distribution is also proposed. To this end, experiments were conducted to investigate the turbulence characteristics in developing and fully developed narrow open-channel flows over a fixed continuous rough bed. Instantaneous 3-D velocities were acquired using a Nortek VectrinoPlus down looking acoustic Doppler velocimeter at streamwise intervals of 0.5 m along the centerline of the flume. Using the experimental data, shape factors and frictional characteristics in narrow open-channel flow were analysed. In addition, this paper investigated changes in the turbulence anisotropy and the decay of turbulence along the mid-depth of the hydraulically rough flow. The new analytical power-law equation developed in this paper accurately describes the fully developed velocity profile in narrow open-channel flow with the dip phenomenon. Reynolds number and aspect ratio significantly influenced the turbulence characteristics in the upper half of the flow field. The decreasing turbulence intensities at mid-depth along the flow developing region demonstrate the attainment of isotropic turbulence as the flow development occurs in the narrow open-channel.     

Numerical investigation of drag force on flow through circular array of cylinders

A 2-D model for flow through a circular patch with an array of vertical circular cylinders in a channel is established using the Navier-Stokes equations with a hybrid RANS/LES turbulence model–the Scale Adaptive Simulation (SAS) model. The applicability of the model is first validated by test cases where experimental data are available for comparison with the computed results. It is verified that the present model can predict well the average velocity and turbulence structure. The drag force and drag coefficient are then calculated using the present model for a number of cases with different solid volume fractions, cylinder Reynolds numbers and patch diameters. It is shown that the drag coefficient increases with increasing solid volume fraction, but decreases with increasing Reynolds number. However, the drag coefficient is independent of the diameter of circular batch when the solid volume fraction and Reynolds number are kept constant.     

模拟敞开式UV反应槽的湍流模型比较研究

应用商用CFD软件F luent模拟敞开式UV反应槽内的流速分布及停留时间分布。分别采用4种湍流模型来评估湍动效应对流场分布的影响,其中包括标准k-ε模型、低雷诺数k-ε模型、雷诺应力模型(RSM)及分离涡(DES)模型。将模拟结果与实验数据进行对比发现,DES模型对于回流区的模拟好于其他模型;而在回流区域外,低雷诺数k-ε模型则更为合适。     

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

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

宽浅式尾水渠体型设计

采用一维和三维k-ε紊流数值模拟方法相结合,对宽深比达到500∶1的尾水渠水流流场进行了模拟。利用一维计算速度快的优点对30个尾水渠体型共60组工况进行了水面线计算,对壁面糙率进行了敏感性分析,提出了满足设计要求、工程量最省的体型。利用三维k-ε数值计算能得到详细流场的优点,分析了不同机组开启时尾水渠内流速分布和水面特性。研究结果表明：推荐方案在小流量时能满足尾水管出口的埋深,在大流量时保证设计有效发电水头且工程量最省。     

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…     

《南水北调与水利科技》优先数字出版声明

The turbulent flows through the channels with abrupt cross-sectional changes are common and importantphysical process in nature.For a better prediction of the mean flow and turbulent characteristics for this problem,atwo-dimensional depth-averaged numerical model is developed.The model is robust and accurate in reproducing therecirculation flow behind a groyne and turbulent flows in channels with abrupt cross-sectional changes,when com-pared to the available experimental data of mean velocities and turbulence kinetic energy.Our results reveal that theabrupt cross-sectional change of a channel can affect the flow pattern significantly and introduces the complex turbu-lence characteristics.In particular,when the channel has an abrupt expansion,the mean flow pattern is mainly in lon-gitudinal direction with rather small transverse component.Meanwhile,a recirculating region forms behind the expan-sion position and the turbulence has very strong intensity within this region.For the flow in the channel with an ab-rupt contraction,the longitudinal component of the flow is decreased by the obstruction on one side and accelerated onthe other side,whereas the transverse velocity is small.The turbulence is extraordinarily strong in the regions adja-cent to the contraction wall in the narrow channel.In both cases of abrupt cross-sectional changes,the TKE is genera-ted dominantly by the shear of the longitudinal velocities.     

FLOW STRUCTURE OF PARTLY VEGETATED OPEN-CHANNEL FLOWS WITH EELGRASS

Aquatic vegetation can influence the transport of sediment and contaminants by changing the mean velocity and turbulent flow structure in channels. It is important to understand the hydraulics of the flows over vegetation in order to manage fluvial processes. Experiments in an open-channel flume with natural vegetation were carried out to study the influence of vegetation on the flows. In a half channel with two different densities of vegetation, the flow velocity, Reynolds stresses, and turbulence intensities were measured using an Acoustic Doppler Velocimeter (ADV). We obtained velocity profiles in the lateral direction, Reynolds stresses in the vertical direction, and the flow transition between the vegetated and non-vegetated zones in different flow regimes. The results show that the streamwise velocity in the vegetated zone with higher density is almost entirely blocked. Reynolds stress distribution distinguishes with two different regions: inside and above the vegetation canopies. The turbulence intensities increase with increasing Reynolds number. The coherent vortices dominate the vertical transport of momentum and are advected clockwise between the vegetated zone and non-vegetated zone by secondary currents (a relatively minor flow superimposed on the primary flow, with significantly different speed and direction), generated by the anisotropy of the turbulence.