The kinetic model is implemented in Python 3.11 as well and the mass balance equations are solved using the initial value problem solver scipy.integrate.solve_ivp. In the following simulations, parameters as shown in Table are used unless otherwise stated. The reactant ratio and the sorbent ratio are calculated as given in eqs. ( ) and (