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.     

Uniform flow of viscoelastic fluids past a confined falling cylinder

Uniform steady flow of viscoelastic fluids past a cylinder placed between two moving parallel plates is investigated numerically with a finite-volume method. This configuration is equivalent to the steady settling of a cylinder in a viscoelastic fluid, and here, a 50% blockage ratio is considered. Five constitutive models are employed (UCM, Oldroyd-B, FENE-CR, PTT and Giesekus) to assess the effect of rheological properties on the flow kinematics and wake patterns. Simulations were carried out under creeping flow conditions, using very fine meshes, especially in the wake of the cylinder where large normal stresses are observed at high Deborah numbers. Some of the results are compared with numerical data from the literature, mainly in terms of a drag coefficient, and significant discrepancies are found, especially for the constant-viscosity constitutive models. Accurate solutions could be obtained up to maximum Deborah numbers clearly in excess of those reported in the literature, especially with the PTT and FENE-CR models. The existence or not of a negative wake is identified for each set of model parameters.     

Very slow flow of Bingham viscoplastic fluid around a circular cylinder

Numerical simulations have been used to study the flow of a Bingham viscoplastic fluid around a circular cylinder in an infinite medium with negligible inertia effects. Papanastasiou's regularisation technique has been adopted to approximate the model. The case corresponding to preponderant plasticity effects has been particularly studied and convergence of the solutions examined in detail. The flow kinematics and stresses have been determined. The rigid zones have been identified and characterised. At large Oldroyd numbers, when plasticity effects become preponderant, a viscoplastic boundary layer appears around the cylinder. The characteristics of this viscoplastic boundary layer are quantified. The results are compared with existing theoretical results, concerning particularly the predictions of the viscoplastic boundary layer theory and the plasticity theory.     

Experimental study on the flow regimes past a confined prism undergoing self-sustained oscillations

An experimental study based on Particle Image Velocimetry (PIV) is presented with the objective of studying the flow regimes that appear in the flow past a confined prism undergoing self-sustained oscillations at low Reynolds numbers (Re). The square-section prism, placed inside a 3D square cross-section vertical channel with a confinement ratio of 1/2.5, was tethered to the channel walls and, therefore, it was allowed to move freely transverse to the incoming flow. Re (based on the prism cross-section height) was varied in the range from 100 to 700. Three different prism to fluid density ratios (m¹) were considered: 0.56, 0.70, and 0.91. These two parameters, Re and m¹, were used to map the results obtained. In particular, it was found that five different regimes appear: (1) steady prism with steady recirculation bubble, (2) steady prism with unsteady vortex shedding wake, (3) large amplitude low frequency oscillating prism with unsteady vortex shedding wake, (4) small amplitude high frequency oscillating prism with unsteady vortex shedding wake, and (5) irregular/chaotic motion of both the prism and the wake. The PIV results and associated numerical simulations were used to analyze the different prism and wake states.     

Series solutions of annular axisymmetric stagnation flow and heat transfer on moving cylinder

The steady, laminar, incompressible flow and heat transfer of a viscous fluid between two circular cylinders for two different types of thermal boundary conditions are investigated. The governing Navier-Stokes and thermal equations of the flow are reduced to a nonlinear system of ordinary differential equations. The equations are solved analyt- ically using the homotopy analysis method （HAM）. Convergence of the HAM solutions is discussed in detail. These solutions are then compared with recently obtained numericM and perturbative solutions. Plots of the velocity and temperature profiles are provided for various values of the relevant parameters.     

New contributions on laminar flow of inelastic non-Newtonian fluid in the two-dimensional symmetric expansion: Creeping and slowly moving flow conditions

