Packed beds find extensive application in a wide variety of industries to cany out a large number of diverse processes. The main objective of the present work is to develop models to predict the arrangement of particles and based on them, to determine and evaluate the structural characteristics of packed beds. These problems have received only a limited attention in the literature. As a first attempt, spheres of uniform size are considered. Beds of aspect ratio up to 2 (referred to as low aspect ratio beds) are analyzed by application of principles of analytical geometry. Expressions are derived for the location of particles and for the structural characteristics of the beds, both of which show periodicity. This leads to the concept of a unit cell which is the repetitive section of the bed whose characteristics are the same as those of the complete bed. The beds fall into three distinct groups — those with aspect ratio between 1 and l√3⁄2, between 1√3⁄2 and 2, and with aspect ratio 2. Equations are distinct for each group. The aspect ratio shows marked influence on the structural characteristics of the beds. Agreement of the predictions on the overall void fraction with the available experimental data is excellent. Radial void fraction profiles are estimated by defining a concentric cylindrical channel (CCC) of an arbitrary thickness and with the cylindrical surface through the radial position of interest located at the middle of the CCC, and by accounting for the solid volumes of all the segments (in this CCC) of spheres with centers lying within a distance of a particle radius on either side of the cylindrical surface. The curved boundaries of the sphere segments are rigorously accounted for. The results show that the entire bed is filled with variations in the void fraction, starting from a value of unity at the wall and zero (or close to zero) towards the axis of the bed. Monte Carlo model for the simulation of high aspect ratio beds has not proved successful even with any of a wide variety of distribution functions for the coordinates of the sphere dropping point. With uniform distribution, the only distribution used in all the reports so far, and with normal distribution, there is not even a qualitative agreement with the reported data on void fraction variations. Distributions with asymmetric density functions such as exponential, Weibull, gamma and beta, show considerable improvement; beta distribution being the best. However even the best results with beta distribution show satisfactory agreement with the experimental data only up to about 2dp from the wall. Simulations with the cluster growth model, modified to account for the confining nature of the wall, lead to more satisfactory results. The proposed algorithm consists of building up the cluster, sphere by sphere, by calculating all possible interior and wall sites for placing an incoming sphere in a stable and non-overlapping position on the current cluster. A preference parameter is defined to place the new sphere at locations along the cross section of the column at which the experimental void fraction profiles show prominent minima, that is, locations around which the bed has relatively high solid volume. Void fraction profiles in beds of various aspect ratios simulated by this model show good agreement with the corresponding experimental data. The structural characteristics of the high aspect ratio beds thus simulated are evaluated. The number of spheres per unit length, Ni is correlated with the aspect ratio. It becomes proportional to the square of the aspect ratio, with the proportionality constant being close to 0.9, for aspect ratios greater than about 10. This follows since in these beds the overall void fraction becomes constant at 0.4. Majority of the spheres have contacts (with neighboring spheres) between 4 and 7, with the lower and upper limits for the coordination number being 2 and 9. The radial profile of the average coordination number (averaged over the height of the bed at the given radial position) shows small oscillations about a mean value of about 6 over almost the entire bed cross section starting from a distance of about ldp from the wall. At a distance of 0.5dp from the wall the predominant number of contacts is four while the mean value is about 4.3. The overall coordination number (averaged over the entire bed) shows inverse dependence on the aspect ratio. For random packings, that is, as the aspect ratio becomes infinity, the overall coordination number tends to six which corresponds to regular cubic arrangement. Cumulative number fraction, CNf is a global measure of the arrangement of spheres in beds of high aspect ratio. Its radial variation shows four distinct regions whose locations are independent of the aspect ratio The CNf values in each region are correlated with aspect ratio The correlations combined with that of NL lead to a very useful and effective model for predicting void fraction profiles in a bed of any specified aspect ratio The validity of the predictive model is demonstrated

Packed beds find extensive application in a wide variety of industries. The objective of the present work is to analyze and evaluate the effects of the wall on structural characteristics, hydrodynamics and heat transfer in packed beds of spheres. As a first attempt, spheres of uniform size are considered. The cylindrical wall of the bed confines the location of the particles thus leading to significant radial variations in void fraction and specific lateral surface area. The two characteristics at any given radial position r* are estimated by defining a concentric cylindrical channel (CCC) of an arbitrary thickness such that its boundaries are equidistant from the cylindrical surface passing through r* and accounting for the solid volumes or lateral surface areas of the segments of spheres (cap, slice, rod and annular ring) contained in the CCC and with centers lying within a distance of a particle radius from r*.The curved boundaries of the sphere segments are rigorously accounted for. The low aspect ratio beds (aspect ratio less than or equal to 2) show three distinct types of behavior. In beds of aspect ratio 2, the void fraction starts from a value of unity at the wall and decreases to a minimum and then increases to unity at the center of the bed. In beds with aspect ratio between l\/¯3/2, there is a continuous decrease in void fraction from unity at the wall to a fairly low value towards the axis and then a slight increase followed by another decrease. The profiles for aspect ratio less than l\/¯3/2 show a continuous decrease from a value of unity at the wall to zero towards the axis. In contrast, beds of high aspect ratio show heavily damped oscillations in the void fraction up to about five particle diameters from the wall and then a constant value. The lateral surface area variations in low aspect ratio beds show a steep fall from a very high value near the wall, and in high aspect ratio beds an oscillatory nature, though not as strong as in the corresponding void fraction profiles. The distribution of flow in packed beds for steady flow of an incompressible Newtonian fluid under isothermal conditions is modeled by using Ergun equation with Brinkman-type correction to account for the viscous effects in the region close to the wall. The confining effect of the wall is incorporated through the radial variations in void fraction and specific lateral surface area. The hydraulic radius in the region next to the wall is modified to take into account the resistance of the wall surface to flow. The resulting model equations with appropriate boundary conditions are solved numerically by collocation technique. The influence of aspect ratio in the range 1.25 to 20.3 and Reynolds number from 0.1 to 1000, the two most important factors affecting the flow behavior, is evaluated. The velocity profiles show a peak in the region close to the wall thus indicating severe channeling effect in this region. The magnitude and location of the peak depend on aspect ratio and Reynolds number. The model predictions agree remarkably with reported experimental data on velocity profiles in a bed of aspect ratio 10.7, and on the effect of Reynolds number on friction factors in beds of low aspect ratio. The radial variations in void fraction, velocity and effective thermal conductivity are incorporated in the two-dimensional pseudo-homogeneous steady-state model to analyze the wall effects on heat transfer in packed beds. Both constant wall temperature and constant wall flux boundary conditions are adopted. The equations are solved numerically using finite difference technique. The radial temperature profiles are seen to be fairly uniform in beds of low aspect ratio thus showing that the often made assumption of complete radial thermal mixing in low aspect ratio beds is valid. Beds of high aspect ratio show strong radial gradients. For constant heat flux condition the slope of the temperature profile remains constant after a small distance from the Inlet thus leading to thermally fully-developed flow. For this condition the heat transfer equations are solved analytically to obtain expressions for Nusselt number and the radial temperature profiles. There is a significant difference in the temperature profiles evaluated in the presence and absence of wall effects. Good agreement is found between the Nusselt numbers obtained from the model and reported experimental data.