Application of the γ-Reθ Transition Model to Simulations of the Flow Past a Circular Cylinder

Based on the finite volume method, the flow past a two-dimensional circular cylinder at a critical Reynolds number (Re = 8.5 × 10⁵) was simulated using the Navier-Stokes equations and the γ-Re_θ transition model coupled with the SST k ? ω turbulence model (hereinafter abbreviated as γ-Re_θ model). Considering the effect of free-stream turbulence intensity decay, the SST k ? ω turbulence model was modified according to the ambient source term method proposed by Spalart and Rumsey, and then the modified SST k ? ω turbulence model is coupled with the γ-Re_θ transition model (hereinafter abbreviated as γ-Re_θ-SR model). The flow past a circular cylinder at different inlet turbulence intensities were simulated by the γ-Re_θ-SR model. At last, the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers were each simulated by the γ-Re_θ-SR model, and the three flow states were analyzed. It was found that compared with the SST k ? ω turbulence model, the γ-Re_θ model could simulate the transition of laminar to turbulent, resulting in better consistency with experimental result. Compared with the γ-Re_θ model, for relatively high inlet turbulence intensities, the γ-Re_θ-SR model could better simulate the flow past a circular cylinder; however the improvement almost diminished for relatively low inlet turbulence intensities The γ-Re_θ-SR model could well simulate the flow past a circular cylinder at subcritical, critical and supercritical Reynolds numbers.     

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.     

Effects of Synthetic jets on a D-Shaped Cylinder wake at a Subcritical Reynolds Number

Effects of synthetic jets on the wake of a D-shaped cylinder is investigated experimentally at a Reynolds number Re_H= 47,000, based on incoming free-stream velocity and the cylinder height (H). The synthetic jets are introduced immediately from the upper and lower trailing edges of the cylinder. The upper and lower synthetic jets are operated in an in-phase or anti-phase mode, and at a momentum ratio C_μ= 1.0% and perturbation frequency St_A= 0.11 ?0.37. The cylinder wake with perturbation is examined in detail and compared with that without, based on smoke-wire flow visualization, pressure transducer and hotwire rake measurements, and data analyses of spectra, tempo-spatial cross-correlation and proper orthogonal decomposition (POD). Large-scale vortical structures in the cylinder wake are significantly modified by the synthetic jets perturbations, exhibiting symmetric or asymmetric patterns, depending on the perturbation frequency and phase relationship of the synthetic jets. These observations are internally correlated with the drag force variations.     

Scrutinizing URANS in Shedding Flows: The Case of Cylinder in Cross-Flow in the Subcritical Regime

To unravel the widespread perception that the RANS (Reynolds-averaged Navier-Stokes) concept is unreliable in predicting the dynamics of separated flows, we assessed the performance of two RANS closure levels, the linear eddy-viscosity (LEVM) and the second-moment (Reynolds stress, RSM) approaches in a massively separated generic flow over a bluff body. Considered is the canonical, zero-turbulence, cross-flow over an infinite cylinder with reference to our LES and the available DNS and experiments at two Reynolds numbers, Re = 3.9 × 10³ and 1.4 × 10⁵, both within the sub-critical regime with laminar separation. Both models capture successfully the vortex shedding frequency, but the low frequency modulations are detected only by the RSM. At high Reynolds numbers the RSM is markedly superior to the LEVM showing very good agreement with the LES and experimental data. The RSM, accounting naturally for the stress anisotropy and phase lag between the stress and strain eigenvectors, is especially successful in reproducing the growth rate of the turbulent kinetic energy in the initial shear layer which proved to be crucial for accurate prediction of the separation-induced transition. A scrutiny of the unsteady RANS (URANS) stress terms based on the conditional phase-averaged LES data shows a remarkable similarity of the normalized coherent and stochastic (modeled) stress components for the two Reynolds numbers considered. The mixed (cross) correlations, while non-negligible at the low Re number, diminish fast relative to the stochastic ones with increasing Reynolds number and, in the whole, are not significant to undermine the URANS concept and its applicability to high Re flows of industrial relevance.     

Separation zones in the flow near a small-angle convex corner

