Mass transfer coefficient for drying of moist particulate in a bubbling fluidized bed

Similar documents
Modeling and Control of Moisture Content in a Batch Fluidized Bed Dryer Using Tomographic Sensor

Modeling and Control of a Fluidised Bed Dryer

Inter-particle force and stress models for wet and dry particulate flow at the intermediate flow regime

Minimum fluidization velocity, bubble behaviour and pressure drop in fluidized beds with a range of particle sizes

IDENTIFICATION OF DEFLUIDIZATION REGION IN A GAS-SOLID FLUIDIZED BED USING A METHOD BASED ON PRESSURE FLUCTUATION MEASUREMENTS

TOPICAL PROBLEMS OF FLUID MECHANICS 17 ONE-DIMENSIONAL TEMPERATURE DISTRIBUTION OF CONDENSING ANNULAR FINS OF DIFFERENT PROFILES

Published in Powder Technology, 2005

cen29305_ch08.qxd 11/30/05 3:05 PM Page 451 INTERNAL FORCED CONVECTION CHAPTER 8 Liquid or gas flow through pipes or ducts is commonly used in heating

Travel Grouping of Evaporating Polydisperse Droplets in Oscillating Flow- Theoretical Analysis

Estimation of Hottest Spot Temperature in Power Transformer Windings with NDOF and DOF Cooling

Engineering Thermodynamics. Chapter 3. Energy Transport by Heat, Work and Mass

On Temporal Instability of Electrically Forced Axisymmetric Jets with Variable Applied Field and Nonzero Basic State Velocity

Investigation of Packing Effect on Mass Transfer Coefficient in a Single Drop Liquid Extraction Column

Modeling and Simulation of an Air-cooling Condenser under Transient Conditions

Section 8.5. z(t) = be ix(t). (8.5.1) Figure A pendulum. ż = ibẋe ix (8.5.2) (8.5.3) = ( bẋ 2 cos(x) bẍ sin(x)) + i( bẋ 2 sin(x) + bẍ cos(x)).

Simultaneous heat and mass transfer studies in drying ammonium chloride in a batch-fluidized bed dryer

Residence time distribution and dispersion of gas phase in a wet gas scrubbing system

Mixing Process of Binary Polymer Particles in Different Type of Mixers

Gas-side mass transfer coefficient of a laboratory column equipped with one sieve tray

Experimental and Numerical Investigation of Two- Phase Flow through Enlarging Singularity

Experimental Study of Temperature Variations in Fluidized Beds

Fluidization Regimes Characterized Using a Fast X- Ray Tomography Setup

2014 International Conference on Computer Science and Electronic Technology (ICCSET 2014)

Luis Manuel Santana Gallego 100 Investigation and simulation of the clock skew in modern integrated circuits. Clock Skew Model

A population balance approach for continuous fluidized bed dryers

The Use of COMSOL for Integrated Hydrological Modeling

Available online at Energy Procedia 100 (2009) (2008) GHGT-9

Throughflow Velocity Crossing the Dome of Erupting Bubbles in 2-D Fluidized Beds

Title: Fast X-Ray Tomography for the Quantification of the Bubbling-, Turbulent- and Fast Fluidization- Flow Regimes and Void Structures.

ME 354 MECHANICS OF MATERIALS LABORATORY STRESSES IN STRAIGHT AND CURVED BEAMS

P = ρ{ g a } + µ 2 V II. FLUID STATICS

Bureau International des Poids et Mesures

IMPLEMENTATION OF FULLY COUPLED HEAT AND MASS TRANSPORT MODEL TO DETERMINE TEMPERTATURE AND MOISTURE STATE AT ELEVATED TEMPERATURES.

Prediction of Minimum Fluidisation Velocity Using a CFD-PBM Coupled Model in an Industrial Gas Phase Polymerisation Reactor

A SIMPLIFIED APPROACH TO THE DRYING OF SOLIDS IN A BATCH FLUIDISED BED

COMPRESSIBILITY AND FLOW PROPERTIES OF A COHESIVE LIMESTONE POWDER IN A MEDIUM PRESSURE RANGE

The Slug Flow Behavior of Polyethylene Particles Polymerized by Ziegler-Natta and Metallocene Catalysts

MODIFYING GELDART CLASSIFICATION FOR VARIOUS COHESION FORCES

ERASMUS UNIVERSITY ROTTERDAM Information concerning the Entrance examination Mathematics level 2 for International Business Administration (IBA)

Modeling of dispersed phase by Lagrangian approach in Fluent

Fig.8-1 Scheme of the fluidization column

Heat Transfer Analysis of a Space Radiating Fin with Variable Thermal Conductivity

Polynomial Degree and Finite Differences

SIMILARITY METHODS IN ELASTO-PLASTIC BEAM BENDING

NATURAL CONVECTION IN HEAT EXCHANGERS - A NEW INCENTIVE FOR MORE COMPACT HEAT EXCHANGERS?

STUDY ON BUBBLE BEHAVIOR OF INERT GASES AT ENTRANCE NOZZLE IN SODIUM-COOLED FAST REACTOR

BALKANTRIB O5 5 th INTERNATIONAL CONFERENCE ON TRIBOLOGY JUNE Kragujevac, Serbia and Montenegro

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Psychrometrics

