Unsteady characteristics of the cavity flow around a rough circular cylinder

Effects of surface roughness on the unsteady cavitating flow around a two-dimensional circular cylinder were experimentally investigated at Reynolds numbers from 1.36 × 10⁵ to 1.78 × 10⁵. Two patterns of surface roughness were investigated, a double-cut pattern and a single-cut pattern. The cavity elongates with an increase of the surface roughness, especially in supercavitating flow. However, for some roughness parameters tested, the cavity length exhibits an extreme decrease. In a particular case of the double-cut pattern, there exists the minimum cavity behind the cylinder.     

基于重叠网格与结构网格的圆柱绕流数值模拟

为研究重叠网格与结构网格在圆柱绕流数值模拟中的区别,以二维圆柱为例,利用有限元分析软件ANSYS 19.2中的DM与Mesh建立模型并划分重叠网格,利用ANSYS 19.2中的ICEM建立模型并划分结构网格。采用FLUENT 19.2中laminar模型模拟分析系统中的平均升力系数、平均阻力系数、斯特劳哈尔数St等流体动力特性。通过改变流体流速得到两种不同网格下各6组雷诺数Re,这6组雷诺数在60~160之间。结果表明:结构网格与重叠网格的St都随着Re的增加而增加,但相同雷诺数下重叠网格对应的St数值更大,St的增长速度更快;重叠网格与结构网格的平均升力系数与平均阻力系数随着Re的增加趋于稳定的速度都加快,但结构网格的平均升力系数与平均阻力系数趋于稳定的速度更快,且两种网格的平均升力系数与平均阻力系数趋于稳定速度的差距逐渐缩小,当Re=160时,两种网格的平均升力系数与平均阻力系数趋于稳定的速度几乎相同;当雷诺数在60~160之间时,采用重叠网格计算出来的斯特劳哈尔数比结构网格更加接近理论值;从升力功率谱密度分布曲线中可以看出,随着雷诺数的增加,两种网格下的频率逐渐变大,并且相同雷诺数下重叠网格的频率比结构网格大。     

Mixed convection boundary layer flow of a viscoelastic fluid over a horizontal circular cylinder

The steady mixed convection boundary layer flow of a viscoelastic fluid over a horizontal circular cylinder in a stream flowing vertically upwards is numerically studied for both cases of heated and cooled cylinders. The governing partial differential equations are transformed into dimensionless forms using an appropriate transformation and then solved numerically using the Keller-box method. The comparison between the solutions obtained and those for a Newtonian fluid is found to be very good. Effects of the mixed convection and elasticity parameters on the skin friction and heat transfer coefficients for a fluid having the Prandtl number equal to one are also discussed. It is found that for some values of the viscoelastic parameter and some negative values of the mixed convection parameter (opposing flow) the boundary layer separates from the cylinder. Heating the cylinder delays separation and can, if the cylinder is warm enough, suppress the separation completely. Similar to the case of a Newtonian fluid, cooling the cylinder brings the separation point nearer to the lower stagnation point. However, for a sufficiently cold cylinder there will not be a boundary layer.     

The effect of particles on fluid turbulence in a turbulent boundary layer over a cylinder

The turbulent fluid and particle interaction in the turbulent boundary layer for cross flow over a cylinder has been experimentally studied. A phase-Doppler anemometer was used to measure the mean and fluctuating velocities of both phases. Two size ranges of particles (30μm–60μm and 80μm–150μm) at certain concentrations were used for considering the effects of particle sizes on the mean velocity profiles and on the turbulent intensity levels. The measurements clearly demonstrated that the larger particles damped fluid turbulence. For the smaller particles, this damping effect was less noticeable. The measurements further showed a delay in the separation point for two phase turbulent cross flow over a cylinder. The project supported by the National Natural Science Foundation of China     

Stress boundary layers in the viscoelastic flow past a cylinder in a channel： limiting solutions

The flow past a cylinder in a channel with the aspect ratio of 2:1 for the upper convected Maxwell (UCM) fluid and the Oldroyd-B fluid with the viscosity ratio of 0.59 is studied by using the Galerkin/Least-square finite element method and a p-adaptive refinement algorithm. A posteriori error estimation indicates that the stress-gradient error dominates the total error. As the Deborah number, D_e, approaches 0.8 for the UCM fluid and 0.9 for the Oldroyd-B fluid, strong stress boundary layers near the rear stagnation point are forming, which are characterized by jumps of the stress-profiles on the cylinder wall and plane of symmetry, huge stress gradients and rapid decay of the gradients across narrow thicknesses. The origin of the huge stress-gradients can be traced to the purely elongational flow behind the rear stagnation point, where the position at which the elongation rate is of 1/2D_e approaches the rear stagnation point as the Deborah number approaches the critical values. These observations imply that the cylinder problem for the UCM and Oldroyd-B fluids may have physical limiting Deborah numbers of 0.8 and 0.9, respectively.The project supported by the National Natural Science Foundation of China (50335010 and 20274041) and the MOLDFLOW Comp. Australia.     

