The calculations of dispersion coefficients inside two‐dimensional randomly packed beds of circular particles

The longitudinal (D_L) and transverse (D_T) dispersion coefficients in two‐dimensional (2‐D) randomly packed beds of circular particles in a laminar flow regime are derived. A 2‐D discrete system of particles is divided into cells using modified Voronoi diagrams. The relationship between the variation of the stream function and the averaged vorticity is obtained from computational fluid dynamics (CFD) simulations. The whole flow pattern is then obtained by using the principle of energy dissipation rate minimization. The obtained values of D_L agree well with 3‐D experimental data for all velocities investigated. At very high velocities, D_T in 2‐D appears to be higher than 3‐D experimental data. In addition, the effects of particle‐size distributions, packing structure, and porosity on the D_L and D_T were studied. One result was that an increase in the width of the particle‐size distribution resulted in higher values of D_L and D_T at high velocities. © 2012 American Institute of Chemical Engineers AIChE J, 59: 1002–1011, 2013     

Near-wall porosity characteristics of fixed beds packed with wood chips

Packed beds of fuel wood chips are commonly found in thermal conversion processes such as combustion or gasification. Wood chips in particular are mostly used as fuel for small-scale domestic heating boilers but also for commercial-scale combustion units. The characterization of spatial voidage distribution inside the wood chip beds is of great importance for flow and reactor modelling. This study focuses on the radial porosity variations of cylindrical beds of three different types of commercially available wood chips including chips classified as G30 size class. The conventional technique of consolidating packed beds with a resin was chosen as the experimental procedure. The radial voidage distribution in different cylindrical beds is determined by image analysis of sections of the solidified packings. Additionally, a packing of monosized spheres was investigated in order to assess the selected procedure in comparison with widely available literature data for spheres. The results are discussed and summarized in a mathematical expression correlating the radial voidage distribution depending on average wood chip size, packing core porosity and dimensionless distance from the tube wall.     

Deposition of Brownian particles in packed beds

An analysis is presented to describe the deposition of Brownian particles from suspensions flowing through granular media. The constricted tube porous media model recently proposed by Payatakes et al.[8] was used to characterize the media and the deposition process was assumed to be controlled by convective diffusion as well as the nature of the surface interaction between the media grains and the Brownian particles. Three different tube configurations were considered and the differences among the results were found to be relatively unimportant. Agreement between the analysis and available experimental data was found to be reasonably good for the limiting case when the surface effect is favorable.     

Effective thermal conductivity of packed beds of silicon-copper particles

The effective thermal conductivity of a bed of silicon-copper particles, as used in the technical synthesis of silicones was measured with a non-stationary method. The particles have an irregular shape and an average diameter of 0·111 mm.From the measured temperature profiles and from literature it was concluded that a mechanism, in which the energy transfer occurs through a series path consisting of an effective solid-path length and gas-path length, plays the main part in the model for heat transfer suggested by Masamune and Smith.The values for the effective thermal conductivity were in the region 0·192–0·440 kcal m^?1 h^?1 deg^?1.     

Mixing of randomly moving particles in liquid—solid fluidized beds

The deterministic diffusion equation of the Fickian type has been widely employed in modeling solids mixing in fluidized beds. In this work, a stochastic or statistical model has been proposed for the solids mixing in a liquid—solid fluidized bed of relatively large particles. Specifically, the particle diffusivity under the steady fluidized state has been derived by combining the well-known ‘Ornstein-Uhlenbeck’ equation of stochastic processes and the so-called partition function of particle velocities. The values of particle diffusivity calculated from the model have been shown to agree well with the available experimental data. Furthermore, it has been shown that the solids-mixing process in the bed can be related to the time-dependent nature of particle motion.     

A new pseudo-continuous model for the fluid flow in packed beds

A new pseudo-continuous model, describing the fluid flow through packed bed reactors, is developed by formulating the Navier-Stokes equations for a statistically described domain geometry. A coupling term, occurring in the momentum equation, represents the fluid-solid interaction and is modelled by comparison with the exact Navier-Stokes equations. Static and dynamic, inviscid and viscous contributions are separated, and stringent requirements for the dependence on the local voidage are found which are met by a simple model. Good agreement with correlations and experiments is found for the predictions of the pressure drop and the radial distribution of the axial velocity.     

