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

Traceback (most recent call last):

"/home/Fenics/Code_Python_cluster/Sintering/", 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/", line 220, in solve
    _solve_varproblem(*args, **kwargs)
  File "/home/anaconda3/envs/fenics11/lib/python3.11/site-packages/dolfin/fem/", line 247, in _solve_varproblem

*** -------------------------------------------------------------------------
*** DOLFIN encountered an error. If you are not able to resolve this issue
*** using the information listed below, you can ask for help at
*** 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'}

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.