Geometrical optimization of helical flow in grooved micromixers

Owing to the enhancement of surface effects at the micro-scale, patterned grooves on a micro-channel floor remain a powerful method to induce helical flows within a pressure driven system. Although there have been a number of numerical studies on geometrical effects concerning fluid mixing within the staggered herringbone mixer, all have focused mainly on the groove angle and depth, two factors that contribute greatly to the magnitude of helical flow. Here we present a new geometrical factor that significantly affects the generation of helical flow over patterned grooves. By varying the ratio of the length of the grooves to the neighboring ridges, helical flow can be optimized for a given groove depth and channel aspect ratio, with up to 50% increases in transverse flow possible. A thorough numerical study of over 700 cases details the magnitude of helical flow over unsymmetrical patterned grooves in a slanted groove micro-mixer, where the optimized parameters for the slanted groove mixer can be translated to the staggered herringbone mixer. The optimized groove geometries are shown to have a large dependence on the channel aspect ratio, the groove depth ratio, and the ridge length.     

A practical guide to the staggered herringbone mixer

An analytical model of mixing in the staggered herringbone mixer (SHM) was derived to estimate mixing parameters and provide practical expressions to guide mixer design and operation for a wide range of possible solutes and flow conditions. Mixing in microfluidic systems has historically been characterized by the mixing of a specific solute system or by the redistribution of flow streams; this approach does not give any insight into the ideal operational parameters of the mixer with an arbitrary real system. For Stokes-flow mixers, mixing can be computed from a relationship between solute diffusivity, flow rate, and mixer length. Confocal microscopy and computational fluid dynamics (CFD) modeling were used to directly determine the extent of mixing for several solutes in the staggered herringbone mixer over a range of Reynolds numbers (Re) and Péclet numbers (Pe); the results were used to develop and evaluate an analytical model of its behavior. Mixing was found to be a function of only Pe and downstream position in the mixer. Required mixer length was proportional to log(Pe); this analytical model matched well with the confocal data and CFD model for Pe<5 x 10(4), at which point the experiments reached the limit of resolution. For particular solutes, required length and mixing time depend upon Re and diffusivity. This analytical model is applicable to other solute systems, and possibly to other embodiments of the mixer, to enable optimal design, operation, and estimation of performance.     

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.     

Effects of geometrical characteristics of surface roughness on droplet wetting

Surface roughness is known to alter the wettability on a solid substrate. In general, either Wenzel or Cassie-Baxter theory is adopted to describe the apparent contact angle. Following the minimum free energy pathway associated with the imbibition process, we have derived a generalized expression for the apparent contact angle on a textured surface and the liquid-gas contact area within the groove that plays a key role. Depending on the geometrical characteristics of the grooves, the surface wetting falls into three regimes: (i) single stable state which is either Wenzel (completely wetted roughness) or Cassie-Baxter (completely nonwetted roughness) state, (ii) two stable states (Wenzel and Cassie-Baxter) separated by an energy barrier, and (iii) single stable state with partially wetted roughness. The sufficient condition for each regime is derived and several groove geometries are given to show the free energy path. Alteration in the geometric parameters may lead to the wetting crossover. We also show that the Cassie-Baxter can occur at a hydrophilic surface for particular pore shapes.     

Toolbox for the design of optimized microfluidic components

A computational "toolbox" for the a priori design of optimized microfluidic components is presented. These components consist of a microchannel under low-Reynolds number, pressure-driven flow, with an arrangement of grooves cut into the top and bottom to generate a tailored cross-channel flow. An advection map for each feature (i.e., groove of a particular shape and orientation) predicts the lateral transport of fluid within the channel due to that feature. We show that applying these maps in sequence generates an excellent representation of the outflow distribution for complex designs that combine these basic features. The effect of the complex three-dimensional flow field can therefore be predicted without solving the governing flow equations through the composite geometry, and the resulting distribution of fluids in the channel is used to evaluate how well a component performs a specified task. The generation and use of advection maps is described, and the toolbox is applied to determine optimal combinations of features for specified mixer sizes and mixing metrics.     

Rapid droplet mixers for digital microfluidic systems

