明渠柱体上游马蹄涡的运动学特征研究

马蹄涡是水下柱体结构周围床面冲刷的主要动力。为掌握马蹄涡的运动学特征,引入了基于旋转强度的马蹄涡识别方法,提出了基于椭圆涡拟合的马蹄涡特征提取方法,其中,椭圆涡由Oseen涡和倾斜单向剪切流叠加而成。基于上述方法,通过开展明渠圆柱绕流PIV实验,对柱体上游对称面时均流场中马蹄涡的运动学特征进行了研究。结果表明,在柱体雷诺数大于104的紊流柱体绕流中,马蹄涡上游的流动分离点随柱体雷诺数的增大而缓慢向下游移动,但马蹄涡的位置保持不变,涡中心距柱体中轴约0.67D(D为柱体直径或宽度),距床面约0.06D;随着柱体雷诺数的增大,马蹄涡的半径不变而强度增加,使得马蹄涡作用于床面的剪切应力增大。研究结果建立了客观提取马蹄涡运动学特征的数学方法,得到了流动分离点、马蹄涡位置等运动学特征随水流条件的变化规律,揭示了明渠柱体绕流与其它流动中柱体绕流的差异,初步构建了马蹄涡的动力作用概念模型。     

Three-dimensional numerical simulation of flow around a circular cylinder under combined steady and oscillatory flow

Combined steady and oscillatory flow past a circular cylinder is investigated numerically by three-dimensional Direct Numerical Simulation (DNS). The incompressible Navier-Stokes equations are solved by finite element method (FEM). The aim of this study is to investigate influence of the existence of steady current on the flow regime and hydrodynamics forces. The computational results of pure oscillatory flow past a circular cylinder agree well with the experimental data. The flow ratios (ratio of current velocity to oscillatory velocity amplitude) are 0.0, 0.2 0.5 and 1.0 in the study. The influence of the flow ratio on the vortex shedding regime and hydrodynamic forces on the cylinder are investigated numerically.

     

Flexural-gravity waves due to unsteady motion of point source under a floating plate in fluid of finite depth

The paper deals with the theoretical investigation of the straight unsteady motion of a point source submerged into the liquid of finite depth below the floating ice plate. The formulae describing the plate deflection and obtained with the help of integral and asymptotic methods are numerically analyzed with respect to velocity and submergence depth of point source, basin depth, ice plate thickness and Young’s modulus.

     

LIQUID PHASE FLOW ESTIMATION IN GAS-LIQUID TWO-PHASE FLOW USING INVERSE ANALYSIS AND PARTICLE TRACKING VELOCIMETRY

An inverse analysis algorithm is proposed for estimating liquid phase flow field from measurement data of bubble motion. This kind of technology will be applied in future for various estimation of fluid flow in rivers, lakes, sea surface flow, and also microscopic channel flow as the problem-handling in civil, mechanical, electronic, and chemical engineering. The relationship between the dispersion motion and the carrier phase flow is governed and expressed by the translational motion equation of spherical dispersion. The equation consists of all the force components including inertia, added inertia, drag, lift, pressure gradient force and gravity force.Using this equation enables us to estimate the carrier phase flow structure using only the data of the dispersion motion.Whole field liquid flow structure is also estimated using spatial or temporal interpolation method. In order to verify this principle, the Taylor-Green vortex flow, and the Karman vortex shedding from a square cylinder have been chosen. The results show that the combination of the inverse analysis and Particle Tracking Velocimetry (PTV) with the spatio-temporal postprocessing algorithm could reconstruct well the carrier phase flow of the gas-liquid two-phase flow.     

浅尾水深水平底板上表面射流水力特性的数值模拟研究

采用VOF(volume of fluid)方法追踪自由液面,辅以Realizable(可实现化)湍流模型封闭两相流时均方程,对浅尾水深水平底板上的表面射流水力特性进行了数值模拟。微分方程的离散采用有限体积法,速度与压力耦合求解使用了压力隐式算子分裂PISO(Pressure-Implicit with Splitting of Operators)算法。研究并分析了入射速度、尾水深和进口宽度等对流场流线分布、速度场及横断面最大流速沿程衰减规律等水力特性的影响。研究发现:流场内存在一个大的漩涡,入射流速越大,漩涡长度越长;横断面上水平流速和纵向流速在漩涡区和明渠流区分布规律不同;横断面最大流速沿程衰减主要受进口流速大小的影响。     

三维薄翼入水水动力的计算

