How to simply plot velocity profil

2

As you have created a numpy array you can use matplotlib.

However, you can export the function to paraview with XDMFFile.write_checkpoint you can visualize the velocity field at every degree of freedom.