dirichlet p boundary conditions applied
projection method implemented