本文提供了一个计算三维薄翼高速入水水动力性能的非定常数值升力面方法。假设薄翼高速入水时,翼的背面为全空泡状态,当翼的随边进入液面后,此空泡域按一定的方向拖至自由液面。在翼的拱弧面上布置线涡和线源,在翼后空泡面上布置线源,所有的奇点系在自由液面上方存在映象以满足自由液面条件。根据翼迎流面的法向不可穿透条件和翼背面及翼后空泡面上的空泡动力学条件确定所有待定线涡和线源的强度,并由此可以计算各时刻翼上压力分布。利用本文提供的方法首先计算了二维薄翼入水的水动力性能,将结果与理论解进行比较,两者结果较为一致。然后,研究了三维薄翼的入水问题。     

A numerical study on dispersion of particles from the surface of a circular cylinder placed in a gas flow using discrete vortex method

The dispersion of particles emitted from the surface of a circular cylinder placed in a gas flow at the Reynolds number of 200 000 is numerically investigated using the discrete vortex method coupled with a Lagrangian approach for solid particle tracking. The wake vortex patterns, the temporal-spatial distributions and trajectories as well as the dispersion functions for particles with various Stokes numbers（St） ranging from 0.001 to 1.0 are obtained. The numerical results reveal that：（1） Solid particles on the cylinder surface are picked up and then transported away from the cylinder by the wake vortex flow.（2） Solid particles emitted from the cylinder surface always follow the vortices in the cylinder wake, and the response of particles to wake vortices is directly related to their Stokes numbers（particles with St= 0.001, 0.0038, 0.01 can distribute both in the vortex core and around the vortex periphery, whereas those with St= 0.1, 1.0 can not enter the vortex core and congregate mainly around the vortex periphery）.（3） The particles move in rolling state in the wake region, and the dispersion intensity of particles in the lateral direction decreases remarkably as the Stokes number of particles is increased from 0.001 to 1.0.     

Numerical investigation of a swirling flow under the optimal perturbation

The nonlinear evolution of 3-D instability of a viscous swirling flow, namely, the Oseen vortex, is addressed by direct numerical simulation with a Reynolds number equal to 5000. The global optimal perturbation is considered as the initial perturbation. In axisymmetric cases, three flow regimes are found: (1) the linear growth; (2) the decay of perturbation energy; (3) secondary energy growth. The linear energy growth, which is characterized by the amplification of radial perturbations, is arrested by the interaction between the vortex ring and the Oseen vortex core. The development of the vortex structure is also investigated for non-symmetric flows.

     

侧向进水口立轴旋涡数值模拟

为了研究侧向进水口立轴旋涡的运动特性,以RNG模型为基础,引入VOF方法捕捉自由液面的变化,对立轴旋涡的形态及其运动轨迹进行了三维数值模拟,并对旋涡水面线和三维速度分布进行了分析。结果表明,侧向进水口前的旋涡是非对称的,其涡轴是一条空间曲线,旋涡下部在进水口的抽吸拉伸作用下发生扭曲。旋涡运动可分为三部分,第一部分水面线呈漏斗状,与对称旋涡类似,第二部分旋涡运动逐渐从竖向转为横向,第三部分为横向运动。在旋涡的上半部,轴向速度随着水深的增加而增大,当旋涡逐渐由竖向转为横向运动时,轴向速度逐渐减小直至为零。最大切向速度发生在淹没深度的52%处。     

THE BUBMAC METHOD-A NUMERICAL METHOD BUBBLE DYNAMICS IN VORTEX CORE

This paper presents a numerical method called BUBMAC method,whichis based on the Marker And Cell(MAC) technique,to numerically simulate the problemof bubble dynamics in vortex core. With the incorporation of the azimuthal velocity intothe momentum equation, the model takes into account the complete interaction betweenthe bubble dynamics and the vortex flow field.The three momentum equations (N-Seqs.) are solved numerically by the finite-difference methed, and the motion of bubblesurface is described by tracing massless marker particles distributed only on the surface ofthe bubble. With some important medifications to the original MAC method, the numeri-cal accuracy of the method is conisderably increased. The validation of the proposed BUB-MAC method is checked by comparing the numerical results with some available analyticalsolutions in the cases of spherical bubble evolution and with numerical results in the casesof nonspherical collapse of bubble and bubble evolution in tip-vortex core. The compar-iso     

水平进水口自由表面旋涡数值模拟研究

