Reduction in drag and vortex shedding frequency through porous sheath around a circular cylinder

A numerical study on the laminar vortex shedding and wake flow due to a porous‐wrapped solid circular cylinder has been made in this paper. The cylinder is horizontally placed, and is subjected to a uniform cross flow. The aim is to control the vortex shedding and drag force through a thin porous wrapper around a solid cylinder. The flow field is investigated for a wide range of Reynolds number in the laminar regime. The flow in the porous zone is governed by the Darcy–Brinkman–Forchheimer extended model and the Navier–Stokes equations in the fluid region. A control volume approach is adopted for computation of the governing equations along with a second‐order upwind scheme, which is used to discretize the convective terms inside the fluid region. The inclusion of a thin porous wrapper produces a significant reduction in drag and damps the oscillation compared with a solid cylinder. Dependence of Strouhal number and drag coefficient on porous layer thickness at different Reynolds number is analyzed. The dependence of Strouhal number and drag on the permeability of the medium is also examined. Copyright © 2009 John Wiley & Sons, Ltd.     

Flow about yawed,stranded cables

The development of a steady lift force on a stranded cable, which is yawed with respect to a flow, is a unique characteristic of a cable when compared to a circular cylinder. Comparisons of lift and normal drag coefficients and wake characteristics were made between stranded cable models and the cylinder. These were based upon surface pressure and hot-wire measurements and flow visualization studies conducted in a low speed wind tunnel on rigid cables and cylinders. The models were yawed to four different yaw angles and tested within the Reynolds number range of 5,000 and 50,000. Pressure profiles for the yawed cables indicated that the lift force is directed towards the side where the primary strands are more nearly aligned with the flow. The pressure profiles also indicated that the lift force is generated by asymmetric separation. The small scale irregularities associated with wires within individual strands also appeared to have an effect on the cable's lift and drag characteristics. Results show that cables have significantly different shedding characteristics and near-wake shear layer structure when compared to the circular cylinder. For the flow regime tested, the Strouhal number showed no dependence on Reynolds number nor spanwise position along the cable.List of symbols C _dn normal drag coefficient - C _l lift coefficient - C _p pressure coefficient - D actual diameter, based on circumscribing circle for the cable - f _v shedding frequency - L/D length to actual diameter ratio - ppd peak-to-peak distance, unit span - Re Reynolds number based on actual diameter - S Strouhal number, - V free stream velocity - cable angle - azimuthal angle     

Numerical Simulation of Forced Convective Heat Transfer Past a Square Diamond-Shaped Porous Cylinder

Fluid flow and heat transfer around and through a porous cylinder is an important issue in engineering applications. In this paper a numerical study is carried out for simulating the fluid flow and forced convection heat transfer around and through a square diamond-shaped porous cylinder. The flow is two-dimensional, steady, and laminar. Conservation laws of mass, momentum, and heat transport equations are applied in the clear region and Darcy–Brinkman–Forchheimer model for simulating the flow in the porous medium has been used. Equations with the relevant boundary conditions are numerically solved using a finite volume approach. In this study, Reynolds and Darcy numbers are varied within the ranges of $1<Re<45$ and $10^{-6}<Da<10^{- 2}$ , respectively. The porosity $(\varepsilon )$ is 0.5. This paper presents the effect of Reynolds and Darcy numbers on the flow structure and heat transfer characteristics. Finally, these parameters are compared among solid and porous cylinder. It was found that the drag coefficient decreases and flow separation from the cylinder is delayed with increasing Darcy number. Also the size of the thermal plume decreases by decreasing Darcy number.     

Effect of magnetic Reynolds number on the two‐dimensional hydromagnetic flow around a cylinder

Numerical experiments have been conducted to study the effect of magnetic Reynolds number on the steady, two‐dimensional, viscous, incompressible and electrically conducting flow around a circular cylinder. Besides usual Reynolds number Re, the flow is governed by the magnetic Reynolds number R_m and Alfvén number β. The flow and magnetic field are uniform and parallel at large distances from the cylinder. The pressure Poisson equation is solved to find the pressure fields in the entire flow region. The effects of the magnetic field and electrical conductivity on the recirculation bubble, drag coefficient, standing vortex and pressure are presented and discussed. For low interaction parameter (N<1), the suppression of the flow‐separation is nearly independent of the conductivity of the fluid, whereas for large interaction parameters, the conductivity of the fluid strongly influences the control of flow‐separation. Copyright © 2008 John Wiley & Sons, Ltd.     

