Two‐dimensional prediction of time dependent,turbulent flow around a square cylinder confined in a channel

This paper presents two‐dimensional and unsteady RANS computations of time dependent, periodic, turbulent flow around a square block. Two turbulence models are used: the Launder–Sharma low‐Reynolds number k–ε model and a non‐linear extension sensitive to the anisotropy of turbulence. The Reynolds number based on the free stream velocity and obstacle side is Re=2.2×10⁴. The present numerical results have been obtained using a finite volume code that solves the governing equations in a vertical plane, located at the lateral mid‐point of the channel. The pressure field is obtained with the SIMPLE algorithm. A bounded version of the third‐order QUICK scheme is used for the convective terms. Comparisons of the numerical results with the experimental data indicate that a preliminary steady solution of the governing equations using the linear k–ε does not lead to correct flow field predictions in the wake region downstream of the square cylinder. Consequently, the time derivatives of dependent variables are included in the transport equations and are discretized using the second‐order Crank–Nicolson scheme. The unsteady computations using the linear and non‐linear k–ε models significantly improve the velocity field predictions. However, the linear k–ε shows a number of predictive deficiencies, even in unsteady flow computations, especially in the prediction of the turbulence field. The introduction of a non‐linear k–ε model brings the two‐dimensional unsteady predictions of the time‐averaged velocity and turbulence fields and also the predicted values of the global parameters such as the Strouhal number and the drag coefficient to close agreement with the data. Copyright © 2009 John Wiley & Sons, Ltd.     

High-fidelity simulations of moving and flexible airfoils at low Reynolds numbers

The present paper highlights results derived from the application of a high-fidelity simulation technique to the analysis of low-Reynolds-number transitional flows over moving and flexible canonical configurations motivated by small natural and man-made flyers. This effort addresses three separate fluid dynamic phenomena relevant to small fliers, including: laminar separation and transition over a stationary airfoil, transition effects on the dynamic stall vortex generated by a plunging airfoil, and the effect of flexibility on the flow structure above a membrane airfoil. The specific cases were also selected to permit comparison with available experimental measurements. First, the process of transition on a stationary SD7003 airfoil section over a range of Reynolds numbers and angles of attack is considered. Prior to stall, the flow exhibits a separated shear layer which rolls up into spanwise vortices. These vortices subsequently undergo spanwise instabilities, and ultimately breakdown into fine-scale turbulent structures as the boundary layer reattaches to the airfoil surface. In a time-averaged sense, the flow displays a closed laminar separation bubble which moves upstream and contracts in size with increasing angle of attack for a fixed Reynolds number. For a fixed angle of attack, as the Reynolds number decreases, the laminar separation bubble grows in vertical extent producing a significant increase in drag. For the lowest Reynolds number considered (Re _c  = 10⁴), transition does not occur over the airfoil at moderate angles of attack prior to stall. Next, the impact of a prescribed high-frequency small-amplitude plunging motion on the transitional flow over the SD7003 airfoil is investigated. The motion-induced high angle of attack results in unsteady separation in the leading edge and in the formation of dynamic-stall-like vortices which convect downstream close to the airfoil. At the lowest value of Reynolds number (Re _c  = 10⁴), transition effects are observed to be minor and the dynamic stall vortex system remains fairly coherent. For Re _c  = 4 × 10⁴, the dynamic-stall vortex system is laminar at is inception, however shortly afterwards, it experiences an abrupt breakdown associated with the onset of spanwise instability effects. The computed phased-averaged structures for both values of Reynolds number are found to be in good agreement with the experimental data. Finally, the effect of structural compliance on the unsteady flow past a membrane airfoil is investigated. The membrane deformation results in mean camber and large fluctuations which improve aerodynamic performance. Larger values of lift and a delay in stall are achieved relative to a rigid airfoil configuration. For Re _c = 4.85 × 10⁴, it is shown that correct prediction of the transitional process is critical to capturing the proper membrane structural response.     

