## Abstract

The high packing fractions of dense granular flows make them an attractive option as a heat transfer fluid or thermal energy storage medium for high temperature applications. Previous works studying the heat transfer to dense flows have identified an increased thermal resistance adjacent to the heated surface as a limiting factor in the heat transfer to a discrete particle flow. While models exist to estimate the heat transfer to dense flows, no physics-based model describing the heat transfer in the near-wall layer is found; this is the focus of the present study. Discrete element method (DEM) simulations were used to examine the near-wall flow characteristics, identifying how parameters such as the near-wall packing fraction and number of particle-wall contacts may affect the heat transfer from the wall. A correlation to describe the effective thermal conductivity (ETC) of the wall-adjacent layer (with thickness of a particle radius) was derived based on parallel thermal resistances representing the heat transfer to particles in contact with the wall, particles not in contact with the wall, and void spaces. Empirical correlations based on DEM results were developed to estimate the near-wall packing fraction and number of particle-wall contacts. The contribution from radiation was also incorporated using a simple enclosure analysis. The ETC correlation was validated by incorporating it into dense flow models for chute flows and cylindrical flows and comparing with the experimental data for each.

## Introduction

The heat transfer to particulate flows has been studied for decades due to its various industrial applications, including drying and heating of particulates such as food grains and coal. More recently, the heat transfer characteristics of particle-based heat transfer fluids are of interest for concentrated solar power tower applications, as particles offer many intriguing characteristics over the current fluids. The present work studies the heat transfer behavior of internal gravity-driven dense granular flows.

Internal granular flows are characterized by shearing at the wall boundaries and are driven by particle collisions. Dense flows demonstrate high packing fractions of approximately 60%. Studies of the heat transfer in a static particle bed indicate that the packing fraction plays a primary role in the heat transfer capabilities of the bed, as the primary pathway for heat transfer is through the interstitial gas [1]. Since the interstitial gas typically exhibits low thermal conductivity compared with the particle material, the effective thermal conductivity of a bed decreases with decreasing packing fraction. Therefore, the high packing fraction of dense flows makes them a viable option as a heat transfer fluid.

While the heat transfer characteristics of static particle beds have been studied quite extensively both experimentally (for example, Refs. [2–4]) and theoretically (as summarized in Ref. [1]), the heat transfer to dense flows has been less widely studied. Sullivan and Sabersky [5] provided one of the pioneering works in the field with their experimental and theoretical analysis of the heat transfer to a vertical chute flow. Using a small rectangular cross section with a single heated wall, they examined the variation in the heat transfer measured for different particles types, particles sizes, and flowrates. In an attempt to understand the significance of the discrete nature of the flow, they compared their experimental results with an analytic solution for a continuum plug flow in their setup. Their experimental results did not agree with the results predicted by the continuum solution; instead, the granular flows demonstrated inferior heat transfer performance, which was attributed to the discrete nature of the flow and the discrete interaction of the particles with the heated wall. The analytic solution was modified to incorporate an increased thermal contact resistance at the heated wall, which accounted for reduced heat transfer due to a more structured arrangement of particles adjacent to a surface. The contact resistance was represented as a thin air gap a fraction of a particle diameter thick (0.085 in their case) and was assumed constant for all test cases. The resulting semi-empirical correlation relating the Nusselt number with the Péclet number showed better agreement with their experimental results. Natarajan and Hunt [6] expanded upon the work of Sullivan and Sabersky, studying longer heated lengths and higher flowrates. Their experimental results demonstrated similar trends to those observed by Sullivan and Sabersky for low velocity, high density flows; namely, an increasing Nusselt number with increasing Péclet number. At higher flowrates, however, Natarajan and Hunt observed a critical point after which the Nusselt number decreased and plateaued. The semi-empirical model developed by Sullivan and Sabersky captured the results observed for low velocity flows but was unable to capture the decrease observed at higher flowrates. The deviation from the model at higher flowrates was explained by recognizing that the density of the flow adjacent to the wall varies with flowrate, rendering the constant contact resistance assumption of Sullivan and Sabersky invalid when considering a wide range of flowrates.

