Investigation of slip boundary conditions in the T-shaped microchannel

Flow regimes and mixing performance in a T-type micromixer at high Reynolds numbers were studied by numerical solution of the Navier–Stokes equations. The Reynolds number was varied from 1 to 1000. The cross section of the mixing channel was 100 μm × 200 μm, and its length was 1400 μm. The transverse inlet channels were symmetric about the mixing channel, and their cross-section was 100 μm × 100 μm, and the total length was 800 μm. Five different flow regimes were identified: (i) stationary vortex-free flow (Re < 5); (ii) stationary symmetric vortex flow with two horseshoe vortices (5 < Re < 150); (iii) stationary asymmetric vortex flow (150 < Re < 240); (iv) non-stationary periodic flow (240 < Re < 400); and (v) stochastic flow (Re > 400). Maximum mixing efficiency is obtained for nonstationary asymmetric vortex flow. In this case, an S-shaped vortex structure is formed in the flow field. The effect of the slip conditions on the flow pattern and mixing efficiency is studied. The slip length varied from 1 to 70 μm in the calculations. It was shown that the mixing can be controlled by hydrophobic coating.     

Influence of aspect ratio on the mean flow field of a surface-mounted finite-height square prism

The flow around surface-mounted, finite-height square prisms at a Reynolds number of Re = 4.2 × 10⁴ was investigated experimentally in a low-speed wind tunnel using particle image velocimetry. The thickness of the boundary layer on the ground plane relative to the width of the prism was δ/D = 1.5. Four prism aspect ratios were tested, AR = 9, 7, 5, and 3, to study how the aspect ratio influences the flow field close to the prism. Upstream of the prism, lowering the aspect ratio from AR = 9 to AR = 3 causes the stagnation point on the upstream face to move closer to the free end, but there is no influence on the location and strength of the horseshoe vortex. Lowering the aspect ratio from AR = 9 to AR = 3 causes the cross-stream vortices in the upper and lower halves of the wake to move downstream and upstream, respectively; the latter vortex is absent for AR = 3, suggesting this prism sits below the critical aspect ratio. Above the free end of the prism, within the region of separated flow, lowering the aspect ratio from AR = 9 to AR = 3 shifts the location of the cross-stream vortex farther downstream. For the prism of AR = 3, reverse flow above the free end is stronger yet more unsteady compared to the more slender prisms, while the streamwise edge vortices are smaller and weaker.     

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.     

Separated flow structures around a cylindrical obstacle in a narrow channel

A detailed experimental study is performed on the separated flow structures around a low aspect-ratio circular cylinder (pin-fin) in a practical configuration of liquid cooling channel. Distinctive features of the present arrangement are the confinement of the cylinder at both ends, water flow at low Reynolds numbers (Re = 800, 1800, 2800), very high core flow turbulence and undeveloped boundary layers at the position of the obstacle. The horseshoe vortex system at the junctions between the cylinder and the confining walls and the near wake region behind the obstacle are deeply investigated by means of Particle Image Velocimetry (PIV). Upstream of the cylinder, the horseshoe vortex system turns out to be perturbed by vorticity bursts from the incoming boundary layers, leading to aperiodical vortex oscillations at Re = 800 or to break-away and secondary vorticity eruptions at the higher Reynolds numbers. The flow structures in the near wake show a complex three-dimensional behaviour associated with a peculiar mechanism of spanwise mass transport. High levels of free-stream turbulence trigger an early instabilization of the shear layers and strong Bloor–Gerrard vortices are observed even at Re = 800. Coalescence of these vortices and intense spanwise flow inhibit the alternate primary vortex shedding for time periods whose length and frequency increase as the Reynolds number is reduced. The inhibition of alternate vortex shedding for long time periods is finally related to the very large wake characteristic lengths and to the low velocity fluctuations observed especially at the lowest Reynolds number.     

Experimental study on turbulent mixing of spray droplets in crossflow

