Runtime Error, fenics.2019.**. How can I solve this problem?

Errors
1

Traceback (most recent call last):
File

"/home/Fenics/Code_Python_cluster/Sintering/PolyParticle.py", line 725, in <module>
    solve(CH_WeakForm_L == CH_WeakForm_R, u_new, solver_parameters=PhaseFields_params)
  File "/home/anaconda3/envs/fenics11/lib/python3.11/site-packages/dolfin/fem/solving.py", line 220, in solve
    _solve_varproblem(*args, **kwargs)
  File "/home/anaconda3/envs/fenics11/lib/python3.11/site-packages/dolfin/fem/solving.py", line 247, in _solve_varproblem
    solver.solve()
RuntimeError: 

*** -------------------------------------------------------------------------
*** DOLFIN encountered an error. If you are not able to resolve this issue
*** using the information listed below, you can ask for help at
***
***     fenics-support@googlegroups.com
***
*** Remember to include the error message listed below and, if possible,
*** include a *minimal* running example to reproduce the error.
***
*** -------------------------------------------------------------------------
*** Error:   Unable to solve linear system using PETSc Krylov solver.
*** Reason:  Solution failed to converge in 10000 iterations (PETSc reason DIVERGED_ITS, residual norm ||r|| = 7.672363e+80).
*** Where:   This error was encountered inside PETScKrylovSolver.cpp.
*** Process: 0
*** 
*** DOLFIN version: 2019.1.0
*** -------------------------------------------------------------------------

The solver I choose for this equataion is:

PhaseFields_params = {'linear_solver': 'tfqmr', 'preconditioner': 'ilu'}
2

Please read Read before posting: How do I get my question answered?

As a first step, you may want to change your PhaseFields_params to use a direct solver, and see if that works or not.