Agent-based simulations of complex droplet pattern formation in a two-branch microfluidic network

Non-Newtonian droplet-based microfluidics logic gates

Droplet-based microfluidic logic gates have many applications in diagnostic assays and biosciences due to their automation and the ability to be cascaded. In spite of many bio-fluids, such as blood exhibit non-Newtonian characteristics, all the previous studies have been concerned with the Newtonian fluids. Moreover, none of the previous studies has investigated the operating regions of the logic gates. In this research, we consider a typical AND/OR logic gate with a power-law fluid. We study the effects of important parameters such as the power-law index, the droplet length, the capillary number, and the geometrical parameters of the microfluidic system on the operating regions of the system. The results indicate that AND/OR states mechanism function in opposite directions. By increasing the droplet length, the capillary number and the power-law index, the operating region of AND state increases while the operating region of OR state reduces. Increasing the channel width will decrease the operating region of AND state while it increases the operating region of OR state. For proper operation of the logic gate, it should work in both AND/OR states appropriately. By combining the operating regions of these two states, the overall operating region of the logic gate is achieved.

Simulation before fabrication: a case study on the utilization of simulators for the design of droplet microfluidic networks

The functional performance of passively operated droplet microfluidics is sensitive with respect to the dimensions of the channel network, the fabrication precision as well as the applied pressure because the entire network is coupled together. Especially, the local and global hydrodynamic resistance changes caused by droplets make the task to develop a robust microfluidic design challenging as plenty of interdependencies which all affect the intended behavior have to be considered by the designer. After the design, its functionality is usually validated by fabricating a prototype and testing it with physical experiments. In case that the functionality is not implemented as desired, the designer has to go back, revise the design, and repeat the fabrication as well as experiments. This current design process based on multiple iterations of refining and testing the design produces high costs (financially as well as in terms of time). In this work, we show how a significant amount of those costs can be avoided when applying simulation before fabrication. To this end, we demonstrate how simulations on the 1D circuit analysis model can help in the design process by means of a case study. Therefore, we compare the design process with and without using simulation. As a case study, we use a microfluidic network which is capable of trapping and merging droplets with different content on demand. The case study demonstrates how simulation can help to validate the derived design by considering all local and global hydrodynamic resistance changes. Moreover, the simulations even allow further exploration of different designs which have not been considered before due to the high costs.

Simulating microfluidic networks allows to check a design even before first prototypes are realized.

Microscale distribution and dynamic surface tension of pulmonary surfactant normalize the recruitment of asymmetric bifurcating airways

We investigate the influence of bifurcation geometry, asymmetry of daughter airways, surfactant distribution, and physicochemical properties on the uniformity of airway recruitment of asymmetric bifurcating airways. To do so, we developed microfluidic idealized in vitro models of bifurcating airways, through which we can independently evaluate the impact of carina location and daughter airway width and length. We explore the uniformity of recruitment and its relationship to the dynamic surface tension of the lining fluid and relate this behavior to the hydraulic (PHyd) and capillary (PCap) pressure drops. These studies demonstrate the extraordinary importance of PCap in stabilizing reopening, even in highly asymmetric systems. The dynamic surface tension of pulmonary surfactant is integral to this stability because it modulates PCap in a velocity-dependent manner. Furthermore, the surfactant distribution at the propagating interface can have a very large influence on recruitment stability by focusing surfactant preferentially to specific daughter airways. This implies that modification of the surfactant distribution through novel modes of ventilation could be useful in inducing uniformly recruited lungs, aiding in gas exchange, and reducing ventilator-induced lung injury.

NEW & NOTEWORTHY The dynamic surface tension of pulmonary surfactant is integral to the uniformity of asymmetric bifurcation airway recruitments because it modulates capillary pressure drop in a velocity-dependent manner. Also, the surfactant distribution at the propagating interface can have a very large influence on recruitment stability by focusing surfactant preferentially to specific daughter airways. This implies that modification of the surfactant distribution through novel modes of ventilation could be useful in inducing uniformly recruited lungs, reducing ventilator-induced lung injury.

Between giant oscillations and uniform distribution of droplets: The role of varying lumen of channels in microfluidic networks

