Stefano
Rebughini‡
^{ab},
Mauro
Bracconi‡
^{a},
Anthony G.
Dixon
*^{b} and
Matteo
Maestri
*^{a}
^{a}Laboratory of Catalysis and Catalytic Processes, Dipartimento di Energia, Politecnico di Milano, via La Masa 34, 20156, Milano, Italy. E-mail: matteo.maestri@polimi.it
^{b}Department of Chemical Engineering, Worcester Polytechnic Institute, Worcester, MA, 01609 USA. E-mail: agdixon@wpi.edu

Received
29th November 2017
, Accepted 20th December 2017

First published on 8th January 2018

Hierarchical modeling is applied for the investigation of micro packed bed reactors. This method allows for the use of Computational Fluid Dynamics (CFD) simulations in the analysis of representative complex geometries, where a full-scale CFD simulation of the entire reactor is not possible. Detailed and computationally demanding analyses are used to study a selected number of conditions and phenomena. Then, lumped parameters are derived from CFD results by means of engineering correlations. These parameters are incorporated in simplified reactor models based on macroscopic conservation equations. We provide evidence for the potential of the approach by using as a show-case a micro packed bed reactor in the context of highly exothermic selective oxidation processes. This reactor configuration consists of catalytic particles packed in the channels of a honeycomb matrix, which is expected to strongly enhance the radial heat transfer. In particular, we first focus on the analysis of energy transfer mechanisms by CFD and their interpretation via a 1D model and we provide an assessment of existing correlations with respect to the unconventional configuration (2 mm channel equivalent diameter and 0.8 mm sphere diameter). These correlations are then implemented in a pseudo-continuous (i.e. macroscopic) 2D model to allow for a systematic investigation of the capabilities of the micro packed bed reactor in dealing with the selective oxidation of o-xylene to phthalic anhydride. We found that due to the enhanced radial heat transfer micro packed bed reactors allow for quasi-isothermal operations, thus extending the range of operating conditions possible without occurring in adverse thermal behavior of the reactor. On a more general basis, we prove that the hierarchical approach to chemical reactor engineering is an effective tool to bring the application of fundamental modeling at a level of complexity relevant to full-scale applications, otherwise not possible because of the impractical computational costs.

In previous work,^{3} Maestri and co-workers proposed a hierarchical modeling approach as an effective route to overcome this problem, thus enabling a practical implementation of the fundamental multiscale modeling approach. It consists of using detailed computationally demanding methods – based on computational fluid dynamics simulations (CFD) of the reactor geometry – to first study a selected and limited number of conditions. Then, CFD analysis is used for the derivation of lumped parameters and correlations to be adopted in macroscopic reactor modeling. In particular this approach has been successfully applied to assess the application of gas-to-solid heat and mass transport correlations derived for packed beds in the context of micro packed bed reactors.^{3} In this work, we extend the application of the hierarchical approach to the full-scale simulation of a micro-channel packed bed reactor with external heat transport in order to rate its performance with respect to that of a traditional multi-tubular reactor.

A micro packed bed reactor is a configuration suggested by Vervloet et al.^{4} to enhance the heat transfer properties of a conventional packed bed reactor and to increase the catalytic load per reactor volume of a structured reactor (e.g., honeycomb reactor).

Fig. 1 shows the schematic representation of this reactor configuration, which consists of small spherical catalytic particles packed into the channel of a metallic honeycomb matrix.^{5} The advantage of this reactor configuration would be two-fold. On one side, the honeycomb matrix is expected to show the enhanced heat transfer performance.^{6} In fact, as demonstrated by Groppi and Tronconi,^{7} the metallic connected solid matrix allows for quasi-isothermal operation of highly exothermic processes (e.g. oxidation of methanol to formaldehyde and epoxidation of ethylene) because of the improved effective radial conductivity. On the other side, the catalyst load per reactor volume is similar to conventional packed bed reactors. Indeed, the overall load of active catalyst per unit reactor volume in a micro packed bed reactor is 0.3–0.4 v/v, which is more similar to typical catalytic loads of industrial packed bed reactors (0.4–0.5 v/v) than those of structured reactors (typically ∼0.1 v/v for honeycomb reactors). As such, this reactor configuration, in principle, could enable quasi-isothermal operation of selective oxidation reactions while keeping a catalytic load per reactor volume similar to conventional packed bed reactors. Here, we first extend the previous analysis^{3} used for particle-to-gas energy and mass transport, to the assessment of existing correlations for the description of the external heat transfer phenomena. All these hierarchically derived correlations are then implemented in a steady-state pseudo-continuous 2D model for the full-scale simulation of the micro packed bed reactor. This hierarchical model, which is less computationally demanding than CFD simulations, is then used to investigate the performances of the reactor configuration for the selective oxidation of o-xylene to phthalic anhydride.^{2}

