We propose and analyse finite volume Godunov type methods based on discontinuous flux for a 2×2 system of non-linear partial differential equations proposed by Hadeler and Kuttler to model the dynamics of growing sandpiles generated by a vertical source on a flat bounded rectangular table. The problem considered here is the so-called partially open table problem where sand is blocked by a wall (of infinite height) on some part of the boundary of the table. The novelty here is the corresponding modification of boundary conditions for the standing and the rolling layers and generalization of the techniques of the well-balancedness proposed in [1]. Presence of walls may lead to unbounded or discontinuous surface flow density at equilibrium resulting in solutions with singularities propagating from the extreme points of the walls. A scheme has been proposed to approximate efficiently the Hamiltonians with the coefficients which can be unbounded and discontinuous. Numerical experiments are presented to illustrate that the proposed schemes detect these singularities in the equilibrium solutions efficiently and comparisons are made with the previously studied finite difference and Semi-Lagrangian approaches by Finzi Vita et al.

We consider a dense granular shear flow in a two-dimensional system. Granular systems(composed of a large number of macroscopic particles) are far from equilibrium due toinelastic collisions between particles: an external driving is needed to maintain themotion of particles. Theoretical description of driven granular media is especiallychallenging for dense granular flows. This paper focuses on a gravity-driven densegranular Poiseuille flow in a channel. A special focus here is on the intriguingphenomenon of fluid-solid coexistence: a solid plug in the center of the system,surrounded by fluid layers. To find and analyze various flow regimes, a multi-scaleapproach is taken. On macro scale, granular hydrodynamics is employed. On micro scale,event-driven molecular dynamics simulations are performed. The entire phase diagram ofparameters is explored, in order to determine which flow regime occurs in various regionsin the parameter space.

An evolution equation for the Cauchy stress tensor is proposed, taking into account the elastic, viscous and plastic characteristics of complex fluids. Hypoplasticity, in particular, is incorporated to model the plastic features. The model is applied to investigate time-dependent Poiseuille flows between two parallel plates to simulate non-Newtonian behavior of complex rheological fluids. Results show that while different degrees of elastic and viscous features can be indexed by varying the values of the model parameters, hypoplasticity is capable of simulating the plastic characteristics. The fluid tends to be divided into different parallel layers as hypoplastic effects enhance gradually. Applications of the model may be found in geomorphic fluid motions involved granular materials.