Transient growth in the asymptotic suction boundary layer

In the present experimental setup, the transient disturbance growth in a spatially invariant boundary layer flow, i.e., the asymptotic suction boundary layer (ASBL), has been investigated. The choice of the ASBL brings along several advantages compared with an ordinary spatially growing boundary layer. A unique feature of the ASBL is that the Reynolds number (Re) can be varied without changing the boundary layer thickness, which in turn allows for parameter variations not possible to carry out in traditional boundary layer flows. A spanwise array of discrete surface roughness elements was mounted on the surface to trigger modes with different spanwise wavenumbers (β). It is concluded that for each mode there exists a threshold roughness Reynolds number (Re _k ), below which no significant transient growth is present. The experimental data suggests that this threshold Re _k is both a function of β and Re. An interesting result is that the energy growth curves respond differently to a change in Re _k when caused by a change in roughness height k, implying that Re remains constant, compared with a change in the free-stream velocity U_￥U_\infty, which also affects the Re. The scaling of the energy growth curves both in level and the downstream direction is treated and appropriate scalings are found. The result shows a complex non-linear receptivity mechanism. Optimal perturbation theory, which has failed to predict the energy evolution in growing boundary layers, is tested for the ASBL and shows that it may satisfactorily predict the evolution of all transiently growing modes that are triggered by the roughness elements.     

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.     

Time‐dependent turbulent cavitating flow computations with interfacial transport and filter‐based models

Turbulent cavitating flow computations need to address both cavitation and turbulence modelling issues. A recently developed interfacial dynamics‐based cavitation model (IDCM) incorporates the interfacial transport into the computational modelling of cavitation dynamics. For time‐dependent flows, it is known that the engineering turbulence closure such as the original k–ε model often over‐predicts the eddy viscosity values reducing the unsteadiness. A recently proposed filter‐based modification has shown that it can effectively modulate the eddy viscosity, rendering better simulation capabilities for time‐dependent flow computations in term of the unsteady characteristics. In the present study, the IDCM along with the filter‐based k–ε turbulence model is adopted to simulate 2‐D cavitating flows over the Clark‐Y airfoil. The chord Reynolds number is Re=7.0 × 10⁵. Two angles‐of‐attack of 5 and 8° associated with several cavitation numbers covering different flow regimes are conducted. The simulation results are assessed with the experimental data including lift, drag and velocity profiles. The interplay between cavitation and turbulence models reveals substantial differences in time‐dependent flow results even though the time‐averaged characteristics are similar. Copyright © 2005 John Wiley & Sons, Ltd.     

Non-unique turbulent boundary layer flows having a moderately large velocity defect: a rational extension of the classical asymptotic theory

The classical analysis of turbulent boundary layers in the limit of large Reynolds number Re is characterised by an asymptotically small velocity defect with respect to the external irrotational flow. As an extension of the classical theory, it is shown in the present work that the defect may become moderately large and, in the most general case, independent of Re but still remain small compared to the external streamwise velocity for non-zero pressure gradient boundary layers. That wake-type flow turns out to be characterised by large values of the Rotta–Clauser parameter, serving as an appropriate measure for the defect and hence as a second perturbation parameter besides Re. Most important, it is demonstrated that also this case can be addressed by rigorous asymptotic analysis, which is essentially independent of the choice of a specific Reynolds stress closure. As a salient result of this procedure, transition from the classical small defect to a pronounced wake flow is found to be accompanied by quasi-equilibrium flow, described by a distinguished limit that involves the wall shear stress. This situation is associated with double-valued solutions of the boundary layer equations and an unconventional weak Re-dependence of the external bulk flow—a phenomenon seen to agree well with previous semi-empirical studies and early experimental observations. Numerical computations of the boundary layer flow for various values of Re reproduce these analytical findings with satisfactory agreement.     

LES and RANS for Turbulent Flow over Arrays of Wall-Mounted Obstacles

