[maybe]

Hi, I am maybe and I am trying to simulate the combustion of an Aluminum droplet using FLUENT 2024.
I developed a chemical kinetic mechanism including chemistry, thermodynamics, transportation, and surface reactions files in CHEMKIN format. And I try to see how the mechanism performs in a 2D simulation.
Firstly, I tried to simulate with the reaction option off, the computation went successfully and converged in 200 steps.

However, when I turn on the reaction option, errors came out in only a few steps. Here is the information.

  iter  continuity  x-velocity  y-velocity      energy           h          h2           o          o2           c          co         co2          oh         ho2         h2o        h2o2       al<g>         alo        alo2        al2o       al2o2    al2o3<g>         alh        alh2        alh3        aloh     al<oh>2     al<oh>3       alo2h    al2o3<l>       al<l>           n          n2     time/iter
     1  1.0000e+00  6.9505e-03  1.2445e-04  6.9018e-05  0.0000e+00  0.0000e+00  0.0000e+00  1.7061e-06  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0:01:10    4
     2  7.1443e-01  3.7648e-03  8.9145e-05  9.1665e-05  0.0000e+00  0.0000e+00  0.0000e+00  9.6128e-01  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0:00:47    3
     3  4.8746e-01  2.0723e-03  7.6055e-05  1.0010e-04  0.0000e+00  0.0000e+00  0.0000e+00  4.7635e-01  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0:00:28    2
     4  3.2844e-01  1.1713e-03  6.2949e-05  1.0750e-04  0.0000e+00  0.0000e+00  0.0000e+00  3.1517e-01  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0.0000e+00  0:00:13    1
 bad sensitivities detected, kisat =            7
 mga =               -Infinity
 need =           1           1           1
 nx, nf, nh, iusr =           31          31           0           1
 xe =   23308800.7207318       1.075652459236753E-016  9.601353022704363E-018
  9.917178022035544E-017  2.241449632615759E-003  1.649842102830128E-017
  4.599166856980210E-017  4.042684444173737E-014  2.437783169727147E-018
  1.215205181698718E-016  4.356661446808399E-018  2.038071443603131E-021
  0.000000000000000E+000  3.392791192505031E-014  1.348725289701485E-016
  2.435644795460029E-017  7.088183581774156E-016  5.407546209926491E-017
  3.467964024796940E-016  5.059852344389265E-020  4.731723358178852E-021
  1.780154937894211E-019  2.065252034350235E-016  1.881015166757336E-017
  2.829687906474808E-016  6.481443523899239E-016  1.206888015566587E-017
  3.684724508601841E-016  3.454129479019592E-017  2.323711578783470E-002
  1.266347171814820E-005
 fe =   94.0000000000000       1.977884983624166E-011   101325.000000000     
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  3.952525166729972E-323
  3.952525166729972E-323  0.000000000000000E+000  0.000000000000000E+000
  2.121995790965272E-314  4.047879836377333E-320  0.000000000000000E+000
  0.000000000000000E+000  0.000000000000000E+000  0.000000000000000E+000
  0.000000000000000E+000  5.928787750094959E-323  7.905050333459945E-323
  0.000000000000000E+000  0.000000000000000E+000  1.482196937523740E-323
  9.881312916824931E-324  9.881312916824931E-324  0.000000000000000E+000
  0.000000000000000E+000
fa =                     NaN                     NaN                     NaN
                   ...(a lot of NaN)
 ga =                     NaN                     NaN                     NaN
                    ..(a lot of NaN again)
        Stabilizing species-0 to enhance linear solver robustness.
        Stabilizing species-0 using GMRES to enhance linear solver robustness.

Divergence detected in AMG solver: species-0        Stabilizing species-1 to enhance linear solver robustness.
        Stabilizing species-1 using GMRES to enhance linear solver robustness.

Divergence detected in AMG solver: species-1        Stabilizing species-2 to enhance linear solver robustness.
        Stabilizing species-2 using GMRES to enhance linear solver robustness.

........similar infomation............

Divergence detected in AMG solver: species-27        Stabilizing temperature to enhance linear solver robustness.
        Stabilizing temperature using GMRES to enhance linear solver robustness.

Divergence detected in AMG solver: temperature
Error at host: floating point exception

Error at Node 0: floating point exception

Error: floating point exception
Error Object: ()

Error: floating point exception
Error Object: #f

I couldn't find out where the problem is. I have tried to lower the URF, change a denser mesh (all mesh quality is good).

Could anyone help me with this problem?
Thanks!