Flow behavior of wakes in a three-phase slurry bubble column with viscous liquid medium

Similar documents
Multiple effects of operating variables on the bubble properties in three-phase slurry bubble columns

World Academy of Science, Engineering and Technology International Journal of Chemical and Molecular Engineering Vol:7, No:12, 2013

Hydrodynamic and Mass Transfer Characteristics of External-Loop Airlift Reactors without an Extension Tube above the Downcomer

CHAPTER 2 EXPERIMENTAL SETUP AND PROCEDURE

Numerical Simulations of a Train of Air Bubbles Rising Through Stagnant Water

Simulation of Gas Holdup in a Bubble Column with a Draft Tube for Gas Dispersion into an Annulus

Experimental Studies on the Instabilities of Viscous Fingering in a Hele-Shaw Cell

EXPERIMENTAL STUDY ON THE HYDRODYNAMIC BEHAVIORS OF TWO CONCENTRIC CYLINDERS

EFFECTS OF CHEMICAL ADDITIVES ON THE PRESSURE DROP IN THE PIPES

Development of High-speed Gas Dissolution Device

Air entrainment in Dip coating under vacuum

Determination of Bubble Size Distribution in Pressurized 3D Bubble Columns Using Borescope-Based Imaging Technique

Instruction Manual. Pipe Friction Training Panel

Journal of Chemical Engineering of Japan, Advance Publication doi: /jcej.12we164; published online on July 25, 2012

Visual Observation of Nucleate Boiling and Sliding Phenomena of Boiling Bubbles on a Horizontal Tube Heater

Measurement of both gas and liquid velocity profiles for bubble-induced turbulent flow

EFFECT OF CORNER CUTOFFS ON FLOW CHARACTERISTICS AROUND A SQUARE CYLINDER

Irrigation &Hydraulics Department lb / ft to kg/lit.

MODELING AND SIMULATION OF VALVE COEFFICIENTS AND CAVITATION CHARACTERISTICS IN A BALL VALVE

LOW PRESSURE EFFUSION OF GASES adapted by Luke Hanley and Mike Trenary

The water supply for a hydroelectric plant is a reservoir with a large surface area. An outlet pipe takes the water to a turbine.

EWGAE 2010 Vienna, 8th to 10th September

Applied Fluid Mechanics

The Effect of Von Karman Vortex Street on Building Ventilation

Experimental Investigation of the Rise Behavior of Live-Oil Droplets during Deep-Sea Oil Spills

Effect of Nozzle Twisted Lance on Jet Behavior and Spitting Rate in Top Blown Process

Analysis of Pressure Rise During Internal Arc Faults in Switchgear

Effect of the cross sectional shape of the recirculation channel on expulsion of air bubbles from FDBs used in HDD spindle motors

Bioreactor System ERT 314. Sidang /2011

DIAGNOSTICS OF IMPULSE LINE BLOCKAGE WITH A MULTI-SENSING DIFFERENTIAL PRESSURE TRANSMITTER AT THE AIR LINE

ANNUAL REPORT UIUC, August 16, Bubble Formation, Breakup and Coalescence in Stopper-rod Nozzle Flow and Effect on Multiphase Mold Flow

LOW PRESSURE EFFUSION OF GASES revised by Igor Bolotin 03/05/12

WATER HYDRAULIC HIGH SPEED SOLENOID VALVE AND ITS APPLICATION

Numerical Investigation of Air Bubbles Evolution and Coalescence from Submerged Orifices Based on OpenFOAM

Questions. theonlinephysicstutor.com. facebook.com/theonlinephysicstutor. Name: Edexcel Drag Viscosity. Questions. Date: Time: Total marks available:

Preliminary Assessments of Combined Effects of Surface Tension and Viscosity on Bubble Column Hydrodynamics

For Multi-Parameter Meters see mvx

Gas-liquid two-phase flow in a downward facing open channel

ON THE EFFECT OF LIFT FORCES IN BUBBLE PLUMES