A note on the diffusional deposition of aerosol particles in packed beds

Aerosol penetration through packed beds of spheres under diffusional conditions is given by the empirical equation

$\text{C}{}_{\mathrm{L}}\text{C}{}_{\mathrm{L}}\text{=}\text{exp}\text{?6·39}\text{D}{}^{\mathrm{<mtext>2</mtext><mtext>3</mtext>}}\text{ν}{}^{\mathrm{<mtext>2</mtext><mtext>3</mtext>}}\text{R}{}^{\mathrm{<mtext>5</mtext><mtext>3</mtext>}}\text{L}$

recently proposed by Gebhart et al. (1973). The present note shows that this equation may be derived independently from Wilson and Geankoplis' (1966) correlation for mass transfer coefficients. It therefore follows that standard expressions for such coefficients may be used to predict aerosol collection under diffusional conditions.     

Numerical derivation of dispersion coefficients for flow through three‐dimensional randomly packed beds of monodisperse spheres

The longitudinal (D_L) and transverse (D_T) dispersion coefficients for flow through randomly packed beds of discrete monosized spherical particles are studied. The three‐dimensional (3‐D) porous‐medium model consists of thousands of spherical particles that are divided into cells using Voronoi diagrams. The relationship between the variation of the dual stream function and the vorticity between neighboring particles is derived using Laurent series. The whole flow pattern at low particle Reynolds number is then obtained by minimization of the dissipation rate of energy with respect to the dual stream function. The D_L is obtained by fitting the resulting effluent curve to a 1‐D solution of a continuous model. The D_T is obtained by fitting the numerical concentration profile to an approximate 2‐D solution. The derived D_L and D_T values are in agreement with 3‐D experimental data from the literature enabling a study of the effects of pore structure and porosity on D_L and D_T. © 2013 The Authors AIChE Journal published by Wiley Periodicals, Inc. on behalf of American Institute of Chemical Engineers AIChE J 60: 749–761, 2014     

Numerical simulations of transfer and transport properties inside packed beds of spherical particles

In this study, we investigate the transport and transfer properties inside packed beds of spherical particles by means of CFD simulations. Heat and mass transfer properties have been computed in packing configurations of increasing complexity at low to moderate Reynolds numbers (1<R_e<80). Only liquid-phase flows are studied (300<S_c<1000). The problem of contact points between particles, which is inherent to finite-volume numerical methods, is solved by applying a contraction of 2% on all the particles of the bed. We show that this treatment has very little influence on the results when analyzed with dimensionless numbers as N_u=f(R_e, P_r) or S_h=f(R_e, S_c). Finally, a very dense packing of spheres was built using a Discrete Element Method and used to represent the real granular media. Transfer predictions by the model are very realistic. Longitudinal and transverse dispersion coefficients are determined inside geometries containing hundreds of particles. Predictions of dispersion coefficients are consistent with literature, but a correction is applied to improve results, because the bed contraction leads to the underestimation of the transverse dispersion coefficient. The model is found to be very promising to study the “near wall channelling” phenomena inside small packed columns induced by the heterogeneity of the porosity profile close to the wall.     

Flow of generalized newtonian liquids through fixed beds of nonspherical particles

The flow of generalized Newtonian liquids through a random fixed bed of particles has been investigated and a universal method of calculation of the creeping flow was suggested. The usefulness of this method has been verified experimentally for the flow of power law, Ellis and Sutterby liquids through fixed beds of different nonspherical particles.     

Mass transfer and holdup in gas-liquid cocurrent flow packed beds containing water-repellent particles

Water-repellent particles were prepared by spraying polytetrafluoroethylene (PTFE) on activated carbon. Gas-liquid volumetric mass transfer coefficient and holdup were determined in gas-liquid cocurrent upflow and downflow packed beds from the measurements of gas desorption and volume, respectively. As the PTFE loading increased, the gas-liquid volumetric mass transfer coefficients in both upflow and downflow modes were enhanced. Gas holdup in the upflow mode increased with the PTFE loading, while the dynamic liquid holdup in the downflow mode decreased. The enhancement of the mass transfer rate from gas to liquid can be explained by the effect of water-repellency on the surface of activated carbon.     

Experimental investigation of pressure drop in packed beds of irregular shaped wood particles