The steady flow of generalized Newtonian fluid in a two-dimensional 1:3 sudden expansion was studied numerically. Finite volume method was applied to solve the momentum equations along with the continuity equation and the Power law rheological model within the laminar flow regime for a range of Reynolds number and Power law index values. The values of generalized Reynolds number, based on physical and rheological properties, upstream channel height and bulk velocity, were varied between 0.0001 ≤ Re_gen ≤ 10, while the Power law index values mapped the 0.60 ≤ n ≤ 1.40 range, allowing for the investigation of both shear-thinning and shear-thickening effects at creeping as well as slowly moving fluid flow conditions. We report accurate results of a systematic study with a focus on most important characteristics of recirculating fluid flow in the downstream section of sudden expansion geometry. It is shown that for the creeping flow regime there exist finite sized redevelopment length, extra pressure drop (Couette correction) and recirculation zones (also called as Moffatt vortices) that are influenced by the non-Newtonian viscous behaviour.     

A note on the numerical simulations of flow past a wavy square-section cylinder

The flow past a square-section cylinder with a geometric disturbance is investigated by numerical simulations. The extra terms, due to the introduction of mapping transformation simulating the effect of disturbance into the transformed Navier-Stokes equations, are correctly derived, and the incorrect ones in the previous literature are pointed out and analyzed. Furthermore, the relationship between the vorticity, especially on the cylinder surface, and the disturbance is derived and explained theoretically. The computations are performed at two Reynolds numbers of 100 and 180 and three amplitudes of waviness of 0.006, 0.025 and 0.167 with another aim to explore the effects of different Reynolds numbers and disturbance on the vortex dynamics in the wake and forces on the body. Numerical results have shown that, at the mild waviness of 0.025, the Kairmain vortex shedding is suppressed completely for Re = 100, while the forced vortex dislocation is appeared in the near wake at the Reynolds number of 180. The drag reduction is up to 21.6% at Re = 100 and 25.7% at Re = 180 for the high waviness of 0.167 compared with the non-wavy cylinder. The lift and the Strouhal number varied with different Reynolds numbers and the wave steepness are also obtained.     

An immersed interface method for flow past circular cylinder in the vicinity of a plane moving wall

Two‐dimensional flows past a stationary circular cylinder near a plane boundary are numerically simulated using an immersed interface method with second‐order accuracy. Instead of a fixed wall, a moving wall with no‐slip boundary is considered to avoid the complex involvement of the boundary layer and to focus only on the shear‐free wall proximity effects for investigating the force dynamics and flow fields. To analyze the convergence and accuracy of our implementation, numerical studies have been first performed on a simple test problem of rotational flow, where the second order of convergence is confirmed through numerical experiments and an optimal range of relative grid‐match ratio of Lagrangian to Eulerian grid sizes has been recommended. By comparing the force quantities and the Strouhal number, the accuracy of this method has been demonstrated on the flow past a stationary isolated cylinder. The cylinder is then put in proximity to the wall to investigate the shear‐free wall proximity effects in the low Reynolds number regime (20≤Re≤200). The gap ratio, e/D, where e denotes the gap between the cylinder and the moving wall and D denotes the diameter of the cylinder, is taken from 0.10 to 2.00 to determine the critical gap ratio, (e/D)_critical, for the alternate vortex shedding, where the fluid forces, flow fields and the streamwise velocity profiles are studied. One of the key findings is that the (e/D)_critical for the alternate vortex shedding decreases as the Reynolds number increases. We also find that, in this low Reynolds number regime, the mean drag coefficient increases and peaks at e/D = 0.5 with the increase of e/D and keeps decreasing gently from e/D = 0.5 to e/D = 2.0, while the mean lift coefficient decreases monotonically with the increase of e/D. New correlations are then proposed for computing force coefficients as a function of Re and e/D for a cylinder in the vicinity of a moving plane wall. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical solutions for impulsively started and decelerated viscous flow past a circular cylinder

A finite difference study of the unsteady two-dimensional flow past a circular cylinder has been conducted using vorticity and streamfunction as the dependent variables. The two cases considered were impulsively started and decelerated flows. The impulsively started problem was considered to validate the method and has yielded results which agree quite closely with existing results from both calculations and experiments. The decelerated flow analysis produced results which can be explained in terms of induced velocity effects from existing wake vortices for both suddenly stopped and uniformly decelerated flows.     

