A serpentine laminating micromixer combining splitting/recombination and advection

Mixing enhancement has drawn great attention from designers of micromixers, since the flow in a microchannel is usually characterized by a low Reynolds number (Re) which makes the mixing quite a difficult task to accomplish. In this paper, a novel integrated efficient micromixer named serpentine laminating micromixer (SLM) has been designed, simulated, fabricated and fully characterized. In the SLM, a high level of efficient mixing can be achieved by combining two general chaotic mixing mechanisms: splitting/recombination and chaotic advection. The splitting and recombination (in other terms, lamination) mechanism is obtained by the successive arrangement of "F"-shape mixing units in two layers. The advection is induced by the overall three-dimensional serpentine path of the microchannel. The SLM was realized by SU-8 photolithography, nickel electroplating, injection molding and thermal bonding. Mixing performance of the SLM was fully characterized numerically and experimentally. The numerical mixing simulations show that the advection acts favorably to realize the ideal vertical lamination of fluid flow. The mixing experiments based on an average mixing color intensity change of phenolphthalein show a high level of mixing performance was obtained with the SLM. Numerical and experimental results confirm that efficient mixing is successfully achieved from the SLM over the wide range of Re. Due to the simple and mass producible geometry of the efficient micromixer, SLM proposed in this study, the SLM can be easily applied to integrated microfluidic systems, such as micro-total-analysis-systems or lab-on-a-chip systems.     

Design and evaluation of a Dean vortex-based micromixer

A mixer, based on the Dean vortex, is fabricated and tested in an on-chip format. When fluid is directed around a curve under pressure driven flow, the high velocity streams in the center of the channel experience a greater centripetal force and so are deflected outward. This creates a pair of counter-rotating vortices moving fluid toward the inner wall at the top and bottom of the channel and toward the outer wall in the center. For the geometries studied, the vortices were first seen at Reynolds numbers between 1 and 10 and became stronger as the flow velocity is increased. Vortex formation was monitored in channels with depth/width ratios of 0.5, 1.0, and 2.0. The lowest aspect ratio strongly suppressed vortex formation. Increasing the aspect ratio above 1 appeared to provide improved mixing. This design has the advantages of easy fabrication and low surface area.     

CFD simulation and analysis of emulsion droplet formation from straight-through microchannels

We recently proposed a technique for preparing monodisperse emulsions with a coefficient of variation below 5% from a silicon array of micrometer-sized channels perpendicular to the plate surface, named a straight-through microchannel (MC). This study involved three-dimensional computational fluid dynamics (CFD) simulations to calculate the formation of an oil-in-water (O/W) emulsion droplet from straight-through MCs with circular and elliptic cross sections. The CFD results demonstrated that the oil phase that passed through the elliptic MCs exceeding a threshold aspect ratio between 3 and 3.5 was cut off spontaneously into a small droplet with a diameter of approximately 40 microm. Sufficient space for water at the channel exit had to be maintained for successful droplet formation. The formation and shrinkage of a neck inside the channel caused an increased pressure difference inside the channel and an increased velocity value near the neck. The pressure and velocity values at the neck drastically changed, and the neck was cut off instantaneously just before the completion of droplet formation. This process was triggered by a gradually increased pressure difference between the circular neck and inflating oil phase. The findings obtained in this paper provide useful numerical and visual information about the droplet formation phenomena from the straight-through MCs. The CFD results were verified by the experimental results, showing that the CFD approach can help design a suitable channel structure.     

格子Boltzman方法数值模拟微通道内不互溶液液两相流动

利用基于场介子格子Boltzman(LB)方法对T形微通道内油-水不互溶两相流动进行了数值模拟.通过静止液滴算例改进该两相流LB模型且验证了模型的有效性,成功数值模拟了T形微通道内油水两相流的五种不同流型.数值模拟结果和实验数据定性和定量吻合很好,表明本文数值方法的有效性.此外,采用改进后的LB两相流模型模拟了不同界面张力和接触角条件下的液滴形成过程,讨论了这两个变量对流型和液滴形状的影响.结果表明LB方法是研究微通道内不互溶两相流的一种有效的数值工具,能够对各种涉及不互溶多相流微流体设备的设计提供指导.     

Numerical simulation of immiscible liquid-liquid flow in microchannels using lattice Boltzmann method

