Experimental study on swirling flow of dilute surfactant solution with deformed free-surface

An experimental investigation was performed on a swirling flow of dilute surfactant solution with deformed free-surface in a cylindrical container driven by the constantly rotating bottom wall. The purpose of the experiment was to estimate weak viscoelasticity in the tested surfactant solutions as well as to investigate the flow characteristics. The tested fluid was an aqueous solution of CTAC (CTAC: cetyltrimethyl ammonium chloride), which is a cationic surfactant. Water, 40 ppm, 60 ppm and 200 ppm CTAC solution flows were tested at Froude numbers ranging from 2.59 to 16.3. Particle image velocimetry (PIV) was used to measure the secondary velocity field in the meridional plane. The deformed free-surface level was extracted from the PIV images. At a similar Froude number, the depth of the dip formed at the center region of the free surface was decreased for CTAC solution flow compared with water flow. The inertia-driven vortex at the up-right corner in the meridional plane becomes more and more weakened with increase of the solution concentration or viscoelasticity. Through analyzing the overall force balance compared with water flow, the first normal stress difference characterizing the viscoelasticity was estimated for the dilute CTAC solution flows. The result supports the viscoelasticity-based turbulent drag-reduction mechanism of surfactant solution flow.     

Flow of drag-reducing surfactant solutions in rough pipes

The paper reports the results of experimental study of the flow of hexadecyltrimethylammonium chloride (CTAC) solutions with addition of sodium salicylate (NaSal) in the rough pipes. Measurements were performed in the range of the surfactant concentration from 200 to 400 ppm at a constant molar ratio CTAC/NaSal of 1:2. Five pipes of the relative roughness k/D varying from 1.2 × 10^?2 to 5.6 × 10^?2, obtained by the covering of inner surface of the pipes with glued silicon carbide particles of different size, were studied. The roughness was observed to increase the drag of flow of CTAC/NaSal solutions already at Reynolds numbers higher than 800. With increasing relative roughness k/D, the critical value of Reynolds number, at which the drag reduction disappears, was found to decrease. However, no influence of the roughness on the critical shear stress was noted. The ratio of the critical Reynolds number for rough pipes to that of hydraulically smooth pipes was independent of the surfactant concentration. The degree of drag reduction by the flow of surfactants was greater in rough pipes than in smooth pipes.     

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.     

Instabilities in viscoelastic flow through an axisymmetric sudden contraction

Kinematics and dynamics of the viscoelastic flow in an axisymmetric 4 : 1 sudden contraction geometry are studied for a highly elastic polyisobutylene (PIB) based polymer solution (referred to as PIB-Boger fluid). The critical conditions for the onset of the elastic instabilities and the dynamics of the resulting secondary flows are measured for various flow rates. The spatio-temporal characteristics of the flow are determined by instantaneous pressure measurements and streakline photography. The nonlinear dynamics of the global flow field both upstream and downstream of the contraction plane are systematically examined. New dynamic flow behavior and elastic instabilities downstream of the contraction plane are reported. It is shown that the instantaneous pressure measurements along with flow visualization can be used as an effective tool to characterize viscoelastic flows in complex geometries.     

DNS of viscoelastic turbulent channel flow with rectangular orifice at low Reynolds number

Direct numerical simulations of turbulent viscoelastic-fluid flow in a channel with a rectangular orifice were performed to investigate the influence of viscoelasticity on turbulence statistics and turbulent structures downstream of the orifice. The geometry considered is periodic rectangular orifices with 1:2 expansion. The constitutive equation follows the Giesekus model, valid for polymer (or surfactant) solutions, which are generally capable of reducing the turbulent frictional drag in a smooth channel. The friction Reynolds number and the Weissenberg number were set to 100 and 20-30, respectively. A drag reduction of about 20% was achieved in the viscoelastic flows. The onset Reynolds number for the transition from a symmetric to an asymmetric state was found to be shifted to higher values than that for the Newtonian flow. In the viscoelastic flow, the turbulent kinetic energy was decreased and fewer turbulent eddies were observed, as the Kelvin-Helmholtz vortices were quickly damped. Away from the orifice, quasi-streamwise vortices in the viscoelastic flow were sustained for a longer period, accompanied by energy exchange from elastic energy of the viscoelastic fluid to kinetic energy.     

