No Arabic abstract
We report a novel technique to passively create strong secondary flows at moderate to high flow rates in microchannels, accurately control them and finally, due to their deterministic nature, program them into microfluidic platforms. Based on the flow conditions and due to the presence of the pillars in the channel, the flow streamlines will lose their fore-aft symmetry. As a result of this broken symmetry the fluid is pushed away from the pillar at the center of the channel (i.e. central z-plane). As the flow needs to maintain conservation of mass, the fluid will laterally travel in the opposite direction near the top and bottom walls. Therefore, a NET secondary flow will be created in the channel cross-section which is depicted in this video. The main platform is a simple straight channel with posts (i.e. cylindrical pillars - although other pillar cross-sections should also function) placed along the channel. Channel measures were 200 mumtimes50 mum, with pillars of 100 mum in diameter. Positioning the pillars in different locations within the cross-section of the channel will result in induction of different secondary flow patterns, which can be carefully engineered. The longitudinal spacing of the pillars is another design parameter (600 mum spacing was used for this video). The device works over a wide range (moderate to high) flow rates. We used 150 muL/min in this experiment. The device has 3 inlets where a dye stream is co-flowed between two water streams. In this video, one can see the effect of the net secondary flow created by inertia in the microchannel by visualizing the cross-section of the fluorescently labeled stream. Confocal images are sequentially taken at the inlet and after 49 consecutive pillars.
Experimental and numerical investigations are performed to provide an assessment of the transport behavior of an ultrasonic oscillatory two-phase flow in a microchannel. The work is inspired by the flow observed in an innovative ultrasonic fabric drying device using a piezoelectric bimorph transducer with microchannels, where a water-air two-phase flow is transported by harmonically oscillating microchannels. The flow exhibits highly unsteady behavior as the water and air interact with each other during the vibration cycles, making it significantly different from the well-studied steady flow in microchannels. The computational fluid dynamics (CFD) modeling is realized by combing the turbulence Reynolds-averaged Navier-Stokes (RANS) k-${omega}$ model with the phase-field method to resolve the dynamics of the two-phase flow. The numerical results are qualitatively validated by the experiment. Through parametric studies, we specifically examined the effects of vibration conditions (i.e., frequency and amplitude), microchannel taper angle, and wall surface contact angle (i.e., wettability) on the flow rate through the microchannel. The results will advance the potential applications where oscillatory or general unsteady microchannel two-phase flows may be present.
Long, shallow microchannels embedded in thick soft materials are widely used in microfluidic devices for lab-on-a-chip applications. However, the bulging effect caused by fluid--structure interactions between the internal viscous flow and the soft walls has not been completely understood. Previous models either contain a fitting parameter or are specialized to channels with plate-like walls. This work is a theoretical study of the steady-state response of a compliant microchannel with a thick wall. Using lubrication theory for low-Reynolds-number flows and the theory for linearly elastic isotropic solids, we obtain perturbative solutions for the flow and deformation. Specifically, only the channels top wall deformation is considered, and the ratio between its thickness $t$ and width $w$ is assumed to be $(t/w)^2 gg 1$. We show that the deformation at each stream-wise cross-section can be considered independently, and that the top wall can be regarded as a simply supported rectangle subject to uniform pressure at its bottom. The stress and displacement fields are found using Fourier series, based on which the channel shape and the hydrodynamic resistance are calculated, yielding a new flow rate--pressure drop relation without fitting parameters. Our results agree favorably with, and thus rationalize, previous experiments.
In fully-developed pressure-driven flow, the spreading of a dissolved solute is enhanced in the flow direction due to transverse velocity variations in a phenomenon now commonly referred to as Taylor-Aris dispersion. It is well understood that the characteristics of the dispersion are sensitive to the channels cross-sectional geometry. Here we demonstrate a method for manipulation of dispersion in a single rectangular microchannel via controlled deformation of its upper wall. Using a rapidly prototyped multi-layer microchip, the channel wall is deformed by a controlled pressure source allowing us to characterize the dependence of the dispersion on the deflection of the channel wall and overall channel aspect ratio. For a given channel aspect ratio, an optimal deformation to minimize dispersion is found, consistent with prior numerical and theoretical predictions. Our experimental measurements are also compared directly to numerical predictions using an idealized geometry.
We simulated two particle-based fluid models, namely multiparticle collision dynamics and dissipative particle dynamics, under shear using reverse nonequilibrium simulations (RNES). In cubic periodic simulation boxes, the expected shear flow profile for a Newtonian fluid developed, consistent with the fluid viscosities. However, unexpected secondary flows along the shear gradient formed when the simulation box was elongated in the flow direction. The standard shear flow profile was obtained when the simulation box was longer in the shear-gradient dimension than the flow dimension, while the secondary flows were always present when the flow dimension was at least 25% larger than the shear-gradient dimension. The secondary flows satisfy the boundary conditions imposed by the RNES and have a lower rate of viscous dissipation in the fluid than the corresponding unidirectional flows. This work highlights a previously unappreciated limitation of RNES for generating shear flow in simulation boxes that are elongated in the flow dimension, an important consideration when applying RNES to complex fluids like polymer solutions.
The aim of the present work is to investigate the role of coherent structures in the generation of the secondary flow in a turbulent square duct. The coherent structures are defined as connected regions of flow where the product of the instantaneous fluctuations of two velocity components is higher than a threshold based on the long-time turbulence statistics, in the spirit of the three-dimensional quadrant analysis proposed by Lozano-Duran et al. (J. Fluid Mech., vol. 694, 2012, pp. 100-130). We consider both the direct contribution of the structures to the mean in-plane velocity components and their geometrical properties. The instantaneous phenomena taking place in the turbulent duct are compared with turbulent channel flow at Reynolds numbers of $Re_tau=180$ and $360$, based on friction velocity at the center-plane and channel half height. In the core region of the duct, the fractional contribution of intense events to the wall-normal component of the mean velocity is in very good agreement with that in the channel, despite the presence of the secondary flow in the former. Additionally, the shapes of the three-dimensional objects do not differ significantly in both flows. On the other hand, in the corner region of the duct, the proximity of the walls affects both the geometrical properties of the coherent structures and the contribution to the mean component of the vertical velocity, which is less relevant than that of the complementary portion of the flow not included in such objects. Our results show however that strong Reynolds shear-stress events, despite the differences observed between channel and duct, do not contribute directly to the secondary motion, and thus other phenomena need to be considered instead.