Three-dimensional numerical simulation of the wake flow of an afterbody at subsonic speeds

We numerically investigate the wake flow of an afterbody at low Reynolds number in the incompressible and compressible regimes. We found that, with increasing Reynolds number, the initially stable and axisymmetric base flow undergoes a first stationary bifurcation which breaks the axisymmetry and develops two parallel steady counter-rotating vortices. The critical Reynolds number (Re _cs) for the loss of the flow axisymmetry reported here is in excellent agreement with previous axisymmetric BiGlobal linear stability (BiGLS) results. As the Reynolds number increases above a second threshold, Re _co, we report a second instability defined as a three-dimensional peristaltic oscillation which modulates the vortices, similar to the sphere wake, sharing many points in common with long-wavelength symmetric Crow instability. Both the critical Reynolds number for the onset of oscillation, Re _co, and the Strouhal number of the time-periodic limit cycle, St_sat, are substantially shifted with respect to previous axisymmetric BiGLS predictions neglecting the first bifurcation. For slightly larger Reynolds numbers, the wake oscillations are stronger and vortices are shed close to the afterbody base. In the compressible regime, no fundamental changes are observed in the bifurcation process. It is shown that the steady state planar-symmetric solution is almost equal to the incompressible case and that the break of planar symmetry in the vortex shedding regime is retarded due to compressibility effects. Finally, we report the developments of a low frequency which depends on the afterbody aspect ratio, as well as on the Reynolds and on the Mach number, prior to the loss of the planar symmetry of the wake.     

Two- and three-dimensional laminar flows between disks co-rotating in a fixed cylindrical enclosure

A numerical investigation is performed for the constant property laminar flow of air in the space between a pair of disks clamped co-axially on a central hub and co-rotating in a stationary cylindrical enclosure. Both two- and three-dimensional flow conditions are examined in relation to the interdisk spacing, H, and the disk angular velocity, Ω. Two interdisk spacings are considered, corresponding to aspect ratios Γ = 0.186 and 0.279 (with Γ = H/(R₂+a−R), where R₂ is the disk radius, a is the disk rim–enclosure wall clearance, and R is the hub radius). A range of rotational speeds encompassing the transition from axisymmetric two-dimensional steady flow to non-axisymmetric three-dimensional unsteady flow are considered for various values of the Reynolds number, Re (with $ Re=\Omega R_2^2/v $, where v is the kinematic viscosity of air). Axisymmetric calculations are first performed for both aspect ratios in the range 3858≤Re≤23 150. Fully three-dimensional calculations are then performed for the configuration with Γ = 0.186 and Re = 23 150, and for the configuration with Γ = 0.279 and Re = 7715, 15 430 and 23 150. The axisymmetric calculations performed with Γ = 0.186 confirm many known features of the flow, including the transition from a steady flow to an oscillatory periodic regime. This occurs at ≈Re = 23 150 for a configuration with a/H = 0, and at ≈Re = 14 670 for one with a/H = 0.28 and a finite disk thickness (b/H = 0.2). Three-dimensional calculations performed for Γ = 0.186 with a/H = 0 and Re = 23 150 reveal a circumferentially periodic flow pattern with eight foci of intensified axial component of vorticity. The axisymmetric calculations performed with Γ = 0.279 and Re > 7715 yield a novel, non-unique steady solution for the velocity field that is asymmetric with respect to the interdisk mid-plane. No experimental verification of this finding exists to date, but similar situations are known to arise in the context of anomalous modes of the Taylor–Couette flow. Relaxing the axisymmetry constraint allows this flow to evolve to an oscillatory three-dimensional regime of increasing irregularity with increasing rotational speed. In this case, the number of foci of intensified axial vorticity varies with time, ranging from six at Re = 7715 to between six and eight at Re = 23 150. © 1998 John Wiley & Sons, Ltd.     

Multiplicity of steady solutions in a two-sided lid-driven cavity with different aspect ratios

