Reply To: ContactAngle2d.cpp
The issue is that the densities of the two SuperLattices are not the densities of the fluid components. The density of sLattice1 is the total fluid density (= rho_1 + rho_2), while the “density” of sLattice2 is an order parameter to define whether it is fluid 1 or fluid 2 (= rho_1 – rho_2).
Therefore to plot the total density you only need to plot “density1”. It should be slightly higher inside the droplet due to the Laplace pressure. The reason you currently see a smaller value in the droplet is because “density1+density2” is two times the density of fluid 2, which is zero inside the droplet.
I hope this helps.