Fluid forces on a square cylinder in oscillating flows with non‐zero‐mean velocities

The unsteady forces on a square cylinder in sinusoidally oscillating flows with non‐zero‐mean velocities are investigated numerically by using a weakly compressible‐flow method with three‐dimensional large eddy simulations. The major parameters in the analysis are Keulegan–Carpenter number (KC) and the ratio between the amplitude and the mean velocities of the approaching flow (A_R). By varying the values of KC and A_R the resulting drag and lift of the cylinders are analyzed systematically at two selected approaching‐flow attack angles (0 and 22.5^°). In the case of the non‐zero attack angle, results show that both the drag and lift histories can be adequately described by Morison equations. However, Morison equations fail to correctly describing the lift history as the attack angle is zero. In addition, when the ratio of A_R/KC is near the Strouhal number of the bluff‐body flow, the resulting drag is promoted due to the occurrence of resonance. Based on the results of systematic analyses, finally, the mean and inertia force coefficients at the two selected attack angles are presented as functions of KC and A_R based on the Morison relationships. Copyright © 2008 John Wiley & Sons, Ltd.     

Flow past a square cylinder at low Reynolds numbers

Results are presented for the flow past a stationary square cylinder at zero incidence for Reynolds number, Re ? 150. A stabilized finite‐element formulation is employed to discretize the equations of incompressible fluid flow in two‐dimensions. For the first time, values of the laminar separation Reynolds number, Re_s, and separation angle, θ_s, at Re_s are predicted. Also, the variation of θ_s with Re is presented. It is found that the steady separation initiates at Re = 1.15. Contrary to the popular belief that separation originates at the rear sharp corners, it is found to originate from the base point, i.e. θ_s=180^° at Re = Re_s. For Re > 5, θ_s approaches the limit of 135 ^°. The length of the separation bubble increases approximately linearly with increasing Re. The drag coefficient varies as Re^?0.66. Flow characteristics at Re ? 40 are also presented for elliptical cylinders of aspect ratios 0.2, 0.5, 0.8 and 1 (circle) having the same characteristic dimension as the square and major axis oriented normal to the free‐stream. Compared with a circular cylinder, the flow separates at a much lower Re from a square cylinder leading to the formation of a bigger wake (larger bubble length and width). Consequently, at a given Re, the drag on a square cylinder is more than the drag of a circular cylinder. This suggests that a cylinder with square section is more bluff than the one with circular section. Among all the cylinder shapes studied, the square cylinder with sharp corners generates the largest amount of drag. Copyright © 2010 John Wiley & Sons, Ltd.     

Passive control of wake flow by two small control cylinders at Reynolds number 80

Passive control of the wake behind a circular cylinder in uniform flow is studied by numerical simulation at Re_D=80. Two small control cylinders are placed symmetrically along the separating shear layers at various stream locations. In the present study, the detailed flow mechanisms that lead to a significant reduction in the fluctuating lift but maintain the shedding vortex street are clearly revealed. When the stream locations lie within 0.8≤X_C/D≤3.0, the alternate shedding vortex street remains behind the control cylinders. In this case, the symmetric standing eddies immediately behind the main cylinder and the downstream delay of the shedding vortex street are the two primary mechanisms that lead to a 70–80% reduction of the fluctuating lift on the main cylinder. Furthermore, the total drag of all the cylinders still has a maximum 5% reduction. This benefit is primarily attributed to the significant reduction of the pressure drag on the main cylinder. Within X_C/D>3.0, the symmetry of the standing eddy breaks down and the staggered vortex street is similar to that behind a single cylinder at the same Reynolds number. In the latter case, the mean pressure drag and the fluctuating lift coefficients on the main cylinder will recover to the values of a single cylinder.     

Drag and Separation Flow Past Solid Sphere with Porous Shell at Moderate Reynolds Number