Spelt et al. [7], Patton et al. [8], and Ahn [9] studied the heat transfer to inclined chute flows over a wide range of flowrates. Similar to Natarajan and Hunt, their results displayed a critical point after which the Nusselt number began to decrease with increasing Péclet number. Once again, the model developed by Sullivan and Sabersky captured the trends observed at lower flowrates but was unable to capture the decrease observed at higher flowrates. It is important to note that the best agreement with the model was found when using an air gap thickness of approximately 1/40th of a particle diameter. The thinner air gap was attributed to the difference in system configuration; gravity pushes particles toward the heated wall in an inclined chute, whereas gravity has no effect on the compactness of particles at the wall in a vertical chute. Ahn attempted to modify the model developed by Sullivan and Sabersky to incorporate the variation in the contact resistance with packing fraction; their model, however, was successful only for certain particle materials.

Denloye and Botterill [10] also studied the heat transfer to vertical chute flows, comparing their experimental results with the “packet renewal” model originally developed by Mickley and Fairbanks [11] for fluidized beds. The modified version of the model developed by Baskakov [12], which incorporates a contact resistance at the wall, was also compared with the experimental data. The thermal contact resistance was approximated using the near-wall thermal conductivity model developed by Yagi and Kunii [13] for static beds. The modified model showed better agreement with the experimental data, once again indicating the importance of the wall-adjacent layer on the overall heat transfer to the flow.

More recently, Morris et al. [14] sought to characterize the thermal contact resistance for inclined chute flows using the same physics applied in discrete element method (DEM) simulations. Examination of the heat transfer from a heated wall to a single particle indicated that the heat transfer decreases quickly when the particle loses contact with the wall, as the small layer of air between the wall and the particle acts as an insulation layer. This phenomenon was incorporated into their model to estimate the contact resistance using a particle-wall distribution function. Using DEM simulations to examine the particle-wall distribution function for different flow parameters and angles of incline, an empirical correlation to describe the Nusselt number in the wall-adjacent layer was developed in the form of a seventh-order polynomial. Incorporation of this model into a two-phase continuum model as a boundary condition showed reasonable agreement with results from DEM flow simulations [15]. The model has not yet been compared with the experimental results.

Our previous work [16] experimentally investigated the heat transfer to gravity-driven dense granular flows of ceramic particles in cylindrical tubes. The effect of flowrate and other flow parameters on the heat transfer to the flow at low operating temperatures was studied, along with the effect of temperature on the heat transfer for a single flowrate. Comparison of the experimental results with a continuum solution for the system studied emphasized the importance of the discrete nature of the flow, as the granular flows demonstrated inferior heat transfer performance. This work seeks to develop a model to predict the heat transfer characteristics observed experimentally and to gain a better understanding of the near-wall heat transfer behavior, examining how the contact resistance varies with different flow properties. DEM simulations were used to examine different flow characteristics in the wall-adjacent layer, assessing how near-wall particle structure may influence the heat transfer behavior. These observations were used to develop a physics-based model to predict the heat transfer through the wall-adjacent layer in the form of an effective thermal conductivity (ETC). The ETC correlation developed was incorporated into a generalized heat transfer model derived for the cylindrical flow system studied experimentally and compared with the experimental results. The ETC correlation was also incorporated into the model developed by Sullivan and Sabersky [5] and compared with the experimental results of both Sullivan and Sabersky and Natarajan and Hunt [6] to emphasize the applicability and versatility of the model developed.

## Modeling Approach

*j*th layer of the composite. Note that Eq. (1) models the flow as a plug flow; this is a common assumption employed for internal granular flows, as such flows display an essentially uniform velocity profile with a thin shear layer [17,18]. Assuming constant thermal properties in each layer, Eq. (1) was solved subjected to the following boundary conditions:

Here, $h$ represents the convective heat transfer coefficient of the flow, $Di$ represents the inner tube diameter, and $W$ is a new parameter defined to represent the ratio of the outer radii of the two layers. Note that $W$ depends solely on the system configuration; therefore, according to Eq. (7), the Nusselt number depends solely on the ratio of the bulk-to-wall-adjacent layer thermal conductivities for a given system configuration.

Examination of Eq. (7) for a single layer (i.e., when $kb=kwa$) indicates that the Nusselt number converges to a constant value of 8, as expected for a thermally developed plug flow subjected to a constant wall heat flux [19]. It is interesting to note that Eq. (7) suggests that the Nusselt number varies with the system configuration (i.e., $Di/dp$), increasing with an increase in the tube-to-particle diameter ratio. As $Di/dp\u2192\u221e\u2009i.e.,\u2009\u2009W\u21921$, the Nusselt number again approaches the continuum solution of 8. This is reasonable, as the relative size of the bulk layer ($W$) approaches a value of 1. This observation suggests that the relative size of the bulk layer should be considered when designing dense flow systems.

### Effective Thermal Conductivity of Wall-Adjacent Layer.

