We simulated reversible equilibrium reactions among these chemical species ... Due to the high level of non-linearity and large number of variables, the Newton-Raphson (NR) method was used to solve ...

In absence of reaction terms, each equation has an equilibrium point: The time constants τ are derived from the literature (see Table 1), and we assume them to be typical values for the specific kind ...