Combined Effect of Surface Roughness and Heterogeneity of Wall Potential on Electroosmosis in Microfluidic/Nanofuidic Channels

2010 ◽  
Vol 132 (4) ◽  
Author(s):  
S. Bhattacharyya ◽  
A. K. Nayak

The motivation of the present study is to generate vortical flow by introducing channel wall roughness in the form of a wall mounted block that has a step-jump in ζ-potential on the upper face. The characteristics for the electrokinetic flow are obtained by numerically solving the Poisson equation, the Nernst–Planck equation, and the Navier–Stokes equations, simultaneously. A numerical method based on the pressure correction iterative algorithm (SIMPLE) is adopted to compute the flow field and mole fraction of the ions. The potential patch induces a strong recirculation vortex, which in turn generates a strong pressure gradient. The strength of the vortex, which appears adjacent to the potential patch, increases almost linearly with the increase in ζ-potential. The streamlines follow a tortuous path near the wall roughness. The average axial flow rate over the block is enhanced significantly. We found that the ionic distribution follow the equilibrium Boltzmann distribution away from the wall roughness. The solutions based on the Poisson–Boltzmann distribution and the Nernst–Planck model are different when the inertial effect is significant. The combined effects due to geometrical modulation of the channel wall and heterogeneity in ζ-potential is found to produce a stronger vortex, and hence a stronger mixing, compared with either of these. Increase in ζ-potential increases both the transport rate and mixing efficiency. A novelty of the present configuration is that the vortex forms above the obstacle even when the patch potential is negative.

Micromachines ◽  
2021 ◽  
Vol 12 (4) ◽  
pp. 431
Author(s):  
Lanju Mei ◽  
Defu Cui ◽  
Jiayue Shen ◽  
Diganta Dutta ◽  
Willie Brown ◽  
...  

This paper investigates the electroosmotic micromixing of non-Newtonian fluid in a microchannel with wall-mounted obstacles and surface potential heterogeneity on the obstacle surface. In the numerical simulation, the full model consisting of the Navier–Stokes equations and the Poisson–Nernst–Plank equations are solved for the electroosmotic fluid field, ion transport, and electric field, and the power law model is used to characterize the rheological behavior of the aqueous solution. The mixing performance is investigated under different parameters, such as electric double layer thickness, flow behavior index, obstacle surface zeta potential, obstacle dimension. Due to the zeta potential heterogeneity at the obstacle surface, vortical flow is formed near the obstacle surface, which can significantly improve the mixing efficiency. The results show that, the mixing efficiency can be improved by increasing the obstacle surface zeta potential, the flow behavior index, the obstacle height, the EDL thickness.


2000 ◽  
Vol 411 ◽  
pp. 1-38 ◽  
Author(s):  
C. ROSS ETHIER ◽  
SUJATA PRAKASH ◽  
DAVID A. STEINMAN ◽  
RICHARD L. LEASK ◽  
GREGORY G. COUCH ◽  
...  

Numerical and experimental techniques were used to study the physics of flow separation for steady internal flow in a 45° junction geometry, such as that observed between two pipes or between the downstream end of a bypass graft and an artery. The three-dimensional Navier–Stokes equations were solved using a validated finite element code, and complementary experiments were performed using the photochromic dye tracer technique. Inlet Reynolds numbers in the range 250 to 1650 were considered. An adaptive mesh refinement approach was adopted to ensure grid-independent solutions. Good agreement was observed between the numerical results and the experimentally measured velocity fields; however, the wall shear stress agreement was less satisfactory. Just distal to the ‘toe’ of the junction, axial flow separation was observed for all Reynolds numbers greater than 250. Further downstream (approximately 1.3 diameters from the toe), the axial flow again separated for Re [ges ] 450. The location and structure of axial flow separation in this geometry is controlled by secondary flows, which at sufficiently high Re create free stagnation points on the model symmetry plane. In fact, separation in this flow is best explained by a secondary flow boundary layer collision model, analogous to that proposed for flow in the entry region of a curved tube. Novel features of this flow include axial flow separation at modest Re (as compared to flow in a curved tube, where separation occurs only at much higher Re), and the existence and interaction of two distinct three-dimensional separation zones.


Author(s):  
Rubén Bruno Díaz ◽  
Jesuino Takachi Tomita ◽  
Cleverson Bringhenti ◽  
Francisco Carlos Elizio de Paula ◽  
Luiz Henrique Lindquist Whitacker