Stability of flow past a cylinder: Energy budget of eigenmodes

Global linear stability analysis of the flow past a circular cylinder at the onset of primary wake instability is carried out. The real and imaginary parts of the most unstable eigenmode, responsible for vortex shedding, are very similar but associated with a spatial shift in the vortex structures. This shift results in the convection of vortices that are observed in the unsteady flow, which is actually a consequence of global absolute instability. The kinetic energy density, associated with the most unstable eigenmode, is studied. At the onset of the instability the energy density of the disturbance field is found to be stronger in the far wake compared with the near wake. With increase in Re the region where the disturbance is strong moves upstream closer to the cylinder. However, the maximum value of the kinetic energy density of the disturbance lies outside the recirculation zone even for Re upto 100. A linearized mechanical energy equation for the time evolution of the kinetic energy density of the disturbance is utilized to examine the energy budget of the most unstable eigenmode at various Re. It is found that the most significant contribution to the growth rate of the disturbance arises from the transfer of the energy due to the strain rate of the base flow to the perturbation. The stabilizing effect of the viscous dissipation increases with increase in Re, but saturates for Re beyond ～70. Copyright © 2009 John Wiley & Sons, Ltd.     

Stress boundary layers in the viscoelastic flow past a cylinder in a channel： limiting solutions

The flow past a cylinder in a channel with the aspect ratio of 2:1 for the upper convected Maxwell (UCM) fluid and the Oldroyd-B fluid with the viscosity ratio of 0.59 is studied by using the Galerkin/Least-square finite element method and a p-adaptive refinement algorithm. A posteriori error estimation indicates that the stress-gradient error dominates the total error. As the Deborah number, D_e, approaches 0.8 for the UCM fluid and 0.9 for the Oldroyd-B fluid, strong stress boundary layers near the rear stagnation point are forming, which are characterized by jumps of the stress-profiles on the cylinder wall and plane of symmetry, huge stress gradients and rapid decay of the gradients across narrow thicknesses. The origin of the huge stress-gradients can be traced to the purely elongational flow behind the rear stagnation point, where the position at which the elongation rate is of 1/2D_e approaches the rear stagnation point as the Deborah number approaches the critical values. These observations imply that the cylinder problem for the UCM and Oldroyd-B fluids may have physical limiting Deborah numbers of 0.8 and 0.9, respectively.The project supported by the National Natural Science Foundation of China (50335010 and 20274041) and the MOLDFLOW Comp. Australia.     

An experimental study of the dynamic characteristics in boundary layer in the flow past a two-dimensional circular cylinder at subcritical reynolds numbers

An experimental study of the dynamic characteristics of flow past a two-dimensional circular cylinder is described. The fluctuationsoof wall shear stress, surface-pressure and velocity of the flow are measured with hot-film, hot-wire and pressure transducer. The frequency feature of fluctuations of wall shear stress is given. The cross-correlation functions of these fluctuations at any two points are calculated. The experimental results reveal that there is an overall syncronous fluctuation, at the shedding frequency, in boundary layer in the flow past a two-dimensional circular cylinder at subcritical Reynolds number.     

Numerical study on the dynamics of droplet passing through a cylinder obstruction in confined microchannel flow

The droplet dynamics passing through a cylinder obstruction was investigated with direct numerical simulations with FE-FTM (Finite Element-Front Tracking Method). The effect of droplet size and capillary number (Ca) was studied for both Newtonian and viscoelastic fluids. In the case of Newtonian droplet immersed in Newtonian medium, the droplet breakup induced by the geometric hindrance depends on the droplet size. As Ca increases, the short droplets (1.3 times longer than the channel width) break up while passing through the obstruction. However, the breakup does not occur for longer droplets (1.8 times longer than the channel width). When the viscoelastic fluid characterized by the Oldroyd-B model is considered, the Newtonian droplet immersed in viscoelastic medium breaks up into two smaller droplets while passing through the cylinder obstruction with increasing De_m (Deborah number of the medium). We also show that the normal stress difference plays a key role on the droplet breakup and the droplet extension. The normal stress difference is enhanced in the negative wake region due to the droplet flow, which also promotes droplet extension in that region. This numerical study provides information not only on underlying physics of the droplet flows passing through a cylinder obstruction but also on the useful guidelines for microfluidic applications.     