This study demonstrates how the hierarchical modeling approach can be applied to investigate in detail the capabilities of novel reactor geometries, by enabling the use of CFD insights at a level of complexity otherwise not accessible because of the unbearable computational cost. It also proves that CFD can be employed to assess and screen existing correlations, thus enabling a rational derivation or selection of the most adequate one with respect to the particular set of operating conditions and reactor geometry. In particular, a CFD resolved-particle model of a single packed channel of a micro packed bed reactor is used to provide simulations of heat transfer, against which existing heat transfer correlations are evaluated using a 1D heat transfer model. The best results are then incorporated into a 2D effective continuum reactor model for o-xylene partial oxidation to phthalic anhydride.

Catalyst properties | ||
---|---|---|

Thermal conductivity | 1.5 | W m^{−1} K^{−1} |

Density | 2100 | kg m^{−3} |

Specific heat at constant pressure | 925 | J kg^{−1} K^{−1} |

Fig. 2 Insight of the computational grid used for CFD simulations with the mesh discretization in both fluid and solid domain. |

Given the low Reynolds numbers considered in the simulations, prism layers were not created near the walls.^{13} At the fluid–solid interface the boundary conditions for the heat diffusion fluxes suggested by Murthy and Mathur^{14} are adopted. Non-slip and Neumann boundary conditions are applied at the walls for velocity and pressure, respectively. The CFD simulations are performed by keeping constant the temperature of the channel walls and with the operating conditions and reactor geometries reported in Tables 2 and 3, respectively.

Operating conditions | ||
---|---|---|

Feed temperature | 293 | K |

External temperature | 295 | K |

Outlet pressure | 1 | atm |

Feed gas species | Air | |

Reynolds numbers | 70–220 |

D
_{p} (mm) |
D
_{c} (mm) |
L (mm) | N |
---|---|---|---|

0.6 | 2.0 | 8.0 | 3.33 |

0.6 | 2.5 | 8.0 | 4.17 |

0.6 | 3.0 | 8.0 | 5.00 |

(1) |

In eqn (1)ρ^{G} is the gas density, u the flow velocity, C^{G}_{p} the specific heat at constant pressure for the gas phase, T^{G} the temperature of the gas phase, T_{wall} the honeycomb matrix temperature, S_{V} the specific heat exchange area between the packed channel and the honeycomb matrix and U is the overall heat transfer coefficient. The operating conditions have been selected to ensure constant gas properties in the channel, thus eqn (1) can be re-written as:

(2) |

The temperature at the reactor inlet (T^{0}) and at the reactor outlet (T^{G}) in eqn (2) are estimated directly from the CFD simulations as mixing-cup average:^{15}

(3) |

(4) |

(5) |

In eqn (4) and (5)G is the specific mass flow rate, ρ^{G} the gas density, ε^{H} the void fraction of the honeycomb matrix, and ε^{P} the void fraction of the packed channel. In eqn (4)ω^{G}_{i} is the mass fraction of the i-th species in the gas phase and R_{i} the reaction rate for the i-th species. In eqn (5)T^{G} is the temperature of the gas phase, T^{H} the temperature of the honeycomb matrix, C^{G}_{p} is the specific heat at constant pressure for the gas phase, S_{V} specific the heat exchange area between the packed channel and the honeycomb matrix, Q^{REACT} the reaction heat and U the overall heat transfer coefficient. The size of the pellets adopted in the micro packed bed reactors results in negligible internal transport limitations due to the very small characteristic diffusion length.

