# Computing Roe's average density for General Equation of State

I am solving a 1d Shock tube problem of compressible fluid obeying Euler equations(Hyperbolic pde). I am trying to simulate it using Finite Volume Method using Roe's Scheme. Half of the tube contains highly compressed fluid and other half is not compressed with density ratio of ~1000.

Euler Equation: $$∂Q/∂t + ∂F/∂x = 0; A = ∂F/∂Q$$

In Roe's Scheme, the parameters of velocity 'u' and Enthalpy 'H', at the interface are need for determining the 'A' a 3x3 matrix(for 1D Euler equation) . They are obtained by satisfying the following condition.

$$A (Q_\text{Left} - Q_\text{Right}) = \text{F}_\text{Left} - \text{F}_\text{Right}\ .$$

But the density, $\rho$ remains undetermined. While reading through literature, it was mentioned that it is 'Natural' to consider

$$\rho = \sqrt{\rho_\text{Left}.\rho_\text{Right}}\ .$$

My question is why is it so and what happens if we consider density to be a arithmetic mean(or some other mean) of the two density values.

I am trying to obtain Matrix 'A' for a different Equation of State(EoS) other than ideal gas EoS. If I choose a particular value of the density(calculated from density on either side of the face), the Roe's average values are greatly simplified.

EoS considered: $$Pressure, P(\rho, E) = (\gamma-1).\rho E+g(\rho)$$

I was obtaining the density from following condition for simplification of equations for 'u' and 'H'. $$dg(\rho)/d\rho = (g(\rho_\text{Left})-g(\rho_\text{Right}))/(\rho_\text{Left}-\rho_\text{Right})$$

The reason for my question is enquire about the mathematical implications in choosing a density, other than 'Natural' square root value.

• I assume that EoS means Equation of State, but it might not be obvious for everybody. – nicoguaro Aug 4 '17 at 13:41
• Thank you @ nicoguaro for point it out this. By EoS, I meant Equation of State. I made the necessary edit in the question – Rakesh Srr Aug 4 '17 at 14:37
• Welcome to Scicomp.SE! I have to agree with nicoguaro, your question as it is currently formulated is near impossible to understand (I literally have no idea what you are asking). Can you add more details such as what state are you actually modeling, how the density enters into the equation, and what your goal here is? (Note: I'm not voting to close for now, because I believe we should give some time -- say, 24 hours -- for simple clarifications, as well as some guidance, before calling for a time-out.) – Christian Clason Aug 4 '17 at 14:39
• I have edited the question to include more information – Rakesh Srr Aug 4 '17 at 15:33
• As far as I know, any matrix is suitable as long as it fulfills the conditions to be a Roe matrix. The impact of the choice of an average state is minimal, so you should feel free chose the most convenient. I'll see if I can find some references to back this before I post it as an answer. – Matthieu Aug 8 '17 at 19:00

• Thank you Ramanathan. I have revisited the paper you mention. My current understanding is apart from the following equation:$$A (Q_\text{Left} - Q_\text{Right}) = \text{F}_\text{Left} - \text{F}_\text{Right}\ .$$. It has to also satify: $\Delta$($\rho$u) = $\Delta$($\rho).$u+$\rho$ $\Delta.$(u).. Let me know if my understanding is right. – Rakesh Srr Aug 15 '17 at 16:47