Abstract

Gas holdup is an important parameter in the carbonation reaction for processing red mud. In this paper, for estimating the gas holdup of Venturi carbonation reactor, water model experiments were performed. The effects of superficial gas velocities, superficial liquid velocities, and liquid-to-solid and height-to-diameter ratio on the gas holdup are investigated. Through dimensional analysis using Buckingham π principle, a derivation of the empirical correlation is also proposed. The results indicate that gas holdup increases with the increase in gas velocity, with the decrease in liquid velocity, with the increase in liquid-to-solid ratio, and with the increase in height-to-diameter ratio. The calculated results from the empirical correlation agree well with the experimental data, which is important for designing a highly efficient Venturi reactor with a high temperature, high pressure, and three phases of gas, liquid, and solid.

1. Introduction

Red mud is a solid residue, which is generated during alumina extraction from bauxite. According to the properties of bauxite and alumina production conditions, 0.6–2.5 tons of residue were generated for every ton of Al2O3 during alumina production [1]. Approximately 120 million tons of Bayer red mud were generated every year worldwide, and substantial reserves of over 2.7 billion tons have accumulated [2]. Every alumina plant allocates a large yard for stacking the Bayer red mud. The industrial solid waste encroaches a great amount of land, causes air pollution, contaminates soil or groundwater, and results in serious environmental problems because of the tiny particles and alkaline slurry. The utilization and recycling of red mud is currently a critical issue and needs to be addressed as soon as possible.

Researchers have made significant efforts to explore the utilization of red mud in aluminum industry. Some of them involve the extraction of valuable materials from red mud, such as the recovery of alumina, soda, ferric oxide, and titanium oxide [35]. Most of the red mud is still stacked in yard, although a significant amount of red mud is utilized as subgrade materials [6, 7] and additive for cement [8] and ceramics [9], and for various other industrial uses [10, 11]. The biggest hurdle in the use of red mud is the composition of soda in red mud slurry. Soda needs to be removed prior to utilizing the red mud in industrial processes.

Although many processes have been proposed for the utilization of red mud, only a few have been applied in industrial processes due to unfavorable local economic and market conditions. It is worth mentioning that most of these processes suffer from high costs and low yields. For fully using the Bayer process residue, the team of Northeastern University, China, proposed a novel approach, called the calcification-carbonation method [12, 13]. The method mainly involves two steps of calcification and carbonation. In the first step of calcification, lime is added to red mud slurry in an autoclave at the temperature of 180–240°C, whereas the calcification reaction occurs in the solution, during which the silicon phase is transformed into hydrogarnet (3CaO·Al2O3·xSiO2·(6 − 2x)H2O). As a result, soda in red mud is extracted:

In the second step of carbonation, carbon dioxide (CO2) is bubbled into the calcified red mud slurry in the carbonation reactor at 80–120°C under a CO2 partial pressure of 0.5–1.2 MPa. The hydrogarnet is decomposed into calcium silicate, calcium carbonate, and aluminum hydroxide, whereas alumina is recovered during the process.

Based on this method, 49% of alumina and 96.8% of soda are extracted from red mud. The new-phase red mud, which mainly consists of calcium carbonate, calcium silicate, and trace amounts of Na2O, is obtained. The new-phase red mud satisfies the requirements for the production of Portland cement and, therefore, can be used fully.

The main unit operation for the process is the carbonation reactor. A Venturi carbonation reactor was adopted for the carbonation reaction, which is a type of self-stirring reactor with a cylindrical vertical cross section and a conical bottom equipped with a Venturi tube (Figure 1). The carbonation reaction is a three-phase reaction consisting of gas, liquid, and solid. The ascending CO2 gas, which is laterally introduced near the bottom, disperses into the red mud slurry within the reactor. The carbonation reaction is performed during the rising of gas and slurry. The efficiency of the carbonation reaction is governed by the concentration of gas in liquid and the contact surface of bubbles with solid. Generally speaking, in gas-liquid-solid flow systems, this contact mainly depends on several factors, including the bubble sizes and their distributions [14, 15], gas holdup [16, 17], and bubble rise velocity [18, 19].