In order to approximate the Nusselt number using Eq. (7), an understanding of the effective thermal conductivities of the bulk and wall-adjacent layers is required. As summarized by Van Antwerpen et al. [1], numerous correlations exist in the literature to estimate the bulk effective thermal conductivity of a static particle bed, while fewer exist for estimating the effective thermal conductivity of the wall-adjacent layer. While the existing correlations vary in complexity, all correlations depend on three primary parameters: the packing fraction of the bed, the interstitial gas thermal conductivity, and the solid material thermal conductivity. As noted earlier, previous works modeling the heat transfer to dense granular flows relied on empirical constants and static bed correlations to approximate the increased thermal resistance at the wall. Morris et al. [14] noted, however, that static bed correlations, which are derived for densely packed beds, do not capture the intermittent nature of particle contacts in a flow. Our work, therefore, sought to characterize the effective thermal conductivities of the bulk and wall-adjacent layers for a flowing case.

Discrete element method simulations were used to examine the particle flow mechanics and their variation with flowrate, inferring how the observed characteristics may affect the heat transfer to the flow. A description of the DEM simulation methodology is beyond the scope of this work; however, a summary of the DEM simulation setup and validation studies performed can be found in Refs. [20–22]. The DEM simulations were used to calculate the packing fraction and instantaneous number of particle-wall contacts for various flowrates (a description of calculation procedure is found in Ref. [23]).

An example of the variation in the packing fraction observed in both the bulk and wall-adjacent layers is displayed in Fig. 2 for 300 *μ*m spherical particles in a 7.75 mm ID tube. The packing fraction of the bulk layer varies minimally from the static bed value of 62%, decreasing by less than 4% for the range of flowrates tested. Static bed values for the bulk layer effective thermal conductivity were therefore deemed a reasonable approximation for the flowing case. The packing fraction of the wall-adjacent layer displays a slightly greater variation with flowrate, decreasing by up to 14%. A more interesting observation comes from examining the variation in the number of particle-wall contacts. On average, 10–60% of the particles in the wall-adjacent layer are not in contact with the wall, depending on the flowrate, particle diameter, and tube diameter. This variation is depicted pictorially in Fig. 3 using DEM images of the particle arrangements in the wall-adjacent layer for different flowrates. The red particles in Fig. 3 represent particles in contact with the wall, while the blue particles represent particles not in contact with the wall, and the yellow represents void spaces where no particle is present. These images indicate that the number of particle-wall contacts decrease with increasing flowrate. The variation in packing fraction with flowrate is also highlighted by the increased amount of yellow space at faster flowrates (Fig. 3(b)).

The significance of this decrease in particle-wall contacts on the heat transfer through the wall-adjacent layer is demonstrated by Morris et al. [14]. Looking at a single particle, they noted that the heat transfer from a heated wall to the particle decreases significantly as soon as the particle loses contact with the wall, as the thin gas layer that forms between the wall and the particle acts as an insulation layer. Consequently, when estimating the heat transfer through the wall-adjacent layer, it is important to consider the heat transfer to both particles in contact with the wall and those not in contact with the wall. This work sought to incorporate these observations into a simple correlation to estimate the effective thermal conductivity of the wall-adjacent layer for a dense particle flow.

#### Conduction in the Wall-Adjacent Layer.

The work of Legawiec and Ziolkowski [24] was used as a guide in the derivation of an effective thermal conductivity for the wall-adjacent layer. They first characterized the heat transfer to a single particle, considering the heat transfer through both the interstitial gas and the particle, then calculated the total ETC of the wall-adjacent layer based on the number of particles present in the wall-adjacent layer. Since they were considering a static particle bed, they assumed that all particles were in contact with the wall and that the amount of void spaces where no particle was present was small enough that the heat transfer through these spaces was negligible.

*n*th heat transfer pathway, $kn$ represents the thermal conductivity of the

*n*th pathway, $An$ represents the area of the tube wall associated with the

*n*th pathway, the subscripts $wa$, $C$, $NC$, and $void$ denote the heat transfer pathway associated with the wall-adjacent layer, contact particles, no-contact particles, and void spaces, respectively, and $NC$ and $NNC$ represent the number of particles in contact with the wall and not in contact with the wall, respectively. The thermal conductivity of the void spaces is equivalent to the thermal conductivity of the interstitial gas $kg$. The “contact” and “no-contact” thermal conductivities are determined for a single particle. Similarly, the contact and no-contact areas represent the tube wall area considered in the calculation of the heat transfer to a single particle of each type; they are calculated based on the projected area of a single particle on the wall. The area associated with the wall-adjacent layer ($Awa$) represents the total tube wall area of the tube length under consideration $Awa=Avoid+NcAC+NNCANC$. Note that the model was derived for axial flow through a cylindrical tube to mimic our experimental studies. Equations (9)–(13) were rearranged to give the ETC of the wall-adjacent layer, yielding

