Dear all,
We are a grouo of students,we tried to implement a 2D channel multhiphase flow, we based on the matlab examples.After we stopped, there were two problems:

Once you have chosen to simulate liquids, which are the values attributed to the constant G (amplitude of the iteraction Strength) in ShanChen.m ?

regarding the initialization of the two fluids, we have:
fIn = reshape( tNS’ * zeros(1,lxly), 9, lx, ly);
gIn = reshape( tNS’ * ones(1,lxly), 9, lx, ly);
So Hypothesis is that the channel is filled with a fluid and the other is vacuum , but numerical simulations does uncorrect results.
how can we initialize these two fluids?
anyone have any ideas to solve these two questions?
Thanks for help