A numerical study of steady and unsteady viscoelastic flow past bounded cylinders

We consider two-dimensional, inertia-free, flow of a constant-viscosity viscoelastic fluid obeying the FENE-CR equation past a cylinder placed symmetrically in a channel, with a blockage ratio of 0.5. Through numerical simulations we show that the flow becomes unsteady when the Deborah number (using the usual definition) is greater than De ≈ 1.3, for an extensibility parameter of the model of L² = 144. The transition from steady to unsteady flow is characterised by a small pulsating recirculation zone of size approximately equal to 0.15 cylinder radius attached to the downstream face of the cylinder. There is also a rise in drag coefficient, which shows a sinusoidal variation with time. The results suggest a possible triggering mechanism leading to the steady three-dimensional Gortler-type vortical structures, which have been observed in experiments of the flow of a viscoelastic fluid around cylinders. The results reveal that the reason for failure of the search for steady numerical solutions at relatively high Deborah numbers is that the two-dimensional flow separates and eventually becomes unsteady. For a lower extensibility parameter, L² = 100, a similar recirculation is formed given rise to a small standing eddy behind the cylinder which becomes unsteady and pulsates in time for Deborah numbers larger than De ≈ 4.0–4.5.     

The flow field around a micropillar confined in a microchannel

The flow field over a low aspect ratio (AR) circular pillar (L/D = 1.5) in a microchannel was studied experimentally. Microparticle image velocimetry (μPIV) was employed to quantify flow parameters such as flow field, spanwise vorticity, and turbulent kinetic energy (TKE) in the microchannel. Flow regimes of cylinder-diameter-based Reynolds number at 100 ⩽ Re_D ⩽ 700 (i.e., steady, transition from quasi-steady to unsteady, and unsteady flow) were elucidated at the microscale. In addition, active flow control (AFC), via a steady control jet (issued from the pillar itself in the downstream direction), was implemented to induce favorable disturbances to the flow in order to alter the flow field, promote turbulence, and increase mixing. Together with passive flow control (i.e., a circular pillar), turbulent kinetic energy was significantly increased in a controllable manner throughout the flow field.     

Influence of aspect ratio on the flow above the free end of a surface-mounted finite cylinder

The flow above the free end of a surface-mounted finite-height cylinder was studied in a low-speed wind tunnel using particle image velocimetry (PIV). Velocity measurements were made in vertical and horizontal measurement planes above the free end of finite cylinders of aspect ratios AR = 9, 7, 5 and 3, at a Reynolds number of Re = 4.2 × 10⁴. The relative thickness of the boundary layer on the ground plane was δ/D = 1.7. Flow separating from the leading edge formed a prominent recirculation zone on the free-end surface. The legs of the mean arch vortex contained within the recirculation zone terminate on the free-end surface on either side of the centreline. Separated flow from the leading edge attaches onto the upper surface of the cylinder along a prominent attachment line. Local separation downstream of the leading edge is also induced by the reverse flow and arch vortex circulation within the recirculation zone. As the cylinder aspect ratio is lowered from AR = 9 to AR = 3, the thickness of the recirculation zone increases, the arch vortex centre moves downstream and higher above the free-end surface, the attachment position moves downstream, and the termination points of the arch vortex move upstream. A lowering of the aspect ratio therefore results in accentuated curvature of the arch vortex line. Changes in aspect ratio also influence the vorticity generation in the near-wake region and the shape of the attachment line.     

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.     

Static wall layers in plane channel displacement flows