Here, $\u2200C$ and $\u2200NC$ represent the volume of a single contact and no-contact particle intersected by the wall-adjacent layer (i.e., partial volume of a sphere), respectively, and $A$ represents the ratio of the tube-to-particle diameters $A=Di/dp$. Note that contact particles are assumed to have a point contact and the packing fraction is estimated from the DEM simulations. Since the correlation was derived to describe experiments using ceramic particles, heat transfer through the small contact areas is negligible. This assumption becomes invalid for high particle-to-gas thermal conductivity ratios (>1000) [1]. The volume and the wall projected area for a single particle of both groups (i.e., contact and no-contact) were determined based on geometry.

Here, $K$ represents the ratio of the solid-to-gas thermal conductivities $K=ks/kg$. The thermal conductivities calculated using Eq. (18) show minimal variation from the curved wall solution (<1.5%). The largest variation is observed for small tube-to-particle radii ratios, where the effect of curvature is greatest.

*a*from all system configurations ($a=0.22)$ was implemented in our correlations. Substituting Eqs. (15) and (18) into Eq. (14) and letting $a=0$ for $kC$ and $a=0.22$ for $kNC$, the ETC of the wall-adjacent layer can be written in its simplified form

#### Radiation in the Wall-Adjacent Layer.

Experimental and theoretical studies of the heat transfer through static particle beds indicate that the heat transfer increases with temperature due to the increased radiative heat transfer contribution [1,2,26]. Similar observations are expected for dense granular flows at high operating temperatures. While numerous works have sought to model the radiation within a static particle bed, a clear modeling approach has not been determined due to the complexity of radiation heat transfer and the random packing arrangement of particles in a bed. Furthermore, the radiation heat transfer in the wall-adjacent layer has received minimal attention. The present work therefore sought to develop a simple model for the radiation contribution in the wall-adjacent layer to provide a first-order approximation for the heat transfer to dense flows at high operating temperatures.

An enclosure analysis was used as the basis for the model development. The particle flow near the wall was classified into different particle layers, as defined in Fig. 5. The first layer (i.e., the wall-adjacent layer) has a thickness of a particle radius; the remaining two layers are each one particle diameter thick. Each layer of particles was considered a surface of the enclosure, with the surface areas of all particles in the layer summed to find the total surface area (as indicated by the different colors in Fig. 5). According to Pitso [27], the radiation from a single sphere in a randomly packed bed penetrates no more than 2–3 particle diameters. A similar assumption was employed in the present model; the enclosure was therefore assumed to consist of four layers: the wall and three particle layers.

*y*th particle layer, $Ty$ represents the average temperature of the

*y*th layer, and $\sigma $ represents the Stefan–Boltzmann constant. If $Tw\u2212Ty\u226a4Twy\xaf2$ (where $Twy\xaf$ represents the average temperature of $Tw$ and $Ty$), then $Tw4\u2212Ty4\u22484Twy\xaf3Tw\u2212Ty$. This simplification seems appropriate for the flows studied since the temperature drop across a few particle diameters is significantly less than the wall temperature. A linear temperature distribution was assumed across layers 2 and 3 (i.e., $T=\Delta T,T2\u2212T3=\Delta T),$ while the temperature gradient across layer 1 (i.e., the wall-adjacent layer) was assumed twice that of the other layers due to the increased thermal resistance adjacent to the wall (i.e., $T=\Delta T).$ These temperature gradient assumptions were guided by our experimental results, which showed a 2× to 3× temperature gradient near the wall [16]. Finally, assuming $\Delta T\u226aTw,\u2009Twy\xaf\u2248Tw$. Incorporating the above assumptions, Eq. (22) can be rewritten as