An underwater explosion is an explosion where the point of detonation is below the surface of the water.

LUGB Vortex Flowmeter

Gas Vapor Injection on Refrigerant Cycle Using Piston Technology

Effect of Coiled Capillary Tube Pitch on Vapour Compression Refrigeration System Performance

Integral type Differential pressure flowmeter VNT Series

The Study of Bubbly Gas-Water Flow Using

ZIN Technologies PHi Engineering Support. PHi-RPT CFD Analysis of Large Bubble Mixing. June 26, 2006

Experimental Study of an Air Lift Pump

EXPERIMENTAL STUDY ON SNOW BEHAVIOR AROUND FENCES INSTALLED ALONG ELEVATED HIGHWAY

Measurement of Bubble Velocity using Spatial Filter Velocimetry

Wind Tunnel Study on the Structural Stability of a Container Crane According to the Boom Shape

NEW PROGRAM! OIL AND GAS TECHNOLOGY PROGRAM

CFD SIMULATIONS IN AN INTERNAL CIRCULATION AIRLIFT OPERATING UNDER HOMOGENEOUS REGIME

Using sea bed roughness as a wave energy dissipater

Vortex flowmeters. Product family introduction Principle of operation Product review Applications Key product features

For Multi-Parameter Meters see mvx

Effect of Perforated Plate Open Area on Gas Holdup in Rayon Fiber Suspensions

Experiment 8: Minor Losses

Workshop 1: Bubbly Flow in a Rectangular Bubble Column. Multiphase Flow Modeling In ANSYS CFX Release ANSYS, Inc. WS1-1 Release 14.

Vortex Flow Meter Wafer or Flange Connection. - Steam - Liquid - Gas

VORTEX SHEDDING AND VORTEX FORMATION FROM A PAIR OF IN-LINE FORCED OSCILLATING TANDEM ARRANGED CIRCULAR CYINDERS IN A UNIFORM FLOW

Dynamic Surface Tensiometer Bubble Pressure Method BP-D5, BP-D5L. Kyowa Interface Science Co., Ltd.

Fluid Flow. Link. Flow» P 1 P 2 Figure 1. Flow Model

QUARTERLY REPORT NO. 3 DE-FG22-94PC /01/95

Numerical simulation of an intermediate sized bubble rising in a vertical pipe

CFD Analysis of the Anti-Surge Effects by Water Hammering

of air-lift reactors based on bubble dynamics

A Study on the Human Impulse Characteristics Of the Standing Shooting Posture

New Approaches for Theoretical Estimation of Mass Transfer Parameters in Both Gas-Liquid and Slurry Bubble Columns

Dynamics of bubble rising at small Reynolds numbers

The effect of two inclined circular plunging jets on air entrainment in an aeration tank

Variation in Pressure in Liquid-Filled Plastic Film Bags Subjected to Drop Impact

Coriolis Mass Flow Meter

CRYSTALLIZATION FOULING IN PACKED COLUMNS

Aerodynamic Measures for the Vortex-induced Vibration of π-shape Composite Girder in Cable-stayed Bridge

Level MEASUREMENT 1/2016

Two-stage pulse tube refrigerator with step displacer as phase shifter

Use a Controlled Vibration to Mixing and Separation of a Gas Bubbles and a Liquid Under Reduced and Microgravity Conditions

UNIT 15 WATER HAMMER AND SURGE TANKS

International Journal of Advanced Engineering and Management Research Vol. 1 Issue Pressure analysis of a Multi-tube Cyclone Separator

EXPERIMENTAL EVALUATIONS OF LOCAL BUBBLE PARAMETERS OF SUBCOOLED BOILING FLOW IN A PRESSURIZED VERTICAL ANNULUS CHANNEL

Akasison Flow phenomena of a siphonic roof outlet

AN ANALYSIS ON HIGH PRESSURE DYNAMIC CALIBRATORS USED IN THE DEFENSE AREAS

