Someone in [this page][1] has claimed to successfully implement [Brent solver][2] in [thermoI.H
][3] instead of Newton solver to resolve the infamous
FOAM FATAL ERROR: Maximum number of iterations exceeded
issue. The part of code he is referring to is:
template<class Thermo, template<class> class Type>
inline Foam::scalar Foam::species::thermo<Thermo, Type>::T
(