Abstract Numerical simulations were carried out with the purpose of investigating the effect of applying circumferential grooves at axial compressor casing passive wall treatment to enhance the stall margin and change the tip leakage flow. The tip leakage flow is pointed out as one of the main contributors to stall inception in axial compressors. Hence, it is of major importance to treat appropriately the flow in this region. Circumferential grooves have shown a good performance in enhancing the stall margin in previous researches by changing the flow path in the tip clearance region. In this work, a passive wall treatment with four circumferential grooves was applied in the transonic axial compressor NASA Rotor 37. Its effect on the axial compressor performance and the flow in the tip clearance region was analyzed and set against the results attained for the smooth wall case. A 2.63% increase in the operational range of the axial compressor running at 100%N, was achieved, when compared with the original smooth wall casing configuration. The grooves installed at compressor casing, causes an increase in the flow entropy generation due to the high viscous effects in this gap region, between the rotor tip surface and casing with grooves. These viscous effects cause a drop in the turbomachine efficiency. For the grooves configurations used in this work, an efficiency drop of 0.7% was observed, compared with the original smooth wall. All the simulations were performed based on 3D turbulent flow calculations using Reynolds Averaged Navier-Stokes equations, and the flow eddy viscosity was determined using the two-equation SST turbulence model. The details of the grooves geometrical dimensions and its implementation are described in the paper.


Author(s):  
Katsuhisa Fujita ◽  
Makoto Katou

The unstable phenomena of thin cylindrical shells subjected to annular axial flow are investigated. In this paper, the analytical model is composed of an elastic axisymmetric shell and a rigid one which are arranged co-axially. Considering the fluid structure interaction between shells and fluid flowing through an annular narrow passage, the coupled equation of motion is derived using Flu¨gge’s shell theory and Navier-Stokes equations. The unstable phenomena of thin cylindrical shells are clarified by using the root locus based on the complex eigenvalue analysis. The numerical parameter studies on the shells with a freely supported end and a rigid one, and with both simply supported ends, are performed taking the dimensins of shells, the characteristics of flowing fluid so on as parameters. The influence of these parameters on the threshold of instability of the coupled vibration between thin cylindrical shells and annular axial flowing fluid are investigated and discussed.


Author(s):  
R. Capanna ◽  
G. Ricciardi ◽  
C. Eloy ◽  
E. Sarrouy

Efficient modelling and accurate knowledge of the mechanical behaviour of the reactor core are needed to estimate the effects of seismic excitation on a nuclear power plant. The fuel assemblies (in the reactor core) are subjected to an axial water flow which modifies their dynamical behaviour. Several fluid-structure models simulating the response of the core to a seismic load has been developed in recent years; most of them require high computational costs. The work which is presented here is a first step in order to simplify the fluid forces modelling, and thus to be able to catch the main features of the mechanical behaviour of reactor core with low computational costs. The main assumption made in this work is to consider the fluid flow as an inviscid potential flow. Thus, the flow can be described only using one scalar function (velocity potential) instead of a vector field and strongly simplifies the fluid mechanics equations, avoiding the necessity to solve Navier-Stokes equations. The pressure distribution around a cylinder is first solved in Fourier space for different values of the parameters (wavenumber, confinement size). The method is applied to a simple geometry (cylinder in an axial flow with a variable confinement) in order to test its effectiveness. The empirical model is then compared to simulations and reference works in literature. The configuration with large confinement has been solved, and results were in agreement with Slender Body Theory. The dependency on the confinement size strongly depends on the wavenumber, but in any case added mass increases as the confinement size decreases. Finally, future perspectives to extend the model to a group of cylinders and to improve the model are discussed (i.e. add viscosity to the model).


2013 ◽  
Vol 444-445 ◽  
pp. 286-292
Author(s):  
Bing Han ◽  
Min Xu ◽  
Xi Pei ◽  
Xiao Min An

The effect of slender body on the rolling characteristics of a double delta wing is found by comparing the numerical simulation results of the double delta wing and wing-body configuration. The coupled computation system solving the Navier-Stokes equations and the rolling motion equation alternatively to obtain the unsteady vortical flow around the two configurations while rolling. The results conclusively showed the upwash effect of the slender body enhanced the energy of strake vortex and merged vortex.The aerodynamic lag of double delta wing is weak, contrarily, the time lag effect of the wing-body configuration is significant. The asymmetry vortices structure nearby the trailing edge are believed to be the main reason for the unsteady time lag effect.


2008 ◽  
Author(s):  
Abdus Samad ◽  
Ki-Sang Lee ◽  
Kwang-Yong Kim