The Discussion of this exercise covers the following points:

Section 2 Multiphase Flow, Flowing Well Performance

Characteristics of Decompression Tank Internally Pressurized With Water Using OpenFOAM Syamsuri 1, a

Study on Fire Plume in Large Spaces Using Ground Heating

Air Bubble Defects in Dispensing Nanoimprint Lithography

To plot the following performance characteristics; A pump is a device, which lifts water from a lower level to a higher

Vortex Meters for Liquids, Gas, and Steam

DIGITAL IMAGE ANALYIS OF BUBBLE FLOW DISTRIBUTION Influence of operational parameters

Correction of Pressure Drop in Steam and Water System in Performance Test of Boiler

Investigation of Suction Process of Scroll Compressors

SPH 4C Unit 4 Hydraulics and Pneumatic Systems

The Estimation Of Compressor Performance Using A Theoretical Analysis Of The Gas Flow Through the Muffler Combined With Valve Motion

IBERIAN SPH 2015 SPH MODELLING FOR AIR ENTRAINMENT IN WAVE BREAKING

An innovative technology for Coriolis metering under entrained gas conditions

Old-Exam.Questions-Ch-14 T072 T071

CFD Analysis ofwind Turbine Airfoil at Various Angles of Attack

Experimental investigation of the check valve behaviour when the flow is reversing

Determination of the wind pressure distribution on the facade of the triangularly shaped high-rise building structure

Transcription:

Korean J. Chem. Eng., 28(3), 974-979 (2011) DOI: 10.1007/s11814-010-0403-4 INVITED REVIEW PAPER Flow behavior of wakes in a three-phase slurry bubble column with viscous liquid medium Dae Ho Lim*, Ji Hwa Jang*, Yong Kang*,, and Ki Won Jun** *School of Bio and Applied Chemical Engineering, Chungnam National University, Daejeon 305-764, Korea **Green Chemical Technology Division, Korea Research Institute of Chemical Technology, Daejeon 305-600, Korea (Received 22 June 2010 accepted 16 August 2010) Abstract Flow behavior of wakes has been investigated in a three-phase slurry bubble column of 0.102 m ID and 1.5 m in height. The dependence of wake characteristics such as rising velocity, frequency, holdup and equivalent size on the operating variables was examined by employing an electric resistivity probe method. The gas velocity, liquid viscosity and solid content in the slurry phase were chosen as independent parameters. The rising velocity of wake region increased with an increase in the gas velocity (4.0-12.0 cm/s), liquid viscosity (1.0-50.0 mpa s) or solid content (0-25 wt%) in the slurry phase. The frequency and holdup of wake phase increased with increasing gas velocity, but decreased with increasing liquid viscosity or solid content in the slurry phase. The equivalent size of wake phase increased with increasing gas velocity, liquid viscosity or solid content in the slurry phase. The wake properties and holdup were well correlated with operation variables within these experimental conditions. Key words: Wake Behavior, Three-phase, Slurry, Bubble Column, Viscous Medium INTRODUCTION To whom correspondence should be addressed. E-mail: kangyong@cnu.ac.kr Taking advantage of several unique features such as simplicity, low operating cost, suitability for continuous operation and effective contacting among multi-phases, the bubble columns have been utilized as reactors and contactors for continuous multi-phase processes [1,2]. One of them is the CTL (coal-to-liquid) or GTL (gasto-liquid) process for the production of synthetic liquid fuels [3-5]. Investigations have been conducted on the bubble columns and slurry bubble columns for industrial applications and for the new emerging fields [3-7]. A bubble column is composed of rising bubbles in a continuous liquid medium. Although it looks simple, the flow behavior of rising bubbles in the column is highly irregular, since the bubbles rise randomly as a dispersed phase, propagating force and energy into the surroundings by means of coalescence and splitting. In addition, the rising bubbles are tailed by the wakes behind their rising path, due to the differences of pressures and thus streamlines between the front and the back sides of the rising bubbles. Owing to the vertical motions in the wake, the contacting intensity and efficiency among multi-phases have been extremely high in the wakes. Therefore, the energy dissipation rate, which is closely related to the transport phenomena such as momentum heat and mass transfer, in the wake can be very high [7-9]. For the analysis of the wake phenomenon which is induced by rising bubbles, the types and formation - shedding mechanism of wake phase have been studied [10,11]. Although extensive investigations have been conducted on the flow behavior of dispersed phase such as bubbles and wakes, there still remains uncertainty especially on the wake characteristics in three-phase slurry bubble columns. In addition, most of the previous investigations have focused on the wake phenomenon or phase which is tailing at the rear of a single bubble. The wake flow behavior induced by a single bubble seems to be quite important for the analysis of the mechanism; however, it has been highly sophisticated to be utilized for the actual slurry bubble column reactors. Fig. 1. Schematic diagram of experimental apparatus. 1. Main column 09. Needle valve 2. Gas distributor 10. Valve 3. Pressure probe 11. Liquid flowmeter 4. Resistivity probe 12. Gas flowmeter 5. Pressure sensor 13. Liquid pump 6. A/D converter 14. Liquid reservoir 7. Computer 15. Resistivity circuit 8. Compressor 16. Date acquisition system 974

