Acessibilidade / Reportar erro

Oxidation of zinc sulphide concentrate in a fluidised bed reactor - Part 1: Characterization of the fluid dynamics of the particle bed

Abstract

The oxidation kinetics of a blend of zinc sulphide concentrates in a bench-scale fluidised bed reactor operated in batch regime was studied. This work involves the fluid dynamics characterisation of the bed during the reaction. Its purpose is to determine bed behaviour and to assess the experimental conditions used in the roasting experiments. The characterisation is derived from both fluid dynamics maps and an analysis based on the dynamics of gas bubble growth in the bed. Bed behaviour under any experimental condition is evaluated with the fluid dynamics maps. The transitions fixed/fluidised bed, fluidised bed/pneumatic transport and between several fluidisation regimes are obtained from curves of dimensionless gas velocity versus dimensionless particle diameters. The possibility of fluidisation in the slugging regime is verified by the dynamics of gas bubble growth in the bed.

Zinc sulphide oxidation; Fluidised bed; Fluid dynamics


TRANSPORT PHENOMENA AND UNIT OPERATION

Oxidation of zinc sulphide concentrate in a fluidised bed reactor – Part 1: Characterization of the fluid dynamics of the particle bed

C. A. R. Queiroz* * To whom correspondence should be addressed ; R. J. Carvalho; F. J. Moura

Departamento de Ciência dos Materiais e Metalurgia, Pontifícia Universidade Católica do Rio de Janeiro, Phone: +(55) (021) 3114-1563, Fax: +(55) (021) 3114-1236, Rua Marquês de São Vicente 225, CEP 22453-900, Gávea, Rio de Janeiro - RJ, Brazil. E-mail: rjcar@dcmm.puc-rio.br

ABSTRACT

The oxidation kinetics of a blend of zinc sulphide concentrates in a bench-scale fluidised bed reactor operated in batch regime was studied. This work involves the fluid dynamics characterisation of the bed during the reaction. Its purpose is to determine bed behaviour and to assess the experimental conditions used in the roasting experiments. The characterisation is derived from both fluid dynamics maps and an analysis based on the dynamics of gas bubble growth in the bed. Bed behaviour under any experimental condition is evaluated with the fluid dynamics maps. The transitions fixed/fluidised bed, fluidised bed/pneumatic transport and between several fluidisation regimes are obtained from curves of dimensionless gas velocity versus dimensionless particle diameters. The possibility of fluidisation in the slugging regime is verified by the dynamics of gas bubble growth in the bed.

Keywords: Zinc sulphide oxidation; Fluidised bed; Fluid dynamics.

INTRODUCTION

Considering the growing competitiveness of the world market and the ever more restrictive and specific environmental regulations, the importance of optimising industrial processes and minimising harmful emissions in the environment is paramount.

With the objective of evaluating the oxidising roasting process, unit operation of zinc production, a study of the oxidation kinetics of a blend of zinc sulphide concentrates in a bench-scale fluidised bed reactor operated in batch regime was carried out at DCMM-PUC-Rio.

It is well known that the fluidised bed reactors used in pyrometallurgical processes involving gas-solid reactions are strongly affected by the fluid dynamics of the system. Representative experimental results and mathematical modelling are necessary for scaling up the reactor. Therefore, a detailed characterisation of the fluid dynamics of the particle bed is extremely important.

In this context, the fluid dynamics of the particle bed can be characterised by both the fluid-solid and Geldart diagrams and fluid dynamics maps and an analysis based on the dynamics of gas bubble growth in the bed.

The fluid-solid diagrams show regions where the bed remains fixed or fluidised or when pneumatic transport occurs. These diagrams are obtained by setting the boundaries between these states from relationships involving Reynolds and Froude dimensionless numbers. Within these boundaries, the operational conditions, particle diameter and superficial gas velocity, are plotted, which allows establishment of the fluidisation regime.

The Geldart diagram shows the relationships of the difference between the solid and gas densities and the particle diameter. Bed behaviour is obtained according to a practical classification, suggested by Geldart (Abrahamsen and Geldart, 1980). This diagram gives qualitative information about bed behaviour during incipient fluidisation. More details on the theoretical construction of this diagram can be found in Queiroz (1997).

The fluid-solid and Geldart diagrams are practical but insufficient for determining bed behaviour when the superficial gas velocity used is higher than that necessary to obtain the incipient fluidisation of the particles. In this case, several fluidisation regimes occur. These regimes are particulate, bubbling, slugging, turbulent and fast fluidisation. The fluid dynamics maps are useful for determining the fluidisation behaviour under any operational condition. In this diagram, the boundaries between the several fluidisation regimes are established from curves of dimensionless gas velocity as a function of the dimensionless particle diameter. Finally, if the slugging regime is possible, the dynamics of gas bubble growth in the bed should be analysed.

