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.     

Experiments on elongational flow of dilute polymer solutions Part I: Jet reaction and excess pressure drop for the flow through small apertures

Experiments have been made with dilute polymer solutions on the reaction of jets issuing from small orifices and the excess pressure drop for orifice and capillary flows.Under the flow conditions with vortices occurring upstream of the aperture, the jet reaction is nearly zero below some mean velocity for PEO solutions and similarly zero below some generalized Reynolds number for Separan solutions. The normalized jet reactions, when they possess positive values, are correlated with the generalized Reynolds number irrespective of the aperture diameters for both kinds of solution.In most cases, the pressure is higher than in the corresponding water flow, but for some flows with no vortex it is lower. For the vortex flow of PEO solutions the normalized excess pressure drop is inversely proportional to the Reynolds number for both orifices and capillaries, while for Separan solutions this quantity is not correlated with the generalized Reynolds number for orifice flow but is correlated with it for capillary flow.     

On the use of flow through a contraction in estimating the extensional viscosity of mobile polymer solutions

We consider the use of pressure measurements in contraction flows in the determination of the extensional viscosity behaviour of polymer solutions. The experimental data are interpreted on the basis of the recent theory of Binding. The resulting extensional viscosities are compared with those obtained from a commercial Spin Line Rheometer.We conclude that contraction flows provide a convenient means of determining the extensional viscosity of shear-thinning polymer solutions. The case is not so clear for constant viscosity Boger fluids.In the course of the experiments, it is shown that excess pressure losses in the contractions can be brought about by two distinct flow mechanisms in the case of Boger fluids. In the axisymmetric case, both vortex enhancement and excess pressure loss are observed, although there is not a strict one-to-one correlation between these phenomena. In the planar case, vortex enhancement is not conspicuously present, although there is still a substantial excess pressure loss at high flow rates. This excess must be associated with the ‘bulb’ flow field which essentially replaces the vortex-enhancement regime of the axisymmetric case.     

Extensional flow resistance of dilute polyacrylamide and surfactant solutions

The extensional flow behaviour of dilute aqueous solutions of a partiallyhy-drolyzed polyacrylamide and a surfactant were investigated in an extensional flow cell. The cell was designed such that fluids were subjected to steady shear before undergoing extensional motion in a converging channel. Extensional resistance was monitored by measuring the pressure drop through the channel. Such measurements were made over a range of extensional rates at fixed values of the upstream shear rate. Solutions of different concentrations were tested — up to 40 ppm of polyacrylamide and 450 ppm of surfactant — at various temperatures in the case of surfactant and for different types and amounts of salt in the case of polyacrylamide. Of the results, the more notable are that the extensional resistance of polyacrylamide solutions is affected much more by CaCl₂ than by NaCl and that surfactant solutions do not exhibit extensional resistance unless they are pre-sheared.     

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.     

Two-dimensional viscoelastic flows: experimentation and modeling

Extensive experimental data on the birefringence in converging and diverging flows of a polymeric melt have been obtained. The birefringence and pressure drop measurements were carried out in working cells of planar geometry having different contraction angles and contraction ratios. For investigation of diverging or abrupt expansion flow, the direction of flow in the cells was reversed. The theoretical predictions are based upon the Leonov constitutive equation and a finite element scheme with streamwise integration.In contrast to Newtonian and second-order fluids, viscoelastic fluids at high shear rates show significant differences in pressure drop and birefringence (i.e. stresses) in converging and diverging flows. For a constant flow rate, the pressure drop is higher and the birefringence smaller in diverging flows than in converging flows. This difference increases with increasing flow rate. Further, for the same contraction ratio but different contraction angles, the birefringence maximum increases considerably with contraction angle. In addition, an increase in contraction ratio has the same effect.The viscoelastic constitutive equation of Leonov has been shown to describe all the above viscoelastic effects observed in the experiments. In general, a reasonable agreement between theory and experiment has been obtained, which shows the usefulness of the Leonov model in describing actual flows.     

Creation of very-low-Reynolds-number chaotic fluid motions in microchannels using viscoelastic surfactant solution

