Hi,

Now, I have a problem about the user subroutine in defining the density of the compressible fluid.

My goal is to define the background material as incompressible air and then specie as compressible fluid, but I can not find the variable to define the density for specie. If I use den(i), it will affect the whole domain including the specie itself.

Does anyones know how to define the density of specie in densit.f?

Uups... I you have success in defining a compressible fluid in a VOF simulation when using COMET, let me know. Because the VOF-algorithm for interface capturing works only if both fluids are incompressible. Unless you got one of the latest code developments...