Gas holdup (ε) is an important parameter in hydrodynamics, which is widely used to determine the characteristics of gas-liquid-solid phase in fluidized beds and bubble columns [20, 21]. The gas holdup is the volume fraction of gas in the reactor, indicating the size of the contact area of gas with liquid. The gas holdup significantly affects the efficiencies of mass transfer, heat transfer, and microrate of reaction. The gas holdup distribution is mainly related to the superficial gas velocity, solids’ concentration, pressure distribution, density of each phase, and the physical properties of the liquid [22].

This is particularly important for the reactor proposed in this study, where the gas holdup characteristics are not yet clear, especially at a quantitative level. Various measuring techniques, such as pressure difference techniques [23, 24], high-speed camera [25], double optical probe [26, 27], particle image velocimetry (PIV) technology [28, 29], and tomography techniques [30, 31], have been proposed to study the hydrodynamic behaviors of gas-liquid flows. Numerical modeling and simulations were proposed to predict the distribution of a radial gas holdup, global gas holdup, and the two-phase flow field [3234]. These techniques provide useful information to better understand the gas-liquid flows under different conditions. However, these techniques may have certain limitations. For example, the measurement of gas holdup in the transient flow, which is related to fluid velocities, dimensions of the equipment, and the properties of the fluid, is generally difficult and, in many cases, involves high costs. However, such detailed information is highly desirable for the present new reactor.

In this paper, based on the water model experiments for the Venturi carbonation reactor, the gas holdup and the relationship between the gas holdup and various parameters in the carbonation process are studied, and an empirical correlation is derived for gas holdup as a function of superficial gas velocities, superficial liquid velocities, and liquid-to-solid and height-to-diameter ratio. The applicability of the correlation is validated, following which, the effects of various parameters on the gas holdup of the Venturi carbonation reactor can be quantified that benefits the performance of the novel reactor.

2. Experimental

2.1. Principles

For the simulation of carbonation reactor, the geometric and dynamic similarities are necessary. A water model with geometric similarity, in which the main geometrical sizes are proportional to the prototype on a scale of 1 : 6 is established. Meanwhile, the dynamic similarity follows the rule that the modified Froud numbers (Fr) for the model and prototype are equal. The modified Froud number is defined usingwhere u is the characteristic velocity (m/s), H is the depth of pool (m), ρl and are the densities of liquid and gas (kg/m3), respectively, and is the acceleration due to gravity (m/s2).

The characteristic velocity can be expressed usingwhere Q is the flow rate of fluid (m3/h), d is the throat diameter of the jet reactor (m), and π is the circular constant. Combining equations (2) and (3), the modified Froud number is obtained, which is given by

From , the quantity of fluid injected into the model can be calculated:where the subscripts p, m, l, and stand for prototype, model, liquid, and gas, respectively.

2.2. Apparatus

The water model used in the experiment refers to the design for industrial implementation of red mud process, which includes a main reactor, a fluid circulation system, and an observation system. The newly designed Venturi jet flow device used for the calcified slag carbonation is a prototype, in which the Venturi tube is used as a liquid entry point at the bottom. This part helps in generating a negative pressure relative to the ambient atmospheric pressure, which sucks gas into the reactor through gas inlet. The gas mainly escapes from the outlet at the top, whereas liquid and bubbles are removed from the right side of the outlet (Figure 1). The geometrical parameters for the reactor are listed in Table 1.

In the water model system, materials made up of organic glass were used to ensure a steady movement of the fluid and observe the mixing of fluid, circulation stream, and bubble behavior (Figure 2). The size and the distribution of gas bubbles and their movements were directly observable.

2.3. Methodology

The pressure-difference method was used to measure the gas holdup in the reactor. According to Bernoulli principle, the pressure difference between the two cross sections of fluidized bed in the reactor was measured. For the gas-liquid system, the density of the mixture of gas and liquid is given by

The pressure difference between the two cross sections is given by