This study presents a continuation method to calculate flow bifurcation in a two-sided lid-driven cavity with different aspect ratios for anti-parallel motion. In anti-parallel motion, the top and bottom walls of the cavity move in opposite directions simultaneously, while the two walls both moving to the right give parallel motion at the same speed. Comprehensive bifurcation diagrams of the cavity flows with different aspect ratios of the cavities are derived via Keller’s continuation method, and linear- stability analysis is used to identify the nature of the various flow solutions. The Reynolds number (1 ≤ Re ≤ 1,200) is used as the continuation parameter to trace the solution curves. In anti-parallel motion, the evolution of the bifurcation diagrams in cases with different aspect ratios (1 ≤ AR ≤ 2.5) is illustrated. Two stable symmetric flows and one stable asymmetric flow are identified, and the existent regions of the stable flows in the aspect ratios and Reynolds numbers are distinguished. The newly found asymmetric flow state can be obtained at a high aspect ratio and a low Reynolds number.     

Simulations of Turbulent Flow in a Plane Asymmetric Diffuser

Large-eddy simulations (LES) of a planar, asymmetric diffuser flow have been performed. The diverging angle of the inclined wall of the diffuser is chosen as 8.5°, a case for which recent experimental data are available. Reasonable agreement between the LES and the experiments is obtained. The numerical method is further validated for diffuser flow with the diffuser wall inclined at a diverging angle of 10°, which has served as a test case for a number of experimental as well as numerical studies in the literature (LES, RANS). For the present results, the subgrid-scale stresses have been closed using the dynamic Smagorinsky model. A resolution study has been performed, highlighting the disparity of the relevant temporal and spatial scales and thus the sensitivity of the simulation results to the specific numerical grids used. The effect of different Reynolds numbers of the inflowing, fully turbulent channel flow has been studied, in particular, Re _b  = 4,500, Re _b  = 9,000 and Re _b  = 20,000 with Re _b being the Reynolds number based on the bulk velocity and channel half width. The results consistently show that by increasing the Reynolds number a clear trend towards a larger separated region is evident; at least for the studied, comparably low Reynolds-number regime. It is further shown that the small separated region occurring at the diffuser throat shows the opposite behaviour as the main separation region, i.e. the flow is separating less with higher Re _b . Moreover, the influence of the Reynolds number on the internal layer occurring at the non-inclined wall described in a recent study has also been assessed. It can be concluded that this region close to the upper, straight wall, is more distinct for larger Re _b . Additionally, the influence of temporal correlations arising from the commonly used periodic turbulent channel flow as inflow condition (similar to a precursor simulation) for the diffuser is assessed.     

Heat transfer enhancement in grooved channels due to flow bifurcations

The flow bifurcation scenario and heat transfer characteristics in grooved channels, are investigated by direct numerical simulations of the mass, momentum and energy equations, using the spectral element methods for increasing Reynolds numbers in the laminar and transitional regimes. The Eulerian flow characteristics show a transition scenario of two Hopf bifurcations when the flow evolves from a laminar to a time-dependent periodic and then to a quasi-periodic flow. The first Hopf bifurcation occurs to a critical Reynolds number Rec₁ that is significantly lower than the critical Reynolds number for a plane-channel flow. The periodic and quasi-periodic flows are characterized by fundamental frequencies ω₁ and m· ω₁+n·ω₂, respectively, with m and n integers. Friction factor and pumping power evaluations demonstrate that these parameters are much higher than the plane channel values. The time-average mean Nusselt number remains mostly constant in the laminar regime and continuously increases in the transitional regime. The rate of increase of this Nusselt number is higher for a quasi-periodic than for a periodic flow regime. This higher rate originates because better flow mixing develops in quasi-periodic flow regimes. The flow bifurcation scenario occurring in grooved channels is similar to the Ruelle-Takens-Newhouse transition scenario of Eulerian chaos, observed in symmetric and asymmetric wavy channels.     

Numerical simulations of laminar flow over a 3D backward-facing step

A numerical investigation of laminar flow over a three-dimensional backward-facing step is presented with comparisons with detailed experimental data, available in the literature, serving to validate the numerical results. The continuity constraint method, implemented via a finite element weak statement, was employed to solve the unsteady three-dimensional Navier–Stokes equations for incompressible laminar isothermal flow. Two-dimensional numerical simulations of this step geometry underestimate the experimentally determined extent of the primary separation region for Reynolds numbers Re greater than 400. It has been postulated that this disagreement between physical and computational experiments is due to the onset of three-dimensional flow near Re ≈ 400. This paper presents a full three-dimensional simulation of the step geometry for 100⩽ Re⩽ 800 and correctly predicts the primary reattachment lengths, thus confirming the influence of three-dimensionality. Previous numerical studies have discussed possible instability modes which could induce a sudden onset of three-dimensional flow at certain critical Reynolds numbers. The current study explores the influence of the sidewall on the development of three-dimensional flow for Re greater than 400. Of particular interest is the characterization of three-dimensional vortices in the primary separation region immediately downstream of the step. The complex interaction of a wall jet, located at the step plane near the sidewall, with the mainstream flow reveals a mechanism for the increasing penetration (with increasing Reynolds number) of three-dimensional flow structures into a region of essentially two-dimensional flow near the midplane of the channel. The character and extent of the sidewall-induced flow are investigated for 100⩽Re⩽ 800. © 1997 John Wiley & Sons, Ltd.     

