Inlet Velocity Boundary Condition for Shan-Chen Model

Lattice Boltzmann Method LBM: Theory
#1

Dear Friends,

I have a question on implementing inlet velocity BC in the Shan-Chen multiphase model. For calculating the interaction force, the effective mass (e.g. psi(x)) for the neighbor nodes (e.g. psi(x+e) & psi(x-e)) should be known. In the case of inlet velocity, the function psi(x-e) is out of the computational domain. How can I compute the interaction force for inlet boundary nodes? I really appreciate your kind suggestions and helps.

Regards,
Hesam

1 Like
#2

Dear all,

I’ve found a hint on this problem in Sukop’s books. I implemented it in my code and it seems work.

Thanks,
Hesam

1 Like
#3

Hi Hesam,

Can you please mention the page numbers and lines in which you found the hints? It would be of great help for me as I am working on similar thing.

Thanks in advance,
Narender Koosukuntla

#4

You can define a constant psi_inlet according to the phase inlet boundary. That means you define an inlet rho (or pressure).

Two papers,

PHYSICAL REVIEW E 71, 046703 (2005)
Computers & Mathematics with Applications, 61(12) 3577-3584, 2011

might be helpful.

#5

Dear Narender,

You can find it in the first paragraph of page 80.

Regards,
Hesam

1 Like
#6

Hi Hesam,

Thanks a lot for the information.

Best
Narender