The knowledge of the pressure drop across a packed bed of irregular shaped wood particles is of great importance for achieving optimal control and maximum efficiency in many applications, such as wood drying, pyrolysis, gasification and combustion. In this work the effect of porosity, average particle size and main particle orientation on the pressure drop in a packed bed is investigated. To this end, particle size distributions and porosities are determined experimentally.Based on the experimental results obtained in this study, the form coefficient C and the permeability K of the Forcheimer equation are calculated for different packed beds. The Ergun equation requires an average equivalent particle diameter that is derived from the measured particle size distribution. This equivalent diameter and the corresponding bed porosity are used in the well known Ergun equation in order to derive adapted shape factors A and B.Since a change in bed porosity and particle size, caused by the degradation of the wood particles and gravity, can be expected in a reacting packed bed, a set of shape factors for use with the Ergun equation is determined that are independent of porosity and particle diameter and fit the experimental data very well.     

A macroscopic model for shear-thinning flow in packed beds based on network modeling

The flow of non-Newtonian fluids in packed beds and other porous media is important in several applications such as polymer processing, filtration, and enhanced oil recovery. Expressions for flowrate versus pressure gradient are desirable for a-priori prediction and for substitution into continuum models. In this work, physically representative network models are used to model the flow of shear-thinning fluids, including power-law and Ellis fluids. The networks are used to investigate the effects of fluid rheology and bed morphology on flow.A simple macroscopic model is developed for the flow of power-law and Ellis fluids in packed beds using results from the network model. The model has the same general functionality as those developed using the popular bundle-of-tubes approach. The constant β, which appears in these models, is often directly derived from the tortuosity and a simple representation of the porous media. It is shown here that this can lead to incorrect and ambiguous values of the constant. Furthermore, the constant is a weak function of the shear-thinning index, indicating that no single bundle-of-tubes could ever properly model flow for a wide variety of shear-thinning fluids.The macroscopic model is compared to experimental data for shear-thinning fluids available in the literature. The model fits the data well when β is treated as an experimental parameter. The best-fit values of β vary, which is expected because even the constant C in the Blake-Kozeny equation varies depending on the source consulted. Additionally, physical effects, such as adsorption and filtration, as well as rheological effects such as viscoelasticity may affect the value of β. We believe that in the absence of these effects, β equals approximately 1.46 for packed beds of uniform spheres at relatively moderate values of the shear-thinning index (>0.3).     

Computational study of forced convective heat transfer in structured packed beds with spherical or ellipsoidal particles

Randomly packed bed reactors are widely used in chemical process industries, because of their low cost and ease of use compared to other packing methods. However, the pressure drops in such packed beds are usually much higher than those in other packed beds, and the overall heat transfer performances may be greatly lowered. In order to reduce the pressure drops and improve the overall heat transfer performances of packed beds, structured packed beds are considered to be promising choices. In this paper, the flow and heat transfer inside small pores of some novel structured packed beds are numerically studied, where the packed beds with ellipsoidal or non-uniform spherical particles are investigated for the first time and some new transport phenomena are obtained. Three-dimensional Navier–Stokes equations and RNG k–ε turbulence model with scalable wall function are adopted for present computations. The effects of packing form and particle shape are studied in detail and the flow and heat transfer performances in uniform and non-uniform packed beds are also compared with each other. Firstly, it is found that, with proper selection of packing form and particle shape, the pressure drops in structured packed beds can be greatly reduced and the overall heat transfer performances will be improved. The traditional correlations of flow and heat transfer extracted from randomly packings are found to overpredict the pressure drops and Nusselt number for all these structured packings, and new correlations of flow and heat transfer are obtained. Secondly, it is also revealed that, both the effects of packing form and particle shape are significant on the flow and heat transfer in structured packed beds. With the same particle shape (sphere), the overall heat transfer efficiency of simple cubic (SC) packing is the highest. With the same packing form, such as face center cubic (FCC) packing, the overall heat transfer performance of long ellipsoidal particle model is the best. Furthermore, with the same particle shape and packing form, such as body center cubic (BCC) packing with spheres, the overall heat transfer performance of uniform packing model is higher than that of non-uniform packing model. The models and results presented in this paper would be useful for the optimum design of packed bed reactors.     

Modification of Ergun equation for application in trickle bed reactors randomly packed with trilobe particles using computational fluid dynamics technique

