棱柱绕流流动的数值模拟

利用数值方法对长宽比为1／3, 1和3的棱柱绕流在雷诺数为100的非稳态流动特性进行了分析和研究。采用有限体积法对棱柱绕流的二维流动N-S方程进行离散求解,分析和研究了非稳态的棱柱绕流流场,升力系数,阻力系数和涡动特性,数值模拟的结果与相关文献的数据比较吻合。通过上述研究能够为了解棱柱绕流的非稳态流动特性提供有力的帮助。而对棱柱三维流动的模拟分析和对雷诺数的变化对棱柱流动特性的影响进行研究,将为掌握棱柱绕流的工程特性打下基础。     

Numerical investigation of flow and combined natural-forced convection from an isothermal square cylinder in cross flow

Unsteady flow and heat transfer from a horizontal isothermal square cylinder is studied numerically using a three-dimensional computational model to investigate the influence of buoyancy on the forced flow and heat transfer characteristics. The numerical model is based on a horizontal square cylinder subjected to laminar fluid flow in an unconfined channel. The governing equations in 3D form are solved using a fractional step method based on the finite difference discretization in addition to a Crank–Nicholson scheme employed to the convective and the viscous terms. Two working fluids–air (Pr = 0.7) and water (Pr = 7)–are considered, and the flow and heat transfer simulations were carried out for the Reynolds and Richardson numbers in the intervals 55 ≤ Re ≤ 250 and 0 ≤ Ri ≤ 2, respectively. The flow characteristics such as time-averaged drag/lift, rms drag/rms lift coefficients as well as Strouhal number were computed. The heat transfer from the cylinder is assessed by mean Nusselt number (and rms Nusselt number) over the total heated cylinder walls. As the buoyancy increases, the mass and the velocity of the fluid flowing underneath the cylinder increases. The fluid is injected into the near wake region with an upward motion which significantly alters the flow field in the downstream as well as upstream regions. The effects of Reynolds, Richardson and Prandtl numbers on the flow field and temperature distributions are discussed in detail. It is shown that the flow and heat transfer characteristics are influenced more for air than water. To fill the void in the literature, useful empirical correlations of practical importance are derived for pure forced and pure natural as well as mixed convection. The mixed convection correlations, in terms of the ratio of pure forced convection, are also developed, and their implications are discussed.     

Milu-CG method and the numerical study on the flow around a rotating circular cylinder

I.IntroductionTilenowaroundarotatitlgcircularcylinderisacomplexunsteadyone.ItincludesmanycomplicatedtlowphenomenaSuchastheunsteadyboundarylayerseparation,thegenerationandsheddingofvorticesandtheinteractionwitllwakesetc..Therotationofacircularcylillderarounditsaxiswilldecreaseandsuppresstheflowseparationandvortexsheddingononesideofthecylinder,whileincreasinganddevelopingonanothel,side.Atransverseliftforcewillactonthecylinder,andthisphenomenoniscalledtheMagnuseffect.Themost.importantparameterf…     

A coupled continuous and discontinuous finite element method for the incompressible flows

In this paper, we develop a coupled continuous Galerkin and discontinuous Galerkin finite element method based on a split scheme to solve the incompressible Navier–Stokes equations. In order to use the equal order interpolation functions for velocity and pressure, we decouple the original Navier–Stokes equations and obtain three distinct equations through the split method, which are nonlinear hyperbolic, elliptic, and Helmholtz equations, respectively. The hybrid method combines the merits of discontinuous Galerkin (DG) and finite element method (FEM). Therefore, DG is concerned to accomplish the spatial discretization of the nonlinear hyperbolic equation to avoid using the stabilization approaches that appeared in FEM. Moreover, FEM is utilized to deal with the Poisson and Helmholtz equations to reduce the computational cost compared with DG. As for the temporal discretization, a second‐order stiffly stable approach is employed. Several typical benchmarks, namely, the Poiseuille flow, the backward‐facing step flow, and the flow around the cylinder with a wide range of Reynolds numbers, are considered to demonstrate and validate the feasibility, accuracy, and efficiency of this coupled method. Copyright © 2016 John Wiley & Sons, Ltd.     