Solutions of flexible high-molecular-weight polymers or some kinds of surfactant are viscoelastic fluids. The elastic stress is induced in such viscoelastic fluid flows and grows nonlinearly with the flow-rate resulting in many particular flow phenomena, including purely elastic instability. The purely elastic instability can even result in a kind of chaotic fluid motion, the so-called elastic turbulence, which is a recently discovered flow phenomenon and arises at arbitrarily small Reynolds number. By using viscoelastic surfactant solution, we attempted to create the peculiar chaotic fluid motions in several specially designed microchannels in which flows with curvilinear streamlines can be generated. The viscoelastic working fluids were aqueous solutions of surfactant, CTAC/NaSal (cetyltrimethyl ammonium chloride/sodium salicylate). CTAC solutions with weight concentration of 200 ppm (part per million) and 1000 ppm, respectively, at room temperature were tested. For comparison, water flows in the same microchannels were also visualized. The Reynolds numbers for all the microchannel flows were quite small (for solution flows, the Reynolds numbers were the order of or smaller than one) and the flow should be definitely laminar for Newtonian fluid. It was found that the regular laminar flow patterns for low-Reynolds-number Newtonian fluid flow in different microchannels were strongly deformed in solution flows: either asymmetrical flow structures or time-dependent vortical fluid motions appeared. These chaotic flow phenomena were considered to be induced by the viscoelasticity of the CTAC solutions. Discussions about the potential applications using such kind of chaotic fluid motions were also made.     

Drag reduction and heat transfer in surfactant solutions with excess counterion

We present the results of a study of turbulent drag reduction in a small circulating loop using surfactant solutions with excess counterion. In addition, these solutions were used in measurements of heat transfer, both in pipe flow and in an impinging jet. Both frictional drag and heat transfer were reduced in the pipe flow experiments. Measurements of heat transfer in the impinging jet revealed a dependence on the molar concentration ratio of the counterion. When the counterion was added at a molar concentration 30 times higher than that of the surfactant, the resulting surfactant solution did not reduce the rate of heat transfer in the impinging jet. By using this surfactant system in an impinging jet, we show both a reduction in pipe friction and normal heat transfer potential in a circulating heat exchange system. In order to investigate this difference in heat transfer between pipe flows and impinging jet flows, a comparison was made of the wall shear stress between these two flow regimes. The estimated wall shear stress was of the same order in both flows, and thus was not considered to be the primary cause of the difference in heat transfer. It is instead suggested that the micellar structure of the surfactant is influenced by a compressive deformation of the impinging flow in a manner that is different from the shear deformation observed in pipe flow.     

Roles of interfacial dynamics in the interaction behaviours between deformable oil droplets

Drainage and deformation of the intervening film can arguably represent the dynamic nature of colliding soft matters. The development of an interaction force analysis between soft interfaces helps to probe the drop deformation and the interfacial properties. Based on the SRYL model, the fluid flow inside the droplet and the convection-diffusion of the surfactant at the oil-water interface is coupled to model the distribution of a non-ionic surfactant (Span80) during drop deformation using AFM. This study quantifies the in situ interfacial concentration with a trace amount of surfactant at the interface and indicates its effect on the interaction forces between two immersed oil droplets in an aqueous solution.     

The effects of chain conformation in the microfluidic entry flow of polymer–surfactant systems

An associative polymer–surfactant system has been used to observe the effects of chain conformation in the entry flow through a microfabricated planar 16:1:16 contraction–expansion geometry. The well-studied system of the flexible polymer poly(ethylene oxide) (PEO) and anionic surfactant sodium dodecyl sulfate (SDS) was used. Dilute polymer solutions with increasing SDS concentration were characterized in steady and dynamic shear, as well as capillary breakup extensional rheology. Based on this characterization, the primary quantitative difference is an increase in zero-shear viscosity as a result of the PEO chain expansion brought on by association of SDS surfactant micelles. However, these quantitatively similar solutions were observed to exhibit much more qualitatively different flow patterns via fluorescent streak imaging in the entry flow. In contrast to previous work on PEO solutions, the PEO–SDS systems were observed to transition to a steady viscoelastic flow regime characterized by stable lip vortices at much lower elasticity and Weissenberg numbers. The resulting insight gained regarding the utility of microfluidic flows in elucidating effects of subtle conformational changes further illustrates the potential for using microfabricated devices as rheometric tools for measuring the properties of dilute and weakly viscoelastic fluids.     

Liquid distribution and electrical conductivity in foam

Experiments were conducted with aqueous foam generated by bubbling nitrogen through anionic, cationic and nonionic surfactant solutions. The ratio of the electrical conductivity of the foam to that of the liquid was found to increase monotonically with the volumetric liquid fraction in the foam. The choice of surfactant as well as the degree of inhomogeneity in bubble size were found to be without effect on the relationship. However, for a fixed liquid fraction, it was found that decreasing the mean bubble size can decrease the conductivity ratio somewhat, as well as accelerate the approach toward Lemlich's limit for low density foam. This effect was attributed to increased suction within the Plateau borders.     

