Skip to content

Additional force term

Viewing 3 posts - 1 through 3 (of 3 total)
  • Author
    Posts
  • #1762
    Wolf
    Member

    Hello everyone!rnI have a simple 3D Poiseuille flow with pressure inflow and outflow boundary conditions. I need to impose an additional force term to a set of the cells inside the domain. The term is simply added to the post-collision distributions: f* = f* + a, where a is a force term. What would be the best way to implement this with OpenLB? Thank you!

    #2150
    mathias
    Keymaster

    Hi Wolf,rnrnThere is already a model to include external forces inplemeted. Please look at the “”forcedPoiseuille”” example. You can also add any term to f*. Just like rn lattice.get(iX,iY)[iPop] += a;rnrnBut mind that in OpenLB f = f(standard)-f^eq(1,0)!rnrnBest regardsrnMathias

    #2151
    Wolf
    Member

    Hi Mathias,rnThanks a lot! That was really helpful!rnrnBest Regards,rnWolf

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