On the basis of an asymptotic analysis of the Navier-Stokes system of equations for large Reynolds numbers (Re → ∞), the plane incompressible fluid flow near a surface having a convex corner with a small angle 2θ* is investigated. It is shown that for θ* = O(Re^?1/4), in addition to the known solution that describes a separated flow completely localized in a thin “viscous” sublayer of the interaction region near the corner point, another solution corresponding to a flow with a developed separation zone is possible. For θ ₀ = Re^1/4 θ* = O(1), the longitudinal dimension of this zone varies from finite values up to values of the order of Re^?3/8. The nonuniqueness of the solution is established on a certain range of variation of the parameter θ ₀. The dependence of the drag coefficient on the angle θ* is found.     

Two-dimensional interaction between an incident shock and a turbulent boundary layer in the presence of an entropy layer

The results of an experimental and numerical investigation of flow and heat transfer in the region of the interaction between an incident oblique shock and turbulent boundary layers on sharp and blunt plates are presented for the Mach numbers M_∞ = 5 and 6 and the Reynolds numbers Re_∞L = 27×10⁶ and 14×10⁶. The plate bluntness and the incident shock position were varied. It is shown that the maximum Stanton number St _m in the shock incidence zone decreases with increase in the plate bluntness radius r to a certain value and then varies only slightly with further increase in r. In the case of a turbulent undisturbed boundary layer heat transfer is diminished with increase in r more slowly than in the case of a laminar undisturbed flow. In the presence of an incident shock the bluntness of the leading edge of the flat plate results in a greater decrease in the Stanton number than in the absence of the shock. With increase in the bluntness of the leading edge of the plate the separation zone first sharply lengthens and then decreases in size or remains constant.     

Low Reynolds number effects in open-channel turbulent boundary layers

The present study reports measurements of a turbulent boundary layer in an open-channel flow using fiber-optic laser Doppler anemometry. The Reynolds numbers based on momentum thickness and depth of flow are in the range 750≤Re _θ ≤2,400 and 15,300≤Re _h ≤54,200, respectively. It is shown that an accurate estimate of the wall shear stress can be made by fitting a fifth-order polynomial to the near-wall data. The effect of Reynolds number on the mean turbulence intensity and triple correlation is examined using both conventional scaling laws and the recent scaling laws proposed by George and Castillo. The present results show that different scaling laws lead to different conclusions on low Reynolds number effects.     

Drag reduction using surfactants in a rotating cylinder geometry

Turbulent Couette flow between two circular cylinders has been used for drag reduction experiments using surfactants. In the experiments presented here, only the outer cylinder rotates, the inner cylinder remains at rest and accurate measurements of the torque at the inner cylinder are measured. Water is used as a reference fluid. A drag reducing surfactant called Arquad S-50 (Akzo Nobel Surface Chemistry LLC, Chicago, Ill., USA) (5 mM)+NaSal (12.5 mM) was used as the drag reduction agent. This surfactant can reduce the drag up to 70% (a Reynolds number of about 70,000–150,000) as measured by pressure drop in a pipe flow. Experiments in Couette flow also show drag reduction in the turbulent range. Two arrangements were used, (1) one small trip-wire on the inner cylinder, and (2) four larger trip-wires on the outer cylinder. These trips reduce the critical Reynolds number for transition from laminar to turbulent flow. In case (1), we obtained 18% drag reduction at 5,000<Re<15,000 and in case (2), we obtained an average reduction of about 20% at 2,000<Re<10,000, increasing up to 30% at Re=15,000. The paper also discusses two important problems. First, the shear rate is not constant in the radial gap in circular Couette flow. For non-Newtonian fluids, where the molecular viscosity is a function of the shear rate, this effect must be considered. Second, which viscosity should be used in the Reynolds number? For pipe flow measurements, most authors use the viscosity of the solvent (generally water and Newtonian). For measurements in the Couette flow, we use a different approach, which is described in this paper. We conclude that Couette flow is a useful method for drag reduction investigations. Its advantage is the much smaller geometry in comparison to those of conventional test facilities such as wind tunnels, water, or oil channels or in tubes.     

Enhanced Dissipation and Inviscid Damping in the Inviscid Limit of the Navier–Stokes Equations Near the Two Dimensional Couette Flow