EXPERIMENTAL RESEARCH OF THE INFLUENCE OF PARTICLE SIZE AND FLUIDIZATION VELOCITY ON ZEOLITE DRYING IN A TWO-COMPONENT FLUIDIZED BED

Effect of Uniform Horizontal Magnetic Field on Thermal Instability in A Rotating Micropolar Fluid Saturating A Porous Medium

Principles of Food and Bioprocess Engineering (FS 231) Problems on Heat Transfer

Air and Heat Flow through Large Vertical Openings

Determinants of generalized binary band matrices

Experimental Study on Port to Channel Flow Distribution of Plate Heat Exchangers

Comments on A Time Delay Controller for Systems with Uncertain Dynamics

FLOW AND HEAT-TRANSFER MODELLING OF THREE-DIMENSIONAL JET IMPINGEMENT ON A CONCAVE SURFACE

Buckling Behavior of Long Symmetrically Laminated Plates Subjected to Shear and Linearly Varying Axial Edge Loads

Solving Systems of Linear Equations Symbolically

BASIC DESIGN EQUATIONS FOR MULTIPHASE REACTORS

PRESENTATION SLIDES: Hydrodynamic Scale-Up of Circulating Fluidized Beds

Exploring the relationship between a fluid container s geometry and when it will balance on edge

Design and Analysis of Silicon Resonant Accelerometer

FACTORS INFLUENCING MEASUREMENTS OF THE THERMAL CONDUCTIVITY BY HOT BALL METHOD

Section 2.1: Reduce Rational Expressions

Experimental Investigation on Segregation of Binary Mixture of Solids by Continuous Liquid Fluidization

Full Papers. Modular Simulation of Fluidized Bed Reactors. 1 Introduction. By Rouzbeh Jafari, Rahmat Sotudeh-Gharebagh, and Navid Mostoufi*

OPTIMIZATION OF A BRAYTON CYCLE AT MAXIMUM POWER OPERATING CONDITION AND AT MAXIMUM THERMAL EFFICIENCY OPERATING CONDITION

Part I.

UNIVERSITY OF BOLTON WESTERN INTERNATIONAL COLLEGE FZE. BEng (HONS) IN CIVIL ENGINEERING SEMESTER ONE EXAMINATION 2016/2017 GROUND AND WATER STUDIES 1

MODELING THE FIBER SLIPPAGE DURING PREFORMING OF WOVEN FABRICS

The Fundamentals of Duct System Design

Studies on flow through and around a porous permeable sphere: II. Heat Transfer

School of Business. Blank Page

AN EXPERIMENTAL STUDY OF THE FROST FORMATION ON A COLD SURFACE IN FREE CONVECTIVE FLOW

An Analytical Electrothermal Model of a 1-D Electrothermal MEMS Micromirror

MAGNETIC RESONANCE (MR) MEASUREMENTS OF THE MASS FLUX IN GAS-SOLID FLUIDIZED BEDS

Lecture 3: DESIGN CONSIDERATION OF DRIERS

VISUALISATION OF TWO-PHASE GAS-LIQUID PIPE FLOWS USING ELECTRICAL CAPACITANCE TOMOGRAPHY

Effect of column diameter on dynamics of gas-solid fluidized bed: A statistical approach

PARAMETER IDENTIFICATION, MODELING, AND SIMULATION OF A CART AND PENDULUM

Comparison of Transmission Loss Allocation Methods in Deregulated Power System

Experiment (4): Flow measurement

A NEW ANALYTICAL BRIDGE PIER SCOUR EQUATION

A comparative study between Electrical Capacitance Tomography and Time-resolved X-ray tomography.

CHAPTER 3 MODELLING AND ANALYSIS OF THE PACKED COLUMN

High-order approximations for unsteady-state diffusion and reaction in slab, cylinder and sphere catalyst

The Mean Version One way to write the One True Regression Line is: Equation 1 - The One True Line

PREDICTION OF MASS FLOW RATE AND PRESSURE DROP IN THE COOLANT CHANNEL OF THE TRIGA 2000 REACTOR CORE

Chem(bio) Week 2: Determining the Equilibrium Constant of Bromothymol Blue. Keywords: Equilibrium Constant, ph, indicator, spectroscopy

Simple Examples. Let s look at a few simple examples of OI analysis.

Modeling the bond of GFRP and concrete based on a damage evolution approach

EE2351 POWER SYSTEM OPERATION AND CONTROL UNIT I THE POWER SYSTEM AN OVERVIEW AND MODELLING PART A

INTERPRETATION OF THE SHEAR STRENGTH OF UNSATURATED SOILS IN UNDRAINED LOADING CONDITIONS

Applied Fluid Mechanics

BOUSSINESQ-TYPE MOMENTUM EQUATIONS SOLUTIONS FOR STEADY RAPIDLY VARIED FLOWS. Yebegaeshet T. Zerihun 1 and John D. Fenton 2

Applied Fluid Mechanics

Design Criteria for a Packed-Fluidized Bed: Homogeneous Fluidization of Geldart s Class B Solids

Suggested Answers to 2017 O Level Physics 5059 Paper 2

Available online at ScienceDirect. Procedia Engineering 121 (2015 ) 19 26

