Hi,
First of all, as you haven’t disclosed the mesh, the code isn’t reproducible. Couldn’t you just use dolfinx.mesh.create_rectangle for now?
You are also not setting any solver options.
Note that you could also compare with: Bad convergence of stokes problem
that solves a similar problem