I am working on heat pipe two phase simulation.In this density variation with temperature is an important parameter. I do not encounter any problem in my simulation if i opt constant vapor density, however when i go for ideal gas or polynomial (density variation) i am getting floating point exception error.. How to avoid this problem without changing density variation(ideal gas).