Follow:

Factors Affecting Mass Transfer Rates in Water-Lean Solvents

https://doi.org/10.1021/acs.iecr.0c00940

“The most straightforward way of evaluating the impact that each parameter has on mass transfer rates is by considering the penetration of CO2 into the solvent upon contact between vapor and liquid phases. (21,22) By doing this, we are assured that no significant phenomena (e.g., solvent depletion in the interface, shifting in reaction orders) will be overlooked.

Following this approach, CO2 concentration in the vapor–liquid interface is constrained by the thermodynamical equilibrium dictated by Henry’s law for dilute components. In the liquid bulk, diffusion of CO2 is followed by its reaction with amine molecules, forming carbamate and a protonated base according to Reactions R1 and R2. These reactions constitute the so-called zwitterion mechanism. (23) There is an ongoing debate in the literature regarding whether the zwitterion mechanism correctly represents the real phenomena behind the reaction between CO2 and amines. (24) Regardless, the rate equations obtained by employing the zwitterion mechanism can be shown to give results consistent with those derived by consideration of the one-step termolecular mechanism. (25,26) Furthermore, the reasons for choosing the zwitterion mechanism in our approach will be evidenced briefly in Section 2.3.

(R1)
(R2)

One can consider reactions R1 and R2 to derive an expression for the reaction rate of CO2. Typically, the strongest base in a solution will be the amine itself, and so this will be the preferable reactant, B, in reaction R2. Furthermore, the zwitterion is an unstable molecule, and its concentration can be approximated by zero at all times. This results in the reaction rate expression given by eq 1, where A is CO2, B is an amine, C is a carbamate, and D is a protonated amine (notice that CC = CD, so that CD will always be omitted in the following equations):

(1)

Equation 1 gives an expression for the irreversible reaction between CO2 and the amine forming carbamate and a protonated amine. For more general uses, however, it is interesting to consider that the conversion might be bounded by chemical equilibrium. This can easily be taken into account by modifying eq 1 into eq 2.

(2)

The reaction rate expressions for amine and carbamate come directly from stoichiometric relationships, so that rB = 2·rA and rC = – rA. The penetration equation for an arbitrary component i is then given by eq 3, where ri is the reaction rate of component i.

(3)

Substituting for A, B, and C results in eq 4a–c.

(4a)
(4b)
(4c)
For our algorithm, the expression related to the transport of the protonated amine is omitted. This is because no other charged product is created in our approach other than the carbamate and the protonated amine itself. Since both products are formed simultaneously (CC = CD), and since one has to impose electrostatic equilibrium to the solvent at all times, the diffusivity of the protonated amine and that of the carbamate must be equal, and an expression for one shall have the same form as the expression for the other. Therefore, eq 4c accounts for the transport of both products. Furthermore, an approximation taken in this work is assuming that the diffusivity of these products is similar to that of the amine itself (DC ≈ DB). This is in accordance with the approach adopted by several other studies. (21,27)

The resolution of eq 4a–c requires two sets of boundary conditions, one for the interface and another for the liquid bulk, along with one set of initial conditions. The initial conditions are ordinarily given by setting the initial concentrations of all components before absorption homogeneously throughout the whole liquid solvent, as seen in the following equations:

(5a)
(5b)
(5c)

In the interface, one can assume that the concentration of CO2 is in equilibrium with the partial pressure of CO2 in the vapor phase by Henry’s law and that there is no vaporization either of the amine or of the products. This results in the following equations:

(6a)
(6b)
(6c)

In the liquid bulk, far from the interface at x = δ, all of the concentrations are the same as those before the absorption took place. This can be seen in the following equations:

(7a)
(7b)
(7c)
In conclusion, the system of partial differential equations outlined in eq 4a–c subject to the initial conditions (eq 5a–c) and the boundary conditions (eqs 6a–c and 7a–c) has to be integrated from t = 0 to τ and from x = 0 to δ so that one can obtain the concentration profiles of CO2, amine, and carbamate in the solvent. The integration limits are the absorption time τ and the penetration depth δ.

As one can see, with the aforementioned assumptions and approximations, the CO2 mass transfer rates can be somewhat characterized by fixing two transport properties, two kinetic properties, and two thermodynamical properties. These are as follows.

  • Transport properties: CO2 diffusivity DA and amine diffusivity DB.

  • Kinetic properties: kinetic rate constant k2 and kinetic rate constant ratio kb/k–1.

  • Thermodynamical properties: equilibrium coefficient K and Henry’s coefficient HA.

These six parameters will be discussed further in the next sections. Before that, however, an explanation should be given about the equilibrium coefficient K. This property can be defined as by the following equation:

(8)

As mentioned previously, in this work, we are assuming that the only mechanisms through which CO2 is absorbed are physical absorption and carbamate formation, i.e., no other CO2-consuming reactions such as bicarbonate formation are taken into account. With this in mind, it is practical to discuss the equilibrium coefficient in terms of the loading α, which accounts for mols of CO2 absorbed per moles of amine in solution. As a function of α, the concentrations of CO2, amine, and carbamate are given, respectively, by eqs 910, and 11.

(9)
(10)
(11)

Substituting eqs 911 in eq 8, one ends up with eq 12. This expression shows how to obtain the equilibrium coefficient K by employing α.

(12)
Finally, two of the properties shown previously are not characteristic of the solvent but of the process or the algorithm itself. One of them is the penetration time τ. This property is equivalent to the surface renewal time scale, and in this study, it has been fixed at a constant value obtained using the Rocha et al. model (28,29) assuming an absorber working with structured packing of the type Intalox 2T (S = 0.0221 m, CE = 0.9) and an effective liquid-phase velocity of Ue = 0.5 m·s–1. This results in a surface renewal time of τ ≈ 0.05 s. The second property is the penetration depth δ. This is the depth in which no more concentration gradient can be seen as one goes deeper into the solvent, i.e., all concentrations are those of the liquid bulk. It is important to fix this parameter at a value big enough so that the whole penetration profile can be obtained, but not so big as to jeopardize the speed and convergence of the calculations.”

Leave a Comment