Search

We analyze the steady viscoelastic fluid flow in slowly varying contracting channels of arbitrary shape and present a theory based on the lubrication approximation for calculating the flow rate-pressure drop relation at low and high Deborah ($De$) numbers. Unlike most prior theoretical studies leveraging the Oldroyd-B model, we describe the fluid viscoelasticity using a FENE-CR model and examine how the polymer chains' finite extensibility impacts the pressure drop. We employ the low-Deborah-number lubrication analysis to provide analytical expressions for the pressure drop up to $O(De^4)$. We further consider the ultra-dilute limit and exploit a one-way coupling between the parabolic velocity and elastic stresses to calculate the pressure drop of the FENE-CR fluid for arbitrary values of the Deborah number. Such an approach allows us to elucidate elastic stress contributions governing the pressure drop variations and the effect of finite extensibility for all $De$. We validate our theoretical predictions with two-dimensional numerical simulations and find excellent agreement. We show that, at low Deborah numbers, the pressure drop of the FENE-CR fluid monotonically decreases with $De$, similar to the previous results for the Oldroyd-B and FENE-P fluids. However, at high Deborah numbers, in contrast to a linear decrease for the Oldroyd-B fluid, the pressure drop of the FENE-CR fluid exhibits a non-monotonic variation due to finite extensibility, first decreasing and then increasing with $De$. Nevertheless, even at sufficiently high Deborah numbers, the pressure drop of the FENE-CR fluid in the ultra-dilute and lubrication limits is lower than the corresponding Newtonian pressure drop.

We provide an experimental framework to measure the flow rate--pressure drop relation for Newtonian and shear-thinning fluids in two common deformable configurations: (\textit{i}) a rectangular channel and (\textit{ii}) an axisymmetric tube. Using the Carreau model to describe the shear-dependent viscosity, we identify the key dimensionless rheological number, $Cu$, which characterizes shear thinning, and we show that our experiments lie within the power-law regime of shear rates. To rationalize the experimental data, we derive the flow rate-pressure drop relation taking into account the two-way-coupled fluid-structure interaction between the flow and its compliant confining boundaries. We thus identify the second key dimensionless number, $\alpha$, which characterizes the compliance of the conduit. We then compare the theoretical flow rate-pressure drop relation to our experimental measurements, finding excellent agreement between the two. We further contrast our results for shear-thinning and Newtonian fluids to highlight the influence of $Cu$ on the flow rate-pressure drop relation. Finally, we delineate four distinct physical regimes of flow and deformation by mapping our experimental flow rate-pressure drop data for Newtonian and shear-thinning fluids into a $Cu-\alpha$ plane., Comment: 10 pages, 4 figures; v2 accepted for publication in Physical Review Fluids

We demonstrate the use of the Lorentz reciprocal theorem in obtaining corrections to the flow rate for oscillatory flow in channels. We start from the unsteady Stokes equations and derive the suitable reciprocity relations for both rigid and deformable channels, assuming all quantities can be expressed as time-harmonic phasors. For the case of a rigid channel, the auxiliary problem is a steady flow solution, and the reciprocal theorem allows us to calculate the first-order correction in the Womersley number to the steady flow rate. For the case of a deformable channel, we account for the fluid--structure interaction using domain perturbation under the assumption of a small compliance number. Using the oscillatory flow solution in a rigid channel as the auxiliary problem, we obtain the first-order correction in the compliance number to the rigid channel flow rate. This correction, validated against numerical simulations, provides the leading-order effect of the interplay between the oscillations of the fluid flow and the compliance of the channel for arbitrary values of the Womersley number., Comment: 13 pages, 4 figures

We present a direct derivation of the typical time derivatives used in a continuum description of complex fluid flows, harnessing the principles of the kinematics of line elements. The evolution of the microstructural conformation tensor in a flow and the physical interpretation of different derivatives then follow naturally., Comment: 1 figure

We analyze the steady non-Newtonian fluid-structure interaction between the flow of an Oldroyd-B fluid and a deformable channel. Specifically, we provide a theoretical framework for calculating the leading-order effect of the fluid's viscoelasticity on the flow rate-pressure drop relation and on the deformation of the channel's elastic wall. We first identify the characteristic scales and dimensionless parameters governing the fluid-structure interaction in slender and shallow channels. Applying the lubrication approximation for the flow and employing a perturbation expansion in powers of the Deborah number $De$, we derive a closed-form expression for the pressure as a function of the non-uniform shape of the channel in the weakly viscoelastic limit up to $\mathrm{O}(De)$. Coupling the hydrodynamic pressure to the elastic deformation, we provide the leading-order effect of the interplay between the viscoelasticity of the fluid and the compliance of the channel on the pressure and deformation fields, as well as on the flow rate-pressure drop relation. For the flow-rate-controlled regime and in the weakly viscoelastic limit, we show analytically that both the compliance of the deforming top wall and the viscoelasticity of the fluid decrease the pressure drop. Furthermore, we reveal a trade-off between the influence of compliance of the channel and the fluid's viscoelasticity on the deformation. While the channel's compliance increases the deformation, the fluid's viscoelasticity decreases it., Comment: 5 figures

