Please present your code that has the convergence study included.
I’ve shown in other examples how to use ufl to set up the analytical solution. See for instance: Poisson problem with Neumann boundary condition - #6 by dokken
or
Error control: Computing convergence rates — FEniCSx tutorial