Direct numerical simulation and global stability analysis of three-dimensional instabilities in a lid-driven cavity

The first bifurcation in a lid-driven cavity characterized by three-dimensional Taylor–Görtler-Like instabilities is investigated for a cubical cavity with spanwise periodic boundary conditions at Re=1000. The modes predicted by a global linear stability analysis are compared to the results of a direct numerical simulation. The amplification rate, and the shape of the three-dimensional perturbation fields from the direct numerical simulation are in very good agreement with the characteristics of the steady S1 mode from the stability analysis, showing that this mode dominates the other unstable unsteady modes. To cite this article: J. Chicheportiche et al., C. R. Mecanique 336 (2008).     

Finite element analysis of vortex shedding using equal order interpolations

An operator splitting and element‐by‐element conjugated gradient solver, and equal order interpolations are applied for solving time dependent Navier–Stokes (NS) equations to simulate flow induced vortex shedding in the present study. In addition, the convection term is corrected by balanced tensor diffusivity, which can stabilize the numerical simulation and overcome the numerical oscillations. The evolution of the interested flowing properties with time is analyzed by using spectral analysis. The developed code has been validated by the application of two examples: a driven cavity flow and a flow induced vortex vibration. Results from the first example for Reynolds number Re=10³ and Re=10⁴ are compared with other numerical simulations. Results from the second example, uniform flow past a square rod over a wide range of high Reynolds numbers from Re=10³～10⁵, are compared with experimental data and other numerical studies. Copyright © 2002 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     

Spatial structure of the flow through an axisymmetric sudden expansion

Spatio-temporal velocity fields of an axisymmetric sudden expansion were measured using an ultrasonic velocity profiler and analyzed to investigate the transitional scheme of the spatial structure using two-dimensional Fourier transform and proper orthogonal decomposition techniques. The variation of the zero-crossing point, the fluctuation energy directed upstream and the eigenmode spectrum all have the same transitional scheme as a function of the Reynolds number. The transitional scheme can be classified Re _d<1,000 for the laminar regime, Re _d=1,000–3,000 for the transitional regime and Re _d>3,000 for the turbulent regime. Especially, in the transitional regime, we found large changes in the flow structure at Re _d=1,500 and 2,000. The jump at Re _d=2,000 is caused by the change in the flow condition upstream. The jump at Re _d=1,500 clearly shows a change in the spatial structure of the flow.     

The wake of falling disks at low Reynolds numbers

We visualized the wake structure of circular disks falling vertically in quiescent water.The evolution of the wake was shown to be similar to the flow patterns behind a fixed disk.The Reynolds number,Re = Ud/ν,is in the range of 40 200.With the ascension of Reynolds numbers,a regular bifurcation occurred at the first critical Reynolds number Re c 1,leading to a transition from an axisymmetric wake structure to a plane symmetric one;A Hopf bifurcation took place at the second critical Reynolds number Re c 2,as the wake structure became unsteady.Plane symmetry of the wake structure was first lost as periodic vortex shedding appeared,but recovered at higher Reynolds number.The difference between the two critical Reynolds numbers was found to be shape-dependent,as we compared our results for thin discs with those for other falling bodies,such as spheres and cones.This observation could be understood in terms of the instability mechanism of the vortical structure.     

Numerical Simulation of Aerodynamic Characteristics of Two Rectangular Cylinders in Side-by-side Arrangement

We present numerical results of flows around two rectangular cylinders in side-by-side arrangement by three-dimensional computations. The two rectangular cylinders are arranged with various distances between the cylinders. The three-dimensional flow structures around two rectangular cylinders denote significant features depending on the distance between the cylinders. In our computations, the Reynolds number Re is set as 10,000, and both the aspect ratio of section of two rectangular cylinders and the distance between the cylinders are considered as parameters to calculate the flow around two rectangular cylinders. The obtained numerical results are also compared with experimental data.     

