RANS and LES computations of the tip-leakage vortex for different gap widths

In hydraulic turbines, the tip-leakage vortex is responsible for flow instabilities and for promoting erosion due to cavitation. To better understand the tip vortex flow, Reynolds-averaged Navier–Stokes (RANS) and large eddy simulation (LES) computations are carried out to simulate the flow around a NACA0009 blade including the gap between the tip and the wall. The main focus of the study is to understand the influence of the gap width on the development of the tip vortex, as for instance its trajectory. The RANS computations are performed using the open source solver OpenFOAM 2.1.0, two incidences and five gaps are considered. The LESs are achieved using the YALES2 solver for one incidence and two gaps.

The validation of the results is performed by comparisons with experimental data available downstream the trailing edge. The position of the vortex core, the mean velocity and the mean axial vorticity fields are compared at three different downstream locations. The results show that the mean behaviour of the tip vortex is well captured by the RANS and LES computations compared to the experiment. The LES results are also analysed to bring out the influence of the gap width on the development of the tip-leakage vortex. Finally, a law that matches the vortex trajectory from the leading edge to the mid-chord is proposed. Such a law can be helpful to determine, in case of cavitation, if the tip vortex will interact with the walls and cause erosion.     

On the use of immersed boundary methods for shock/obstacle interactions

This paper describes the implementation of immersed boundary method using the direct-forcing concept to investigate complex shock–obstacle interactions. An interpolation algorithm is developed for more stable boundary conditions with easier implementation procedure. The values of the fluid variables at the embedded ghost-cells are obtained using a local quadratic scheme which involves the neighboring fluid nodes. Detailed discussions of the method are presented on the interpolation of flow variables, direct-forcing of ghost cells, resolution of immersed-boundary points and internal treatment. The method is then applied to a high-order WENO scheme to simulate the complex fluid–solid interactions. The developed solver is first validated against the theoretical solutions of supersonic flow past triangular prism and circular cylinder. Simulated results for test cases with moving shocks are further compared with the previous experimental results of literature in terms of triple-point trajectory and vortex evolution. Excellent agreement is obtained showing the accuracy and the capability of the proposed method for solving complex strong-shock/obstacle interactions for both stationary and moving shock waves.     

湍流边界层中下扫流与“反发卡涡”

用氢气泡法观测湍流边界层的下扫流和有关的流动结构.实验中发现一种新型涡结构，它的特征与典型的发卡涡正好相反.发卡涡的头部指向下游，而它的头部指向上游； 发卡涡的两腿之间，由于涡的诱导产生上升流，而它则在其两腿之间，由于涡的诱导产生下扫流. 关键词： 湍流边界层 流动显示 流动结构 发卡涡     

An improved penalty immersed boundary method for fluid–flexible body interaction

An improved penalty immersed boundary (pIB) method has been proposed for simulation of fluid–flexible body interaction problems. In the proposed method, the fluid motion is defined on the Eulerian domain, while the solid motion is described by the Lagrangian variables. To account for the interaction, the flexible body is assumed to be composed of two parts: massive material points and massless material points, which are assumed to be linked closely by a stiff spring with damping. The massive material points are subjected to the elastic force of solid deformation but do not interact with the fluid directly, while the massless material points interact with the fluid by moving with the local fluid velocity. The flow solver and the solid solver are coupled in this framework and are developed separately by different methods. The fractional step method is adopted to solve the incompressible fluid motion on a staggered Cartesian grid, while the finite element method is developed to simulate the solid motion using an unstructured triangular mesh. The interaction force is just the restoring force of the stiff spring with damping, and is spread from the Lagrangian coordinates to the Eulerian grids by a smoothed approximation of the Dirac delta function. In the numerical simulations, we first validate the solid solver by using a vibrating circular ring in vacuum, and a second-order spatial accuracy is observed. Then both two- and three-dimensional simulations of fluid–flexible body interaction are carried out, including a circular disk in a linear shear flow, an elastic circular disk moving through a constricted channel, a spherical capsule in a linear shear flow, and a windsock in a uniform flow. The spatial accuracy is shown to be between first-order and second-order for both the fluid velocities and the solid positions. Comparisons between the numerical results and the theoretical solutions are also presented.     

激波冲击V形界面重气体导致的壁面与旋涡作用及其对湍流混合的影响