幂律流体圆柱绕流的格子波尔兹曼模拟简

基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明：插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中.     

Simulation of a Viscous Flow Around an Oscillating Cylinder at Low Keulegan-Carpenter Numbers

A finite difference simulation method for a viscous flow around a circular cylinder sinusoidally oscillating at low Keulegan-Carpenter numbers is presented. Navier-Stokes equations in finite difference form are solved on a moving grid system, based on a time dependent coordinate transformation. Evolution with time of the flow structures induced by a circular cylinder performing sinusoidal oscillations in a fluid at rest, by means of stream lines, pressure contours and vortex shedding is studied in detail at Keulegan-Carpenter numbers, Kc = 9.4 and 14. The time dependent drag and lift are also explained.     

Development of a hybrid finite volume/element solver for incompressible flows

In this paper, we report our development of an implicit hybrid flow solver for the incompressible Navier–Stokes equations. The methodology is based on the pressure correction or projection method. A fractional step approach is used to obtain an intermediate velocity field by solving the original momentum equations with the matrix‐free implicit cell‐centred finite volume method. The Poisson equation derived from the fractional step approach is solved by the node‐based Galerkin finite element method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centres and the auxiliary variable at cell vertices, making the current solver a staggered‐mesh scheme. Numerical examples demonstrate the performance of the resulting hybrid scheme, such as the correct temporal convergence rates for both velocity and pressure, absence of unphysical pressure boundary layer, good convergence in steady‐state simulations and capability in predicting accurate drag, lift and Strouhal number in the flow around a circular cylinder. Copyright © 2007 John Wiley & Sons, Ltd.     

幂律流体圆柱绕流的格子波尔兹曼模拟

基于插值补充格子波尔兹曼方法和幂律流体的本构方程,建立了贴体坐标系下适用于幂律流体的格子波尔兹曼模型,模拟了幂律流体的圆柱绕流问题,采用非平衡外推格式处理圆柱表面的速度无滑移边界,利用应力积分法确定曳力系数和升力系数,并与基于标准的格子波尔兹曼方法和有限容积法获得的数值数据进行对比,吻合良好. 进行了网格无关性验证之后,分析了稳态流动时,不同雷诺数下幂律指数对于尾迹长度、分离角、圆柱表面黏度分布、表面压力系数及曳力系数的影响,以及非定常流动中,幂律指数对于流场、曳力系数、升力系数和斯特劳哈尔数的影响. 获得的变化规律与基于其他数值模拟方法得到的结果相一致,充分验证了模型的有效性和正确性. 结果表明:插值补充格子波尔兹曼方法可以用来模拟幂律流体在具有复杂边界流场内的流动问题,通过引入不同的非牛顿流体本构方程,该方法还可以进一步应用于其他类型的非牛顿流体研究中.      

定常流作用于波浪形边界附近的圆柱上的水动力

本文介绍在水洞中做的作用于波浪形边界附近的水动力特性的实验研究。在基于圆柱模型直径的 R_c=10~4～1.9×10~4 范围内,测量了圆柱在波谷、波峰和不同距离上的阻力、升力脉动变化的频率。流谱显示实验揭示了尾流结构随距离-直径比 G/D 的变化情况及圆柱与边界相互作用的机制。     

流向振荡圆柱绕流的格子Boltzmann方法模拟

用一种新近发展起来的格子Boltzmann方法(LBM)在相对较小的雷诺数(Re \le 200)条件下模拟了不可压缩的流向振荡圆柱绕流问题, 考查了涡脱落模态和升阻力特性. 通过模拟, 在近尾流区发现了实验研究中已经发现的对称/反对称的涡脱落模态, 包括有些传统数值方法未发现的模态. 研究了频率锁定区域的范围及其与振幅的关系, 发现振幅越大, 发生锁定的频率区域越宽. 此外还对升阻力进行了定量意义的模拟,研究了振荡频率和振幅与升阻力的关系.      

On the development of lift and drag in a rotating and translating cylinder

