Abstract:In this work, we implement the local <math>μ(I)</math> rheology law and its non-local modification in a PISO-VOF numerical scheme with a sharp-interface-capturing THINC method for two-dimensional dense granular flows in both steady and transient states. The non-local effect on the constitutive relation is addressed by the gradient expansion model to account for additional momentum transport in view of the spatial variation of flow inertial number and the stress model singularity near the static state, <math>I=0</math>, is handled by the Bercovier-Engelman regularization scheme. Our simulations provide the first numerical evidence for the crucial role of non-local momentum transport on degrading the velocity profiles predicted with a local rheology model as those reported from two-dimensional discrete element simulations. For example, we successfully reproduce the sub-Bagnold velocity profile for steady inclined flows of height comparable to the <math>Hstop</math> and capture how the linear velocity profile in a simple shear creeping flow is degraded to a S-shape. More importantly, we exploit the current solver to study the layer-accumulated deviation, ε, between the local and the non-local velocity profile (Bagnold to its degradation) for inclined flows. For steady flows, ε decays monotonically with flow Froude number Fr and a critical <math>Frc≈</math>0.25 is detected below which ε rises rapidly necessitating a non-local constitutive relation. For a layer developing from rest to its steady state, the time evolution of <math>ε(t)</math> shows further dependence on flow height and can be collapsed into a monotonically growing trend with local instantaneous inertial number <math>I(t)</math>. A critical <math>Icr</math> is also determined and the local model prediction can be erroneous whenever <math>I(t)</math>. Interestingly, the result for simple shear flows without gravity also confirms that non-local effect becomes non-negligible when the macroscopic inertial number drops to the order of <math>Icr</math>. Finally, we simulate the two-dimensional column collapse flows with both local and non-local rheology and find a shorter run-out distance with the presence of non-local momentum transport.

Continuum Model for Nanoscale Multiphase Flows

Numerical Investigation on Local Non-equilibrium Flows Using a Diatomic Nonlinear Constitutive Model

Numerical Analysis of the Micro-Couette Flow Using a Non-Newton–fourier Model with Enhanced Wall Boundary Conditions

Multi-Scale Gas-Liquid Interfacial Model

A bulk-surface continuum theory for fluid flows and phase segregation with finite surface thickness

Two-phase flows with bulk-surface interaction: thermodynamically consistent Navier-Stokes-Cahn-Hilliard models with dynamic boundary conditions

Nonequilibrium model for compressible two-phase two-pressure flows with surface tension

Continuum description of rarefied gas dynamics. I. Derivation from kinetic theory

A unified continuum and variational multiscale formulation for fluids, solids, and fluid-structure interaction

Derivation of a two-phase flow model with two-scale kinematics, geometric variables and surface tension using variational calculus

Continuum models for the contact line problem

Finite element methods for a class of continuum models for immiscible flows with moving contact lines

A cut-cell method for the numerical simulation of 3D multiphase flows with strong interfacial effects

Upscaling a Navier-Stokes-Cahn-Hilliard model for two-phase porous-media flow with solute-dependent surface tension effects

Continuum simulation for regularized non-local μ(I) model of dense granular flows

Mesoscale modeling of near-contact interactions for complex flowing interfaces

An implicit surface tension model for the analysis of droplet dynamics

A Contoured Continuum Surface Force Model for Particle Methods

Multicontinuum Homogenization for Coupled Flow and Transport Equations

A unified two-scale gas-liquid multi-fluid model with capillarity and interface regularization through a mass transfer between scales