[m.r.soufivand]

As in the article " https://doi.org/10.1016/j.cep.2019.04.015 ", I want to simulate direct condensation of steam to subcooled water spray.

The results are good when I simulate only hydrodynamics, but when I compile the UDF to mass transfer and add source terms, In the first few iterations, it diverges and gives the following error:
"WARNING: Invalid cp (0.000000e+00 J/kgK) for water-vapor at temperature nan K
Error at host: floating point exception
Error: floating point exception
Error Object: #f"



I have attached the UDF to this post. I would be grateful if you could help me.