Numerical Study of Electro-Osmotic Fluid Flow and Vortex Formation

Review on Blood Flow Dynamics in Lab-on-a-Chip Systems: An Engineering Perspective

Under different transport mechanisms, blood flow dynamics, heavily linked to the flow shear rate conditions, in "lab-on-a-chip" (LOC) systems are found to result in varying transport phenomena. This Review examines the blood flow patterns in LOC systems through the role of viscoelastic properties such as dynamic blood viscosity and elastic behavior of the red blood cells. The study of blood transport phenomena in LOC systems through key parameters of capillary and electro-osmotic forces is provided through experimental, theoretical, and numerous numerical approaches. The disturbance triggered by electro-osmotic viscoelastic flow is particularly discussed and applied in the enhancement of the mixing and separating capabilities of LOC devices handling blood and other viscoelastic fluids for future research opportunities. Furthermore, the Review identifies the challenges in the numerical modeling of blood flow dynamics under the LOC systems, such as the call for more accurate and simplified blood flow models and the emphasis on numerical studies of viscoelastic fluid flow under the electrokinetic effect. More practical assumptions for zeta potential conditions while studying the electrokinetic phenomena are also highlighted. This Review aims to provide a comprehensive and interdisciplinary perspective on blood flow dynamics in microfluidic systems driven by capillary and electro-osmotic forces.

Particle Size-Dependent Electrophoresis in Polymer Solutions

It has long been known that the electrophoretic velocity of a charged particle is independent of its size under the thin-Debye-layer limit. This so-called Smoluchowski velocity is, however, valid only for Newtonian fluids. A couple of recent theoretical studies predict the rheology-induced particle size dependence of electrophoresis in non-Newtonian fluids. This work presents the first experimental demonstration of such dependence in viscoelastic poly(ethylene oxide) (PEO) solutions. Three different-sized particles are observed to travel at the same electrophoretic velocity in a Newtonian buffer through a rectangular microchannel. In contrast, their measured electrophoretic velocities in the PEO solution exhibit an increasing trend for larger particles, which is consistent with theoretical prediction. This particle size dependence is found to grow with an increasing concentration or length of the PEO polymer. Both trends are attributed to enhanced fluid elasticity, as characterized by the increasing elasticity number.

Ion transport and current rectification in a charged conical nanopore filled with viscoelastic fluids

The ionic current rectification (ICR) is a non-linear current-voltage response upon switching the polarity of the potential across nanopore which is similar to the I–V response in the semiconductor diode. The ICR phenomenon finds several potential applications in micro/nano-fluidics (e.g., Bio-sensors and Lab-on-Chip applications). From a biological application viewpoint, most biological fluids (e.g., blood, saliva, mucus, etc.) exhibit non-Newtonian visco-elastic behavior; their rheological properties differ from Newtonian fluids. Therefore, the resultant flow-field should show an additional dependence on the rheological material properties of viscoelastic fluids such as fluid relaxation time \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\lambda )$$\end{document}(λ) and fluid extensibility \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\varepsilon )$$\end{document}(ε). Despite numerous potential applications, the comprehensive investigation of the viscoelastic behavior of the fluid on ionic concentration profile and ICR phenomena has not been attempted. ICR phenomena occur when the length scale and Debye layer thickness approaches to the same order. Therefore, this work extensively investigates the effect of visco-elasticity on the flow and ionic mass transfer along with the ICR phenomena in a single conical nanopore. The Poisson–Nernst–Planck (P–N–P) model coupled with momentum equations have been solved for a wide range of conditions such as, Deborah number, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\le De \le 100$$\end{document}1≤De≤100, Debye length parameter, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1\le \kappa R_t \le 50$$\end{document}1≤κRt≤50, fluid extensibility parameter, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0.05\le \varepsilon \le 0.25$$\end{document}0.05≤ε≤0.25, applied electric potential, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-40\le V \le 40$$\end{document}-40≤V≤40, and surface charge density \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma = -10$$\end{document}σ=-10 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$-50$$\end{document}-50. Limited results for Newtonian fluid (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$De = 0$$\end{document}De=0, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varepsilon = 0$$\end{document}ε=0) have also been shown in order to demonstrate the effectiveness of non-Newtonian fluid behaviour over the Newtonian fluid behaviour. Four distinct novel characteristics of electro-osmotic flow (EOF) in a conical nanopore have been investigated here, namely (1) detailed structure of flow field and velocity distribution in viscoelastic fluids (2) influence of Deborah number and fluid extensibility parameter on ionic current rectification (ICR) (3) volumetric flow rate calculation as a function of Deborah number and fluid extensibility parameter (4) effect of viscoelastic parameters on concentration distribution of ions in the nanopore. At high applied voltage, both the extensibility parameter and Deborah number facilitate the ICR phenomena. In addition, the ICR phenomena are observed to be more pronounced at low values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa R_t$$\end{document}κRt than the high values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa R_t$$\end{document}κRt. This effect is due to the overlapping of the electric double layer at low values of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa R_t$$\end{document}κRt.