This work presents a numerical optimization procedure for a low-speed axial flow fan blade with weighted average surrogate model. Reynolds-averaged Navier-Stokes equations with SST turbulence model are discretized by finite volume approximations and solved on hexahedral grids for flow analyses. The blade profile as well as stacking line is modified to enhance blade total efficiency, i.e., the objective function. Six design variables related to blade lean and blade profile are selected, and a design of experiments technique produces design points where flow analyses are performed to obtain values of the objective function. PBA model is employed as a surrogate model for optimization. A search algorithm is used to find the optimal design in the design space from the constructed surrogate model for the objective function. As a main result, the efficiency is increased effectively by the present optimization procedure.


Author(s):  
Youn-Sung Kim ◽  
Hyeon-Seok Shim ◽  
Kwang-Yong Kim

This study aims to evaluate effects of blade pitch and inlet guide vane (IGV) angle on the performance characteristics of a submersible axial-flow pump. According to the results of the previous study, the efficiency at the design and over-load conditions were significantly affected by the angle of IGV due to change in the incidence angle. To investigate the interactional effects of IGV and blade angle are analyzed using three-dimensional Reynolds-averaged Navier-Stokes equations with shear stress transport turbulence model. The hexahedral grids are used in the computational domain and a grid-dependency test is performed to obtain an optimal number of the grids. In this study, combinations of three different blade angles and two different IGV angles are tested. Adjusting angle of IGV increases the total pressure of the pump with a blade pitch increase, which can increase the efficiency of the pump in operating range.


2018 ◽  
Vol 848 ◽  
pp. 1040-1072 ◽  
Author(s):  
Chenguang Zhang ◽  
Harris Wong ◽  
Krishnaswamy Nandakumar

Fluid flow in microchannels has wide industrial and scientific applications. Hence, it is important to explore different driving mechanisms. In this paper, we study the net transport or fluid pumping in a two-dimensional channel induced by a travelling temperature wave applied at the bottom wall. The Navier–Stokes equations with the Boussinesq approximation and the convection–diffusion heat equation are made dimensionless by the height of the channel and a velocity scale obtained by a dominant balance between buoyancy and viscous resistance in the momentum equation. The system of equations is transformed to an axial coordinate that moves with the travelling temperature wave, and we seek steady solutions in this moving frame. Four dimensionless numbers emerge from the governing equations and boundary conditions: the Reynolds number $Re$, a Reynolds number $Rc$ based on the wave speed, the Prandtl number $Pr$ and the dimensionless wavenumber $K$. The system of equations is solved by a finite-volume method and by a perturbation method in the limit $Re\rightarrow 0$. Surprisingly, the leading and first-order perturbation solutions agree well with the computed axial flow for $Re\leqslant 10^{3}$. Thus, the analytic perturbation solutions are used to study systematically the effects of $Re$, $Rc$, $Pr$ and $K$ on the dimensionless induced axial flow $Q$. We find that $Q$ varies linearly with $Re$, and $Q/Re$ versus any of the three remaining dimensionless groups always exhibits a maximum. The global maximum of $Q/Re$ in the parameter space is subsequently determined for the first time. This induced axial flow is driven solely by the Reynolds stress.


Author(s):  
Bakhtier Farouk ◽  
Murat K. Aktas

Formation of vortical flow structures in a rectangular enclosure due to acoustic streaming is investigated numerically. The oscillatory flow field in the enclosure is created by the vibration of a vertical side wall of the enclosure. The frequency of the wall vibration is chosen such that a standing wave forms in the enclosure. The interaction of this standing wave with the horizontal solid walls leads to the production of Rayleigh type acoustic streaming flow patterns in the enclosure. All four walls of the enclosure considered are thermally insulated. The fully compressible form of the Navier-Stokes equations is considered and an explicit time-marching algorithm is used to explicitly track the acoustic waves. Numerical solutions are obtained by employing a highly accurate flux corrected transport (FCT) algorithm for the convection terms. A time-splitting technique is used to couple the viscous and diffusion terms of the full Navier-Stokes equations. Non-uniform grid structure is employed in the computations. The simulation of the primary oscillatory flow and the secondary (steady) streaming flows in the enclosure is performed. Streaming flow patterns are obtained by time averaging the primary oscillatory flow velocity distributions. The effect of the amount of wall displacement on the formation of the oscillatory flow field and the streaming structures are studied. Computations indicate that the nonlinearity of the acoustic field increases with increasing amount of the vibration amplitude. The form and the strength of the secondary flow associated with the oscillatory flow field and viscous effects are found to be strongly correlated to the maximum displacement of the vibrating wall. Total number of acoustic streaming cells per wavelength is also determined by the strength and the level of the nonlinearity of the sound field in the resonator.


Sign in / Sign up

Export Citation Format

Share Document