Instability of the separated shear layer in flow past a cylinder: Forced excitation

The receptivity of the separated shear layer for Re = 300 flow past a cylinder is investigated by forced excitation via an unsteady inflow. In order to isolate the shear layer instability, a numerical experiment is set up that suppresses the primary wake instability. Computations are carried out for one half of the cylinder, in two dimensions. The flow past half a cylinder with steady inflow is found to be stable for Re = 300. However, an inlet flow with pulsatile perturbations, of amplitude 1% of the mean, results in the excitation of the shear layer mode. The frequency of the perturbation of the inlet flow determines the frequency associated with the shear layer vortices. For a certain range of forced frequencies the recirculation region undergoes a low‐frequency longitudinal contraction and expansion. An attempt is made to relate this instability to a global mode of the wake determined from a linear stability analysis. Interestingly, this phenomenon disappears when the outflow boundary of the computational domain is shifted sufficiently downstream. This study demonstrates the need of carefully investigating the effect of the location of outflow boundaries if the computational results indicate the presence of low‐frequency fluctuations. The effect of Re and amplitude of unsteadiness at the inlet are also presented. All computations have been carried out using a stabilized finite element formulation of the incompressible flow equations. Copyright © 2007 John Wiley & Sons, Ltd.     

Laminar-turbulent transition in the boundary layer on cones in a hypersonic flow at high Reynolds numbers per meter

The laminar-turbulent transition is experimentally studied in boundary-layer flows on cones with a rectangular axisymmetric step in the base part of the cone and without the step. The experiments are performed in an A-1 two-step piston-driven gas-dynamic facility with adiabatic compression of the working gas with Mach numbers at the nozzle exit M _∞ = 12–14 and pressures in the settling chamber P₀ = 60–600 MPa. These values of parameters allow obtaining Reynolds numbers per meter near the cone surface equal to Re _1e = (53–200) · 10⁶ m ⁻¹. The transition occurs at Reynolds numbers Re _tr = (2.3–5.7) · 10⁶. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 76–83, May–June, 2007.     

Numerical and experimental study of the fine structure of a stratified fluid flow over a circular cylinder

The fine structure of the flow field of a continuously stratified fluid around a circular cylinder for small values of the Froude number was investigated in laboratory and numerical experiments. The parameters of the leading perturbation, the internal-wave field, and the cylinder wake were calculated using a two-dimensional model. The existence of the previously experimentally observed high-gradient density layers in the wake that are parallel to the flow axis was for the first time confirmed by numerical calculations. Results of the numerical and experimental studies are in good agreement with each other and with analytical models for small values of the Froude number. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 43–54, January–February, 2007.     

Numerical modeling of the turbulent viscous-gas flow past a cylinder with a circular vortex cell in the presence of suction from the central body surface

The compressibility effect on the cylinder drag reduction due to air suction through the surface of a central body in a circular vortex cell is estimated on the basis of the solution of the steady Reynolds equations closed by the shear stress transfer model, together with the continuity, energy, and state equations.     

Efficient computation of mean drag for the subcritical flow past a circular cylinder using general Galerkin G2

General Galerkin (G2) is a new computational method for turbulent flow, where a stabilized Galerkin finite element method is used to compute approximate weak solutions to the Navier–Stokes equations directly, without any filtering of the equations as in a standard approach to turbulence simulation, such as large eddy simulation, and thus no Reynolds stresses are introduced, which need modelling. In this paper, G2 is used to compute the drag coefficient c_D for the flow past a circular cylinder at Reynolds number Re=3900, for which the flow is turbulent. It is found that it is possible to approximate c_D to an accuracy of a few percent, corresponding to the accuracy in experimental results for this problem, using less than 10⁵ mesh points, which makes the simulations possible using a standard PC. The mesh is adaptively refined until a stopping criterion is reached with respect to the error in a chosen output of interest, which in this paper is c_D. Both the stopping criterion and the mesh‐refinement strategy are based on a posteriori error estimates, in the form of a space–time integral of residuals times derivatives of the solution of a dual problem, linearized at the approximate solution, and with data coupling to the output of interest. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical investigation on vortex-induced vibration of an elastically mounted circular cylinder at low Reynolds number using the fictitious domain method