In order to calculate the radiation ETC, the view factors from the wall to each layer must be approximated. The ray-tracing algorithm of starccm+ was employed to examine the radiation heat transfer between the surfaces and extract the view factors to each surface. The centroid data from the DEM flows discussed in the Effective Thermal Conductivity of Wall-Adjacent Layer section were used to create a solid model of the particles in the near-wall region of the cylinder using solidworks. The surface of each particle was split according to the cylindrical surface “cutting” the particle at the boundary of each layer; the total surface area of each layer was then assigned to the sum of the individual surfaces in the associated layer. The process of splitting particle surfaces was a very time intensive and difficult process, as each particle had to be manually selected; consequently, ensuring that all particles were selected was a challenge due to the random packing. A small section of the near-wall layers was therefore used in the modeling process. As a result, there was a small amount of radiation “leakage” to the top, bottom, and side surfaces. Table 1 summarizes the system configuration modeled and the resulting view factor calculations. A slow flowrate was modeled to match the flowrates tested experimentally at high operating temperatures in our previous work [16]. The view factors may change with increasing flowrate due to the variation in the wall-adjacent layer packing fraction; however, due to the time-intensive nature of view factor calculation, further examination into the effect of flowrate on view factors was beyond the scope of the present study. Note that the leakage was incorporated into the primary view factors (proportionately to each of the layers based on the view factors) such that the view factors of the enclosure summed to one.

## Effective Thermal Conductivity Results

### Flow Parameters.

The effective thermal conductivity of the wall-adjacent layer due to conduction can be determined using Eqs. (19)–(21); however, an understanding of the wall-adjacent layer packing fraction and the number of particle-wall contacts is required. A series of DEM simulations were performed using starccm+ to study these parameters and how they varied with tube diameter, particle diameter, and flowrate. Six different configurations, using different combinations of tube and particle diameter, were tested at 15 different flowrates, ranging from creeping flow (*U* = 0.002 m/s) to the upper limit of choked flow (which ensures a dense flow and varies for each system configuration). Two different tube diameters were employed (7.75 and 10.92 mm ID), along with three different particle diameters (300, 450, and 600 *μ*m). These values were selected as they reflect the systems tested experimentally in our previous work. Future work may examine a broader range of parameters to verify the trends observed here. The centroids of all particles were tracked at each time-step, for a total of 1000 time steps, and used to calculate packing fraction and number of particle-wall contacts at a given instant; the parameter values reported are time-averaged values.

where $U$ represents the average flow velocity and $g$ represents the gravitational acceleration. The correlation agrees with the simulation results for all test cases to within 5%.

The correlation agrees with the simulation results to within 16%. Although this error may seem large, the error observed in $kwa$ as a result of error in both the packing fraction and the number of particle-wall contacts is less than 7%.

It is important to note that the simulations were performed using particle material properties representative of the ceramic zirconia–silica particles used in our experiments. We recognize that the behavior of granular flows depends on many structural and contact properties. We suspect that the general trends discussed above are true regardless of particle material; however, we are uncertain whether the relative magnitudes will change. Since no universal rheological model for granular flows exists, it is difficult to find an all-encompassing model. Future work may seek to test a wider range of system parameters.

### Effect of Flowrate and Temperature.

Incorporating the empirical correlations discussed in the Flow Parameters section into the effective thermal conductivity correlation (Eqs. (19)–(21)), the variation in the heat transfer at the wall with flowrate at low operating temperatures was examined (Fig. 8). As expected, the ETC of the wall-adjacent layer decreases with increasing flowrate due to the decreasing packing fraction and number of particle-wall contacts. The relative contribution of the packing fraction and the number of contacts to the decrease observed was examined by varying one parameter while keeping the other constant. Each case is plotted in Fig. 8, using the slow flowrate value for the parameter held constant. It is clear that the decrease in the number of particle-wall contacts is the driving force behind the decrease in the wall-adjacent layer ETC. This makes sense, as the number of particle-wall contacts decreases more significantly with flowrate than the packing fraction (see Figs. 6(a) and 7). Finally, the significance of including the heat transfer to both particles with and without contact is also examined in Fig. 8 by studying the ETC assuming all particles contact the wall (i.e., when $a=0$). The thermal conductivity is greater than that predicted with NC particles included due to the superior heat transfer exhibited by contact particles.

The effect of temperature on the heat transfer through the wall-adjacent layer is displayed in Fig. 9 for various particle diameters (assuming $Di/dp$ is constant). As expected, the heat transfer in the wall-adjacent layer increases with increasing temperature due to enhanced thermal properties and an increased radiation contribution at higher temperatures. Larger particle diameters exhibit superior heat transfer at elevated temperatures due to increased radiation contributions resulting from larger void spaces and therefore penetration distances. This is an encouraging result for the validity of the ETC correlation developed, as the same phenomenon was observed in studies of the ETC in the bulk of a static bed. A similar result is observed when considering the heat transfer to a dense flow, as seen on the right-hand axes (i.e., dashed lines) of Fig. 9. Here, the two-layer model (Eq. (7)) and the wall-adjacent layer ETC correlation (Eqs. (19) and (24)) were used to predict the heat transfer coefficients for a dense flow through a cylindrical tube for the same particle diameters. Note that the *y*-axis plots $hDi$, such that the variation due to tube diameter is excluded. This result suggests that larger particle diameters may be more beneficial for high temperature applications.