Fluid rheological effects on streaming dielectrophoresis in a post-array microchannel

Recent studies have demonstrated the strong influences of fluid rheological properties on insulator-based dielectrophoresis (iDEP) in single-constriction microchannels. However, it is yet to be understood how iDEP in non-Newtonian fluids depends on the geometry of insulating structures. We report in this work an experimental study of fluid rheological effects on streaming DEP in a post-array microchannel that presents multiple contractions and expansions. The iDEP focusing and trapping of particles in a viscoelastic polyethylene oxide solution are comparable to those in a Newtonian buffer, which is consistent with the observations in a single-constriction microchannel. Similarly, the insignificant iDEP effects in a shear-thinning xanthan gum solution also agree with those in the single-constriction channel except that gel-like structures are observed to only form in the post-array microchannel under large DC electric fields. In contrast, the iDEP effects in both viscoelastic and shear-thinning polyacrylamide solution are significantly weaker than in the single-constriction channel. Moreover, instabilities occur in the electroosmotic flow and appear to be only dependent on the DC electric field. These phenomena may be associated with the dynamics of polymers as they are electrokinetically advected around and through the posts.

A Hierarchical Grid Solver for Simulation of Flows of Complex Fluids

Tree-based grids bring the advantage of using fast Cartesian discretizations, such as finite differences, and the flexibility and accuracy of local mesh refinement. The main challenge is how to adapt the discretization stencil near the interfaces between grid elements of different sizes, which is usually solved by local high-order geometrical interpolations. Most methods usually avoid this by limiting the mesh configuration (usually to graded quadtree/octree grids), reducing the number of cases to be treated locally. In this work, we employ a moving least squares meshless interpolation technique, allowing for more complex mesh configurations, still keeping the overall order of accuracy. This technique was implemented in the HiG-Flow code to simulate Newtonian, generalized Newtonian and viscoelastic fluids flows. Numerical tests and application to viscoelastic fluid flow simulations were performed to illustrate the flexibility and robustness of this new approach.

Insulator-based dielectrophoretic focusing and trapping of particles in non-Newtonian fluids

Insulator-based dielectrophoretic (iDEP) microdevices have been limited to work with Newtonian fluids. We report an experimental study of the fluid rheological effects on iDEP focusing and trapping of polystyrene particles in polyethylene oxide, xanthan gum, and polyacrylamide solutions through a constricted microchannel. Particle focusing and trapping in the mildly viscoelastic polyethylene oxide solution are slightly weaker than in the Newtonian buffer. They are, however, significantly improved in the strongly viscoelastic and shear thinning polyacrylamide solution. These observed particle focusing behaviors exhibit a similar trend with respect to electric field, consistent with a revised theoretical analysis for iDEP focusing in non-Newtonian fluids. No apparent focusing of particles is achieved in the xanthan gum solution, though the iDEP trapping can take place under a much larger electric field than the other fluids. This is attributed to the strong shear thinning-induced influences on both the electroosmotic flow and electrokinetic/dielectrophoretic motions.

Influence of Non-Structural Parameters on Dual Parallel Jet Characteristics of Porous Nozzles

As an important actuator of the dual parallel jet, the porous nozzle has some non-structural parameters (such as inlet pressure, nozzle spacing ratio, etc.) which have a significant influence on energy transport, chemical combustion and pollutant generation. The research on the microfluidic state of the porous nozzle dual parallel jet, however, remains insufficient because of its microjet pattern and complex intersection process. In this paper, the authors used numerical simulation and an experimental method to clarify the influence of porous nozzles’ non-structural parameters on dual parallel jet characteristics. The results show that the inlet pressure only changes the pressure peak value on the parallel jet axis; the starting point (SP) and peak point (PP) on the parallel jet axis, which are located at X_sp = 22 mm and X_pp = 75 mm, respectively, are not changed; and with the increase in the nozzle spacing ratio, the merging points (MPs) on the parallel jet axis are X_mp = 25 mm, 32 mm and 59 mm, respectively. The merging point and the combined point move to a farther distance and the inner deflection angle of the jet is weakened.