Chemical Engineers' Diffusion Study
Chemical Engineers' Diffusion Study
com
Abstract
Ternary diffusion coefficients of monoethanolamine (MEA) and N-methyldiethanolamine (MDEA) in aqueous solutions have been measured at
303.2, 313.2, and 323.2 K. The systems studied are aqueous solutions containing total amine concentrations of 2.5 and 4.0 kmol/m3 and each
solution prepared with four different amine molar ratios. The main diffusion coefficients (D11 and D22) and cross-diffusion coefficients (D12 and
D21) and the density and viscosity of the aqueous amine solutions are discussed and analyzed as a function of temperature and their concentration.
# 2008 Taiwan Institute of Chemical Engineers. Published by Elsevier B.V. All rights reserved.
Fig. 1. Mutual diffusion coefficients of aqueous MEA (2.00 kmol/m3) solution: Fig. 2. Ternary diffusion coefficients of Na2SO4 (1) + MgSO4 (2) + H2O at a
(*) data of Hikita et al. (1980); (&) data of Snijder et al. (1993); (~) this total salt concentration of 0.40 kmol/m3 at 298.2 K; lines, smoothed values.
study; line, calculated using equation of Snijder et al. (1993).
Table 2
method developed by Deng and Leaist (1991). For ternary Ternary diffusion coefficient of MEA (1) + MDEA (2) + H2O for a total amine
mutual diffusion coefficient measurement, the ternary mutual concentration of 2.50 kmol/m3
diffusion coefficients of Na2SO4 (1) + MgSO4 (2) + H2O for a T (K) Concentration Dij (109 m2/s)
total salt concentration of 0.4 kmol/m3 at 298.2 K were (kmol/m3)
measured. The results are presented in Table 1 and a c̄1 c̄2 D11 D12 D21 D22
comparison of the results is shown in Fig. 2. The x-axis
303.2 0.00 2.50 0.449
coordinate in Fig. 2 is the relative molar ratio of amines and was 0.50 2.00 0.506 0.004 0.070 0.464
defined as 1.00 1.50 0.604 0.032 0.061 0.477
1.50 1.00 0.696 0.073 0.029 0.511
c̄1
f1 ¼ (3) 2.00 0.50 0.795 0.113 0.014 0.535
c̄1 þ c̄2 2.50 0.00 0.975
where c̄1 and c̄2 are concentrations in the carrier stream. In this 313.2 0.00 2.50 0.576
figure, the four mutual diffusion coefficients (D11, D22, D12, and 0.50 2.00 0.681 0.151 0.057 0.555
1.00 1.50 0.807 0.154 0.088 0.614
D21), obtained from this study are generally in good agreement 1.50 1.00 0.908 0.165 0.061 0.642
with the data of Deng and Leaist (1991). 2.00 0.50 1.073 0.194 0.010 0.672
In Table 1 and Fig. 2, both D12 and D21 have negative values. 2.50 0.00 1.227
Deng and Leaist (1991) discussed these coupled phenomena for 323.2 0.00 2.50 0.727
MgCl2 + MgSO4 + H2O system using the limiting diffusion 0.50 2.00 0.956 0.009 0.121 0.729
coefficients of ions at infinite dilution in H2O, D0ik , and extended 1.00 1.50 1.137 0.097 0.029 0.729
it for Na2SO4 + MgSO4 + H2O system. The diffusion coeffi- 1.50 1.00 1.248 0.117 0.044 0.817
2.00 0.50 1.451 0.229 0.025 0.876
cient of Na+ ions is larger than that of SO42 ions. As a result,
2.50 0.00 1.536
an electric field is generated by the gradient in Na2SO4. To
648 C.-C. Ko et al. / Journal of the Chinese Institute of Chemical Engineers 39 (2008) 645–651
Fig. 3. Ternary diffusion coefficients of MEA (1) + MDEA (2) + H2O solutions Fig. 5. Main diffusion coefficients D22 of MEA (1) + MDEA (2) + H2O solu-
at a total amine concentration of 2.50 kmol/m3 at 303.2 K: (&) D11; (b) D12; tions at a total amine concentration of 2.50 kmol/m3: (&) 303.2 K; (b)
(") D21; (*) D22; lines, smoothed values. 313.2 K; (") 323.2 K; lines, smoothed values.
presented in Table 2. The corresponding mutual diffusion Also, at this limit, D21, the diffusion of MDEA due to the
coefficients in binary systems were also measured and are also concentration gradient of MEA, becomes zero since the
presented in Table 2. A plot of ternary mutual diffusion concentration gradient in MEA cannot produce a coupled flow
coefficients of MEA (1) + MDEA (2) + H2O versus the relative of MDEA, i.e., the solution is free of MDEA. Similarly, as
concentration of solutes at 303.2 K is shown in Fig. 3. The zero f 1 ! 0 (i.e., c̄1 ! 0), D22 is approaching the binary diffusion
values of Dik were also shown for easy observation of the coefficient in MDEA + H2O and D12 becomes zero since there
limiting case as c̄i ! 0. Fig. 3 shows that both D11 and D22 is no MEA in the solution.
increase as f 1 increases. The ratio of D12–D11 varies from 0.01 to 0.285 while the
In Fig. 3, it shows that at the limiting condition of f 1 ! 1 ratio of D21–D22 varies from 0.05 to 0.16. It can be concluded
(i.e., c̄2 ! 0), the main mutual diffusion coefficient, D11 that for the diffusion of MEA, the cross-diffusion effect due to
approaches the binary diffusion coefficient for MEA + H2O. the concentration gradient of MDEA is larger than that for the
diffusion of MDEA due to the concentration gradient of MEA.
In Fig. 4, the main diffusion coefficients, D11 of MEA
(1) + MDEA (2) + H2O solutions are shown as a function of
temperature at a total amine concentration of 2.5 kmol/m3. This
figure shows that D11 increases as the temperature increases at a
constant molar ratio and it also increases as the molar ratio
increases at a constant temperature. The main diffusion
coefficients, D22 of MEA (1) + MDEA (2) + H2O solutions
are shown in Fig. 5 as a function of temperature at a total
alkanolamine concentration of 2.5 kmol/m3. Similarly, D22
increases as the temperature increases at the same molar ratio
and also increases as the molar ratio increases at a constant
temperature. Fig. 6 shows the cross-diffusion coefficients D12
of MEA (1) + MDEA (2) + H2O solutions as a function of
temperature at a total amine concentration of 2.5 kmol/m3. In
this figure, as f 1 ! 0 (i.e., c̄1 ! 0), D12 becomes zero since
there is no MEA in the solution. In Fig. 6, D12 increases as the
temperature increases at the same molar ratio and also increases
as the molar ratio increases. Thus, at a constant temperature,
D12 increases when the amount of surrounding MDEA (a larger
Fig. 4. Main diffusion coefficients D11 of MEA (1) + MDEA (2) + H2O at a
molecule compared to MEA) molecules decreases.
total amine concentration of 2.50 kmol/m3: (&) 303.2 K; (b) 313.2 K; (") The ternary diffusion coefficients of MEA (1) + MDEA
323.2 K; lines, smoothed values. (2) + H2O for a total amine concentration of 4.00 kmol/m3 are
C.-C. Ko et al. / Journal of the Chinese Institute of Chemical Engineers 39 (2008) 645–651 649
Fig. 6. Comparison of cross-diffusion coefficients D12 of MEA (1) + MDEA Fig. 7. Comparison of main diffusion coefficients D11 of MEA (1) + MDEA
(2) + H2O for a total amine concentration of 2.50 kmol/m3: (&) 303.2 K; (v) (2) + H2O solutions at 303.2 K: (&) 2.50 kmol/m3; (*) 4.00 kmol/m3.
313.2 K; (") 323.2 K.
presented in Table 3. The dependence of the ternary diffusion to MEA + H2O, and the binary diffusion coefficient, D1, at
coefficients on the temperature and molar ratio of amine is 2.5 kmol/m3 is 0.975 109 m2/s and D1 at 4.0 kmol/m3 is
similar as in the total amine concentration of 2.5 kmol/m3 0.878 109 m2/s. This corresponds to the general behavior of
solution. the mutual diffusion coefficient in binary system, i.e., the
At 303.2 K, a plot of D11 as a function of various total amine diffusion coefficient increases as the concentration of solute
concentrations of 2.5 and 4.0 kmol/m3 is shown in Fig. 7. The decreases. Due to the presence of MDEA in solution, D11
main diffusion coefficient, D11 increases as f 1 increases for both becomes smaller than D1. This effect increases as the
total amine concentrations and D11 in a total amine concentration of MDEA increases. Thus, at a constant molar
concentration of 2.5 kmol/m3 is higher than that in a total ratio, D11 in a total amine concentration of 2.5 kmol/m3 is
amine concentration of 4.0 kmol/m3. At f 1 = 1, system reduces higher than that in a total amine concentration of 4.0 kmol/m3.
The plot of D22 as a function of different total amine
Table 3 concentrations at 303.2 K is shown in Fig. 8. This figure shows
Ternary diffusion coefficient of MEA (1) + MDEA (2) + H2O for a total amine
concentration of 4.00 kmol/m3
T (K) Concentration Dij (109 m2/s)
(kmol/m3)
c̄1 c̄2 D11 D12 D21 D22
303.2 0.00 4.00 0.309
1.00 3.00 0.311 0.026 0.092 0.346
1.50 2.50 0.361 0.065 0.054 0.344
2.50 1.50 0.514 0.108 0.064 0.368
3.00 1.00 0.593 0.117 0.033 0.424
4.00 0.00 0.878
313.2 0.00 4.00 0.398
1.00 3.00 0.432 0.017 0.150 0.449
1.50 2.50 0.565 0.043 0.099 0.464
2.50 1.50 0.711 0.112 0.069 0.496
3.00 1.00 0.828 0.130 0.044 0.542
4.00 0.00 1.108
323.2 0.00 4.00 0.509
1.00 3.00 0.694 0.051 0.199 0.566
1.50 2.50 0.727 0.072 0.101 0.553
2.50 1.50 0.921 0.083 0.065 0.618
3.00 1.00 1.038 0.113 0.059 0.721
Fig. 8. Main diffusion coefficients D22 of MEA (1) + MDEA (2) + H2O at
4.00 0.00 1.412
303.2 K: (&) 2.50 kmol/m3; (*) 4.00 kmol/m3.
650 C.-C. Ko et al. / Journal of the Chinese Institute of Chemical Engineers 39 (2008) 645–651
Fig. 10. Main diffusion coefficients of MEA (1) + MDEA (2) + H2O as func-
Fig. 9. Comparison of cross-diffusion coefficients D21 of MEA (1) + MDEA tion of viscosity of solution at a total amine concentration of 2.50 kmol/m3 and
(2) + H2O solutions at 303.2 K: (&) 2.50 kmol/m3; (5) 4.00 kmol/m3.
303.2 K: (4) D11; (&) D22; lines, smoothed values.
Table 4
that at the same amine molar ratio, D22 decreases as the total
Densities and viscosities of MEA (1) + MDEA (2) + H2O solutions
amine concentration increases—a similar behavior as D11. The
T (K) c̄1 (kmol/m3) c̄2 (kmol/m3) r (g/cm3) m (mPa s) plot of D21 as a function of different total amine concentrations
303.2 0.00 2.50 1.0249 2.6959 at 303.2 K is shown in Fig. 9. In this figure, the values of D21
0.50 2.00 1.0214 2.3469 show only slight difference between the total amine concentra-
1.00 1.50 1.0173 2.0050 tions of 2.5 and 4.0 kmol/m3; D21 in a total amine concentration
1.50 1.00 1.0132 1.7424
2.00 0.50 1.0093 1.5725
of 4.0 kmol/m3 is slightly higher than that in a total amine
2.50 0.00 1.0048 1.3291 concentration of 2.5 kmol/m3. Usually, D21, the diffusion
313.2 0.00 2.50 1.0204 2.0370 coefficient of MDEA due to the concentration gradient of MEA
0.50 2.00 1.0169 1.8037 (a smaller molecule compared to MDEA), is normally small
1.00 1.50 1.0131 1.5557 because of the size difference and do not vary appreciably with
1.50 1.00 1.0091 1.3884 the total concentration of amines.
2.00 0.50 1.0052 1.2636
In the Stokes–Einstein equation, the diffusion coefficient is
2.50 0.00 1.0014 1.1409
related to the viscosity of the solvent. In order to find out the
323.2 0.00 2.50 1.0210 1.6283
0.50 2.00 1.0119 1.4187
dependence of the diffusion coefficient on the viscosity of
1.00 1.50 1.0079 1.2782 solution, the viscosity and density of solutions have also been
1.50 1.00 1.0042 1.1654 measured, and the results are presented in Table 4. At a constant
2.00 0.50 1.0007 1.1245 temperature and total concentration of amines, both the density
2.50 0.00 0.9974 1.0058 and viscosity of the solution decrease as the concentration of
303.2 0.00 4.00 1.0405 6.6817 MEA increases. In Fig. 10, the main diffusion coefficients, D11
1.00 3.00 1.0335 4.8453
and D22, are plotted as function of viscosity of solutions. As
1.50 2.50 1.0300 4.1669
2.50 1.50 1.0222 3.0756 shown in this figure, both D11 and D22 decrease as the viscosity
3.00 1.00 1.0179 2.6597 of solutions increases; also D11 shows stronger dependence on
4.00 0.00 1.0095 2.0928 the viscosity than that of D22.
313.2 0.00 4.00 1.0344 4.6996
1.00 3.00 1.0278 3.4950 4. Conclusion
1.50 2.50 1.0244 3.0946
2.50 1.50 1.0169 2.3346
3.00 1.00 1.0130 2.0722
Ternary diffusion coefficients of MEA and MDEA in
4.00 0.00 1.0052 1.7321 aqueous solutions at a total amine concentrations of 2.5 and
323.2 0.00 4.00 1.0280 3.5042 4.0 kmol/m3 with four different molar amine ratios have been
1.00 3.00 1.0218 2.6434 measured at 303.2, 313.2, and 323.2 K using the Taylor
1.50 2.50 1.0186 2.3507 dispersion technique. For the total amine concentration, both
2.50 1.50 1.0120 1.8960 D11 and D22 increase as the temperature increases at a constant
3.00 1.00 1.0086 1.6192
amine molar ratio and also increases as the amine molar ratio
4.00 0.00 1.0005 1.3513
increases at a constant temperature. At a constant temperature
C.-C. Ko et al. / Journal of the Chinese Institute of Chemical Engineers 39 (2008) 645–651 651
and a total amine concentration, both D11 and D22 decrease as Hikita, H., H. Ishihawa, K. Uku, and T. Murakami, ‘‘Diffusivities of Mono-, Di-
, and Tri-ethanolamines in Aqueous Solutions,’’ J. Chem. Eng. Data, 25,
the viscosity of solution increases. From the ratios of Dik–Dii, it
324 (1980).
was found that for the diffusion of MEA, the cross-diffusion Hikita, H., H. Ishikawe, T. Murakaml, and T. J. Ishii, ‘‘Densities and Viscosities
effect due to the concentration gradient of MDEA is larger than and Amine Diffusivities of Aqueous MIPA, DIPA, DGA, and EDA
that for the cross-diffusion of MDEA due to the concentration Solutions,’’ J. Chem. Eng. Jpn., 14, 411 (1981).
gradient of MEA. The results of this study can be used in the Horng, S.-Y. and M.-H. Li, ‘‘Kinetics of Absorption of Carbon Dioxide into
kinetic study for the absorption of CO2 using the MEA Aqueous Solutions of Monoethanolamine + Triethanolamine,’’ Ind. Eng.
Chem. Res., 41, 257 (2002).
(1) + MDEA (2) + H2O as an absorbent. Kao, C.-Z. and M.-H. Li, ‘‘Diffusion Coefficients of Aqueous N,N-Diethy-
lethanolamine, N,N-Dimethylethanolamine, Monoisopropanolamine,
Acknowledgement Piperazine, and Sulfolane Solutions,’’ J. Chin. Inst. Chem. Engrs., 37,
341 (2006).
Kohl, A. L. and R. B. Nielsen, Gas Purification, 5th Ed., p. 40, Gulf, Houston,
This research was supported by the National Science
U.S.A. (1997).
Council of the Republic of China under grant NSC 96-2221-E- Leaist, D. G. and S. M. Abdu, ‘‘Ternary Mutual Diffusion Coefficients and
033-020. Critical Micelle Concentrations of Aqueous Sodium Dodecyl Sulfate + -
Lithium Dodecyl Sulfate Solutions at 25 8C,’’ J. Chem. Eng. Data, 46, 922
References (2001).
Leaist, D. G., Y. Li, and R. Poissant, ‘‘Coupled Diffusion in Aqueous Weak
Acid + Alkanolamine Absorbents,’’ J. Chem. Eng. Data, 43, 1048 (1998).
Alizadeh, A., C. A. Nieto de Castro, and W. A. Wakeham, ‘‘The Theory of the
Taylor Dispersion Technique for Liquid Diffusivity Measurements,’’ Int. J. Li, M.-H. and M.-D. Lai, ‘‘Solubility and Diffusivity of N2O and CO2 in
Thermophys., 1, 243 (1980). (Monoethanolamine + N-Methyldiethanolamine + Water) and in (Mono-
Austgen, D. M., G. T. Rochelle, and C. C. Chen, ‘‘Model of Vapor–Liquid ethanolamine + 2-Amino-2-Methyl-1-Propanol + Water),’’ J. Chem. Eng.
Data, 40, 486 (1995).
Equilibria for Aqueous Acid Gas–Alkanolamine Systems. 2. Representation
of H2S and CO2 Solubility in Aqueous MDEA and CO2 Solubility in Liao, C. H. and M. H. Li, ‘‘Kinetics of Absorption of Carbon Dioxide into
Aqueous Mixtures of MDEA with MEA or DEA,’’ Ind. Eng. Chem. Res., Aqueous Solutions of Monoethanolamine + N-Methyldiethanolamine,’’
Chem. Eng. Sci., 57, 4569 (2002).
30, 543 (1991).
Baldauf, W. and H. Knapp, ‘‘Measurements of Diffusivities in Liquids by the Rowley, R. L., M. E. Adams, T. L. Marshall, J. L. Oscarson, W. V. Wilding, and
Dispersion Method,’’ Chem. Eng. Sci., 38, 1031 (1983). D. J. Anderson, ‘‘Measurement of Diffusion Coefficients Important in
Chakravarty, T., U. K. Phukan, and R. H. Weiland, ‘‘Reaction of Acid Gases Modeling the Absorption Rate of Carbon Dioxide into Aqueous N-
Methyldiethanolamine,’’ J. Chem. Eng. Data, 42, 310 (1997).
with Mixtures of Amines,’’ Chem. Eng. Prog., 81, 32 (1985).
Chang, L.-C., T.-I. Lin, and M.-H. Li, ‘‘Mutual Diffusion Coefficients of Some Snijder, E. D., M.J.M.T. Riele, G. F. Versteeg, and W.P.M.V. Swaaij, ‘‘Diffusion
Aqueous Alkanolamines Solutions,’’ J. Chem. Eng. Data, 50, 77 (2005). Coefficients of Several Aqueous Alkanolamine Solutions,’’ J. Chem. Eng.
Chen, Y.-J., T.-W. Shih, and M.-H. Li, ‘‘Heat Capacity of Aqueous Mixtures of Data, 38, 475 (1993).
Taylor, F. R. S., ‘‘Dispersion of Soluble Matter in Solvent Flowing Slowly
Monoethanolamine with N-Methyldiethanolamine,’’ J. Chem. Eng. Data,
46, 51 (2001). through a Tube,’’ Proc. Roy. Soc. Lond. A, 219, 186 (1953).
Cussler, E. L., Multicomponent Diffusion, p. 66, Elsevier, Amsterdam, The Tyrrell, H. J. V. and K. R. Harris, Diffusion in Liquids, p. 21, Butterworth,
London, U.K. (1984).
Netherlands (1976).
Deng, Z. and D. G. Leaist, ‘‘Ternary Mutual Coefficients of Xiao, J., C.-W. Li, and M.-H. Li, ‘‘Kinetics of Absorption of Carbon Dioxide into
MgCl2 + MgSO4 + H2O and Na2SO4 + MgSO4 + H2O from Taylor Disper- Aqueous Solutions of 2-Amino-2-methyl-1-propanal + Monoethanolamine,’’
sion Profile,’’ Can. J. Chem., 69, 1548 (1991). Chem. Eng. Sci., 55, 161 (2000).