Large-eddy simulation (LES) has been applied to calculate the turbulent flow over staggered wall-mounted cubes and staggered random arrays of obstacles with area density 25%, at Reynolds numbers between 5 × 10³ and 5 10⁶, based on the free stream velocity and the obstacle height. Re = 5 × 10³ data were intensively validated against direct numerical simulation (DNS) results at the same Re and experimental data obtained in a boundary layer developing over an identical roughness and at a rather higher Re. The results collectively confirm that Reynolds number dependency is very weak, principally because the surface drag is predominantly form drag and the turbulence production process is at scales comparable to the roughness element sizes. LES is thus able to simulate turbulent flow over the urban-like obstacles at high Re with grids that would be far too coarse for adequate computation of corresponding smooth-wall flows. Comparison between LES and steady Reynolds-averaged Navier-Stokes (RANS) results are included, emphasising that the latter are inadequate, especially within the canopy region.     

Transition from vortex to wall driven turbulence production in the Taylor–Couette system with a rotating inner cylinder

Axisymmetrically stable turbulent Taylor vortices between two concentric cylinders are studied with respect to the transition from vortex to wall driven turbulent production. The outer cylinder is stationary and the inner cylinder rotates. A low Reynolds number turbulence model using the k‐ω formulation, facilitates an analysis of the velocity gradients in the Taylor–Couette flow. For a fixed inner radius, three radius ratios 0.734, 0.941 and 0.985 are employed to identify the Reynolds number range at which this transition occurs. At relatively low Reynolds numbers, turbulent production is shown to be dominated by the outflowing boundary of the Taylor vortex. As the Reynolds number increases, shear driven turbulence (due to the rotating cylinder) becomes the dominating factor. For relatively small gaps turbulent flow is shown to occur at Taylor numbers lower than previously reported. Copyright © 2002 John Wiley & Sons, Ltd.     

Transitionally rough zero pressure gradient turbulent boundary layers

Near-wall measurements are performed to study the effects of surface roughness and viscous shear stresses on the transitionally rough regime (5 < k ⁺ < 70) of a zero pressure gradient turbulent boundary layer. The x-dependence is known from the eleven consecutive measurements in the streamwise direction, which allows for the computation of the streamwise gradients in the boundary layer equations. Thus, the skin friction is computed from the integrated boundary layer equation with errors of 3 and 5% for smooth and rough, respectively. It is found that roughness destroys the viscous layer near the wall, thus, reducing the contribution of the viscous stress in the wall region. As a result, the contribution in the wall shear stress due to form drag increases, while the viscous stress decreases. This yields Reynolds number invariance in the skin friction as k ⁺ increases into the fully rough regime. Furthermore, the roughness at the wall reduces the high peak of the streamwise component of the Reynolds stress in the near-wall region. However, for the Reynolds wall-normal and shear stress components, its contribution is not significantly altered for sand grain roughness.     

Computation of two‐dimensional flows past ram‐air parachutes

Computational results for flow past a two‐dimensional model of a ram‐air parachute with leading edge cut are presented. Both laminar (Re=10⁴) and turbulent (Re=10⁶) flows are computed. A well‐proven stabilized finite element method (FEM), which has been applied to various flow problems earlier, is utilized to solve the incompressible Navier–Stokes equations in the primitive variables formulation. The Baldwin–Lomax model is employed for turbulence closure. Turbulent flow computations past a Clarck‐Y airfoil without a leading edge cut, for α=7.5°, result in an attached flow. The leading edge cut causes the flow to become unsteady and leads to a significant loss in lift and an increase in drag. The flow inside the parafoil cell remains almost stagnant, resulting in a high value of pressure, which is responsible for giving the parafoil its shape. The value of the lift‐to‐drag ratio obtained with the present computations is in good agreement with those reported in the literature. The effect of the size and location of the leading edge cut is studied. It is found that the flow on the upper surface of the parafoil is fairly insensitive to the configuration of the cut. However, the flow quality on the lower surface improves as the leading edge cut becomes smaller. The lift‐to‐drag ratio for various configurations of the leading edge cut varies between 3.4 and 5.8. It is observed that even though the time histories of the aerodynamic coefficients from the laminar and turbulent flow computations are quite different, their time‐averaged values are quite similar. Copyright © 2001 John Wiley & Sons, Ltd.     

Effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics of a sphere

Experiments were conducted in water and wind tunnels on spheres in the Reynolds number range 6 × 10³ to 6.5 × 10⁵ to study the effect of natural ventilation on the boundary layer separation and near-wake vortex shedding characteristics. In the subcritical range of Re (<2 × 10⁵), ventilation caused a marginal downstream shift in the location of laminar boundary layer separation; there was only a small change in the vortex shedding frequency. In the supercritical range (Re > 4 × 10⁵), ventilation caused a downstream shift in the mean locations of boundary layer separation and reattachment; these lines showed significant axisymmetry in the presence of venting. No distinct vortex shedding frequency was found. Instead, a dramatic reduction occurred in the wake unsteadiness at all frequencies. The reduction of wake unsteadiness is consistent with the reduction in total drag already reported. Based on the present results and those reported earlier, the effects of natural ventilation on the flow past a sphere can be categorized in two broad regimes, viz., weak and strong interaction regimes. In the weak interaction regime (subcritical Re), the broad features of the basic sphere are largely unaltered despite the large addition of mass in the near wake. Strong interaction is promoted by the closer proximity of the inner and outer shear layers at supercritical Re. This results in a modified and steady near-wake flow, characterized by reduced unsteadiness and small drag. Received: 8 September 1998 / Accepted: 1 January 2000     

LES-based Study of the Roughness Effects on the Wake of a Circular Cylinder from Subcritical to Transcritical Reynolds Numbers

This paper investigates the effects of surface roughness on the flow past a circular cylinder at subcritical to transcritical Reynolds numbers. Large eddy simulations of the flow for sand grain roughness of size k/D = 0.02 are performed (D is the cylinder diameter). Results show that surface roughness triggers the transition to turbulence in the boundary layer at all Reynolds numbers, thus leading to an early separation caused by the increased momentum deficit, especially at transcritical Reynolds numbers. Even at subcritical Reynolds numbers, boundary layer instabilities are triggered in the roughness sublayer and eventually lead to the transition to turbulence. The early separation at transcritical Reynolds numbers leads to a wake topology similar to that of the subcritical regime, resulting in an increased drag coefficient and lower Strouhal number. Turbulent statistics in the wake are also affected by roughness; the Reynolds stresses are larger due to the increased turbulent kinetic energy production in the boundary layer and separated shear layers close to the cylinder shoulders.     

Vortex shedding in cylinder flow of shear-thinning fluids: I. Identification and demarcation of flow regimes

An experimental study on the flow of non-Newtonian fluids around a cylinder was undertaken to identify and delimit the various shedding flow regimes as a function of adequate non-dimensional numbers. The measurements of vortex shedding frequency and formation length (l_f) were carried out by laser-Doppler anemometry in Newtonian fluids and in aqueous polymer solutions of CMC and tylose. These were shear thinning and elastic at weight concentrations ranging from 0.1 to 0.6%. The 10 and 20 mm diameter cylinders (D) used in the experiments had aspect ratios of 12 and 6 and blockage ratios of 5 and 10%, respectively. The Reynolds number (Re^*) was based on a characteristic shear rate of U_∞/(2D) and ranged from 50 to 9×10³ thus encompassing the laminar shedding, the transition and shear-layer transition regimes. Increasing fluid elasticity reduced the various critical Reynolds numbers (Re_etr^*, Re_lf^*, Re_bbp^*) and narrowed the extent of the transition regime. For the 0.6% tylose solution the transition regime was even suppressed. On the other end, pseudoplasticity was found to be indirectly responsible for the observed reduction in Re_otr^*: it increases the Strouhal number which in turn increases the vortex filaments, precursors of the transition regime. Elasticity was better quantified by the elasticity number Re′/We than by the Weissenberg number. This elasticity number involves the calculation of the viscosity at a high characteristic shear rate, typical of the boundary layer, rather than at the average value (U_∞/(2D)) used for the Reynolds number, Re^*.     