基于多组分混合物质量分数模型,采用色散最小耗散可控的高分辨率有限体积方法,数值模拟了弱激波冲击V形空气/SF_6界面后,界面不稳定性生成的旋涡与固体壁面作用问题.激波冲击V形界面之后,因斜压效应诱导涡量沉积在界面附近,形成沿界面规则排列的多个涡对结构.旋涡的诱导作用使界面不断变形和卷起,同时旋涡之间不断发生相互并对,诱导更多更小尺度的旋涡产生.旋涡诱导作用的叠加效应,使界面尖端处的初始涡对向上下壁面发展.随后,涡结构开始与壁面发生复杂的相互作用.旋涡与壁面作用后沿壁面加速,使得物质界面沿壁面伸展,随后,旋涡从壁面回弹,并诱导二次旋涡产生.旋涡与壁面相互作用的过程,能够明显加剧物质混合.本文从物质混合的角度研究了该过程的机理,分析了旋涡与壁面作用对物质混合的影响.     

Comparison between the 3D numerical simulation and experiment of the bubble near different boundaries

Based on the potential flow theory, the vortex ring is introduced to simulate the toroidal bubble, and the boundary element method is applied to simulate the evolution of the bubble. Elastic-plasticity of structure being taken into account, the interaction between the bubble and the elastic-plastic structure is computed by combining the boundary element method (BEM) and the finite element method (FEM), and a corresponding 3D computing program is developed. This program is used to simulate the three-dimensional bubble dynamics in free field, near wall and near the elastic-plastic structure, and the numerical results are compared with the existing experimental results. The error is within 10%. The effects of different boundaries upon the bubble dynamics are presented by studying the bubble dynamics near different boundaries. Supported by the National Natural Science Foundation of China (Grant No. 50779007), the National Science Foundation for Young Scientists of China (Grant No. 50809018), the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20070217074), the Defense Advanced Research Program of Science and Technology of Ship Industry (Grant No. 07J1.1.6), and Harbin Engineering University Foundation (Grant No. HEUFT07069)     

湍流热对流大尺度环流反转时的角涡特性

本文采用DNS方法计算二维方腔Rayleigh–Bénard热对流.在软湍流区热对流场呈现大尺度环流和两个反向转动的角涡,并出现了大尺度环流的反转现象.连续的温度等值线和流线图清晰地描述了反转现象的全过程.在反转过程中,角涡的大小尺度变化起到重要的作用.对角涡大小尺度变化的分析发现,在反转现象中其角涡尺度随时间的变化出现剧烈的振荡,而没有反转现象的热对流场中角涡尺度变化只有小幅的脉动.对反转过程前后的角涡大小尺度、典型位置速度及角点附近温度等流动特性进行了探讨和分析,发现反转是在瞬间完成的,角涡内速度脉动较小、温度较高,反转前角涡尺度与角涡侧壁垂向速度变化具有同步性.     

Development of Lattice Boltzmann Flux Solver for Simulation of Incompressible Flows

A lattice Boltzmann flux solver (LBFS) is presented in this work for simulation of incompressible viscous and inviscid flows. The new solver is based on Chapman-Enskog expansion analysis, which is the bridge to link Navier-Stokes (N-S) equations and lattice Boltzmann equation (LBE). The macroscopic differential equations are discretized by the finite volume method, where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers. The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. LBFS is validated by its application to simulate the viscous decaying vortex flow, the driven cavity flow, the viscous flow past a circular cylinder, and the inviscid flow past a circular cylinder. The obtained numerical results compare very well with available data in the literature, which show that LBFS has the second order of accuracy in space, and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.     

二维管道充填过程的修正SPH模拟

通过施加一种密度初始化方法对传统光滑粒子动力学(SPH)方法进行修正,提出一种修正SPH方法.同时,为了提高边界上数值计算的准确性,提出一种新的固壁边界处理方法.通过修正SPH方法模拟液滴拉伸问题和溃坝问题,验证修正SPH方法的准确性和可靠性.随后,对研究很少的管道充填过程进行修正SPH模拟,并讨论Re对流场及涡的影响.数值结果表明,修正SPH方法能够准确模拟牛顿流体管道充填过程,且流动受Re的影响较大.     

Numerical prediction of tip vortex cavitation behavior and noise considering nuclei size and distribution