Convective and diffusive surfactant transfer in multiphase liquid systems

Laser interferometry was used to investigate diffusive and convective mass transfer in a multicomponent fluid mixture with a liquid–liquid or liquid–gas interface. For this purpose, an immobile gas bubble or insoluble fluid droplet, having the shape of a short cylinder with a free lateral surface, was inserted into a thin liquid layer. In the case of non-uniform distribution of the dissolved surfactant component, the Marangoni convection near the drop/bubble was initiated by the surface tension inhomogeneities, depending on the surfactant concentration. The applied experimental techniques allowed us to study the structure and evolution of the convective flows and concentration fields in a liquid layer, which due to its small thickness were nearly two-dimensional. Making use of both the vertical and horizontal orientation of the liquid layer, we investigated the mass transfer process at different levels of the interaction between gravity and capillary forces. During the experiments, we detected new solutocapillary phenomena, which were found to be caused by oscillatory regimes of solutal convection occurring around air bubbles and chlorobenzene drops in heterogeneous aqueous solutions of alcohol with a vertical surfactant concentration gradient. The role of the oscillatory instability in the processes of drop saturation by the surfactant from its water solution and an inverse process of surfactant extraction from the drop into the surrounding homogeneous fluid (water) was determined. A reasonable explanation for the driving mechanisms of the discovered effects has been proposed.     

Anomalous motion of viscous fingers in surfactant solutions in a Hele–Shaw cell

Viscous fingering in surfactant solutions in a rectangular Hele–Shaw cell was investigated. Test fluids were aqueous solutions of cetyltrimethylammonium bromide (CTAB) with sodium salicylate (NaSal) as a counter ion, and the ratio of mole concentration of CTAB and that of NaSal was 1–7.7. Two fluids that had a mole concentration different from that of CTAB were used. Air was injected into the cell and the growth of the interface between air and a CTAB/NaSal solution was observed. The fingertip grew similar to the finger growth in shear-thinning fluids at low pressure gradients. It took a cuspidate shape at the intermediate pressure gradient, and a sudden protrusion at a critical shear rate occurred. In high shear rate regions, the finger behaved as in a less shear-thinning fluid. These phenomena relate to rheological properties of the test fluids. Comparison with flow curves for CTAB/NaSal systems showed that the critical shear rate related to the shear rate at which a bending point appeared in the flow curve.     

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.     

Heat transfer characteristics of water and APG surfactant solution in a micro-channel heat sink

The steady increase in internal heat production of cost and high performance electronic components has lead researchers to seek improved ways to remove the heat generated. Single-phase liquid flow has been considered as a potential solution for solving this cooling problem. However, when considering that any solution needs to be of low cost and low mass fluxes and yet retain low temperature gradients across the electronic components, it seems that two-phase boiling flow is preferred. Surfactant solutions have been introduced in connection with enhancement of the boiling processes. We investigated the effects of surfactant solution flows through a micro-channel heat sink. The experimental setup included a high-speed IR radiometer and a CCD camera that were used to characterize the test module. The module consisted of inlet and outlet manifolds that distributed surfactant solutions through an array of 26 parallel micro-channels. The experimental results have shown that there exists an optimal solution concentration and mass flux for enhancing heat removal. Surfactant solution boiling flows were also found to stabilize the maximum and average surface temperatures for a wide range of applied heat fluxes. In addition, the use of surfactant solutions at low mass fluxes has led to CHF enhancement when compared to regular water flows. In the last part of this work, possible explanations for the observed non-ionic surfactant effects are presented.     

Extracting extensional properties through excess pressure drop estimation in axisymmetric contraction and expansion flows for constant shear viscosity,extension strain-hardening fluids

In this study, hyperbolic contraction–expansion flow (HCF) devices have been investigated with the specific aim of devising new experimental measuring systems for extensional rheological properties. To this end, a hyperbolic contraction–expansion configuration has been designed to minimize the influence of shear in the flow. Experiments have been conducted using well-characterized model fluids, alongside simulations using a viscoelastic White–Metzner/FENE-CR model and finite element/finite volume analysis. Here, the application of appropriate rheological models to reproduce quantitative pressure drop predictions for constant shear viscosity fluids has been investigated, in order to extract the relevant extensional properties for the various test fluids in question. Accordingly, experimental evaluation of the hyperbolic contraction–expansion configuration has shown rising corrected pressure drops with increasing elastic behaviour (De=0～16), evidence which has been corroborated through numerical prediction. Moreover, theoretical to predicted solution correspondence has been established between extensional viscosity and first normal stress difference. This leads to a practical means to measure extensional viscosity for elastic fluids, obtained through the derived pressure drop data in these HCF devices.     

