divergent result

[lollipopcorn]

Hi,

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_test.cpp: http://text-share.com/view/275cf34f
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.

Thanks,
John

[mathias]

Unstabilities can have many sources, the boundary or initial conditions, the model some geometric source.. so it is hard to tell you in a forum. You need to start form a case where it is stable and start altering things step by step to find the source. Further, you may come to our spring school next year. Best Mathias

[Author]

Posts