An official website of the United States government
Lubrication theory is adapted to incorporate the large normal stresses that occur for order-one Deborah numbers,$$De$$, the ratio of the relaxation time to the residence time. Comparing with the pressure drop for a Newtonian viscous fluid with a viscosity equal to that of an Oldroyd-B fluid in steady simple shear, we find numerically a reduced pressure drop through a contraction and an increased pressure drop through an expansion, both changing linearly with$$De$$at high$$De$$. For a constriction, there is a smaller pressure drop that plateaus at high$$De$$. For a contraction, much of the change in pressure drop occurs in the stress relaxation in a long exit channel. An asymptotic analysis for high$$De$$, based on the idea that normal stresses are stretched by an accelerating flow in proportion to the square of the velocity, reveals that the large linear changes in pressure drop are due to higher normal stresses pulling the fluid through the narrowest gap. A secondary cause of the reduction is that the elastic shear stresses do not have time to build up to their steady-state equilibrium value while they accelerate through a contraction. We find for a contraction or expansion that the high$$De$$analysis works well for$$De>0.4$$.
more »
« less
Flow of an Oldroyd-B fluid in a slowly varying contraction: theoretical results for arbitrary values of Deborah number in the ultra-dilute limit
Pressure-driven flows of viscoelastic fluids in narrow non-uniform geometries are common in physiological flows and various industrial applications. For such flows, one of the main interests is understanding the relationship between the flow rate$$q$$and the pressure drop$$\Delta p$$, which, to date, is studied primarily using numerical simulations. We analyse the flow of the Oldroyd-B fluid in slowly varying arbitrarily shaped, contracting channels and present a theoretical framework for calculating the$$q-\Delta p$$relation. We apply lubrication theory and consider the ultra-dilute limit, in which the velocity profile remains parabolic and Newtonian, resulting in a one-way coupling between the velocity and polymer conformation tensor. This one-way coupling enables us to derive closed-form expressions for the conformation tensor and the flow rate–pressure drop relation for arbitrary values of the Deborah number ($$De$$). Furthermore, we provide analytical expressions for the conformation tensor and the$$q-\Delta p$$relation in the high-Deborah-number limit, complementing our previous low-Deborah-number lubrication analysis. We reveal that the pressure drop in the contraction monotonically decreases with$$De$$, having linear scaling at high Deborah numbers, and identify the physical mechanisms governing the pressure drop reduction. We further elucidate the spatial relaxation of elastic stresses and pressure gradient in the exit channel following the contraction and show that the downstream distance required for such relaxation scales linearly with$$De$$.
more »
« less
- Award ID(s):
- 2246791
- PAR ID:
- 10569263
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 988
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The extensional rheology of dilute suspensions of spheres in viscoelastic/polymeric liquids is studied computationally. At low polymer concentration$$c$$and Deborah number$$\textit{De}$$(imposed extension rate times polymer relaxation time), a wake of highly stretched polymers forms downstream of the particles due to larger local velocity gradients than the imposed flow, indicated by$$\Delta \textit{De}_{\textit{local}}\gt 0$$. This increases the suspension’s extensional viscosity with time and$$\textit{De}$$for$$De \lt 0.5$$. When$$\textit{De}$$exceeds 0.5, the coil-stretch transition value, the fully stretched polymers from the far-field collapse in regions with$$\Delta \textit{De}_{\textit{local}} \lt 0$$(lower velocity gradient) around the particle’s stagnation points, reducing suspension viscosity relative to the particle-free liquid. The interaction between local flow and polymers intensifies with increasing$$c$$. Highly stretched polymers impede local flow, reducing$$\Delta \textit{De}_{\textit{local}}$$, while$$\Delta \textit{De}_{\textit{local}}$$increases in regions with collapsed polymers. Initially, increasing$$c$$aligns$$\Delta \textit{De}_{\textit{local}}$$and local polymer stretch with far-field values, diminishing particle–polymer interaction effects. However, beyond a certain$$c$$, a new mechanism emerges. At low$$c$$, fluid three particle radii upstream exhibits$$\Delta \textit{De}_{\textit{local}} \gt 0$$, stretching polymers beyond their undisturbed state. As$$c$$increases, however,$$\Delta \textit{De}_{\textit{local}}$$in this region becomes negative, collapsing polymers and resulting in increasingly negative stress from particle–polymer interactions at large$$\textit{De}$$and time. At high$$c$$, this negative interaction stress scales as$$c^2$$, surpassing the linear increase of particle-free polymer stress, making dilute sphere concentrations more effective at reducing the viscosity of viscoelastic liquids at larger$$\textit{De}$$and$$c$$.more » « less
-
Free surface flows driven by boundary undulations are observed in many biological phenomena, including the feeding and locomotion of water snails. To simulate the feeding strategy of apple snails, we develop a centimetric robotic undulator that drives a thin viscous film of liquid with the wave speed$$V_w$$. Our experimental results demonstrate that the behaviour of the net fluid flux$$Q$$strongly depends on the Reynolds number$$Re$$. Specifically, in the limit of vanishing$$Re$$, we observe that$$Q$$varies non-monotonically with$$V_w$$, which has been successfully rationalised by Pandeyet al.(Nat. Commun., vol. 14, no. 1, 2023, p. 7735) with the lubrication model. By contrast, in the regime of finite inertia ($${Re} \sim O(1)$$), the fluid flux continues to increase with$$V_w$$and completely deviates from the prediction of lubrication theory. To explain the inertia-enhanced pumping rate, we build a thin-film, two-dimensional model via the asymptotic expansion in which we linearise the effects of inertia. Our model results match the experimental data with no fitting parameters and also show the connection to the corresponding free surface shapes$$h_2$$. Going beyond the experimental data, we derive analytical expressions of$$Q$$and$$h_2$$, which allow us to decouple the effects of inertia, gravity, viscosity and surface tension on free surface pumping over a wide range of parameter space.more » « less
-
An analysis is presented for the axisymmetric lubrication resistance between permeable spherical particles. Darcy's law is used to describe the flow in the permeable medium and a slip boundary condition is applied at the interface. The pressure in the near-contact region is governed by a non-local integral equation. The asymptotic limit$$K=k/a^{2} \ll 1$$is considered, where$$k$$is the arithmetic mean permeability, and$$a^{-1}=a^{-1}_{1}+a^{-1}_{2}$$is the reduced radius, and$$a_1$$and$$a_2$$are the particle radii. The formulation allows for particles with distinct particle radii, permeabilities and slip coefficients, including permeable and impermeable particles and spherical drops. Non-zero particle permeability qualitatively affects the axisymmetric near-contact motion, removing the classical lubrication singularity for impermeable particles, resulting in finite contact times under the action of a constant force. The lubrication resistance becomes independent of gap and attains a maximum value at contact$$F=6{\rm \pi} \mu a W K^{-2/5}\tilde {f}_c$$, where$$\mu$$is the fluid viscosity,$$W$$is the relative velocity and$$\tilde {f}_c$$depends on slip coefficients and weakly on permeabilities; for two permeable particles with no-slip boundary conditions,$$\tilde {f}_c=0.7507$$; for a permeable particle in near contact with a spherical drop,$$\tilde {f}_c$$is reduced by a factor of$$2^{-6/5}$$.more » « less
-
A lubrication analysis is presented for the resistances between permeable spherical particles in near contact,$$h_0/a\ll 1$$, where$$h_0$$is the minimum separation between the particles, and$$a=a_1 a_2/(a_1+a_2)$$is the reduced radius. Darcy's law is used to describe the flow inside the permeable particles and no-slip boundary conditions are applied at the particle surfaces. The weak permeability regime$$K=k/a^{2} \ll 1$$is considered, where$$k=\frac {1}{2}(k_1+k_2)$$is the mean permeability. Particle permeability enters the lubrication resistances through two functions of$$q=K^{-2/5}h_0/a$$, one describing axisymmetric motions, the other transverse. These functions are obtained by solving an integral equation for the pressure in the near-contact region. The set of resistance functions thus obtained provide the complete set of near-contact resistance functions for permeable spheres and match asymptotically to the standard hard-sphere resistances that describe pairwise hydrodynamic interactions away from the near-contact region. The results show that permeability removes the contact singularity for non-shearing particle motions, allowing rolling without slip and finite separation velocities between touching particles. Axisymmetric and transverse mobility functions are presented that describe relative particle motion under the action of prescribed forces and in linear flows. At contact, the axisymmetric mobility under the action of oppositely directed forces is$$U/U_0=d_0K^{2/5}$$, where$$U$$is the relative velocity,$$U_0$$is the velocity in the absence of hydrodynamic interactions and$$d_0=1.332$$. Under the action of a constant tangential force, a particle in contact with a permeable half-space rolls without slipping with velocity$$U/U_0=d_1(d_2+\log K^{-1})^{-1}$$, where$$d_1=3.125$$and$$d_2=6.666$$; in shear flow, the same expression holds with$$d_1=7.280$$.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2024.223
