Abstract

Heat transfer in packed beds and their thermal response have been of great interest for scientists and engineers for the last several years, since they play a crucial role in determining design and operation of reactors. Heat transfer of a packed bed is characterised through lumped parameter, namely, effective thermal conductivity. In the present studies, experiments were performed to investigate the thermal conductivity of a packed bed in radial direction. The packed bed was formed using iron ore particles. To determine the effective thermal conductivity a new transient methodology is proposed. The results obtained were compared with the models proposed by ZBS and Kunii and Smith.

1. Introduction

Packed beds have extensive applications in industrial thermal systems: metallurgical process units, chemical reactors, thermal storage units, heat exchangers, combustors, and so forth. Heat transfer in packed beds and their thermal response have been of great interest for scientists and engineers. In this context, authors have been involved in development of a process in which the sensible heat of off-gas is being utilised to produce directly reduced iron from a packed bed of iron ore fines and coal fines. The heat is being transferred indirectly to the bed via reactor wall. Since the reduction of iron ore with coal is highly endothermic, the production rate of such a reactor will be primarily controlled by the amount of heat that can be transferred through the reactor wall and subsequently through the packed bed itself. In view of this, study of heat transfer from the reactor wall to the particle bed and that through the bed itself gains considerable relevance.

Yagi and Kunii [1] in their classical paper discuss all the possible mechanisms that occur during heat transfer through packed bed. It involves mechanisms such as thermal conduction through solids, heat transfer through the contact surface area, radiant heat transfer between surfaces of neighbouring particles and beyond through voids, and heat transfer through convection and mixing of fluid.

Yagi and Kunii [1] also showed that, in packed beds with no imposed fluid flow, significant portion of the heat flows through the fluid film near the contact point of two particles. In general, heat transfer through packed beds has been characterized by an effective thermal conductivity in the bulk of the bed and a heat transfer coefficient or an effective thermal conductivity at the wall. Subsequently, researchers have carried out experiments exploring variation of void fraction near the packed bed container wall [24], variation of thermal conductivity due the effect of radiation [5, 6], variation of transport properties within a bed under moving and stagnant fluid conditions [7], and so forth. Some of the studies are also being made with radial [816] and axial [1, 1724] heating.

The following salient points are deduced from these studies.(1)The variation in effective thermal conductivity and heat transfer coefficient reported by various researchers goes as high as 100%. Often these variations have been attributed to inherent nature of packed bed heat transfer.(2)A relatively less number of experiments have been reported at higher temperatures where radiation phenomena can contribute significantly to the heat transfer. It is important to note that various metallurgical processes such as extraction, melting, and heat treatment are carried out at higher temperatures where radiation plays a significant role.(3)These studies showed difference in heat transfer characteristics along vertical (gravitational) and horizontal directions. These variations can be attributed to differences in convective flows of fluid phase in the voids as well as in contact resistances arising out of bed weight along the vertical and horizontal directions.

Large varieties of methods have been used to measure the thermal conductivity of packed beds. Tsotsas and Martin [25] classified the methods to determine thermal conductivity of packed beds as steady state and transient methods [8, 10, 13, 16, 26], radial and axial methods [824], and absolute and comparative methods [7, 27]. In steady state methods, the required boundary conditions are imposed and temperature profile within the bed is allowed to come to the steady state. Heat flux is measured and thermal conductivity of the bed is obtained by solution of Fourier or Laplace equation. In transient methods, thermal response of the bed is recorded by giving a pulse input in the form of temperature change or heat flux change. The bed conductivity is obtained by fitting the solution of Fourier equation to the measured temperature-time curve. In absolute methods, the bed conductivity is calculated from experimental results without any additional information. In comparative methods, the heat transfer in the packed bed is in series with that through a material with known conductivity. The bed conductivity is obtained by using ratios of temperature gradients within the materials.

In the present work, a new transient experimental methodology is proposed. It should be noted that thermal conductivity along the radial direction in a cylindrical bed can only be determined using transient experiments, unless the bed is annular and a sink for heat is provided [19]. Heat transfer measurements in packed beds generally show large scatter in data. Many researchers [21, 2830] have explained this scatter as various possibilities of errors during the experiments of heat transfer in a packed bed. Authors opine that, in addition to the inherent uncertainties involved in packed bed experiments, possibilities also exist that some additional factors such as the stress distribution within the bed due its own weight as well as external load and the heat flow in the gravitational direction or in the opposite direction (changes in free convective contribution) also need to be considered. In this regard, authors have conducted some interesting experiments and those results will be communicated in their future correspondence.

Many metallurgical processes involving packed beds operate at high temperatures, that is, greater than 300°C, wherein radiative contribution to the overall heat transfer plays an important role. However, conducting experiments at higher temperatures is more difficult than that at lower temperatures and therefore a relatively less number of experiments have been reported at higher temperatures.

