No abstract

Abstract. We consider two physically and mathematically distinct regularization mechanisms of scalar hyperbolic conservation laws. When the flux is convex, the combination of diffusion and dispersion are known to give rise to monotonic and oscillatory traveling waves that approximate shock waves. The zero-diffusion limits of these traveling waves are dynamically expanding dispersive shock waves (DSWs). A richer set of wave solutions can be found when the flux is non-convex. This review compares the structure of solutions of Riemann problems for a conservation law with non-convex, cubic flux regularized by two different mechanisms: 1) dispersion in the modified Korteweg-de Vries (mKdV) equation; and 2) a combination of diffusion and dispersion in the mKdV-Burgers equation. In the first case, the possible dynamics involve two qualitatively different types of DSWs, rarefaction waves (RWs) and kinks (monotonic fronts). In the second case, in addition to RWs, there are traveling wave solutions approximating both classical (Lax) and non-classical (undercompressive) shock waves. Despite the singular nature of the zero-diffusion limit and rather differing analytical approaches employed in the descriptions of dispersive and diffusive-dispersive regularization, the resulting comparison of the two cases reveals a number of striking parallels. In contrast to the case of convex flux, the mKdVB to mKdV mapping is not one-to-one. The mKdV kink solution is identified as an undercompressive DSW. Other prominent features, such as shock-rarefactions, also find their purely dispersive counterparts involving special contact DSWs, which exhibit features analogous to contact discontinuities. This review describes an important link between two major areas of applied mathematics, hyperbolic conservation laws and nonlinear dispersive waves.

No abstract

No abstract

A partially-wetting liquid can deform the underlying elastic substrate upon which it rests. This situation requires the development of theoretical models to describe the wetting forces imparted by the drop onto the solid substrate, particularly those at the contact-line. We construct a general solution using a displacement potential function for the elastic deformations within a finite elastic substrate associated with these wetting forces, and compare the results for several different contact-line models. Our work incorporates internal contributions to the surface stress from both liquid/solid Σls and Σsg solid/gas solid surface tensions (surface stress), which results in a non-standard boundary-value problem that we solve using a dual integral equation. We compare our results to relevant experiments and conclude that the generalization of solid surface tension Σls ≠ Σsg is an essential feature in any model of partial-wetting. The comparisons also allow us to systematically eliminate some proposed contact-line models.

Continuum modelling of granular flow has been plagued with the issue of ill-posed dynamic equations for a long time. Equations for incompressible, two-dimensional flow based on the Coulomb friction law are ill-posed regardless of the deformation, whereas the rate-dependent μ(I)-rheology is ill-posed when the non-dimensional inertial number I is too high or too low. Here, incorporating ideas from critical-state soil mechanics, we derive conditions for well-posedness of partial differential equations that combine compressibility with I-dependent rheology. When the I-dependence comes from a specific friction coefficient μ(I), our results show that, with compressibility, the equations are well-posed for all deformation rates provided that μ(I) satisfies certain minimal, physically natural, inequalities.

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 µ(I), Φ(I)-rheology, which postulates that the bulk friction coefficient µ (i.e. the ratio of the shear stress to the pressure) and the solids volume fraction φ are functions of the inertial number I only. Although the µ(I), Φ(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 µ(I), Φ(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 µ(I) and Φ(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.

Granular materials segregate by size under shear, and the ability to quantitatively predict the time required to achieve complete segregation is a key test of our understanding of the segregation process. In this paper, we apply the Gray-Thornton model of segregation (developed for linear shear profiles) to a granular flow with an exponential profile, and evaluate its ability to describe the observed segregation dynamics. Our experiment is conducted in an annular Couette cell with a moving lower boundary. The granular material is initially prepared in an unstable configuration with a layer of small particles above a layer of large particles. Under shear, the sample mixes and then re-segregates so that the large particles are located in the top half of the system in the final state. During this segregation process, we measure the velocity profile and use the resulting exponential fit as input parameters to the model. To make a direct comparison between the continuum model and the observed segregation dynamics, we locally map the measured height of the experimental sample (which indicates the degree of segregation) to the local packing density. We observe that the model successfully captures the presence of a fast mixing process and relatively slower re-segregation process, but the model predicts a finite re-segregation time, while in the experiment re-segregation occurs only exponentially in time.