Liquid Feed Injection in a High Density Riser

Transcription:

Mass transfer coefficient for drying of moist particulate in a uling fluidized ed Dr. Yassir T. Makkawi 1, and Prof. Raffaella Ocone 2 1 Chemical Engineering & Applied Chemistry, Aston University, Birmingham B4 7ET, UK 2 Chemical Engineering, Heriot-Watt University, Edinurgh EH14 4AS, UK Astract Experiments on drying of moist particles y amient air were carried out to measure the mass transfer coefficient in a uling fluidized ed. Fine glass eads of mean diameter 125 m were used as the ed material. Throughout the drying process, the dynamic material distriution was recorded y Electrical Capacitance Tomography (ECT) and the exit air condition was recorded y a temp/humidity proe. The ECT data were used to otain qualitative and quantitative information on the ule characteristics. The exit air moisture content was used to determine the water content in the ed. The measured overall mass transfer coefficient was in the range of 0.0145-0.021 m/s. A simple model ased on the availale correlations for ule-cloud and cloud-dense interchange (two-region model) was used to predict the overall mass transfer coefficient. Comparison etween the measured and predicted mass transfer coefficient have shown reasonale agreement with the measurements. The results were also used to determine the relative importance of the two transfer regions. Keywords: Buling fluidized ed; Electrical Capacitance Tomography; Mass Transfer Coefficient; Particulate drying Corresponding author: Tel +44 0)121 204 3398, E-mail: y.makkawi@aston.ac.uk 1

1. Introduction Buling fluidized ed technology is one of the most effective mean for interaction etween solid and gas flow, mainly due to its good mixing and high heat and mass transfer rate. When applied to drying of non-porurs wetted solid particles, the water is drawn-off driven y the difference in water concentration etween the ule phase and the surrounding phases. Different mechanisms can regulate this process, depending on the ed hydrodynamics (i.e. ules characteristics) and the water content in the ed. Therefore, the design of uling fluidized ed dryer requires an understanding of the comined complexity in hydrodynamics and mass transfer mechanism. Mass transfer in fluidized eds is generally known to e more complex than heat transfer. In a gas-solid fluidized ed drier there are different phases which all contriute to the removal of moisture from the wet particles. These are the ule phase, its surrounding cloud phase and the dense annular phase. Thus, in experimental determination of the overall mass transfer coefficient, one must have in hand detailed quantitative data on the ule characteristics, as well as the drying rate. The most widely used mass transfer model of Kunii and Levenspiel [1] expresses the overall mass transfer in a uling ed in terms of the cloud-ule interchange and dense-cloud interchange. The first mechanism is assumed to arise from the contriution of ule throughflow (i.e. circulating gas from the cloud phase and in and out of the ule) in addition to the diffusion from a thin cloud layer into the ule. The second mechanism is assumed to arise only from diffusion etween the dense phase and the cloud oundary. For particles of aout 500 m, some researchers assume that the transfer is of a purely diffusional nature, and thus neglect the contriution of ule throughflow. However, Walker [2] and Sit and Grace [3] pointed out that, pure diffusional model may significantly underestimate, the overall mass transfer coefficient. Kunii and Levenspiel [1] reported that the true overall mass transfer coefficient may fall closer to either of the acting mechanisms depending on the operating conditions (particle size, gas velocity, etc.). They suggested accounting for the first mechanism y summing the diffusional and throughflow, and adding those to the second mechanism in a similar fashion as for additive resistance. Because of the growing interest on modelling as a tool for effective research and design, researchers on uling fluidized ed drying or mass transfer in general are nowadays seeking to validate or develop new mass transfer coefficient equations required for accurate prediction of the process kinetics. Ciesielczyk and Iwanowski [4] presented a fluidized ed drying model ased on the interphase mass transfer coefficient. They incorporated Mori and Wen [5] model 2

for the ule diameter into the mass transfer coefficient formulation of Kunii and Levenspiel [1] to predict a generalized drying curve for solid particles in a uling ed. The results have shown good agreement with the experimental measurement for group B particles of Geldart classification. Kerkhof [6] discussed some modeling aspect of atch fluidized ed drying, where Kunii and Levenspiel [1] mass transfer coefficient comined with particles internal diffusion model was used to simulate the thermal degradation of life-science products. The study suggested that the exchange etween the ule phase and the particles inside the ule play an important role in the drying process. Recently, Scala [7] experimentally studied the mass transfer around a freely active particle in a dense fluidized ed of inert particles. He concluded that the active particle only reside in the dense phase and never enters the ule phase, hence it has no direct contriution to the ule-dense phase interchanges. This contradicting oservation confirms that despite the considerale effort on developing fluidized ed mass transfer coefficients, there still remain uncertainties with respect to the assumptions used in developing these coefficients. The aim of this study was to measure the overall mass transfer coefficient in a conventional uling fluidized ed dryer. The measurement was used for the validation of a simple predictive model ased on the availale correlations and also to assess the various assumptions used in developing these correlations. For this purpose, the drying rate and the corresponding mass alance were calculated from the measurement of the drying air at outlet conditions. Information on the ed hydrodynamics were otained from Electrical Capacitance Tomography (ECT) images. 2. Experimental The primary ojective of the experiment was to measure the mass transfer coefficient for a drying process in a conventional uling fluidized ed. This required detailed knowledge of the fluidized ed hydrodynamics and drying rate. For this purpose, non-porous wet solid particles of glass eads were contained in a vertical column and fluidized using air at amient temperature. The fluidising air was virtually dry and otained from a high-pressure compressor. An advanced imaging ECT sensor was used to provide dynamic information on the fluidized ed material distriution. The sensor was connected to a data acquisition unit and a computer. The air outlet temperature and its relative humidity were recorded using a temperature/humidity proe. Since the air condition at the inlet of the fluidization column was constant and completely independent 3

