A robust PDE solver for the 3D Stokes/Navier-Stokes systems on the grid environment