Rectangular Microfluidic Channels Research Articles

Homebrew Photolithography for the Rapid and Low-Cost, "Do It Yourself" Prototyping of Microfluidic Devices.

Photolithography is the foundational process at the root of micro-electromechanical (MEMS) and microfluidic systems manufacture. The process is descendant from the semiconductor industry, originating from printed circuit board and microprocessor fabrication, itself historically performed in a cleanroom environment utilizing expensive, specialist microfabrication equipment. Consequently, these conditions prove cost-prohibitive and pose a large barrier to entry. We present a novel homebrew, "do-it-yourself" method for performing photolithography to produce master mold wafers using only household appliances and homemade equipment at the bench side, outside of a cleanroom, producing a range of designs including spiral, serpentine, rectangular, and circulatory. Our homebrew processes result in the production of microfluidic channels with feature resolution of ∼85 μm width and 50 μm height utilizing inkjet-printed photomasks on transparency film to expose dry-film photoresist. From start to finish, the entire process takes under <90 min and costs <£300. With SU8 epoxy negative photoresist and a chrome photomask, our low-cost UV exposure apparatus and homemade spincoater could be used to produce PDMS devices containing large arrays of identical microwells measuring 4.4 μm in diameter. We show that our homebrew method produces both rectangular and spiral microfluidic channels with better results than can be achieved by SLA 3D printing by comparison, and amenable to bonding into multilayer functional microfluidic devices. As these methods are fundamental to microfluidics manufacture, we envision that this work will be of value to researchers across a broad range of disciplines, such as those working in resource-constrained countries or conditions, with many and widely varying applications.

Spatio-temporal programming of lyotropic phase transition in nanoporous microfluidic confinements

HypothesisThe nanoporous polydimethylsiloxane (PDMS) surfaces of a rectangular microfluidic channel, selectively uptakes water molecules, concentrating the solute molecules in an aqueous phase, that could drive phase transitions. Factors such as surface wettability, channel geometry, the surface-to-volume ratio, and surface topography of the confinements could play a key role in tuning the phase transitions spatio-temporally. ExperimentsUsing a lyotropic chromonic liquid crystal as model biological material, confined within nanoporous microfluidic environments, we study molecular assembly driven by nanoporous substrates. By combining timelapse polarized imaging, quantitative image processing, and a simple mathematical model, we analyze the phase transitions and construct a master diagram capturing the role of surface wettability, channel geometry and embedded topography on programmable lyotropic phase transitions. FindingsIntrinsic PDMS nanoporosity and confinement cross-section, together with the imposed wettability regulate the rate of the N-M phase transition; whereas the microfluidic geometry and embedded topography enable phase transition at targeted locations. We harness the emergent long-range order during N-M transition to actuate elasto-advective transport of embedded micro-cargo, demonstrating particle manipulation concepts governed by tunable phase transitions. Our results present a programmable physical route to material assembly in microfluidic environment, and offer a new paradigm for assembling genetic components, biological cargo, and minimal synthetic cells.

Marangoni destabilization of bidimensional-confined gas–liquid co-flowing streams in rectangular microfluidic channels

In microchannels, the stability of a fluid jet injected into another immiscible fluid strongly depends on its degree of geometric confinement. When the width of the jet, w, is larger than the channel height, H, the surface tension driven Rayleigh–Plateau instability is suppressed so that the 2D (bidimensional)-confined jet is absolutely stable and never collapses into bubbles (or drops) in contrast to what occurs when w ≤ H [Dollet et al., “Role of the channel geometry on the bubble pinch-off in flow-focusing,” Phys. Rev. Lett. 100(3), 034504 (2008); Guillot et al., “Stability of a jet in confined pressure-driven biphasic flows at low Reynolds number in various geometries,” Phys. Rev. E 78(1), 016307 (2008)]. We here demonstrate both experimentally and theoretically that this picture is, indeed, no longer valid when Marangoni effects are considered. We experimentally show that the addition of small length alcohol molecules into the liquid phase destabilizes a 2D-confined gas–water microfluidic stream ( w &amp;gt; H), leading to the generation of steady non-linear waves and further to the production of bubbles. Using a simple hydrodynamic model, we show through a linear analysis that the destabilization of the gas stream may result from a Marangoni instability due to the fast adsorption of the alcohol molecules, which occurs on a timescale comparable to that of the microfluidic flow.

