Classical rescaling of phase field distribution to 1/3 rho as common in Lee-literature. The source terms were removed from the moments and returned to the distributions based on the original code of Safari.

fix the issue in FullVectorConnector derived classes: add initialization function of base class in derived classes

MultiphaseVelocityBCAlgorithm::applyBC(): vleocity in phase field taken from getBoundaryVelocityX1(). There is still something wrong as the actual velocity in the simulation is much higher than set in the input file!