of the uling ed operating conditions, only one proe was installed at the freeoard (air exit). The detailed experimental set-up is shown in Figure 1. 2.1. Equipment and materials The fluidization was carried out in a cast acrylic column, 13.8 cm diameter and 150 cm high. The column was transparent, thus allowing for direct visual oservation. A PVC perforated gas distriutor with a total of 150 holes (~1.8% free area), was placed 24 cm aove the column ase. The upstream piping was fitted with pressure regulator, moisture trap, valve and three parallel rotamaters. A one-step valve was connected efore the moisture trap and was used as the upstream main flow controller. The particles used were allotini (non-porous glass eads) with a mean diameter of 125 m and a density of 2500 kg/m 3 (Geldart A/B mixture). The detailed physical properties of the particles are given in Ta. 1. Distilled water at amient condition was used to wet the particles. A variale speed granule shaker was utilised to produce the final wetted mixture. The Electrical Capacitance Tomography imaging system used (ECT from Process Tomography Limited, Manchester, UK), consisted of two adjacent sensor rings each containing 8 electrode of 3.8 cm length. All electrodes were connected to the computer through a data acquisition module (DAM 2000). The PC was equipped with custom communication hardware and software (ECT32) that allow for online and off-line dynamic image display. The system is capale of taking crosssectional images of the ed at two adjacent levels simultaneously at 100 frames per second. Further details aout the ECT system used in this study and its application to fluidization analysis can e found in Makkawi et al. [8] and Makkawi and Ocone [9]. The exit air quality was measured using a temperature and humidity proe (Type: Vaisala HMI 31, Vantaa, Finland, measuring range: 0-100% RH, -40-115 C ). The proe was hung y a connecting wire inside the fluidized ed freeoard approximately 10 cm aove the maximum expanded ed height. 2.2. Procedure The experimental procedure employed was completely non-intrusive. This is descried in the following steps in the order of their occurrence: (a) A total weight of 4.5 kg of a dry allotini mixture was placed in a granule shaker after eing wetted y distilled water. The shaker was firmly clamped and operated 4

continuously for at least 25 minutes to ensure an even distriution of water content. Distilled water was used to eliminate any possile interference with the ECT signal (ECT works for non-conducting materials only) () The wetted particles were then loaded into the fluidization column. Prior to commencement of drying, the ECT sensor was calirated for two extreme cases. This was carried out y sliding the ECT sensor up to the freeoard to calirate for the empty ed case, and down to the static ed area to calirate for the packed ed case. It should e mentioned that, ecause the water content was limited to a maximum of 45 ml (1% moisture on dry solid weight asis), the possile changes in the particle/air permittivity during the drying process would e negligile. Further details on the sensitivity of the ECT system to moisture content in the fluidized ed dryer can e found in Chaplin and Pugsley [10] and Chaplin et al. [11]. (c) The wet ed material was fluidized at the required air flow rate. This was carefully adjusted to ensure the ed operation at the single ule regime. The temperature and relative humidity were recorded at 2 minutes intervals. Simultaneously, and at the 5 minutes intervals, a segment of 60 seconds ECT data were recorded. At the same time, the expanded ed height during fluidization was otained from visual oservations. (d) Finally, the drying rate was otained from the measured air flow rate and temperature/humidity data at inlet and outlet using psychometric charts and mass alance calculations. The recorded ECT data were further processed off-line and loaded into in-house developed MATLAB algorithm to estimate the ule characteristics. The aove descried procedure was repeated for the three different operating conditions summarized in Ta. 2. To ensure data reproduciility, each operating condition was repeated three times, making a total of nine experiment tests. 2.3. Measurement of mass transfer coefficient Considering a section of the ed as shown in Figure 2, the overall mass transfer coefficient etween the ule phase and the surrounding dense phase, k d, can e defined y the following rate equation: dc S u kd ( Cd C) dz V (1) 5