Table 1. Physical and rheological properties of liquid phase Dynamic viscosity (mpa s) Flow behavior of wakes in three-phase slurry bubble column 975 Surface tension (mn/m) Density (kg/m 3 ) K (Pa s n ) n Diffusivity (cm 2 /s) Kinematic viscosity (m 2 /s) Water 1 72.9 1000 0.001 1 2.22 10 5 9.61 10 7 CMC 0.1 wt% 11 73.2 1001 21.69 10 3 0.882 0.48 10 5 1.10 10 5 CMC 0.2 wt% 24 73.3 1002 43.82 10 3 0.847 0.26 10 5 2.40 10 5 CMC 0.3 wt% 38 73.6 1003 71.69 10 3 0.825 0.19 10 5 3.79 10 5 CMC 0.4 wt% 50 73.9 1004 0.102 10 3 0.802 0.15 10 5 4.98 10 5 Therefore, in the present study, characteristics of wake phase were examined in a three-phase slurry bubble column, where various bubbles with different sized are generated and their coalescence and breakage could occur simultaneously and continuously. Because, the resultant flow behavior of wake phase could be utilized for the practical applications in the various fields. EXPERIMENTS Fig. 1 shows the schematic diagram of experimental apparatus which was made of stainless-steel column (0.102 m 1.5 m high). A perforated plate, which was served as a gas distributor, was installed between the main column section and a 0.2 m high stainless-steel distributor box. The slurry phase was composed of aqueous solutions of carboxymethyl cellulose (CMC) and glass bead. The viscosity of liquid phase was in the range of 1.0-50 mpa s, with a variation of CMC amount added in the tap water. The physical properties of liquid phase are listed in Table 1. The size of glass bead was in the range of 40-70 µm and solid content in the slurry phase was 0-20 wt%. Filtered, oil-free compressed air, which was used as a gas phase, was admitted at the bottom of the column through the gas distributor [12-14]. After steady state was attained in an adjusted experimental condition, the bubbles as well as wakes were measured. The wake characteristics such as rising velocity, frequency, holdup and size were measured by means of a dual electrical resistivity probe method. The probe, which was applied by 1.75 V DC, was installed at 0.4 m from the gas distributor. The probe consisted of two 7 mm diameter stainless-steel pipes coated with epoxy resin. The vertical and horizontal distances between the tips of the two probes were 2 mm and 3 mm, respectively. The probes were located at the center between the wall and the center of the column. The tips of the probes were made of 0.2 mm platinum wire to detect the difference in conductivity of gas and liquid [5,7]. The detected analogy signals, which were obtained at the preselected sampling rate, were stored on a personal computer with the DT2805 Lab card. The sampling rate and time were 1,000 HZ and 30 sec, respectively, comprising a total sample length of 30,000, which was enough to detect the fuel spectrum of the signals [4,5]. The wake characteristics were determined from the relationship be- Fig. 2. Typical examples of bubble and wake in a viscous bubble column. U G 10 2 [m/s] µ L 10 3 [Pa s] S C [wt%] (a) : 10 11 5 (b) : 12 11 5 Fig. 3. Effects of gas velocity on the rising velocity and frequency of wake in a three-phase slurry bubble column. µ L 10 3 [Pa s] : 1 05 11 38 S C [wt%] : 5 25 25 25 Korean J. Chem. Eng.(Vol. 28, No. 3)