Centrifugal spray injected at various angles in gas crossflow has been studied experimentally using PIV visualization system and image-processing techniques. Experiments were carried out inside a rectangular duct (95 mm × 95 mm in cross-section) at ambient temperature and pressure, with different gas Reynolds numbers (vary from 12,900 to 45,000) and three injection angles (60°, 90° and 120°). The spray angle of the centrifugal nozzle is 80°, with D₃₂ of 80 μm. The instantaneous images of droplets distribution and the values of the line-averaged D₃₂ at different positions on the cross-sections along the flow field for each condition were obtained, and their flow field configurations were achieved. Quantitative assessments of mixing degree between two phases for different injection angles were determined using a spatial unmixedness parameter. It is found that the addition of droplets into the gas crossflow enhanced the turbulence intensity of the gas crossflow and caused different-scale vortices. The flow field structure, to a great extent, is dependent on the injection angle. The entrainment and centrifugal force of large vortex lead to uneven droplet distribution and moreover influence the mixing of droplets and gas crossflow. A better mixing result can be obtained with the injection angle of nozzles of 60°.     

Tip gap height effects on flow structure and heat/mass transfer over plane tip of a high-turning turbine rotor blade

The effects of tip gap height-to-chord ratio, h/c, on the flow structure and heat/mass transfer over the plane tip surface of a large-scale high-turning turbine rotor blade have been investigated for h/c = 1.0%, 2.0%, 3.0% and 4.0%. For near-wall tip gap flow visualizations, a high-resolution oil film method is employed, and the naphthalene sublimation technique is used for local heat/mass transfer rate measurements. From the tip surface visualizations, a pair of vortices named “tip gap vortices” is identified in the leading edge region within the tip gap. The overall tip gap flow is characterized not only by the tip gap vortices but also by the flow separation/re-attachment process along the pressure-side tip edge. Within the separation bubble, there exist complicated near-wall flows moving toward a mid-chord flow converging area. With increasing h/c, the tip gap vortices, the flow separation/re-attachment, and the converging flows within the separation bubble tend to be intensified. In general, higher thermal load is found along the loci of the tip gap vortices and along the re-attachment line, while lower thermal load is observed behind the tip gap vortex system and near the mid-chord flow converging area. Heat/mass transfer characteristics with the variation of h/c are discussed in detail in conjunction with the tip gap flow features. Based on the flow visualizations and heat/mass transfer data, new realistic tip gap flow models have been proposed for h/c = 1.0 and 4.0%.     

Pulsating flow in a planar diffuser upstream of automotive catalyst monoliths

The flow distribution across automotive exhaust catalysts has a significant effect on their conversion efficiency. The exhaust gas is pulsating and flow distribution is a function of engine operating condition, namely speed (frequency) and load (flow rate). This study reports on flow measurements made across catalyst monoliths placed downstream of a wide-angled planar diffuser presented with pulsating flow. Cycle-resolved particle image velocimetry (PIV) measurements were made in the diffuser and hot wire anemometry (HWA) downstream of the monoliths. The ratio of pulse period to residence time within the diffuser (defined as the J factor) characterises the flow distribution. During acceleration the flow remained attached to the diffuser walls for some distance before separating near the diffuser inlet later in the cycle. Two cases with J ∼ 3.5 resulted in very similar flow fields with the flow able to reattach downstream of the separation bubbles. With J = 6.8 separation occurred earlier with the flow field resembling, at the time of deceleration, the steady flow field. Increasing J from 3.5 to 6.8 resulted in greater flow maldistribution within the monoliths; steady flow producing the highest maldistribution in all cases for the same Re.     

Effect on the flow and heat transfer characteristics for sinusoidal pulsating laminar flow in a heated square cylinder

Two-dimensional numerical simulation is performed to understand the effect of flow pulsation on the flow and heat transfer from a heated square cylinder at Re = 100. Numerical calculations are carried out by using a finite volume method based on the pressure-implicit with splitting of operators algorithm in a collocated grid. The effects of flow pulsation amplitude (0.2 ≤ A ≤ 0.8) and frequency (0 ≤ f _p  ≤ 20 Hz) on the detailed kinematics of flow (streamlines, vorticity patterns), the macroscopic parameters (drag coefficient, vortex shedding frequency) and heat transfer enhancement are presented in detail. The Strouhal number of vortices shedding, drag coefficient for non-pulsating flow are compared with the previously published data, and good agreement is found. The lock-on phenomenon is observed for a square cylinder in the present flow pulsation. When the pulsating frequency is within the lock-on regime, time averaged drag coefficient and heat transfer from the square cylinder is substantially augmented, and when the pulsating frequency in about the natural vortex shedding frequency, the heat transfer is also substantially enhanced. In addition, the influence of the pulsating amplitude on the time averaged drag coefficient, heat transfer enhancement and lock-on occurrence is discussed in detail.     

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.     

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.     