## Model Validation

The effective thermal conductivity correlation developed for the wall-adjacent layer was validated by implementing it into dense flow models and comparing with the experimental data. In particular, the experimental data from our low and high temperature studies in cylindrical tubes and data published in the literature for vertical chute flows were considered. The models were found to have good agreement with the experimental results, as discussed in the Experimental Results and Published Data sections.

### Experimental Results.

The heat transfer to the dense flows through cylindrical tubes studied in our previous work [16] was modeled using the two-layer Nusselt number correlation presented in Eq. (7). Our experimental studies were subdivided into two parts: investigation into the effect of flowrate and other system parameters, such as tube and particle diameter, on the heat transfer to the flow at low operating temperatures, and investigation into the effect of temperature on the heat transfer to the flow. Results from the low temperature studies using the 7.6 mm and 10.8 mm ID tubes with 270 *μ*m zirconia–silica particles (as presented in Fig. 8 of Ref. [16]) are compared with the two-layer model in Fig. 10. The two-layer model results were calculated using thermal properties at the mean flow temperature of the tube (*T* = 105 °C). The error bars and dashed lines denote the 95% confidence interval for the experimental data and two-layer model, respectively. Note that the radiation term in the wall-adjacent layer ETC was neglected for low-temperature results presented in Fig. 10. A similar comparison for the remaining three data sets displayed in Fig. 8 of Ref. [16] can be found in Ref. [25].

The experimental results, apart from a couple stray values, show good agreement with the two-layer model within the experimental uncertainty. The worst agreement is observed at the two extremes of the flowrates tested. The slow flowrates (large inverse Graetz numbers) demonstrate the largest experimental uncertainty due to small temperature differences across the flow; the fast flowrates at shorter heated lengths (smallest inverse Graetz numbers) approach the thermal entrance region, where the thermally developed solution plotted is no longer appropriate. Close examination of Fig. 10 indicates that the 10.8 mm ID tube yields slightly larger Nusselt numbers than the 7.6 mm ID tube (direct comparison can be seen in Fig. 7 of Ref. [16]). This variation may potentially be explained by the two-layer model, which suggests that the heat transfer decreases slightly with a decreasing tube-to-particle diameter ratio (decreasing $W$).

Results from the high temperature studies, again using the 7.6 mm and 10.8 mm ID tubes with 270 *μ*m zirconia–silica particles, are presented in Fig. 11, along with the corresponding values predicted by the two-layer model. Note that the results are plotted as a function of the wall temperature here. The effective thermal conductivity of both the bulk and wall-adjacent layers is required to estimate the heat transfer using the two-layer model. Ideally, the bulk layer ETC would be calculated using the mean temperature and the wall-adjacent layer ETC would be calculated using the wall temperature. The model, however, was developed to be used as a predictive tool; as such, knowledge of both the mean and wall temperatures is unlikely. A relationship between the two is not straight forward; therefore, a single temperature to approximate the ETC of both layers is required. Using a single temperature will lead to minor inaccuracies in the value calculated; implementation of the wall temperature resulted in smaller errors and was therefore selected as the reference temperature. From a design standpoint, the expected heat transfer to a flow at a desired wall operating temperature may be more valuable; for example, if designing based on the maximum operating temperature of a tube. Due to the smaller temperature differences present in the low temperature studies and absence of radiation, accurate selection of a reference temperature was less of a concern.

The two-layer model shows good agreement with the experimental data within experimental uncertainty, with the largest variation observed at the highest temperatures. The relative contribution from radiation is also displayed in Fig. 11 by plotting the model prediction without radiation included. As expected, the contribution from radiation is negligible at low temperatures and begins to grow at higher temperatures. For the systems tested, the inclusion of radiation increases the heat transfer coefficient by up to 29%.

### Published Data.

The ETC correlation developed for the wall-adjacent layer was further validated using experimental results published by Sullivan and Sabersky [5] and Natarajan and Hunt [6] for vertical chute flows. Sullivan and Sabersky and Natarajan and Hunt studied the heat transfer to chute flows (rectangular cross section) with a single heated wall at a constant temperature, each using different heated lengths. Both found good agreement between their experimental results and the single-resistance model developed by Sullivan and Sabersky for low flowrates.