The main objective of this work is to determine the fluid dynamics behaviour of particle beds of a blend of zinc sulphide concentrates using air as reactant gas. The concepts used to plot the fluid dynamics maps are presented and discussed. The dynamics of gas bubble growth in the bed are analysed when there is the possibility of slugging of the bed.

Fluid Dynamics Maps

The fluid dynamics maps consist of diagrams that provide the bed behaviour under any operating condition (Grace, 1990). In these diagrams, the boundaries between the fluid dynamics regimes are obtained from curves of dimensionless gas velocity (U*) versus dimensionless particle diameter (dp*), defined as

where dp is the mean particle diameter (m), rs is the particle density (kg/m3), rg is the gas density (kg/m3), m is the gas viscosity (Pa.s), g is the acceleration due to gravity (m/s2) and U is the superficial gas velocity (m/s).

The relationships between U* and dp* given below provide the boundaries between the fluid dynamics regimes of the beds.

Fixed/Fluidised Bed Boundary

The expression for the dimensionless minimum fluidisation velocity (U*mf) in terms of the dimensionless particle diameter gives the fixed/fluidised bed boundary:

where C1 and C2 are the dimensionless empirical constants for particle shape factor and bed porosity, C1 = 25.7 and C2 = 0.0365, as suggested by Richardson and Grace (Grace, 1982).

Fluidised Bed/Pneumatic Transport Boundary

The dimensionless terminal velocities () that define the fluidised bed/pneumatic transport boundary are given by

where Re is the Reynolds number (-).

To draw the lines that define the transition between the particulate, bubbling, slugging, turbulent and spouting fluidisation regimes, the following equations were used:

Onset of Bubbles in the Bed

Geldart and Abrahamsen (1978) proposed the following equation for the minimum bubbling velocity in the bed (Umb):

Onset of Slugs in the Bed

Stewart and Davidson (1967) give the minimum slugging velocity (Ums) according to the following expression:

where Umf is the minimum fluidisation velocity (m/s) and D is the bed diameter (m).

The conditions necessary for the onset of slugs in the bed are not always satisfied. Furthermore, it should be observed that the location of the slugging regime in the diagram is specific for a particular bed since Equation 8 indicates that the minimum slugging velocity depends on bed diameter.

Turbulent Fluidisation

Yerushalmi and Conkurt (1967) observed that for a certain superficial gas velocity denoted by Uc, the pressure fluctuations that occurred in the slugging bed are reduced and the large voids tend to disappear. Increasing even further the gas velocity to Urt, the pressure fluctuations fall to a value much lower than that corresponding to Uc. The frequency of the pressure oscillations becomes much higher and normally the large voids are no longer observed. For these velocities, the following correlations were obtained:

Spouting Regime

Geldart (1973) observed that this regime depends only on particle diameter and the difference between the solids and gas densities. The transition is given by

Equations 7 to 11 written in dimensionless form become

For fast fluidisation, there are no equations available. This regime depends on the rate of feeding the solids into the bed.

Dynamics of Gas Bubble Growth in the Bed

The fluid dynamics maps show that slugging of the bed may occur if the superficial gas velocity is higher than the minimum slugging velocity and lower than the value at which the turbulent regime occurs. However, two additional conditions are necessary for the onset of the slugging regime. The maximum stable bubble size for the system must be on the order of or greater than the inside diameter of the bed, and the height-to-diameter ratio of the bed must be large enough that slugs can form (Grace, 1982).

In this context, a study based on the dynamics of gas bubble growth in the bed was conducted, not only to determine the maximum stable bubble size, but also to evaluate the conditions under which these bubbles were stable.

The study of the dynamics of gas bubble growth was based on the upward velocity of the bubbles in the bed. The procedure described by Grace (1982) uses the following equations:

where Ub(H/2) is the bubble velocity at the midsection of the bed (m/s), db(H/2) is the bubble diameter at the midsection of the bed (m), dbo is the initial bubble diameter (m), dbm is the maximum stable bubble diameter (m), emf is the voidage of the bed under minimum fluidisation condition (-), fs is the particle shape factor or sphericity defined as (-), d is the mean volumetric fraction of the bed occupied by bubbles (-), H is the height of the bed (m), Hmf is the height of the bed under minimum fluidisation condition (m) and ms is the mass of solids (kg).

The calculation procedure consists of the following steps:

1. Mean volumetric fraction of the bed occupied by bubbles (d) is assumed;

