I tried to apply what I learned from example cylinder2D, bstep2D, and multicomponent2D to create a test simulation where water is intruding into the air within a channel filled with circular obstacles.
I changed the density of air to 900 in order to decrease the density ratio RhoWater/RhoAir = 900/1000.
But the density of water becomes ‘nan’ during the running process.
My code is in the link below:
Water_Air dat file: http://text-share.com/view/62367840
It would be greatly appreciated if anyone could give me some hint on where could possibly go wrong.