Immiscible kerosene-water two-phase flows in microchannels connected by a T-junction were numerically studied by a Lattice Boltzmann (LB) method based on field mediators.The two-phase flow lattice Boltzmann model was first validated and improved by several test cases of a still droplet.The five distinct flow regimes of the kerosene-water system,previously identified in the experiments from Zhao et al.,were reproduced.The quantitative and qualitative agreement between the simulations and the experimental dat...     

Numerical simulation of two-phase partition chromatography in microchannels for moderated log P measurements

A finite element simulation has been used in order to study the partition chromatography process of one species between an aqueous mobile phase and an organic stationary phase located at the bottom of a rectangular microchannel. The transient model incorporates convection--diffusion of the species in the water phase coupled to the diffusion in the stationary organic phase by the way of the partition kinetics at the interface. The time evolution of the injected species concentration is analyzed versus the velocity of the mobile phase, the detecting position and the thickness of the stationary phase. The comparison of simulation results with both experimental data and analytical model confirm its validity. These simulations show that thin channels can be used to measure log P of molecules from their retention time. Finally, we have shown how the sample velocity can be optimized for a given geometry of the channel and diffusion coefficient of the species.     

Design, fabrication and characterization of monolithic embedded parylene microchannels in silicon substrate

This paper presents a novel channel fabrication technology of bulk-micromachined monolithic embedded polymer channels in silicon substrate. The fabrication process favorably obviates the need for sacrificial materials in surface-micromachined channels and wafer-bonding in conventional bulk-micromachined channels. Single-layer-deposited parylene C (poly-para-xylylene C) is selected as a structural material in the microfabricated channels/columns to conduct life science research. High pressure capacity can be obtained in these channels by the assistance of silicon substrate support to meet the needs of high-pressure loading conditions in microfluidic applications. The fabrication technology is completely compatible with further lithographic CMOS/MEMS processes, which enables the fabricated embedded structures to be totally integrated with on-chip micro/nano-sensors/actuators/structures for miniaturized lab-on-a-chip systems. An exemplary process was described to show the feasibility of combining bulk micromachining and surface micromachining techniques in process integration. Embedded channels in versatile cross-section profile designs have been fabricated and characterized to demonstrate their capabilities for various applications. A quasi-hemi-circular-shaped embedded parylene channel has been fabricated and verified to withstand inner pressure loadings higher than 1000 psi without failure for micro-high performance liquid chromatography (microHPLC) analysis. Fabrication of a high-aspect-ratio (internal channel height/internal channel width, greater than 20) quasi-rectangular-shaped embedded parylene channel has also been presented and characterized. Its implementation in a single-mask spiral parylene column longer than 1.1 m in a 3.3 mm x 3.3 mm square size on a chip has been demonstrated for prospective micro-gas chromatography (microGC) and high-density, high-efficiency separations. This proposed monolithic embedded channel technology can be extensively implemented to fabricate microchannels/columns in high-pressure microfluidics and high-performance/high-throughput chip-based micro total analysis systems (microTAS).     

Computer simulation and theory of the diffusion- and flow-induced concentration dispersion in microfluidic devices and HPLC systems based on rectangular microchannels

The dynamics of formation of solute peaks in microfluidic systems are investigated by computer simulation. A finite-element numerical procedure is applied to analyze the diffusion- and flow-controlled concentration dispersion in a 40 microm-high rectangular flow-through channel. Two-dimensional concentration profiles are shown for channels with cross sections of large aspect ratio. The final shapes of the peaks are formed during a very short time period, ranging from a few milliseconds to about 1s for low and high flow velocities, respectively. The observed standard half-width sigma of the peaks is found to strictly follow a linear function of t(1/2) over the whole time range. The extrapolated long-term peak characteristics are in perfect agreement with theoretical predictions. For comparison, theoretical results on the concentration dispersion for solute peaks in open-channel liquid-chromatography (HPLC) are re-examined and applied.     

Electrokinetically-driven flow mixing in microchannels with wavy surface

This paper investigates the mixing characteristics of electrokinetically-driven flow in microchannels with different wavy surface configurations. Numerical simulations are performed to analyze the influence of the wave amplitude and the length of the wavy section on the mixing efficiency within the microchannel. Typically, straight channels have a poor mixing performance because the fluid flow is restricted to the low Reynolds number regime, and hence mixing takes place primarily as a result of diffusion effects. However, the wavy surfaces employed in the current microchannels increase the interfacial contact area between the two species in the microchannel and therefore improve the mixing efficiency. The mixing performance is further enhanced by the application of a heterogeneous charge pattern on the wavy surfaces. The numerical results show that the heterogeneous charge pattern generates flow circulations near the microchannel walls. These circulations are shown to provide an effective enhancement in the mixing performance. Overall, the present results show that the mixing performance is improved by increasing the magnitude of the heterogeneous surface zeta potential upon the wavy surface or by increasing the wave amplitude or the length of the wavy section in the microchannel.     