where C is the water concentration in the ule phase, C d is the concentration in the surrounding dense phase, u, S and V are characteristic features of the ule representing the rising velocity, the interphase area and the volume, respectively. For moisture-free inlet air, Eq. 1 is suject to the following oundary conditions: C ( ) 0 at 0 C in air z and C out C at z H (2) out The ule moisture content at the outlet C can e given y: (drying rate) C out = (ule mass flow rate) m air ( Cout Cin) m air (3) where m air and m are the mass flow rate of the fluidising air and ules respectively. Because of the assumption that the ules rise much faster than the gas through the dense phase and the inlet air was virtually dry, Eq. 3 reduces to: C out ( Cout ) air where surface. (4) ( C ) is otained from the measured temperature and humidity at the uling ed out air For a spherical ule, S V ratio appearing in Eq. 1 reduces to d 6, where d is the ule diameter. It should e mentioned that for a perforated distriutor (such as the one used in this experiment), coalescence of ules mainly takes place at a few centimetres aove the distriutor, therefore, the entrance effects are neglected and the ule characteristics are assumed to e independent of height (this was confirmed from the ECT images). Finally, assuming that the water concentration in the dense phase is uniform and remains unchanged during the ule rise ( C w w ) and integrating Eq. 1 from z 0 to z H, d water the mass transfer coefficient is otained as follows: d u Cd C out kd ln 6H Cd where respectively. ed d, u and H are the ule diameter, the ule velocity and the expanded ed height, (5) 6

2.4. Measurement of ule characteristics Experimental determination of the overall mass transfer requires knowledge of the ule diameter and velocity (see Eq. 5). Using the ECT, the size and velocity of the ules (or voids ) in a gas-solid fluidized ed can e otained. The distinct lowering of the solid fraction when the ule passes across the sensor area allows for identification of the ule events in a given time and space. The ule velocity was then calculated from the delay time determined from a detailed analysis of the signal produced y the two adjacent sensors, such that: u t where t t2 t 1, t 1 and t 2 represent the time when the ule peak passes through the lower and upper level sensors respectively, and represents the distance etween the centre of the two sensors, which is 3.8 cm. The method is demonstrated for a typical ECT data in Figure 3. (6) The ule diameter was otained from the ECT data of relative solid fraction at the moment of ule peak across the sensor cross-section. From this, the ed voidage fraction (the fraction occupied y ules) was calculated as follows: d D (7) is the ule fraction, P is the relative solid fraction (i.e. packed ed: P 1; where 1 P empty ed: P 0 ) and D is the ed/column diameter. This procedure is demonstrated for a typical ECT data in Figure 4. Further details on the application of twin-plane ECT for the measurements of ule characteristics in fluidized ed can e found in Makkawi and Wright [12]. 3. Theoretical Models 3.1. Estimation of mass transfer coefficient In a uling fluidized ed, it is widely elieved that mass transfer occurs in two distinct regions: at the dense-cloud interface, k dc, and at the cloud-ule interface, k c. The overall mass transfer may e dense-cloud controlled; cloud-ule controlled or equally controlled y the two mechanisms depending on the operating conditions. The following theoretical formulations of these acting mechanisms are mainly ased on the following assumptions: i. The fluidized ed operates at single ule regime. 7

