Unsteady flow structure and vorticity convection over the airfoil oscillating at high reduced frequency

Unsteady vortex structures and vorticity convection over the airfoil (NACA 0012), oscillating in the uniform inflow, are studied by flow visualization and velocity measurements. The airfoil, pivoting at one-third of the chord, oscillates periodically near the static stalling angle of attack (AOA) at high reduced-frequency. The phase-triggering and modified phase-averaged techniques are employed to reconstruct the pseudo instantaneous velocity field over the airfoil. During the down stroke cycle, the leading-edge separation vortex is growing and the vortex near the trailing edge begins to shed into the wake. During the upstroke cycle, the leading-edge separation vortex is matured and moves downstream, and the counter clockwise vortex is forming near the trailing edge. Convection speeds and wavelength of the unsteady vortex structure over the airfoil equal to that of the counter clockwise vortex shed into the wake. This kind of vortex structure is termed as “synchronized shedding” type. The wavelength of unsteady vortex structure over the airfoil is significantly different from that at low reduced-frequency. Consistent convection speeds of the leading-edge separation vortex are acquired from the spatial-temporal variations of local circulation and local surface vorticity generation, and equals that predicted from flow visualization. Spatial-temporal variations of the local surface vorticity generation clearly reveal the formation and passage of the leading-edge separation vortex only in the region where the flow does not separate completely from the surface. Significant amounts of the surface vorticity are generated within the leading-edge region of the airfoil during the upstroke cycle. Only negligible amount of surface vorticity is produced within the region of complete flow separation. During the down stroke cycle, the surface vorticity generation is mild along the airfoil surface, except the leading-edge region where a small scale leading-edge separation vortex is forming and growing.     

Drag reduction and thrust generation by tangential surface motion in flow past a cylinder

Sensitivity of drag to tangential surface motion is calculated in flow past a circular cylinder in both two- and three-dimensional conditions at Reynolds number \(\textit{Re} \le 1000\). The magnitude of the sensitivity maximises in the region slightly upstream of the separation points where the contour lines of spanwise vorticity are normal to the cylinder surface. A control to reduce drag can be obtained by (negatively) scaling the sensitivity. The high correlation of sensitivities of controlled and uncontrolled flow indicates that the scaled sensitivity is a good approximation of the nonlinear optimal control. It is validated through direct numerical simulations that the linear range of the steady control is much higher than the unsteady control, which synchronises the vortex shedding and induces lock-in effects. The steady control injects angular momentum into the separating boundary layer, stabilises the flow and increases the base pressure significantly. At \(\textit{Re}=100\), when the maximum tangential motion reaches 50% of the free-stream velocity, the vortex shedding, boundary-layer separation and recirculation bubbles are eliminated and 32% of the drag is reduced. When the maximum tangential motion reaches 2.5 times of the free-stream velocity, thrust is generated and the power savings ratio, defined as the ratio of the reduced drag power to the control input power, reaches 19.6. The mechanism of drag reduction is attributed to the change of the radial gradient of spanwise vorticity (\(\partial _{r} \hat{\zeta }\)) and the subsequent accelerated pressure recovery from the uncontrolled separation points to the rear stagnation point.     

Modeling the flow past an oscillating airfoil by the method of viscous vortex domains

The Lagrangian vortex method for solving the Navier-Stokes equations is applied for numerically modeling the unsteady flow past a wing airfoil executing angular oscillations in a viscous incompressible flow. Formulas relating the unsteady forces on the airfoil and the vorticity field are derived. The calculated results are compared with the experimental data for the NACA-0012 airfoil executing harmonic oscillations in an air flow at the Reynolds number Re = 4.4 × 10⁴.     

INFLUENCE OF SEPARATION ON SOUND GENERATED BY VORTEX-STEP INTERACTION

