Dear Marcin,
I have two questions about multiphase simulation in Sailfish, which I hope you could kindly give me some guidance.
1. In sc_laplace_2d.py the LBBinaryFluidShanChen model is applied. The pressure is calculated in this formula:p1 = (phi1 + rho1) + G * phi1 * rho1. I suppose rho and phi is the density of the two different phase/ fluids, right? But why the pressure is calculated this way? Can you give me a hint of related papers? Thanks!
2. Can I apply Density or Velocity boundary condition under multi-phase models using the latest software?
I have checked all multi-phase examples using LBSingleFluidShanChen, LBBinaryFluidShanChen, LBBinaryFluidFreeEnergy, and find that the boundary conditions are either wall or periodic condition.
I have tried to use Density or Velocity boundary condition(not in
initial condition) under those multi-phase models as mentioned, but it
failed with the following: Logic Error: cuStreamSynchronize failed: an illegal instruction was encountered.
Is there any example which can successfully apply density or velocity condition under multi-phase models? If not, can you tell me which part of Sailfish needs adjustment? Thanks!