Influence of Three-Dimensional Roughness on Pressure-Driven Flow Through Microchannels

2003 ◽  
Vol 125 (5) ◽  
pp. 871-879 ◽  
Author(s):  
Yandong Hu ◽  
Carsten Werner ◽  
Dongqing Li

Surface roughness is present in most of the microfluidic devices due to the microfabrication techniques or particle adhesion. It is highly desirable to understand the roughness effect on microscale flow. In this study, we developed a three-dimensional finite-volume-based numerical model to simulate pressure-driven liquid flow in microchannels with rectangular prism rough elements on the surfaces. Both symmetrical and asymmetric roughness element arrangements were considered, and the influence of the roughness on pressure drop was examined. The three-dimensional numerical solution shows significant effects of surface roughness in terms of the rough elements’ height, size, spacing, and the channel height on both the velocity distribution and the pressure drop. The compression-expansion flow around the three-dimensional roughness elements and the flow blockage caused by the roughness in the microchannel were discussed. An expression of the relative channel height reduction due to roughness effect was presented.

2017 ◽  
Vol 832 ◽  
pp. 793-824 ◽  
Author(s):  
Kiran Bhaganagar ◽  
Narasimha Rao Pillalamarri

A fundamental study has been conducted to understand the front characteristics and the mixing in the flow of density currents over rough surfaces. A large-eddy simulation (LES) has been performed for lock-exchange release density currents over rough walls to shed light on the unsteady mixing processes. A volume-penalization method, which is a special case of the immersed-boundary method, has been implemented to realize the bottom-mounted rough topology. In this study, the LES has been conducted in a channel with a lower wall covered with three-dimensional cube- and pyramid-shaped roughness elements, such that all cases have the same base area, but differences in the roughness solidity and volume fraction of roughness. Both cases of identical roughness elements and those with randomness in height have been considered. The maximum roughness height for all cases is kept at a constant fraction (10 %) of the total channel height. The study focuses on the instantaneous mixing processes in lock-exchange release currents over rough surfaces. An important contribution of the work is that qualitative and quantitative analysis has been conducted to demonstrate additional mixing mechanisms due to the presence of surface roughness that enhances dilution of the current. Enhanced mixing due to roughness is related to the strength of the shear layer resulting from the roughness, and hence depends on friction Reynolds number ($Re_{\unicode[STIX]{x1D70F}}$). The combined role of current characteristics and $Re_{\unicode[STIX]{x1D70F}}$ together dictate the mixing processes and extent of dilution in density currents over surface roughness.


Author(s):  
Ainul Haque ◽  
Ameeya Kumar Nayak

In this paper, a mathematical model has been developed to analyze the combined electroosmotic and pressure driven flow of power law fluids in a micro channel in the presence of Joule heating effects. The effects of Navier slip boundary condition and thermal radiation is examined for effective heat transfer in a hydrophobic microchannel. The analytical treatment has been performed for fluid flow and heat transfer effects in terms of flow governing parameters. This study highlights the effect of channel height to the electric double layer thickness and observed the flow variation due to heat transfer effect with the available scientific data. For a pure EOF, velocity slip have more significant role to get a maximum flow rate as expected. For both pseudo-plastic and dilatent fluids Nusselt number is decreased with the increment of the hydrophobic parameter and dimensionless pressure gradient where as increment in Joule heating effect enhance the heat transfer rate.


Author(s):  
Stepan Tolkachev ◽  
Victor Kozlov ◽  
Valeriya Kaprilevskaya

In this article, the results of research about stationary and secondary disturbances development behind the localized and two-dimensional roughness elements are presented. It is shown that the two-dimensional roughness element has a destabilizing effect on the disturbances induced by the three-dimensional roughness element lying upstream. In this case, the two-dimensional roughness element causes the appearance of stationary structures, and then secondary perturbations, whose frequency range lies lower than in the case of the stationary vortices excited by a three-dimensional roughness element.


2019 ◽  
Vol 141 (12) ◽  
Author(s):  
Bayode E. Owolabi ◽  
David J. C. Dennis ◽  
Robert J. Poole