In this work we study the long time inviscid limit of the two dimensional Navier–Stokes equations near the periodic Couette flow. In particular, we confirm at the nonlinear level the qualitative behavior predicted by Kelvin’s 1887 linear analysis. At high Reynolds number Re, we prove that the solution behaves qualitatively like two dimensional Euler for times \({{t \lesssim Re^{1/3}}}\), and in particular exhibits inviscid damping (for example the vorticity weakly approaches a shear flow). For times \({{t \gtrsim Re^{1/3}}}\), which is sooner than the natural dissipative time scale O(Re), the viscosity becomes dominant and the streamwise dependence of the vorticity is rapidly eliminated by an enhanced dissipation effect. Afterwards, the remaining shear flow decays on very long time scales \({{t \gtrsim Re}}\) back to the Couette flow. When properly defined, the dissipative length-scale in this setting is \({{\ell_D \sim Re^{-1/3}}}\), larger than the scale \({{\ell_D \sim Re^{-1/2}}}\) predicted in classical Batchelor–Kraichnan two dimensional turbulence theory. The class of initial data we study is the sum of a sufficiently smooth function and a small (with respect to Re^?1) L² function.     

Passive Boundary Layer Flow Control Using Porous Lamination

The flow over a porous laminated flat plate is investigated from a flow control perspective through experiments and computations. A square array of circular cylinders is used to model the porous lamination. We determine the velocities at the fluid–porous interface by solving the two-dimensional Navier–Stokes and the continuity equations using a staggered flow solver and using LDV in experiments. The control parameters for the porous region are porosity, \(\phi \) and Reynolds number, Re, based on the diameter of the circular cylinders used to model the porous lamination. Computations are conducted for \(0.4< \phi < 0.9\) and \(25< Re < 1000\), and the experiments are conducted for \(\phi = 0.65\) and 0.8 at \(Re \approx 391,\ 497\) and 803. The permeability of the porous lamination is observed to induce a slip velocity at the interface, effectively making it a slip wall. The slip velocity is seen to be increasing functions of \(\phi \) and Re. For higher porosities at higher Re, the slip velocity shows non-uniform and unsteady behavior and a breakdown Reynolds number is defined based on this characteristic. A map demarcating the two regimes of flow is drawn from the computational and experimental data. We observe that the boundary layer over the porous lamination is thinner than the Blasius boundary layer and the shear stress is higher at locations over the porous lamination. We note that the porous lamination helps maintain a favorable pressure gradient at the interface which delays separation. The suitable range of porosities for effective passive separation control is deduced from the results.     

Unsteady characteristics of low-Re flow past two tandem cylinders

This study investigated the two-dimensional flow past two tandem circular or square cylinders at Re = 100 and D / d = 4–10, where D is the center-to-center distance and d is the cylinder diameter. Numerical simulation was performed to comparably study the effect of cylinder geometry and spacing on the aerodynamic characteristics, unsteady flow patterns, time-averaged flow characteristics and flow unsteadiness. We also provided the first global linear stability analysis and sensitivity analysis on the physical problem for the potential application of flow control. The objective of this work is to quantitatively identify the effect of the cylinder geometry and spacing on the characteristic quantities. Numerical results reveal that there is wake flow transition for both geometries depending on the spacing. The characteristic quantities, including the time-averaged and fluctuating streamwise velocity and pressure coefficient, are quite similar to that of the single cylinder case for the upstream cylinder, while an entirely different variation pattern is observed for the downstream cylinder. The global linear stability analysis shows that the spatial structure of perturbation is mainly observed in the wake of the downstream cylinder for small spacing, while moves upstream with reduced size and is also observed after the upstream cylinder for large spacing. The sensitivity analysis reflects that the temporal growth rate of perturbation is the most sensitive to the near-wake flow of downstream cylinder for small spacing and upstream cylinder for large spacing.     

Taguchi design of three dimensional simulations for optimization of turbulent mixed convection in a cavity

This study discusses the application of Taguchi method in assessing maximum heat transfer rate for the turbulent mixed convection in an enclosure embedded with rotating isothermal cylinder. The simulations were planned based on Taguchi’s L16 orthogonal array with each trial performed under different conditions of position of the cylinder, Reynolds number (Re) and Rayleigh number (Ra). The thermal lattice Boltzmann based on D3Q19 methods without any turbulent submodels was purposed to simulate the flow and thermal fields. A relaxation time method with the stability constants is introduced to solve turbulent natural convection problems. Signal-to-noise ratios (S/N) analysis were carried out in order to determine the effects of process parameters and optimal factor settings. Finally, confirmation tests verified that Taguchi method achieved optimization of heat transfer rate with sufficient accuracy.     