LES of incompressible heat and fluid flows past a square cylinder at high Reynolds numbers

This paper presents large eddy simulation (LES) results of incompressible heat and fluid flows around a square cylinder (SC) at zero incident angle at high Reynolds numbers (Re) in the range from 1.25×10⁵ to 3.5×10⁵. LES results are obtained on the basis of swirling strength based sub-grid model, and a higher order upwind scheme developed with respect to the Taylor expansion. It was found that, for the zero incident SC wake flows at a Reynolds number in the range {Re₅ = Re/10⁵ ∈ [1.25, 3.5]}, the Strouhal number equals to 0.1079, completely independent of the Reynolds number; the coefficient of drag is around 1.835 with an uncertainty of about 1.9%, almost non-sensitive to the Re. When Re is beyond 3.0×10⁵, the time-averaged peak value of sub-grid viscosity is over 340, implying that the role of sub-grid model is crucial in some regions where vortex motion is active and vortex interaction is intense. The time–spanwise (t-z) averaged sub-grid viscosity ratio profiles and the profiles of fluctuations of the sub-grid viscosity ratio and velocity components at four locations downstream of the SC are presented. The fields of the t-z averaged sub-grid viscosity ratio, and the instantaneous fields of streamwise and spanwise vorticities are also reported and discussed. The predicted mean Nusselt number is compared with empirical correlations, revealing that swirling strength based LES has its potential in predicting natural and industrial flows.     

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.     

Unsteady characteristics of the cavity flow around a rough circular cylinder

Effects of surface roughness on the unsteady cavitating flow around a two-dimensional circular cylinder were experimentally investigated at Reynolds numbers from 1.36 × 10⁵ to 1.78 × 10⁵. Two patterns of surface roughness were investigated, a double-cut pattern and a single-cut pattern. The cavity elongates with an increase of the surface roughness, especially in supercavitating flow. However, for some roughness parameters tested, the cavity length exhibits an extreme decrease. In a particular case of the double-cut pattern, there exists the minimum cavity behind the cylinder.     

Influence of drag laws on pressure and bed material recirculation rate in a cold flow model of an 8 MW dual fluidized bed system by means of CPFD

A cold flow model of an 8 MW dual fluidized bed (DFB) system is simulated using the commercial computational particle fluid dynamics (CPFD) software package Barracuda. The DFB system comprises a bubbling bed connected to a fast fluidized bed with the bed material circulating between them. As the hydrodynamics in hot DFB plants are complex because of high temperatures and many chemical reaction processes, cold flow models are used. Performing numerical simulations of cold flows enables a focus on the hydrodynamics as the chemistry and heat and mass transfer processes can be put aside. The drag law has a major influence on the hydrodynamics, and therefore its influence on pressure, particle distribution, and bed material recirculation rate is calculated using Barracuda and its results are compared with experimental results. The drag laws used were energy-minimization multiscale (EMMS), Ganser, Turton–Levenspiel, and a combination of Wen–Yu/Ergun. Eleven operating points were chosen for that study and each was calculated with the aforementioned drag laws. The EMMS drag law best predicted the pressure and distribution of the bed material in the different parts of the DFB system. For predicting the bed material recirculation rate, the Ganser drag law showed the best results. However, the drag laws often were not able to predict the experimentally found trends of the bed material recirculation rate. Indeed, the drag law significantly influences the hydrodynamic outcomes in a DFB system and must be chosen carefully to obtain meaningful simulation results. More research may enable recommendations as to which drag law is useful in simulations of a DFB system with CPFD.