In this study, we examine the development length requirements for laminar Couette–Poiseuille flows in a two-dimensional (2D) channel as well as in the three-dimensional (3D) case of flow through a square duct, using a combination of numerical and experimental approaches. The parameter space investigated covers wall to bulk velocity ratios, r, spanning from 0 (purely pressure-driven flow) to 2 (purely wall driven-flow; 4 in the case of a square duct) and a wide range of Reynolds numbers (Re). The results indicate an increase in the development length (L) with r. Consistent with the findings of Durst et al. (2005, “The Development Lengths of Laminar Pipe and Channel Flows,” ASME J. Fluids Eng., 127(6), pp. 1154–1160), L was observed to be of the order of the channel height in the limit as Re→0, irrespective of the condition at the inlet. This, however, changes at high Reynolds numbers, with L increasing linearly with Re. In all the cases considered, a uniform velocity profile at the inlet was found to result in longer entry lengths than in a flow developing from a parabolic inlet profile. We show that this inlet effect becomes less important as the limit of purely wall-driven flow is approached. Finally, we develop correlations for predicting L in these flows and, for the first time, also present laser Doppler velocimetry (LDV) measurements of the developing as well as fully-developed velocity profiles, and observe good agreement between experiment, analytical solution, and numerical simulation results in the 3D case.


2011 ◽  
Vol 679 ◽  
pp. 77-100 ◽  
Author(s):  
QIANLONG LIU ◽  
ANDREA PROSPERETTI

The finite-Reynolds-number three-dimensional flow in a channel bounded by one and two parallel porous walls is studied numerically. The porous medium is modelled by spheres in a simple cubic arrangement. Detailed results on the flow structure and the hydrodynamic forces and couple acting on the sphere layer bounding the porous medium are reported and their dependence on the Reynolds number illustrated. It is shown that, at finite Reynolds numbers, a lift force acts on the spheres, which may be expected to contribute to the mobilization of bottom sediments. The results for the slip velocity at the surface of the porous layers are compared with the phenomenological Beavers–Joseph model. It is found that the values of the slip coefficient for pressure-driven and shear-driven flow are somewhat different, and also depend on the Reynolds number. A modification of the relation is suggested to deal with these features. The Appendix provides an alternative derivation of this modified relation.


2010 ◽  
Vol 660 ◽  
pp. 527-537 ◽  
Author(s):  
NILS TILTON ◽  
DENIS MARTINAND ◽  
ERIC SERRE ◽  
RICHARD M. LUEPTOW

A generalized solution for pressure-driven flow through a permeable rotating inner cylinder in an impermeable concentric outer cylinder, a situation used commercially in rotating filtration, is challenging due to the interdependence between the pressure drop in the axial direction and that across the permeable inner cylinder. Most previous approaches required either an imposed radial velocity at the inner cylinder or radial throughflow with both the inner and outer cylinders being permeable. We provide an analytical solution for rotating Couette–Poiseuille flow with Darcy's law at the inner cylinder by using a small parameter related to the permeability of the inner cylinder. The theory works for suction, injection and even combined suction/injection, when the axial pressure drop in the annulus is such that the transmembrane pressure difference reverses sign along the axial extent of the system. Corresponding numerical simulations for finite-length systems match the theory very well.


2012 ◽  
Vol 1 (2) ◽  
pp. 41 ◽  
Author(s):  
D. Barchiesi ◽  
S. Kessentini

The fabrication process of nanodevices is continually improved. However, most of the nanodevices, such as biosensors present rough surfaces with mean roughness of some nanometers even if the deposition rate of material is more controlled. The effect of roughness on performance of biosensors was fully addressed for plane biosensors and gratings, but rarely addressed for biosensors based on Local Plasmon Resonance. The purpose of this paper is to evaluate numerically the influence of nanometric roughness on the efficiency of a dimer nano-biosensor (two levels of roughness are considered). Therefore, we propose a general numerical method, that can be applied to any other nanometric shape, to take into account the roughness in a three dimensional model. The study focuses on both the far-field, which corresponds to the experimental detected data, and the near-field, responsible for exciting and then detecting biological molecules. The results suggest that the biosensor efficiency is highly sensitive to the surface roughness. The roughness can produce important shifts of the extinction efficiency peak and a decrease of its amplitude resulting from changes in the distribution of near-field and absorbed electric field intensities.


2013 ◽  
Vol 724 ◽  
pp. 642-670 ◽  
Author(s):  
S. Cherubini ◽  
M. D. De Tullio ◽  
P. De Palma ◽  
G. Pascazio

