串列双圆柱绕流问题的数值模拟

本文运用有限体积方法,对绕串列放置的双圆柱的二维不可压缩流动进行了数值计算。为研究两圆柱不同间距对圆柱相互作用和尾流特征的影响,选取间距比Ｌ／Ｄ（Ｌ为两圆柱中心间的距离,Ｄ为圆柱直径）在１．５～５．０之间每隔０．５共八个有代表性的间距进行了计算模拟。计算均在Ｒｅ＝２００条件下进行。计算结果表明：对该绕流问题,流动特征在很大程度上取决于间距的大小。且间距存在一临界值,间距比从小于临界值变化到大于临界     

Flow-induced oscillation of two circular cylinders in tandem arrangement at low Re

Results are presented for flow-induced vibrations of a pair of equal-sized circular cylinders of low nondimensional mass (m^*=10) in a tandem arrangement. The cylinders are free to oscillate both in streamwise and transverse directions. The Reynolds number, based on the free-stream speed and the diameter of the cylinders, D is 100 and the centre-to-centre distance between the cylinders is 5.5D. The computations are carried out for reduced velocities in the range 2≤U^*≤15. The structural damping is set to zero for enabling maximum amplitudes of oscillation. A stabilized finite element method is utilized to carry out the computations in two dimensions. Even though the response of the upstream cylinder is found to be qualitatively similar to that of an isolated cylinder, the presence of a downstream cylinder is found to have significant effect on the behaviour of the upstream cylinder. The downstream cylinder undergoes very large amplitude of oscillations in both transverse and streamwise directions. The maximum amplitude of transverse response of the downstream cylinder is quite similar to that of a single cylinder at higher Re beyond the laminar regime. Lock-in and hysteresis are observed for both upstream and downstream cylinders. The downstream cylinder undergoes large amplitude oscillations even beyond the lock-in state. The phase between transverse oscillations and lift force suffers a 180 jump for both the cylinders almost in the middle of the synchronization regime. The phase between the transverse response of the two cylinders is also studied. Complex flow patterns are observed in the wake of the freely vibrating cylinders. Based on the phase difference and the flow patterns, the entire flow range is divided into five sub-regions.     

Flow past a sphere in density-stratified fluid

Stratified flow past a three-dimensional obstacle such as a sphere has been a long-lasting subject of geophysical, environmental and engineering fluid dynamics. In order to investigate the effect of the stratification on the near wake, in particular, the unsteady vortex formation behind a sphere, numerical simulations of stratified flows past a sphere are conducted. The time-dependent Navier–Stokes equations are solved using a three-dimensional finite element method and a modified explicit time integration scheme. Laminar flow regime is considered, and linear stratification of density is assumed under Boussinesq approximation. The effects of stratification is implemented by density transport without diffusion. The computed results include the characteristics of the near wake as well as the effects of stratification on the separation angle. Under increased stratification, the separation on the sphere is suppressed and the wake structure behind the sphere becomes planar, resembling that behind a vertical cylinder. With further increase in stratification, the wake becomes unsteady, and consists of planar vortex shedding similar to von Karman vortex streets.     

Vortex shedding and acoustic resonance of single and tandem finned cylinders

The effect of fins on vortex shedding and acoustic resonance is investigated for isolated and two tandem cylinders exposed to cross-flow in a rectangular duct. Three spacing ratios between the tandem cylinders (S/D_e=1.5, 2 and 3) are tested for a Reynolds number range from 1.6×10⁴ to 1.1×10⁵. Measurements of sound pressure as well as mean and fluctuating velocities are performed for bare and finned cylinders with three different fin densities. The effect of fins on the sound pressure generated before the onset of acoustic resonance as well as during the pre-coincidence and coincidence resonance is found to be rather complex and depends on the spacing ratio between cylinders, the fin density and the nature of the flow-sound interaction mechanism.For isolated cylinders, the fins reduce the strength of vortex shedding only slightly, but strongly attenuate the radiated sound before and during the occurrence of acoustic resonance. This suggests that the influence of the fins on correlation length is stronger than on velocity fluctuations. In contrast to isolated cylinders, the fins in the tandem cylinder case enhance the vortex shedding process at off-resonant conditions, except for the large spacing case which exhibits a reversed effect at high Reynolds numbers. Regarding the acoustic resonance of the tandem cylinders, the fins promote the onset of the coincidence resonance, but increasing the fin density drastically weakens the intensity of this resonance. The fins are also found to suppress the pre-coincidence resonance for the tandem cylinders with small spacing ratios (S/D_e=1.5, 2 and 2), but for the largest spacing case (S/D_e=3), they are found to have minor effects on the sound pressure and the lock-in range of the pre-coincidence resonance.     

Vortex-induced vibration of two elastically coupled cylinders in side-by-side arrangement