NUMERICAL MODELLING AND PARTICLE IMAGE VELOCIMETRY MEASUREMENT OF THE LAMINAR FLOW FIELD INDUCED BY AN ENCLOSED ROTATING DISC

The fluid flow field within an enclosed cylindrical chamber with a rotating flat disc was calculated using a finite volume computational fluid dynamics (CFD) model and compared with particle image velocimetry (PIV) measurements. Two particular laminar cases near the Transitional flow regime were investigated: Reynolds number Re=2.5×1 ⁴, chamber aspect ratio G (h/R_d)=0.2 and Re=4.2×10⁴, G (h/R_d)=0.217. This enabled direct comparison with the numerical and experimental results reported by other researchers. The computational details and some major factors that affect the computed accuracy and convergence speed are also discussed in detail. PIV results containing some 4300 velocity vector points in each of seven planes for each case were obtained from the flow field parallel to the rotating disc. It was found that PIV results could be obtained in planes within the boundary layers as well as the core flow by careful use of a thin laser illumination sheet and correct choice of laser pulse separation. There was close agreement between numerical results, the present PIV measurements and other reported experimental and numerical results.     

Two-dimensional unsteady laminar flow of a power law fluid across a square cylinder

The two-dimensional and unsteady free stream flow of power law fluids past a long square cylinder has been investigated numerically in the range of conditions 60≤Re≤160 and 0.5≤n≤2.0. Over this range of Reynolds numbers, the flow is periodic in time. A semi-explicit finite volume method has been used on a non-uniform collocated grid arrangement to solve the governing equations. The global quantities such as drag coefficients, Strouhal number and the detailed kinematic variables like stream function, vorticity and so on, have been obtained for the above range of conditions. While, over this range of Reynolds number, the flow is known to be periodic in time for Newtonian fluids, a pseudo-periodic flow regime displaying more than one dominant frequency in the lift is observed for shear-thinning fluids. This seems to occur at Reynolds numbers of 120 and 140 for n=0.5 and 0.6, respectively. Broadly speaking, the smaller the value of the power law index, lower is the Reynolds number of the onset of the pseudo-periodic regime. This work is concerned only with the fully periodic regime and, therefore, the range of Reynolds numbers studied varies with the value of the power law index. Not withstanding this aspect, in particular here, the effects of Reynolds number and of the power law index have been elucidated in the unsteady laminar flow regime. The leading edge separation in shear-thinning fluids produces an increase in drag values with the increasing Reynolds number, while shear-thickening fluid behaviour delays this separation and shows the lowering of the drag coefficient with the Reynolds number. Also, the preliminary results suggest the transition from the steady to unsteady flow conditions to occur at lower Reynolds numbers in shear-thinning fluids than that in Newtonian fluids.     

Theoretical analysis of three-dimensional bifurcated flow inside a diagonally lid-driven cavity

The instability mechanism of fully three-dimensional, highly separated, shear-driven confined flow inside a diagonally lid-driven cavity was investigated. The analysis was conducted on 100³ and 200³ stretched grids by a series of direct numerical simulations utilizing a standard second-order accuracy finite volume code, openFoam. The observed oscillatory instability was found to set in via a subcritical symmetry breaking Hopf bifurcation. Critical values of the Reynolds number Re _cr = 2320 and the non-dimensional angular oscillating frequency \({\omega_{\rm cr}=0.249}\) for the transition from steady to oscillatory flow were accurately determined. An oscillatory regime of the bifurcated flow was analyzed in depth, revealing and characterizing the spontaneous symmetry breaking mechanism. Characteristic spatial patterns of the base flow and the main flow harmonic were determined for the velocity, vorticity and helicity fields. Lagrangian particle tracers were utilized to visualize the mixing phenomenon of the flow from both sides of the diagonal symmetry plane.     

Effect of acceleration on the wavy Taylor vortex flow

In this paper, we use a laser optical technique to investigate the characteristics of a wavy Taylor vortex flow between two concentric cylinders, with the inner cylinder subjected to a wide range of predetermined acceleration and the outer one at rest. We focus on the inner/outer radius ratio of 0.894, with an acceleration (dRe/dt*) from 0.1123 to 2,247, and Reynolds number from Re/Re _c =1.0 to 36. The results show that, with increasing Reynolds number, there is an initial increase in the wavelength of the wavy vortex flow (λ), and a decrease in the wave speed (c) before they asymptote to a constant value, which is a function of the acceleration. As for the wave amplitude (A), it is found that the effect of acceleration is significant only in a very narrow range of Reynolds numbers. Received: 21 August 2001 / Accepted: 22 November 2001     

