I have some problems with boundary condition in Palabos. Only periodic and bounceback boundary condition works fine with the SCMP Model, but the others like velocity or outflow boundaries produce mistakes. It is possible to use a flow boundaries with this model?
You’re right: except for bounce-back and periodicity, the boundary conditions in Palabos are made for the “incompressible, single-component Navier-Stokes” models. This seems to be a delicate topic, not only because implementing an adequate boundary conditions in Palabos requires some work, but also because, to my knowledge, the literature on how to implement boundary conditions with the Shan/Chen model is very vague. Nevertheless, all of this is part of the “under development” section of the Palabos project.
In the meantime you might find that you will get a stable (if not very accurate) simulation by imposing an equilibrium distribution on the boundaries (in the two-component model, you need to do this for both fluids, and it is your responsibility to come up with a reasonable value for the density). This can be done with the “createEquilibriumBoundaryConditoinXD” boundary condition, or even more brutally, by manually imposing the equilibrium with a given velocity and density value after each iteration through a “initializeAtEquilibrium” command applied to the boundary nodes.
I am trying to run multiphase flow simulation in porous media. I am injecting one phase to displace another one.
I started with rayleighTaylor and then, I tried to change boundary conditions.
In my simulation, I tried using inlet and outlet boundary conditions as:
but this gives me Nan answer!
Could anybody help me with this?
Moreover, I have a question about rayleighTaylor problem.
For setup of rayleighTaylor problem, we define two regions with different densities.
What these densities really mean? Are they physical densities?