Granular flow velocity model

Due to the hyperbolicity and nonlinearity of the model equations, associated with possible shocks in granular flows over non-trivial topography, numerical solutions with the traditional high-order accuracy methods are often accompanied with numerical oscillations of the depth profile and velocity field. This usually leads to numerical instabilities unless these are properly counteracted by a sufficient amount of artificial numerical diffusion. Here, a non-oscillatory central (NOC) difference scheme with a total variation diminishing (TVD) limiter for the cell reconstruction is employed, see e.g. [4], [12] we obtain numerical solutions without spurious oscillations. In order to test the model equations, we consider an ideal mountain subregion in which the talweg is defined by the slope function... 

The theoretical studies of rapid granular flows are generally based on the assumption that the energy dissipation in a binary particle collision is determined by a constant coefficient of restitution e, the ratio of the relative approach to recoil velocities normal to the point of impact on the particle. However, measurements show that the coefficient of restitution is a strong function of the relative impact velocity [10]. Physically, the energy dissipation relates to the plastic deformation of the particle s surface. Thus, a realistic microscopic model should include the deformation history of the particle s surface. However, such a model might become computationally demanding and thus not feasible. 

To model the particle velocity fluctuation covariances caused by particle-particle collisions and particle interactions with the interstitial gas phase, the concept of kinetic theory of granular flows is adapted (see chap 4). This theory is based on an analogy between the particles and the molecules of dense gases. The particulate phase is thus represented as a population of identical, smooth and inelastic spheres. In order to predict the form of the transport equations for a granular material the classical framework from the kinetic theory of... 

Several early and recent investigators have derived variations of the throughput, axial velocity, and residence time expressions (Hogg et al., 1973 Perron and Bui, 1990) but Seaman s expressions given here have found most practical applications and are recommended for use in industrial kilns so long as there are no true granular flow models to predict the bed behavior as are available for fluids of isotropic materials. 

In the literature numerous two-fluid models of different complexity have been proposed to predict the fluidized bed reactor cold flow and reactive flow behaviors. Four decades ago emphasis was placed on the modeling of the velocity fluctuation co-variance terms in the dispersed particle fluid phase momentum equations. The early one-dimensional models were normally closed by an elasticity modulus parameterization for the particle phase collisional pressure and a constant viscosity parameter for the corresponding shear stresses. Later, with the improved computer memory and speed capacities, multi-dimensional flow models and more advanced model closures were developed based on the kinetic theory of granular flow (KTGF). Moreover, the... 

---