2. The height of the bed (H) is calculated;

3. The bubble velocity (Ub(H/2)) is calculated;

4. The value d is verified by the equation,

If equation (23) is not satisfied, a new value of d is assumed and the procedure is repeated until it is satisfied. At this point, the values of d and H give estimates for the fraction of the bed occupied by bubbles and the expanded bed height, respectively.

RESULTS AND DISCUSSION

Analysis of the Fluid Dynamics Maps

Fluid dynamics maps of ZnS and ZnO particle beds fluidised respectively with air and a gas mixture containing 15% SO2 and 85% N2 at 883 and 1213 K were plotted. The temperatures are the extremes used in the kinetics study. The density of the blend of zinc sulphide concentrates is 4130 kg/m3.

The effect of superficial gas velocity on bed behaviour for 58mm ZnS particles fluidised with air at 883 K is shown in Figure 1. Particulate fluidisation occurs when the superficial gas velocity (U) is equal to the minimum fluidisation velocity (Umf) and bubbling fluidisation takes place for U = 2 Umf. When U = 20 Umf, the possibility of slugging of the bed exists. For U > 80 Umf, the particles enter the pneumatic transport regime.


The effect of temperature on the behaviour of beds of different particle sizes fluidised at a superficial gas velocity of 0.097 m/s is depicted in Figure 2. With increasing temperature, the positions of the beds shifted downwards and to the left with relation to the fixed/fluidised bed boundary. The same displacement occurs on the lines that define the several fluidisation regimes with the exception of the spouting bed boundary, which is not affected by temperature. It can be observed that the increase in temperature and the possibility of slugging concur, as evidenced by the 272 mm particles. Nevertheless, the influence of temperature is small and is restricted to beds that operate under conditions near those required for slugging flow.


The effects of gas composition and concentrate conversion at 1213 K and a superficial gas velocity of 0.097 m/s are illustrated in Figure 3 for beds of different particle sizes. This figure shows a comparison between the fluid dynamic behaviour of beds consisting of unconverted concentrate, fluidised by air, and fully converted concentrate, fluidised by a gas containing 15% SO2 and 85% N2. It is shown that displacement of the positions of the beds in relation to the fixed/fluidised bed boundary is proportional to displacement of the lines that define the fluidisation regimes. This result is justified since the changes in the gas physical properties are followed by the variation in concentrate density during the oxidation reaction. Thus, the modifications in gas and concentrate compositions do not significantly affect bed behaviour.


In order to ensure adequate conditions for mass transfer and good solids mixing, the bed should operate in the bubbling mode and not too far from the boundary with the particulate regime. This will cause the formation of small bubbles that rise in the bed, increasing gas-particle contact and minimising the bypassing of the reactant gas, improving solids mixing and enhancing mass transport. With this in mind, suitable superficial gas velocities and particle diameters for the kinetics study were obtained from the dashed area displayed in Figure 4. This area corresponds to superficial gas velocities between 0.0158 (Umf) and 0.0631 m/s and particle diameters in the 162-324 mm range. For particles smaller than 162 mm fluidised at a gas velocity of 0.0631 m/s, there is the possibility of slugging of the bed. Particles larger than 324 mm remain static (fixed bed).


Analysis Based on the Dynamics of Gas Bubble Growth in the Bed

The dynamics of gas bubble growth were evaluated for two initial masses of the zinc sulphide concentrate blend and a mean particle diameter of 165 mm. Figure 5 shows the results obtained for initial bubble diameter, bubble diameter at the midsection of the bed, maximum stable bubble diameter and height of the bed as a function of superficial gas velocity.


It can be observed in this figure that initial bubble diameter, bubble diameter at the midsection of the bed and maximum stable bubble diameter increase with increasing superficial gas velocity. The bubble diameters can become greater than the diameter of the bed. It is also observed that the increase in bed height is not proportional to the increase in bubble diameter. It can be seen that for an initial mass of 0.075 kg and a superficial gas velocity of 0.299 m/s there is the possibility of slugging, because the bubble diameter at the midsection of the bed is 0.034 m, while the diameter of the reactor is 0.045 m. For an initial mass of 0.045 kg and a superficial gas velocity of 0.235 m/s, the bubble diameter at the midsection of the bed is 0.020 m. For superficial gas velocities above 0.250 and 0.345 m/s and initial masses of 0.045 and 0.075 kg, respectively, the initial bubble diameter is greater than the bed height, which makes the bubble unstable.

