Reply To: Implement Gravity in Shen Chen model

[zshi6193]

Dear Max,

Thanks for your reply. Much appreciated.
I am sorry that I have one more question regarding the pressure boundary used in Shan-Chen two compartment multiphase model. To set the pressure inlet and pressure outlet, I used the following methods:
slatticeWater.defineRho(supergeometry,3(inlet),rhowater)
slatticeOil.defineRho(supergeometry,3(inlet),zero)
slatticeWater.defineRho(supergeometry,4(outlet),rhowater)
slatticeOil.defineRho(supergeometry,4(outlet),zero).
(I applied above four codes every iteraction to keep the constant pressure)
Alternatively, I also try:
Bcwater.addPressureBoundary(supergeometry,3(inlet),omega)
BcWater.addPressureBoundary(supergeometry,4(outlet),omega)
slatticeOil.defineRho(supergeometry,3(inlet),zero).
slatticeOil.defineRho(supergeometry,4(outlet),zero).

The middle part of the domain is filled with oil and the rest is filled with water. But under above set up, the equilibrium state can not be achieved. The pre-filled oil layer moves downwards.
I am not sure that I set up the boundary conditions correctly or not? .Could you help me to check it? Ideally, at the same pressure at the top and bottom (no pressure gradient), the oil in the middle part should be stable.

Thank you very much.
Simon