旋涡运动是一种复杂的水流运动,目前尚难得出理论上的精确解。自由表面旋涡是水工建筑物进水口前的常见水力学现象,其发生常会对水工建筑物造成不利影响。旋涡数值模拟难度较大,尤其是吸气旋涡较难模拟出。建立了进水口旋涡三维数值模拟的模型,模拟出吸气旋涡现象。分析了计算结果与物理模型试验结果的差异,探讨了VOF法在水气两项流中的运用,为进一步进行吸气旋涡的三维数值模拟提供了一定经验。     

方型人工鱼礁周围水流运动的数值模拟研究

采用计算流体动力学(CFD)技术,模拟得到了5种不同来流速度(0.1、0.2、0.4、0.6和0.8 m/s)下方型人工鱼礁(边长为3 m)周围的水流场。铅垂平面上的计算结果表明:当水流贴近礁体迎流面时,水流抬升而形成上升流;不同来流速度工况下,上升流最大速度均约为来流速度的0.64倍,上升流平均速度均为来流速度的0.12倍,而上升流的最大高度为礁体高度的2.62～2.65倍;由于流动分离,在礁体顶面形成一小旋涡区,而在礁体背水面后端形成背涡区;5种来流速度下的背涡区长度均为礁体高度的3.0～3.5倍,而背涡区高度为礁体高度的1.1～1.2倍。水平面上的计算结果表明:紧贴礁体四周为水流减速区,在礁体两侧形成小旋涡区(缓流区),而在礁体背部形成一大旋涡区(背涡区);背涡区内的水流旋涡结构随来流速度大小而变,但背涡区的最大宽度基本不随来流速度的改变而改变,约为礁体宽度的2倍。     

环状气泡动力学数值模型

本文假设气泡运动过程中流场是无旋、不可压缩的理想流体,运用势流理论导出气泡边界面运动的控制方程,建立气泡动力学模型,用边界积分法求解得到气泡等边界的变形及位置.在环状气泡阶段,本文在前人研究成果的基础上分别建立了轴对称涡面模型、轴对称涡环模型以及三维涡环气泡模型.对比不同模型的优缺点及其适用范围,总结了环状气泡模型建立的详细步骤及其必须注意的关键技术,并开发了相应的程序.运用不同模型模拟重力场中气泡的运动,并进行讨论.最后计算了流场中滞后流的速度以及压力随时间的变化,得到一些有规律性的曲线,旨在为相关环状气泡动态特性研究提供参考.     

Laboratory Investigation of Stilling Basin Slope Effect on Bed Scour at Downstream of Stepped Spillway: Physical Modeling of Javeh RCC Dam

Stepped spillway and stilling basin are one of the most important energy dissipation structures. Eventhough, most of energy dissipated by these structures, but in skimming flow, the upstream flow motion is nonaerated and the residual energy capable to destroyed structures during floods. In this study, effect of stilling basin slope on bed scour, downstream of Javeh dam was investigating. Experiments performed in hydraulic structures laboratory of the University of Kerman with six different discharges (5, 7, 13, 17, 25 and 30 l/s.m) and five various stilling basin slope (0.02, 0.01, 0, ?0.01 and???0.02). The parameters such as maximum scour depth (d_s), flow velocity (in three point), water depth on upstream and downstream of stepped spillway and stilling basin, the distance of the maximum scour depth to sill (L_s) and the gheometery of scour hole measured. Result shown that when stilling basin slopes was 0.02, the average of maximum relative scour depth, 47% Increased and in ?0.02, 52.2% Decreased. In addition, the distance of maximum scour depth until stilling basin increased by increasing and decreased by decreasing the stilling basin slope.

     

Physical and numerical modeling of swirling flow in a scroll vortex intake

Scroll vortex intakes are hydraulic structures commonly used in water supply, drainage and sewerage systems. Water flows into the intake via an eccentric approach channel and a vortex chamber imparts swirl to the flow, leading to a stable air core vortex along the dropshaft. Although much effort has been devoted to investigate the scroll vortex flow, yet the understanding on its hydraulic characteristics is still far from complete due to a lack of detailed velocity field measurements. This paper presents the first comprehensive measurement of the three-dimensional (3D) flow field of a scroll vortex intake using non-intrusive Laser Doppler Anemometry (LDA). A validated 3D Computational Fluid Dynamics (CFD) model with the Volume-Of-Fluid (VOF) method is developed for interpreting the measurement. It is found that the vortex flow in the scroll chamber resembles a free vortex and the circulation is approximately equal to that at the chamber inlet, with a thin bottom boundary layer. The vortex flow at the bellmouth outlet possesses a circulation constant smaller than that in the chamber, and its vertical velocity component is approximately constant across the flow thickness.     