Effect of Inflow Turbulence on an Airfoil Flow with Laminar Separation Bubble: An LES Study

The present paper is concerned with numerical investigations on the effect of inflow turbulence on the flow around a SD7003 airfoil. At a Reynolds number Re_c =?60,000, an angle of attack α =?4^° and a low or zero turbulence intensity of the oncoming flow, the flow past the airfoil is known to be dominated by early separation, subsequent transition and reattachment leading to a laminar separation bubble with a distinctive pressure plateau. The objective of the study is to investigate the effect of inflow turbulence on the flow behavior. For this purpose, a numerical methodology relying on a wall-resolved large-eddy simulation, a synthetic turbulence inflow generator and a specific source term concept for introducing the turbulence fluctuations within the computational domain is used. The numerical technique applied allows the variation of the free-stream turbulence intensity (TI) in a wide range. In order to analyze the influence of TI on the arising instantaneous and time-averaged flow field past the airfoil, the present study evaluates the range 0% ≤ TI ≤?11.2%, which covers typical values found in atmospheric boundary layers. In accordance with experimental studies it is shown that the laminar separation bubble first shrinks and finally completely vanishes for increasing inflow turbulence. Consequently, the aerodynamic performance in terms of the lift-to-drag ratio increases. Furthermore, the effect of the time and length scales of the isotropic inflow turbulence on the development of the flow field around the airfoil is analyzed and a perceptible influence is found. Within the range of inflow scales studied decreasing scales augment the receptivity of the boundary layer promoting an earlier transition.     

Turbulent Drag Reduction by a Near Wall Surface Tension Active Interface

In this work we study the turbulence modulation in a viscosity-stratified two-phase flow using Direct Numerical Simulation (DNS) of turbulence and the Phase Field Method (PFM) to simulate the interfacial phenomena. Specifically we consider the case of two immiscible fluid layers driven in a closed rectangular channel by an imposed mean pressure gradient. The present problem, which may mimic the behaviour of an oil flowing under a thin layer of different oil, thickness ratio h₂/h₁ =?9, is described by three main flow parameters: the shear Reynolds number Re _τ (which quantifies the importance of inertia compared to viscous effects), the Weber number We (which quantifies surface tension effects) and the viscosity ratio λ = ν₁/ν₂ between the two fluids. For this first study, the density ratio of the two fluid layers is the same (ρ₂ = ρ₁), we keep Re _τ and We constant, but we consider three different values for the viscosity ratio: λ =?1, λ =?0.875 and λ =?0.75. Compared to a single phase flow at the same shear Reynolds number (Re _τ =?100), in the two phase flow case we observe a decrease of the wall-shear stress and a strong turbulence modulation in particular in the proximity of the interface. Interestingly, we observe that the modulation of turbulence by the liquid-liquid interface extends up to the top wall (i.e. the closest to the interface) and produces local shear stress inversions and flow recirculation regions. The observed results depend primarily on the interface deformability and on the viscosity ratio between the two fluids (λ).     

Experimental Investigation of the Transverse Self-Oscillations of Circular Cylinders Mounted with a Narrow Clearance in a Plane Channel

The results of an investigation of the flow past and the behavior of free bluff bodies mounted in pipes and channels with a narrow clearance, conducted in the Institute of Mechanics of Moscow State University, are presented. The drag of circular cylinders of different size and mass in a circulation-free water flow in a plane channel of rectangular cross-section was studied in the transverse self-oscillation regime. The experiments were conducted for Reynolds numbers based on the cylinder diameter 1.7 ? 10⁴ ≤ Re ≤ 7.2 ? 10⁴, relative clearances \(\bar S\) based on a cross-sectional area ranging from 0.76 to 0.9, and cylinder-to-water density ratios ρ _c /ρ ranging from 1.29 to 8.2. Only the case of intense transverse self-oscillations accompanied by impact interaction with the channel wall was considered. The dependence of the period-average cylinder drag coefficient C _x on the basic dimensionless relevant parameters is obtained. The dependence of the dimensionless self-oscillation frequency determined in [1] is refined. The kinematic and dynamic features of the flows past spheres in cylindrical pipes and cylinders in plane channels are compared in the transverse self-oscillation regime.     

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall

Hydrodynamic properties of squirmer swimming in power-law fluid near a wall considering the interaction between squirmer and wall are numerically studied with an immersed boundary-lattice Boltzmann method. The power-law index, Reynolds number, initial orientation angle of squirmer, and initial distance of squirmer from the wall are all taken into account to investigate the swimming characteristics for pusher (β?<?0), neutral squirmer (β?=?0), and puller (β?>?0) (three kinds of swimmer types) near the no-slip boundary. Four new kinds of swimming modes are found. Results show that, for the pushers and pullers, the wall displays an increasing attraction with increasing power-law index n, which differs from the neutral squirmer who always departs from the wall after the first collision with the wall. Both the initial orientation angle and initial distance from the wall only affect the moving situations rather than the moving modes of the squirmers. However, the squirmers depart from the wall as the Reynolds number increases and chaotic orbits appear for some squirmers at Re?=?5. Several typical flow fields are analyzed and the power consumption and torque for different kinds of flows are also studied. It is found that, as the absolute value of β increases, the power consumption generally increases in shear-thinning (n?=?0.4), Newtonian (n?=?1), and shear-thickening (n?=?1.6) fluids. Moreover, the pushers (β?<?0) and the pullers (β?>?0) expend almost the same power if the absolute value of β remains the same. In addition, the power consumption of the squirmers is highly dependent on the power-law index n.     

Model-based control of vortex shedding at low Reynolds numbers

Model-based feedback control of vortex shedding at low Reynolds numbers is considered. The feedback signal is provided by velocity measurements in the wake, and actuation is achieved using blowing and suction on the cylinder’s surface. Using two-dimensional direct numerical simulations and reduced-order modelling techniques, linear models of the wake are formed at Reynolds numbers between 45 and 110. These models are used to design feedback controllers using \(\mathcal {H}_\infty \) loop-shaping. Complete suppression of shedding is demonstrated up to Re \(=\) 110—both for a single-sensor arrangement and for a three-sensor arrangement. The robustness of the feedback controllers is also investigated by applying them over a range of off-design Reynolds numbers, and good robustness properties are seen. It is also observed that it becomes increasingly difficult to achieve acceptable control performance—measured in a suitable way—as Reynolds number increases.     

A modified experimental system for measurement of Magnus loads

A modified and improved system for the measurement of Magnus loads on spinning models, at various spin rates and at various angles of attack is presented. The objective of the modifications was to increase the measuring-system sensitivity to the Magnus loads, while increasing its stiffness relative to the other loads. The spin mechanism was designed to minimize interference with the outer flow. The system was designed to operate in conjunction with any standard balance (of the appropriate diameter), so that it could be used at all wind-tunnel speeds, from low subsonic to supersonic. The versatility of the new system was demonstrated in successful measurements at various Mach numbers. Test Mach numbers varied from M = 0.025 to M = 0.4, with spin rates from zero to 30,000 rpm, and angles of attack from = 0° to = 90°. The test results show that Magnus loads can be measured accurately, both at low and high angles of attack and spin rates.List of symbols C _y side force coefficient - C _y/sin² Magnus force coefficient - D body diameter - f spin rate (dimensional) - L model length - M Mach number - P nondimensional spin rate - Re Reynolds number based on body diameter Re = VD/ - Re _cf cross flow Reynolds number Re _cf = Re sin - rpm revolutions per minute - V free stream velocity - X _cp center of pressure location, measured from the nose - angle of attack     

Evaluation of the Spectral Element Dynamic Model for Large-Eddy Simulation on Unstructured,Deformed Meshes

In this paper, the Spectral-Element Dynamic Model (SEDM), suited for Large-Eddy Simulation (LES) using Discontinuous Finite Element Methods (DFEM), is assessed using unstructured meshes. Five test cases of increasing complexity are considered, namely, the Taylor-Green vortex at Re =?5000, the turbulent channel flow at Re_τ =?587, the circular cylinder in cross-flow at Re_D =?3900, the square cylinder in cross-flow at Re_D =?22400 and the channel with periodic constrictions at Re_h =?10595. Various discretization parameters such as the grid spacing, polynomial degree and numerical flux are assessed and very accurate results are reported in all cases. This consistency in the results demonstrates the versatility of the SEDM approach and its ability to gage the actual resolution and quality of the mesh and, accordingly, to introduce an amount of sub-grid dissipation which is adapted to the spatial discretization considered.     

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.