The following boundary conditions at the reactor inlet (z = 0) are applied for eqn (4) and (5), respectively:

ω^{G}_{i}(z = 0) = ω^{G,0}_{i} | (6) |

T^{G}(z = 0) = T^{G,0} | (7) |

The axial dispersion is neglected due to the high Péclet number. The gas-phase homogeneous reactions and radiative heat transfer are negligible, as reported in ref. 2 and 16.

The honeycomb matrix is described as a continuum^{17} consisting of a static, thermally connected solid phase.^{7} The heat balance for the honeycomb matrix is given by the following equation:

(8) |

In eqn (8)ρ^{H} and C^{H}_{p} are the density and the specific heat at constant pressure of the honeycomb matrix, respectively. k^{H}_{eff,ax} and k^{H}_{eff,r} are the axial and radial thermal conductivity of the honeycomb matrix. The boundary conditions applied to eqn (8) are:

Inlet conditions at z = 0:

(9) |

Outlet conditions at z = L:

(10) |

Symmetry condition at r = 0:

(11) |

Wall conditions at r = R:

(12) |

In eqn (12)T_{cool} is the coolant temperature and h_{w} is the wall heat transfer coefficient between the honeycomb matrix and the external coolant.

In this 2D model, the solid properties are considered constant and reported in Table 4. The thermo-physical properties of the ideal gas are evaluated by using the OpenSMOKE++ library.^{18}

Packed bed^{2} |
Micro packed bed | |
---|---|---|

Tube diameter (cm) | 2.54 | 2.54 |

Particle diameter (cm) | 0.4 | 0.08 |

CPSI/wall thickness | (not required) | 50/25 |

Length (m) | 3.0 | 1.42 |

Specific mass flow rate (kg m^{−2} s^{−1}) |
1.3 | 0.45 |

Catalyst conductivity (W m^{−1} K^{−1}) |
1.5 | 1.5 |

Support conductivity (W m^{−1} K^{−1}) |
(not required) | 100 |

The partial differential equations (PDEs) characterizing this model are solved by using the method of lines.^{19} The axial and radial reactor lengths are split into a suitable number of grid points and in each point the spatial derivatives are approximated by the Euler backward differentiation method. Thus, the PDEs problem is converted into a Differential Algebraic Equations (DAEs) system. The DAE system is solved by using the OpenSMOKE++ library, whose integrator can exploit the tridiagonal block structure of the Jacobian and which resulted in very efficient and robust convergence.^{18}

a) Heat exchange between the single packed channel and the honeycomb matrix (Fig. 1a), accounted for by an overall heat transfer coefficient (U) in eqn (5) and (8);

b) Heat transfer within the honeycomb matrix (Fig. 1b), represented by axial (k^{H}_{eff,ax}) and radial (k^{H}_{eff,r}) effective thermal conductivities in eqn (8);

First, we apply the hierarchical modeling approach to derive these lumped parameters. Then, these parameters are implemented in the 2D model, which is used to analyze the behavior of micro packed bed reactors during selective oxidation processes.

For the effective thermal conductivities, we rely on the previous investigation of Groppi and Tronconi for a honeycomb matrix.^{20} Thus, the following correlations for the axial and radial effective thermal conductivity are employed:

k^{H}_{eff,ax} = k^{H}(1 − ε^{H}) | (13) |

(14) |

The correlations above result from a theoretical analysis of the heat conduction in the unit square cell of the honeycomb matrix and they have been applied to investigate the heat transfer phenomena in externally cooled monolith reactors. They have been validated by using experimental results^{7,17,21,22} and with 3D simulations of the honeycomb structure.^{23} Therefore, these correlations are a reliable description of the heat transfer within the honeycomb matrix and we implement them in the 2D model of micro packed bed reactors.