AbstractThis work provides a global optimization analysis, looking for perturbations inducing the largest energy growth at a finite time in a boundary-layer flow in the presence of smooth three-dimensional roughness elements. Amplification mechanisms are described which can bypass the asymptotical growth of Tollmien–Schlichting waves. Smooth axisymmetric roughness elements of different height have been studied, at different Reynolds numbers. The results show that even very small roughness elements, inducing only a weak deformation of the base flow, can localize the optimal disturbance characterizing the Blasius boundary-layer flow. Moreover, for large enough bump heights and Reynolds numbers, a strong amplification mechanism has been recovered, inducing an increase of several orders of magnitude of the energy gain with respect to the Blasius case. In particular, the highest value of the energy gain is obtained for an initial varicose perturbation, differently to what found for a streaky parallel flow. Optimal varicose perturbations grow very rapidly by transporting the strong wall-normal shear of the base flow, which is localized in the wake of the bump. Such optimal disturbances are found to lead to transition for initial energies and amplitudes considerably smaller than sinuous optimal ones, inducing hairpin vortices downstream of the roughness element.


2016 ◽  
Vol 796 ◽  
pp. 158-194 ◽  
Author(s):  
Holger B. E. Kurz ◽  
Markus J. Kloker

The effects of a spanwise row of finite-size cylindrical roughness elements in a laminar, compressible, three-dimensional boundary layer on a wing profile are investigated by direct numerical simulations (DNS). Large elements are capable of immediately tripping turbulent flow by either a strong, purely convective or an absolute/global instability in the near wake. First we focus on an understanding of the steady near-field past a finite-size roughness element in the swept-wing flow, comparing it to a respective case in unswept flow. Then, the mechanisms leading to immediate turbulence tripping are elaborated by gradually increasing the roughness height and varying the disturbance background level. The quasi-critical roughness Reynolds number above which turbulence sets in rapidly is found to be $Re_{kk,qcrit}\approx 560$ and global instability is found only for values well above 600 using nonlinear DNS; therefore the values do not differ significantly from two-dimensional boundary layers if the full velocity vector at the roughness height is taken to build $Re_{kk}$. A detailed simulation study of elements in the critical range indicates a changeover from a purely convective to a global instability near the critical height. Finally, we perform a three-dimensional global stability analysis of the flow field to gain insight into the early stages of the temporal disturbance growth in the quasi-critical and over-critical cases, starting from a steady state enforced by damping of unsteady disturbances.


Author(s):  
Hongjun Song ◽  
Yi Wang ◽  
Kapil Pant

The utilization of cross-stream diffusion under laminar flow for precise analyte handling plays a critical role in microfluidic biochemical assays such as sample preparation, concentration gradient generation, and molecular interactions. The non-uniform velocity profile along the cross-section of a rectangular microchannel with arbitrary aspect ratio under pressure-driven flow results in unique, heterogeneous species transport including Taylor dispersion and position-dependent diffusion scaling law. Although numerical methods such as finite difference method, finite element method, the method of lines and lattice Boltzmann (LB) method have been used for quantitative study of the phenomena, they inherently suffer from several limitations, such as difficulty to provide direct, physical insight into the underlying transport mechanism and prohibitive computational cost to suppress the artificial numerical diffusion (ND). To address these issues, several analytical models have been proposed, which share several common assumptions such as large aspect ratio and neglecting depth-wise diffusion due to the non-uniform axial velocity in the 3D convection-diffusion equation, markedly limiting their utility. In this paper, we present a three dimensional (3D) analytical model to investigate the diffusion of analyte between two cross streams in rectangular microchannels with arbitrary aspect ratios under pressure-driven flow. The 3D convection-diffusion equation is solved in a Fourier series form using a double integral transformation method and associated eigensystem calculation. Therefore, the model for the first time is capable of capturing the non-uniform transport rate (i.e., the ‘butterfly effect’) and the position-dependent scaling-law of diffusion (1/3-power at the channel wall and 1/2-pwer at the half-depth plane) through an analytical solution. Our analytical model was extensively validated against both experimental and numerical data in terms of the concentration distribution, diffusion scaling law and the mixing efficiency with excellent agreement (the relative error is much less than 0.5% in various benchmark test cases.) Quantitative comparison between our analytical model and other prior analytical models in extensive parameter space was also performed, which convincingly demonstrates that our model accommodates much broader transport regimes and more practical microfluidic applications.


Sign in / Sign up

Export Citation Format

Share Document