Direct numerical simulation of a NACA0012 in full stall

This work aims at investigating the mechanisms of separation and the transition to turbulence in the separated shear-layer of aerodynamic profiles, while at the same time to gain insight into coherent structures formed in the separated zone at low-to-moderate Reynolds numbers. To do this, direct numerical simulations of the flow past a NACA0012 airfoil at Reynolds numbers Re = 50,000 (based on the free-stream velocity and the airfoil chord) and angles of attack AOA = 9.25° and AOA = 12° have been carried out. At low-to-moderate Reynolds numbers, NACA0012 exhibits a combination of leading-edge/trailing-edge stall which causes the massive separation of the flow on the suction side of the airfoil. The initially laminar shear layer undergoes transition to turbulence and vortices formed are shed forming a von Kármán like vortex street in the airfoil wake. The main characteristics of this flow together with its main features, including power spectra of a set of selected monitoring probes at different positions on the suction side and in the wake of the airfoil are provided and discussed in detail.     

The effect of incidence angle on the mean wake of surface-mounted finite-height square prisms

The mean wake of a surface-mounted finite-height square prism was studied experimentally in a low-speed wind tunnel to explore the combined effects of incidence angle (α) and aspect ratio (AR). Measurements of the mean wake velocity field were made with a seven-hole pressure probe for finite square prisms of AR = 9, 7, 5 and 3, at a Reynolds number of Re = 3.7 × 10⁴, for incidence angles from α = 0° to 45°. The relative thickness of the boundary layer on the ground plane, compared to the prism width, was δ/D = 1.5. As the incidence angle increases from α = 0° to 15°, the mean recirculation zone shortens and the mean wake shifts in the direction opposite to that of the mean lift force. The downwash is also deflected to this side of the wake and the mean streamwise vortex structures in the upper part of the wake become strongly asymmetric. The shortest mean recirculation zone, and the greatest asymmetry in the mean wake, is found at the critical incidence angle of α_critical ≈ 15°. As the incidence angle increases from α = 15° to 45°, the mean recirculation zone lengthens and the mean streamwise vortex structures regain their symmetry. These vortices also elongate in the wall-normal direction and become contiguous with the horseshoe vortex trailing arms. The mean wake of the prism of AR = 3 has some differences, such as an absence of induced streamwise vorticity near the ground plane, which support its classification as lying below the critical aspect ratio for the present flow conditions.     

Flow pattern and molecular visualization of DNA solutions through a 4:1 planar micro-contraction

The present work explores unusual flow behavior of entangled fluids in an abrupt contraction flow device. Fluorescent imaging was carried out on four different entangled DNA solutions with concentrations ranging from 0.1 to 1.0% (with a wide range of entanglements per chain Z = 7–55). For weakly entangled solutions (Z < 30), vortex flow was dominant at high flow rates. However, for well-entangled DNA solutions (Z ≥ 30), unusual time dependant shear banding was observed at the contraction entrance. Upon reducing the slip length by adding sucrose to the well-entangled DNA solution, vortex flow became dominant again. In vortex flow, most DNA chains remained coiled at the corner in regular recirculation. However, when jerky-shear-banding flow developed, significant stable stretching of DNA chains occurred at the center-line, with quasi-periodic switching between stretching and recoil at the corner.     

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.     

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.     

URANS vs. experiments of flow induced motions of multiple circular cylinders with passive turbulence control

Two-dimensional Unsteady Reynolds-Average Navier–Stokes equations with the Spalart–Allmaras turbulence model are used to simulate the flow induced motions of multiple circular cylinders with passive turbulence control (PTC) in steady uniform flow. Four configurations with 1, 2, 3, and 4 cylinders in tandem are simulated and studied at a series of Reynolds numbers in the range of 30 000<Re<120 000. Simulation results are verified by experimental data measured in the Marine Renewable Energy Laboratory. Good agreement was observed between the values of vorticity, amplitude ratio, and frequency ratio predicted by numerical simulations and experimental measurements. The amplitude and frequency response show the initial and upper branches in vortex induced vibration (VIV), transition from VIV to galloping, and galloping branch for all PTC-cylinders. The maximum amplitude of 2.9 diameters for the first cylinder is achieved at Re=104 356 in the numerical results. Compared with the first cylinder, the VIV initial branch starts at higher Re for the downstream cylinders due to the presence of the upstream cylinder(s). 2P and 2P+2S vortex patterns are observed at Re=62 049 and Re=90 254 for the single PTC-cylinder. Furthermore, the shed vortices of the downstream cylinders are strongly disrupted and modified by the vortices shed from the upstream one in the cases of multiple PTC-cylinders.     