Droplet transition from non-axisymmetric to axisymmetric shape: Dynamic role of lubrication film in a rectangular microfluidic channel

In the past 20 years, droplet microfluidics is burgeoning in many chemical and biological applications due to the unique capability of droplets to act as confined containers. Confinement is ensured even in the case of squeezed droplets within microchannels much smaller than droplet volumes due to the presence of a lubrication thin film that prevents contact between droplets and the channel walls. The thickness of the lubrication film depends on the dynamics of the entire microfluidic system, affecting the actual droplet's shape and velocity. Therefore, this film is extensively studied to obtain insight into the dynamics of flowing droplets, especially when confined in small channels. Circular cross section channels are the most studied for their axial symmetry, but practical applications present most likely non-axisymmetric channels, as a result of fabrication processes, such as soft lithographic rectangular channels. The latter showed unique transitional morphological behavior of droplets, which assumes an axisymmetric or non-axisymmetric shape during their flow inside a non-axisymmetric channel, depending on the lubrication film. This work gives a comprehensive experimental characterization of the dynamics of the lubrication film during the droplet shape transition. We settled on a novel approach based on the optical diffraction of a localized light beam provided by two-facing optical waveguides integrated with the microfluidics circuit. The technique allows for studying the dynamics of flowing droplets and their relationship with the lubrication film thickness. Additionally, this experimental system enables a precise definition of two regimes of lubrication film, and the critical capillary number at which the transition occurs.

Elastohydrodynamic Relaxation of Soft and Deformable Microchannels.

Hydrodynamic flows in compliant channels are of great interest in physiology and microfluidics. In these situations, elastohydrodynamic coupling leads to (i)a nonlinear pressure-vs-flow-rate relation, strongly affecting the hydraulic resistance; and (ii), because of the compliance-enabled volume storage, a finite relaxation time under a stepwise change in pressure. This latter effect remains relatively unexplored, even while the timescale can vary over a decade in typical situations. In this study we provide time-resolved measurements of the relaxation dynamics for thin and soft, rectangular microfluidic channels. We describe our data using a perturbative lubrication approximation of the Stokes equation coupled to linear elasticity, while taking into account the effect of compliance and resistance of the entrance. The modeling allows us to completely describe all of the experimental results. Our Letter is relevant for any microfluidic scenario wherein a time-dependent driving is applied and provides a first step in the dynamical description of compliant channel networks.

Effect of Surfactant Dynamics on Flow Patterns Inside Drops Moving in Rectangular Microfluidic Channels

Drops contained in an immiscible liquid phase are attractive as microreactors, enabling sound statistical analysis of reactions performed on ensembles of samples in a microfluidic device. Many applications have specific requirements for the values of local shear stress inside the drops and, thus, knowledge of the flow field is required. This is complicated in commonly used rectangular channels by the flow of the continuous phase in the corners, which also affects the flow inside the drops. In addition, a number of chemical species are present inside the drops, of which some may be surface-active. This work presents a novel experimental study of the flow fields of drops moving in a rectangular microfluidic channel when a surfactant is added to the dispersed phase. Four surfactants with different surface activities are used. Flow fields are measured using Ghost Particle Velocimetry, carried out at different channel depths to account for the 3-D flow structure. It is shown that the effect of the surfactant depends on the characteristic adsorption time. For fast-equilibrating surfactants with a characteristic time scale of adsorption that is much smaller than the characteristic time of surface deformation, this effect is related only to the decrease in interfacial tension, and can be accounted for by the change in capillary number. For slowly equilibrating surfactants, Marangoni stresses accelerate the corner flow, which changes the flow patterns inside the drop considerably.