Axisymmetric viscous, two-dimensional steady and incompressible fluid flow past a solid sphere with porous shell at moderate Reynolds numbers is investigated numerically. There are two dimensionless parameters that govern the flow in this study: the Reynolds number based on the free stream fluid velocity and the diameter of the solid core, and the ratio of the porous shell thickness to the square root of its permeability. The flow in the free fluid region outside the shell is governed by the Navier–Stokes equation. The flow within the porous annulus region of the shell is governed by a Darcy model. Using a commercially available computational fluid dynamics (CFD) package, drag coefficient and separation angle have been computed for flow past a solid sphere with a porous shell for Reynolds numbers of 50, 100, and 200, and for porous parameter in the range of 0.025–2.5. In all simulation cases, the ratio of b/a was fixed at 1.5; i.e., the ratio of outer shell radius to the inner core radius. A parametric equation relating the drag coefficient and separation point with the Reynolds number and porosity parameter were obtained by multiple linear regression. In the limit of very high permeability, the computed drag coefficient as well as the separation angle approaches that for a solid sphere of radius a, as expected. In the limit of very low permeability, the computed total drag coefficient approaches that for a solid sphere of radius b, as expected. The simulation results are presented in terms of viscous drag coefficient, separation angles and total drag coefficient. It was found that the total drag coefficient around the solid sphere as well as the separation angle are strongly governed by the porous shell permeability as well as the Reynolds number. The separation point shifts toward the rear stagnation point as the shell permeability is increased. Separation angle and drag coefficient for the special case of a solid sphere of radius r = a was found to be in good agreement with previous experimental results and with the standard drag curve.     

带有柔性薄板三维方柱流固耦合的数值模拟研究

采用双向流固耦合方法,对带有柔性薄板三维方柱的流场变化特性进行了研究。通过对比单方柱,分析了带有柔性薄板三维方柱阻力系数、升力系数以及斯特劳哈尔数的变化规律。研究表明,在方柱尾流区域附加一柔性薄板可以使其阻力系数降低34.6%,同时其变化幅值大大减小;其升力系数的均方根减小84.8%,流场脉动大幅度减小;斯特劳哈尔数降低79.5%。研究结果表明,在三维方柱后设置柔性薄板可有效抑制涡脱落,从而改善三维方柱的尾流特性。     

Reference values for drag and lift of a two‐dimensional time‐dependent flow around a cylinder

This paper presents a numerical study of a two‐dimensional time‐dependent flow around a cylinder. Its main objective is to provide accurate reference values for the maximal drag and lift coefficient at the cylinder and for the pressure difference between the front and the back of the cylinder at the final time. In addition, the accuracy of these values obtained with different time stepping schemes and different finite element methods is studied. Copyright © 2004 John Wiley & Sons, Ltd.     

Experimental study of sedimentation characteristics of spheroidal particles

The drag of non-spherical particles is a basic, important parameter for multi-phase flow. As the first step in research in this area, the terminal velocities, Ut, of hemispherical and spherical segment particles with maximal diameters of 6-21 mm were measured in static fluids by using a high-speed video camera. The drag coefficient, CD, measured for Reynolds number, Re of 10^1-10^5, has been obtained and compared with those for a sphere. The Re based on the terminal velocity has a logarithmic linear relationship with Ar number for both the facet facing upwards or downwards for the two experimental spheroidal particles, and their Co values are greater than those of spheres. A shape function that depends on the initial orientation of the particle facet is presented to correct for the shape effects.     

Effects of aspect ratio on the drag of a wall-mounted finite-length cylinder in subcritical and critical regimes

The effect of cylinder aspect ratio (??H/d, where H is the cylinder height or length, and d is the cylinder diameter) on the drag of a wall-mounted finite-length circular cylinder in both subcritical and critical regimes is experimentally investigated. Two cases are considered: a smooth cylinder submerged in a turbulent boundary layer and a roughened cylinder immersed in a laminar uniform flow. In the former case, the Reynolds number Re _d (??U _?? d/??, with U _?? being the free-stream velocity and ?? the fluid viscosity) was varied from 2.61?×?10⁴ to 2.87?×?10⁵, and two values of H/d (2.65 and 5) were examined; in the latter case, Re _d ?=?1.24?×?10⁴?C1.73?×?10⁵ and H/d?=?3, 5 and 7. In the subcritical regime, both the drag coefficient C _D and the Strouhal number St are smaller than their counterparts for a two-dimensional cylinder and reduce monotonously with decreasing H/d. The presence of a turbulent boundary layer causes an early transition from the subcritical to critical regime and considerably enlarges the Re _d range of the critical regime. No laminar separation bubble occurs on the finite-length cylinder immersed in the turbulent boundary layer, and consequently, the discontinuity is not observed in the C _D?CRe _d and St?CRe _d curves. In the roughened cylinder case, the Re _d range of the critical regime grows gradually with decreasing H/d, while the C _D crisis becomes less obvious. In both cases, H/d has a negligible effect on the critical value of Re _d at which transition occurs from the subcritical to critical regime.     

