Skip to content

divergent result

Viewing 2 posts - 1 through 2 (of 2 total)
  • Author
  • #1932


    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:

    It would be greatly appreciated if anyone could give me some hint on where could possibly go wrong.



    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

Viewing 2 posts - 1 through 2 (of 2 total)
  • You must be logged in to reply to this topic.