We study steady and pulsating displacement flows of a Bingham fluid by a Newtonian fluid, along a plane channel. For sufficiently large yield stress a static residual wall layer can result during the displacement. The flow is parameterised by the Reynolds number (Re), the Bingham number (B) and the viscosity ratio (M). Perhaps intuitively, thicker layers are found with larger M and at lower Re. The residual layer is formed on the advective timescale of the displacement but drains on a slower timescale governed by M. For larger M truly stationary layers are only found for large t when the layer has thinned sufficiently to become static. Increased Re results in increased energy production locally around the finger. For large enough Re the energy production can play a significant role in yielding the fluid. As the energy production rate increases it also becomes focused around the corner or shoulder region of the front, and spreads axially along the initial part of the residual layer. This causes fluid to yield increasingly far behind the front and allows for the layer to thin. As B increases the static layer tends to decrease (see also [1], [2]). At small Re the static layer thickness appears to be independent of M. At large Re the layer thickness is dependent on M and decreases asymptotically to a constant value as B → ∞.For pulsating displacement flow rates, Q(t) = 2(1 + Asin  ωt) : A ∈ [0, 1) we study two ranges: ωRe ? 2π and ωRe ? 2π. For the viscous regime (ωRe ? 2π) a pseudo-steady 1D model predicts that the residual layer should remain static for 3(1 + Asin  ωt) < MB. In practice we find that partial mobilisation of the residual layer occurs even when this inequality is satisfied, but not if MB becomes significantly larger than 3(1 + A). For ωRe ? 2π we mobilise the layer for significantly larger values of MB and at smaller A, than in the viscous regime. This effect is traced to the occurrence of out-of-phase velocity fluctuations in the displacing fluid within a wall layer close to the interface.     

Flow structures around an equilateral triangle arrangement of three spheres

This paper represents the results of an experimental study on the flow structure around a single sphere and three spheres in an equilateral-triangular arrangement. Flow field measurements were performed using a Particle Image Velocimetry (PIV) technique and dye visualization in an open water channel for a Reynolds number of Re = 5 × 10³ based on the sphere diameter. The distributions and flow features at the critical locations of the contours of the velocity fluctuations, the patterns of sectional streamlines, the vorticity contours, the turbulent kinetic energy, the Reynolds stress correlations and shedding frequency are discussed. The gap ratios (G/D) of the three spheres were varied in the range of 1.0 ⩽ G/D ⩽ 2.5 where G was the distance between the sphere centers, and D was the sphere diameter which was taken as 30 mm. Due to the interference of the shedding shear layers and the wakes, more complex features of the flow patterns can be found in the wake region of the two downstream spheres behind the leading sphere. For G/D = 1.25, a jet-like flow around the leading sphere through the gap between the two downstream spheres occurred, which significantly enhanced the wake region. It was observed that a continuous flow development involving shearing phenomena and the interactions of shedding vortices caused a high rate of fluctuations over the whole flow field although most of the time-averaged flow patterns were almost symmetric about the two downstream spheres.     

Two-phase flow behavior inside a header connected to multiple parallel channels

The main objective of this work is to examine the flow distribution of two-phase mixture to parallel channels and to investigate the flow behavior at header-channel junctions simulating the corresponding parts of compact heat exchangers. The cross-section of the header and the channels were fixed to 14 mm × 14 mm and 12 mm × 1.6 mm, respectively. The mass flux and the mass quality ranges were 70–165 kg/m² s and 0.3–0.7, respectively. Air and water were used as the test fluids. The flow distribution at the fore part of the header (region A) is affected only by the upstream flow configuration and the rate of liquid flow separation decreased a flowing downwards. On the other hand, in the rear part, the downstream effect predominates over the upstream effect due to strong flow recirculation near the end plate. In this part, the liquid separation increased (region B) and then decreased (region C) as the mixture proceeds downwards. The validity of the existing models for branching flows at parallel T-junction was tested, and turned out to be appropriate for region A. However, the models were not applicable to the rear part due to a strong flow recirculation. Moreover, the effect of the membranes in channels was investigated, but that was minor.     

Experimental investigation of a developing two-phase bubbly flow in horizontal pipe

