The solution speed is very slow

I would suggest looking at: How to choose the optimal solver for a PDE problem? - #2 by nate
and
Default absolute tolerance and relative tolerance - #4 by nate
as it gives a lot of general guidance.

If you want to understand what happens in a “Newton solve” i would suggest reading Custom Newton solvers — FEniCSx tutorial