Thermally-actuated microfluidic membrane valve for point-of-care applications

Microfluidics has enabled low volume biochemistry reactions to be carried out at the point-of-care. A key component in microfluidics is the microfluidic valve. Microfluidic valves are not only useful for directing flow at intersections but also allow mixtures/dilutions to be tuned real-time and even provide peristaltic pumping capabilities. In the transition from chip-in-a-lab to lab-on-a-chip, it is essential to ensure that microfluidic valves are designed to require less peripheral equipment and that they are transportable. In this paper, a thermally-actuated microfluidic valve is presented. The valve itself is fabricated with off-the-shelf components without the need for sophisticated cleanroom techniques. It is shown that multiple valves can be controlled and operated via a power supply and an Arduino microcontroller; an important step towards transportable microfluidic devices capable of carrying out analytical assays at the point-of-care. It is been calculated that a single actuator costs less than $1, this highlights the potential of the presented valve for scaling out. The valve operation is demonstrated by adjusting the ratio of a water/dye mixture in a continuous flow microfluidic chip with Y-junction channel geometry. The power required to operate one microfluidic valve has been characterised both theoretically and experimentally. Cyclical operation of the valve has been demonstrated for 65 h with 585 actuations. The presented valve is capable of actuating rectangular microfluidic channels of 500 μm × 50 μm with an expected temperature increase of up to 5 °C. The fastest actuation times achieved were 2 s for valve closing (heating) and 9 s for valve opening (cooling).

Reducing Taylor-Aris dispersion by exploiting lateral convection associated with acoustic streaming

Flow physicochemical processes, like adsorption/desorption-based separations and chemical conversions in open rectangular microfluidic channels are fundamentally limited by axial dispersion. In this contribution, the implementation of acoustics is investigated to reduce axial dispersion in microchannels. The extent of dispersion reduction is studied by means of plug injections. It was observed that the long range vortex flow of acoustic streaming can efficiently reduce dispersion by more than a factor of two. The observations made in the present study pave the way towards the development of novel, highly efficient analytical or preparative separations and microreactors with uniform residence time distributions. We elaborate on the practical implementations of both areas.

Surface tension driven flow of blood in a rectangular microfluidic channel: Effect of erythrocyte aggregation

Microfluidic platforms have increasingly been explored for in vitro blood diagnostics and for studying complex microvascular processes. The perfusion of blood in such devices is typically achieved through pressure-driven setups. Surface tension driven blood flow provides an alternative flow delivery option, and various studies in the literature have examined the behavior of blood flow in such fluidic devices. In such flows, the influence of red blood cell (RBC) aggregation, the phenomenon majorly responsible for the non-Newtonian nature of blood, requires particular attention. In the present work, we examine differences in the surface tension driven flow of aggregating and non-aggregating RBC and Newtonian suspensions, in a rectangular microchannel. The velocity fields were obtained using micro-PIV techniques. The analytical solution for blood velocity in the channel is developed utilizing the power law model for blood viscosity. The results showed that RBC aggregation has an impact at the late stages of the flow, observed mainly in the bluntness of the velocity profiles. At the initial stages of the flow, the shearing conditions are found moderately elevated, preventing intense RBC aggregate formation. As the flow decelerates in the channel, RBC aggregation increases, affecting the flow characteristics.

Numerical Analysis of Heat Transfer Enhancement in a Micro-Channel Due to Mechanical Stirrers