Here, $BC$ and $BNC$ represent the correction terms for the contact and no-contact particles for a cylindrical system, respectively $BC=1+(1/A);BNC=1+(1.35/A)$, $Am$ represents a modified tube-to-particle diameter ratio, $P$ represents the perimeter of the cross section, and $\beta $ represents a shape-dependent correction factor. Note that the flat plate solution is applicable for chute flows $Di/dp\u2192\u221e$, which is captured by the correction terms evaluating to 1.

While the correlations developed in the present study for the packing fraction and number of particle-wall contacts are based on data obtained for flows through cylindrical tubes, they were used to estimate values for the chute flows under consideration by utilizing the modified tube-to-particle diameter ratio. This may not be the best parameter to compare the flow physics in a cylinder and chute; however, without any chute data, the modified tube-to-particle diameter ratio was used in an attempt to stay consistent. The hydraulic diameter was also discussed as a potential controlling parameter; however, incorporation of the hydraulic diameter yielded negative values for the number of NC particles in a chute. A wider range of DEM data for different cross-sectional shapes would help clarify the uncertainties present in extrapolating the empirical correlations to different system configurations. An understanding of the physics controlling the packing fraction and number of particle-wall contacts in the wall-adjacent layer would also help explain the trends observed.

The experimental results for the glass particles tested by Sullivan and Sabersky are plotted in Fig. 12, along with the Nusselt numbers predicted by their single-resistance model, using their value of $\chi =0.085$ and the value calculated using $kwa$. Note that the bulk ETC was calculated using the correlation developed by Zehner and Schlünder [28], which has been shown to provide reasonably accurate predictions over a wide range of solid–gas thermal conductivity ratios [1]. The Nusselt numbers calculated using $kwa$ show good agreement with the experimental data, and therefore, the empirical value of $\chi $ implemented by Sullivan and Sabersky. This suggests that for the small range of slow flowrates tested, it is reasonable to assume that the thermal contact resistance remains essentially constant.

Natarajan and Hunt expanded upon the work of Sullivan and Sabersky by testing a wider range of flowrates and a longer heated length. Their experimental results for the glass particles tested are displayed in Fig. 13. Once again, the Nusselt numbers predicted using a constant $\chi =0.085$ and a variable $\chi $ calculated using $kwa$ are also plotted. Both model predictions show good agreement with the experimental data for slow flowrates (*U* < 9.0 cm/s). At faster flowrates, however, a constant resistance does not capture the decrease in Nusselt number observed. Incorporating the wall-adjacent layer ETC correlation, which shows a decrease in ETC with increasing flowrate (corresponding to an increase in $\chi $), captures the decrease in Nusselt number observed experimentally quite well.

## Conclusion

This work sought to gain a better understanding of the near-wall heat transfer characteristics of internal dense granular flows and develop a physics-based model to predict the heat transfer to such a flow. Experimental results from our companion study [16] corroborated the observation of numerous studies in the literature that dense granular flows demonstrate inferior heat transfer compared with a continuum flow. It is well-accepted that the decrease in heat transfer observed in discrete flows results from an increased thermal resistance adjacent to the heated surface due to a more structured arrangement of particles. While previous works have modeled the reduced heat transfer to dense flows, they rely on empirical parameters or correlations derived for static beds to estimate the near-wall contact resistance. The present work sought to overcome this deficiency by deriving an ETC correlation for the wall-adjacent layer based on the fundamental physics, considering potential variation for a flowing bed.

Discrete element method simulations were employed to examine the near-wall flow physics and understand how parameters such as the near-wall packing fraction and number of particle-wall contacts may affect the heat transfer from the wall; both parameters were shown to decrease with increasing flowrate. Since the heat transfer to a single particle decreases drastically upon losing contact with the wall, it is important to consider heat transfer to both particles in contact with the wall and particles that are not. Using a parallel circuit of thermal resistances, an ETC correlation for the wall-adjacent layer was derived. The resulting correlation depends on the tube diameter, particle diameter, thermal conductivities of the solid and gas materials, the packing fraction of the wall-adjacent layer, the number of particles in contact with the wall, and the wall temperature. Empirical correlations were developed based on DEM simulation results to estimate the near-wall packing fraction and number of particle-wall contacts. The wall temperature allows the radiation contribution to be incorporated.