A low cost laser-speckle photographic technique for velocity measurement in slow flows

A low cost, low power laser-speckle photographic technique has been developed and is duscussed for the measurement of point velocities in slow laminar flows. The technique is particularly suitable for axisymmetric flows where the two velocity components can be easily measured. The accuracy of the technique is established by measurement of the velocity distribution for Poiseuille flow and from data obtained for acceleration of an inelastic Newtonian fluid through a four-to-one circular contraction. Preliminary results are also presented in the contracting flow field for a non shear-thinning highly elastic fluid. These data are particularly significant for verification of finite element numerical solutions currently being developed for viscoelastic fluids in circular entry flows.     

Pressure drop and heat transfer characteristics of boiling water in sub-hundred micron channel

The current work focuses on the pressure drop, heat transfer and stability in two phase flow in microchannels with hydraulic diameter of less than one hundred microns. Experiments were conducted in smooth microchannels of hydraulic diameter of 45, 65 μm, and a rough microchannel of hydraulic diameter of 70 μm, with deionised water as the working fluid. The local saturation pressure and temperature vary substantially over the length of the channel. In order to correctly predict the local saturation temperature and subsequently the heat transfer characteristics, numerical techniques have been used in conjunction with the conventional two phase pressure drop models. The Lockhart–Martinelli (liquid–laminar, vapour–laminar) model is found to predict the two phase pressure drop data within 20%. The instability in two phase flow is quantified; it is found that microchannels of smaller hydraulic diameter have lesser instabilities as compared to their larger counterparts. The experiments also suggest that surface characteristics strongly affect flow stability in the two phase flow regime. The effect of hydraulic diameter and surface characteristics on the flow characteristics and stability in two phase flow is seldom reported, and is of considerable practical relevance.     

Turbulent channel flow of dilute polymeric solutions: Drag reduction scaling and an eddy viscosity model

Direct numerical simulation of viscoelastic turbulent channel flows up to the maximum drag reduction (MDR) limit has been performed. The simulation results in turn have been used to develop relationships between the flow and fluid rheological parameters, i.e. maximum chain extensibility, Reynolds number, Re_τ, and Weissenberg number, We_τ and percent drag reduction (%DR) as well as the slope increment of the mean velocity profile. Moreover, based on the trends observed in the mean velocity profile and the overall momentum balance three different regimes of drag reduction (DR), namely, low drag reduction (LDR; 0 ≤ %DR ≤ 20), high drag reduction (HDR; 20 ≤ %DR ≤ 52) and MDR (52 ≤ %DR ≤ 74) have been identified and mathematical expressions for the eddy viscosity in these regimes are presented. It is found that both in LDR and HDR regimes the eddy viscosity varies with the distance from the channel wall. However, in the MDR regime the ratio of the eddy viscosity to the Newtonian one tends to a very small value around 0.1 within the channel. Based on these expressions a procedure that relies on the DNS predictions of the budgets of momentum and viscoelastic shear stress is developed for evaluating the mean velocity profile.     

Efficient microfluidic rectifiers for viscoelastic fluid flow

In this work we propose a new type of microfluidic rectifier, which is able to operate efficiently under creeping flow conditions. The flow of Newtonian and non-Newtonian fluids was investigated experimentally in different microchannels with triangular (nozzle/diffuser) and hyperbolic shapes in order to achieve high anisotropic flow resistance between the two flow directions. The Newtonian fluid used was de-ionized water and the viscoelastic fluids were aqueous solutions of polyacrylamide and polyethylene oxide with different molecular weights. Pressure drop measurements were performed in addition to visualizations of the flow patterns by streak line photography for a wide range of flow rates. For the Newtonian flows, inertia leads to the appearance of recirculations for both flow directions, but no significant rectification effects appear. For the viscoelastic fluids, two distinct behaviors are identified: at low flow rates, the pressure drops are similar in both flow directions; above a critical flow rate (or Deborah number), the flow patterns become quite different, leading to different flow rates in the forward and backward flow directions for the same pressure drop, i.e., rectification effects emerge. In particular, the viscoelastic fluid flow becomes unsteady in the forward direction, due to the presence of elastic instabilities, which leads to a significant increase in the flow resistance. Flow resistance ratios greater than three were achieved for the hyperbolic rectifier, clearly in excess of the value for the triangular-shaped rectifier and for other geometries proposed in the literature for operation in creeping flow conditions. This high diodicity is associated with the distinct nature of the extensional flows in the forward and backward directions of the hyperbolic-type microgeometry.     