Experimental results for various water and air superficial velocities in developing adiabatic horizontal two-phase pipe flow are presented. Flow pattern maps derived from videos exhibit a new boundary line in intermittent regime. This transition from water dominant to water–gas coordinated regimes corresponds to a new transition criterion C_T = 2, derived from a generalized representation with the dimensionless coordinates of Taitel and Dukler.Velocity, turbulent kinetic energy and dissipation rate, void fraction and bubble size radial profiles measured at 40 pipe diameters for J_L = 4.42 m/s by hot film velocimetry and optical probes confirm this transition: the gas influence is not continuous but strongly increases beyond J_G = 0.06 m/s. The maximum dissipation rate, derived from spectra, is increased in two-phase flow by a factor 5 with respect to the single phase case.The axial evolution of the bubble intercept length histograms also reveal the flow organization in horizontal layers, driven by buoyancy effects. Bubble coalescence is attested by a maximum bubble intercept evolving from 2.5 to 4.5 mm along the pipe. Turbulence generated by the bubbles is also manifest by the 4-fold increase of the maximum turbulent dissipation rate along the pipe.     

Computational and experimental characterization of a liquid jet plunging into a quiescent pool at shallow inclination

A circular water jet (Re = 1.6 × 10⁵; We = 8.8 × 10³) plunging at shallow angles (θ ≈ 12.5°) into a quiescent pool is investigated computationally and experimentally. A surprising finding from the work is that cavities, of the order of jet diameter, are formed periodically in the impact location, even though the impinging flow is smooth and completely devoid of such a periodicity. Computational prediction of these frequencies was compared with experimental findings, yielding excellent agreement. The region in the vicinity of the impact is characterized by strong churning due to splashing and formation of air cavities. Measured velocity profiles indicate a concentration of momentum beneath the free surface slightly beyond the impact location (X/D_j ≈ 14), with a subsequent shift towards the free surface further downstream of this point (X/D_j ≈ 30). This shift is due primarily to the action of buoyancy on the cavity/bubble population. Comparisons of the mean velocity profile between simulations and experiments are performed, yielding good agreement, with the exception of the relatively small churning flow region. Further downstream (X/D_j ≳ 40), the flow develops mostly due to diffusion and the location of peak velocity coincides with the free surface. In this region, the free surface acts as an adiabatic boundary and restricts momentum diffusion, causing the peak velocity to occur at the free surface.     

Low Reynolds number flow over a square cylinder with a detached flat plate

A numerical study of the alteration of a square cylinder wake using a detached downstream thin flat plate is presented. The wake is generated by a uniform flow of Reynolds number 150 based on the side length of the cylinder, D. The sensitivity of the near wake structure to the downstream position of the plate is investigated by varying the gap distance (G) along the wake centerline in the range 0 ⩽ G ⩽ 7D for a constant plate length of L = D. A critical gap distance is observed to occur at G_c ≈ 2.3D that indicates the existence of two flow regimes. Regime I is characterised by vortex formation occurring downstream of the gap while for regime II, formation occurs within the gap. By varying the plate length and gap distance, a condition is found where significant unsteady total lift reduction can occur. The root mean square lift reduction is limited by an unsteady stall process on the plate.     

Effect of side ratio on fluid flow and heat transfer from rectangular cylinders using the PANS method

A computational study of heat transfer from rectangular cylinders is carried out. Rectangular cylinders are distinguished based on the ratio of the length of streamwise face to the height of the cross-stream face (side ratio, R). The simulations were performed to understand the heat transfer in a flow field comprising separation, reattachment, vortex shedding and stagnation. The Partially-Averaged Navier–Stokes (PANS) modeling approach is used to solve the turbulent flow physics associated and the wall resolve approach is used for the near wall treatment because of the flow separation involved. The simulations were performed using a finite volume based opensource software, OpenFOAM, at Reynolds number (Re) = 22,000 for rectangular cylinder at constant temperature kept in an air stream. Two critical side ratios were obtained, R = 0.62 and 3.0. At R = 0.62, the maximum value of the drag coefficient (C_d) = 2.681 was observed which gradually reduced by 54% at R = 4.0. The base pressure coefficient and global Nusselt number also attained the maximum value at R = 0.62 and from R = 2.5 to 3.0 a sharp discontinuous increase by 140% in the Strouhal number was observed. At R = 0.62, it was observed that the separated flow reattaches at the trailing edge after rolling over the side face and therefore increases the overall Nusselt number. The phase averaging was also performed to analyze the unsteady behavior of heat transfer.     