The gas density is very small as compared to the liquid density ρl and, therefore, can be ignored. From equations (6) and (7), the gas holdup can be calculated usingwhere is the gas holdup, is the acceleration due to gravity (m/s2), ρmix is the density of gas-liquid mixture in the reactor (kg/m3), ρl is the density of liquid in the reactor (kg/m3), h is the altitude difference between the two measuring sections of a fluidized bed in a reactor (m), and P is the pressure difference between the two cross sections (kPa).

2.4. Dimensional Analysis

Dimensional analysis is a method for reducing complex physical problems to their simplest form by studying the dimensions of variables that describe the physical phenomenon. The principal use of dimensional analysis is to make the physical equation dimensionally homogeneous. The method is of great generality and mathematical simplicity. In engineering, empirical correlation among various parameters can be derived from experimental data using the method of dimensional analysis, which also enables weighing the magnitude of each operating parameter.

3. Dimensional Analysis of Gas Holdup

Based upon the water model experiments, the gas holdup () in the reactor is mainly related to superficial gas velocity (), superficial liquid velocity (Ul), liquid-to-solid ratio (Rls), and height-to-diameter ratio (Rhd) [23]. Based on the dimensional analysis, the influence of various factors on the gas holdup in the reactor is studied. According to Buckingham’s theorem (π - principle), a general functional form of the independent variables was obtained and is given in the following equations [3537]:orwhere is the gas holdup, and Ul are the superficial velocities of gas and liquid (m/s), respectively, Rls is the liquid-to-solid ratio, Rhd is the height-to-diameter ratio, and d1, d2, and d5 are the diameter of gas inlet of the jet reactor (m), diameter of the throat of the jet reactor (m), and diameter of the reactor (m), respectively. Furthermore, H is the height of liquid level in the reactor (m), and μl are the viscosities of gas and liquid (Pa·s), respectively, and and ρl are the densities of gas and liquid (kg/m3), respectively. The dimensions of the variables are listed in Table 2.

Based on π-principle, the total number of variables (n) is 13 and the number of independent variables (k) is 5, so that n − k = 8, which is the number of dimensionless quantities that can be set up. Furthermore, ,,ρl, μl, and d5 were selected as the independent variables, whereas , Rls, and Rhd were selected as the dimensionless values. Then, according to the π-principle, the dimensionless forms of the remaining n − k independent variables were defined as follows [3537]:

This way, a dimensionless functional form is obtained as follows:

An explicit functional form is given by

In addition, d1, d2, and H are quantitative in nature. Therefore, equation (20) can be expressed as

4. Derivation of the Empirical Correlation

Generally, the criterion equation can be expressed in the form of power function with the independent variable being under certain conditions. Therefore, equation (21) can be rewritten aswhere a0, a1, a2, a3, and a4 are the fitting coefficients. After taking the logarithm on both sides of equations (22), equation (23) is obtained:

Based on the linear relationship of equation (23), the fitting coefficients a1, a2, a3, and a4 are obtained by plotting the experimental data in Origin 8.6 (See Figure S1 in the Supplementary Material), whereas a0 is optimized using the Quasi-Newton method:

The values of the coefficients are the mean values of the slopes in their respective graphs. The fitting equation can be expressed using

The results obtained using equation (25) are consistent with the experimental results. The mean statistical variance is 0.97 (See Table S2 in the Supplementary Material), while the mean relative error is within 8.7% (See Table S3 in the Supplementary Material). Furthermore, equation (25) is valid within a range of the following parameters: of 7.08–21.23 m·s−1; Ul of 0.89–1.77 m·s−1; Rls of 10–20; Rhd of 3–5. Figure 3 shows that the calculated results obtained from the fitting equation agree well with the experiment data. Equation (25) is significant for designing a highly efficient Venturi reactor involving high temperature, high pressure, and three phases of gas, liquid, and solid.

5. Discussion

5.1. Compatibility of the Similarity Criteria

