Multiphase flow in microfluidic systems --control and applications of droplets and interfaces

Micro- and nanotechnology can provide us with many tools for the production, study and detection of colloidal and interfacial systems. In multiphase flow in micro- and nanochannels several immiscible fluids will be separated from each other by flexible fluidic interfaces. The multiphase coexistence and the small-volume confinement provide many attractive characteristics. Multiphase flow in microfluidic systems shows a complicated behavior but has many practical uses compared to a single-phase flow. In this paper, we discuss the methods of controlling multiphase flow to generate either micro- or nano-droplets (or bubbles) or stable stratified interfaces between fluidic phases. Furthermore, applications of the droplets and interfaces in microchannels are summarized.     

Deformation of isotropic and anisotropic liquid droplets dispersed in a cellulose liquid crystalline derivative

In this work we study the strain-induced deformation of both isotropic and anisotropic liquid droplets dispersed in a liquid crystalline cellulose matrix. We have produced two types of acetoxypropylcellulose (APC) solid films one with a droplet dispersion of the commercial liquid crystal E7 from Merck, and another with a droplet dispersion of silicone oil. To produce the solid films a solution of APC (60%wt) in dimethylacetamide (DMAc) with 15%wt of either the commercial nematic liquid crystal E7 or the silicone oil was prepared. After homogenization the phase separated solutions were submitted to a shear flow mechanical field and casted onto a Teflon plate. We performed mechanical uniaxial stress-strain tests in the free standing films recording continuously the strain and images of the deformed droplets. The mechanical tests were carried out using a mini stress-strain testing machine apparatus and an Olympus optical polarizing microscope with an attached CCD camera. The images obtained from the mechanical tests for each value of the imposed strain were then analyzed comparing the images of deformed droplets with those of the undeformed ones, extracting in this way the local strain field. The droplet deformation data obtained show, as expected, significant differences in the local strain field when stretching parallel and perpendicular to the initial shear direction. No significant differences were found in the local strain fields obtained from the silicone oil and the E7 droplets films. The local strain fields variation with the imposed strain are compared with the predictions of the theory developed for nematic elastomers by Warner and Terentjev (Liquid crystal elastomers. Clarendon Oxford Press, Oxford, 2003).     

Burning rates of turbulent iso-octane aerosol mixtures in spherical flame explosions

Experimental studies of aerosol combustion under quiescent and turbulence conditions have been conducted to quantify the differences in the flame structure and burning rates between aerosol and gaseous mixtures. Turbulence was generated by variable speed fans to yield rms turbulence velocities between 0.5 and 4.0 m/s and this was uniform and isotropic. Homogeneously distributed and near monodispersed iso-octane-air aerosol clouds were generated using a thermodynamic condensation method. Spherically expanding flames, following central ignition, at near atmospheric pressures were employed to quantify the flame structure and propagation rate. The effects of the diameter of fine fuel droplets on flame propagation were investigated. It is suggested that the inertia of fuel droplets is an important cause of flame enhancement during early flame development. During later stages, cellular flame instability and the effective, gaseous phase, equivalence ratio becomes important. The latter effect leads has increases the flame speed of rich mixtures, but decreases that of lean ones. Droplet enhancement of burning velocity can be significant at low turbulence but is negligible at high turbulence.     

Al2O3:C optically stimulated luminescence droplets: Characterization and applications in medical beams

Optically stimulated luminescence (OSL) detectors, which are widely used in radiation protection, offer a number of potential advantages for radiotherapy dosimetry. In this study we characterized 1-μl of OSL droplets consisting of a mixture of Al₂O₃:C powder and a photo-curable polymer, in addition to results described in a previous work (Nascimento et al., 2013). The concentration test showed that droplets have a higher spatial resolution than other commonly used Al₂O₃:C-based detectors. Our results from the dose response, reproducibility and dependence with accumulative dose were obtained for droplets with a powder/polymer concentration that showed a high Signal to Noise Ratio (SNR) without compromising the droplet malleability. Additional test results show the response of such droplets in percentage depth dose curves and dose profiles of clinical beams.     

Several surfaces with special wettability: Influence on spreading and motion of W/O emulsion droplets

Physical and chemical modifications were made on the surface of the aluminum sheet to change the surface properties and superhydrophobic–hydrophilic wettability gradient surface was made on the perspex surface by using microstructure-pattering technique and self-assembled-monolayer method. By using high-speed video camera system and optical tensiometer, this paper discusses the influence of special surfaces with different wettability on spreading and motion of water, oil, and W/O emulsion droplets both experimentally and theoretically. In addition, the paper also discusses the influence of the superhydrophobic–hydrophilic wettability gradient on fluidity of W/O emulsion droplets and the coalescence process of droplets. The results showed that the contact angle of W/O emulsion droplets on the modified surfaces was related to the water and oil distribution at the three-phase line. On the wettability gradient surface, the droplet moved spontaneously when the droplet was located at the junction of the gradient. A quasi-steady theoretical model was used to analyze the driving and resistant forces acting on a droplet to improve the understanding of the self-transport behavior of the droplets.     