Correlations for heat transfer coefficient and pressure drop in the annulus of concentric helical coils

Heat transfer and pressure drop characteristics in the annulus of concentric helical coils heat exchangers were experimentally investigated. The effects of coil curvature ratio, flow configuration, number of turns and addition of surfactant were investigated. Five test coils were designed and manufactured to study the effect of different parameters on heat transfer and pressure drop. The liquids used in the present study were water and oleyl-dihydroxy-etheyl-amine-oxide (ODEAO, C₂₂H₄₅NO₃ = 371) non-ionic aqua surfactant solution flowing through the annulus side. The inner side Reynolds number range 11,000–27,000 and the annulus side range 5,000–19,000. The results showed that the annulus Nusselt number decreases as annulus curvature ratio increases and increases when number of turns decrease. Moreover, the friction factor increases with the curvature ratio and also increases as number of turns decreases. Both Nusselt number and friction factor decrease when ODEAO concentration increases.     

An incompressible SPH method for simulation of unsteady viscoelastic free-surface flows

In this paper, an incompressible smoothed particle hydrodynamics (SPH) method is presented to solve unsteady free-surface flows. Both Newtonian and viscoelastic fluids are considered. In the case of viscoelastic fluids, both the Maxwell and Oldroyd-B models are investigated. The proposed SPH method uses a Poisson pressure equation to satisfy the incompressibility constraints. The solution algorithm is an explicit predictor-corrector scheme and employs an adaptive smoothing length based on density variations. To alleviate the numerical difficulties encountered when fluid is highly stretched, an artificial stress term is incorporated into the momentum equation which reduces the risk of unrealistic fractures in the material. Two challenging test cases, the impacting drop and the jet buckling problems, are solved to demonstrate the capability of the proposed scheme in handling viscoelastic flows with complex free surfaces. The jet buckling test case was solved for a wide range of Weissenberg numbers. It was shown that in all cases the method is stable and fairly accurate and agrees well with the available data.     

The flow behavior of fiber suspensions in Newtonian fluids and polymer solutions.

This is the second part of a study examining the mechanical properties and capillary flow of fiber suspensions in Newtonian fluids and in polymer solutions. In part I results for the viscous and elastic properties of the fiber suspensions were presented and it was shown that the fiber suspensions exhibited normal stresses in Newtonian as well as in viscoelastic suspending media. It was thus expected that circulating secondary flows would occur near the entrance to a capillary. Four types of fillers (glass, carbon, nylon and vinylon fibers) suspended in glycerin, HEC solutions and Separan solutions were investigated. The entrance flow patterns were visualized and the pressure fluctuations measured. The visualization enabled the eddies occurring in the fiber suspensions in Newtonian fluids to be analysed and classified into two tpyes. The results from the flow visualization were correlated with the pressure fluctuations. Empirical equations for the tube length correction factor due to the elasticity were obtained.     

Bagley correction: the effect of contraction angle and its prediction

The excess pressure losses due to end effects (mainly entrance) in the capillary flow of a branched polypropylene melt were studied both experimentally and theoretically. These losses were first determined experimentally as a function of the contraction angle ranging from 10° to 150°. It was found that the excess pressure loss function decreases for the same apparent shear rate with increasing contraction angle from 10° to about 45°, and consequently slightly increases from 45° up to contraction angles of 150°. Numerical simulations using a multimode K-BKZ viscoelastic and a purely viscous (Cross) model were used to predict the end pressures. It was found that the numerical predictions do agree well with the experimental results for small contraction angles up to 30°. However, the numerical simulations under-predict the end pressure for larger contraction angles. The effects of viscoelasticity, shear, and elongation on the numerical predictions are also assessed in detail. Shear is the dominant factor controlling the overall pressure drop in flows through small contraction angles. Elongation becomes important at higher contraction angles (greater than 45°). It is demonstrated in abrupt contractions (angle of 180°) that both the entrance pressure loss and the vortex size are strongly dependent on the extensional viscosity for this branched polymer. It is suggested that such an experiment (visualisation of entrance flow) can be useful in evaluating the validity of constitutive equations and it can also be used to fitting parameters of rheological models that control the elongational viscosity.