The mixing of analytes and reagents for a biological or chemical lab-on-a-chip is an important, yet difficult, microfluidic operation. As volumes approach the sub-nanoliter regime, the mixing of liquids is hindered by laminar flow conditions. An electrowetting-based linear-array droplet mixer has previously been reported. However, fixed geometric parameters and the presence of flow reversibility have prevented even faster droplet mixing times. In this paper, we study the effects of varying droplet aspect ratios (height:diameter) on linear-array droplet mixers, and propose mixing strategies applicable for both high and low aspect ratio systems. An optimal aspect ratio for four electrode linear-array mixing was found to be 0.4, with a mixing time of 4.6 seconds. Mixing times were further reduced at this ratio to less than three seconds using a two-dimensional array mixer, which eliminates the effects of flow reversibility. For lower aspect ratio (相似文献   

Mixing of non-Newtonian fluids in wavy serpentine microchannel using electrokinetically driven flow

A numerical investigation is performed into the mixing performance of electrokinetically driven non-Newtonian fluids in a wavy serpentine microchannel. The flow behavior of the non-Newtonian fluids is described using a power-law model. The simulations examine the effects of the flow behavior index, the wave amplitude, the wavy-wall section length, and the applied electric field strength on the mixing performance. The results show that the volumetric flow rate of shear-thinning fluids is higher than that of shear-thickening fluids, and therefore results in a poorer mixing performance. It is shown that for both types of fluid, the mixing performance can be enhanced by increasing the wave amplitude, extending the length of the wavy-wall section, and reducing the strength of the electric field. Thus, although the mixing efficiency of shear-thinning fluids is lower than that of shear-thickening fluids, the mixing performance can be improved through an appropriate specification of the flow and geometry parameters. For example, given a shear-thinning fluid with a flow behavior index of 0.8, a mixing efficiency of 87% can be obtained by specifying the wave amplitude as 0.7, the wavy-wall section length as five times the characteristic length, the nondimensional Debye-Huckel parameter as 100, and the applied electric field strength as 43.5 V/cm.     

Microcirculation within grooved substrates regulates cell positioning and cell docking inside microfluidic channels

Immobilization of cells inside microfluidic devices is a promising approach for enabling studies related to drug screening and cell biology. Despite extensive studies in using grooved substrates for immobilizing cells inside channels, a systematic study of the effects of various parameters that influence cell docking and retention within grooved substrates has not been performed. We demonstrate using computational simulations that the fluid dynamic environment within microgrooves significantly varies with groove width, generating microcirculation areas in smaller microgrooves. Wall shear stress simulation predicted that shear stresses were in the opposite direction in smaller grooves (25 and 50 microm wide) in comparison to those in wider grooves (75 and 100 microm wide). To validate the simulations, cells were seeded within microfluidic devices, where microgrooves of different widths were aligned perpendicularly to the direction of the flow. Experimental results showed that, as predicted, the inversion of the local direction of shear stress within the smaller grooves resulted in alignment of cells on two opposite sides of the grooves under the same flow conditions. Also, the amplitude of shear stress within microgrooved channels significantly influenced cell retainment in the channels. Therefore, our studies suggest that microscale shear stresses greatly influence cellular docking, immobilization, and retention in fluidic systems and should be considered for the design of cell-based microdevices.     

Scale-up of the vortex wave microfilter using the power ratio

The vortex wave technique enhances microfiltration performance by combining the mixing characteristics of oscillatory flow and flow deflectors. The ability to reproduce these highly convective mixing patterns on a scale applicable to pilot plant studies has been investigated in this paper. The scale-up of a single channel membrane unit (18200 mm²) to a double channel element (106080 mm²) has been investigated in terms of geometric, kinematic and dynamic similarity. A non-dimensional approach has been adopted, wherever possible, to characterise the performance of both microfilters. The kinetics of each system are dominated by the Reynolds number, and the dynamics have been classified in terms of a new non-dimensional number — the power ratio. The power ratio has been employed to predict accurately the performance of the scaled-up filter, and also to provide an indication of transitional flow patterns. The inhibition of a gel-layer within the vortex wave modules has been demonstrated by the non-linear relationship between flux and transmembrane pressure. Experimental results and a theory which includes the effects of osmotic-pressure show that the flux is directly proportional to the cube root of the transmembrane pressure. Non-dimensional analysis of the experimental results indicates a similar trend, whereby the flux is directly proportional to the cube root of power dissipation under two-dimensional, laminar flow conditions.     

Exploring DNA groove water dynamics through hydrogen bond lifetime and orientational relaxation

Dynamics of water molecules in the grooves of DNA are of great interest both for practical (functionality of DNA) and fundamental (as examples of confined systems) interest. Here the authors employ atomistic molecular dynamics simulations to understand varying water dynamics at the minor and the major grooves of a 38 base-pair long DNA duplex in water. In order to understand and quantify the diversity in the nature of hydrogen bond due to many hydrogen bond donors and acceptors present in the four bases, they have undertaken study of hydrogen bond lifetime (HBLT) correlation functions of all the specific hydrogen bonds between the base atoms and water molecules. They find that the HBLT correlation functions are in general multiexponential, with the average lifetime depending significantly on the specificity and may thus be biologically relevant. The average hydrogen bond lifetime is longer in the minor groove than that in the major groove by almost a factor of 2. Analysis further shows that water hydrogen bonds with phosphate oxygen have substantially shorter lifetimes than those with the groove atoms. They also compute two different orientational time correlation functions (OTCFs) of the water molecules present at the major and the minor grooves and attempt to correlate OTCF with HBLT correlation function. The OTCFs in the minor groove exhibit three time scales, with the time constant of the slowest component one to two orders of magnitude longer than what is observed for bulk water. A slow component is also present for the major groove water but with shorter time constant. Interestingly, correlation between reformations allowed HBLT correlation function [C(HB)(t)] and the OTCF markedly deviates from each other in the grooves, indicating enhanced rigidity of water molecules in the grooves.     

Electrostatic and dispersion interactions during protein adsorption on topographic nanostructures

Recently, biomaterials research has focused on developing functional implant surfaces with well-defined topographic nanostructures in order to influence protein adsorption and cellular behavior. To enhance our understanding of how proteins interact with such surfaces, we analyze the adsorption of lysozyme on an oppositely charged nanostructure using a computer simulation. We present an algorithm that combines simulated Brownian dynamics with numerical field calculation methods to predict the preferred adsorption sites for arbitrarily shaped substrates. Either proteins can be immobilized at their initial adsorption sites or surface diffusion can be considered. Interactions are analyzed on the basis of Derjaguin-Landau-Verway-Overbeek (DLVO) theory, including electrostatic and London dispersion forces, and numerical solutions are derived using the Poisson-Boltzmann and Hamaker equations. Our calculations show that for a grooved nanostructure (i.e., groove and plateau width 8 nm, height 4 nm), proteins first contact the substrate primarily near convex edges because of better geometric accessibility and increased electric field strengths. Subsequently, molecules migrate by surface diffusion into grooves and concave corners, where short-range dispersion interactions are maximized. In equilibrium, this mechanism leads to an increased surface protein concentration in the grooves, demonstrating that the total amount of protein per surface area can be increased if substrates have concave nanostructures.     

Changing DNA grooves--a 1,4,5,8-naphthalene tetracarboxylic diimide bis-intercalator with the linker (beta-Ala)(3)-Lys in the minor groove

We have been investigating a modular, threading DNA polyintercalator design based upon the 1,4,5,8-naphthalene tetracarboxylic diimide (NDI) intercalating unit. Previously, we have reported the NMR analysis of a bis-intercalator-DNA complex in which the peptide linker between NDI units was found to occupy the DNA major groove (Guelev, Lee, Sorey, Hoffman, Iverson, Chem. Biol. 2001, 8, 415-425). Here we describe the NMR analysis of a complex between a related bis-intercalator known to display altered DNA sequence specificity. In this case, the linker resides in the DNA minor groove. We have thus shown that within this set of sequence specific bis-intercalators, both DNA grooves can be accessed, setting the stage for longer threading polyintercalators designed to have linkers occupying both grooves in an alternating fashion.     

Passive and Active Mixing in Microfluidic Devices

In microfluidics the Reynolds number is small, preventing turbulence as a tool for mixing, while diffusion is that slow that time does not yield an alternative. Mixing in microfluidics therefore must rely on chaotic advection, as well-known from polymer technology practice where on macroscale the high viscosity makes the Reynolds numbers low and diffusion slow. The mapping method is used to analyze and optimize mixing also in microfluidic devices. We investigate passive mixers like the staggered herringbone micromixer (SHM), the barrier embedded micromixer (BEM) and a three-dimensional serpentine channel (3D-SC). Active mixing is obtained via incorporating particles that introduce a hyperbolic flow in e.g. two dimensional serpentine channels. Magnetic beads chains-up in a flow after switching on a magnetic field. Rotating the field creates a physical rotor moving the flow field. The Mason number represents the ratio of viscous forces to the magnetic field strength and its value determines the fate of the rotor: a single, an alternating single and double, or a multiple part chain-rotor results. The type of rotor determines the mixing quality with best results in the alternating case where crossing streamlines introduce chaotic advection. Finally, an active mixing device is proposed that mimics the cilia in nature. The transverse flow induced by their motion indeed enhances mixing at the microscale.     

Anisotropic wetting characteristics on submicrometer-scale periodic grooved surface

Submicrometer-scale periodic structures consisting of parallel grooves were prepared on azobenzene-containing multiarm star polymer films by laser interference. The wetting characteristics on the patterned surfaces were studied by contact angle measurements. Macroscopic distortion of water drops was found on such small-scale surface structures, and the contact angles measured from the direction parallel to the grooves were larger than those measured from the perpendicular direction. A thermodynamic model was developed to calculate the change in the surface free energy as a function of the instantaneous contact angle when the three-phase contact line (TPCL) moves along the two orthogonal directions. It was found that the fluctuations, i.e., energy barriers, on the energy versus contact angle curves are crucial to the analysis of wetting anisotropy and contact angle hysteresis. The calculated advancing and receding contact angles from the energy versus contact angle curves were in good agreement with those measured experimentally. Furthermore, with the groove depth increasing, both the degree of wetting anisotropy and the contact angle hysteresis perpendicular to the grooves increased as a result of the increase in the energy barrier. The theoretical critical value of the groove depth, above which the anisotropic wetting appears, was determined to be 16 nm for the grooved surface with a wavelength of 396 nm. On the other hand, the effect of the groove wavelength on the contact angle hysteresis perpendicular to the grooves was also interpreted on the basis of the thermodynamic model. That is, with the wavelength decreasing, the contact angle hysteresis increased due to the increase in the number of energy barriers. These results may provide theoretical evidence for the design and application of anisotropic wetting surface.     

Entropy of water in the hydration layer of major and minor grooves of DNA

Transport properties (translational and rotational) of water in the two grooves of the B-DNA duplex are known to be different from those in the bulk. Here, we use a recently developed theoretical scheme to compute the entropies of water molecules in both of the grooves of DNA and compare them with that in the bulk. The scheme requires as input both translational and rotational velocity autocorrelation function (C(V)(t) and C(omega)(t), respectively) data. These velocity autocorrelation functions were computed from an atomistic MD simulation of a B-DNA duplex (36 base pairs long) in explicit water (TIP3P). The average values of the entropy of water at 300 K in both of the grooves of DNA (the TS value in the major groove is 6.71 kcal/mol and that in the minor groove is 6.41 kcal/mol) are found to be significantly lower than that in bulk water (the TS value is 7.27 kcal/mol). Thus, the entropic contribution to the free energy change (TDeltaS) of transferring a minor groove water molecule to the bulk is 0.86 kcal/mol and of transferring a major groove water to the bulk is 0.56 kcal/mol at 300 K, which is to be compared with 1.44 kcal/mol for melting of ice at 273 K. We also calculate the energy of interaction of each water molecule with the rest of the atoms in the system and hence calculate the chemical potential (Helmholtz free energy per water molecule, A = E - TS) in the different domains. The identical free energy value of water molecules in the different domains proves the robustness of the scheme. We propose that the configurational entropy of water in the grooves can be used as a measure of the mobility (or microviscosity) of water molecules in a given domain.     

Yeast cell microfiltration: Flux enhancement in baffled and pulsatile flow systems

The filtration performance of yeast cell harvesting was greatly improved by using an oscillatory flow mixing technique in both tubular and flat sheet membrane systems; flux increases of several fold were achieved. The effects of various operating parameters such as frequency, amplitude and ratio of net forward flow rate to oscillatory flow rate were investigated. A simple alternative to the dual head pump, which involved the use of solenoid valves for producing flow reversal, was developed. The flux was significantly increased using this new operating mode. As the only increase of power consumption was from inertia loss due to change of flow direction the additional power consumption is minimal.     

Numerical analysis of a rapid magnetic microfluidic mixer

This paper presents a detailed numerical investigation of the novel active microfluidic mixer proposed by Wen et al. (Electrophoresis 2009, 30, 4179-4186). This mixer uses an electromagnet driven by DC or AC power to induce transient interactive flows between a water-based ferrofluid and DI water. Experimental results clearly demonstrate the mixing mechanism. In the presence of the electromagnet's magnetic field, the magnetic nanoparticles create a body force vector that acts on the mixed fluid. Numerical simulations show that this magnetic body force causes the ferrofluid to expand significantly and uniformly toward miscible water. The magnetic force also produces many extremely fine finger structures along the direction of local magnetic field lines at the interface in both upstream and downstream regions of the microchannel when the external steady magnetic strength (DC power actuation) exceeds 30 Oe (critical magnetic Peclet number Pe(m),cr = 2870). This study is the first to analyze these pronounced finger patterns numerically, and the results are in good agreement with the experimental visualization of Wen et al. (Electrophoresis 2009, 30, 4179-4186). The large interfacial area that accompanies these fine finger structures and the dominant diffusion effects occurring around the circumferential regions of fingers significantly enhance the mixing performance. The mixing ratio can be as high as 95% within 2.0 s. at a distance of 3.0 mm from the mixing channel inlet when the applied peak magnetic field supplied by the DC power source exceeds 60 Oe. This study also presents a sample implementation of AC power actuation in a numerical simulation, an experimental benchmark, and a simulation of DC power actuation with the same peak magnetic strength. The simulated flow structures of the AC power actuation agree well with the experimental visualization, and are similar to those produced by DC power. The AC and DC power actuated flow fields exhibited no significant differences. This numerical study suggests approaches to maximize the performance of the proposed rapid magnetic microfluidic mixer, and confirms its exciting potential for use in lab-on-a-chip systems.     

Differential solvation and tautomer stability of a model base pair within the minor and major grooves of DNA

2-(2'-Hydroxyphenyl)benzoxazole (HBO) may be used as a model base pair to study solvation, duplex environment, and tautomerization within the major and minor groves of DNA duplexes. In its ground state, HBO possesses an enol moiety which may be oriented syn or anti relative to the imino nitrogen of the benzoxazole ring. In the absence of external hydrogen-bond donors and acceptors HBO exists as the internally hydrogen-bonded syn-enol, a mimic of the rare base pair tautomer found in DNA, which may be photoinduced to tautomerize and form the keto tautomer, a mimic of the dominant base pair tautomer. Previously, we demonstrated that when incorporated into DNA such that the enol moiety is positioned in the major groove, HBO is not solvated, exists exclusively as the internally hydrogen-bonded syn-enol which is efficiently photoinduced to tautomerize, and the corresponding keto tautomer is preferentially stabilized. In stark contrast, we now show that when HBO is incorporated in DNA such that the enol moiety is positioned in the minor groove, the enol tautomer is preferentially stabilized. Molecular dynamics simulations suggest that this results from the formation of a stable hydrogen-bond between the HBO enol and the O4' atom of an adjacent nucleotide, an H-bond acceptor that is only available in the minor groove. The differential stabilization of the enol and keto tautomers in the major and minor grooves may reflect the functions for which these environments evolved, including duplex replication, stability, and recognition.     

Erosion resistance of bionic functional surfaces inspired from desert scorpions

In this paper, a bionic method is presented to improve the erosion resistance of machine components. Desert scorpion (Androctonus australis) is a typical animal living in sandy deserts, and may face erosive action of blowing sand at a high speed. Based on the idea of bionics and biologic experimental techniques, the mechanisms of the sand erosion resistance of desert scorpion were investigated. Results showed that the desert scorpions used special microtextures such as bumps and grooves to construct the functional surfaces to achieve the erosion resistance. In order to understand the erosion resistance mechanisms of such functional surfaces, the combination of computational and experimental research were carried out in this paper. The Computational Fluid Dynamics (CFD) method was applied to predict the erosion performance of the bionic functional surfaces. The result demonstrated that the microtextured surfaces exhibited better erosion resistance than the smooth surfaces. The further erosion tests indicated that the groove surfaces exhibited better erosion performance at 30° injection angle. In order to determine the effect of the groove dimensions on the erosion resistance, regression analysis of orthogonal multinomials was also performed under a certain erosion condition, and the regression equation between the erosion rate and groove distance, width, and height was established.     

A microfluidic mixer with grooves placed on the top and bottom of the channel

A new microfluidic mixer is presented consisting of a rectangular channel with grooves placed in the top and bottom. This not only increases the driving force behind the lateral flow, but allows for the formation of advection patterns that cannot be created with structures on the bottom alone. Chevrons, pointing in opposite directions on the top and bottom, are used to create a pair of vortices positioned side by side. Stripes running the width of the channel generate a pair of vertically stacked vortices. Computational fluid dynamics (CFD) simulations are used to model the behavior of the systems and provide velocity maps at cross-sections within the mixer. Experiments demonstrate the mixing that results when two segregated species enter the mixer side-by-side and pass through two cycles of the mixer (i.e., two alternating sets of four stripes and four chevrons).