Drag reduction for liquid flow through micro-apertures

Pressure drops in the flow through micro-orifices and capillaries were measured for silicone oils, aqueous solutions of polyethylene glycol (PEG), and surfactant aqueous solutions. The diameter of micro-orifices ranged from 5 μm to 400 μm. The corresponding length/diameter ratio was from 4 to 0.05 and capillary diameters were 105 μm and 450 μm. The following results were obtained: silicone oils of 10^?6 m²/s and 10^?5 m²/s in kinematic viscosity generated a reduction of pressure drop (RPD), that is, drag reduction, similar to the RPD of water and a glycerol/water mixture reported in the previous paper by the present authors. When RPD occurred, the pressure drop (PD) of silicone oils of 10^?6 m²/s and 10^?5 m²/s had nearly the same magnitude. Namely, the difference in viscosity did not influence RPD. A 10³ ppm aqueous solution of PEG20000 provided almost the same PD as that of PEG8000 for the 400 μm to 15 μm orifices, but a greater PD than that of PEG8000 for the 10 μm to 5 μm orifices. A non-ionic surfactant and a cationic surfactant were highly effective in RPD compared with anionic surfactants: the non-ionic and cationic surfactant solutions had PD one order of magnitude lower than that of water under some flow conditions in the concentration range from 1 ppm to 10⁴ ppm, but the anionic surfactant solutions did not generate RPD except in the case of the smallest orifice of 5 μm in diameter. The PD of the non-ionic surfactant solution showed a steep rise at a Reynolds number (Re_t) for 400 μm to 15 μm orifices. The Re_t provides the relationship Re_t = K/D, where D is the orifice diameter, and K is a constant of 2 × 10^?2 m for the 100–20 μm orifices irrespective of liquid concentration. Capillary flow experiment revealed that the PEG, non-ionic and cationic surfactant solutions generated RPD also in a laminar flow through the capillary of 105 μm in diameter, but not in the flow through the capillary of 450 μm in diameter. In order to clarify the cause of RPD, an additional experiment was carried out by changing the orifice material from metal to acrylic resin. The result gave a different appearance of RPD, suggesting that RPD is related to an interfacial phenomenon between the liquid and wall. The large RPDs found in the present experiment are very interesting from both academic and practical viewpoints.     

Laminar,transitional and turbulent annular flow of drag-reducing polymer solutions

Mean and rms axial velocity-profile data obtained using laser Doppler anemometry are presented together with pressure-drop data for the flow through a concentric annulus (radius ratio κ = 0.506) of a Newtonian (a glycerine–water mixture) and non-Newtonian fluids—a semi-rigid shear-thinning polymer (a xanthan gum) and a polymer known to exhibit a yield stress (carbopol). A wider range of Reynolds numbers for the transitional flow regime is observed for the more shear-thinning fluids. In marked contrast to the Newtonian fluid, the higher shear stress on the inner wall compared to the outer wall does not lead to earlier transition for the non-Newtonian fluids where more turbulent activity is observed in the outer wall region. The mean axial velocity profiles show a slight shift (～5%) of the location of the maximum velocity towards the outer pipe wall within the transitional regime only for the Newtonian fluid.     

Super-stable parallel flows of multiple visco-plastic fluids

We consider the stability of a multi-layer plane Poiseuille flow of two Bingham fluids. It is shown that this two-fluid flow is frequently more stable than the equivalent flow of either fluid alone. This phenomenon of super-stability results only when the yield stress of the fluid next to the channel wall is larger than that of the fluid in the centre of the channel, which need not have a yield stress. Our result is in direct contrast to the stability of analogous flows of purely viscous generalised Newtonian fluids, for which short wavelength interfacial instabilities can be found at relatively low Reynolds numbers. The results imply the existence of parameter regimes where visco-plastic lubrication is possible, permitting transport of an inelastic generalised Newtonian fluid in the centre of a channel, lubricated at the walls by a visco-plastic fluid, travelling in a stable laminar flow at higher flow rates than would be possible for the single fluid alone.     

