The application of microfluidic devices for DNA amplification has recently been extensively studied. Here, we review the important development of microfluidic polymerase chain reaction (PCR) devices and discuss the underlying physical principles for the optimal design and operation of the device. In particular, we focus on continuous-flow microfluidic PCR on-chip, which can be readily implemented as an integrated function of a micro-total-analysis system. To overcome sample carryover contamination and surface adsorption associated with microfluidic PCR, microdroplet technology has recently been utilized to perform PCR in droplets, which can eliminate the synthesis of short chimeric products, shorten thermal-cycling time, and offers great potential for single DNA molecule and single-cell amplification. The work on chip-based PCR in droplets is highlighted.

Using a phase-field model to describe fluid/fluid interfacial dynamics and a lattice Boltzmann model to address hydrodynamics, two dimensional (2D) numerical simulations have been performed to understand the mechanisms of droplet formation in microfluidic T-juntion. Although 2D simulations may not capture underlying physics quantitatively, our findings will help to clarify controversial experimental observations and identify new physical mechanisms. We have systematically examined the influence of capillary number, flow rate ratio, viscosity ratio, and contact angle in the droplet generation process. We clearly observe that the transition from the squeezing regime to the dripping regime occurs at a critical capillary number of 0.018, which is independent of flow rate ratio, viscosity ratio, and contact angle. In the squeezing regime, the squeezing pressure plays a dominant role in the droplet breakup process, which arises when the emerging interface obstructs the main channel. The droplet size depends on both the capillary number and the flow rate ratio, but is independent of the viscosity ratio under completely hydrophobic wetting conditions. In the dripping regime, the droplet size will be significantly influenced by the viscosity ratio as well as the built-up squeezing pressure. When the capillary number increases, the droplet size becomes less dependent on the flow rate ratio. The contact angle also affects the droplet shape, size, and detachment point, especially at small capillary numbers. More hydrophobic wetting properties are expected to produce smaller droplets. Interestingly, the droplet size is dependent on the viscosity ratio in the squeezing regime for less hydrophobic wetting conditions

The Boltzmann equation describes the dynamics of rarefied gas flows, but the multidimensional nature of its collision operator poses a real challenge for its numerical solution. In this paper, the fast spectral method [36], originally developed by Mouhot and Pareschi for the numerical approximation of the collision operator, is extended to deal with other collision kernels, such as those corresponding to the soft, Lennard–Jones, and rigid attracting potentials. The accuracy of the fast spectral method is checked by comparing our numerical solutions of the space-homogeneous Boltzmann equation with the exact Bobylev–Krook–Wu solutions for a gas of Maxwell molecules. It is found that the accuracy is improved by replacing the trapezoidal rule with Gauss–Legendre quadrature in the calculation of the kernel mode, and the conservation of momentum and energy are ensured by the Lagrangian multiplier method without loss of spectral accuracy. The relax-to-equilibrium processes of different collision kernels with the same value of shear viscosity are then compared; the numerical results indicate that different forms of the collision kernels can be used as long as the shear viscosity (not only the value, but also its temperature dependence) is recovered. An iteration scheme is employed to obtain stationary solutions of the space-inhomogeneous Boltzmann equation, where the numerical errors decay exponentially. Four classical benchmarking problems are investigated: the normal shock wave, and the planar Fourier/Couette/force-driven Poiseuille flows. For normal shock waves, our numerical results are compared with a finite difference solution of the Boltzmann equation for hard sphere molecules, experimental data, and molecular dynamics simulation of argon using the realistic Lennard–Jones potential. For planar Fourier/Couette/force-driven Poiseuille flows, our results are compared with the direct simulation Monte Carlo method. Excellent agreements are observed in all test cases, demonstrating the merit of the fast spectral method as a computationally efficient method for rarefied gas dynamics

Based on the fast spectral approximation to the Boltzmann collision operator, we present an accurate and efficient deterministic numerical method for solving the Boltzmann equation. First, the linearized Boltzmann equation is solved for Poiseuille and thermal creep flows, where the influence of different molecular models on the mass and heat flow rates is assessed, and the Onsager-Casimir relation at the microscopic level for large Knudsen numbers is demonstrated. Recent experimental measurements of mass flow rates along a rectangular tube with large aspect ratio are compared with numerical results for the linearized Boltzmann equation. Then, a number of two-dimensional microflows in the transition and free-molecular flow regimes are simulated using the nonlinear Boltzmann equation. The influence of the molecular model is discussed, as well as the applicability of the linearized Boltzmann equation. For thermally driven flows in the free-molecular regime, it is found that the magnitudes of the flow velocity are inversely proportional to the Knudsen number. The streamline patterns of thermal creep flow inside a closed rectangular channel are analysed in detail: when the Knudsen number is smaller than a critical value, the flow pattern can be predicted based on a linear superposition of the velocity profiles of linearized Poiseuille and thermal creep flows between parallel plates. For large Knudsen numbers, the flow pattern can be determined using the linearized Poiseuille and thermal creep velocity profiles at the critical Knudsen number. The critical Knudsen number is found to be related to the aspect ratio of the rectangular channel.

Using a lattice Boltzmann multiphase model, three-dimensional numerical simulations have been performed to understand droplet formation in microfluidic cross-junctions at low capillary numbers. Flow regimes, consequence of interaction between two immiscible fluids, are found to be dependent on the capillary number and flow rates of the continuous and dispersed phases. A regime map is created to describe the transition from droplets formation at a cross-junction (DCJ), downstream of cross-junction to stable parallel flows. The influence of flow rate ratio, capillary number, and channel geometry is then systematically studied in the squeezing-pressure-dominated DCJ regime. The plug length is found to exhibit a linear dependence on the flow rate ratio and obey power-law behavior on the capillary number. The channel geometry plays an important role in droplet breakup process. A scaling model is proposed to predict the plug length in the DCJ regime with the fitting constants depending on the geometrical parameters.