Normal flux method at the boundary for SPH