The tip vortex cavitation behavior and sound generation were numerically analyzed. A numerical scheme combining Eulerian flow field computation and Lagrangian particle trace approach was applied to simulate tip vortex cavitation. Flow field was computed by using hybrid method which combines Reynolds-Averaged Navier-Stokes solver with Dissipation Vortex Model. The trajectory and behavior of each cavitation bubble were computed by Newton’s second law and Rayleigh-Plesset equation, respectively. According to nuclei population data, the cavitation nuclei were distributed and convected into the tip vortex flow. Calculated volume of the cavitation bubble and the trajectory were used as the input of cavitation bubble noise analysis. The relationship of cavitation inception, sound pressure level, and cavitation nuclei size was studied at several cavitation numbers. It was found that cavitation inception of smaller nuclei is more sensitive to the change of cavitation number and cavitation noise due to the cavitated smallest nuclei has the most influence on overall tip vortex cavitation noise.     

Intrinsic and extrinsic vortex nucleation mechanisms in the flow

We propose very general vortex nucleation mechanisms[1] analogous to a hydrodynamic instability and calculate associated critical velocity in agreement with experiments. The creation of vortices via extrinsic mechanism is driven by a formation of the surface vorticity sheet created by the flow, which reaches a critical size. Such a sheet screens an attraction of a half-vortex ring to the wall, the barrier for the vortex nucleation disappears and the vortex nucleation is started. In the intrinsic mechanism the creation of a big vortex ring, which transforms into the vortex, is driven by a fluctuative generation of small vortex rings.     

Vortex sheet in superfluid<Superscript>3</Superscript>He-A

The vortex sheet (VS) is the most unexpected discovery in rotating superfluids during the last ten years. Usually superfluids respond to rotation by creating an array of vortex lines, which are parallel to the rotation axis, and the circulation around them is quantized. In the VS the vorticity is located on a 2 dimensional sheet that folds to equidistant layers in a rotating container. The VS is one out of five stable vortex structures in³He-A. The stability of the VS in³He-A arises from a special structure, which consists of a nonsingular vorticity bound to a topologically stable domain wall. The vortex sheet forms the equilibrium state of³He-A at rotation velocities larger than ～3 rad/s, but it is also created as metastable state at lower velocities. Experimentally the vortex sheet is distinguished from its NMR response.     

Eulerian simulation of sedimentation flows in vertical and inclined vessels

Sedimentation of particles in inclined and vertical vessels is numerically simulated using a finite volume method where the Eulerian multiphase model is applied. The particulate phase as well as the fluid phase is regarded as a continuum while the viscosity and solid stress of the particulate phase are modelled by the kinetic theory of granular flows. The numerical results show an interesting phenomenon of the emergence of two circulation vortices of the sedimentation flow in a vertical vessel but only one in the inclined vessel. Several sensitivity tests are simulated to understand the factors that influence the dual-vortex flow structure in vertical sedimentation. Results show that a larger fluid viscosity makes the two vortex centres much closer to each other and the boundary layer effect at lateral walls is the key factor to induce this phenomenon. In the fluid boundary layer particles settle down more rapidly and drag the local carrier fluid to flow downward near the lateral walls and thus form the dual-vortex flow pattern.     

Lattice-Boltzmann Simulations of Fluid Flows in MEMS

The lattice Boltzmann model is a simplified kinetic method based on the particle distribution function. We use this method to simulate problems in MEMS, in which the velocity slip near the wall plays an important role. It is demonstrated that the lattice Boltzmann method can capture the fundamental behaviors in micro-channel flow, including velocity slip, nonlinear pressure drop along the channel and mass flow rate variation with Knudsen number. The Knudsen number dependence of the position of the vortex center and the pressure contour in micro-cavity flows is also demonstrated.     

Intrinsic and extrinsic vortex nucleation mechanisms in the flow

We propose very general vortex nucleation mechanisms[1] analogous to a hydrodynamic instability and calculate associated critical velocity in agreement with experiments. The creation of vortices via extrinsic mechanism is driven by a formation of the surface vorticity sheet created by the flow, which reaches a critical size. Such a sheet screens an attraction of a half-vortex ring to the wall, the barrier for the vortex nucleation disappears and the vortex nucleation is started. In the intrinsic mechanism the creation of a big vortex ring, which transforms into the vortex, is driven by a fluctuative generation of small vortex rings. Work supported by NORDITA and Landau Institute     

Thermal fluctuations of vortex matter in trapped bose-einstein condensates