The establishment of the criteria model is very simple. However, the absence of the physical analysis leads to incorrect results [38]. In this study, thirteen physical quantities involved with the gas holdup were analyzed (equation (9)). These quantities are the determinants, which significantly affect the size, amount, and residence time of gas bubbles in the fluid. Some physical quantities involved with the gas holdup include temperature and surface tension. These parameters are relatively stable in fluid flow and do not have a considerable effect on the gas holdup and, therefore, are not included in the model. In fact, the densities of gas and liquid ( and ρl) and viscosities of gas and liquid ( and μl) are stable. Their dimensions can be used to make dimensionless complexes according to π-theorem which requires that both sides of the physical equation must be dimensionally homogeneous. Eight dimensionless complexes were made (see equations (11)–(18)), wherein Π3 and Π4 represent the Reynolds numbers, which often appear in the analysis of heat and mass transfer in fluid flow:

The Reynolds numbers are the main similarity criteria characterizing the similarities of inertia force and viscosity force in fluid flow. Π2, Π5, and Π6 characterize the geometric similarities of apparatus, and Π0, Π1, and Π7 assure the similarities of magnitude of gas, solid, and liquid. The criterion model for estimating the gas holdup of the jet reactor is constructed in the case of a serious physical analysis and understanding of its physical background. Incompatible similarity criteria are not found in the dimensionless complexes. Therefore, it is available in the process of multiphase flow.

The data used in the dimensional analysis was obtained from water model experiments, which were based on the geometric and dynamic similarities. Froude number was used as the dynamic similarity criteria to construct the water model. Using the throat diameter of the jet reactor as a length scale, the Froude number characterizes the similarities of inertia force and gravity in fluid, which assure the dynamic similarity of the water model with the prototype (see Section 2.1).

Figure 4 shows a comparison of the results of the proposed correlation with those of the water experiments. It can be seen that the calculated results are in good agreement with the experimental results. Therefore, the proposed correlation can be used for predicting the gas holdup in the novel reactor.

5.2. Effect of Various Parameters on Gas Holdup

The water model experiments showed that the gas holdup in the Venturi reactor is mainly affected by the superficial gas velocity, superficial liquid velocity, liquid-to-solid ratio (Rls), and height-to-diameter ratio (Rhd). Among them, the superficial gas velocity showed an obvious impact on gas holdup. As shown in Figure 3(a), when Rls = 20, Rhd = 4, and Ul = 1.33 m·s−1, gas holdup increased with the increase in superficial gas velocity. Notably, the increase in gas velocity increased the number of gas bubbles. With the increase in the superficial gas velocity, the strength of stirring also intensified. Larger gas bubbles were disintegrated into smaller bubbles, whereas the rising velocity of small gas bubbles became slower, due to which, the gas holdup increased.

The superficial liquid velocities showed a weaker effect on gas holdup. With the increase in the superficial liquid velocity, the gas holdup decreased (Figure 3(b)). Fan et al. [39] reported that the bubble rise velocity () can be expressed in terms of a bubble slip velocity (Us) and the linear liquid velocity (Ul) using

Equation (27) shows that the gas holdup has an inverse relationship with the liquid velocity. In air-water systems, the slip velocity is the major contributor to bubble rise velocity (/) [39]. However, the bubble rise velocity increases with the increase in liquid velocity, which decreases the gas holdup.

When solid particles are added in the solution, the properties of liquid, such as viscosity, change. The existence of solid particles and the increase in viscosity prevented the disintegration of large bubbles into small bubbles and increased the coalescence of smaller bubbles into larger ones. Furthermore, the residence time of bubbles in the solution decreased. The increase in liquid-to-solid ratio means that the solid particles decreased, due to which the gas holdup increased (Figure 3(c)).

The ratio of the height to diameter of the reactor (Rhd) is an important parameter for the equipment. The experimental conditions were set to be as follows: Rhd = 3, 4, and 5, Ul = 1.33 m·s−1, and Rls = 20. As shown in Figure 3(d), for low superficial gas velocity and high height-to-diameter ratio, the results show a significant difference between the calculated and experimental values. However, with the increase in the value of Rhd, the gas holdup increased.