Impact of aspect ratio on flow boiling of water in rectangular microchannels

In this paper we focus on the impact of varying the aspect ratio of rectangular microchannels, on the overall pressure drop involving water boiling. An integrated system comprising micro-heaters, sensors and microchannels has been realized on (1 1 0) silicon wafers, following CMOS compatible process steps. Rectangular microchannels were fabricated with varying aspect ratios (width [W] to depth [H]) but constant hydraulic diameter of 142 ± 2 μm and length of 20 mm. The invariant nature of the hydraulic diameter is confirmed through two independent means: physical measurements using profilometer and by measuring the pressure drop in single-phase fluid flow. The experimental results show that the pressure drop for two-phase flow in rectangular microchannels experiences minima at an aspect ratio of about 1.6. The minimum is possibly due to opposing trends of frictional and acceleration pressure drops, with respect to aspect ratio. In a certain heat flux and mass flux range, it is observed that the two-phase pressure drop is lower than the corresponding single-phase value. This is the first study to investigate the effect of aspect ratio in two-phase flow in microchannels, to the best of our knowledge. The results are in qualitative agreement with annular flow model predictions. These results improve the possibility of designing effective heat-sinks based on two-phase fluid flow in microchannels.     

Analysis of an organized turbulent structure using a pattern recognition technique in a drag-reducing surfactant solution flow

This work presents the investigation for an organized turbulent structure in a drag-reducing flow of dilute surfactant solution by utilizing a particle image velocimetry system to perform the pattern recognition technique on a trajectory in four quadrants of streamwise and wall-normal velocity fluctuations. The pattern recognition is added to a new algorithm in order to directly capture the spatial rotation motion. The Reynolds number based on the channel height and bulk mean velocity was set to 1.5 × 10⁴. Surfactant solution with a weight concentration of 150 ppm was employed and the drag reduction rate was 65%. In the drag-reducing flow, we observe increased frequencies of occurrence of the flow events that correspond to a meandering motion in the wall-normal direction of the high-and low-speed regions. Three findings from investigation of the ensemble-averaged Reynolds shear stress and vortex structure are as follows: (i) the Reynolds shear stress in the large fluctuation range occurs in the narrow region; (ii) Size, strength, arrangement and inclination in the spatial vortex structure in the drag-reducing flow differ from those of the water; and (iii) all trajectory contributions for the wall friction coefficient decrease. Finally, we interpreted that the viscoelasticity characterizing the viscoelastic stress and relaxation time in rheological properties of the flow changes specific elementary vortex for the drag-reducing flow, and the trajectories of each flow pattern change drastically.     

A 1-D theory for extensional deformation of a viscoelastic filament under exponential stretching

We consider a viscoelastic filament placed between two coaxial discs, with the bottom plate fixed and the top plate pulled at an exponential rate. Using a slender rod approximation, we derive a one-dimensional (1-D) model which describes the deformation of a viscoelastic filament governed by the Oldroyd-B constitutive model. It is assumed that the flow is axisymmetric and that inertia and gravity are negligible. One solution of the model equations corresponds to ideal uniaxial elongation. A linear stability analysis shows that this solution is unstable for a Newtonian fluid and for viscoelastic filaments with small Deborah number (De ≤ 0.5). For Deborah number greater than 0.5, ideal uniaxial elongation is linearly stable. Numerical solution of the nonlinear equations confirms the result of the linear stability analysis. For initial conditions close to ideal uniaxial flow, our results show that if De > 0.5, the central portion of the filament undergoes considerable strain hardening. As a result, the sample remains almost cylindrical and the deformation approaches pure uniaxial extension as the Hencky strain increases. For De ≤ 0.5, the Trouton ratio based on the effective extension rate at the mid-plane radius gives a much better approximation to the true extensional viscosity than that based on the imposed stretch rate.     

Purely elastic instabilities in three-dimensional cross-slot geometries