976 D. H. Lim et al. tween the reformed and digitized probe signals and the wake dwell and lag time. RESULTS AND DISCUSSION Typical examples of bubble and wake phase detected by the probes in the three-phase slurry bubble columns can be seen in Fig. 2. The wake region can be distinguished from the bubble phase at the back of the bubble. Effects of gas velocity on the frequency and rising velocity of wake phase, which were obtained from the relationship between the reformed and digitized probe signals, can be seen in Fig. 3. The frequency and rising velocity increase with increasing gas velocity. The increase of gas velocity can lead to the increase of gas input into the slurry bubble column, which results in the increase of the number of bubbles and thus of bubble coalescence in the column. Since the bubbles are tailed by wakes as a dispersed phase, the number of wakes could also increase in line with that of bubbles, with an increase in the gas velocity. In addition, the wake regions could also coalesce with each other whenever bubble coalescence could occur. It has been reported that the bubble rising velocity increases with increasing gas velocity [14-16]; thus, the rising velocity of the wake region can also increase with increasing gas velocity, since the rising bubbles are tailed by wakes. Effects of liquid viscosity on the frequency and rising velocity of wake phase can be seen in Fig. 4. The frequency decreases but the rising velocity increases, with an increase in the liquid viscosity. It is noted that the bubble frequency decreases but the rising velocity of bubbles increases, with increasing liquid viscosity. These can be due to the fact that the bubble coalescence could increase with increasing liquid viscosity [14-16]. Since the buoyancy force acting on the bubbles, which is due to the difference of density between the bubble and slurry phase, increases with increasing bubble size, the rising velocity of wakes as well as bubbles could increase with increasing liquid viscosity. Because, the bubble coalescence increases and thus the frequency of bubbles decreases, with an increase in the liquid viscosity in three-phase slurry bubble columns [6-9]. Effects of solid content on the frequency and rising velocity of wake region can be seen in Fig. 5. The frequency of wake phase decreases gradually but the rising velocity increases, with increasing solid content in the slurry phase. The amount of solid content in the slurry phase could act on the bubble coalescence in that the solid particle could inhibit bubbles from contacting with slurry phase, which results in the increase of probability of bubbles contacting each other, that is, the increase of bubble coalescence. Since the coalescence and breaking or disintegration of wakes are quite similar to those of bubbles in the three-phase slurry bubble columns, the effects of solid content on the wake motion and flow behavior could be also explained as motioned above. Effects of gas velocity on the wake holdup can be seen in Fig. 6. The holdup of wake was determined from the relationship between the reformed and digitized probe signals of wake dwell and lag time [7,11]. As expected, the wake holdup increases with increasing gas velocity in three-phase slurry columns. As indicated, the wake phase Fig. 4. Effects of liquid viscosity on the rising velocity and frequency of wake in a three-phase slurry bubble column. U G 10 2 [m/s] : 4 04 10 12 S C [wt%] : 5 15 25 25 Fig. 5. Effects of solid content in the slurry phase on the rising velocity and frequency of wake in a three-phase slurry bubble column. U G 10 2 [m/s] : 4 08 10 12 µ L 10 3 [Pa s] : 1 11 11 01 March, 2011