An analysis of the sound produced when a line vortex interacts at low Mach number with forward or backward facing steps is made. The radiation is dominated by an aeroacoustic dipole whose strength is equal to the unsteady drag on the step. The drag is determined by the vorticity distribution, and a correct estimate of the sound must therefore include contributions from vorticity in the separated flow induced by the vortex. The separation is modelled by assuming that the shed vorticity rolls up into a concentrated core, fed by a connecting sheet from the edge of the step of negligible circulation. The motion everywhere is irrotational except at the impinging vortex and the separation core, and the trajectory of the core is governed by an emended Brown & Michael equation. For large steps it is found that estimates of the generated sound that neglect separation are typically an order of magnitude too large. The sound levels predicted for small steps with and without separation are of comparable magnitudes, although the respectivephasesare different.Turbulentflow over a step frequently involves separation and large surface pressure fluctuations at reattachment zones. The results of this paper suggest that numerical schemes for determining the noise generated by turbulent flow over a step must take proper account of “forcing” of the separation region by the impinging turbulence and of vorticity production via the no-slip condition.     

A Parallel 3D Unsteady Flow Analysis in an Asymmetrically Constricted Vessel

In this study, parallel computation of unsteady incompressible flow in an asymmetrically constricted 3D vessel has been presented. A time accurate cell centered finite volume method (FVM) in conjunction with pseudo-compressibility technique and Roe's flux difference splitting of nonlinear terms has been employed for solving the Navier-Stokes (NS) equations on the multiple instruction multiple data (MIMD) machine VPP700. The influence of Reynolds' number ( Re ) and the Strouhal number ( St ) on flow dynamic factors like wall pressure (WP), wall shear stress (WSS), central axis velocity (CAV), etc., have been analyzed. Three-dimensional (3D) features in the formation and detachment of separation zones, which are sensitive to both Re and St have been noticed on the diverging wall of the constriction.     

Interaction of zones of flow separation in a centrifugal impeller-stationary vane system

In a radial flow pump operating in off-design conditions, regions of stall can exist on the rotating impeller blade and on the downstream diffuser blade, vane or tongue. Interaction of these stall zones can generate complex patterns of vorticity concentrations. In turn, these vorticity concentrations are related to sources of unsteady stagnation enthalpy. The form of these patterns is strongly dependent on the instantaneous location of the impeller trailing-edge relative to the leading-edge of the vane.Comparison of instantaneous with ensemble-averaged images shows that the flow structure in the gap region between the impeller and the vane is highly repetitive. Away from this region, in particular in the separated shear layer from the vane, the nonrepetitive nature of the vorticity field is manifested in substantial reduction of peak levels of vorticity in the ensemble-averaged image, relative to the instantaneous image.The three-dimensional flow structure resulting from these separation zone interactions was characterized via end views of the flow patterns. Particularly pronounced concentrations of vorticity can occur in this plane. They tend to be located in the shear layer at the outer edge of the large-scale separation zone. These vorticity concentrations are, however, highly non-stationary for successive passages of the impeller blade. Ensemble-averaging reveals that they persist primarily on the endwalls of the diffuser.The authors are grateful to the Office of Naval Research for support of this research program     

Parametric study of unsteady laminar flows

Results of a parametric study of unsteady laminar flows are analyzed. Three-dimensional unsteady equations of hydromechanics for a compressible medium are solved. The range of the characteristic Reynolds number Re = 400–900 is considered. It is demonstrated that the laminar flow in a plane channel ceases to be steady at Re = 415. As the Reynolds number increases, the unsteady processes become more intense, disturbances penetrate inward the channel, and separation zones lose their stability. In the vicinity of the channel exit, however, the flow tends to stabilize, though it remains unsteady. No transition to a turbulent flow occurs in the examined range of Reynolds numbers.     

Effect of homogeneous stratification of the fluid on the dynamics of a kármán vortex street behind a circular cylinder

