1. Introduction
Monolithic honeycomb structures are used extensively as supports in catalytic processes, especially in internal combustion engine exhaust aftertreatment applications [
1,
2]. Tubular monoliths possess higher open facial areas and greater mechanical strength relative to conventional pellets or packed beds. These features also allow high specific surface areas, low differential pressure, and relatively homogeneous thermal conditions even at high throughput [
3,
4,
5,
6,
7]. The open structure of monoliths also decreases weight and thermal mass relative to pellet media, which is important in mobile applications where both rapid heating and light-off improve conversion efficiency during transient operation. Ceramic monoliths are more commonly used than metallic types for automotive applications due to easier wash-coat adhesion [
6]. Despite this, corrugated metal monoliths offer thermal and design benefits, which may outweigh their disadvantages.
Metallic monolith performance surpasses that of ceramics in several ways. Thin foils used in their construction increase the open area, decreasing the pressure drop [
8] and weight, while maintaining a mechanically robust structure. Metals are more thermally conductive than ceramics, which decreases thermal stratification in the reactor [
6,
9]. Finally, metal monoliths can be formed using rolling techniques, which offer additional geometries not available through extrusion [
10]. Monolith cell profiles can be sinusoidal, square, or triangular, each of which offers unique flow characteristics. Exceptionally small cell densities are possible, which creates large specific surfaces area for high catalyst loading. Open-cell monoliths, which behave like metal foams, can be built using perforated sheet stock [
10,
11], offering inter-channel flow mixing at a relatively lower cost. These thermal and fluid factors offer additional degrees of freedom in catalyst manufacturing not available to conventional ceramics. Full comprehension of these factors is of great value and will inform the next generation of novel catalytic technologies.
Previous work in monolith catalyst heat transfer modeling has primarily taken two approaches. One approach is to resolve an individual channel of the monolith in high detail, capturing turbulence and thermal stratification throughout. Heat transfer between channels is then inferred by expanding the single channel result to the full scale of the reactor. Studies by Cornejo et al. [
12,
13] examined the heat transfer effects within single cells of metallic monoliths. Their investigations used computational fluid dynamics (CFD) to investigate the relation between the Nusselt number (Nu) and the flow characteristics in the channel. CFD is more computationally cost effective when used to resolve smaller portions of a larger system.
The second approach to model monolith heat transfer is to consider the entire reactor using CFD or quasi-dimensional simulations without resolving individual channels. For example, an automotive catalytic converter was modeled by Hayes et al. [
14] using CFD and 2D symmetry approximation. This model could solve simple flow cases in a matter of seconds; however, simulating catalyst light-off took up to 2 h. Modeling the entire reactor can be computationally less expensive under certain conditions but loses some of the finer physics and chemistry present at the channel level. This tradeoff is well illustrated by Sadeghi et al. [
15], where the same reactor conditions are modeled using increasing dimensionality and different mathematical approaches. Despite modeling the same conditions, adding dimensions to the approach changed nitric oxide (NO) concentration and conversion rates over the length of the reactor. In this work, we present a quasi-dimensional modeling approach of the second variety, where the whole heat exchanger was resolved using quasi-dimensional methods. At the same time, our approach resolved cellular channel effects using geometric and thermal assumptions developed in prior works for single channel analyses. The net benefit of this approach is a computationally lean model, which still resolves dimensional stratification of the entire reactor.
Although comprehensive studies are sparse, metal monolith heat transfer has been of interest for many decades. An early foundational study by R.K. Shah [
16] developed an initial understanding of thermal transfer in arbitrary shaped ducts. Here, steady-state Nusselt number (
Nu) values were determined for developed laminar flow in numerous duct geometries. These
Nu values are critical for monolith models, as
Nu is a constant, and can be applied to every repeated tubular cell in the monolith. Simple models can be developed using accurate estimation of gas properties and proper selection of
Nu. While fully developed
Nu values differ in the literature for different cellular profiles, the constant and uniform
Nu assumption can apply to a fixed profile and greatly simplifies modeling. For metal monoliths where the channel hydraulic diameter (
Dh) is much smaller than their length (
Dh << L), this assumption is accurate.
A later study by Boger and Heibel [
17] calculated the apparent convective heat transfer coefficients in copper, aluminum, and cordierite monoliths. Heated gas was flowed through the samples, which were cooled with a circulating water jacket. Results showed apparent heat transfer coefficients of 1000 W/m
2K or greater in metallic monoliths, which exceeded that of cordierite by an order of magnitude. They concluded that monolith thermal conductivity caused the higher measured heat transfer rates. Similar apparent rates were measured in a reactor study by Roh et al. [
18] using coated FeCrAl foil monoliths. The FeCrAl monolith demonstrated superior heat exchange when compared to a packed bed reactor operating under identical temperatures and flows.
More recent modeling studies have also explored coupling of dissimilar monoliths in both CFD and using genetic algorithm optimization. Regarding Cornejo et al. [
19] and Reinao et al. [
20], the use of concentric cylinder/ring paired monoliths with a single inflow of reactants was showed to affect catalyst conversion. The variation and fine-tuning of the relative fraction of the cylinder to ring can distribute heat optimally for a reaction. Regarding Cherif et al. [
21], application of a cheaper base-metal catalyst and conversion of methane within a cylindrical bed was optimized through a genetic algorithm approach. Careful sizing of concentric ring sections and varying catalyst type can create a more thermally uniform reactor, which performs similarly to a catalyst bearing only platinum materials. Together, these recent works demonstrate the importance of understanding concentric monolith heat exchange and the value of an increasing number of modeling procedures available for its study.
Catalyst substrate selection also has an important role in heat exchange within a reactor. Highly endothermic reactors such as steam reformers require high heating rates to maintain efficient conversion. Understanding and enhancing heat exchange allows both efficient operation and reactor intensification. This is especially crucial for waste-heat-driven compact reforming due to lower possible temperature differentials than in electrically heated or flame-heated reactors. For example, compact methane steam reforming (MSR) was achieved by Tonkovich et al. [
22] at temperatures near 850 °C. In contrast, normal waste heat temperatures available from diesel engines and turbines are only 600 °C [
23] and 650 °C [
24], respectively. The importance of support thermal conductivity was demonstrated in a study by Ryu et al. [
25], where a metallic monolith increased conversion and activity of MSR over that of a powdered catalyst bed. The metal monolith reactor also achieved this conversion while using only 18% of the powdered catalyst metal loading, significantly reducing the cost of such a reactor.
In this work, heat exchange between two metal monoliths brazed to a shared central mantle was experimentally measured and then modeled using a quasi-dimensional approach. Once model–experimental agreement was established, the model was then used to parametrically investigate heat exchange behavior of monoliths under varied thermal and physical constraints and determine the effects of each on the heat exchange rate and efficiency. The model is more computationally efficient than previous multi-dimensional approaches while providing dimensional insights not found in correlations and simpler models. Through its speed and general applicability, the developed approach expedites the design and construction process of experimental metal monolith reactors that include heat transfer by design. Furthermore, the model can be re-calibrated when new experimental data are available, increasing specificity and accuracy. By retaining only essential thermal and physical properties in a quasi-dimensional model, this work demonstrates a dimensionally resolved yet fast method for modeling thermal behavior in a metallic monolith reactor.
4. Results and Discussion
To calibrate the heat transfer model, a constant fully developed Nusselt number (Nu) was sought. A single Nu was found by minimizing root mean square error (RMSE) across the entire SRM experimental dataset. Once a satisfactory Nu was established, the calibrated model was then subjected to a new full-factorial study where geometric parameters of the modeled monolith were varied. Effects of various construction and flow parameters were analyzed, and their outcomes fit using simple linear equations to provide useful trends in how geometry affects heat exchange.
A total of 420 unique experimental test conditions were measured under the conditions designated in the SRM, with three replicates for each gas and flow orientation. Range-finding and experimental preparations yielded two additional complete replicates for the co-flow air condition, and the data were included with the three original replicates to create five distinct measurements of each condition. Runtime to model all 420 conditions using baseline model divisions (100) and subdivisions (100) was 3528.4 s or 8.4 s per condition. A stock clock rate Intel i7 5820K was used to run the model on a single thread with no added acceleration. The original goal of computationally cheap and accurate monolith modeling was thus proven, and the results are discussed in following passages.
4.1. Determination of Optimal Nusselt Number
With constant geometry, the value of
Nu should remain constant across different gases and flow orientations. As it is also critical to convective heat transfer, it is necessary to determine to model the heat exchange process. A binary search algorithm was used to determine an optimized
Nu value from the data. In his examination of arbitrary ducts, Shah [
16] established bounds of 2.0 ≤
Nu ≤ 4.0 for sinusoidal profiles. These bounds were used to constrain the search algorithm. The modeled reactor bulk outlet temperature was used to determine model validity, with the inherent stratification of the model weighted with each channel’s relative mass flow at the outlet to determine a mean value. The RMSE in the model outlet temperature was used as the algorithm loss function and was minimized using
Nu = 3.12 for the entire dataset. This value is within the range of values established by Shah and Cornejo et al. [
12,
13]. Fitting the model to experimental data using this value resulted in a near-linear fit.
Figure 11 shows these data as a scatter plot, with distinction between the flow pattern of the experiment and the inert gas used. Modeling trends show slight overprediction of the outlet temperature for co-flow and underprediction in counterflow. The slopes of both flow patterns follow that of the dashed prediction line, indicating that the error between the model and measurement is offset by a fixed value. Factors such as measurement uncertainty and day-to-day variability could induce such error. The thermocouples used to measure inlet and outlet temperatures are reversed between co- and counterflow conditions. Thermocouple aging, drift, and the asymmetrical insulation and flow characteristics between the two ends of the reactor could induce additional thermal effects, which the model does not capture. The co-flow outlet thermocouple is in indirect thermal communication with the cooled exhaust outlet, while the counterflow outlet thermocouple is in indirect thermal communication with the exhaust inlet. Minor radiative heat transfer from the heated inlet/outlet tubes would influence the thermocouple measurement in these cases, increasing the counterflow temperature measurement and decreasing the co-flow measurement under otherwise identical conditions. This is evident in the consistent overprediction and underprediction of the model, which increases at higher temperatures.
The minimized root mean square error was determined to be 20.16 K, which is acceptable considering the high flow rates and temperatures of the heat exchange assembly. For the case of reaction modeling, a change in temperature of 20 K would not change the reaction yield or activity under high throughput conditions. As shown in
Figure 12, the error was centered around an average of 0 K and the relative modeling error (%) was generally 5% or less. This is acceptable from a heat transfer standpoint, considering the high-speed solution offered with this approach. After design screening, further accuracy can be obtained through more accurate modeling methods like CFD.
4.2. Parametric Study
After
Nu calibration, the geometry of the modeled monolith was varied in a full-factorial analysis focused on the heat exchange effectiveness and integral heat transfer rate as responses. This model examined just the cylindrical inner monolith and assumed a constant elevated wall temperature to drive the heat exchange process. Input variables included the cell density of the monolith (CPSI), the flow rate through the monolith, the wall temperature (
Twall), and the aspect ratio of the monolith (
α). The monolith aspect ratio is defined with Equation (10). The flow rate was reported as a gas-hourly space velocity (GHSV), which is a typical metric for catalyst systems. Response data from the full factorial model were processed in JMP 15 ™ to determine the relationships between the four independent variables on the thermal responses. The simple model generated using JMP for both dependent variables is shown in the scatter plots of
Figure 13. Perfect linear agreement is plotted using a solid red line. Linear regression of a model output is an unconventional but useful approach in this analysis. By performing a regression, a simple and relatively accurate relationship between each variable to the response can be obtained. Fit Equations (11) and (12) are formulated such that each independent variable is normalized to vary from +1 to −1, making the slopes of all variables directly comparable. The relative importance of each variable can thus be discerned with the magnitude of its slope.
Heat exchanger effectiveness defined earlier with Equation (9) was fit using multiple linear regression, yielding Equation (11). Minimization of the sample RMSE was used to determine best fit. Effectiveness was found to vary as a function of the four independent variables. RMSE was 4.07%, with an R
2 of 0.92. Most important to the heat exchange effectiveness is an aspect ratio, α. Monolith geometries with smaller diameters and longer channels result in higher heat exchange effectiveness. GHSV has the second-highest impact, with low velocities leading to the highest efficiencies. Physically, increasing flow to a heat exchanger will decrease efficiency under normal circumstances due to decreased residence time. This assumes that the flow conditions within the heat exchanger do not drastically change the convective heat transfer through turbulent transition or other means. Wall temperature does not strongly affect the heat transfer efficiency due to the temperature differential being captured in the heat exchange effectiveness term. Minor impact results from an increased thermal differential near the outer radii of the cylinder. Due to radial conduction, which scales with the natural log of the radial distance, the apparent thermal resistance along the radius increases as the radial position approaches zero. The higher absolute temperature differential increases the overall penetration of heat due to the more complete heating of these outer radial positions. Cell density was shown to have almost no effect on the heat exchange effectiveness. In the three examined densities, flow was wholly laminar and modeled as such with the terminal Nusselt number from experimental fitting. Without physical modeling of the wall structures, the effects of catalyst wash coating, or capturing minute variation in individual cells of the monolith, it is impossible to determine if cell density more strongly affects heat exchange effectiveness or plays no role whatsoever. From this regression, it is clear, however, that its effect is minor.
Values of
hint were also modeled in JMP 15 ™ using a simple linear fit to analyze the individual contributions of each independent variable to the overall system. The fit is reported in Equation (12), with similar scaling coefficients to indicate relative importance to the final value of
hint. The best fit for the integral heat transfer coefficient was found as presented, with an R
2 of 0.98743 and an RMSE 0f 0.481 W m
−2K
−1. Like the heat exchange effectiveness, the cell density was relatively unimportant to the overall heat exchange rate. This makes sense, considering that some of the physical phenomena, which would change experimentally with cell density, such as the void fraction and conduction, were held constant in the model. These values were not varied as inputs to the model as they depend on the foil thickness chosen during monolith construction. As this variable could vary widely based on engineering requirements, the added complexity and expense of its inclusion would not yield comparable information useful for this study. Aspect ratio and temperature both show a similar scale of importance in affecting
hint. As shown previously by Boger and Heibel [
17], higher wall temperatures lead to higher log-mean temperature differential (LMTD) values in a monolith, which, in turn, create higher integral heat transfer coefficients. The same effect is shown in this modeling. Integral heat transfer rates are 1–2 orders of magnitude lower than those measured experimentally by Boger and Heibel. This is due to the larger diameter monolith used in this study, the lower-conductivity monolith material, and the convective condition at the monolith shell. This study used heated gas at the shell, whereas the previous study used circulating liquid water, which would provide at least an order of magnitude larger convective rate at the surface.
Two considerations must still be made regarding cell density. The pressure drop across the assembly is affected by cell density regardless of heat transfer effects. Higher cell densities have lower open cross-sectional areas, leading to higher restrictions. Higher cell densities also have higher surface area to volume ratios, leading to higher wash-coat loading and greater reactor intensification. The finding that cell density does not affect heat exchange effectiveness makes the design process simpler. The trade-off becomes one of catalytic performance required versus wash-coating cost under constant heat exchange capabilities.
Modeled values for h
int were already shown to vary as a function of GHSV, the aspect ratio, and the initial temperature differential (Δ
Tinit). This effect is best illustrated using continuous contour plots for various temperature differentials.
Figure 14 shows a series of these contour plots generated from the model, which further illustrates this relationship. Increasing Δ
Tinit consistently increases h
int regardless of other input parameters. At aspect ratios below unity, the slope of contours shows a decreasing trend. The lower limit of the aspect ratio, which approaches zero, implies an infinitely flat disc of an infinite diameter. In this case, it would be expected that thermal resistance would go to infinity. It follows that these curves will rapidly collapse upon each other as the aspect ratio approaches zero. It is intuitive that a flat disc heated from the rim is not conducive to heat transfer. Moving in the other direction, towards high aspect ratios, it is shown that the contour slopes become flat, reflecting the linear fit shown earlier. An increased aspect ratio increases the surface area to volume ratio, which decreases the thermal path length from the heated outer surface to the internal gases, thereby increasing the overall heat transfer rate. For a cylinder of constant volume whose geometry varies with the aspect ratio, the surface area to volume ratio is given with Equation (13). Aspect ratio in the numerator confirms this observation. The physical significance of the heat exchange coefficient is an indication of overall resistance to heat flow.
Heat exchanger effectiveness was calculated for the same conditions described above and arranged in four contour plots, as shown in
Figure 15. As can be seen, heat exchanger effectiveness varied from 21% to 81% and was highly dependent on input conditions. Low GHSV resulted in higher efficiency, which makes intuitive sense as gas was exposed to the heated wall for a longer duration of time. Increasing the aspect ratio was also shown to increase efficiency. As with integral heat transfer rates, the surface area to volume ratio increases with the aspect ratio and increases exposure of monolith gases to the heated wall. The diameter of the monolith decreases as well, leading to lower thermal resistance between the bulk of monolith gas and the surrounding wall. A larger initial temperature differential showed a marginal increase in efficiency, which was more pronounced when the aspect ratio was large. Because a low aspect ratio creates such a large thermal resistance between the wall and the bulk of flowing gases, changing the initial temperature differential here showed little effect and efficiency was low in all cases. The opposite occurs in high-aspect-ratio conditions, where a higher temperature differential drives greater heat exchange rates, leading to the bulk of gases nearly equalizing with wall temperature via the outlet of the monolith.
To summarize, while parametric investigation confirmed an intuitive relation between thermal inputs, flow rates, and their effect on heat exchange outcomes, it also revealed surprising relationships between some of the parameters. Cell density showed little effect on thermal performance of theoretical monoliths, while a varying aspect ratio showed a significant effect due to the inherent change in the surface area to volume ratios. This means that heat exchange effectiveness would be maximized at minimum GHSV and at high aspect ratios, whereas heat exchange rates are maximized with GHSV, the temperature differential, and the aspect ratio. Because of the competing effect of the gas flow rate on the two parameters, a compromise in the design must be made. For example, in a steam reformer, the monolith design could be optimized to provide high absolute conversion at a low flow, high molar conversion at a high flow, or some middle point between those two. The optimum operation condition depends on the end use, but any of these cases can be targeted through this modeling approach.
5. Conclusions
A computationally inexpensive, yet sufficiently accurate, heat transfer model was developed for inter-monolith heat transfer in a concentric configuration using a control volume approach and discretized solving. The model assumed radial symmetry and resolved the monolith along axial and radial dimensions. The model was validated using an experimental monolith assembly as a heat exchanger between heated air and inert gases. Heat exchange rates were calibrated by solving for a fully developed Nusselt number, which was determined through a binary seek algorithm. This algorithm was used to fit the model to experimentally measured values. Best fit was determined with minimization of the RMSE between the predicted outlet temperature of the model and the observed outlet temperature for all experiments using a single value for Nu. The Nu value found through this method was 3.12, which was within the range reported in prior literature and reasonable for generic internal flow.
After Nusselt number calibration, the geometry of the modeled monolith was varied in a full-factorial analysis focused on the integral heat transfer rate and heat exchange effectiveness responses. Relative effects of all input criteria were analyzed using JMP 15 ™ software and multiple linear regression of the model output was used to simply describe the individual effects of each on the response. Aspect ratio and GHSV were the most important factors in the net response. The initial temperature differential showed a small impact on heat exchanger effectiveness, with little variation between contour plots. However, the initial temperature differential was shown to be nearly as significant as the aspect ratio regarding the integral heat transfer rate, generating a significant shift in the contour plot’s upper bounds. Cell density of the monolith was shown to be unimportant for both output criteria.
Contour plots were generated using regression parameters to show a continuous response using the three significant input conditions. Optimum heat exchange efficiencies were achieved by minimizing gas flow rate GHSV and maximizing the aspect ratio and initial temperature differential. In contrast, optimum heat transfer rates occurred by maximizing gas flow rate GHSV, with all else held constant. Using a baseline condition, α = 1.0 and GHSV = 6000, this corresponded to a heat exchange efficiency increase of 43.2%, and a heat transfer increase of 44.8% at respective optimum conditions. These findings indicated that for a given reactor design, one of these two operation modes must be selected to maximize performance. Maximum conversion will occur where temperatures are the highest, under the high heat exchange condition. Maximum heat recovery will occur in the maximum heat exchange coefficient condition. While one mode or the other will offer a greater benefit to a specific system, increasing the aspect ratio of monoliths was shown to improve heat transfer rates and efficiency regardless of the desired operation mode.
The model presented here easily incorporates changes in monolith geometry, foil material, and process gases, making it well suited for design screening or algorithmic optimization of round catalyst prototypes. While only one geometry was experimentally examined, the full factorial modeling approach showed that rapid screening of monolith parameters reveals immediate design improvements. A future work could further increase accuracy and model robustness through incorporating radiation modeling or mass transfer effects. Nevertheless, the accuracy shown using basic heat transfer and conservation laws is sufficient for early design stage investigation. The modeling procedure of this current work is recommended for identification of promising geometries and flow conditions, which can be later examined thoroughly in computationally expensive CFD models or physical prototype construction.