skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Viscoplastic boundary layers
In the limit of a large yield stress, or equivalently at the initiation of motion, viscoplastic flows can develop narrow boundary layers that provide either surfaces of failure between rigid plugs, the lubrication between a plugged flow and a wall or buffers for regions of predominantly plastic deformation. Oldroyd ( Proc. Camb. Phil. Soc. , vol. 43, 1947, pp. 383–395) presented the first theoretical discussion of these viscoplastic boundary layers, offering an asymptotic reduction of the governing equations and a discussion of some model flow problems. However, the complicated nonlinear form of Oldroyd’s boundary-layer equations has evidently precluded further discussion of them. In the current paper, we revisit Oldroyd’s viscoplastic boundary-layer analysis and his canonical examples of a jet-like intrusion and flow past a thin plate. We also consider flow down channels with either sudden expansions or wavy walls. In all these examples, we verify that viscoplastic boundary layers form as envisioned by Oldroyd. For each example, we extract the dependence of the boundary-layer thickness and flow profiles on the dimensionless yield-stress parameter (Bingham number). We find that, while Oldroyd’s boundary-layer theory applies to free viscoplastic shear layers, it does not apply when the boundary layer is adjacent to a wall, as has been observed previously for two-dimensional flow around circular obstructions. Instead, the boundary-layer thickness scales in a different fashion with the Bingham number, as suggested by classical solutions for plane-parallel flows, lubrication theory and, for flow around a plate, by Piau ( J. Non-Newtonian Fluid Mech. , vol. 102, 2002, pp. 193–218); we rationalize this second scaling and provide an alternative boundary-layer theory.  more » « less
Award ID(s):
1554044 2050396
PAR ID:
10427402
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
813
ISSN:
0022-1120
Page Range / eLocation ID:
929 to 954
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Journal of Fluid Mechanics)
    null (Ed.)
    Prior mathematical work of Constantin & Iyer ( Commun. Pure Appl. Maths , vol. 61, 2008, pp. 330–345; Ann. Appl. Probab. , vol. 21, 2011, pp. 1466–1492) has shown that incompressible Navier–Stokes solutions possess infinitely many stochastic Lagrangian conservation laws for vorticity, backward in time, which generalize the invariants of Cauchy ( Sciences mathématiques et physique , vol. I, 1815, pp. 33–73) for smooth Euler solutions. We reformulate this theory for the case of wall-bounded flows by appealing to the Kuz'min ( Phys. Lett. A , vol. 96, 1983, pp. 88–90)–Oseledets ( Russ. Math. Surv. , vol. 44, 1989, p. 210) representation of Navier–Stokes dynamics, in terms of the vortex-momentum density associated to a continuous distribution of infinitesimal vortex rings. The Constantin–Iyer theory provides an exact representation for vorticity at any interior point as an average over stochastic vorticity contributions transported from the wall. We point out relations of this Lagrangian formulation with the Eulerian theory of Lighthill (Boundary layer theory. In Laminar Boundary Layers (ed. L. Rosenhead), 1963, pp. 46–113)–Morton ( Geophys. Astrophys. Fluid Dyn. , vol. 28, 1984, pp. 277–308) for vorticity generation at solid walls, and also with a statistical result of Taylor ( Proc. R. Soc. Lond. A , vol. 135, 1932, pp. 685–702)–Huggins ( J. Low Temp. Phys. , vol. 96, 1994, pp. 317–346), which connects dissipative drag with organized cross-stream motion of vorticity and which is closely analogous to the ‘Josephson–Anderson relation’ for quantum superfluids. We elaborate a Monte Carlo numerical Lagrangian scheme to calculate the stochastic Cauchy invariants and their statistics, given the Eulerian space–time velocity field. The method is validated using an online database of a turbulent channel-flow simulation (Graham et al. , J. Turbul. , vol. 17, 2016, pp. 181–215), where conservation of the mean Cauchy invariant is verified for two selected buffer-layer events corresponding to an ‘ejection’ and a ‘sweep’. The variances of the stochastic Cauchy invariants grow exponentially backward in time, however, revealing Lagrangian chaos of the stochastic trajectories undergoing both fluid advection and viscous diffusion. 
    more » « less
  2. ; ; (, Journal of Fluid Mechanics)
    null (Ed.)
    We use an online database of a turbulent channel-flow simulation at $$Re_\tau =1000$$ (Graham et al. J. Turbul. , vol. 17, issue 2, 2016, pp. 181–215) to determine the origin of vorticity in the near-wall buffer layer. Following an experimental study of Sheng et al. ( J. Fluid Mech. , vol. 633, 2009, pp.17–60), we identify typical ‘ejection’ and ‘sweep’ events in the buffer layer by local minima/maxima of the wall stress. In contrast to their conjecture, however, we find that vortex lifting from the wall is not a discrete event requiring $$\sim$$ 1 viscous time and $$\sim$$ 10 wall units, but is instead a distributed process over a space–time region at least $$1\sim 2$$ orders of magnitude larger in extent. To reach this conclusion, we exploit a rigorous mathematical theory of vorticity dynamics for Navier–Stokes solutions, in terms of stochastic Lagrangian flows and stochastic Cauchy invariants, conserved on average backward in time. This theory yields exact expressions for vorticity inside the flow domain in terms of vorticity at the wall, as transported by viscous diffusion and by nonlinear advection, stretching and rotation. We show that Lagrangian chaos observed in the buffer layer can be reconciled with saturated vorticity magnitude by ‘virtual reconnection’: although the Eulerian vorticity field in the viscous sublayer has a single sign of spanwise component, opposite signs of Lagrangian vorticity evolve by rotation and cancel by viscous destruction. Our analysis reveals many unifying features of classical fluids and quantum superfluids. We argue that ‘bundles’ of quantized vortices in superfluid turbulence will also exhibit stochastic Lagrangian dynamics and satisfy stochastic conservation laws resulting from particle relabelling symmetry. 
    more » « less
  3. ; ; ; ; (, Journal of Fluid Mechanics)
    We present numerical simulations of laminar and turbulent channel flow of an elastoviscoplastic fluid. The non-Newtonian flow is simulated by solving the full incompressible Navier–Stokes equations coupled with the evolution equation for the elastoviscoplastic stress tensor. The laminar simulations are carried out for a wide range of Reynolds numbers, Bingham numbers and ratios of the fluid and total viscosity, while the turbulent flow simulations are performed at a fixed bulk Reynolds number equal to 2800 and weak elasticity. We show that in the laminar flow regime the friction factor increases monotonically with the Bingham number (yield stress) and decreases with the viscosity ratio, while in the turbulent regime the friction factor is almost independent of the viscosity ratio and decreases with the Bingham number, until the flow eventually returns to a fully laminar condition for large enough yield stresses. Three main regimes are found in the turbulent case, depending on the Bingham number: for low values, the friction Reynolds number and the turbulent flow statistics only slightly differ from those of a Newtonian fluid; for intermediate values of the Bingham number, the fluctuations increase and the inertial equilibrium range is lost. Finally, for higher values the flow completely laminarizes. These different behaviours are associated with a progressive increases of the volume where the fluid is not yielded, growing from the centreline towards the walls as the Bingham number increases. The unyielded region interacts with the near-wall structures, forming preferentially above the high-speed streaks. In particular, the near-wall streaks and the associated quasi-streamwise vortices are strongly enhanced in an highly elastoviscoplastic fluid and the flow becomes more correlated in the streamwise direction. 
    more » « less
  4. ; ; ; ; ; ; ; ; (, Proceedings of the National Academy of Sciences)
    Debris flows are dense and fast-moving complex suspensions of soil and water that threaten lives and infrastructure. Assessing the hazard potential of debris flows requires predicting yield and flow behavior. Reported measurements of rheology for debris flow slurries are highly variable and sometimes contradictory due to heterogeneity in particle composition and volume fraction ( ϕ ) and also inconsistent measurement methods. Here we examine the composition and flow behavior of source materials that formed the postwildfire debris flows in Montecito, CA, in 2018, for a wide range of ϕ that encapsulates debris flow formation by overland flow. We find that shear viscosity and yield stress are controlled by the distance from jamming, Δ ϕ = ϕ m − ϕ , where the jamming fraction ϕ m is a material parameter that depends on grain size polydispersity and friction. By rescaling shear and viscous stresses to account for these effects, the data collapse onto a simple nondimensional flow curve indicative of a Bingham plastic (viscoplastic) fluid. Given the highly nonlinear dependence of rheology on Δ ϕ , our findings suggest that determining the jamming fraction for natural materials will significantly improve flow models for geophysical suspensions such as hyperconcentrated flows and debris flows. 
    more » « less
  5. ; ; ; (, Journal of Fluid Mechanics)
    We consider the time evolution in two spatial dimensions of a double vorticity layer consisting of two contiguous, infinite material fluid strips, each with uniform but generally differing vorticity, embedded in an otherwise infinite, irrotational, inviscid incompressible fluid. The potential application is to the wake dynamics formed by two boundary layers separating from a splitter plate. A thin-layer approximation is constructed where each layer thickness, measured normal to the common centre curve, is small in comparison with the local radius of curvature of the centre curve. The three-curve equations of contour dynamics that fully describe the double-layer dynamics are expanded in the small thickness parameter. At leading order, closed nonlinear initial-value evolution equations are obtained that describe the motion of the centre curve together with the time and spatial variation of each layer thickness. In the special case where the layer vorticities are equal, these equations reduce to the single-layer equation of Moore ( Stud. Appl. Math. , vol. 58, 1978, pp. 119–140). Analysis of the linear stability of the first-order equations to small-amplitude perturbations shows Kelvin–Helmholtz instability when the far-field fluid velocities on either side of the double layer are unequal. Equal velocities define a circulation-free double vorticity layer, for which solution of the initial-value problem using the Laplace transform reveals a double pole in transform space leading to linear algebraic growth in general, but there is a class of interesting initial conditions with no linear growth. This is shown to agree with the long-wavelength limit of the full linearized, three-curve stability equations. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email