In this work we created functional microfluidic chips without actually designing them. We accomplished this by first generating a library of thousands of different random microfluidic chip designs, then simulating the behavior of each design on a computer using automated finite element analysis. The simulation results were then saved to a database which a user can query via to find chip designs suitable for a specific task. To demonstrate this functionality, we used our library to select chip designs that generate any three desired concentrations of a solute. We also fabricated and tested 16 chips from the library, confirmed that they function as predicted, and used these chips to perform a cell growth rate assay. This is one of many different applications for randomly-designed microfluidics; in principle, any microfluidic chip that can be simulated could be designed automatically using our method. Using this approach, individuals with no training in microfluidics can obtain custom chip designs for their own unique needs in just a few seconds.

Computer simulation plays a growing role in the design of microfluidic chips. However, the particle tracers in some existing commercial computational fluid dynamics software are not well suited for accurately simulating the trajectories of particles such as cells, microbeads, and droplets in microfluidic systems. To address this issue, we present a microfluidics-optimized particle simulation algorithm (MOPSA) that simulates the trajectories of cells, droplets, and other particles in microfluidic chips with more lifelike results than particle tracers in existing commercial software. When calculating the velocity of a particle, MOPSA treats the particle as a two-dimensional rigid circular object instead of a single point. MOPSA also checks for unrealistic interactions between particles and channel walls and applies an empirical correcting function to eliminate these errors. To validate the performance of MOPSA, we used it to simulate a variety of important features of microfluidic devices like channel intersections and deterministic lateral displacement (DLD) particle sorter chips. MOPSA successfully predicted that different particle sizes will have different trajectories in six published DLD experiments from three research groups; these DLD chips were used to sort a variety of different cells, particles, and droplets. While some of these particles are not actually rigid or spherical, MOPSA's approximation of these particles as rigid spheres nonetheless resulted in lifelike simulations of the behaviors of these particles (at least for the particle sizes and types shown here). In contrast, existing commercial software failed to replicate these experiments. Finally, to demonstrate that MOPSA can be extended to simulate other properties of particles, we added support for simulating particle density to MOPSA and then used MOPSA to simulate the operation of a microfluidic chip capable of sorting cells by their density. By enabling researchers to accurately simulate the behavior of some types of particles in microfluidic chips before fabricating the chips, MOPSA should accelerate the development of new microfluidic devices for important applications.

Microfluidics researchers are increasingly using computer simulation in many different aspects of their research. However, these simulations are often computationally intensive: simulating the behavior of a simple microfluidic chip can take hours to complete on typical computing hardware, and even powerful workstations can lack the computational capabilities needed to simulate more complex chips. This slows the development of new microfluidic chips for new applications. To address this issue, we present a microfluidic simulation method that can simulate the behavior of fluids and particles in some typical microfluidic chips instantaneously (in around one second). Our method decomposes the chip into its primary components: channels and intersections. The behavior of fluid in each channel is determined by leveraging analogies with electronic circuits, and the behavior of fluid and particles in each intersection is determined by querying a database containing nearly 100,000 pre-simulated channel intersections. While constructing this database takes a nontrivial amount of computation time, once built, this database can be queried to determine the behavior of fluids and particles in a given intersection in a fraction of a second. Using this approach, the behavior of a microfluidic chip can be simulated in just one second on a standard laptop computer, without any noticeable degradation in the accuracy of the simulation. While our current technique has some constraints on the designs of the chips it can simulate (namely, T- or cross-shaped intersections, 90 degree channel turns, a fixed channel width, fluid flow rates between 0 and 2 cm/s, and particles with diameters between 1 and 20 microns), we provide several strategies for increasing the range of possible chip designs that can be simulated using our technique. As a proof of concept, we show that our simulation method can instantaneously simulate the paths followed by particles in both simple and complex microfluidic chips, with results that are essentially indistinguishable from simulations that took hours or even days to complete using conventional approaches.

The ability to thoroughly mix two fluids is a fundamental need in microfluidics. While a variety of different microfluidic mixers have been designed by researchers, it remains unknown which (if any) of these mixers are optimal (that is, which designs provide the most thorough mixing with the smallest possible fluidic resistance across the mixer). In this work, we automatically designed and rationally optimized a microfluidic mixer. We accomplished this by first generating a library of thousands of different randomly designed mixers, then using the non-dominated sorting genetic algorithm II (NSGA-II) to optimize the random chips in order to achieve Pareto efficiency. Pareto efficiency is a state of allocation of resources (e.g. driving force) from which it is impossible to reallocate so as to make any one individual criterion better off (e.g. pressure drop) without making at least one individual criterion (e.g. mixing performance) worse off. After 200 generations of evolution, Pareto efficiency was achieved and the Pareto-optimal front was found. We examined designs at the Pareto-optimal front and found several design criteria that enhance the mixing performance of a mixer while minimizing its fluidic resistance; these observations provide new criteria on how to design optimal microfluidic mixers. Additionally, we compared the designs from NSGA-II with some popular microfluidic mixer designs from the literature and found that designs from NSGA-II have lower fluidic resistance with similar mixing performance. As a proof of concept, we fabricated three mixer designs from 200 generations of evolution and one conventional popular mixer design and tested the performance of these four mixers. Using this approach, an optimal design of a passive microfluidic mixer is found and the criteria of designing a passive microfluidic mixer are established.