In this work, heat transfer characteristics of the bed are studied at higher temperature with 1D radial heating of a bed of iron ore particles. In these experiments, no heat sink is used. Temperature profile inside the bed has been obtained by inserting thermocouples at various radial positions. In order to determine the thermal conductivity one needs the rate of heat input and the temperature profile. In the present methodology, the heat input to the packed bed is determined from the transient temperature inside the bed and not from the electrical power supplied to the furnace. Effective thermal conductivity of the bed is determined using the temperature profiles obtained from the experiments. The values of the conductivity are compared with the two widely used models by Kunii and Smith [31] and Schlunder and coworkers [32, 33]. Authors have conducted experiments in both directions, that is, axial and radial. Experiments in both directions showed interesting facts, which will be communicated later. In the present work, discussion has been made only for the experiments in radial direction.

2. Experimental Setup and Procedure

Figure 1 depicts the experimental setup for the study of radial heat transfer under transient conditions with stagnant fluid. Experiments were carried out for the bed heights of 200 mm and 300 mm. The bed is formed with iron ore particles with the average size of 4.3 mm. As shown in Figure 1(a), a cylindrical electrical resistance furnace was designed and fabricated. In order to ensure a sufficiently long hot zone wherein the temperature is uniform, three independently controlled electrical heating coils are placed, one at the center and the other two on each end of the furnace. By adjusting the heat input to the three segments independently one can achieve a large constant temperature zone. A thin wall cylinder of stainless steel 316L containing the packed bed is inserted vertically within the furnace. The packed bed is placed almost entirely within the constant temperature zone. All the temperature measurements are done in the middle section of the packed bed on a horizontal plane. Thermocouples are placed along a diameter of the packed bed. Temperature in the vicinity of any particle and within it has been verified to be within experimental error of ±0.1%.

Cylinder wall temperatures are obtained by embedding thermocouples into the wall, by first drilling a small blind hole and then closing it after inserting the thermocouple tip. This is to ensure good contact and, specifically, to avoid direct radiation from the furnace. To avoid the fin effect, this thermocouple is led along the wall for some distance before being taken out. Similarly the thermocouples measuring bed temperatures are led axially, since axial variation of temperature is expected to be small due to the large constant temperature zone. To ensure proper positioning of the thermocouple tips, the wires coming out through the bottom flange are weighed down so that the wire remains taut within the bed. Figure 1(b) depicts the thermocouple positions and the respective distances. Instrumentation for the experiments is shown in Figure 2. Calibrated chromel-alumel -type thermocouples of gauge 0.25 mm have been used for temperature measurement. The temperature measuring accuracy was ±1 K. The sheath of the shielded thermocouple wire is properly grounded. Data from the thermocouples are recorded continuously using a high-gain data acquisition card (PCL 818HG, Advantech, USA) and a personal computer. The maximum sampling rate of the card is 100 kS/s. Thermocouples are connected to the data acquisition system card through isolation modules (ADAM 3011, Advantech, USA). The selected temperature range for the present work is 400 K to 900 K.

3. Analysis of the Experimental Data and Discussion

The experiments were conducted in a vertical cylindrical bed kept inside a vertical 3-segment tubular furnace, heated so as to give a large constant temperature zone. The packed bed is contained in a stainless steel tube, and the measurements are made in the radial direction at a midway cross-section. Care is taken to ensure 1D heat transfer in the radial direction. Figure 3 depicts the typical temperature evolution obtained for the bed using thermocouples numbers 1 to 9. Thermocouple 6 shows the lowest temperatures as it is at the centre of the bed.

Thermocouples numbers 1 and 2, 3 and 9, 4 and 8, and 5 and 7 are placed diametrically opposite to each other, distance from the center being almost the same; hence, their temperatures are similar. Thermocouples numbers 10, 11, and 12 are placed to verify that the heat flow is one-dimensional. Thermocouples numbers 4, 10, and 12 are at the same distance of 18.5 mm from the wall but at different axial and azimuthal positions. Similarly thermocouples 11 and 8 are placed 17.5 mm from the wall but have different azimuth angles. It was observed that the variations in temperatures in the axial as well as the azimuthal directions were as low as ±0.15%. Figure 4 depicts the typical radial temperature distributions as time progresses.

Experimentally obtained temperatures are shown using the “” symbol. The lines shown are best fit curves using a quadratic function. The present system is cylindrical in nature; considering the same, the authors also used Bessel function [34] of first and second kind to obtain a best fit. The variation in temperature obtained by quadratic function and Bessel function was observed to be ±0.25%. A schematic of the fitted curve is shown in Figures 5 and 6 for times and . The measured wall temperature is shown by and . The temperature at the wall obtained by extrapolation of the curves up to the wall is shown by and .

For the given system at unit length, heat content within the material at particular time is calculated by integration using the fitted temperature profile as is the best fit polynomial for temperature as a function of radius of a packed bed, is the radius of the packed bed under consideration, is the specific heat of the material, is the bulk density of the bed, and is the heat content within the material at time .