The ETC correlation was validated by incorporating it into dense flow models for different system configurations. A two-layer Nusselt number correlation was developed in this study to estimate the heat transfer to dense flows in vertical cylindrical tubes; the resulting correlation was compared with the experimental results from our companion study. The two-layer model showed good agreement with the experimental results for different tube diameters, particle diameters, flowrates, and temperatures. The wall-adjacent layer ETC correlation was also used to estimate the near-wall contact resistance described in the single-resistance model developed by Sullivan and Sabersky [5] for vertical chute flows. The modified model showed good agreement with the experimental results of both Sullivan and Sabersky and Natarajan and Hunt [6]. Unlike the original model, which assumed that the contact resistance did not vary with flowrate, the modified version was able to capture the decrease in heat transfer observed at faster flowrates.

The agreement between the model developed and the experimental results shows promise for the ETC correlation developed. In its current form, the ETC correlation is limited by an understanding of the variation in the near-wall packing fraction and the number of particle-wall contacts with system parameters; extrapolation of the empirical correlations outside of the regime of parameters for which they were developed presents an uncertainty currently and must be considered with care. An understanding of the mechanics that control the behavior of each parameter as well as a study into the variation across a wider range of system parameters would greatly expand the applicability of the ETC correlation derived. Ultimately, we anticipate our model to serve as a gateway to analyzing the overall heat transfer behavior of dense granular flows and to serve as a means to identify other parameters that can be modified to improve the overall performance of dense flow systems.

## Funding Data

Advanced Research Projects Agency—Energy (ARPA-E), U.S. Department of Energy (Grant No. DE-AR0000414; Funder ID: 10.13039/100000015).

## Nomenclature

*a*=fraction of a particle radius that “no-contact” particles are from wall

*A*=tube-to-particle diameter ratio

*A*=_{n}tube wall area of the

*n*th heat transfer pathway, m^{2}*B*=_{C}correction term for contact particles in cylindrical system

*B*=_{NC}correction term for “no-contact” particles in cylindrical system

*d*=_{p}particle diameter, mm

*D*=_{i}inner tube diameter, m

*dR*=infinitesimal thermal resistance, K/W

*F*=_{w−y}view factor from wall to

*y*th layer*g*=gravitational acceleration, m/s

^{2}*h*=heat transfer coefficient, W/m

^{2}-K*k*=thermal conductivity, W/m-K

*K*=ratio of solid-to-gas thermal conductivities

*k*_{particle}=thermal conductivity of single particle in wall-adjacent layer, W/m-K

*k*_{wa,rad}=thermal conductivity of wall-adjacent layer due to radiation, W/m-K

*L*=length of heated section, m

*n*=_{C}number of contact particles per unit wall area, m

^{–2}*n*_{C,}_{max}=Max. number of contact particles per unit wall area, m

^{–2}*n*=_{NC}number of “no-contact” particles per unit wall area, m

^{−2}*N*=_{C}number of particles in contact with the wall

*N*_{NC}=number of particles

*not*in contact with the wall- $Nukb=$ =
Nusselt number using bulk thermal conductivity

- $NuL\xaf=$ =
length-averaged Nusselt number

*P*=perimeter of flow cross section, m

- Pe
=_{L}Péclet number based on length

- $qwall\u2009\u2033$ =
wall heat flux, W/m

^{2} - $qwall,rad\u2033$ =
radiation heat flux leaving the wall, W/m

^{2} *R*=_{i}inner tube radius, m

*R*=_{n}thermal resistance of

*n*th heat transfer pathway, K/W*r*=radial position, m

*r*=_{p}particle radius, mm

*T*=temperature, °C or K

- $Twy\xaf\u2009$ =
average temperature of wall and

*y*th layer, K *U*=mean axial flow velocity, m/s

*W*=ratio of bulk-to-wall-adjacent layer radii

*z*=axial position, m

### Symbols

- $\alpha $ =
thermal diffusivity, m

^{2}/s - $\beta $ =
shape-dependent correction factor

- $\delta $ =
radius of bulk layer, m

- $\Delta T\u2009$ =
temperature difference, K

- $\sigma $ =
Stefan–Boltzmann constant, W/m

^{2}-K^{4} - $\varphi wa\u2009$ =
packing fraction of wall-adjacent layer

- $\Phi C$ =
intermediate parameter for contact particles

- $\Phi NC\u2009$ =
intermediate parameter for “no-contact” particles

- $\chi \u2009$ =
relative thickness of air gap adjacent to wall in Sullivan and Sabersky model

- $\u2200$ =
volume of single particle in wall-adjacent layer, m

^{3}