“…The steady-state 2D CFD model employed in our study is based on the previously reported CFD models [26,27,34]. We assume a pseudo-homogeneous description of the catalytic bed [22,29], therefore the computational domain representing the catalytic packed bed is solved as a continuum.…”