Vortex-induced vibration (VIV) of two elastically coupled circular cylinders in side-by-side arrangement is investigated numerically. The Reynolds-averaged Navier–Stokes equations are solved by the finite element method for simulating the flow and the equation of motion is solved for calculating the vibration. The mass ratio (the ratio of the mass of the cylinder to the displaced fluid mass) is 2 and the Reynolds number is 5000 in the simulations. Simulations are carried out for one symmetric configuration (referred to be Case A) and one asymmetric configuration (referred to be Case B). In both Case A and Case B, the primary response frequencies of the two cylinders are found to be the same both inside and outside the lock-in regimes. Five response regimes are found in both cases and they are the first-mode lock-in regime, the second-mode lock-in regime, the sum-frequency lock-in regime and two transition regimes. When the vibration is transiting from the first- to the second-mode lock-in regimes, the vibration of each cylinder contains both first- and the second-mode natural frequencies, and the vibrations are usually irregular. In the transition regime between the second-mode lock-in and the sum-frequency lock-in regimes, the response frequencies of both cylinders increases with an increase in the reduced velocity until they are close to the sum of the two natural frequencies. In both cases, the lower boundary reduced velocity of the total lock-in regime (the sum of the five lock-in regimes) is about 3 and the upper boundary reduced velocity is about 11 times the first-to-second-mode natural frequency ratio.     

Numerical simulations of the vortex-shedding in flows around rectangular cylinders

Flows around rectangular cylinders with a series of width-to-height ratios are calculated by means of the Improved Finite Analytic Method (IFAM) and the formation, development and shedding of vortices from the cylinders are simulated successfully. According to these results of time-dependent processes the physical phenomena in the flows are investigated in detail, and the discontinuity of Strouhal number is explained in the case of the width-to-height ratio equal to 3.0. The numerical solutions here show good agreement with the experimental results. In addition, based on several hundreds of the calculated flow patterns a moving picture is made by the computer image processing technology and recorded on a video tape, and then the vivid pictures of the physical process of vortex-shedding can be replayed later and analysed in detail. The project supported by the National Natural Science Foundation of China     

Large eddy simulations of flow interference between two unequal sized square cylinders

A uniform flow past two unequal sized square cylinders arranged in a side-by-side pattern and at a Reynolds number of 50,000 has been investigated using large eddy simulation (LES) technique. The modelling of sub-grid scales of turbulence is done using the Smagorinsky model. The effect of the transverse gap ratio (T/D) on the flow characteristics has been studied. Numerical simulations are carried out for five different transverse gap ratios (T/D), namely 1.120, 1.250, 1.375, 1.750 and 2.500. Results in terms of the aerodynamic forces, Strouhal number, mean base pressure coefficient, streamlines, vorticity, surface pressure distribution, normal and shear stresses are presented. A shift in the stagnation point for the small square cylinder from the centre of its front face towards its gap side is seen at smaller T/D ratios. The presence of a jet-like flow seen in the gap side is more pronounced at T/D = 1.12. A biased gap side flow towards the near wake of the small square cylinder is seen at smaller T/D ratios. No interference effect is seen at T/D = 2.5. The flow behaviour is similar to that of the isolated square cylinder at this gap ratio.     

LES of incompressible heat and fluid flows past a square cylinder at high Reynolds numbers

This paper presents large eddy simulation (LES) results of incompressible heat and fluid flows around a square cylinder (SC) at zero incident angle at high Reynolds numbers (Re) in the range from 1.25×10⁵ to 3.5×10⁵. LES results are obtained on the basis of swirling strength based sub-grid model, and a higher order upwind scheme developed with respect to the Taylor expansion. It was found that, for the zero incident SC wake flows at a Reynolds number in the range {Re₅ = Re/10⁵ ∈ [1.25, 3.5]}, the Strouhal number equals to 0.1079, completely independent of the Reynolds number; the coefficient of drag is around 1.835 with an uncertainty of about 1.9%, almost non-sensitive to the Re. When Re is beyond 3.0×10⁵, the time-averaged peak value of sub-grid viscosity is over 340, implying that the role of sub-grid model is crucial in some regions where vortex motion is active and vortex interaction is intense. The time–spanwise (t-z) averaged sub-grid viscosity ratio profiles and the profiles of fluctuations of the sub-grid viscosity ratio and velocity components at four locations downstream of the SC are presented. The fields of the t-z averaged sub-grid viscosity ratio, and the instantaneous fields of streamwise and spanwise vorticities are also reported and discussed. The predicted mean Nusselt number is compared with empirical correlations, revealing that swirling strength based LES has its potential in predicting natural and industrial flows.     

Unsteady flows of a generalized second grade fluid with the fractional derivative model between two parallel plates

The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced. Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus. The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates. In addition, the solutions of the shear stresses at the plates are also determined. The project supported by the National Natural Science Foundation of China (10372007, 10002003) and CNPC Innovation Fund     

The evolution of starting vortex and secondary separation on an elliptic cylinder