Extending the bounds of ‘steady’ RANS closures: Toward an instability-sensitive Reynolds stress model

The incapability of the conventional Unsteady RANS (Reynolds–Averaged Navier Stokes) models to adequately capture turbulence unsteadiness presents the prime motivation of the present work, which focuses on formulating an instability-sensitive, eddy-resolving turbulence model on the Second-Moment Closure level. The model scheme adopted, functioning as a ‘sub-scale’ model in the Unsteady RANS framework, represents a differential near-wall Reynolds stress model formulated in conjunction with the scale-supplying equation governing the homogeneous part of the inverse turbulent time scale ω_h (ω_h = ɛ_h/k). The latter equation was straightforwardly obtained from the model equation describing the dynamics of the homogeneous part of the total viscous dissipation rate ɛ, defined as ɛ_h = ɛ − 0.5ν∂²k/(∂x_j∂x_j) (Jakirlic and Hanjalic, 2002), by applying the derivation rules to the expression for ω_h. The model capability to account for vortex length and time scales variability was enabled through an additional term in the corresponding length-scale determining equation, providing a selective enhancement of its production, pertinent particularly to the highly unsteady separated shear layer region, modeled in terms of the von Karman length scale (comprising the second derivative of the velocity field) in line with the SAS (Scale-Adaptive Simulation) proposal (Menter and Egorov, 2010). The present model formulation, termed as SRANS model (Sensitized RANS), does not comprise any parameter depending explicitly on grid spacing. The predictive capabilities of the newly proposed length-scale determining model equation, solved in conjunction with Jakirlic and Hanjalic’s (2002) Reynolds stress model equation, are presently demonstrated by computing the flow configurations of increasing complexity featured by boundary layer separation from sharp-edged and continuous curved surfaces: backward-facing step flow, flow over a wall-mounted fence, flow over smoothly contoured periodically arranged hills and flow in a 3-D diffuser. The model performances are also assessed in capturing the natural decay of the homogeneous isotropic turbulence and the near-wall Reynolds stress anisotropy in a plane channel. In most cases considered the fluctuating velocity field was obtained starting from steady RANS results.     

Influence of single rectangular groove on the flow past a circular cylinder

In the present study, flow control mechanism of single groove on a circular cylinder surface is presented experimentally using Particle image velocimetry (PIV). A square shaped groove is patterned longitudinally on the surface of the cylinder with a diameter of 50 mm. The flow characteristics are studied as a function of angular position of the groove from the forward stagnation point of the cylinder within 0° ≤ θ ≤ 150°. In the current work, instantaneous and time-averaged flow data such as vorticity, ω streamline, Ψ streamwise, u/U_o and transverse, v/U_o velocity components, turbulent kinetic energy, TKE and RMS of streamwise, u_rms and transverse, v_rms velocity components are utilized in order to present the results of quantitative analyses. Furthermore, Strouhal numbers are calculated using Karman vortex shedding frequency, f_k obtained from single point spectral analysis. It is concluded that a critical angular position of the groove, θ = 80° is observed. The flow separation is controlled within 0° ≤ θ < 80°. At θ = 80°, the flow separation starts to occur in the upstream direction. The instability within the shear layer is also induced on grooved side of the cylinder with frequencies different than Karman vortex shedding frequency, f_k.     

Flow instability in baffled channel flow

