Absorption of sulfur dioxide in citric acid-sodium citrate buffer solutions

Theoretical models for the absorption process have been developed based on the film theory and on an ap- proximate surface renewal theory. It has been...
0 downloads 0 Views 633KB Size
Ind. Eng. Chem. Res. 1987,26, 1291-1296

1291

Absorption of SO2 in Citric Acid-Sodium Citrate Buffer Solutions Binay K. Dutta,* Ranjan K. Basu, Amit Pandit,+and Parthasarathi Ray Department of Chemical Engineering, Calcutta University, Calcutta 700009, India

Absorption of SO2 in citric acid-sodium citrate buffer solutions has been carried out by using a continuous flow stirred cell with a flat gas-liquid interface. The effects of gas concentration, liquid concentration, pH, and stirring speed on the absorption rate have been investigated. Theoretical models for the absorption process have been developed based on the film theory and on an approximate surface renewal theory. It has been observed that the absorption data could be satisfactorily interpreted using the model based on film theory of mass transfer. Removal of sulfur dioxide from stack gases is extremely important in the control of atmospheric pollution. The most widely used processes for this purpose are dry absorption in lime and limestone injected into the hot flue gas and the wet process where the flue gas is scrubbed in lime slurries, sometimes with additives that increase the rate of absorption. These processes, besides being nonregenerative in nature, create the problem of disposal of large amounts of waste. One regenerative process of sulfur dioxide removal involving absorption in citric acid-sodium citrate buffer solution has received considerable attention in recent times. The process was developed by Erga (1975) and later on by Farrington and Bengtsson (1979). The pH range used in this process is 3-5 so that the dissolved SO2 can subsequently be recovered from the loaded solution by steam stripping. Although some aspects of the process and the nature of solution equilibrium have been discussed by Erga (1980), the details are not available in open literature. Recently some studies on this absorption process have been reported by Dutta et al. (1986). The present work is intended to investigate the absorption characteristics of sulfur dioxide gas in sodium citrate-citric acid buffer solution using a continuous flow stirred vessel with a flat gas-liquid interface.

Experimental Setup and Procedure The absorption cell was made of a 10-cm section of 7.5-cm-i.d. Corning glass pipe with the ends closed by rubber corks. It was fitted with stirrers for both gas and liquid phases. The stirrer shafts were concentric but driven by different motors provided with speed changing devices. The rpm of the stirrers could thus be changed independently so that the effect of rpm in both gas and liquid phases could be studied individually. The liquid holdup in the absorption cell was 137 mL, and the liquid level was maintained constant by a liquid seal. The gas-phase stirrer was placed about half an inch over the liquid surface. Liquid samples for analysis were collected in two sampling tubes of known volumes provided with stopcocks and placed alternately in the exit liquid line. The stirred cell with the sampling device is shown schematically in Figure 1.

The gas phase was a mixture of SO2 and air, prepared by mixing these gases at desired proportions. SO2 was available in cylinder, and air was supplied by a compressor. The individual gas flow rates were measured by calibrated capillary flow meters. The absorbent liquid was taken in a large Morriotte bottle, and the flow rate was measured Present address: Department of Chemical Technology, Bombay 400019, India. 0888-5885/87/2626-1291$01.50/0

by another capillary flow meter. The flow rates of the various streams were adjusted by fine-control needle valves. The system was continuous with respect to both gas and liquid phases. The gas flow was started after the liquid level could be adjusted and maintained at the desired position. The attainment of steady state was ascertained by taking liquid samples at different times during a run and analyzing the SO2content in the samples. It usually required about 15 min for attainment of steady state. However, in all the runs, absorption rates were found out by taking liquid samples after the system operated under uniform conditions for more than 25 min. Liquid flow rate was 45 mL/min and gas flow rate was maintained at 8 L/min. It was found that a maximum of only 2% of the SO2entering was absorbed in the liquid at steady operating conditions. It was also found by analysis of several gas samples that the inlet and exit gas concentrations were virtually equal. Absorption rates were measured in water as well as in sodium citrate-citric acid buffer solutions at different gas and liquid concentrations and pH. Effect of stirrer speed in both of the phases was also studied independently. All the runs were taken at a constant temperature of 35 f 0.5 OC and essentially atmospheric pressure. Fluctuations in flow rates and gas concentration during any run were below 2 % . Several runs were repeated to ensure reproducibility of the experimental data collected. Oxidation of SOz in solution could not be detected. The concentration of dissolved gas was estimated by iodometric titration. pH values of fresh and loaded buffer solutions were measured by Phillips Model PM-9045 digital pH meter.