Aspect ratio effect on flow-induced forces on circular cylinders in a cross-flow

Finite-span circular cylinders with two different aspect ratios, placed in a cross-flow, are investigated experimentally at a cylinder Reynolds number of 46,000. Simultaneous measurements of the flow-induced unsteady forces on the cylinders and the stream velocity in the wake are carried out. These results together with mean drag measurements along the span and available literature data are used to evaluate the flow mechanisms responsible for the induced unsteady forces and the effect of aspect ratio on these forces. The coherence of vortex shedding along the span of the cylinder is partially destroyed by the separated flow emanating from the top and by the recirculating flow behind the cylinder. As a result, the fluctuating lift decreases drastically. Based on the data collected, it is conjectured that the fluctuating recirculating flow behind the cylinder is the flow mechanism responsible for the unsteady drag and causes it to increase beyond the fluctuating lift. The fluctuating recirculating flow is a direct consequence of the unsteady separated flow. The unsteady forces vary along the span, with lift increasing and drag decreasing towards the cylinder base. When the cylinder span is large compared to the wall boundary layer thickness, a submerged two-dimensional region exists near the base. As the span decreases, the submerged two-dimensional region becomes smaller and eventually vanishes. Altogether, these results show that fluctuating drag is the dominant unsteady force in finite-span cylinders placed in a cross-flow. Its characteristic frequency is larger than that of the vortex shedding frequency.List of symbols a span of active element on cylinder, = 2.5 cm - C _D local rms drag coefficient, 2D/ U ² da - C _L local rms lift coefficient, 2l/ U ² da - C _D local mean drag coefficient, 2D/ U ² da - C _D spanwise-averaged C _D for finite-span cylinder - (C _D ) _2D spanwise-averaged mean drag coefficient for two-dimensional cylinder - C _p pressured coefficient - -(C _p ) _b pressure coefficient at = - d diameter of cylinder, = 10.2 cm - D fluctuating component of instantaneous drag - D local rms of fluctuating drag - D local mean drag - E _D power spectrum of fluctuating drag, defined as - E _L power spectra of fluctuating lift, defined as - f _D dominant frequency of drag spectrum - f _L dominant frequency of lift spectrum - f _u dominant frequency of velocity spectrum - h span of cylinder - H height of test section, = 30.5 cm - L fluctuating component of instantaneous lift - L local rms of fluctuating lift - R _Du () cross-correlation function of streamwise velocity and local drag, - R _Lu () cross-correlation function of stream wise velocity and local lift, - Re Reynolds number, U d/y - S _L Strouhal number based on f _L ,f _L d/U - S _D Strouhal number based on f _D ,f _D d/U - S _u Strouhal number based on f _u , f _u d/U - t time - u fluctuating component of instantaneous streamwise velocity - U mean streamwise velocity - mean stream velocity upstream of cylinder - x streamwise distance measured from axis of cylinder - y transverse distance measured from axis of test section - z spanwise distance measured from cylinder base - angular position on cylinder circumference measured from forward stagnation - kinematic viscosity of air - density of air - time lag in cross-correlation function - _D normalized spectrum of fluctuating drag - _L normalized spectrum of fluctuating lift     

Unsteady forces on circular cylinders in a cross-flow