Direct numerical simulation of a particle-laden mixing layer with a chemical reaction

A direct numerical simulation (DNS) is applied to a particle-laden turbulent mixing layer with a chemical reaction, and the effects of particles on turbulence and chemical species’ diffusion and reaction in both zero and finite gravity cases are investigated. Unreactive particles, whose response time, τ_P, is smaller than the Kolmogorov time scale, τ_K [τ_P/τ_K = O(10⁻¹)], are uniformly injected into the high-speed side of the mixing layer. Two reactive chemical species are separately introduced through different sides. The results show that although laden particles generally depress turbulent intensities, they begin to enhance turbulent intensities downstream as the particle size decreases provided that the inlet particle volume fraction is fixed. This is because that the small particles with small particle response time easily accumulate at the circumference of coherent vortices and act to suppress the growth of primitive small-scale coherent vortices upstream but enhance that of relatively developing large-scale ones downstream. Also, since the small-scale turbulence, which promotes the chemical reaction, is suppressed by the laden particles in the entire region, chemical product decreases overall. Furthermore, the presence of finite gravity on the particles acts to depress the turbulent intensities, but its effects on chemical species’ diffusion and reaction are quite small.     

PIV study of large-scale flow organisation in slot jets

The paper reports on particle image velocimetry (PIV) measurements in turbulent slot jets bounded by two solid walls with the separation distance smaller than the jet width (5–40%). In the far-field such jets are known to manifest features of quasi-two dimensional, two component turbulence. Stereoscopic and tomographic PIV systems were used to analyse local flows. Proper orthogonal decomposition (POD) was applied to extract coherent modes of the velocity fluctuations. The measurements were performed both in the initial region close to the nozzle exit and in the far fields of the developed turbulent slot jets for Re ⩾ 10,000. A POD analysis in the initial region indicates a correlation between quasi-2D vortices rolled-up in the shear layer and local flows in cross-stream planes. While the near-field turbulence shows full 3D features, the wall-normal velocity fluctuations day out gradually due to strong wall-damping resulting in an almost two-component turbulence. On the other hand, the longitudinal vortex rolls take over to act as the main agents in wall-normal and spanwise mixing and momentum transfer. The quantitative analysis indicates that the jet meandering amplitude was aperiodically modulated when arrangement of the large-scale quasi-2D vortices changed between asymmetric and symmetric pattern relatively to the jet axis. The paper shows that the dynamics of turbulent slot jets are more complex than those of 2D, plane and rectangular 3D jets. In particular, the detected secondary longitudinal vortex filaments and meandering modulation is expected to be important for turbulent transport and mixing in slot jets. This issue requires further investigations.     

CFD study: Effect of pulsating flow on gas–solid hydrodynamics in FCC riser

Gas–solid flow in a fluid catalytic cracking (FCC) riser exhibits poor mixing in the form of a core–annulus flow pattern and a dense bottom/dilute top distribution of solids. To enhance gas–solid mixing, studies on dense fluidized beds have suggested using a pulsating flow of gas. The present study investigates the effect of pulsating flow on gas–solid hydrodynamics inside the FCC riser employing computational fluid dynamics. Two flow conditions are investigated: a cold flow of air-FCC catalyst in a pilot-scale riser and a reactive flow in an industrial-scale FCC riser. In the cold-flow riser, pulsating flows cause the slug flow of solids and thus increase the average solid accumulation in the flow domain and solid segregation towards the wall. In the industrial FCC riser, pulsating flows produce radial profiles that are more homogeneous. Pulsating flows further improve the conversion and yield in the initial few metres of height. At 7 m, the conversion from pulsating flow is 59%, compared with 44% in without pulsating flow. The results and analysis presented here will help optimize flow conditions in the circulating fluidized bed riser, in not only FCC but also applications such as fast pyrolysis and combustion.