The experimental results of studying the effect of homogeneous stratification of the fluid on the conditions of generation of a Kárman vortex street [1] developing in the wake of a cylinder in steady horizontal motion are described. In a homogeneous medium at Reynolds numbers Re >5 two symmetrical regions of vorticity of opposite sign are formed behind the cylinder and move together with the latter. As the speed of the cylinder increases, the link between the vortices and the cylinder grows weaker, the vortices are stretched out along the flow and at Re > 40 begin to separate alternately, forming a vortex street in the wake. At first, the frequency of vortex separation increases sharply with increase in Re, but then levels off. It is found that in a uniformly stratified fluid the onset of vortex separation from the moving cylinder is delayed. The dependence of the critical Reynolds number (onset of vortex separation) on the internal Froude number is obtained. The effect of stratification of the fluid on the frequency of separation of the vortices in the Kármán street is investigated. The effect of the Froude number on the dependence of the Strouhal number on the Reynolds number is established.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 83–86, January–February, 1986.In conclusion the authors wish to thank A. T. Onufriev for his interest in their work and useful discussions of the results.     

Large-scale vortical structure of turbulent separation bubble affected by unsteady wake

The large-scale vortical structure of a turbulent separation bubble under the influence of an unsteady wake was investigated. The unsteady wake was generated by a spoked-wheel type wake generator installed in front of the separation bubble. This wake generator was rotated either clockwise or counter-clockwise at Re _H=5.600. The mechanism of vortex shedding from the separation bubble was analyzed in detail by taking a conditional average as well as a phase average. Spatial box filtering (SBF) was used to extract the large-scale vortical structure from the turbulent separation bubble affected by the unsteady wake. To elucidate the influence of the unsteady wake on the large-scale vortical structure, conditional averages of the velocity, vorticity and turbulent kinetic energy were calculated. The nature of the convection of the vortical structure under the influence of an unsteady wake was analyzed. The dipole acoustic pressure level was predicted using Curle's integral of wall-pressure fluctuations.     

Unsteady flow calculation in a radial flow centrifugal pump with spiral casing

The prediction of the two-dimensional unsteady flow established in a radial flow centrifugal pump is considered. Assuming the fluid incompressible and inviscid, the velocity field is represented by means of source and vorticity surface distributions as well as a set of point vortices. Using this representation, a grid-free (Lagrangian) numerical method is derived based on the coupling of the boundary element and vortex particle methods. In this context the source and vorticity surface distributions are determined through the non-entry boundary condition together with the unsteady Kutta condition. In order to satisfy Kelvin's theorem, vorticity is shed at the trailing edges of the impeller blades. Then the vortex particle method is used to approximate the convection of the free vorticity distribution. Results are given for a pump configuration experimentally tested by Centre Technique des Industries Mécaniques (CETIM). Comparisons between predictions and experimental data show the capability of the proposed method to reproduce the main features of the flow considered.     

Planar velocity measurements of a two-dimensional compressible wake

The present study describes the application of particle image velocimetry (PIV) to investigate the compressible flow in the wake of a two-dimensional blunt base at a freestream Mach number M_X=2. The first part of the study addresses specific issues related to the application of PIV to supersonic wind tunnel flows, such as the seeding particle flow-tracing fidelity and the measurement spatial resolution. The seeding particle response is assessed through a planar oblique shock wave experiment. The measurement spatial resolution is enhanced by means of an advanced image-interrogation algorithm. In the second part, the experimental results are presented. The PIV measurements yield the spatial distribution of mean velocity and turbulence. The mean velocity distribution clearly reveals the main flow features such as expansion fans, separated shear layers, flow recirculation, reattachment, recompression and wake development. The turbulence distribution shows the growth of turbulent fluctuations in the separated shear layers up to the reattachment location. Increased velocity fluctuations are also present downstream of reattachment outside of the wake due to unsteady flow reattachment and recompression. The instantaneous velocity field is analyzed seeking coherent flow structures in the redeveloping wake. The instantaneous planar velocity and vorticity measurements return evidence of large-scale turbulent structures detected as spatially coherent vorticity fluctuations. The velocity pattern consistently shows large masses of fluid in vortical motion. The overall instantaneous wake flow is organized as a double row of counter-rotating structures. The single structures show vorticity contours of roughly elliptical shape in agreement with previous studies based on spatial correlation of planar light scattering. Peak vorticity is found to be five times higher than the mean vorticity value, suggesting that wake turbulence is dominated by the activity of large-scale structures. The unsteady behavior of the reattachment phenomenon is studied. Based on the instantaneous flow topology, the reattachment is observed to fluctuate mostly in the streamwise direction suggesting that the unsteady separation is dominated by a pumping-like motion.     