We perform Monte Carlo studies of vortices in three dimensions in a cylindrical confinement, with uniform and nonuniform density. The former is relevant to rotating 4He; the latter is relevant to a rotating trapped Bose-Einstein condensate. In the former case, we find dominant angular thermal vortex fluctuations close to the cylinder wall. For the latter case, a novel effect is that at low temperatures the vortex solid close to the center of the trap crosses directly over to a tensionless vortex tangle near the edge of the trap. At higher temperatures an intermediate tensionful vortex liquid located between the vortex solid and the vortex tangle may exist.     

Wave-induced bed flows by a Lagrangian vortex scheme

Nominally 2-dimensional viscous flow induced by gravity waves over a spatially periodic bed is simulated by a Lagrangian vortex scheme. A vortex sheet is introduced on the surface at each time step to satisfy the zero velocity conditions. The sheet is discretised; the vortex-in-cell method is used to convect vorticity and random walks are added to effect viscous diffusion. Good agreement with analytical theory is obtained for velocity profiles in uniform sinusoidal flow and for mass transport due to linear waves. Mass transport for finite amplitude waves is also obtained. For separated flow over rippled beds, which is still liminar, a vortex decay factor is required to produce agreement with experiment and is thought to compensate for large scale 3-dimensional effects.     

A high resolution spatially adaptive vortex method for separating flows. Part I: Two-dimensional domains

A grid-free high-resolution spatially-adaptive vortex method for two-dimensional incompressible flow in bounded domains is presented. The computational algorithm is based on operator splitting in which convection and diffusion are handled separately every time step. In the convection step, computational elements are convected with velocities obtained by fast approximations of the Biot–Savart superposition with second-order Runge–Kutta time integration scheme. Diffusion is performed using the smooth redistribution method that employs a Gaussian basis function for vorticity in the interior. Near solid walls, the core functions are modified to conserve circulation. The no-slip boundary condition is enforced by creating of a vortex sheet that is redistributed to neighboring elements using the redistribution method. The proposed method enables accurate and smooth recovery of the vorticity and does not require explicit use of vortex images or occasional re-meshing. Algorithms for reduction in computational cost by accurately removing elements in overcrowded regions and for spatial adaptivity that allows for variable core sizes and variable element spacing are presented. Computations of flow around an impulsively started cylinder for Reynolds number values of 1000, 3000, and 9500 are preformed to investigate various aspects of the proposed method.     

Accurate computation of Stokes flow driven by an open immersed interface

We present numerical methods for computing two-dimensional Stokes flow driven by forces singularly supported along an open, immersed interface. Two second-order accurate methods are developed: one for accurately evaluating boundary integral solutions at a point, and another for computing Stokes solution values on a rectangular mesh. We first describe a method for computing singular or nearly singular integrals, such as a double layer potential due to sources on a curve in the plane, evaluated at a point on or near the curve. To improve accuracy of the numerical quadrature, we add corrections for the errors arising from discretization, which are found by asymptotic analysis. When used to solve the Stokes equations with sources on an open, immersed interface, the method generates second-order approximations, for both the pressure and the velocity, and preserves the jumps in the solutions and their derivatives across the boundary. We then combine the method with a mesh-based solver to yield a hybrid method for computing Stokes solutions at N² grid points on a rectangular grid. Numerical results are presented which exhibit second-order accuracy. To demonstrate the applicability of the method, we use the method to simulate fluid dynamics induced by the beating motion of a cilium. The method preserves the sharp jumps in the Stokes solution and their derivatives across the immersed boundary. Model results illustrate the distinct hydrodynamic effects generated by the effective stroke and by the recovery stroke of the ciliary beat cycle.     

过渡流下近壁插入圆柱对壁面强化传热的影响

本文利用基于复合网格系统的计算方法,对Re=50～1200的近壁插入圆柱流场进行数值模拟,研究过渡流状态下在壁面附近插入圆柱对下游壁面传热强化的影响。并基于低速循环水槽流动实验台,采用粒子成像测试法(PIV)对Re=100～500的近壁插入圆柱流场进行可视化实验研究,验证了数值模拟方法的可靠性。研究结果表明:近壁插入圆柱流场在Re=100时进入过渡流状态;Re直接影响圆柱尾流中周期性涡脱和壁面涡岛的发生位置及其洗刷效应的大小,随着Re的增大,洗刷效应明显增强,因而,过渡流范围内Re越大,圆柱下游壁面传热强化越大。