Reynolds‐stress modelling of M = 2.25 shock‐wave/turbulent boundary‐layer interaction

M = 2.25 shock‐wave/turbulent‐boundary‐layer interactions over a compression ramp for several angles (8, 13 and 18°) at Reynolds‐number Re=7 × 10³ were simulated with three low‐Reynolds second‐moment closures and a linear low‐Reynolds standard k–ε model. A detailed assessment of the turbulence closures by comparison with both mean‐flow and turbulent experimental quantities is presented. The Reynolds‐stress model which is wall‐topology free and which uses an optimized redistribution closure, is in good agreement with experimental data both for wall‐pressure and mean‐velocity profiles. Detailed analysis of three components of the Reynolds‐stress tensor (comparison with measurements and transport‐equation budgets) provides a critical evaluation of full Reynolds‐stress models for the separated supersonic compression ramp. The discrepancy observed in the shock‐wave foot region, between computations and measurements for the Reynolds‐stresses profiles, could be explained by considering the experimental shock‐wave oscillation and directions for future modelling work are indicated. Copyright © 2007 John Wiley & Sons, Ltd.     

Finite element analysis of two- and three-dimensional flows around square columns in tandem arrangement

Two- (2D) and three-dimensional (3D) finite element analyses for flow around two square columns in tandem arrangement were performed with various column spacings and Reynolds numbers. The computed values were compared with the wind-tunnel results in terms of the aerodynamic characteristics of the leeward column. In most 2D computations, strong vortices were formed behind the windward column, irrespective of widely changed Reynolds numbers. This was different from the experimental phenomena of equivalent spacing, so that the computed time-averaged pressure coefficients were not identical to the experimental values except when the distance between the two columns was adequately wide or narrow. On the other hand, in 3D computation, distinct differences in flow structures behind the column were observed between Reynolds numbers of 10³ and 10⁴ and the pressure coefficient in the 3D analysis with Re=10⁴ agreed well with the experimental value. Thus, the effectiveness of 3D computations and Reynolds number effects on the flow around two square columns have been confirmed. © 1998 John Wiley & Sons, Ltd.     

Flow around a NACA0018 airfoil with a cavity and its dynamical response to acoustic forcing

Trapping of vortices in a cavity has been explored in recent years as a drag reduction measure for thick airfoils. If, however, trapping fails, then oscillation of the cavity flow may couple with elastic vibration modes of the airfoil. To examine this scenario, the effect of small amplitude vertical motion on the oscillation of the shear layer above the cavity is studied by acoustic forcing simulating a vertical translation of a modified NACA0018 profile. At low Reynolds numbers based on the chord (O(10⁴)), natural instability modes of this shear layer are observed for Strouhal numbers based on the cavity width of order unity. Acoustic forcing sufficiently close to the natural instability frequency induces a strong non-linear response due to lock-in of the shear layer. At higher Reynolds numbers (above 10⁵) for Strouhal number 0.6 or lower, no natural instabilities of the shear layer and only a linear response to forcing were observed. The dynamical pressure difference across the airfoil is then dominated by added mass effects, as was confirmed by numerical simulations.     

Effects of a trapped vortex cell on a thick wing airfoil