Viscous flows through configurations manufactured from soft materials apply both pressure and shear stress at the solid-liquid interface, leading to deformation of the cross-section, which affects the flow rate-pressure drop relation. Conventionally, calculating this flow rate-pressure drop relation requires solving the complete elastohydrodynamic problem, which couples the fluid flow and elastic deformation. In this work, we use the reciprocal theorems for Stokes flow and linear elasticity to derive a closed-form expression for the flow rate-pressure drop relation in deformable channels, bypassing the detailed calculation of the solution to the fluid-structure-interaction problem. For small deformations (under a domain perturbation scheme), our theory provides the leading-order effect, of the interplay between the fluid stresses and the compliance of the channel, on the flow rate-pressure drop relation. Our approach uses solely the fluid flow solution and the elastic deformation due to the fluid stress distribution in an undeformed channel, eliminating the need to solve the coupled elastohydrodynamic problem. Unlike previous theoretical studies that neglected the presence of lateral sidewalls and considered shallow geometries of effectively infinite width, our approach allows to determine the influence of confining sidewalls on the flow rate-pressure drop relation. For the flow-rate-controlled situation and the plate-bending theory for the elastic deformation, we show a trade-off between the effect of compliance of the deforming top wall and the drag due to sidewalls on the pressure drop. While increased compliance decreases the pressure drop, the effect of the sidewalls increases it. Our theoretical framework may provide insight into existing experimental data and pave the way for the design of novel optimized soft microfluidic configurations of different cross-sectional shapes., Comment: 9 pages, 2 figures

We analyze the pressure-driven flow of a viscoelastic fluid in arbitrarily shaped, narrow channels and present a theoretical framework for calculating the relationship between the flow rate $q$ and pressure drop $\Delta p$. We utilize the Oldroyd-B model and first identify the characteristic scales and dimensionless parameters governing the flow in the lubrication limit. Employing a perturbation expansion in powers of the Deborah number ($De$), we provide analytical expressions for the velocity, stress, and the $q-\Delta p$ relation in the weakly viscoelastic limit up to $O(De^2)$. Furthermore, we exploit the reciprocal theorem derived by Boyko $\&$ Stone (Phys. Rev. Fluids, vol. 6, 2021, pp. L081301) to obtain the $q-\Delta p$ relation at the next order, $O(De^3)$, using only the velocity and stress fields at the previous orders. We validate our analytical results with two-dimensional numerical simulations of the Oldroyd-B fluid in a hyperbolic, symmetric contracting channel and find excellent agreement. For the flow-rate-controlled situation, both our theory and simulations reveal weak dependence of the velocity field on the Deborah number, so that the velocity can be approximated as Newtonian. In contrast to the velocity, the pressure drop strongly depends on the viscoelastic effects and decreases with $De$. Elucidating the relative importance of different terms in the momentum equation contributing to the pressure drop, we identify that a pressure drop reduction for narrow contracting geometries is primarily due to gradients in the viscoelastic shear stresses, while viscoelastic axial stresses have a minor effect on the pressure drop along the symmetry line., Comment: 28 pages, 9 figures

We present a theoretical model and experimental demonstration of thin liquid film deformations due to a dielectric force distribution established by surface electrodes. We model the spatial electric field produced by a pair of parallel electrodes and use it to evaluate the stress on the interface through Maxwell stresses. By coupling this force with the Young-Laplace equation, we obtain the deformation of the interface. To validate our theory, we design an experimental setup which uses microfabricated electrodes to achieve spatial dielectrophoretic actuation of a thin liquid film, while providing measurements of microscale deformations through digital holographic microscopy. We characterize the deformation as a function of the electrode-pair geometry and film thickness, showing very good agreement with the model. Based on the insights from the characterization of the system, we pattern conductive lines of electrode pairs on the surface of a microfluidic chamber and demonstrate the ability to produce complex two-dimensional deformations. The films can remain in liquid form and be dynamically modulated between different configurations or polymerized to create solid structures with high surface quality., Comment: 14 pages and 7 figures in the main manuscript. 4 sections and 5 figures in the SI

We demonstrate the existence of a fluid-structure instability arising from the interaction of electro-osmotic flow with an elastic substrate. Considering the case of flow within a soft fluidic chamber, we show that above a certain electric field threshold, negative gauge pressure induced by electro-osmotic flow causes the collapse of its elastic walls. We combine experiments and theoretical analysis to elucidate the underlying mechanism for instability and identify several distinct dynamic regimes. The understanding of this instability is important for the design of electrokinetic systems containing soft elements., Comment: 10 pages, 4 figures