Implicit Finite-Volume Scheme to Solve Coupled Saint-Venant and Darcy–Forchheimer Equations for Modeling Flow Through Porous Structures

In hydrodynamic modeling of flow through porous structures, the solution domain might encounter discontinuities. These include, for example, porous structure-open channel interface, porous dam-break, and heterogeneous porous structures. The treatment of discontinuity is challenging within a numerical scheme as it can be a source of instabilities. This study proposes a finite-volume method to solve coupled Saint-Venant and Darcy–Forchheimer equations for simulating free-surface flow through porous structures. For capturing shocks arising at discontinuous regions, an upwind scheme is utilized to maintain the solution monotone. Fully implicit methods can allow the choice of longer time steps. Since the current problem involves two nonlinear systems, namely the open-channel and seepage flow equations, the Picard method is adopted to linearize the system of equations. Unlike typical implicit schemes of seepage flows, herein, both flow depth and velocity matrices appear within the iterative process, threatening the convergence criterion. To converge iteration, the continuity equation's flux term is treated using the dynamic wave equation under the relaxation method. The present model is applicable to simulate gradually and rapidly unsteady flow through homogeneous and heterogeneous porous media under laminar, transitional, and fully developed turbulent flow regimes within various closed and/or open boundary conditions.

     

排沙漏斗悬移质泥沙运动数值模拟

排沙漏斗是一种利用立轴旋涡实现水沙分离的二级排沙设施。为了探明排沙漏斗悬移质泥沙运动特征和排沙机理,采用欧拉-拉格朗日方法模拟了排沙漏斗水气沙三相流动,其中水相和气相采用雷诺应力模型(RSM)和VOF方法描述,悬移质泥沙用颗粒轨道模型描述。通过对计算的悬移质含沙量垂向分布与实验结果的对比,验证了数学模型的可靠性。计算结果表明:在漏斗室外侧区域含沙量分布整体上较均匀,在漏斗室内侧区域的近底部存在一个明显的高浓度带。作者认为,上述排沙漏斗悬移质泥沙含沙量分布特征主要是由排沙漏斗内的二次流引起的。     

HYDRAULIC CHARACTERISTICS OF VERTICAL VORTEX AT HYDRAULIC INTAKES

The trace of vertical vortex flow at hydraulic intakes is of the shape of spiral lines,which was observed in the presented experiments with the tracer technique. It represents the fluid particles flow spirally from the water surface to the underwater and rotate around the vortex-axis multi-cycle. This process is similar to the movement of screw. To describe the multi-circle spiral characteristics under the axisymmetric condition,the vertical vortex would change not only in the radial direction but also in the axial direction. The improved formulae for three velocity components for the vertical vortex flow were deduced by using the method of separation of variables in this article. In the improved formulae,the velocity components are the functions of the radial and axial coordinates,so the multi-circle spiral flow of vertical vortex could be simulated. The calculated and measured results for the vertical vortex flow were compared and the causes of errors were analyzed.     

NUMERICAL STUDY FOR THE NONLINEAR EVOLUTION OF VORTEX SHEET AND ITS APPLICATIONS: PART II. EVOLUTION OF WING SHEDDING VORTEX

1.　 INTRODUCTIONVortex sheet exists widely in a variety of nature and engineering problems,which isusually shed from airplane wing,ship propeller,swim ming fish,and so on.To study thebehaviour of vortex sheetevolution is much helpful to understand the physical mechanism ofthose flows.Therefore,many vortex m ethods[1 - 5 ] ,even involving viscous effects,weredeveloped to simulate the vortex sheet motion.As was discussed in Part I,asuperconvergence scheme based on Multhopp′s discretizati…     

EXPERIMENTAL AND NUMERICAL INVESTIGATION OF FREE SURFACE VORTEX

An experimental model was set up to investigate the formation and evolution of the free surface vortex. A Particle Image Velocimetry （PIV） was used to measure the free surface vortex flow field at different development stages. Flow visualization was used to locate the vortex position and find its structure. Empirical formulas about the critical submergence and the whole field structure were obtained. It is found that the tangential velocity distribution is similar to that of the Rankine vortex and the radial velocity changes little in the vortex functional scope. Vortex starts from the free surface and gradually intensifies to air entrainment vortex. The vortex core moves during the formation and evolution of the free surface vortex. Based on the experimental model, the vortex position and structure were predicted by numerical simulation combined with a vortex model and compared with that of the experiments, which shows satisfactory agreement.