Numerical investigation of the solute dispersion in finite-length microchannels with the interphase transport

This paper utilizes a combined approach of the convection-diffusion theory and the moment analysis to conduct a comprehensive investigation of the solute dispersion under the influence of the interphase transport in finitely long inner coated microchannels. The present work has threefold novel contributions: (1) The 2D solute concentration contours in the stationary phase are calculated for the first time to facilitate the understanding the role of the interphase transport in the solute dispersion in the mobile phase. (2) The skewness of the elution curves is investigated to guide the control of solute band shape at the channel outlet. (3) The 2D diffusion-convection theory and zero-dimensional (0D) moment analysis complement each other to present a characterization of the solute dispersion behaviors more comprehensive than that by either of the two methods alone. Parametric studies are performed to clarify the effects of four major parameters related to the interphase transport (i.e., stationary phase Péclet number, interphase transport rate, partition coefficient, and stationary phase thickness) on the solute dispersion characteristics. The results from this study provide a straightforward understanding of the effects of interphase transport on the solute dispersion in finitely long microchannels and are of potential relevance to the design and operation of the microfluidics-based analytical devices.     

Engineering and analysis of surface interactions in a microfluidic herringbone micromixer

We developed a computational model and theoretical framework to investigate the geometrical optimization of particle-surface interactions in a herringbone micromixer. The enhancement of biomolecule- and particle-surface interactions in microfluidic devices through mixing and streamline disruption holds promise for a variety of applications. This analysis provides guidelines for optimizing the frequency and specific location of surface interactions based on the flow pattern and relative hydraulic resistance between a groove and the effective channel. The channel bottom, i.e., channel surface between grooves, was identified as the dominant location for surface contact. In addition, geometries that decrease the groove-to-channel hydraulic resistance improve contact with the channel top. Thus, herringbone mixers appear useful for a variety of surface-interaction applications, yet they have largely not been employed in an optimized fashion.     

Thermodynamics of twisted DNA with solvent interaction

The imaginary time path integral formalism is applied to a nonlinear Hamiltonian for a short fragment of heterogeneous DNA with a stabilizing solvent interaction term. Torsional effects are modeled by a twist angle between neighboring base pairs stacked along the molecule backbone. The base pair displacements are described by an ensemble of temperature dependent paths thus incorporating those fluctuational effects which shape the multisteps thermal denaturation. By summing over ~10(7)-10(8) base pair paths, a large number of double helix configurations is taken into account consistently with the physical requirements of the model potential. The partition function is computed as a function of the twist. It is found that the equilibrium twist angle, peculiar of B-DNA at room temperature, yields the stablest helicoidal geometry against thermal disruption of the base pair hydrogen bonds. This result is corroborated by the computation of thermodynamical properties such as fractions of open base pairs and specific heat.     

Pressure generation at the junction of two microchannels with different depths

In this study, we report the design of a microchip‐based hydraulic pump that comprises three glass conduits arranged in a T‐geometry, one of which has a 2 mm long segment shallower (0.5–3 μm in depth) than the remaining 15 μm deep microfluidic network. Upon application of an electric field across this microchannel junction, a mismatch in EOF rate is introduced due to a differential in the fluid conductivity across the deep and shallow segments. Using the reported micropump, pressure‐driven velocities up to 3.2 mm/s have been generated in a 15 μm deep separation channel for an applied voltage of 1.75 kV allowing us to operate under separation conditions that yield the minimum plate height. Moreover, we have shown that this flow velocity can be maximized by optimizing the depth in the shallow region of the T‐geometry. Interestingly however, a simple theory accounting for fluid conductivity differences across microchannels of different depths significantly underestimates the pressure‐driven velocities observed in our experiments. The Taylor dispersion coefficient in our system on the other hand compares well with the theoretical predictions reported in the literature. Finally, the functionality of our device has been demonstrated by implementing a reverse‐phase chromatographic separation that was driven by the pressure‐driven flow generated on‐chip.     