Flow instability in baffled channel flow, where thin baffles are mounted on both channel walls periodically in the direction of the main flow, has been numerically investigated. The geometry considered here can be regarded as a simple model for finned heat exchangers. The aim of this investigation is to understand how baffle interval (L) and Reynolds number (Re) influence the flow instability. With a fixed baffle length of one quarter of channel height (H), ratios of baffle interval to channel height (R_B = L/H) between 1 and 4 are considered. The critical Reynolds number of the primary instability, a Hopf bifurcation from steady flow to time-periodic flow, turned out to be minimum when R_B = 3.08. The friction factor (f) is strongly correlated with the critical Reynolds number for R_B ⩽ 2.5. For the particular cases of R_B = 1.456 and R_B = 1.0, we performed Floquet stability analysis in order to study the secondary instability through which time-periodic two-dimensional flow bifurcates into three-dimensional flow. The results obtained in this investigation are in good agreement with those computed from full simulations, and shed light on understanding and controlling flow characteristics in a finned heat exchanger, quite beneficial to its design.     

A study of pulsating flow in automotive catalyst systems

Conversion efficiency, durability and pressure drop of automotive exhaust catalysts are dependent on the flow distribution within the substrate. This study examines the effect of pulsating flow on the flow distribution within these systems. The flow distribution was measured for a range of flow rates at pulsation frequencies of 16, 32, 64 and 100 Hz. It was shown that the flow uniformity at 16 Hz was similar to the steady equivalent whereas improved uniformity was seen at the higher frequencies resulting in a reduced pressure drop. It was further found that flow maldistribution under pulsating conditions was less sensitive to increases in flow rate compared to steady-state flow. Downstream of the monolith strong pulses were observed although the pulse shapes changed across the substrate diameter. Flow maldistribution correlated well with a non-dimensional parameter derived from the inlet flow velocity, pulsation frequency and diffuser length.     

The effect of sidewalls on rectangular jets

An experimental study is presented regarding the influence of sidewalls on the turbulent free jet flow issuing from a smoothly contracting rectangular nozzle of aspect ratio 15. “Sidewalls” are two parallel plates, flush with each of the slots’ short sides, practically establishing bounding walls extending the nozzle sidewalls in the downstream direction. Measurements of the streamwise and lateral velocity mean and turbulent characteristics have been accomplished, with an x-sensor hot wire anemometer, up to an axial distance of 35 nozzle widths, for jets with identical inlet conditions with and without sidewalls. Centreline measurements for both configurations have been collected for three Reynolds numbers, Re_D = 10,000, 20,000 and 30,000. For Re_D = 20,000 measurements in the transverse direction were collected at 13 different downstream locations in the range, x = 0–35 nozzle widths, and in the spanwise direction at three different downstream locations, x = 2, 6 and 25 nozzle widths.Results indicate that, the two jet configurations (with and without sidewalls) produce statistically different flow fields. Sidewalls do not lead to the production of a 2D flow field as undulations in the spanwise mean velocity distribution indicate. They do increase the two-dimensionality of the jet increasing the longevity of 2D spanwise rollers structures formed in the initial stages of entrainment, which are responsible for the convection of longitudinal momentum towards the outer field, establishing larger streamwise mean velocities at the jet edges. In the near field, up to 25 nozzle widths, lower outward lateral velocities in the presence of the sidewalls are held responsible for the decrease of turbulent terms including rms of velocity fluctuations and Reynolds stresses. Skewness factors increase monotonically across the shear layers from negative values to positive forming sharp peaks at the outer edges of the jet, illustrative of the presence of well defined 2D roller structures in the jet with sidewalls.     

A numerical study of the longitudinal thermoconvective rolls in a mixed convection flow in a horizontal channel with a free surface

This paper presents a numerical study of three-dimensional laminar mixed convection within a liquid flowing on a horizontal channel heated uniformly from below. The upper surface is free and assumed to be flat. The coupled Navier–Stokes and energy equations are solved numerically by the finite volume method taking into account the thermocapillary effects (Marangoni effect). When the strength of the buoyancy, thermocapillary effects and forced convective currents are comparable (Ri ⋍ O(1) and Bd = Ra/Ma ⋍ O(1)), the results show that the development of instabilities in the form of steady longitudinal convective rolls is similar to those encountered in the Poiseuille–Rayleigh–Bénard flow. The number and spatial distribution of these rolls along the channel depend on the flow conditions. The objective of this work is to study the influence of parameters, such as the Reynolds, Rayleigh and Biot numbers, on the flow patterns and heat transfer characteristics. The effects of variations in the surface tension with temperature gradients (Marangoni effect) are also considered.     