Modeling anisotropic Reynolds-stress dissipation in particle- or droplet-laden flows

Particles or droplets dispersed in turbulent flows at sufficiently high volume loadings lead to modifications of turbulence characteristics. More specifically, in a detailed experimental investigation by Poelma and coworkers, where the particle phase moves with a non-zero mean velocity relative to the fluid phase, it was found that anisotropic Reynolds-stress dissipation is induced. Recently, we have proposed a model that can account for this effect in RANS-based and PDF method simulations. In our previous work, however, no simulation results of the RANS model formulation were presented. In the present work, a new compact tensorial RANS formulation is presented and the new formulation is validated against the experimental data of Poelma and coworkers.     

Self-aeration and turbulence in a stepped channel: Influence of cavity surface roughness

The strong interactions between free-surface flows and atmospheric surroundings may lead to substantial air–water mixing with void fractions ranging from zero in clear-water to 100%. In this study, the air–water flow properties were studied in a large stepped water channel operating at large Reynolds numbers. Interactions between free-surface and cavity recirculation were systematically investigated in the skimming flow regime. Some surface roughness was introduced on the cavity walls and identical experiments were performed with several configurations. Basic results demonstrated some influence of step surface roughness on the flow properties leading to some counter-intuitive finding. The presence of cavity roughness was associated with higher flow velocities and comparatively lower turbulence levels. Distributions of bubble/droplet chords spanned over several orders of magnitude without significant influence of the cavity roughness. The distributions of turbulence levels and bubble count rates showed some correlation and highlighted strong interactions between entrained particles (bubbles, drops) and the flow turbulence.     

Numerical investigation of condensing steam flow in boundary layers

The paper describes a numerical method for the prediction of condensing steam flow within compressible boundary layers. The method is based on a simple stream function technique, which enables straightforward integration of the nucleation and droplet growth equations in a Lagrangian frame of reference. Calculations show how viscous dissipation and reduced expansion rate within a typical boundary layer influence nucleation and growth, leading to droplet radii and size distributions that differ substantially from those predicted in inviscid flow. The impact of condensation on temperature and velocity profiles, and the implications for thermodynamic loss are also considered.     

A microfluidic chip capable of switching W/O droplets to vertical laminar flow for electrochemical detection of droplet contents

Analysis of droplet contents is a key function involved in droplet-based microfluidic systems. Direct electrochemical detection of droplet contents suffers problems such as relatively poor repeatability, interference of capacitive current and relatively poor detectability. This paper presents a novel hybrid polydimethylsiloxane-glass chip for highly sensitive and reproducible amperometric detection of droplet contents. By wettability-patterning of the channel surface of the hybrid chip, water in oil droplets generated in the upstream part of the central channel can be switched to a two-phase vertical laminar flow (i.e., a continuous oil stream flowing atop a continuous aqueous stream) in the downstream part of the channel. The vertical laminar flow keeps the analyte in the underneath-flowing aqueous stream in direct contact with the sensing electrodes located on the bottom surface of the channel. Therefore, steady-state current signals with high sensitivity (1.2 A M⁻¹ cm⁻² for H₂O₂), low limit of detection (0.12 μM, S/N = 2), and good reproducibility (RSD 1.1% at 0.3 mM H₂O₂) were obtained. The methods for patterning of the inner channel surface are presented, and the behaviors of the microchip in flow profile switching and amperometric detection are discussed. The application of the developed microchip to enzyme kinetics study is also demonstrated.     

Large-Eddy Simulation of Interactions Between a Reacting Jet and Evaporating Droplets

Large-eddy simulation of a turbulent reactive jet with and without evaporating droplets is performed to investigate the interactions among turbulence, combustion, heat transfer and evaporation. A hybrid Eulerian–Lagrangian approach is used for the gas–liquid flow system. Arrhenius-type finite-rate chemistry is employed for the chemical reaction. To capture the highly local interactions, dynamic procedures are used for all the subgrid-scale models, except that the filtered reaction rate is modelled by a scale similarity model. Various representative cases with different initial droplet sizes (St ₀) and mass loading ratios (MLR) have been simulated, along with a case without droplets. It is found that compared with the bigger, slow responding droplets (St ₀ = 16), smaller droplets (St ₀ = 1) are more efficient in suppressing combustion due to their preferential concentration in the reaction zones. The peak temperature and intensity of temperature fluctuations are found to be reduced in all the droplet cases, to a varying extent depending on the droplet properties. Detailed analysis on the contributions of respective terms in a transport equation for grid-scale kinetic energy (GSKE) shows that the droplet evaporation effect on GSKE is small, while the droplet momentum effect depends on St ₀. When the MLR is sufficiently high, the bigger (St ₀ = 16) droplets can have profound influence on GSKE, and consequently on the formation and evolution of large-scale flow structures. On the other hand, the turbulence level is found to be lower in the droplet cases than in the pure flame case, due to the dissipative droplet dynamic effect.