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!