Thermocouples 3 and 9 are placed near the wall with a distance of less than 4 mm from the wall. Achenbach [35] observed that porosity near the wall is higher than that in the bulk (away from the wall). In radial direction, the variation of porosity near the wall levels out at 4 particle diameters. In the present experiments, average particle size is 4.3 mm. As thermocouples 3 and 9 are placed in the near wall region wherein the porosity variation exists, they should not be included in calculation of effective thermal conductivity in bulk. It should be noted that thermocouples numbers 4 and 8 are placed at distance of 22 mm from the wall, which is more than five particle diameters.

Therefore, in order to determine the effective thermal conductivity in the bulk away from the wall region, the measured data of thermocouples 6, 5, 7, 4, and 8 were used. Rate of heat transferred to the bed under consideration, thus, was determined using (1), by fitting a quadratic using temperatures obtained from these thermocouples. Further, the temperature gradients at locations of thermocouples 4 and 8 were obtained from the fitted quadratic equation. Subsequently, the following equation is used to calculate the thermal conductivity in the bulk of the bed along with (1): As discussed in previous section, Figure 7 represents the bed thermal conductivities measured with iron ore particles of average size 4.3 mm for bed heights of 300 mm and 200 mm. Temperature measurements have been done at the bed heights of 150 mm and 100 mm, respectively, below the top of the bed. The data presented correspond to several experiments conducted at varying heat input rates. Considering the variation one expects in packed bed experiments, reproducibility of the results is good. It is observed that effective thermal conductivity of the bed increases with temperature. Figure 7 depicts that the values obtained are in the range of 0.20–0.35 W/m·K. These values compare well with the values reported by several workers in the literature [4, 6, 7, 15, 16, 18, 25, 36]. The data is obtained for the temperature range of 400 to 900 K. The scatter obtained is more for the lower temperature compared to that at higher temperature. The scatter for the temperature range 400 to 600 K is 0.025, for 600 to 700 K is 0.0234, and for 700 to 900 K is 0.0164. It is to be noted that the scatter in effective thermal conductivity is much higher than the error in values arising from error in measurement of temperature as well as thermocouple positions. The conductivity data obtained is compared with two models given by Kunii and Smith [31], Zehner and Schlunder [32], and Bauer and Schluender [33] as depicted in Figure 7.

The model presented by Zehner, Bauer, and Schlunder is popularly known as ZBS model. Both the models are presented as follows as a reference.(1)Model proposed by Kunii and Smith where (2)ZBS model where Kunii and Smith developed an expression for effective thermal conductivity with stagnant fluid based on one-dimensional heat diffusion model for a unit cell of packed spheres. ZBS model considered the heat flux assuming parallel heat flux vectors as a unit cell. ZBS model also accounted for particle shape, radiation effect, fluid pressure dependence, contact conduction, particle flattening, shape and size distribution, and oxidation effects using adjustable parameters for particles. Particle shape factor is considered as 2.5. In the present study, the model is presented by omitting the surface oxide terms and low-pressure terms. Thermal conductivity of the single iron ore particle was measured and observed to be 0.95 to 2.8 for the temperature range of 300 to 1200 K. For the same temperature range, emissivity was considered as 0.81 to 0.71.

The values obtained in the present work were compared with those calculated using the reported models [3133]. They are depicted in Figure 7. It is interesting to note that the thermal conductivity values for packed beds of height 300 mm are higher than those for packed beds of 200 mm height. This can be attributed to the interparticle contact forces generated by the weight of bed. A detailed investigation on this aspect is being carried out and will be communicated in future.

4. Conclusions

A new experimental technique has been developed to measure the thermal conductivity of packed beds along the radial direction in a cylindrical packed bed. The technique involves heating a cylindrical packed bed from the periphery and measuring its transient response. The effective thermal conductivity values obtained were comparable to the values mentioned in the literature. The data obtained from experiments is comparable with ZBS model and the model by Kunii and Smith.

Nomenclature

:Deformation coefficient
:Particle shape factor
:Diameter of particle (m)
:Acceleration due to gravity (m/s2)
:Heat transfer coefficient (W/m2·K)
:Heat transfer coefficient of thermal radiation, solid surface to solid surface (W/m2·K)
:Heat transfer coefficient of thermal radiation, void space to void space (W/m2·K)
:Heat transfer coefficient in packed bed with fluid flowing (W/m2·K)
:Thermal conductivity of packed bed (W/m·K)
:Effective thermal conductivity of packed bed (W/m·K)
:Effective thermal conductivity of packed bed in motionless fluid (W/m·K)
:Effective thermal conductivity in the vicinity of the wall of packed bed (W/m·K)
:Thermal conductivity of stainless steel disk (W/m·K)
:Effective thermal conductivity of packed bed (W/m·K)
:Thermal conductivity of fluid (W/m·K)
:Radiative component of effective thermal conductivity (W/m·K)
:Rate of heat transfer per unit length (W/m)
:Radius of a packed bed (m)
:Measured temperature on the face of the stainless steel disk (K)
:Bed temperature at the interface obtained by extrapolation (K)
:Measure of the effective thickness of the fluid film adjacent to the contact point between two solid particles
:Void fraction
:Emissivity
:Stefan–Boltzmann constant (W/m2·K−4).

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.