A granular flow is normally comprised of a mixture of grain-particles (such as sand, gravel or rocks) of different sizes. In this study, dry granular flows are modeled utilizing a set of equations akin to a two-phase mixture system, in which the interstitial fluid is air. The resultant system of equations for a two-dimensional configuration includes two continuity and two momentum balance equations for the two respective constituents. The density variation is described considering the phenomenon of air entrainment/extrusion at the flow surface, where the entrainment rate is assumed to be dependent on the divergent or convergent behavior of the solid constituent. The density difference between the two constituents is extremely large, so, as a consequence scaling analysis reveals that the flow behavior is dominated by the solid species, yielding small relative velocities between the two constituents. A non-oscillatory central (NOC) scheme with total variation diminishing (TVD) limiters is implemented. Three numerical examples are investigated: the first being related to the flow behaviors on a horizontal plane with an unstable initial condition; the second example is devoted to simulating a dam-break problem with respect to different initial conditions; and in the third one investigates the behavior of a finite mass of granular material flowing down an inclined plane. The key features and the capability of the equations to model the behavior are illustrated in these numerical examples. (C) 2012 The Society of Powder Technology Japan. Published by Elsevier B.V. and The Society of Powder Technology Japan. All rights reserved.