The simplest microfluidic network (a loop) comprises two parallel channels with a common inlet and a common outlet. Recent studies that assumed a constant cross section of the channels along their length have shown that the sequence of droplets entering the left (L) or right (R) arm of the loop can present either a uniform distribution of choices (e.g., RLRLRL...) or long sequences of repeated choices (RRR...LLL), with all the intermediate permutations being dynamically equivalent and virtually equally probable to be observed. We use experiments and computer simulations to show that even small variation of the cross section along channels completely shifts the dynamics either into the strong preference for highly grouped patterns (RRR...LLL) that generate system-size oscillations in flow or just the opposite-to patterns that distribute the droplets homogeneously between the arms of the loop. We also show the importance of noise in the process of self-organization of the spatiotemporal patterns of droplets. Our results provide guidelines for rational design of systems that reproducibly produce either grouped or homogeneous sequences of droplets flowing in microfluidic networks.

Millifluidics as a simple tool to optimize droplet networks: Case study on drop traffic in a bifurcated loop

We report that modular millifluidic networks are simpler, more cost-effective alternatives to traditional microfluidic networks, and they can be rapidly generated and altered to optimize designs. Droplet traffic can also be studied more conveniently and inexpensively at the millimeter scale, as droplets are readily visible to the naked eye. Bifurcated loops, ladder networks, and parking networks were made using only Tygon(®) tubing and plastic T-junction fittings and visualized using an iPod(®) camera. As a case study, droplet traffic experiments through a millifluidic bifurcated loop were conducted, and the periodicity of drop spacing at the outlet was mapped over a wide range of inlet drop spacing. We observed periodic, intermittent, and aperiodic behaviors depending on the inlet drop spacing. The experimentally observed periodic behaviors were in good agreement with numerical simulations based on the simple network model. Our experiments further identified three main sources of intermittency between different periodic and/or aperiodic behaviors: (1) simultaneous entering and exiting events, (2) channel defects, and (3) equal or nearly equal hydrodynamic resistances in both sides of the bifurcated loop. In cases of simultaneous events and/or channel defects, the range of input spacings where intermittent behaviors are observed depends on the degree of inherent variation in input spacing. Finally, using a time scale analysis of syringe pump fluctuations and experiment observation times, we find that in most cases, more consistent results can be generated in experiments conducted at the millimeter scale than those conducted at the micrometer scale. Thus, millifluidic networks offer a simple means to probe collective interactions due to drop traffic and optimize network geometry to engineer passive devices for biological and material analysis.

Origin of periodic and chaotic dynamics due to drops moving in a microfluidic loop device

Droplets moving in a microfluidic loop device exhibit both periodic and chaotic behaviors based on the inlet droplet spacing. We observe that the periodic behavior is an outcome of carrier phase mass conservation principle, which translates into a droplet spacing quantization rule. This rule implies that the summation of exit spacing is equal to an integral multiple of inlet spacing. This principle also enables identification of periodicity in experimental systems with input scatter. We find that the origin of chaotic behavior is through intermittency, which arises when drops enter and leave the junctions at the same time. We derive an analytical expression to estimate the occurrence of these chaotic regions as a function of system parameters. We provide experimental, simulation, and analytical results to validate the origin of periodic and chaotic behavior.

Droplet sorting in a loop of flat microfluidic channels

Motivated by recent experiments, we numerically study the droplet traffic in microfluidic channels forming an asymmetric loop with a long and a short arm. The loop is connected to an inlet and an outlet channel by two right angled T-junctions. Assuming flat channels, we employ the boundary element method (BEM) to numerically solve the two-dimensional Darcy equation that governs two phase flow in the Hele-Shaw limit. The occurrence of different sorting regimes is summarized in sorting diagrams in terms of droplet size, distance between consecutive droplets in the inlet channel, and loop asymmetry for mobility ratios of the liquid phases larger and smaller than one. For large droplet distances, the traffic is regulated by the ratio of the total hydraulic resistances of the long and short arms. At high droplet densities and below a critical droplet size, droplet-droplet collisions are observed for both mobility ratios.

Path selection rules for droplet trains in single-lane microfluidic networks

We investigate the transport of periodic trains of droplets through microfluidic networks having one inlet, one outlet, and nodes consisting of T junctions. Variations of the dilution of the trains, i.e., the distance between drops, reveal the existence of various hydrodynamic regimes characterized by the number of preferential paths taken by the drops. As the dilution increases, this number continuously decreases until only one path remains explored. Building on a continuous approach used to treat droplet traffic through a single asymmetric loop, we determine selection rules for the paths taken by the drops and we predict the variations of the fraction of droplets taking these paths with the parameters at play including the dilution. Our results show that as dilution decreases, the paths are selected according to the ascending order of their hydrodynamic resistance in the absence of droplets. The dynamics of these systems controlled by time-delayed feedback is complex: We observe a succession of periodic regimes separated by a wealth of bifurcations as the dilution is varied. In contrast to droplet traffic in single asymmetric loops, the dynamical behavior in networks of loops is sensitive to initial conditions because of extra degrees of freedom.

