[mmoja]

I have a simple C code
real temp= C_T(c,t);
real press =C_P(c,t);
rho = D*pow(press,B)/pow(temp,Q);
for calculating density real gas.D,B and Q are all constants. if i hook up and iterate i get floating point error and when i add some constants value to temp and press it calculates and returns the density of the constant values and ignores the C_T and C_P... What should I do?

[coglione]

Hello mmoja,

this is just a presumption of what may went wrong, but C_P returns the static gauge pressure which may be a negative value depending on your operating pressure and boundary settings! Thus if constant B is not an even number pow(press,B) may become negative as well as your calculated density. A negative density is definitely a problem for any FV-code.

cheers

[mmoja]

Thanks for the reply, the constant B is 1.01 which is not an even number and my operating pressure is between 1-100bars...not atmospheric. Concerning the boundary conditions I have the gauge pressure as zero at the outlet pressure.... and mass flow for inlet boundary condition...