Granular flows occur in a wide range of situations of practical interest to industry, in our natural environment and in our everyday lives. This paper focuses on granular flow in the so-called inertial regime, when the rheology is independent of the very large particle stiffness. Such flows have been modelled with the $\unicode[STIX]{x1D707}(I),\unicode[STIX]{x1D6F7}(I)$ -rheology, which postulates that the bulk friction coefficient $\unicode[STIX]{x1D707}$ (i.e. the ratio of the shear stress to the pressure) and the solids volume fraction $\unicode[STIX]{x1D719}$ are functions of the inertial number $I$ only. Although the $\unicode[STIX]{x1D707}(I),\unicode[STIX]{x1D6F7}(I)$ -rheology has been validated in steady state against both experiments and discrete particle simulations in several different geometries, it has recently been shown that this theory is mathematically ill-posed in time-dependent problems. As a direct result, computations using this rheology may blow up exponentially, with a growth rate that tends to infinity as the discretization length tends to zero, as explicitly demonstrated in this paper for the first time. Such catastrophic instability due to ill-posedness is a common issue when developing new mathematical models and implies that either some important physics is missing or the model has not been properly formulated. In this paper an alternative to the $\unicode[STIX]{x1D707}(I),\unicode[STIX]{x1D6F7}(I)$ -rheology that does not suffer from such defects is proposed. In the framework of compressible $I$ -dependent rheology (CIDR), new constitutive laws for the inertial regime are introduced; these match the well-established $\unicode[STIX]{x1D707}(I)$ and $\unicode[STIX]{x1D6F7}(I)$ relations in the steady-state limit and at the same time are well-posed for all deformations and all packing densities. Time-dependent numerical solutions of the resultant equations are performed to demonstrate that the new inertial CIDR model leads to numerical convergence towards physically realistic solutions that are supported by discrete element method simulations.
more »
« less
Rotation of a superhydrophobic cylinder in a viscous liquid
The hydrodynamic quantification of superhydrophobic slipperiness has traditionally employed two canonical problems – namely, shear flow about a single surface and pressure-driven channel flow. We here advocate the use of a new class of canonical problems, defined by the motion of a superhydrophobic particle through an otherwise quiescent liquid. In these problems the superhydrophobic effect is naturally measured by the enhancement of the Stokes mobility relative to the corresponding mobility of a homogeneous particle. We focus upon what may be the simplest problem in that class – the rotation of an infinite circular cylinder whose boundary is periodically decorated by a finite number of infinite grooves – with the goal of calculating the rotational mobility (velocity-to-torque ratio). The associated two-dimensional flow problem is defined by two geometric parameters – namely, the number $N$ of grooves and the solid fraction $\unicode[STIX]{x1D719}$ . Using matched asymptotic expansions we analyse the large- $N$ limit, seeking the mobility enhancement from the respective homogeneous-cylinder mobility value. We thus find the two-term approximation, $$\begin{eqnarray}\displaystyle 1+{\displaystyle \frac{2}{N}}\ln \csc {\displaystyle \frac{\unicode[STIX]{x03C0}\unicode[STIX]{x1D719}}{2}}, & & \displaystyle \nonumber\end{eqnarray}$$ for the ratio of the enhanced mobility to the homogeneous-cylinder mobility. Making use of conformal-mapping techniques and inductive arguments we prove that the preceding approximation is actually exact for $N=1,2,4,8,\ldots$ . We conjecture that it is exact for all $N$ .
more »
« less
- Award ID(s):
- 1909407
- NSF-PAR ID:
- 10173301
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 880
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We examine correlations of the Möbius function over $\mathbb{F}_{q}[t]$ with linear or quadratic phases, that is, averages of the form 1 $$\begin{eqnarray}\frac{1}{q^{n}}\mathop{\sum }_{\deg fmore » « less
-
Let $f\in C^{2}(\mathbb{T}^{2})$ have mean value 0 and consider $$\begin{eqnarray}\sup _{\unicode[STIX]{x1D6FE}\,\text{closed geodesic}}\frac{1}{|\unicode[STIX]{x1D6FE}|}\biggl|\int _{\unicode[STIX]{x1D6FE}}f\,d{\mathcal{H}}^{1}\biggr|,\end{eqnarray}$$ where $\unicode[STIX]{x1D6FE}$ ranges over all closed geodesics $\unicode[STIX]{x1D6FE}:\mathbb{S}^{1}\rightarrow \mathbb{T}^{2}$ and $|\unicode[STIX]{x1D6FE}|$ denotes its length. We prove that this supremum is always attained. Moreover, we can bound the length of the geodesic $\unicode[STIX]{x1D6FE}$ attaining the supremum in terms of the smoothness of the function: for all $s\geq 2$ , $$\begin{eqnarray}|\unicode[STIX]{x1D6FE}|^{s}{\lesssim}_{s}\biggl(\max _{|\unicode[STIX]{x1D6FC}|=s}\Vert \unicode[STIX]{x2202}_{\unicode[STIX]{x1D6FC}}f\Vert _{L^{1}(\mathbb{T}^{2})}\biggr)\Vert \unicode[STIX]{x1D6FB}f\Vert _{L^{2}}\Vert f\Vert _{L^{2}}^{-2}.\end{eqnarray}$$more » « less
-
We give an integrability criterion on a real-valued non-increasing function $\unicode[STIX]{x1D713}$ guaranteeing that for almost all (or almost no) pairs $(A,\mathbf{b})$ , where $A$ is a real $m\times n$ matrix and $\mathbf{b}\in \mathbb{R}^{m}$ , the system $$\begin{eqnarray}\Vert A\mathbf{q}+\mathbf{b}-\mathbf{p}\Vert ^{m}<\unicode[STIX]{x1D713}(T),\quad \Vert \mathbf{q}\Vert ^{n}This paper presents a theory to obtain the force per unit length acting on a slender filament with a non-circular cross-section moving in a fluid at low Reynolds number. Using a regular perturbation of the inner solution, we show that the force per unit length has $O(1/\ln (2A))+O(\unicode[STIX]{x1D6FC}/\ln ^{2}(2A))$ contributions driven by the relative motion of the particle and the local fluid velocity and an $O(\unicode[STIX]{x1D6FC}/(\ln (2A)A))$ contribution driven by the gradient in the imposed fluid velocity. Here, the aspect ratio ( $A=l/a_{0}$ ) is defined as the ratio of the particle size ( $l$ ) to the cross-sectional dimension ( $a_{0}$ ) and $\unicode[STIX]{x1D6FC}$ is the amplitude of the non-circular perturbation. Using thought experiments, we show that two-lobed and three-lobed cross-sections affect the response to relative motion and velocity gradients, respectively. A two-dimensional Stokes flow calculation is used to extend the perturbation analysis to cross-sections that deviate significantly from a circle (i.e. $\unicode[STIX]{x1D6FC}\sim O(1)$ ). We demonstrate the ability of our method to accurately compute the resistance to translation and rotation of a slender triaxial ellipsoid. Furthermore, we illustrate novel dynamics of straight rods in a simple shear flow that translate and rotate quasi-periodically if they have two-lobed cross-section, and rotate chaotically and translate diffusively if they have a combination of two- and three-lobed cross-sections. Finally, we show the remarkable ability of our theory to accurately predict the motion of rings, retaining great accuracy for moderate aspect ratios ( ${\sim}10$ ) and cross-sections that deviate significantly from a circle, thereby making our theory a computationally inexpensive alternative to other Stokes flow solvers.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2019.776
Warning: Leaving National Science Foundation Website
You are now leaving the National Science Foundation website to go to a non-government website.
Website:
NSF takes no responsibility for and exercises no control over the views expressed or the accuracy of the information contained on this site. Also be aware that NSF's privacy policy does not apply to this site.