Filtering microfluidic bubble trains at a symmetric junction

We report how a nominally symmetric microfluidic junction can be used to sort all bubbles of an incoming train exclusively into one of its arms. The existence of this "filter" regime is unexpected, given that the junction is symmetric. We analyze this behavior by quantifying how bubbles modulate the hydrodynamic resistance in microchannels and show how speeding up a bubble train whilst preserving its spatial periodicity can lead to filtering at a nominally symmetric junction. We further show how such an asymmetric traffic of bubble trains can be triggered in symmetric geometries by identifying conditions wherein the resistance to flow decreases with an increase in the number of bubbles in the microchannel and derive an exact criterion to predict the same.

Design of pressure-driven microfluidic networks using electric circuit analogy

This article reviews the application of electric circuit methods for the analysis of pressure-driven microfluidic networks with an emphasis on concentration- and flow-dependent systems. The application of circuit methods to microfluidics is based on the analogous behaviour of hydraulic and electric circuits with correlations of pressure to voltage, volumetric flow rate to current, and hydraulic to electric resistance. Circuit analysis enables rapid predictions of pressure-driven laminar flow in microchannels and is very useful for designing complex microfluidic networks in advance of fabrication. This article provides a comprehensive overview of the physics of pressure-driven laminar flow, the formal analogy between electric and hydraulic circuits, applications of circuit theory to microfluidic network-based devices, recent development and applications of concentration- and flow-dependent microfluidic networks, and promising future applications. The lab-on-a-chip (LOC) and microfluidics community will gain insightful ideas and practical design strategies for developing unique microfluidic network-based devices to address a broad range of biological, chemical, pharmaceutical, and other scientific and technical challenges.

Passive droplet trafficking at microfluidic junctions under geometric and flow asymmetries

When droplets enter a junction they sort to the channel with the highest flow rate at that instant. Transport is regulated by a discrete time-delayed feedback that results in a highly periodic behavior where specific patterns can continue to cycle indefinitely. Between these highly ordered regimes are chaotic structures where no pattern is evident. Here we develop a model that describes droplet sorting under various asymmetries: branch geometry (length, cross-section), droplet resistance and pressures. First, a model is developed based on the continuum assumption and then, with the assistance of numerical simulations, a discrete model is derived to predict the length and composition of the sorting pattern. Furthermore we derive all unique sequences that are possible for a given distribution and develop a preliminary estimation of why chaotic regimes form. The model is validated by comparing it to numerical simulations and results from microfluidic experiments in PDMS chips with good agreement.

Sampling from nanoliter plugs via asymmetrical splitting of segmented flow

A microfluidic device is developed to sample from oil-segmented fluid plugs of nanoliter volume by passive splitting. In the device, plugs are pumped into a loop where they split according to the flow resistance of each arm of the loop. The loop structure prevents changes in flow resistance over time, caused by accumulation of plugs in downstream collector capillaries, from interfering with the split ratio. To prevent plugs from recombining at the downstream junction of the loop, a series of posts is fabricated into the junction. This structure allows oil to cross, enabling pressure equalization, but not the aqueous plugs allowing them to be collected into separate capillaries or channels. The split ratio depended on both dimensions of the loop channels and frequency of plugs entering the loop. Split ratios from 1:1 to 34:1 were achieved for samples from 1.7 to 3.3 nL. Long-term stability was demonstrated by splitting over 7000 plugs with 6.3% relative standard deviation of the daughter plug size. The system will have application to high-throughput chemical analysis on a nanoliter scale by enabling several applications including performing multiple assays on a single sample, preserving sample while some is removed for analysis, and splitting reagents or test compounds for use in multiple assays.

Transport of resistance through a long microfluidic channel

We introduce a continuous analytical model of propagation of resistance in pressure-driven flow of two-phase fluid in a single channel. This model can be used to predict and interpret experimental results in droplet microfluidics where the hydrodynamic resistance of a capillary comprises a constant part, specific to the channel and the viscosity of the continuous fluid, and a variable part, related to the presence and distribution of droplets. The continuous model is a convenient generalization of the discrete models as demonstrated via comparisons with discrete simulations and with experiments.