Based on a slit model, a pellet scale model has been developed for calculation of drag force imposed on trilobe catalyst particles in a packed bed reactor. The drag coefficient for single gas phase flow in a porous media has been calculated by CFD simulation and the results compared to the Ergun equation. The results show that the drag coefficient predicted by Ergun equation should be modified for various bed porosities, particle aspect ratio and gas densities. Therefore, a correction factor has been proposed to correct the Ergun equation constants in various conditions for trilobe particles. Comparison between the proposed corrected Ergun equation results and experimental data indicates considerable agreement.     

Experimentally‐based constitutive relations for co‐current gas‐liquid flow in randomly packed beds

Experimental Observations on average pulse velocity and frequency in concurrent gas‐liquid (down) flow through randomly packed beds are used to extract constitutive relations for the gas‐liquid interaction and mean curvature terms that appear in a recently proposed volume‐averaged two‐fluid model for bubbly flow. The proposed closures lead to a reasonably quantitative prediction of the average pressure drop and liquid saturation under bubbly flow conditions and in the near pulse regime. In addition, the proposed closures provide realistic estimates for the location of the bubble‐to‐pulse transition in microgravity and in 1g down‐flow and predict the disappearance of the bubbly flow pattern at low liquid fluxes in 1g down‐flow. © 2016 American Institute of Chemical Engineers AIChE J, 63: 812–822, 2017     

Effective interfacial area in packed beds for absorption with chemical reaction

The extended crossflow model is proposed and developed for describing gas absorption with chemical reaction in a packed bed operated in the trickle flow regime. The model involves the division of the total interfacial area into the interfacial area of the static hold-up and that of the dynamic hold-up. The concept of the effectiveness of the interfacial area of the static hold-up is presented. A procedure is developed for calculating the effectiveness, and consequently, the effective interfacial area for various different absorption processes with chemical reaction, both liquid phase controlled and gas phase controlled. The predictions based on the present model agree well with reported experimental data. The present model is shown to be simpler and more accurate than previous models.     

Experimental and numerical investigation of structure and hydrodynamics in packed beds of spherical particles

In chemical industry, flows often occur in nontransparent equipment, for example in steel pipelines and vessels. Magnetic resonance imaging is a suitable approach to visualize the flow, which cannot be performed with classical optical techniques, and obtain quantitative data in such cases. It is therefore a unique tool to noninvasively study whole‐field porosity and velocity distributions in opaque single‐phase porous media flow. In this article, experimental results obtained with this technique, applied to the study of structure and hydrodynamics in packed beds of spherical particles, are shown and compared with detailed computational fluid dynamics simulations performed with an in‐house numerical code based on an immersed boundary method‐direct numerical simulation approach. Pressure drop and the radial profiles of porosity and axial velocity of the fluid for three packed beds of spheres with different sizes were evaluated, both experimentally and numerically, in order to compare the two approaches. © 2018 The Authors AIChE Journal published by Wiley Periodicals, Inc. on behalf of American Institute of Chemical Engineers AIChE J, 64: 1896–1907, 2018     

A novel coalescence model for binary collision of identical wet particles

A coalescence model that incorporates both capillary and viscous contributions of the liquid binder, along with the liquid bridge volume effects, has been developed. The objective is to predict the coalescence of two spherical particles coated with a thin liquid film approaching with equal initial velocities from opposite directions. The model is based on an overall coefficient of restitution that is determined with the aid of the approximate analytical values of the maximum possible energy dissipation and a critical value that depends on the total initial kinetic energy of particles. The maximum possible energy dissipation accounts for the energy loss due to the viscous and capillary effects and inelastic collision. The proposed simplified method to determine the critical velocity has been compared with the numerical solution of the general equation of motion, and an excellent agreement has been found. The coalescence model has been investigated at the limiting conditions by neglecting either the capillary effect or the viscous effect, respectively. Finally, comparisons have been made with experimental data, and a reasonable agreement has been found.     

Verification of a model for estimating the void fraction in a three-component randomly packed bed

The void fraction in a three-component randomly packed bed was calculated from the authors' model, and the calculated values were compared with published experimental data for spherical and irregularly shaped particles and with results from computer simulations. Results from the model were in good agreement with simulated and published experimental data.