Nonlinear interaction of a vortex pair with clean and surfactant-covered free surfaces

Two-dimensional unsteady viscous-flow problem associated with the normal incidence of a counter-rotating vortex pair on a free surface is analyzed. Effects of surface tension and insoluble surfactants on the generation of free-surface vorticity and surface waves are investigated. A recently developed finite-difference method based on boundary-fitted coordinates is used to solve the fully-nonlinear problem. Results show that in the absence of surfactants and at low Froude number (based on circulation strength and initial separation distance of the vortex pair), waves of short lengths are generated. However, secondary vorticity generated in this case is not strong enough to affect the outward translation of the primary vortices. At intermediate Froude number, a transient wave developing outboard of the primary vortex becomes steep, and eventually breaks because of local instability. Consequently, free-surface vorticity inhibits the outward translation of the primary vortices. Surface tension in a clean free surface dampens the steep short waves, hence also the generation of free-surface vorticity. However, variation in surface tension induced by surfactants intensifies the generation of surface vorticity, thereby causing the primary vortices to rebound. The increase in the rotational part of wave motion results in the dampening of overall free-surface deformations. However, it is found that the shear stress associated with a large gradient of surfactant concentration could cause local steepening of the short wave generated outboard of the primary vortex.     

A three‐dimensional vortex particle‐in‐cell method for vortex motions in the vicinity of a wall

A new vortex particle‐in‐cell method for the simulation of three‐dimensional unsteady incompressible viscous flow is presented. The projection of the vortex strengths onto the mesh is based on volume interpolation. The convection of vorticity is treated as a Lagrangian move operation but one where the velocity of each particle is interpolated from an Eulerian mesh solution of velocity–Poisson equations. The change in vorticity due to diffusion is also computed on the Eulerian mesh and projected back to the particles. Where diffusive fluxes cause vorticity to enter a cell not already containing any particles new particles are created. The surface vorticity and the cancellation of tangential velocity at the plate are related by the Neumann conditions. The basic framework for implementation of the procedure is also introduced where the solution update comprises a sequence of two fractional steps. The method is applied to a problem where an unsteady boundary layer develops under the impact of a vortex ring and comparison is made with the experimental and numerical literature. Copyright © 2001 John Wiley & Sons, Ltd.     

Simulation of cross-flow-induced vibration of cylinder arrays by surface vorticity method

The surface vorticity method (SVM), which is a fast and practical grid-free two-dimensional (2-D) method, and a fluid–structure interaction model incorporating the effects of cylinder motions and displacements is used to simulate the vortex-induced vibration of cylinder arrays at sub-critical Reynolds number Re=2.67×10⁴. The SVM is found to be most suitable for simulating a 2-D cylinder row with large-amplitude vibrations where the vorticity field and the fluid forces of the cylinder row change drastically, and the effect of the stream on the transverse direction vibration is very significant. The fluidelastic instability of a flexible cylinder row at small pitch ratio is also investigated, and the critical reduced velocity of the cylinder row at a reduced damping parameter S_G=1.29 is calculated, which is in good agreement with experimental and analytical results of the unsteady model. Vortex-induced vibration of a staggered cylinder array is simulated using different structural parameters. When the cylinders are relatively more flexible, the flow pattern changes dramatically and the fluid–structure interaction has a dominant impact on the flow field. Compared with grid-based methods, the grid-free SVM is a fast and practical method for the simulation of the FIV of cylinder arrays due to vortex shedding at sub-critical Reynolds numbers.     

Surrogate-based aeroelastic loads prediction in the presence of shock-induced separation