Creeping and low Reynolds number flows of an upper-convected Maxwell (UCM) fluid are investigated numerically in a three-dimensional orthogonal cross-slot geometry. We analyze two different flow configurations corresponding to uniaxial extension and biaxial extension, and assess the effects of extensional flow type, Deborah and Reynolds numbers on flow dynamics near the interior stagnation point. Using these two flow arrangements the amount of stretch and compression near the stagnation point can be varied, providing further insights on the viscoelastic flow instability mechanisms in extensionally dominated flows with an interior stagnation point. The uniaxial extensional flow arrangement leads to the onset of a steady flow asymmetry, followed by a second purely elastic flow instability that generates an unsteady flow at higher flow rates. On the other hand, for the biaxial extension flow configuration a symmetric flow is observed up to the critical Deborah number when the time-dependent purely elastic instability sets in, without going through the steady symmetric to steady asymmetric transition.     

Spherical Couette flow of a viscoelastic fluid Part I: Experimental study of the inner sphere rotation

Flow of viscoelastic fluid contained between a stationary outer sphere and a rotating inner sphere is studied experimentally. In the present investigation, relatively low-concentration polyacrylamide-water solutions are used as viscoelastic fluid, and for the sake of comparison glycerin-water solutions are used as the Newtonian fluid. In experiments, measurements of the rotational torque and the velocity profile in the meridional plane of the spherical gap are made. Various transition phenomena of flow modes that are unique to viscoelastic fluids are investigated by flow visualization for a wide range of rotational Reynolds numbers. Experimental results revealed that a roll-cell-like structure (banded radial structure) caused by elastic instability is generated in the polar region, and that it propagates toward the equatorial region when rotation of the inner sphere is increased, resulting in the formation of two distinct regions: the elastic dominant region and the inertial dominant region. Flow modes are classified and the critical Reynolds numbers are obtained for different gap widths. A correlation is obtained for the torque data in the regime before the onset of instability.     

An experimental investigation of the flow about a sphere rotating in a rivlin-ericksen fluid

Tangential and radial velocity profiles were measured for the flow about a sphere rotating slowly in a Newtonian fluid, contained in a rectangular tank. Velocities were determined from enlarged streak photographs of aluminium particles moving in a collimated “sheet” of light, at several planes through the flow field. Similar velocity profiles were measured for the flow of a 1.50% Natrosol 250 H solution about two spheres of different diameters rotating in two different sized rectangular tanks. A set of velocity distributions were also measured for a sphere rotating in a 0.9% Natrosol 250 H solution. A dye tracer study of the flow about a sphere rotating in this liquid is presented as well. Both Natrosol solutions exhibited viscoelastic behaviour. The Newtonian fluid study was carried out at a Reynolds number of 1.2 and the viscoelastic fluid studies were within the Reynolds number range of 0.05–1.24.The zero shear viscosities of the Natrosol solutions were measured using the falling-sphere method. The non-Newtonian material parameters were obtained by fitting the theoretical curves to the measured velocity data. The values of the elastic and shear thinning parameters for the two liquids obtained in the different geometrical and dynamical situations are compared.     

Simulation of a valveless pump with an elastic tube

A valveless pump consisting of a pumping chamber with an elastic tube was simulated using an immersed boundary (IB) method. The interaction between the motion of the elastic tube and the pumping chamber generated a net flow toward the outlet throughout a full cycle of the pump. The net flow rate of the valveless pump was examined by varying the stretching coefficient (ϕ), bending coefficient (γ), the aspect ratio (l/d) of the elastic tube, and the frequency (f) of the pumping chamber. As the stretching and bending coefficients of the elastic tube increased, the net flow through the valveless pump decreased. Elastic tubes with aspect ratios in the range of 2 ⩽ l/d ⩽ 3 generated a higher flow rate than that generated for tubes with aspect rations of l/d = 1 or 4. As the frequency of the pumping chamber increased, the net flow rate of the pump for l/d = 2 increased. However, the net flow rate for l/d = 3 was nonlinearly related to the pumping frequency due to the complexity of the wave motions. Snapshots of the fluid velocity vectors and the wave motions of the elastic tube were examined over one cycle of the pump to gain a better understanding of the mechanism underlying the valveless pump. The relationship between the average gap in the elastic tube and the average flow rate of the pump was analyzed. A smaller gap in the elastic tube during the expansion mode and a wider gap in the elastic tube during the contraction mode played a dominant role in generating a high average flow rate in the pump, regardless of the stretching coefficient (ϕ), the aspect ratio (l/d) of the elastic tube, or the pumping frequency of the pumping chamber (f).