An efficient temperature control and minimal harmful effects, i.e., slugging and/or particle losses produced by the increase in superficial gas velocity, are achieved under conditions where the bubbles are unstable (Grace, 1982). Therefore, initial diameter must be equal to or greater than bed height. Under this condition the bed is strongly agitated, the particles are projected against the wall by the bursting of the bubbles producing a radial gradient and the solids are concentrated around the wall of the reactor. The centre of the column begins to have a relatively small fraction of solids, which causes the breakup of bubbles with smaller initial diameters and higher frequencies. In this case, despite the low superficial velocity the characteristics of the bed are similar to those observed for the turbulent regime. The strong agitation of the bed provides better heat transfer between the solid and the wall of the reactor and mass transport between small and transitory voids, which minimises the bypassing of the gas, enhancing the gas-solid contact. In the present work this condition is achieved for superficial velocities of 0.256 m/s and 0.384 m/s when the initial masses are 0.045 kg and 0.075 kg, respectively.

Figure 6 shows the range of suitable experimental conditions for the kinetics study based on the dynamics of bubble growth in the bed.


CONCLUSIONS

Fluid dynamics maps and analysis of the dynamics of gas bubble growth allowed prediction of the behaviour of particle beds of a blend of zinc sulphide concentrates according to particle size, solid densities, fluidising gas properties and superficial gas velocity.

The fluid dynamics regime of the bed is mainly a function of particle diameter and superficial gas velocity and is subject to modification by an increase in temperature when the bed operates in the vicinity of the boundaries between the different regimes.

Based on the fluid dynamics maps, the ideal experimental condition to carry out the kinetics study was achieved when the bed was operated in the bubbling fluidising regime, close to the particulate/bubbling boundary. This condition corresponds to superficial gas velocities between 0.0158 (Umf) and 0.0631 m/s and particle diameters in the range of 162-324 mm. However, under this condition, agglomeration of the solids and difficulties in controlling the temperature may occur due to the highly exothermic reactions.

The study based on the dynamics of bubble growth allows determination of superficial gas velocities higher than velocities for which conventional studies predict slugging of the bed. For the present project the optimum experimental condition was achieved for superficial gas velocities of 0.256 m/s and 0.384 m/s when the initial masses fed into the reactor were 0.045 kg and 0.075 kg, respectively. These conditions make it possible to control the temperature without preventing gas-solid contact, thereby avoiding excessive particle losses.

ACKNOWLEDGEMENTS

The authors are grateful to the Brazilian government for providing financial support for this research and C. A. R. Queiroz thanks the Brazilian Research Council (CNPq), for the scholarship awarded.

Received: February 27, 2004

Accepted: August 18, 2004

  • Abrahamsen, A.R. and Geldart, D., Behaviour of Gas-Fluidised Beds of Fine Powders. Part 1. Homogeneous Expansion, Powder Technology, No. 26, 35 (1980).
  • Geldart, D., Types of Gas Fluidisation, Powder Technology, No. 7, 285 (1973).
  • Geldart, D. and Abrahamsen, A.R., Homogeneous Fluidisation of Fines Powders Using Various Gases and Pressures, Powder Technology, No. 19, 133 (1978).
  • Grace, J.R., Fluidised-Bed Hydrodynamics - Handbook of Multiphase Systems, Chap. 8, Hetsroni, G., ed., McGraw Hill, New York (1982).
  • Grace, J.R., High-Velocity Fluidised Bed Reactors, Chem. Eng. Sci., 45, No. 8, 1953 (1990).
  • Queiroz, C.A.R., Ustulação Oxidante de um Concentrado de Sulfeto de Zinco em Reator de Leito Fluidizado. Masters thesis, Pontifical Catholic University of Rio de Janeiro (1997).
  • Stewart, P.S.B. and Davidison, J.F., Slug Flow in Fluidised Beds, Powder Technology, No. 1, 61 (1967).
  • Yerushalmi, J. and Conkurt, N.T., Gas Backmixing in High Velocity Fluidised Beds, Fluidisation, Davidson, J.F. and Keairns, J.L., eds., Cambridge Univ. Press, Cambridge, 387 (1967).
  • *
    To whom correspondence should be addressed
  • Publication Dates

    • Publication in this collection
      14 Mar 2005
    • Date of issue
      Mar 2005

    History

    • Received
      27 Feb 2004
    • Accepted
      18 Aug 2004
    Brazilian Society of Chemical Engineering Rua Líbero Badaró, 152 , 11. and., 01008-903 São Paulo SP Brazil, Tel.: +55 11 3107-8747, Fax.: +55 11 3104-4649, Fax: +55 11 3104-4649 - São Paulo - SP - Brazil
    E-mail: rgiudici@usp.br