ii. The ules are spherically shaped. iii. The ule rise velocity is fast ( u 5 ). U mf iv. The contriution of particles presence within the ule is negligile. v. The contriution of the gas flow through the dense phase is assumed to e negligile. mf As evident from the tomographic analysis of the ule characteristics, assumptions (i)-(iv) are to a great extent a good representation of the actual uling ehaviour. Cloud-ule interchange: According to Davidson and Harrison [13], the mass transfer at the ule-cloud interface arises from two different contriutions: 1. Diffusion across a limited thin layer where the mass transfer coefficient is estimated y: 0.975 g d 0.25 0.5 kc D (8) 2. Convection contriution as a result of ule throughflow, which consists of circulating gas etween the ule and the cloud (leaving the ules from top and re-entering it from the ase as shown in Figure 5. Based on the analysis of Davidson and Harrison [13], Murray [14] reported the following corrected equation for the contriution of ule throughflow: k 0. 25 (9) q U mf Kunii and Levenspiel [1] recommended adding oth acting mechanisms, thus the mass transfer coefficient etween the ule-cloud is given y: 0.25 0.5 g kc 0.975D 0. 25U mf d (10) Dense-cloud interchange: Utilizing the Higie penetration model [15], the mass transfer coefficient at the dense-cloud interface is expressed in terms of the ule-cloude exposure time and the effective diffusivity as follows [16]: k dc 4De t mf 0.5 (11) 8

where d u c t is the exposure time etween the ule and the cloud. From the analysis of tomographic images of the ules and its oundaries, it is oserved that the cloud diameter (the outer ring in Figure 5) is in the range of 1.2-1.8 ule diameters, and that the concentration of the particles within the core of the rising ule is virtually negligile. Therefore, in the current analysis we assumed that d ~ 1. 5d and that the contriution of the particles circulation etween the ule and the surrounding phases to the overall mass transfer process is negligile. The effective diffusivity, D is assumed equal to the molecular diffusivity, D. With these approximations, Eq. 11 reduces to: e c k dc D mf u 0.92 d 0.5 (12) Overall ed exchange: For an equally significant contriution from cloud-ule and dense-cloud interchanges, Kunii and Levenspiel [1] suggested adding oth contriutions in a analogy to parallel resistances, such that the overall ed mass transfer coefficient ( k d ) is given y: 1 k d 1 1 (13) k k dc c Sustituting Eq. 10 and 12 into Eq. 13 yields the overall mass transfer coefficient as follows: k A B 1 1 d (14) 0.5 d A1 B1 where A 0.25 0. 5 d g 0. U d 0.5 1 0.975 25 mf D (15) D 0. 5 B1. 92 mf u (16) 0 4. Results and discussions 4.1. Hydrodynamics Figure 6 shows the measured ule velocity and ule diameter as a function of the water content in the ed. These measurements were taken at different time intervals during the drying process. Each data point represents the average over 60 seconds. Both parameters vary slightly 9

within a limited range. These hydrodynamic oservations suggest that the ule characteristics almost remain independent of the water content, at least within the range of operation conditions considered here. This is due to the fact that the initial water content in the ed was not significant enough to cause considerale hydrodynamic changes. Among the many availale correlations, the following equations have een found to provide good matches with the experimental measurements: Bule velocity [1]: u U U.711gd mf 0 0.5 (17) This a modified form of Davidson and Harrison [13] equation, where 0. 75 and are correction factors suggested y Werther [17] and Hilligardt and Werther [18]. 1 3 3.2Dc Bule diameter [5]: 0.347D d 3 (18) where o D exp 0.3 z D exp 0. z D 0.652 o c 0. 4 no c D o 4D 2 c U U mf 0.4 (19) U mf used in Eqs 17 and 19 was given y: U mf d 2 p p g 150 g 3 mf 1 p mf (20) For the operating condition given in Ta. 2, Eq. 20 gives U mf =0.065 m/s, which closely matches the measured value of 0.062 m/s. Despite the fact that Eqs. 17-20 were all originally developed for dry ed operations; they seem to provide a reasonale match with the experimental measurements made here under wet ed condition (Figure 6). This is not surprising, since the water content in the ed was relatively low as discussed aove. The expanded ed height, used in the experimental estimation of the overall mass transfer coefficient (Eq. 4), is shown in Figure 7. A clear gradual (ut limited) increase in the ed expansion as the water is removed from the ed can e noticed. 10

4.2. Drying rate The drying rate curves for the three experiments conducted are shown in Figure 8. The curve fitting is used to otain the water content in the ed at various times. From Figure 8, it may e concluded that the time of drying is directly proportional to the initial water content, and inversely proportional to the drying air flow rate. For instance, at an air velocity of 0.47 m/s, this time was reduced y half when reducing the initial water content from C o, ed =10% to o ed =5%, while at the initial water content of from 0.47 m/s to 0.33 m/s. C o, ed =10, this time was ~35% longer when reducing the air velocity C, The water concentration in the ed as a function of the drying time is shown in Figure 9. This was otained from the integration of the drying curve function, F, and sutracting it from the initial water content, w o, such that: t w t w o t 0 F t dt (21) 4.3. Mass transfer coefficient 4.3.1 Experimental measurement The measured overall mass transfer coefficient, k d, as a function of the water concentration in the ed is shown in Figure 10. The values of k d falls within the range of 0.0145-0.021 m/s. It is interesting to note that this range is close to the value one can otain from the literature for the mass transfer coefficient from a free water surface to an adjacent slow moving amient air stream (~0.015 m/s) [19]. 4.3.2 Comparison with other experimental data Walker [2] and Sit and Grace [3] measured the mass transfer coefficient in a two-dimensional fluidized ed. The technique employed involves the injection of ozone-rich ule into an airsolid fluidized ed. Their experimental data for various particle sizes are compared with our measurement in Figure 11. Taking into consideration the differences in the experimental set-up and operating conditions, the agreement with our measurement appears satisfactory for the particle size considered in this study. 11

4.3.3 Comparison with theoretical prediction Figure 12 compares the measured mass transfer coefficient with the theoretical predictions otained from the formulations given in section 3.1. The oundaries for the overall mass transfer coefficient are given y: (i) a model accounting for cloud-ule contriution, dense-ule contriution as well as the ule through flow contriution, giving the lower limit (Eq. 14) and (ii) a model accounting for the cloud-ule contriution only, giving the upper limit (Eq. 10). The results also suggest that, within the operating conditions considered here, the drying may well e represented y a purely diffusional model, controlled y either the resistance residing at the dens-cloud interface, or the cloud-ule interface. Finally, Ta. 3 shows the numerical values of the various mass transfer contriutions otained from Eqs 8-16. It is shown that the estimated diffusional resistances, as well as the contriution from the ule throughflow, are all of the same order of magnitude. Previously, Geldart [20] argued that the ule throughflow is not important for small particles and may e neglected. According to our analysis, this may well e the case here. However, generalization of this conclusion should e treated with caution especially when dealing with larger particles. 5. Conclusion Mass transfer coefficient in a uling fluidized ed dryer has een experimentally determined. This work is the first to utilise an ECT system for this purpose. The ECT allowed for quantification of the ule diameter and velocity, as well as providing new insight into the ule-cloud-dense oundaries. The measured overall mass transfer coefficient was in the range of 0.045-0.021 m 2 /s. A simple hydrodynamic and mass transfer model, ased on the availale correlations was used to predict the mass transfer coefficient in a uling fluidized ed. Despite the complexity of the process, and the numer of assumption employed in this analysis, the model ased on pure diffusional mass transfer seems to provide satisfactory agreement with the experimental measurements. This work set the seen for future experimental investigations to otain a generalised correlation for the mass transfer coefficient in fluidized ed dryer, particularly that utilizes the ECT or other similar imaging techniques. Such a correlation is of vital importance for improved fluidized ed dryer design and operation in its widest application. A comprehensive experimental program, covering a wider range of operating conditions (particle size, gas velocity, water content, porous/non-porous particles) is recommended. 12

Nomenclature A [m 2 ] column/ed cross-sectional area A, B 1 1 [m1.5 s -1 ] parameters defined in Eqs. 15, 16 respectively C, [-] water concentration in the dense and ule phases respectively, kg/kg d C d, D [m] diameter D, D e [m 2 s -1 ] molecular and effective diffusivity respectively g [ms -2 ] gravity acceleration constant H [m] expanded ed height k [ms -1 ] overall mass transfer coefficient (etween dense and ule phases) d k [ms -1 ] mass transfer coefficient etween cloud and ule phases c k [ms -1 ] mass transfer coefficient etween dense and cloud phases dc m [kgs -1 ] mass flow rate P [-] relative solid fraction S [m 2 ] surface area U [ms -1 ] superficial gas velocity u [ms -1 ] velocity V [(m 3 ) volume w [g] ed water content z [m] axial coordinate Greek symols [-] ed voidage [-] ule fraction [kg.m -3 ] density [m] distance etween the centre of the two ECT sensors [-] particle sphericity p Suscripts c d mf ule cloud dense minimum fluidization 13

References [1] D. Kunii, and O. Levenspiel, Fluidization Engineering. Second edition, Butterworth- Heinemann, Boston 1991. [2] B. V. Walker, Transaction Institute of Chemical Engineers 1975, 53, 225. [3] S.P. Sit, and J. R. Grace, Chemical Engineering Science 1978, 33, 1115. [4] W. Ciesielczyk and J. Iwanowski, Powder Technology 2006, 24, 1153. [5] S. Mori and C. Y. Wen, AICHE J. 1975, 21, 109. [6] P. J. A. M. Kerkhof, Chemical Engineering and Processing 2000, 39, 69. [7] F. Scala, Chemical Engineering Science 2007, 62, 4159. [8] Y. Makkawi, P. Wright, R. Ocone, Powder Technology 2006, 163, 69. [9] Y. Makkawi and R. Ocone, Chemical Engineering Science 2007, 62, 4304. [10] G. Chaplin and T. Pugsley, Chemical Engineering Science 2005, 60, 7022. [11] G. Chaplin, T. pugsley, L. van der Lee, A. Kantzas, C. Winters, Measurement Science and Technology 2005, 16, 281. [12] Y. Makkawi and P. C. Wright, Powder Technology 2004, 148, 142. [13] J. Davidson and D. Harrison, Fluidized Particles, Camridge University Press, Camridge 1963. [14] J. D. Murray, Journal of Fluid Mechanics 1965, 22, 57. [15] R. Higie, in Fluidization Engineering, Second edition (Eds: D. Kunii and O. Levenspiel), Butterworth-Heinemann, Boston, 1991. 256. [16] T. Chia, and H. Koayashi, Chemical Engineering Science 1970, 25,1375. [17] J. Werther, in Fluidization Engineering, Second edition (Eds: D. Kunii and O. Levenspiel), Butterworth-Heinemann, Boston, 1991. 147 [18] K. Hilligardt and J. Werther, German Chemical Engineering 1986, 9, 215. [19] G. Saravacos and Z. Maroulis, Transport properties of foods, CRC Press, 2001. [20] D. Geldart, Gas Fluidization Technology, John Wiley and Sons, Chichester, UK, 1986. 14

Tale 1 Geldart Group A/B Particle size range (m) 50-180 Mean particle diameter (m) 125 Particle density (kg/m 3 ) 2500 Bulk density (kg/m 3 ) 1300 Hardness according to Mohs 6% Sphericity 80% Pores < 0.02 nm Material Pure soda lime glass allotini. Chemical composition SiO 2 =72%, Na 2 O=13%, CaO=9%, MgO=4%, Al 2 O 3 =1%, K 2 O & Fe 2 O 3 =1% Commercial name Glass eads type S, Art. 4500 Electric permittivity ~3.1 15

Tale 2 Experimental unit Column Distriutor Particles Operating conditions Diameter = 13.8 cm, height = 150 cm, material: cast acrylic Perforated PVC; 150 holes of 2 mm dia. d = 125 m, p = 2500 kg/m 3, Material: glass p Fluidization fluid Air at amient condition (20 C) Static ed height 20 cm Exp. 1 Exp. 2 Exp. 3 Fluidization velocity (m/s) 0.34 0.47 0.47 Initial water content (wt%) 1.0 1.0 0.5 16

Tale 3 Gas velocity Experimental ule characteristics Overall mass transfer coeff. Dense-cloud interchange (diffusion only) Theoretical Cloud-ule interchange (diffusion only) Bule throughflow U (m/s) d (m) u (m/s) k d (m/s) k dc (m/s)- Eq. 12 k c (m/s)- Eq. 8 k q (m/s)- Eq. 9 0.35 0.04 0.99 0.0145 0.0178 0.0194 0.015 17

List of Figs Figure 1. Experimental set up (a) Schematic of the fluidized ed () A photograph of the installation Figure 2. Schematic representation of the method used in experimental calculation of the overall mass transfer coefficient Figure 3. Estimation of ule velocity from ECT data Figure 4. Estimation of ule diameter from ECT measurement (a) ule diameter () ECT solid fraction () ECT slice images Figure 5. Mass transfer mechanism etween dense-cloud and cloud-ule phases demonstrated in a typical ECT image of an isolated rising clouded ule in a fluidized ed Figure 6. Variation of the ule velocity and ule diameter during the drying process Figure 7. Variation of expanded fluidized ed height during the drying process Figure 8. Drying rate curves for the three conducted experiments Figure 9. Variation of water content during drying Figure 10. Experimentally measured overall mass transfer coefficient Figure 11. Experimental overall mass transfer coefficient in comparison with other previously reported results Figure 12. Comparison etween experimental measurement and various theoretical models for mass transfer coefficients List of Tales Tale 1. Physical and chemical properties of the dry particles Tale 2. Summary of experimental operating conditions Tale 3. The measured overall mass transfer coefficient for one selected operating condition in comparison to the theoretical predictions of various contriutions. 18

Summary The drying of moist particles in a uling fluidized ed mainly takes place as a result of mass transfer etween the ules and the surrounding dense solid phase. In this study, the mass transfer coefficient was estimated from experimental measurement of the ule characteristics and drying air condition. The results were used for the validation of empirical correlations and for assessing the various assumptions used in developing these correlations. 19

Air to atmosphere Cast acrylic column Temp/humidity proe 13.8 cm Wet ed material ECT sensor Upper sensor Lower sensor 1500 cm 7.6 cm 3.8 cm 3.8 cm 7.6 cm Gas distriutor 25 cm (a) Fluidized ed ECT sensors () Data acquisition system Drying air Figure 1 20

Wet dense phase u C d H dz Bule C u<<u air Figure 2 21

0.95 upper sensor low er sensor relative solid fraction, P (-) 0.85 0.75 0.65 0.55 ules peak mean 0.45 t 1 t 2 3.05 3.25 3.45 time (second) t Figure 3 22

time (second) 0.85 1.15 1.45 1.75 2.05 2.35 0.08 0.07 ule diameter (m) 0.06 0.05 0.04 0.03 0.02 0.01 0 0 1 2 3 4 5 ule numer (-) relative solid fraction, P (-) 1 0.95 0.9 0.85 0.8 0.75 0.7 0.65 0.6 0.55 0.5 mean 0.85 1.15 1.45 1.75 2.05 2.35 time (seconds) 4 5 3 un-accounted ule un-accounted ule 2 1 Figure 4 23

0.06 m ule dense phase q cloud d k c k dc Figure 5 24

ule velocity, U (m/s) 1.5 1.3 1.1 0.9 0.7 0.5 0.07 (a) U=0.47 m/s U=0.34 m/s U = 0.47 m/s U = 0.35 m/s Eq. 17 0 2 4 6 8 10 water content, C (g/kg solid) ule diameter, d (m) 0.06 0.05 0.04 0.03 0.02 0.01 U = 0.47 m/s U = 0.35 m/s Eq. 18 0 2 4 6 8 10 water content, C (g/kg solid) Figure 6 25

ed expansion, H (m) 0.44 0.40 0.36 0.32 0.28 0.24 U = 0.47 m/s U = 0.35 m/s 0.20 0.0 2.0 4.0 6.0 8.0 10.0 water content, C (g/kg solid) Figure 7 26

drying rate (g/min) 7 6 5 4 3 2 F(t) =-5.1E-11t 6 +7.0E-9t 5-3.6E7t 4 +8.2E-6t 3-6.8E-5t 2-3.2E-4t +0.0059 C o,ed C o,ed C o,ed = 5, U = 0.47 m/s = 10, U = 0.47 m/s = 10, U = 0.34 m/s polynomial fit 1 0 0 5 10 15 20 25 30 35 40 45 50 time (minute) Figure 8 27

water in ed, C ed (g/kg dry solid) 12 10 8 6 4 2 0 C o,ed C o,ed C o,ed = 5, U = 0.47 m/s = 10, U = 0.47 m/s = 10, U = 0.34 m/s polynomial fit 0 5 10 15 20 25 30 35 40 45 50 55 60 time (minute) Figure 9 28

overall mass transfer coefficient, k d (m/s) 0.042 0.036 0.030 0.024 0.018 0.012 0.006 0.000 U=0.35 m/s U=0.47 m/s 0 2 4 6 8 water content (g/kg solid) Figure 10 29

overall mass transfer coefficient, kd (m/s) 0.03 0.025 0.02 0.015 0.01 0.005 0 0 100 200 300 400 particle diameter, d p (m) Sit and Grace (1977) Walker (1975) This study Figure 11 30

mass transfer coefficient (m/s) 0.05 0.04 0.03 0.02 0.01 (a) () (c) (d) () (a) (c) (d) measured cloud-ule transfer coeff (diffusion only, Eq 8) cloud-ule transfer coeff (diffusion-convection, Eq 10) dense-cloud transfer coeff (diffusion only, Eq 12) overall transfer coeff (diffusion-convection, Eq 14) 0.00 0.032 0.036 0.040 0.044 0.048 0.052 ule diameter (m) Figure 12 31