The hierarchical analysis of this heat transfer phenomenon is performed as follows. First, we describe the heat transfer in the micro packing by means of CFD simulations (high hierarchy model) of a single channel of the honeycomb matrix packed with spherical particles. Temperature contours for a tube-to-particle ratio of 4.1 and two different Reynolds numbers are reported Fig. 3. Then, the CFD calculated profiles are analyzed with the 1D pseudo-homogeneous model (low hierarchy model) to estimate the overall heat transfer coefficients by exploiting eqn (1) and (2). These values are now compared with the ones obtained from literature correlations derived for conventional packed bed reactors. We select the correlation proposed by Martin–Zehner,^{24,25} widely adopted for the description of packed beds,^{26} and the correlation by Wellauer et al.,^{27} derived for externally cooled packed bed reactor working at low Reynolds numbers. This condition is close to the one of micro-packed bed reactor, where the small dimension of the pellet employed leads to low Reynolds numbers.

Fig. 3 Contours of temperature on a mid-tube plane through the bed for a tube-to-particle ratio of 4.17 and Reynolds number equal to 110 (a) and 180 (b). |

The Martin–Zehner^{24,25} correlation consists of the wall heat transfer coefficient (h_{w}) estimated with Martin and Nilles^{24} correlation coupled with the formula of Zehner and Schlünder^{25} for the effective radial thermal conductivity (k^{H}_{eff,r}). Wellauer et al.^{27} proposed alternative methods for the evaluation of the wall heat transfer coefficient and of the effective radial thermal conductivity derived in the context of low Reynolds number flows. The equations involved in this correlation are reported in the ESI† (section 1). The combination of these two expressions is based on the Biot number (Bi) as suggested by Dixon^{28} with eqn (15).

(15) |

The prediction of the overall heat transfer coefficient of the correlations of Martin–Zehner^{24,25} and Wellauer et al.^{27} show large deviations, up to 300 W m^{−2} K^{−1}, in the range of operating conditions of micro packed bed reactors, as shown in Fig. 4. In this case, we found out that the values estimated with the correlation of Wellauer et al.^{27} are more similar to the CFD-based results, such as in the range of ±10 W m^{−2} K^{−1}, than those evaluated with the Martin–Zehner^{24,25} correlation. This difference can be ascribed to the flow regime and the operating conditions used to derive the literature correlations. In particular, in terms of flow regime, the correlation of Martin–Zehner^{24,25} has been proposed to describe industrial packed bed reactors, which operate at high Reynolds numbers. Instead, Wellauer et al.^{27} proposed their correlation to describe externally cooled packed bed reactor working at low Reynolds numbers. The micro packed bed reactors are characterized by Reynolds numbers between 30 and 130, which are around the ones investigated in the derivation of the correlation of Wellauer et al.^{27} Moreover, in terms of reactor geometry, the Martin–Zehner^{24,25} correlation describes packed bed with higher tube-to-particle diameter ratio than those represented by the correlation of Wellauer et al.^{27} Thus, the correlation suggested by Wellauer et al.^{27} can better reproduce the CFD-based results because it has been derived from experiments with tube-to-particle diameter ratios similar to those of micro packed bed reactors. The heat exchange between the packed channel and the honeycomb matrix in a micro packed bed reactor can be described by using the correlation suggested by Wellauer et al.^{27} The accuracy of the conventional models is guaranteed since the results of the CFD analysis allows for selecting the empirical correlation which properly accounts for the geometries and operating conditions of the investigated reactor.

This process is industrially carried out in multi-tubular reactors cooled by molten salts. The reaction pathways that lead to the total combustion products are more activated and more exothermic than the ones responsible for phthalic anhydride production. Thus, the temperature control is crucial to ensure high selectivity to the desired product (phthalic anhydride). The kinetic parameters and the rate expressions are taken from Froment et al.,^{2} more details can be found in the ESI† (section 3). Fig. 5 shows the temperature profile obtained in the classical reactor configuration (parameters in Table 4) compared to the one evaluated in a micro-packed bed reactor.

Fig. 5 Temperature axial profile for the packed bed reactor (dashed line) and the micro packed bed reactor (solid line). |

The operating conditions for the micro packed reactor have been selected to guarantee to the micro packed bed the same Gas Hourly Space Velocity (GHSV) and the same pressure drop as the packed bed reactor. Moreover, the comparison between these two reactor configurations is performed by considering an o-xylene-to-air feed molar ratio of 0.009 and a coolant temperature of 625 K, which are the typical industrial conditions as reported by Froment et al.^{2}

