The model is implemented in Aspen Custom Modeler and solved with the standard steady-state solver. The step size of the discretization of the spatial domain in the pellet is chosen as large as possible up to a minimum of 100 discrete points or a maximum step size of 2.5 μm. The CFD4 algorithm is used to calculate the derivatives.