Symmetry Breaking and Turbulence in Perturbed Plane Couette Flow

Perturbed plane Couette flow containing a thin spanwise-oriented ribbon undergoes a subcritical bifurcation at Re≈230 to a steady three-dimensional state containing streamwise vortices. This bifurcation is followed by several others giving rise to a fascinating series of stable and unstable steady states of different symmetries and wavelengths. First, the backwards-bifurcating branch reverses direction and becomes stable near Re≈200. Then the spanwise reflection symmetry is broken, leading to two asymmetric branches which are themselves destabilized at Re≈420. Above this Reynolds number, time evolution leads first to a metastable state whose spanwise wavelength is halved and then to complicated time-dependent behavior. These features are in agreement with experiments. Received 15 December 2001 and accepted 29 March 2002 Published online: 2 October 2002 Communicated by H.J.S. Fernando     

Transient Growth in Exactly Counter-Rotating Couette–Taylor Flow

Transient growth due to non-normality is investigated for the Couette-Taylor problem with counter-rotating cylinders as a function of aspect ratio η and Reynolds number Re. For all Re≤500, transient growth is enhanced by curvature, i.e. is greater for η<1 than for η=1, the plane Couette limit. For fixed Re>130, it is found that the greatest transient growth is achieved for η on the linear stability boundary. Transient growth is approximately 20% higher near the Couette-Taylor linear stability boundary at Re=310, η=0.986 than at Re=310, η=1, near the threshold observed for transition in plane Couette flow. For 106<Re<130, the greatest transient growth occurs for a value of η between the linear stability boundary and one. For Re<106, the flow is linearly stable and the greatest transient growth occurs for a value of η less than one. The energy in the optimal inputs is primarily meridional; that in the optimal outputs is primarily azimuthal. Pseudospectra are calculated for two contrasting cases. Received 5 November 2001 and accepted 29 March 2002 Published online: 2 October 2002 Communicated by H.J.S. Fernando     

Numerical study of buoyancy- and thermocapillary-driven flows in a cavity

Thermocapillary- and buoyancy-driven convection in open cavities with differentially heated endwalls is investigated by numerical solutions of the two-dimensional Navier-Stokes equations coupled with the energy equation. We studied the thermocapillary and buoyancy convection in the cavities, filled with low-Prandtl-number fluids, with two aspect-ratiosA=1 and 4, Grashof number up to 10⁵ and Reynolds number ⋎Re⋎≤10⁴. Our results show that thermocapillary can have a quite significant effect on the stability of a primarily buoyancy-driven flow, as well as on the flow structures and dynamic behavior for both additive effect (i.e., positiveRe) and opposing effect (i.e., negativeRe).     

Three-dimensional flow in a lid-driven cavity with width-to-height ratio of 1.6

The flow in a lid-driven cavity with width-to-height ratio of 1.6 is investigated numerically and experimentally. Experimental investigation use an apparatus with a spanwise length-to-height ratio of $\Uplambda = 10.85.$ Λ = 10.85 . Increasing the Reynolds number, we experimentally find a gradual change from the quasi-two-dimensional basic flow to a three-dimensional flow pattern. The three-dimensional flow has a significant amplitude at considerably low Reynolds numbers. Streak-line photographs and PIV vector maps are presented to illustrate the structure of the finite-amplitude flow pattern. The smooth transition is in contrast to the linear instability predicted by a linear-stability analysis for a cavity with infinite span. LDV measurements confirm the absence of a distinct threshold Reynolds number and indicate an imperfect bifurcation. The deviations between experimental observations and numerical critical Reynolds number for infinite span are explained by conducting three-dimensional simulations for a finite-span geometry. A good agreement between experimental and numerical simulation is obtained. The numerical and experimental data lead to the conjecture of a premature onset of the three-dimensional flow caused by strong secondary flows which are induced by the cavity end walls. Nevertheless, the flow structure in the finite-span cavity carries the same characteristic signatures as the nonlinear flow in the corresponding infinite-length cavity. We conclude that the observed flow can be identified as the continuation of the normal mode C _e ⁴ earlier identified in a linear-stability analysis.