The increase in temperature led to an expansion of the volume of bubbles, due to which the gas density decreased. According to eqaution (6), the decreased gas density is negligible as compared to the liquid density. This change was taken into consideration in the dynamic similarity, where the flow rate of the gas and gas density were included in the model (equation (4)). Therefore, the gas density and liquid density in equation (25) should be considered when designing industrial unit operations. Similarly, the viscosity of red mud slurry is affected by temperature and is also considered in the model.

The surface tension of bubbles may affect the gas holdup in fluid. Liu et al. [40] discussed the applicability of water experimental results to the hot metal ladle in plant and reported that the differences in the surface tension in molten iron were 3.5 times that of water. Mach [41] reported that the surface tension had a significant effect on gas holdup and hydrodynamics in a three-phase fluidized bed, while viscosity had only a slight impact. Lower liquid surface tension facilitates the formation of smaller bubbles and promotes a noncoalescing tendency in the rising gas bubbles, thus reducing the growth of bubbles. Surfactants can increase drag on gas bubbles due to the formation of a surface tension gradient on the bubble surface [39]. A higher drag results in lower bubble rise velocities, which in turn results in a higher gas holdup.

6. Conclusions

Based upon water model experiments, the effects of superficial gas velocities, superficial liquid velocities, and liquid-to-solid and height-to-diameter ratios on gas holdup of the Venturi carbonation reactor were studied. Based upon dimensional analysis and Buckingham’s theorem, an empirical correlation was proposed.

Superficial gas velocity is primarily the major influencing factor for gas holdup of the Venturi carbonation reactor. With the increase in superficial gas velocity, the bubbles in the solution and gas holdup increased. The superficial liquid velocity and the values of Rls and Rhd significantly impacted the gas holdup in the reactor. There is a negative correlation between the superficial liquid velocity and gas holdup of the reactor.

The results obtained using the empirical correlation are consistent with the experimental results and produce a statistical variance of 0.97 and mean relative error of less than 8.7%. The results are significant for designing a highly efficient Venturi reactor having high temperature, high pressure, and three phases of gas, liquid, and solid.

Notations

a0, a1, a2, a3, and a4:Undetermined constants
d:Characteristic diameter of pipe, m
d1:Gas inlet diameter of the reactor, m
d2:Throat diameter of the reactor, m
d5:Diameter of the reactor, m
Fr′:Modified Froud number
:Acceleration due to gravity, m/s2
h:Altitude difference between the two measuring sections of a fluidized bed, m
H:Height of liquid level in the reactor, m
Rls:Liquid-to-solid ratio
Rhd:Height-to-diameter ratio
Q:Flow rate of fluid, cm3/min
u:Characteristic velocity, m/s
U:Superficial velocity, m/s
P:Pressure difference between the two measuring sections, Pa
ρ:Density of fluid, kg/m3
ρmix:Density of the mixture of gas and liquid, kg/m3
μ:Viscosity of the fluid, Pa·s
ε:Gas holdup
π:Circular constant
Subscripts
p, m, l, s, and :The prototype, model, liquid, solid, and gas, respectively.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors are grateful to the National Natural Science Foundation of China and National 863 Plan of China for the financial support (Grant no. U1202274 and Grant no. 51204040).

Supplementary Materials

Figure S1: fitting coefficients a1, a2, a3, and a4 as obtained according to the slopes: (a) a1; (b) a2; (c, d) a3; (e) a4. Table S2: mean statistical variance of the results presented in Figure S1: the first and second columns are obtained from Figure S1 (a); first and third columns are obtained from Figure S1 (b); fourth and fifth columns are obtained from Figure S1 (c); sixth and seventh columns are obtained from Figure S1 (d); eighth column is obtained from Figure S1 (e). Table S3: relative error between the calculated and experimental results: The relative error: δ = (Xi − Xc)/Xt × 100, where Xi is the experimental value and Xc is the calculated value. The mean relative error: , where σ = 8.7%. (Supplementary Materials)