A three-axis piezoelectric load cell was used to measure the local unsteady forces induced on cylinders placed in a cross-flow. In conjunction with this, a single hot-wire was used to traverse the wake at a fixed distance behind the cylinder so that correlations between the induced forces on the cylinder and the wake velocity could be calculated to provide insight into the character of the flow-induced unsteady forces. Experiments were carried out on both two-dimensional and finite-span cylinders at a Reynolds number of 46,000. For the two-dimensional cylinder case, substantial evidence was obtained to demonstrate that the strength of the vortex roll-up along the span was quite uniform. Consequently, the lift-velocity correlation along the span remained unchanged. On the other hand, there was a total lack of correlation between the fluctuating drag and the wake velocity, thus indicating that the drag signal was not quite periodic. In the finite-span cylinder case, the separated flow from the top edge of the cylinder was found to suppress vortex shedding along the span of the cylinder, destroyed its coherence and caused the wake flow to oscillate in the stream direction. This oscillation induced a significant fluctuating drag on the cylinder. Consequently, the fluctuating drag far exceeded the fluctuating lift and the wake velocity was found to correlate well with the drag and not with the lift. This correlation remained intact along the span of the cylinder. Finally, the rms fluctuating lift and drag forces were found to vary along the cylinder span, with the lift increasing and the drag decreasing as the base of the cylinder is approached; thus suggesting that a submerged two-dimensional region exists near the base of the cylinder.List of symbols a span of active element on cylinder - C _D local rms drag coefficient, - C _L local rms lift coefficient, - C _D local mean drag coefficient - (C _D ) _2D spanwise-averaged mean drag coefficient for two dimensional cylinder - d diameter of cylinder (= 10.2 cm) - D fluctuating component of instantaneous drag - D local rms of fluctuating drag - E _D power spectrum of fluctuating drag, defined as - E _L power spectrum of fluctuating lift, defined as - E _U power spectrum of fluctuating streamwise velocity, defined as - f _L dominant frequency of lift spectrum - f _D dominant frequency of drag spectrum - f _u dominant frequency of velocity spectrum - h span of cylinder - H height of test section (= 30.5 cm) - L fluctuating component of instantaneous lift - L local rms of fluctuating lift - R _Du () cross-correlation function of streamwise velocity and local drag - R _Lu () cross-correlation function of streamwise velocity and local lift - Re Reynolds number, - S _L Strouhal number based on f _L , - S _D Strouhal number based on f _D , - S _U Strouhal number based on f _u , - t time - u fluctuating component of instantaneous streamwise velocity - u rms of streamwise fluctuating velocity - u rms of streamwise fluctuating velocity upstream of cylinder - U mean streamwise velocity - U mean stream velocity upstream of cylinder - x streamwise distance measured from axis of cylinder - y transverse distance measured from axis of cylinder - z spanwise distance measured from floor of test section - v kinematic viscosity of air - density of air - time lag in cross-correlation function - _D normalized spectrum of fluctuating drag - _L normalized spectrum of fluctuating lift - _U normalized spectrum of fluctuating streamwise velocity     

Flow past a cylinder: shear layer instability and drag crisis

Flow past a circular cylinder for Re=100 to 10⁷ is studied numerically by solving the unsteady incompressible two‐dimensional Navier–Stokes equations via a stabilized finite element formulation. It is well known that beyond Re ～ 200 the flow develops significant three‐dimensional features. Therefore, two‐dimensional computations are expected to fall well short of predicting the flow accurately at high Re. It is fairly well accepted that the shear layer instability is primarily a two‐dimensional phenomenon. The frequency of the shear layer vortices, from the present computations, agree quite well with the Re^0.67 variation observed by other researchers from experimental measurements. The main objective of this paper is to investigate a possible relationship between the drag crisis (sudden loss of drag at Re ～ 2 × 10⁵) and the instability of the separated shear layer. As Re is increased the transition point of shear layer, beyond which it is unstable, moves upstream. At the critical Reynolds number the transition point is located very close to the point of flow separation. As a result, the shear layer eddies cause mixing of the flow in the boundary layer. This energizes the boundary layer and leads to its reattachment. The delay in flow separation is associated with narrowing of wake, increase in Reynolds shear stress near the shoulder of the cylinder and a significant reduction in the drag and base suction coefficients. The spatial and temporal power spectra for the kinetic energy of the Re=10⁶ flow are computed. As in two‐dimensional isotropic turbulence, E(k) varies as k^?5/3 for wavenumbers higher than energy injection scale and as k^?3 for lower wavenumbers. The present computations suggest that the shear layer vortices play a major role in the transition of boundary layer from laminar to turbulent state. Copyright © 2004 John Wiley & Sons, Ltd.     