Abstract Using two-dimensional numerical simulations of the momentum, mass, and energy conservation equations, we investigate the enhancement of heat transfer in a rectangular micro-fluidic channel. The fluid inside the channel is assumed to be stationary initially and actuated by the motion imparted by mechanical stirrers, which are attached to the bottom of the channel. Based on the direction of the oscillation of the stirrers, the boundary conditions can be classified as either no-slip (when the oscillation is perpendicular to the length of the channel) or periodic (when the oscillation is along the length of the channel). The heat transfer enhancement due to the motion of the stirrers (with respect to the stationary stirrer situation) is analyzed in terms of the Reynolds number (ranging from 0.7 to 1000) and the Peclet number (ranging from 10 to 100). We find that the heat transfer first increases and then decreases with an increase in the Reynolds number for any given Peclet number. The heat transferred is maximum at a Reynolds number of 20 for the no-slip case and at a Reynolds number of 40 for the periodic case. For a given Peclet and Reynolds number, the heat flux for the periodic case is always larger than the no-slip case. We explain the reason for these trends using time-averaged flow velocity profiles induced by the oscillation of the mechanical stirrers.

Modeling and Analysis of SOI Gratings-Based Opto-Fluidic Biosensor for Lab-on-a-Chip Applications

The design, modeling, and analysis of a silicon-on-insulator (SOI) grating coupler integrated with a microfluidic channel for lab-on-a-chip applications are presented. The grating coupler was designed to operate at 1310 nm. The simulated SOI structure consisted of a 220 nm top-Si device layer with an integrated waveguide, grating coupler, and a buried oxide layer of 2 µm. A rectangular microfluidic channel was deposited on the SOI optical grating structure for light and fluid interaction. The fluidic flow through the device was driven by centrifugal and Coriolis forces. The grating structure was designed to achieve a maximum coupling efficiency at the optimized injection angle of the light source. The sensitivity of the grating structure could be analyzed and evaluated using the change in coupled power as a function of the effective refractive index and was found to be 0.928 × 10−6 RIU. The SOI optical grating structure along with the micro fluidic channel on top could be effectively used as an absorbance-based lab-on-a-chip biosensor.

Inertial manipulation of bubbles in rectangular microfluidic channels.

Inertial microfluidics is an active field of research that deals with crossflow positioning of the suspended entities in microflows. Until now, the majority of the studies have focused on the behavior of rigid particles in order to provide guidelines for microfluidic applications such as sorting and filtering. Deformable entities such as bubbles and droplets are considered in fewer studies despite their importance in multiphase microflows. In this paper, we show that the trajectory of bubbles flowing in rectangular and square microchannels can be controlled by tuning the balance of forces acting on them. A T-junction geometry is employed to introduce bubbles into a microchannel and analyze their lateral equilibrium position in a range of Reynolds (1 < Re < 40) and capillary numbers (0.1 < Ca < 1). We find that the Reynolds number (Re), the capillary number (Ca), the diameter of the bubble (D[combining macron]), and the aspect ratio of the channel are the influential parameters in this phenomenon. For instance, at high Re, the flow pushes the bubble towards the wall while large Ca or D[combining macron] moves the bubble towards the center. Moreover, in the shallow channels, having aspect ratios higher than one, the bubble moves towards the narrower sidewalls. One important outcome of this study is that the equilibrium position of bubbles in rectangular channels is different from that of solid particles. The experimental observations are in good agreement with the performed numerical simulations and provide insights into the dynamics of bubbles in laminar flows which can be utilized in the design of flow based multiphase flow reactors.

Pure axial flow of viscoelastic fluids in rectangular microchannels under combined effects of electro-osmosis and hydrodynamics