Theoretical Considerations When SO2is absorbed in a polyfunctional buffer solution such as sodium citrate-citric acid buffer, the following ionization reactions occur. SOp + HzO + H+ + HS03-

(1)

HS03- 6 H+ + S032-

(2)

H,Ci + H+ + H2Ci-

(3)

H2Ci- + H+

+ HCi2-

(4)

HCi2- 6 H+ + Cis(5) Where H3Ci, H2Ci-, HCi2-, and Ci3- represent citric acid and corresponding anions. Reactions 1and 2 only occur if the absorbent liquid is pure water. The hydrolysis of dissolved SO2 (eq 1)is known to be very fast (Eigen et al., 1961),the proton-transfer reactions 2-5 being even faster. Thus, all the above reactions may be regarded as instantaneous. 0 1987 American Chemical Society

1292 Ind. Eng. Chem. Res., Vol. 26, No. 7, 1987

Boundary conditions 11and 12 indicate that there is no net flow of the citric acid or citrate ions and of the charge across the gas-liquid interface. Also, as x S = [SI, (13) In addition, the concentrations of the individual species in eq 6-8 are interrelated by the equilibrium relations which, for eq 1-5, may be expressed as

- Srirred Cell 2 , 3 - G a s inlet 8 outlet 4, 5 - L i q u i d inlet a outlet 6 - Thermometer 7,8- Liquid a G a s stirrers 9 - Static h e a d meter 10,Il- Sampling tubes

-

K,' =

~H+~HSO,-

aH+aS032-

K,' = Figure 1. Experimental setup-stirred

Dso*

d2[S021 ax*

~

d2[HS03-] + DHSOj

+ Dso32-

8x2

a[SO,] at

d2[so32-]

+

d[HS03-] at

d[S03*-] (6)

For H3Ci and corresponding anions d2[H3Ci]

a2[HzCi-]

7

DH,Ci

+DbCI

axL

a[H2Ci-]

YH+YS032-

aH+aH2Ci7

YH3Ci

+

aH+aHCi2-

a2[H+]

___

YHCi2-

uHCiZ- YH'YCi"

a21~0~2-1 -

2DS0,J

d2[H2Ci-] - 2DHC~2

ax2

where K,-K6 are the respective thermodynamic equilibrium constants and K,'-K,' are the effective equilibrium constants in concentration units. It is not possible to obtain analytical solution for the above coupled set of partial differential equations (eq 6-8) subject to prescribed conditions. Numerical solution would also be quite involved. The corresponding film theory equations will,however, be much simpler and would admit of closed-form solutions. The film theory equations are

dX

[SO,] = [SO,],

+ DH2C~

a [H2Ci-] ~

dX

(10)

+

a[HCi2-]

DHC~2-

aw+1

DH+ -- DHS03ax

d[H2Ci-] DHzC~

dX

7 -

-- 0 (11) dX

- 2Dso32

a [ HCi2-]

2 D ~ c ~___ 2 ax

dx2

-

a [so32-] ~

3%

ax

-

d[Ci3-]

7 =o (12)

dx2

= 0 (15)

+ d2[Ci3-]

+ Dci3- ___ dx

= 0 (16)

and d2[H+] DH+

d2[HS03-]

-- DHS03-

dx2

d2[S032-] -

2DSo32-

-

dx2 dx2 d2[H2Ci-] d2[HCi2-] - 2D~ci2DH2C~ dx2 dx2 d2[Ci3-] 3Dci3= 0 (17) dx The above equations assume that there is no potential gradient across the film. The boundary conditions at x = 0 (gas-liquid interface) remain the same as eq 10-12. The boundary condition 13 in this case reduces to ~

x =6

a[Ci3-]

+ Dc,3 7

d[HSO3-]

dx2 d2[ HCi*-]

dX2

[SI = [SI, (9) where [SI is the concentration of any species in eq 6-8. Also

~

dx2

d2[H2Ci-] + DH2Ci?

DHCi2-

t =0

a [H3Ci]

DH3Ci

a* [ HCiL-] -

The appropriate initial and boundary conditions are

DH?C~

d2[S032-]

Dso32d2[H3Ci]

-

d2[Ci3-] d[H+] d[HS03-] d[S03*-] -2-3DCl ____ - ~ax 2 at at at d[H,Ci-] d[HCi2-] a[CiJ-] ~2 ~ 3 (8) at at at

x =0

(14)

YH+YCi3-

d[HCi2-]

a*[HSO3-]

dX'

YHCi2-

= K5-

~

YH2Ci-

YH+YHCi2-

+-+- at

-- DHSOj ax2 DH,Ci

= K4

YH+YHCi2-

aH+aci3-

YH&i

YH+YH2Ci-

YHzCi-

d[Ci3-] (7) at at In addition, the following equation must be satisfied for electrical neutrality of the solution. DH+

YH+YS032-

YH'YH2Ci-

~

K5' =

YHS03-

= K3

~

uH2Ci-

-

ax2

+---- at

K4r =

YH+YH S O ~ -

= K2

QHS03-

aH3Ci

In order to apply the surface renewal theory, the unsteady-state diffusion equations for this system may be expressed as follows. For SO2and corresponding anions

YHS03-

___

cell.

Yso2

= K1

YH+YHSO~-

%02

K2' =

Yso2

s = [SI,

(18)

Concentration distribution of all the species concerned may be obtained by solving eq 15-17 subject to boundary conditions 10-12 and 18 together with the equilibrium relations 14. However, calculation of the interfacial concentrations of HS03- and SO3- are particularly necessary for finding out the theoretical enhancement factor for absorption. The interfacial concentration of hydrogen ion, [H+Ii,is required for this purpose. It can be shown that

Ind. Eng. Chem. Res., Vol. 26, No. 7, 1987 1293 [H+Iimay be calculated by solving the sixth degree algebraic equation (more details are given in the Appendix) [H+]: CU~[H+]: a2[H+]t C Y ~ [ H++ ]a4[H+]i2 ~ + aS[H+]i (Ye = 0 (19)

+

+

+

Liquid i t l r r e r r . p . m .

I-

o

+

where Dd(3/- b3 = (D7K4 - b2D3 - b3Ds)Ki - D1K1’ [SOJi as = (D&,’ - 2b2D4 - b3D7)KiK; (DlD&i + 2D2K2/)KI’ [SO2li = -(DiD7K: + 2020&2/)K1’Ki[S02]i (3b2D5 + b3DB)KiK:K5/ a 5 = -(D1D$(5/ + 2D,D,K,’)K,’K3’K: [SOzli “6 = -2D2DsK1’K2/K3’K4‘K5/[SO2]i (20) b2 and b3 are given by b2 = [H3Ci], D6[H2Ci-I0 D7[HCi2-Io+ DB[Ci3-l0

+

+

b3 = [H+l0- D1[HSO3-lo- 2D2[S032-]o- D3[H2Ci-Io2D4[HCi2-Io- 3[Ci3-], (21)

9

c

a

-.-

I 0

200 Go8

400 phose

600

800

1000

rtiner r,p.m.

Figure 2. Effect of gas-phase stirrer rpm on absorption rate.

c0 x

2

-

Temp: 3 5 ‘ ~ Cos stirrer r,p.m. : 266 0 4 . Liquid stirrrr r.p.m.: a90 V70 A50

The enhancement factor for the absorption process may now be readily calculated by using the film theory.

Ecf = 1 + (DHs~,-([HSO~-]~ - [HSO,-]o) + Dso3~([S032-Ii - [S032-Io))/lD~~,([S021i - [SO2IO)l (22) It was pointed out by Olander (1960) that for single equilibrium reaction the film and surface renewal theory expressions for enhancement factor are identical if the diffusivities of the species are equal. In other cases, an approximate expression for the enhancement factor was suggested by Chang and Rochelle (1981,1982a) which is the same as the film theory expression except that the diffusivity ratios are replaced by the corresponding square roots. For the present system, the enhancement factor for approximate surface renewal theory would be E,, = 1 + (D~so,-’/~([HSo,-]i - [HSO,-]o) + D s o ~ z - ” ~ ( [ S O-~[S032-l~))/(D~,,l/2([sO~li ~-I~ - [SOZIO)) (23) For the case of absorption in pure water, the expression for enhancement factor for the film theory given by Olander (1960) may be used.

The same expression for the surface renewal theory may be approximated by (Chang and Rochelle, 1981)

where the second ionization of the dissolved SO2 is neglected, as it occurs to a negligible extent in solution in pure water. Equations 22-25 may be used to test the chemical absorption mechanism of SO2.

Experimental Results Absorption rates were measured in sodium citratecitric acid buffer solution as also in pure water at different gas and liquid stirrer speeds, gas concentrations, liquid concentrations, and pH. The fraction of the feed gas absorbed at steady state was very small, and the gas concentration could thus be calculated from observed flow rates of air and SO2.

00

h

IO

I2

14

16

1.8

20

Gar c o n c ( v o l % S o 2 1

Figure 3. Effect of gas concentration on absorption rate in water.

Test for Gas-Phase Resistance. For studying gas absorption with chemical reaction, it is convenient to eliminate gas-phase resistance to mass transfer. The significance of gas-phase resistance for the system under study was tested by taking several runs with different gas-phase stirrer speeds while keeping all other system variables fixed for SO2-buffer and SO2-water systems individually. No effect of increasing stirrer rpm from 255 to 997 on the rate of absorption was observed (Figure 2). Thus, even at a speed of 255, the gas-phase stirrer could produce sufficient turbulence in the medium so that the resistance to mass transfer in that phase was negligible. Effect of Different System Variables on Absorption Rate. The specific absorption rates of SO2 in both water and buffer solutions were found to increase linearly with gas concentration (Figures 3 and 4). This is because of increased solubility of the gas (and hence of interfacial concentrations of bisulfite and sulfite ions) with increased partial pressure of SO2in the gas phase. The effect of pH on absorption rate is shown in Figure 5. It is found that the absorption rate increases with pH to some extent at each stirrer speed. The effect of concentration of the feed solution was studied by preparing solutions with different amounts of sodium citrate and citric acid so that the pH remained constant. Figure 6 shows that the absorption rate is nearly independent of feed concentration. This is because at low concentrations of gas used and at the interfacial area per unit volume of liquid available, the observed concentration of the dissolved gas was too small to change the pH of the absorbent substantially even for the case of lowest solution concentration (3.65%). However, if the amount of dissolved sulfur dioxide in the liquid is large, the absorption rate is likely to be more for solutions of higher buffer concentration because of correspondingly larger buffer capacity.

1294 Ind. Eng. Chem. Res., Vol. 26, No. 7, 1987 2.5

T e m p : 35.C Buffer ccnc: 0.178(M) 0 0 s stlrrrr r . p . m . : 2 5 5 o f f e e d solution: 4.42

2.0

Liquid stirrer r . p . m.:

-

e 90 W

001 0.0

I

I

1

I

0.5 1.0 1.5 Gos conc. ( V O I .% so, )

2.0

0.0

Figure 4. Effect of gas concentration on absorption rate in buffer solution. Temp: 35.C

L a

I'

0.0

3.0

1.5

2.0

2.5

(E, /E, )exp Figure 7. Theoretical and experimental enhancement ratio plotfilm theory. SO2 concentration, 1.05-2.00%; pH, 3.58-4.80; solution concentration, 3.65-12.16%;liquid rpm, 50-90.

The physical solubility, however, is likely to be influenced by other species present in the solution, and a method of estimation of the magnitude of such effect was suggested by van Krevelen and Hoftijzer (1948). As the salting out parameters for the different ionic species present were not available, eq 28 was used to estimate the physical solubility of SO2. Diffusivity of SO2 and Other Species. Diffusivity of molecular SOz in water was reported to be 2.00 X cmz/s a t 30 "C (Peaceman, 1951). The diffusivities of H+, HS03-, and SO:- ions were estimated by using the Nernst equation (Robinson and Stokes, 1959) RTXo D=(29)

I 3.5 4.0 4.5 P" of feed Solution

5.0

Figure 5. Effect of pH on absorption rate. Temp: 35'c Goa conc :I.70%SO2 Gas stirrer r.p.m.: 255 of f o r d solullon: 4.42 Liquid stirrer r. p.m. I

1.0

Physical Solubility of SO2. Johnstone and Leppla (1934) reported the ionization constant and solubility in SOz-water system. Rabe and Harris (1963) reported a detailed investigation of the equilibrium and solubility for this system. They recommended the following expression for the Henry's law constant as a function of temperature: H = exp[2851.1/T - 9.37951 (28)

Goa conc: 1.70% SO2 & p s sh'rrer r . p . m . 2 5 5 Liquid stirrer r . p m .

*-

0.5

70

vj32

Y--A

*

a

0.oo

4

8 12 Feed liquid conc.

('/e

16 20 by W t . )

Figure 6. Effect of feed liquid concentration on absorption rate.

Interpretation of Absorption Data The theoretical enhancement factor ( E )for absorption in both SO2-water and SO2-buffer systems may be calculated by using eq 22 and 24 for the film theory and eq 23 and 25 for the approximate surface renewal theory, respectively. By comparison of the specific rates of absorption in S02-buffer (kbuffer) and in SOz-water (kwater) kbuffer kwater

-kbuffer kwater

--

- -Ecf

for film theory

Ewf

for surface renewal theory

(26) (27)

Ews

Equations 26 and 27 provide simple but effective ways of testing the validity of the theories for the system studied. Different physical properties and data required for this purpose have been taken from literature or estimated by using suitable correlations.

where 3 is the Faraday number and vj is the charge of the species. The infinite dilution ionic conductivity, Xo, for the above species was taken to be 350, 50, and 72 cm2/ (ohm-g-mol),respectively (Landolt-Bornstein, 1960). Some experimental data are available on diffusivity of H3Ci and H2Ci-. Muller and Stokes (1957) reported these values a t 25 "C and at different acid concentrations. No data are available on the diffusivities of HCi2- and Ci3-. The values for H&i- and HCi2- were, however, assumed to be equal following Albery et al. (1967) who observed that the diffusivities of dibasic acid anions were only slightly affected by the ionic charge. The same value was also taken for Ci3-, as the third ionization of H3Ci occurs to a very little extent in the pH range used in this study, and even an appreciable error in the diffusivity of Cis will not have any significant effect on the computed results. The diffusivity values were corrected for change in temperature and viscosity wherever necessary by using the Stokes-Einstein relation. The values for H3Ci and HzCiused were 0.839 X and 1.035 X cm2/s, respectively. Ionization Equilibrium Constants. The thermodynamic equilibrium constant for the ionization reaction 1 was reported by Rabe and Harris (1963). K1 = exp[1972.5/T - 10.9670] (30) The equilibrium constant for reaction 2 was reported by at 25 Tarter and Garretson (1941) to be K z = 6.2 x

Ind. Eng. Chem. Res., Vol. 26, No. 7 , 1987 1295

/

Temp: 35'c

,.*. ' '

c

a -

0

1

2

pH and the higher the degree of dissociation attainable. The behavior is also in conformity with the experimental observations.

3

4

5

(Ec/h)e\p

Figure 8. Theoretical and experimental enhancement ratio plotsurface renewal theory. SO2 concentration, 1.05-2.00%; pH, 3.584.80; solutioa concentration, 3.65-12.1670; liquid rpm, 50-90.

OC. As the second ionization of dissolved SO2 was very small in the pH range of this study, this value was used in the calculations. The equilibrium constants for reactions 3-5 are available as functions of temperature (Bates and Pinching, 1949).

PK3 = 1255'6 - 4.5635 + 0.01167T T PK4 = 1585*2- 5.4460

T

+ 0.01640T

PK5 = 1814'9 - 6.3664 + 0.02239T (31) T Activity Coefficients of the Species in Solution. Activity coefficient of unionized H3Ci was obtained from the data provided by Muller and Stokes (1957) at different acid concentrations, and those of citrate ions were given by Conway (1952) at different total ionic concentrations. Activity coefficient of hydrated SOzwas estimated by the relation (Harned and Owen, 1958) log yi = 0.0761 (32) where I is the ionic strength of solution. The ionic activity coefficients of H+, HS03-, and were estimated by using the modified Debye-Huckel limiting law (Lowell et al., 1970). c

The values of CBiand Cqiavailable a t 25 OC were used (Chang and Rochelle, 1981). Comparison of Theoretical and Experimental Enhancement Factor Ratios. The theoretical enhancement factors for both SOz-buffer and SOz-water systems can be calculated by using the above equations and data available. The experimental value of the enhancement ratio is the ratio of the rates of absorption in SOz-buffer and SOz-water systems, respectively, under identical conditions. The theoretical enhancement ratios using both film theory and approximate surface renewal theory are plotted against the experimental values in Figures 7 and 8, respectively. It appears from these parity plots that the experimental data could be better represented by the film theory of mass transfer with multiple equilibrium reactions than by the approximate surface renewal theory. Computed results also indicate higher enhancement at lower feed gas concentration, because the less the interfacial concentration of dissolved sulfur dioxide, the more is the

Conclusions Absorption of SOz in low concentrations using citric acid-sodium citrate buffer solutions exhibits considerable enhancement over absorption in water. Absorption rate is weakly dependent on pH in the range 3-5 but is independent of the concentration of the feed solution. The rate, however, increases linearly with the gas concentration. Absorption data could be better interpreted by the film theory of mass transfer with multiple equilibrium reactions than by the approximate surface renewal theory. Nomenclature u1-a3, b,-b3 = integration constants, eq 34-36 a = activity D = liquid-phase diffusivity DI = Dnso3-ID~+ Dz = DSo32-/DH+ D3 = DH2Ci-/DH+ D4 = DHCiz-/DH+ D, = DCp/DH+ DE= DH2Ci-/DH3Ci 01= DHCi2-/DH3Ci D, = DcPIDH~c~ E = enhancement factor H = Henry's law constant, eq 28 HnCi3-" = citric acid and corresponding anions (n = 0-3) K = thermodynamic equilibrium constant K' = effective equilibrium constant in concentration units kbuffer= rate of absorption in SOz-buffer system k,,,, = rate of absorption in SOz-water system R = universal gas constant S = species in solution T = absolute temperature t = time x = coordinate Greek Symbols y = activity coefficient 6 = film thickness

Subscripts c = citric acid-sodium citrate buffer solution f = film theory i = gas-liquid interface 0 = bulk solution s = refers to surface renewal theory w = water

Appendix Derivation of Equations 19 and 20. Integration of eq 15-17 results in

DH2Ci-

DHCi2-

DH3Ci

DHBCi DCi3---[Ci3-] DH3Ci

[H3Ci] + -[HzCi-] + -[HCiz-] +

Dci33-[Ci3-] DH+

= a2x

+ b2 (35)

= u3x

+ b,

(36)

1296 Ind. Eng. Chem. Res., Vol. 26, No. 7, 1987

Differentiating eq 35 and 36 with respect to x and comparing with eq 11 and 12, a2 = 0 = a,. Then, DH&!-

DHCi2-

[H3Ci] + -[H2Ci-] + -[HCi2-] + DH4Ci DH9Ci Dci3-[Ci3-] DH3Ci

Dei K3’K4’Kj’

= b2 (37)

DHSOL Dso32[H+] - -[HS03-] - 2-[S032-] DH+ DH+ Dci3DH2CiDHCi2[Ci3-] = b, (38) -[H2Ci-] - 2-[HCi2-] - 3DH+ DH+ DH+

Thus, the quantities on the left-hand side in eq 37 and 38 are invariant over the thickness of the film, and it may be noted at this point that b2 and b3 can thus be calculated from the bulk concentrations of the different species involved. At the gas-liquid interface, eq 34-36 may be written as Dso32[HSOa-1, + -[SO~2-ll Dso* [H3Ci],+

DHzC~-

DHCi2-

DHX~

DHICi

= b, (39)

Dci3-

-[Ci3-] = b, (40) DH8Ci

and

Dso32-

DHS03-

DH2Ci-

-[H2Ci-Ii DH+

-

Dci3-

DHCi2-

2-[HCi2-Ii

- 3-[Ci3-Ii

DH+

DH+

= b3 (41)

Now the interfacial concentrations of the various species in solution may be related by using the effective equilibrium constants [eq 141 as

[H2Ci-Ii=

K3’ [H3CiIi [H+Ii

[HCi2-Ii= [ci3-li=

K3’K4’[H3Ci]i

.[H+Ii2

K3K4’K5‘[H3CiIi

W+l?

i3

(43)

Also, from eq 41 and eq 42, [H+]i- DHSO~K1’[S02]i- 2- Dso32-KI’K2’[S02ji DH+ [H+]i DH+ [H+li2 DH2Ci-K,’[H3CiIi DHCiz- K3’K4’[HaCiIi -. . - 2DH+ [H+]i DH+ [H+li2 Dcia-K3’K,’K,’ [H3Cili 3-= b3 (44) [ H+]i3 DH+ On substituting eq 43 in eq 44 and simplifying, we get eq 19 with a’s given by eq 20. The principles underlying the theoretical development are similar to those mentioned by Chang and Rochelle (1982b). Registry No. SOz, 7446-09-5; citric acid, 77-92-9; sodium citrate, 994-36-5.

Literature Cited

-[H,Ci-], + -[HCi2-],+

[H+]i- -[HSO,-]i - 2-[S03’-]i DH’ DH+

----)

DH3Ci3- [ H’]

(42)

On substitution of eq 42 in eq 40 and on simplification,

Albery, W. J.; Greenwood, A. R.; Ribble, R. F. Trans. Faraday S o c . 1967, 63, 360. Bates, R. G.; Pinching, G. D. J . Am. Chem. Soc. 1949, 71, 1274. Chang, C. S.; Rochelle, G. T . AZChE J . 1961, 27, 292. Chang, C. S.; Rochelle, G. T . Ind. Eng. Chem. Fundam. 1962a, 21. 379. Chang, C. S.; Rochelle, G. T . AZChE J . 1962b, 28, 261. Conway, B. E. Electrochemical Data; Elsevier: Amsterdam, 1952. Dutta, B. K.; Basu, R. K.; Pandit, Amit; Ray, P., Paper presented a t the World Congress I11 of Chemical Engineering, Tokyo, Sept 1986. Eigen, M.; Kustin, K.; Maas, G. 2. Phys. Chem. (Frankfurt)1961, 30, 130. Erga, 0. U.S. Patent 3 886060, 1975. Erga, 0. Chem. Eng. Sci. 1980, 35, 162. Farrington, J. F.; Bengtsson, S. 108th AIChE Meeting, 1979, New Orleans; Paper A-79-35. Harned, H. S.; Owen, B. B. The Physical Chemistry o j Nonelectrolytic Solutions; Reinhold: New York, 1958. Johnstone, H. P.; Leppla, P. W. J . Am. Chem. Soc. 1934, 56, 2233. Londolt-Bornstein Physihalish-Chemische T a b r l k n , Bd 11-7, Springer-Verlag: Berlin, 1960. Lowell, P. S.; Ottmers, D. M.; Schwitzgebel, K.; Strange, T. I.; Deberry, D. w . A.P.T.D. 1287, P B 193-029, 1970; USEPA, Washington, DC. Muller, G. T. A,; Stokes, R. H. Trans. Faraday due. 1987, 58, 64%. Olander, D. R. AZChE J . 1960, 6 , 233. Peaceman, D. W. Sc.D. Thesis, MIT, Cambridge, 1951. Rabe, A. E.; Harris, J. F. J . Chem. Eng. Data 1963, 8 , 333. Robinson, R. A.; Stokes, R. H. Electrolytic Solutions; Butterworths: London, 1959. Tarter, H. V.; Garretson, H. H. J . Am. Chem. Soc. 1941, 63, 808. Van Krevelen, D. W.; Hoftijzer, P. J. “Chimie et Industrie, Numero Speciale du XXI e”, Presented a t the Congress International de chiemie Industrielle, Bruxelles, Sept 1948. Received for reuieu J u n e 3 , 1985 Revised manuscript received July 17, 1986 Accepted February 24, 1987