The packed bed configuration is characterized by a hot spot in the temperature profile. Instead, the micro packed bed reactor results in an almost isothermal profile. This different behavior depends on the amount of heat removed in the two reactor configurations. In fact, the heat extracted from the packed bed is initially lower than the heat generated by the reaction. Once the temperature difference between the reactor and the coolant becomes sufficient, the heat is properly removed decreasing the temperature in the reactor. This opposite behavior generates the hot-spot in the packed bed. On the contrary, the metallic honeycomb matrix of the micro packed bed increases the heat removal, which enables the quasi-isothermal operation of this configuration. This is due to the high overall heat transfer coefficient of the micro packed bed which allows for a large heat exchange even in presence of a small (<10 K) temperature difference between the reactor and the coolant. The overall heat transfer coefficient relies on three different contributions in series: the heat transfer between the packed bed and the solid matrix of the honeycomb, the heat conduction in the metallic matrix and the heat exchange between the monolith and the coolant. A tight fitting of the honeycomb in the reactor tube enables a perfect contact between the monolith and the tube wall. The radial heat transfer is strongly promoted by the metallic structure of the honeycomb through the conduction mechanism, when the contribution of the heat transfer between the packed bed and the solid matrix is sufficiently high. The proper design of the micro packed bed is pivotal to provide an intense gas to matrix heat transfer boosting the heat exchange phenomena by means of the conduction in the solid.

This results in an enhanced radial heat transfer in the micro packed bed allowing for a wider range of operating conditions. This is illustrated in Fig. 6 where the maximum temperature for both configurations is reported as a function of the coolant temperature. In the packed bed reactor, a coolant temperature higher than 637 K leads to a thermal loss of the reactor caused by the intensification of the combustion reactions to CO_{x} resulting in a dramatic decrease of the selectivity coupled to an uncontrolled increment of the temperature. On the contrary, the coolant temperature for the micro packed bed reactor can be increased up to 677 K keeping the maximum temperature in the reactor lower than 685 K, when the external heat transfer resistance between the monolith and the coolant is negligible. This behavior depends on the increased radial heat transfer, which allows for better heat removal and enables a quasi-isothermal operation of the reactor, even with the same catalytic load adopted in the packed bed reactor.

The previous analysis has been carried out with an optimal contact between the monolith and the tube. We also analyzed the effect of the presence of a gas gap between the honeycomb and the tube wall or the onset of fouling on the surface between the tube and the coolant. These phenomena introduce additional thermal resistances which reduce the heat performances of the reactor. Despite this, even considering an external heat transfer value of 700 W m^{−2} K^{−1},^{20} the behavior of the reactor is poorly affected as shown in Fig. 6 (full square symbols). The maximum temperature difference between the gas and the coolant remains below 10 K and the quasi-isothermal operating condition in the reactor are kept. The runaway is predicted at 669 K due to the lower heat flux between the monolith and the coolant. However, this is still a major improvement in the thermal behavior with respect to the conventional packed bed. Moreover, the increased heat removal affects the selectivity and yield to the desired product (phthalic anhydride) of the reactor. Fig. 7 reports calculated selectivity, yield, and conversion for the micro packed bed reactor as a function of the coolant temperature in the case of negligible external heat transfer resistance. It is shown that increasing the coolant temperature the conversion of o-xylene increases too but also reduces the selectivity to phthalic anhydride, which is the typical behavior of this reaction scheme (see ESI† section 3). Thus, an optimal temperature for the operation of the reactor is observed leading to a maximum of the phthalic anhydride molar yield. Fig. 7 shows that this value is 0.695 and it corresponds to a constant coolant temperature of 661 K. In particular, this value is slightly higher that the yield of industrial multi-tubular packed bed reactors (0.690)^{2} and quite similar to the value of 0.696 that can be obtained with an isothermal reactor. Moreover, the micro packed bed reactor enables a more stable operativity of the reactor strongly decreasing the hot-spot temperature and increasing the robustness to temperature variations in the reactor. Focusing our attention on the temperature profiles inside the reactors, despite the similar yield values between the multi-tubular packed bed and the micro packed bed rector, the second one is almost isothermal (T_{max}–T_{cool} < 1.5 K). This analysis relies on the adequacy of the correlation used for the description of the heat transfer in the micro packed bed reactor. Fig. 8 shows the effect of two different correlations on the predicted temperature profile. Temperature differences between 2 and 18 K are experienced within the reactor due to the different estimation of the heat transfer coefficients. Above all, the Martin–Zehner correlation^{24,25} usually overpredicts the reactor temperature leading to an earlier runaway of the reactor at 661 K, clearly below the limit of 681 K obtained with the Wellauer method,^{27} which was selected using the hierarchical approach. The results obtained with the Martin–Zehner^{24,25} strongly underestimate the heat exchange masking the potential of micro packed bed reactors. In particular, the runaway of the reactor is predicted at lower temperatures, which actually correspond to the optimum of the yield evaluated using the Wellauer^{27} correlation. Thus, the selection of the correlation has a strong impact on the analysis of the reactor, hindering the operability of the reactor in these convenient conditions.

