Skip to content

divergent result

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

    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

    #2696
    mathias
    Keymaster

    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.