Flow visualization was used to investigate experimentally the evolution process from symmetrical shedding to staggered shedding of the starting vortex and the phenomenon of secondary separation on an elliptic cylinder at moderate Reynolds numbers. The vortex structure of the flow separation was studied. The temporal variation of separation angle and length of wake vortex were given. The photographs and experimental results provided basis for further investigation of the complicated feature of the starting process of unsteady separated flows around an elliptic cylinder. The project supported by the National Natural Science Foundation of China.     

Finite volume computation of incompressible turbulent flows in general co-ordinates on staggered grids

A brief review of the computation of incompressible turbulent flow in complex geometries is given. A 2D finite volume method for the calculation of turbulent flow in general curvilinear co-ordinates is described. This method is based on a staggered grid arrangement and the contravariant flux componets are chosen as primitive variables. Turbulence is modelled either by the standard k–ε model or by a k–ε model based on RNG theory. Convection is approximated with central differences for the mean flow quantities and a TVD-type MUSCL scheme for the turbulence equations. The sensitivity of the method to the grid properties is investigated. An application of this method to a complex turbulent flow is presented. The results of computations are compared with experimental data and other numerical solutions and are found to be satisfactory.     

A least-squares finite element method for time-dependent incompressible flows with thermal convection

The time-dependent Navier–Stokes equations and the energy balance equation for an incompressible, constant property fluid in the Boussinesq approximation are solved by a least-squares finite element method based on a velocity–pressure–vorticity–temperature–heat-flux ( u –P–ω–T– q ) formulation discretized by backward finite differencing in time. The discretization scheme leads to the minimization of the residual in the l²-norm for each time step. Isoparametric bilinear quadrilateral elements and reduced integration are employed. Three examples, thermally driven cavity flow at Rayleigh numbers up to 10⁶, lid-driven cavity flow at Reynolds numbers up to 10⁴ and flow over a square obstacle at Reynolds number 200, are presented to validate the method.     

Shear layer instability and acoustic interaction in solid propellant rocket motors

Segmentation of solid propellant rocket motors has been demonstrated to be a source of unpredicted and undesirable pressure and thrust oscillations. Surface discontinuities are the primary cause of these vortex-shedding-driven oscillations, which result from a strong coupling between the shear layer instability and the acoustic motion in the chamber. The analysis of an axisymmetric geometry corresponding to a {1\over 15} subscale P230 motor of the Ariane 5 rocket is numerically computed. With a suitable mesh for the viscosity value studied, the aeroacoustics in the chamber is fully described. A coupling between the hydrodynamic instability and the organ-pipe acoustic mode is clearly demonstrated. The mechanism for frequency selection is discussed. © 1997 John Wiley & Sons, Ltd.     

分块法研究圆柱绕流升阻力

使用新的分块耦合方法，分别对单圆柱和串列双圆柱绕流进行了数值模拟．对于单圆柱绕流，低Re下计算所得到的定常涡尺寸与实验非常接近．对于串列双圆柱绕流，研究分析了改变双圆柱中心间距对上下游圆柱的升阻力系数和脉动频率所产生的影响，计算结果与实验非常吻合，为进一步研究涡致振动提供了依据．     

Computation of 3D vortex flows past a flat plate at incidence through a variational approach of the full steady euler equations

A variational method for solving directly the full steady Euler equations is presented. This method is based on both Newton's linearization and a least squares formulation. The validity of the Euler model and boundary conditions to capture the vortex sheet is discussed. A finite element approximation of the groups of conservative variables is described and results are given for 3D subsonic flows as well as supersonic flows past a flat plate at high angle of attack.     

Vortex and wake-induced vibrations of a tandem arrangement of two flexible circular cylinders with near wake interference

Results showing the dynamic response of a tandem arrangement of two vertical high aspect ratio (length over diameter) and low mass ratio (mass over mass of displaced fluid) flexible cylinders vibrating at low mode number are presented in this paper. Two circular cylinder models were aligned with the flow, so the downstream or trailing cylinder was immersed in the wake of the leading one. Centre-to-centre distances from 2 to 4 diameters were studied. The models were very similar in design, with external diameters of 16 mm and a total length of 1.5 m. Reynolds numbers up to 12 000 were achieved with reduced velocities, based on the fundamental natural frequency of the downstream cylinder in still water, up to 16. The trailing model had a mass ratio of 1.8 with a combined mass-damping parameter of 0.049, whilst the corresponding figures for the leading cylinder were 1.45 and 0.043, respectively. The dynamic response of the trailing model has been analysed by studying cross-flow and in-line amplitudes, dominant frequencies and modal amplitudes. The dynamic response of the leading one is analysed by means of its cross-flow amplitudes and dominant frequencies and it is also related to the motion of the trailing cylinder by studying the synchronisation between their instantaneous cross-flow motions. Planar digital particle image velocimetry (DPIV) was used to visualise the wake. Different response regimes have been identified based on the type of oscillations exhibited by the cylinders: vortex-induced (VIV), wake-induced (WIV) or combinations of both.