The two-dimensional flow around a rotating cylinder is investigated numerically using a vorticity forces formulation with the aim of analyzing quantitatively the flow structures, and their evolutions, that contribute to the lift and drag forces on the cylinder. The Reynolds number considered, based on the cylinder diameter and steady free stream speed, is Re=200, while the non-dimensional rotation rate (ratio of the surface speed and free stream speed) selected was α=1 and 3. For α=1 the wake behind the cylinder for the fully developed flow is oscillatory due to vortex shedding, and so are the lift and drag forces. For α=3 the fully developed flow is steady with constant (high) lift and (low) drag. Each of these cases is considered in two different transient problems, one with angular acceleration of the cylinder and constant speed, and the other one with translating acceleration of the cylinder and constant rotation. We characterize quantitatively the contributions of individual fluid elements (vortices) to aerodynamic forces, explaining and quantifying the mechanisms by which the lift is generated in each case. In particular, for high rotation (when α=3), we explain the relation between the mechanisms of vortex shedding suppression and those by which the lift is enhanced and the drag is almost suppressed when the fully developed flow is reached.     

Electromagnetic control of seawater flow around circular cylinders

We investigate numerically the electromagnetic control of seawater flows over an infinitely long circular cylinder. Stripes of electrodes and magnets, wrapped around the cylinder surface, produce a tangential body force (Lorentz force) that stabilizes the flow. This mechanism delays flow separation, reduces drag and lift, and finally suppresses the von Kármán vortex street. Results from two-dimensional simulations of the Navier–Stokes equations in a range 10<Re<300 and Lorentz force calculations are presented. Emphasis is placed on the disclosure of physical phenomena as well as a quantitative detection of the flow field and forces. It is shown that the drag strongly depends on the geometry of the electromagnetic actuator and on its location at the cylinder surface. The effect of flow control increases with larger Reynolds numbers, since the boundary layer thickness and the penetration depth of the Lorentz force are closely connected.     

Shape optimisation of an oscillating body in fluid flow by adjoint equation and ALE finite element methods

The purpose of this paper is to determine the shape of an oscillating body by minimising drag and lift forces, located in a transient incompressible viscous fluid flow by means of the Arbitrary Lagrangian Eulerian finite element method and an optimal control theory. A performance function is expressed by the drag and lift forces. The performance function should be minimised satisfying the state equation and the constant volume condition. Therefore, this problem can be transformed into a minimisation problem without constraint by the Lagrange multiplier method. The adjoint equation and the gradient of the performance function are used to update the shape of the body. In this study, as a minimisation technique, the weighted gradient method is applied. The final shape is obtained of which drag and lift forces are reduced by 66.2% and 92.8%, respectively. The final shape obtained by this study is compared with the final shape of the non-oscillating body. The obtained final shape of the oscillating body is significantly different from the non-oscillating body.     

不可压缩粘性流动的CBS有限元解法

对于二维不可压缩粘性流动,首先通过坐标变换的方式得到了的不含对流项的NS方程,并给出了CBS有限元方法求解的一般过程。结合一类同时含有压力和速度的出口边界条件,对方腔顶盖驱动流、后向台阶绕流和圆柱绕流进行了计算。所得结果与基准解符合良好,验证了CBS算法对于定常、非定常粘性不可压缩流动问题的可行性和所用出口边界条件的无反射特性。特别的,对于圆柱绕流,Re=100时非定常升、阻力系数及漩涡脱落等非定常都得到了较好地模拟,为一进步研究自激振动等更加复杂的非定常流动问题奠定了基础。     

Numerical simulation of Reynolds number effects on velocity shear flow around a circular cylinder

Three-dimensional Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are performed to investigate the shear effects on flow around a circular cylinder at Reynolds numbers of Re=60–1000. The shear parameter, β, which is based on the velocity gradient, cylinder diameter and upstream mean velocity at the center plane of the cylinder, varies from 0 to 0.30. Variations of Strouhal number, drag and lift coefficients, and unsteady wake structures with shear parameter are studied, along with their dependence on Reynolds number. The presented simulation provides detailed information for the flow field around a circular cylinder in shear flow. This study shows that the Strouhal number exhibits no significant variation with shear parameter. The stagnation point moves to the high-velocity side almost linearly with shear parameter, and this result mainly influences the aerodynamic forces acting on a circular cylinder in shear flow. Both the Reynolds number and shear parameter influence the movement of the stagnation point and separation point. Mode A wake instability is suppressed into parallel vortex shedding mode at a certain shear parameter. The lift force increases with increasing shear parameter and acts from the high-velocity side to the low-velocity side. In addition, a simple method to estimate the lift force coefficient in shear flow is provided.     