Control and applications of immiscible liquids in microchannels

Photolithography was used in combination with photocleavable self-assembled monolayers to pattern surface free energies inside microchannels enabling the control of the boundary between immiscible liquids. While aqueous solutions are confined to the hydrophilic pathways by surface forces alone, organic liquids are confined to the hydrophobic region only if the aqueous liquid first occupies the hydrophilic region. In this way, stable liquid boundaries between immiscible liquids are possible as long as the pressures are maintained below critical values. The maximum pressures are determined by the interfacial tension of the aqueous solution and organic liquid, channel depth, and advancing contact angle (theta;(a)). Experimental results on maximum pressures are in good agreement with the analytical values. The ability to confine and position the boundary between immiscible liquids inside microchannels leads to a broad range of applications in microfluidic systems, which is exemplified by fabrication of a semipermeable membrane in a surface-patterned channel via interfacial polymerization.     

Particle electrophoresis and dielectrophoresis in curved microchannels

Studies of particle electrophoresis have so far been limited to primarily theoretical or numerical analyses in straight microchannels. Very little work has been done on particle electrophoretic motions in real microchannels that may have one or multiple turns for reducing the devices size or achieving other functions. This article presents an experimental and numerical study of particle electrophoresis in curved microchannels. Polystyrene microparticles are found to migrate across streamlines and flow out of a spiral microchannel in a focused stream near the outer wall. This transverse focusing effect arises from the dielectrophoretic particle motion induced by the nonuniform electric field intrinsic to curved channels. The experimental observations agree quantitatively with the numerical predictions.     

Electroosmotic flow in microchannels with arbitrary geometry and arbitrary distribution of wall charge

General solutions are developed for direct current (DC) and alternating current (AC) electroosmotic flows in microfluidic channels with arbitrary cross-sectional geometry and arbitrary distribution of wall charge (zeta potential). The applied AC electric field can also be of arbitrary waveform. By proposing a nondimensional time scale varpi defined as the ratio of the diffusion time of momentum across the electric double-layer thickness to the period of the applied electric field, we demonstrate analytically that the Helmholtz-Smoluchowski electroosmotic velocity is an appropriate slip condition for AC electroosmotic flows in typical microfluidic applications. With this slip condition approach, electroosmotic flows in rectangular and asymmetric trapezoidal microchannels with nonuniform wall charge, as examples, are investigated. The unknown constants in the proposed general solutions are numerically determined with a least-squares method through matching the boundary conditions. We find that the wall charge affects significantly the electroosmotic flow while the channel geometry does not. Moreover, the flow feature is characterized by another nondimensional time scale Omega defined as the ratio of the diffusion time of momentum across the channel hydraulic radius to the period of the applied electric field. The onset of phase shift between AC electroosmotic velocity and applied electric field is also examined analytically.     

Modeling the interactions between compliant microcapsules and pillars in microchannels

Using a computational model, we investigate the motion of microcapsules inside a microchannel that encompasses a narrow constriction. The microcapsules are composed of a compliant, elastic shell and an encapsulated fluid; these fluid-filled shells model synthetic polymeric microcapsules or biological cells (e.g., leukocytes). Driven by an imposed flow, the capsules are propelled along the microchannel and through the constricted region, which is formed by two pillars that lie in registry, extending from the top and bottom walls of the channels. The tops of these pillars (facing into the microchannel) are modified to exhibit either a neutral or an attractive interaction with the microcapsules. The pillars (and constriction) model topological features that can be introduced into microfluidic devices or the physical and chemical heterogeneities that are inherently present in biological vessels. To simulate the behavior of this complex system, we employ a hybrid method that integrates the lattice Boltzmann model (LBM) for fluid dynamics and the lattice spring model (LSM) for the micromechanics of elastic solids. Through this LBM/LSM technique, we probe how the capsule's stiffness and interaction with the pillars affect its passage through the chambers. The results yield guidelines for regulating the movement of microcarriers in microfluidic systems and provide insight into the flow properties of biological cells in capillaries.     

Nanochannels in SU-8 with floor and ceiling metal electrodes and integrated microchannels

Sacrificially etched 2-D nanofluidic channels and nanospaces with integrated floor and ceiling electrodes and arbitrary channel geometries have been demonstrated with channel heights from 20 nm to 400 nm, widths from 800 nm to 40 microm, and lengths up to 3 mm, using SU-8 as the channel structural material.