Flow behavior of wakes in three-phase slurry bubble column 977 Fig. 6. Effects of gas velocity on the wake holdup in a three-phase slurry bubble column. µ L 10 3 [Pa s] : 11 01 11 38 S C [wt%] : 05 25 25 25 is a tailing at the rear of the rising bubbles, and thus, the amount of tailing could increase with increasing bubble holdup by increasing gas input into the slurry column. Therefore, we can state that the gas-liquid contact could be more effective with increasing gas velocity, since the wake holdup increases with increasing gas velocity. It has been understood that the contacting flow of fluid element in the wake phase is quite vigorous owing to the vortex motion at the Fig. 7. Effects of liquid viscosity on the wake holdup in a threephase slurry bubble column. U G 10 2 [m/s] : 04 08 10 10 S C [wt%] : 15 15 15 05 Fig. 8. Effects of solid content in the slurry phase on the wake holdup in a three-phase slurry bubble column. U G 10 2 [m/s] : 4 08 10 12 µ L 10 3 [Pa s] : 1 11 11 38 rear of the rising bubbles. Since the wake phase has been assumed to be composed mainly of liquid or slurry phase, the contact between the gas bubbles and liquid phase in the wake region could be most effective in the column. This can be the reason why the energy dissipation rate can be high in the wake region [1,2,11-13]. Effects of liquid viscosity and solid content in the slurry phase on the wake holdup in the column can be seen in Figs. 7 and 8. Here the wake holdups decrease gradually with increasing liquid viscosity or solid content in the slurry phase. Since the bubble coalescence and thus wake coalescence could increase with increasing liquid viscosity or solid content in the slurry phase, the holdups of wake region could decrease with increasing liquid viscosity or solid content in the slurry phase, as in the cases of bubble holdups [1,2]. The equivalent size of wakes was determined from the determined values of frequency and holdup, and this value was checked based on the concept that the shape of wake would be ellipsoidal. And then, the equivalent diameter of it was determined as a spherical shape which has the same volume as that of an ellipsoid. Thus obtained equivalent size of wake region is shown in Figs. 9-11, as a function of gas velocity liquid viscosity and solid content in the slurry, respectively. In these figures, the wake size increases with increasing gas velocity, liquid viscosity or solid content in the slurry phase. As indicated, the size of the wake region can be increased in the condition that the coalescence between wake regions could be more frequent. Since the contact and coalescence of the wake region are accompanied by those bubbles, the wake size increases in line with the increase of bubble size. In other words, the coalescence and contact of the wake region can increase with increasing gas velocity, liquid viscosity or solid content in the slurry phase. The frequency, rising velocity, holdup and equivalent size of wake region have been well correlated in terms of operating variables as Korean J. Chem. Eng.(Vol. 28, No. 3)

978 D. H. Lim et al. Fig. 9. Effects of gas velocity on the wake equivalent diameter in a three-phase slurry bubble column. µ L 10 3 [Pa s] : 1 01 11 38 S C [wt%] : 5 15 15 15 Fig. 11. Effects of solid content in the slurry phase on the wake equivalent diameter in a three-phase slurry bubble column. U G 10 2 [m/s] : 08 10 10 12 µ L 10 3 [Pa s] : 11 11 38 38 0.51 0.29 0.25 d W =0.056U G S C µ L The correlation coefficients of Eqs. (1)-(4) are 0.90, 0.92, 0.91 and 0.92, respectively. CONCLUSION (4) Fig. 10. Effects of liquid viscosity on the wake equivalent diameter in a three-phase slurry bubble column. U G 10 2 [m/s] : 4 04 08 10 S C [wt%] : 5 15 15 05 Eqs. (1)-(4): The wake properties such as frequency, rising velocity, and equivalent size and its holdup were successfully examined in a three-phase slurry bubble column by employing the resistivity probe method. The frequency of wake region increased with increasing gas velocity but decreased with increasing liquid viscosity or solid content in the slurry phase. The rising velocity and equivalent size of wake region increased with increasing gas velocity, liquid viscosity or solid content in the slurry phase. The holdup of wake region increased with increasing gas velocity, but it decreased with increasing liquid viscosity or solid content in the slurry phase. The wake properties and its holdup were well correlated in terms of operating variables as: 0.62 0.24 0.21 F W =3.150U G S C µ L 0.78 0.37 0.15 U W =3.174U G S C µ L 0.29 0.36 0.11 ε W =0.070U G S C µ L 0.51 0.29 0.25 d W =0.056U G S C µ L ACKNOWLEDGEMENT 0.62 0.24 0.21 F W =3.150U G S C µ L 0.78 0.37 0.15 U W =3.174U G S C µ L 0.29 0.36 0.11 ε W =0.070U G S C µ L March, 2011 (1) (2) (3) This work was supported by the New & Renewable energy of the Korea Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korea Government Ministry of Knowledge Economy (3A-3-1, 2009).

