compressible flow computation
hello everybody
From the (legendry) book of Suhas Patankar we know the formulation of the flow computation for incompressible flow. I have myself written codes for constant density and variable density incompressible flows but can somebody tell me how to compute "COMPRESSIBLE" flow ? My flow is subsonic so much consideration abt transonic etc is not required. Also can somebody tell me how FLUENT computes COMPRESSIBLE flow ? The only ref I could get is "PhD thesis of KC Karki" which was written under guidance of Suhas Patankar only. THANKING YOU. 
Re: compressible flow computation
There are 2 issues involved to extend to compressible flow: 1. density treatment in continuity equation (density change in other transport equations is the same as variabledensity incompressible flow) and 2. boundary conditions.
When deriving pressure correction equation, you have u_new=u_old+u', and density_new=density_old+density' so mass flux = density_new*u_new*area = mass_flux_old + density_old*u'*area + u_old*density'*area + density'*u'*area The last term can be neglected, and compared to incompressible flow, you have extra term and you can link density' to p' by ideal gas law. Boundary conditions: for example subsonic inlet  you have to have 2 types of inlets: fixed mass and fixed velocities. The pressure at inlet is extrapolated from interior, and since temperature is fixed, you have changing density. If mass fixed, then normal velocity to the inlet surface has to change. If velocity fixed, then mass flux has to change. Just the idea. 
Re: compressible flow computation
Thank you Mr John. Actually with the above concepts I derived the pressure correction equation but it does not satisfy the condition aP= aE + aW + aN + aS . this is because we have rho_e, rho_n etc like terms which bring P'_e , P'_n etc like terms and no matter whether one uses central differencing or upwinding for density , one gets an equation where aP is slightly different from aE + aW + aN + aS . Am I write ? Also in FLUENT , they use some " density equation" but I wonder how do they derive it ?

Re: compressible flow computation
Abhijeet,
I guess the following paper may be helpful: K.C. Karki, S.V. Patankar, "Pressure Based Calculation Procedure for Viscous Flows at All Speeds in Arbitrary Configurations", AIAA J v27 n9, 1989. > Also in FLUENT , they use some " density equation" but I wonder how do they derive it ? I am not familiar with Fluent, but it may refer to the other possibility of solving compressible flows (usually with M>0.3). In this method, in contrast to the pressurebased methods, density (rather than pressure) is one of the solved variables  it is the conserved variable in the continuity equation. This method is straightforward, but difficulties are encountered for incompressible flows, where density is constant, and it is not linked to the pressure via an equation of state. The same applies (to a lesser degree) for lowMach flows as well. This was the motivation for devising the pressurebased methods. 
Thanks to Rami & John

Re: compressible flow computation
AP is NOT equal to AE+AW+AN+AS in compressible flow formula. Indeed using your structured mesh concept, AP is equal to AE_of_WEST_CELL + AW_of_EAST_CELL + AN_of_SOUTH_CELL + AS_of_NORTH_CELL and the matrix is no longer symmetry.

All times are GMT 4. The time now is 13:24. 