I am currently doing a project using multiphase ShanChen to simulate bubble detachment. One of the set-ups we are discussing is a Poiseuille inflow of liquid in a 2D cylinder and a smaller vapour inflow perpendicular to it and located in the middle of the cylinder. I am currently working on top of the cylinder2D example.
In summary, the current boundary conditions are: left inlet: Poiseuille flow; right outlet: Dirichlet pressure condition p=0; bottom inlet: Poiseuille flow.
However I am having trouble assigning a different inlet density (even using a small density ratio) for the bottom one using defineRhoU() and iniEquilibrium() with the associated material number. Is this the correct approach? How can I make it better?