Effect of Thermal Non-Equilibrium on Convective Instability in a Ferromagnetic Fluid-Saturated Porous Medium

The effect of local thermal non-equilibrium (LTNE) on the onset of thermomagnetic convection in a ferromagnetic fluid-saturated horizontal porous layer in the presence of a uniform vertical magnetic field is investigated. A modified Forchheimer-extended Darcy equation is employed to describe the flow in the porous medium, and a two-field model is used for temperature representing the solid and fluid phases separately. It is found that both the critical Darcy–Rayleigh number and the corresponding wave number are modified by the LTNE effects. Asymptotic solutions for both small and large values of scaled interphase heat transfer coefficient H _t are presented and compared with those computed numerically. An excellent agreement is obtained between the asymptotic and the numerical results. Besides, the influence of magnetic parameters on the instability of the system is also discussed. The available results in the literature are recovered as particular cases from the present study.     

Application of the Skeleton Model of a Highly Porous Cellular Material in Modeling Supersonic Flow past a Cylinder with a Forward Gas-Permeable Insert

The supersonic (M_∞ = 4.85) flow past a cylinder with a forward insertmade of a highlyporous cellular material is numerically modeled within the framework of the Reynolds-averaged Navier–Stokes equations. The air flow in the gas-permeable insert is described on the basis of a skeleton model of a highly-porous medium, whose determining parameters are the porosity coefficient (95%) and the pore dimensions (1 mm) of the actual cellular material. The aerodynamic drag coefficients of the model with different lengths of the porous forward insert are calculated on the unit Reynolds number range from 6.9 × 10⁵ to 13.8 × 10⁶ m^?1. They are in agreement with the available experimental data, which indicates the adequacy of the proposed skeleton model in describing the actual properties of highly-porous materials.     

Interference effect of two equal‐sized square cylinders in tandem arrangement: With planar shear flow

Numerical simulations have been performed for flow past two equal‐sized square cylinders in tandem arrangement subjected to incoming planar shear flow. Effect of L/d ratio and the shear parameter has been studied. The range of L/d ratio (ratio of center‐to‐center distance (L) to cylinder width (d)) is varied from 2 to 7 and the non‐dimensional shear parameter (K) is varied from 0.0 to 0.4 in steps of 0.1. For all the cases the Reynolds number (Re) based on centerline velocity and cylinder width is fixed at 100. The results are compared with that of isolated square cylinder with uniform flow. Strouhal number decreases with increasing shear parameter. There are more than one shedding frequency at high shear parameters and L/d ratios. The mean drag coefficient is decreased with shear parameter and lesser than that of the single cylinder. The root mean square (RMS) value of both lift and drag coefficients is higher for the downstream cylinder for all values of shear parameter. With increasing L/d ratio, for both lift and drag, the RMS value increases and then decreases for upstream cylinder, whereas it continuously increases for the downstream cylinder. The stagnation point is moved towards the top leading edge with increasing shear. The critical L/d ratio, which is defined as the distance between two cylinders, beyond which the vortex shedding from the upstream cylinder occurs, decreases with increasing shear parameter. Copyright © 2007 John Wiley & Sons, Ltd.     

The aerodynamics of a cylinder submerged in the wake of another

Flow-induced fluctuating lift (C_Lf) and drag (C_Df) forces and Strouhal numbers (St) of a cylinder submerged in the wake of another cylinder are investigated experimentally for Reynolds number (Re)=9.7×10³–6.5×10⁴. The spacing ratio L^⁎ (=L/D) between the cylinders is varied from 1.1 to 4.5, where L is the spacing between the cylinders and D is the cylinder diameter. The results show that C_Lf, C_Df and St are highly sensitive to Re due to change in the inherent nature of the flow structure. How the flow structure is dependent on Re and L^⁎ is presented in a flow structure map. Zdravkovich and Pridden (1977) observed a ‘kink’ in time-mean drag distribution at L^⁎≈2.5 for Re>3.1×10⁴, but not for Re≤3.1×10⁴. The physics is provided here behind the presence and absence of the ‘kink’ that was left unexplained since then.     

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

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