Flow behavior of wakes in three-phase slurry bubble column 979 d P d W NOMENCLATURE : particle diameter [µm] : equivalent diameter [cm] D : column diameter [m] F w : wake frequency [Hz] K : consistency index in power-law model [PaS n ] L : column height [m] S C : solid content [wt%] t : time [sec] U G : gas velocity [m/s] : wake rising velocity [cm/s] U W Greek Letters a : diffusivity [cm/s] ε w : wake phase holdup [-] µ L : liquid viscosity [mpa s] n : kinematic viscosity [m/s] ρ G : gas density [kg/m 3 ] ρ L : liquid density [kg/m 3 ] ρ S : solid content density [kg/m 3 ] s : surface tension [mn/m] REFERENCES 1. L. S. Fan and K. Tsuchiya, Bubble dynamics in liquids and liquidsolid suspension, Stoneham, MA, Butterworth-Heinemann (1990). 2. W. D. Deckwer, Bubble Column Reactors, John Wiley and Sons, NY (1992). 3. R. Krishna and S. T. Sie, Fuel Processing Technol., 64, 73 (2000). 4. J. H. Jang, M. J. Seo, D. H. Lim, Y. Kang, H. Jung and H. T. Lee, Korean Chem. Eng. Res., 47, 587 (2009). 5. I. S. Shin, S. M. Son, U. Y. Kim, Y. Kang, S. D. Kim and H. Jung, Korean J. Chem. Eng., 26, 587 (2009). 6. H. Li and A. Prakash, Powder Technol., 113, 158 (2000). 7. Y. Kang, Y. J. Cho, K. J. Woo, K. I. Kim and S. D. Kim, Chem. Eng. Sci., 55, 411 (2000). 8. R. C. Chen, F. M. Wang and T. J. Lim, Chem. Eng. Sci., 54, 4831 (1999). 9. I. Komasawa and T. Otake, J. Chem. Japan, 13, 103 (1980). 10. J. Chen, P. Gupta, S. Degaleesan, M. H. Al-Dajjam, M. P. Dudukovic and B. A. Toseland, Flow Measurement and Instrumentation, 9, 91 (1998). 11. K. Tsuchiya and L. S. Fan, Chem. Eng. Sci., 43, 1167 (1988). 12. S. M. Son, K. I. Lee, S. H. Kang and S. D. Kim, AIChE J., 53, 3011 (2007). 13. Y. Kang and S. D. Kim, I&EC Process Des. Dev., 25, 717 (1986). 14. S. H. Kang, S. M. Son, Y. Kang, J. W. Bae and K. W. Jun, Korean J. Chem. Eng., 25, 897 (2008). 15. S. M. Son, I. S. Shin, S. H. Kang, Y. Kang and S. D. Kim, Korean J. Chem. Eng., 24, 866 (2007). 16. U. Y. Kim, S. M. Son, S. H. Kang, Y. Kang and S. D. Kim, Korean J. Chem. Eng., 24, 892 (2007). Korean J. Chem. Eng.(Vol. 28, No. 3)