Fig. 7 Effect of the coolant temperature on the o-xylene conversion, phthalic anhydride selectivity and phthalic anhydride yield for the micro packed bed reactor. |

Fig. 8 Correlation effect on the maximum temperature in the micro packed bed reactor as a function of the coolant temperature. Wellauer et al.^{27} solid line and full circular symbol. Martin–Zehner^{24,25} solid line and empty square symbol. |

Thus, the capabilities of the hierarchical modeling are evident being able to define the most adequate correlation with a reasonable computational effort.

Our analysis demonstrates that the high heat removal of the honeycomb matrix allows a quasi-isothermal operation of the reactor. The connected solid matrix should also reduce the temperature radial gradients by the enhanced radial heat transfer, which allows for a significant increment of the reactor tube diameter. To better clarify this point, we performed a parametric analysis on the reactor geometry by increasing the tube diameter. In particular, Fig. 9 illustrates the maximum temperature and the yield as a function of the coolant temperature for different tube diameters. We observe that increasing the tube diameter slightly affects the maximum temperature in the reactor. In fact, Fig. 9 shows that the maximum temperature in the reactor increases while increasing the tube diameter. This is associated with the increased distance between the center of the reactor and the molten salts. Therefore, even with the enhanced heat transfer properties of the metallic matrix, the center of the reactor reaches higher temperature for higher tube diameter. Despite this trend in the temperature profile, increasing the tube diameter does not remarkably affect the reactor performance in terms of phthalic anhydride yield, which reaches the maximum value at 661 K for all the tube diameters investigated. In fact, even in the worst scenario of a tube diameter of 7.62 cm at a coolant temperature of 661 K, the micro packed bed is quasi-isothermal (T_{max}–T_{cool} < 5.5 K).

On a more general basis, this analysis clearly demonstrates that the hierarchical modeling approach can be adopted to analyze complex reactor configurations, when fundamental CFD modeling is still limited by the impractical computational cost. In particular, the hierarchical modeling has been applied to investigate micro packed bed reactors by comparing them to multi-tubular packed bed reactors. Even if this comparison is based on a simplified kinetic scheme and on a limited number of operating conditions, this preliminary analysis shows that the micro packed bed reactors can be a suitable alternative to multi-tubular packed bed reactors.