This paper presents an analysis of the combined electro-osmotic and pressure-driven axial flows of viscoelastic fluids in a rectangular microchannel with arbitrary aspect ratios. The rheological behavior of the fluid is described by the complete form of Phan-Thien–Tanner (PTT) model with the Gordon–Schowalter convected derivative which covers the upper convected Maxwell, Johnson–Segalman and FENE-P models. Our numerical simulation is based on the computation of 2D Poisson–Boltzmann, Cauchy momentum and PTT constitutive equations. The solution of these governing nonlinear coupled set of equations is obtained by using the second-order central finite difference method in a non-uniform grid system and is verified against 1D analytical solution of the velocity profile with less than 0.06% relative error. Also, a parametric study is carried out to investigate the effect of channel aspect ratio (width to height), wall zeta potential and the Debye–Huckel parameter on 2D velocity profile, volumetric flow rate and the Poiseuille number in the mixed EO/PD flows of viscoelastic fluids with different Weissenberg numbers. Our results show that, for low channel aspect ratios, the previous 1D analytical models underestimate the velocity profile at the channel half-width centerline in the case of favorable pressure gradients and overestimate it in the case of adverse pressure gradients. The results reveal that the inapplicability of the Debye–Huckel approximation at high zeta potentials is more significant for higher Weissenberg number fluids. Also, it is found that, under the specified values of electrokinetic parameters, there is a threshold for velocity scale ratio in which the Poiseuille number is approximately independent of channel aspect ratio.

Jetting of a shear banding fluid in rectangular ducts.

Non-Newtonian fluids are susceptible to flow instabilities such as shear banding, in which the fluid may exhibit a markedly discontinuous viscosity at a critical stress. Here we report the characteristics and causes of a jetting flow instability of shear banding wormlike micelle solutions in microfluidic channels with rectangular cross sections over an intermediate volumetric flow regime. Particle-tracking methods are used to measure the three-dimensional flow field in channels of differing aspect ratios, sizes, and wall materials. When jetting occurs, it is self-contained within a portion of the channel where the flow velocity is greater than the surroundings. We observe that the instability forms in channels with aspect ratio greater than 5, and that the location of the high-velocity jet appears to be sensitive to stress localizations. Jetting is not observed in a lower concentration solution without shear banding. Simulations using the Johnson-Segalman viscoelastic model show a qualitatively similar behavior to the experimental observations and indicate that compressive normal stresses in the cross-stream directions support the development of the jetting flow. Our results show that nonuniform flow of shear thinning fluids can develop across the wide dimension in rectangular microfluidic channels, with implications for microfluidic rheometry.

Computational study of a magnetic design to improve the diagnosis of malaria: 2D model

This paper investigates the feasibility of a cost effective high gradient magnetic separation based device for the detection and identification of malaria parasites in a blood sample. The design utilizes magnetic properties of hemozoin present in malaria-infected red blood cells (mRBCs) in order to separate and concentrate them inside a microfluidic channel slide for easier examination under the microscope. The design consists of a rectangular microfluidic channel with multiple magnetic wires positioned on top of and underneath it along the length of the channel at a small angle with respect to the channel axis. Strong magnetic field gradients, produced by the wires, exert sufficient magnetic forces on the mRBCs in order to separate and concentrate them in a specific region small enough to fit within the microscope field of view at magnifications typically required to identify the malaria parasite type. The feasibility of the device is studied using a model where the trajectories of the mRBCs inside the channel are determined using first-order ordinary differential equations (ODEs) solved numerically using a multistep ODE solver available within MATLAB. The mRBCs trajectories reveal that it is possible to separate and concentrate the mRBCs in less than 5min, even in cases of very low parasitemia (1–10 parasites/µL of blood) using blood sample volumes of around 3µL employed today.

Electrokinetic and aspect ratio effects on secondary flow of viscoelastic fluids in rectangular microchannels