A surrogate-based modeling strategy is presented for robust and efficient prediction of unsteady aeroelastic loads in the presence of shock-induced separation. Enriched piston theory predictions are extended with a data-driven nonlinear autoregressive with exogenous inputs model to account for hysteresis from the interplay of a dynamically deforming surface with the separation bubble in a shock/boundary layer interaction. The approach is evaluated for prescribed surface motion and shock-induced panel flutter responses, with good agreement observed in each scenario relative to unsteady Reynolds-averaged Navier–Stokes simulations. For the latter, excellent agreement is observed in the prediction of the stability boundary and oscillation frequency. In contrast, the oscillation amplitude conservatively deviates from the Reynolds-averaged Navier–Stokes solution with increasing dynamic pressure. The online computational cost of the extended approach is orders of magnitude less than that required for predictions using an unsteady Reynolds-averaged Navier–Stokes model.     

Reynolds number dependence of two-dimensional laminar co-rotating vortex merging

The paper reports unsteady Navier–Stokes calculations of laminar two-dimensional co-rotating vortex merging for various Reynolds numbers. The unsteady, incompressible two-dimensional Navier–Stokes equations were solved with fourth-order Runge–Kutta temporal discretization and fourth-order symmetric compact schemes for spatial discretization. Calculations of the unsteady Taylor vortex benchmark showed that fourth-order accurate solutions for all primitive variables were indeed achieved. Calculations for a pair of equal-strength co-rotating vortices show good agreement with reported direct numerical simulation and experiments for the evolution of the separation distance and core radius. It is found that the time required for merging is inversely proportional to the square root of the Reynolds number. According to previous experimental research, it was also found that complete merging in laminar regime undergoes four stages with physical meaning. The physical mechanism responsible for the merging process is investigated and it is found that the antisymmetric vorticity dynamics plays an important role until full merging.     

Simultaneous quantitative flow measurement using PIV on both sides of the air–water interface for breaking waves

An experimental method for simultaneously measuring the velocity fields on the air and water side of unsteady breaking waves is presented. The method includes a novel technique for seeding the air flow such that the air velocity can be resolved in the absence of wind. Low density particles that have large Stokes drag and ability to respond to high-frequency flow fluctuations are used to seed the air flow. Multi-camera, multi-laser particle image velocimetry setups are applied to small-scale unsteady breaking waves, yielding fully time-resolved velocity fields. The surface tension of the fluid is altered and controlled to form spilling breaking waves. Results for the velocity and vorticity fields of representative spilling breakers, which show shedding of an air-side vortex and well-documented generation of water-side vorticity, are presented and discussed.     

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

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

Interactions of three-dimensional vicous axisymmetric vortex rings with a free surface

The unsteady nonlinear interaction of three-dimensional vortices with a free surface is a great challenge in fluid mechanics, which has deep theoretical significance and important practical background. Applying the three-dimensional VOF method, the interactions of three-dimensional axisymmetric vortex rings with a free surface in an incompressible viscous fluid are numerically simulated. The influence of the Froude number and the surface tension are studied and the evolution of the vorticity, the trajectories of the vortex rings and the baroclinic vorticity on the surface are obtained. The results agreed well with the experiments reported in the literature. The project supported by the National Natural Science Foundation of China     

Three‐dimensional vortex simulation of unsteady flow in hydraulic turbines

On the basis of the Helmholtz decomposition, a grid‐free numerical scheme is provided for the solution of unsteady flow in hydraulic turbines. The Lagrangian vortex method is utilized to evaluate the convection and stretch of the vorticity, and the BEM is used to solve the Neumann problem to define the potential flow. The no‐slip boundary condition is satisfied by generating vortex sticks at the solid surface. A semi‐analytical regularization technique is applied to evaluate the singular boundary surface integrals of the potential velocity and its gradients accurately. The fast multipole method was extended to evaluate the velocity and velocity gradients induced by the discretized vortex blobs in the Lagrangian vortex method. The successful simulation for the unsteady flow through a hydraulic turbine's runner has manifested the effectiveness of the proposed method. Copyright © 2011 John Wiley & Sons, Ltd.