Compressible code problems
 

Hello all

I am trying to modify rhoCentralFoam to introduce godunov finite volume scheme. I have attached my initial attempt in this thread. Basically, I have modified the main code rhoCentralFoam.C now called gFoam.C and createFields.H

I have added new variables rhoFlux, rhoUFlux and rhoEFlux which contain the numerical flux across all the faces.

Also added two new files numericalFlux.H and KFVS.H which contain the flux functions.

Its not working as I am getting nan on the oblique shock case. It would be nice if somebody can look at the code and tell me if there are any mistakes.

I assumed that rPsi = R*T where R=gas constant and T is the temperature.

Thanks
praveen