Oil emulsion separation with fluidic oscillator generated microbubbles

Surfactants stabilise oil droplets in water, forming a dispersed oil–water emulsion. Treatment of oily effluents is a serious challenge owing to the high stability and colloidal nature of the oil droplets. In many applications, microbubbles are employed for separation purposes due to their buoyancy and increased surface area to volume ratio. This property has been exploited in the water treatment industry for separation in a process known as dissolved air flotation (DAF). Though practically efficient, the process is energy intensive operating at >5 bars and consequently consuming ∼90% of the total energy required in water purification plants. In this study microbubbles were produced by fluidic oscillation via a no-moving part diverter valve to cut down the energy consumption considerably. Microbubbles are applied for the separation of emulsified oil in a process known as microflotation. The mean bubble size generated by fluidic oscillation from the 50 μm pore diffuser was ∼100 μm, otherwise coarse bubbles were produced under steady flow. The effect of surfactant concentration on oil droplet size was investigated. It was found that oil droplet size varied inversely proportional to surfactant concentration. In addition, it was found that the oil removal efficiency also depends on the surfactant concentration. The maximum oil removal efficiency by Microflotation was found to be 91% under lowest surfactant concentration tested (0.3 wt%) whilst at highest surfactant concentration used (10 wt%); lowest recovery efficiency (19.4%) was recorded.     

Laminar flow of power-law fluids past a rotating cylinder

In this work, the continuity and momentum equations have been solved numerically to investigate the flow of power-law fluids over a rotating cylinder. In particular, consideration has been given to the prediction of drag and lift coefficients as functions of the pertinent governing dimensionless parameters, namely, power-law index (1 ≥ n ≥ 0.2), dimensionless rotational velocity (0 ≤ α ≤ 6) and the Reynolds number (0.1 ≤ Re ≤ 40). Over the range of Reynolds number, the flow is known to be steady. Detailed streamline and vorticity contours adjacent to the rotating cylinder and surface pressure profiles provide further insights into the nature of flow. Finally, the paper is concluded by comparing the present numerical results with the scant experimental data on velocity profiles in the vicinity of a rotating cylinder available in the literature. The correspondence is seen to be excellent for Newtonian and inelastic fluids.     

PIV, 2D-LIF and 1D-Raman measurements of flow field,composition and temperature in premixed gas turbine flames

Several laser diagnostic measurement techniques have been applied to study the lean premixed natural gas/air flames of an industrial swirl burner. This was made possible by equipping the burner with an optical combustion chamber that was installed in the high-pressure test rig facility at the DLR Institute of Combustion Technology in Stuttgart. The burner was operated with preheated air at various operating conditions with pressures up to p = 6 bar and a maximum thermal power of P = 1 MW.The instantaneous planar flow field inside the combustor was studied with particle image velocimetry (PIV). Planar laser induced fluorescence (PLIF) of OH radicals on a single-shot basis was used to determine the shape and the location of the flame front as well as the spatial distribution of reaction products. 1D laser Raman spectroscopy was successfully applied for the measurement of the temperature and the concentration of major species under realistic gas turbine conditions.Results of the flow field analysis show the shape and the size of the main flow regimes: the inflow region, the inner and the outer recirculation zone. The highly turbulent flow field of the inner shear layer is found to be dominated by small and medium sized vortices. High RMS fluctuations of the flow velocity in the exhaust gas indicate the existence of a rotating exhaust gas swirl. From the PLIF images it is seen that the primary reactions happened in the shear layers between inflow and the recirculation zones and that the appearance of the reaction zones changed with flame parameters. The results of the multiscalar Raman measurements show a strong variation of the local mixture fraction allowing conclusions to be drawn about the premix quality. Furthermore, mixing effects of unburnt fuel and air with fully reacted combustion products are studied giving insights into the processes of the turbulence–chemistry interaction.