The secondary flow of PTT fluids in rectangular cross-sectional plane of microchannels under combined effects of electroosmotic and pressure driving forces is the subject of the present study. Employing second-order central finite difference method in a very refined grid network, we investigate the effect of electrokinetic and geometric parameters on the pattern, strength and the average of the secondary flow. In this regard, we try to illustrate the deformations of recirculating vortices due to change in the dimensionless Debye–Huckel and zeta potential parameters as well as channel aspect ratio. We demonstrate that, in the presence of thick electric double layers, significant alteration occurs in the secondary flow pattern by transition from favorable to adverse pressure gradients. Moreover, it is found that for polymer-electrolyte solutions with large Debye lengths, the secondary flow pattern and the shape of vortices are generally dependent upon the width-to-height ratio of the channel cross section. Also, the inspections of strength and average of secondary flow reveal that the sensitivity of these quantities with respect to the electrokinetic, geometric and rheological parameters increases by increasing the absolute value of velocity scale ratio. In this regard, utilizing the curve fitting of the results, several empirical expressions are presented for the strength and average of the secondary flow under various parametric conditions. The obtained relations with the other predictions for secondary flow are of high practical importance when dealing with the design of microfluidic devices that manipulate viscoelastic fluids.

Quasistatic packings of droplets in flat microfluidic channels

As observed in recent experiments, monodisperse droplets self-assemble spontaneously in different ordered packings. In this work, we present a numerical study of the droplet packings in the flat rectangular microfluidic channels. Employing the boundary element method, we numerically solve the Stokes equation in two-dimension and investigate the appearance of droplet packing and transition between one and two-row packings of monodisperse emulsion droplets. By calculating packing force applied on the droplet interface, we investigate the effect of flow rate, droplet size, and surface tension on the packing configurations of droplets and transition between different topological packings.

Flow patterns of liquid–liquid two-phase flow in non-Newtonian fluids in rectangular microchannels

This paper presents an experimental investigation on the flow patterns of Newtonian/non-Newtonian fluids two-phase flow in T-shaped rectangular microchannels by using a high-speed camera. Three different dimensions of microchannels are employed with 400μm×400μm, 400μm×600μm, 400μm×800μm in depth and width, respectively. The cyclohexane is used as the dispersed phase, and various concentrations of carboxyl methyl cellulose (CMC) aqueous solutions as the continuous phase. Results show that slug flow, droplet flow, parallel flow and jet flow are the main flow patterns for cyclohexane-CMC two-phase flow in rectangular microchannels. The larger the concentration of CMC solutions and the smaller the size of the microchannels are, the wider of the regimes are for the droplet flow and parallel flow. A model is proposed to predict the transitions between adjacent flow regimes for liquid–liquid two-phase flow in non-Newtonian fluids in rectangular microchannels.

Pressure effects on electroosmotic flow of power-law fluids in rectangular microchannels

In this paper, the fully developed electroosmotic flow of power-law fluids in rectangular microchannels in the presence of pressure gradient is analyzed. The electrical potential and momentum equations are numerically solved through a finite difference procedure for a non-uniform grid. A complete parametric study reveals that the pressure effects are more pronounced at higher values of the channel aspect ratio and smaller values of the flow behavior index. The Poiseuille number is found to be an increasing function of the channel aspect ratio for pressure assisted flow and a decreasing function of this parameter for pressure opposed flow. It is also observed that the Poiseuille number is increased by increasing the zeta potential. Furthermore, the results show that an increase in the flow behavior index results in a lower flow rate ratio, defined to be the ratio of the flow rate to that of a Newtonian fluid at the same conditions. Moreover, whereas the flow rate ratio in the presence of an opposed pressure gradient is smaller than that of a favorable pressure force for shear thinnings, the opposite is true for shear-thickening fluids.

Circulation within confined droplets in Hele-Shaw channels

Liquid droplets flowing through a rectangular microfluidic channel develop a vortical flow field due to the presence of shear forces from the surrounding fluid. In this paper, we present an experimental and computational study of droplet velocities and internal flow patterns in a rectangular pressure-driven flow for droplet diameters ranging from 0.1 to 2 times the channel height. Our study shows excellent agreement with asymptotic predictions of droplet and interfacial velocities for infinitesimally small droplets. As the droplet diameter nears the size of the channel height, the droplet velocity slows significantly, and the changing external flow field causes a qualitative change in the location of internal vortices. This behavior is relevant for future studies of mass transfer in microfluidic devices.