You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
136<br />
A Simple Model for Gas-Grain Two Phase Isentropic Flow<br />
and a High Resolution Scheme for the Numerical Solution of<br />
the Governing Equations<br />
J. Hudson andD.Harris ∗<br />
This contribution presents a simple model for the isentropic two phase flow of a<br />
solid granular material dispersed in a gas obeying the perfect gas law. The dispersed<br />
and continuous phases are both treated as continua and an Eulerian description of<br />
the flow is adopted.<br />
The first key physical quantity in the model is the solids volume fraction and the<br />
model comprises continuity equations and balances of linear momentum for the gas<br />
and solid phases.<br />
The second key physical quantity is the fluctuation energy for the solids phase,<br />
otherwise known as the granular temperature. The fluctuation energy arises in the<br />
theory of the statistical mechanics of granular materials and dense gases and was<br />
later introduced into continuum models. The model is complete by the fluctuation<br />
energy equation. The presence of the solids fluctuation energy gives rise to a solids<br />
pressure in the balance of momentum for the solids phase. We assume the simplest<br />
possible equation of state for the solids pressure, namely the analogy of the perfect<br />
gas law. In the absence of viscosity the classical equal pressures model is ill-posed,<br />
but we demonstrate numerically that the presence of the solids pressure due to the<br />
solids fluctuation energy gives rise to regimes in which the model is indeed well-posed.<br />
The hyperbolicity of the system of equations is investigated numerically and the<br />
numerical solution of the equations in conservation form is accomplished using a<br />
high-resolution scheme 1 in the hyperbolic regime. Careful attention is given to the<br />
discretisation of the inhomogeneous terms. Three one-dimensional test cases are considered,<br />
namely a simple advection test problem, a square pulse test problem and a<br />
steady state problem and we use these to obtain quantitative and qualitative insight<br />
into the predictions of the model.<br />
∗ School of Mathematics, University of Manchester.<br />
1 Hudson et al, submitted to J. Comp. Phys.