The effects of a trapped vortex cell (TVC) on the aerodynamic performance of a NACA0024 wing model were investigated experimentally at Re = 10⁶ and 6.67×10⁵6.67\times 10^{5}. The static pressure distributions around the model and the wake velocity profiles were measured to obtain lift and drag coefficients, for both the clean airfoil and the controlled configurations. Suction was applied in the cavity region to stabilize the trapped vortex. For comparison, a classical boundary layer suction configuration was also tested. The drag coefficient curve of the TVC-controlled airfoil showed sharp discontinuities and bifurcative behavior, generating two drag modes. A strong influence of the angle of attack, the suction rate and the Reynolds number on the drag coefficient was observed. With respect to the clean airfoil, the control led to a drag reduction only if the suction was high enough. Compared to the classical boundary layer suction configuration, the drag reduction was higher for the same amount of suction only in a specific range of incidence, i.e., α = −2° to α = 6° and only for the higher Reynolds number. For all the other conditions, the classical boundary layer suction configuration gave better drag performances. Moderate increments of lift were observed for the TVC-controlled airfoil at low incidence, while a 20% lift enhancement was observed in the stall region with respect to the baseline. However, the same lift increments were also observed for the classical boundary layer suction configuration. Pressure fluctuation measurements in the cavity region suggested a very complex interaction of several flow features. The two drag modes were characterized by typical unsteady phenomena observed in rectangular cavity flows, namely the shear layer mode and the wake mode.     

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.     

Some Reynolds number effects on two- and three-dimensional turbulent boundary layers

Experimental data for a two-dimensional (2-D) turbulent boundary layer (TBL) flow and a three-dimensional (3-D) pressure-driven TBL flow outside of a wing/body junction were obtained for an approach Reynolds number based on momentum thickness of Re _θ =23,200. The wing shape had a 3:2 elliptical nose, NACA 0020 profiled tail, and was mounted on a flat wall. Some Reynolds number effects are examined using fine spatial resolution (Δy ⁺=1.8) three-velocity-component laser-Doppler velocimeter measurements of mean velocities and Reynolds stresses at nine stations for Re _θ =23,200 and previously reported data for a much thinner boundary layer at Re _θ =5,940 for the same wing shape. In the 3-D boundary layers, while the stress profiles vary considerably along the flow due to deceleration, acceleration, and skewing, profiles of the parameter correlate well and over available Reynolds numbers. The measured static pressure variations on the flat wall are similar for the two Reynolds numbers, so the vorticity flux and the measured mean velocities scaled on wall variables agree closely near the wall. The stresses vary similarly for both cases, but with higher values in the outer region of the higher Re _θ case. The outer layer turbulence in the thicker high Reynolds number case behaves similarly to a rapid distortion of the flow, since stream-wise vortical effects from the wall have not diffused completely through the boundary layer at all measurement stations. Received: 9 June 2000/Accepted: 26 January 2001     

Numerical study of vortex shedding from a rotating cylinder immersed in a uniform flow field

A numerical study is made of the unsteady two‐dimensional, incompressible flow past an impulsively started translating and rotating circular cylinder. The Reynolds number (Re) and the rotating‐to‐translating speed ratio (α) are two controlled parameters, and the influence of their different combinations on vortex shedding from the cylinder is investigated by the numerical scheme sketched below. Associated with the streamfunction (ψ)–vorticity (ω) formulation of the Navier–Stokes equations, the Poisson equation for ψ is solved by a Fourier/finite‐analytic, separation of variable approach. This approach allows one to attenuate the artificial far‐field boundary, and also yields a global conditioning on the wall vorticity in response to the no‐slip condition. As for the vorticity transport equation, spatial discretization is done by means of finite difference in which the convection terms are handled with the aid of an ENO (essentially non‐oscillatory)‐like data reconstruction process. Finally, the interior vorticity is updated by an explicit, second‐order Runge–Kutta method. Present computations fall into two categories. One with Re=10³ and α≤3; the other with Re=10⁴ and α≤2. Comparisons with other numerical or physical experiments are included. Copyright © 2000 John Wiley & Sons, Ltd.