Skip to content

Why am i keep getting inf velocity when using shan-chen model?

OpenLB – Open Source Lattice Boltzmann Code Forums on OpenLB General Topics Why am i keep getting inf velocity when using shan-chen model?

Viewing 2 posts - 1 through 2 (of 2 total)
  • Author
    Posts
  • #4698
    nv4dll
    Participant

    Dear all!
    I’m using shan-chen two-component model for simulation, and I was keep getting infinite lattice velocity during the simulation, and I tried to change maxvelocity on the Boundary and it doesn’t work.
    The output is like this :
    [LatticeStatistics] step=90; t=0.09; uMax=inf; avEnergy=-nan; avRho=-nan
    Can anyone help me? I’m kind of desperate. I want to late one heavier fluid driven another lighter fluid out of a tube.

    Here is some of my code:
    lattice1.iniEquilibrium( superGeometry.getMaterialIndicator({1,3,4}), rhoF1,uF );
    lattice1.defineRhoU( superGeometry.getMaterialIndicator({1,3,4}), rhoF1,uF );
    Do I need to initialize the second fluid in the Geometry? Since the second fluid isn’t in the tube at the begining.

    I checked the results and it seem that there is no fluid out of the tube.
    I will be VERY APPRECIATE if someone can help me out!
    Thanks for the feedback and advices

    #4699
    mathias
    Keymaster

    Dear nv4dll,

    It seems that your simulation has diverged – it is not stable. You should always start with an existing set-up and change step by step small settings and check after each step the results. I also recomment you to come to our next spring school in March in April. We have a section on multi-component/phase flows.

    Best
    Mathias

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