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
-
-
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
-
Large-eddy simulation was used to model turbulent atmospheric surface layer (ASL) flow over canopies composed of streamwise-aligned rows of synthetic trees of height,$$h$$, and systematically arranged to quantify the response to variable streamwise spacing,$$\delta _1$$, and spanwise spacing,$$\delta _2$$, between adjacent trees. The response to spanwise and streamwise heterogeneity has, indeed, been the topic of a sustained research effort: the former resulting in formation of Reynolds-averaged counter-rotating secondary cells, the latter associated with the$$k$$- and$$d$$-type response. No study has addressed the confluence of both, and results herein show secondary flow polarity reversal across ‘critical’ values of$$\delta _1$$and$$\delta _2$$. For$$\delta _2/\delta \lesssim 1$$and$$\gtrsim 2$$, where$$\delta$$is the flow depth, the counter-rotating secondary cells are aligned such that upwelling and downwelling, respectively, occurs above the elements. The streamwise spacing$$\delta _1$$regulates this transition, with secondary cell reversal occurring first for the largest$$k$$-type cases, as elevated turbulence production within the canopy necessitates entrainment of fluid from aloft. The results are interpreted through the lens of a benchmark prognostic closure for effective aerodynamic roughness,$$z_{0,{Eff.}} = \alpha \sigma _h$$, where$$\alpha$$is a proportionality constant and$$\sigma _h$$is height root mean square. We report$$\alpha \approx 10^{-1}$$, the value reported over many decades for a broad range of rough surfaces, for$$k$$-type cases at small$$\delta _2$$, whereas the transition to$$d$$-type arrangements necessitates larger$$\delta _2$$. Though preliminary, results highlight the non-trivial response to variation of streamwise and spanwise spacing.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2024.223