- M. P. Dudukovic, Science, 2009, 325, 698–701 CrossRef CAS PubMed .
- G. F. Froment, K. B. Bischoff and J. De Wilde, Chemical Reactor Analysis and Design, John Wiley & Sons, Incorporated, 3rd edn, 2010 Search PubMed .
- S. Rebughini, A. Cuoci and M. Maestri, Chem. Eng. J., 2016, 289, 471–478 CrossRef CAS .
- D. Vervloet, F. Kapteijn, J. Nijenhuis and J. R. Van Ommen, Catal. Today, 2013, 216, 111–116 CrossRef CAS .
- M. Iovane, R. Zennaro, P. Forzatti, G. Groppi, L. Lietti, E. Tronconi, C. G. Visconti, S. Rossini and E. Mignone, US Pat., 20120184631, 2012 Search PubMed .
- K. Pangarkar, T. J. Schildhauer, J. R. Van Ommen, J. Nijenhuis, F. Kapteijn and J. A. Moulijn, Structured packings for multiphase catalytic reactors, Ind. Eng. Chem. Res., 2008, 47(10), 3720–3751 CrossRef CAS .
- E. Tronconi and G. Groppi, Simulation of structured catalytic reactors with enhanced thermal conductivity for selective oxidation reactions, Catal. Today, 2001, 69, 63–73 CrossRef .
- H. G. Weller, G. Tabor, H. Jasak and C. Fureby, A tensorial approach to computational continuum mechanics using object-oriented techniques, Comput. Phys., 1998, 12, 620–631 CrossRef .
- M. W. J. Chase, NIST-JANAF Thermochemical Tables, American Inst. of Physics, 1998 Search PubMed .
- C. Kloss, C. Goniva, A. Hager, S. Amberger and S. Pirker, Pcfd, 2012, 12, 140–152 CrossRef .
- S. Rebughini, A. Cuoci and M. Maestri, Chem. Eng. Sci., 2016, 141, 240–249 CrossRef CAS .
- A. G. Dixon, M. Nijemeisland and E. H. Stitt, Comput. Chem. Eng., 2013, 48, 135–153 CrossRef CAS .
- S. A. Logtenberg, M. Nijemeisland and A. G. Dixon, Chem. Eng. Sci., 1999, 54, 2433–2439 CrossRef CAS .
- J. Y. Murthy and S. R. Mathur, J. Heat Transfer, 1998, 120, 583–591 CrossRef .
- R. K. Shah and A. L. London, Laminar flow forced convection in ducts: a source book for compact heat exchanger analytical data, Academic press, 1978 Search PubMed .
- G. Groppi and E. Tronconi, Chem. Eng. Sci., 2000, 55, 2161–2171 CrossRef CAS .
- G. Groppi and E. Tronconi, AIChE J., 1996, 42, 2382–2387 CrossRef CAS .
- A. Cuoci, A. Frassoldati, T. Faravelli and E. Ranzi, Comput. Phys. Commun., 2015, 192, 237–264 CrossRef CAS .
- J. H. Ferziger and M. Peric, Computational Methods for Fluid Dynamics, Springer, Berlin Heidelberg, 2001 Search PubMed .
- E. Tronconi, G. Groppi and C. G. Visconti, Curr. Opin. Chem. Eng., 2014, 5, 55–67 CrossRef .
- G. Groppi, A. Belloli, E. Tronconi and P. Forzatti, Chem. Eng. Sci., 1995, 50, 2705–2715 CrossRef CAS .
- E. Tronconi and G. Groppi, Study on the thermal behavior of structured plate-type catalysts with metallic supports for gas/solid exothermic reactions, Chem. Eng. Sci., 2000, 55, 6021–6036 CrossRef CAS .
- D. Schlereth and O. Hinrichsen, Ind. Eng. Chem. Res., 2014, 53, 11550–11556 CrossRef CAS .
- H. Martin and M. Nilles, Chem. Ing. Tech., 1993, 65, 1468–1477 CrossRef CAS .
- P. Zehner and E. U. Schlünder, Chem. Ing. Tech., 1970, 42, 933–941 CrossRef CAS .
- A. G. Dixon, Can. J. Chem. Eng., 2012, 90, 507–527 CrossRef CAS .
- T. Wellauer, D. L. Cresswell and E. J. Newson, in Chemical Reaction Engineering—Boston, American Chemical Society, 1982, vol. 196, pp. 41–527 Search PubMed .
- A. G. Dixon, Chem. Eng. Process., 1996, 35, 323–331 CrossRef CAS .

## Footnotes |

† Electronic supplementary information (ESI) available: Equations of Wellauer correlation, details on notation and symbols, kinetic scheme. See DOI: 10.1039/c7re00195a |

‡ These authors contributed equally to this work. |

This journal is © The Royal Society of Chemistry 2018 |