We experimentally demonstrate the phenomenon of electroosmotic dipole flow that occurs around a localized surface charge region under the application of an external electric field in a Hele-Shaw cell. We use localized deposition of polyelectrolytes to create well-controlled surface charge variations, and show that for a disk-shaped spot, the internal pressure distribution that arises, results in uniform flow within the spot and dipole flow around it. We further demonstrate the superposition of surface charge spots to create complex flow patterns, without the use of physical walls., Comment: Manuscript 9 pages, 5 figures; Supplementary Material 9 pages, 2 supplementary figures

In a wide range of applications, microfluidic channels are implemented in soft substrates. In such configurations, where fluidic inertia and compressibility are negligible, the propagation of fluids in channels is governed by a balance between fluid viscosity and elasticity of the surrounding solid. The viscous-elastic interactions between elastic substrates and non-Newtonian fluids are particularly of interest due to the dependence of viscosity on the state of the system. In this work, we study the fluid-structure interaction dynamics between an incompressible non-Newtonian fluid and a slender linearly elastic cylinder under the creeping flow regime. Considering power-law fluids and applying the thin shell approximation for the elastic cylinder, we obtain a non-homogeneous p-Laplacian equation governing the viscous-elastic dynamics. We present exact solutions for the pressure and deformation fields for various initial and boundary conditions for both shear-thinning and shear-thickening fluids. We show that in contrast to Stokes' problem where a compactly supported front is obtained for shear-thickening fluids, here the role of viscosity is inversed and such fronts are obtained for shear-thinning fluids. Furthermore, we demonstrate that for the case of a step in inlet pressure, the propagation rate of the front has a $t^{\frac{n}{n+1}}$ dependence on time ($t$), suggesting the ability to indirectly measure the power-law index ($n$) of shear-thinning liquids through measurements of elastic deformation., Comment: 27 pages, 5 figures

The interaction of a thin viscous film with an elastic sheet results in coupling of pressure and deformation, which can be utilized as an actuation mechanism for surface deformations in a wide range of applications, including microfluidics, optics, and soft robotics. Implementation of such configurations inherently takes place over finite domains and often requires some pre-stretching of the sheet. Under the assumptions of strong pre-stretching and small deformations of the lubricated elastic sheet, we use the linearized Reynolds and Foppl-von Karman equations to derive closed-form analytical solutions describing the deformation in a finite domain due to external forces, accounting for both bending and tension effects. We provide a closed-form solution for the case of a square-shaped actuation region and present the effect of pre-stretching on the dynamics of the deformation. We further present the dependence of the deformation magnitude and timescale on the spatial wavenumber, as well as the transition between stretching- and bending-dominant regimes. We also demonstrate the effect of spatial discretization of the forcing (representing practical actuation elements) on the achievable resolution of the deformation. Extending the problem to an axisymmetric domain, we investigate the effects arising from nonlinearity of the Reynolds and Foppl-von Karman equations and present the deformation behavior as it becomes comparable to the initial film thickness and dependent on the induced tension. These results set the theoretical foundation for implementation of microfluidic soft actuators based on elastohydrodynanmics., Comment: 19 pages, 9 figures

We analyze flow of non-Newtonian fluids in a Hele-Shaw cell, subjected to spatially non-uniform electroosmotic slip. Motivated by their potential use for increasing the characteristic pressure fields, we specifically focus on power-law fluids with wall depletion properties. We derive a p-Poisson equation governing the pressure field, as well as a set of linearized equations representing its asymptotic approximation for weakly non-Newtonian behavior. To investigate the effect of non-Newtonian properties on the resulting fluidic pressure and velocity, we consider several configurations in one- and two-dimensions, and calculate both exact and approximate solutions. We show that the asymptotic approximation is in good agreement with exact solutions even for fluids with significant non-Newtonian behavior, allowing its use in the analysis and design of microfluidic systems involving electro-kinetic transport of such fluids., Comment: 20 pages, 7 figures

We present an analytical study, validated by numerical simulations, of electroosmotic flow in a Hele-Shaw cell with non-uniform surface charge patterning. Applying the lubrication approximation and assuming thin electric double layer, we obtain a pair of uncoupled Poisson equations which relate the pressure and the stream function, respectively, to gradients in the zeta potential distribution parallel and perpendicular to the applied electric field. We solve the governing equations for the fundamental case of a disk with uniform zeta potential and show that the flow-field in the outer region takes the form of a pure dipole. We illustrate the ability to generate complex flow-fields around smooth convex regions by superposition of such disks with uniform zeta potential and a uniform pressure driven flow. This method may be useful for future on-chip devices, allowing flow control without the need for mechanical components., Comment: Accepted to Physics of Fluids