A direct numerical simulation of two-dimensional (2D) flow past an elastically mounted circular cylinder at low Reynolds number using the fictitious domain method had been undertaken. The cylinder motion was modelled by a two degree-of-freedom mass–spring–damper system. The computing code was verified against a benchmark problem in which flow past a stationary circular cylinder is simulated. Then, analyses of vortex-induced vibration (VIV) responses, drag and lift forces and the phase and vortex structures were carried out. Results show that the cylinder's non-dimensional cross-flow response amplitude reaches its summit of 0.572 in the ‘lock-in’ regime. The ‘2S’, instead of the ‘2P’, vortex shedding mode is dominated in the ‘lower’ branch for this 2D low-Re VIV. A secondary oscillation is observed in the lift force when ‘lock-in’ occurs. It is shown that this secondary component changes the phase, offset the energy input by the primary component and thus reduces the cylinder responses. Effects of the Skop–Griffin parameter on cylinder responses were also investigated.     

Effect of the boundary layer on the base pressure in a two-dimensional flow with a Mach number M = 5

New data on the base pressure in a two-dimensional ow with a Mach number M = 5 are obtained for a wide range of variation of the normalized boundary-layer thickness in the flow-separation cross section. The test results are compared with Tanner’s theory, and a conclusion is made that this numerical model has to be corrected. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 23–28, May–June, 2005.     

On the effect of the nonparallelism of the boundary layer flow over a flat plate on its stability characteristics

The effect of the nonparallelism of the boundary layer flow over a flat plate on its stability characteristics has been investigated by several authors, and it was claimed that the results of the theoretical calculations are already in good agreement with the experimental observations. However, this is not true. In this paper, this problem is reinvestigated, using two different methods. It is found that within the framework of linear theory, the theoretical results are in fact not in good agreement with the experimental observations. To settle this problem, nonlinear effect must be taken into consideration. Projects Supported by the Science Fund of the Chinese Academy of Sciences     

Numerical investigation of the effect of convex transverse curvature and concave grooves on the turbulent boundary layer along a cylinder in axial flow

Axisymmetric turbulent boundary layers that develop around streamwise oriented long cylinder-like objects can be found in many applications, such as towed array sonars or marine seismic streamers. In many of these applications, turbulent fluctuations within the boundary layer flow can have a negative impact compared with laminar flow conditions. The aim of the present work is to design a surface modification that influences the turbulent boundary layer around a cylinder in axial flow in order to reduce turbulent fluctuations. To design the surface we consider recent findings regarding the turbulence damping effects of groove-like surface structures and combine these insights with the effect of convex transverse curvature on turbulence. We use large-eddy simulations to investigate the flow around a cylinder of modified design and around a reference circular cylinder. Both flows have a radius-based Reynolds number of ${\text{Re}}_a\approx 1.23\cdot 10^4$. The modified design leads to a 20 % decrease in the average wall shear stress and results in local reductions in the turbulent intensities, Reynolds stress, the temporal velocity spectrum, and the turbulent dissipation rate. The analysis within the anisotropy-invariant space reveals a tendency towards flow relaminarization. However, the new design has no effect on turbulent pressure fluctuations. We provide suggestions on how to further improve the surface design to achieve even greater flow stabilization.     

Numerical study on the dynamics of droplet passing through a cylinder obstruction in confined microchannel flow

The droplet dynamics passing through a cylinder obstruction was investigated with direct numerical simulations with FE-FTM (Finite Element-Front Tracking Method). The effect of droplet size and capillary number (Ca) was studied for both Newtonian and viscoelastic fluids. In the case of Newtonian droplet immersed in Newtonian medium, the droplet breakup induced by the geometric hindrance depends on the droplet size. As Ca increases, the short droplets (1.3 times longer than the channel width) break up while passing through the obstruction. However, the breakup does not occur for longer droplets (1.8 times longer than the channel width). When the viscoelastic fluid characterized by the Oldroyd-B model is considered, the Newtonian droplet immersed in viscoelastic medium breaks up into two smaller droplets while passing through the cylinder obstruction with increasing De_m (Deborah number of the medium). We also show that the normal stress difference plays a key role on the droplet breakup and the droplet extension. The normal stress difference is enhanced in the negative wake region due to the droplet flow, which also promotes droplet extension in that region. This numerical study provides information not only on underlying physics of the droplet flows passing through a cylinder obstruction but also on the useful guidelines for microfluidic applications.