Numerical simulation of natural and mixed convection flows by Galerkin‐characteristic method

A numerical investigation is performed to study the solution of natural and mixed convection flows by Galerkin‐characteristic method. The method is based on combining the modified method of characteristics with a Galerkin finite element discretization in primitive variables. It can be interpreted as a fractional step technique where convective part and Stokes/Boussinesq part are treated separately. The main feature of the proposed method is that, due to the Lagrangian treatment of convection, the Courant–Friedrichs–Lewy (CFL) restriction is relaxed and the time truncation errors are reduced in the Stokes/Boussinesq part. Numerical simulations are carried out for a natural convection in squared cavity and for a mixed convection flow past a circular cylinder. The computed results are compared with those obtained using other Eulerian‐based Galerkin finite element solvers, which are used for solving many convective flow models. The Galerkin‐characteristic method has been found to be feasible and satisfactory. Copyright © 2006 John Wiley & Sons, Ltd.     

等边布置三圆柱绕流的数值分析

用有限元方法对于等边三角形布置的三个相同直径的二维圆柱绕流问题进行了数值模拟，分别求得了在不同间距比下的流场分布和各圆柱的升、阻力系数以及斯托哈罗数。计算结果表明，小间距比情形下三圆柱之间的干扰是严重的，流动并不对称于中心轴线，而是偏向下游的某个圆柱。数值计算结果与有关文献和实验进行了对比，结果令人满意。     

Rolling/sliding of a particle on a flat wall in a linear shear flow at finite Re

Recently Lee and Balachandar proposed analytically-based expressions for drag and lift coefficients for a spherical particle moving on a flat wall in a linear shear flow at finite Reynolds number. In order to evaluate the accuracy of these expressions, we have conducted direct numerical simulations of a rolling particle for shear Reynolds number up to 100. We assume that the particle rolls on a horizontal flat wall with a small gap separating the particle from the wall (L = 0.505) and thus avoiding the logarithmic singularity. The influence of the shear Reynolds number and the translational velocity of the particle on the hydrodynamic forces of the particle was investigated under both transient and the final drag-free and torque-free steady state. It is observed that the quasi-steady drag and lift expressions of Lee and Balachandar provide good approximation for the terminal state of the particle motion ranging from perfect sliding to perfect rolling. With regards to transient particle motion in a wall-bounded shear flow it is observed that the above validated quasi-steady drag and lift forces must be supplemented with appropriate wall-corrected added-mass and history forces in order to accurately predict the time-dependent approach to the terminal steady state. Quantitative comparison with the actual particle motion computed in the numerical simulations shows that the theoretical models quite effective in predicting rolling/sliding motion of a particle in a wall-bounded shear flow at moderate Re.     

Numerical study on the suppression of the vortex-induced vibration of an elastically mounted cylinder by a traveling wave wall

In the present paper, the commercial CFD code “Fluent” was employed to perform 2-D simulations of an entire process that included the flow around a fixed circular cylinder, the oscillating cylinder (vortex-induced vibration, VIV) and the oscillating cylinder subjected to shape control by a traveling wave wall (TWW) method. The study mainly focused on using the TWW control method to suppress the VIV of an elastically supported circular cylinder with two degrees of freedom at a low Reynolds number of 200. The cross flow (CF) and the inline flow (IL) displacements, the centroid motion trajectories and the lift and drag forces of the cylinder that changed with the frequency ratios were analyzed in detail. The results indicate that a series of small-scale vortices will be formed in the troughs of the traveling wave located on the rear part of the circular cylinder; these vortices can effectively control the flow separation from the cylinder surface, eliminate the oscillating wake and suppress the VIV of the cylinder. A TWW starting at the initial time or at some time halfway through the time interval can significantly suppress the CF and IL vibrations of the cylinder and can remarkably decrease the fluctuations of the lift coefficients and the average values of the drag coefficients; however, it will simultaneously dramatically increase the fluctuations of the drag coefficients.