CFD Study of Gas Holdup and Frictional Pressure Drop of Vertical Riser Inside IC Reactor

: The internal circulation system in Internal Circulation (IC) reactor plays an important role in increasing volumetric loading rate and promoting the mixing between sludge and wastewater. In order to design the internal circulation system, the ﬂow behaviors of gas-liquid inside vertical riser should be studied in detail. In the present study, the Multiple Flow Regimes model is adopted to capture the phase interface for di ﬀ erent ﬂow conditions. The ﬂow patterns, internal circulation ﬂow rate, gas holdup, and frictional pressure drop of vertical riser are investigated. The results show that the bubble ﬂow inside a vertical riser is in a stable ﬂow condition. There exists a maximum value for internal circulation ﬂow rate with the increasing superﬁcial gas velocity. The parameters of Martinelli models for gas holdup and frictional pressure drop are improved based on Computational Fluid Dynamics (CFD) results. The deviations between the calculated gas holdup and frictional pressure drop by improved model and experimental value are reduced to 14% and 13.2%, respectively. The improved gas holdup and frictional pressure drop model can be used for the optimal design of internal circulation system. riser.


Introduction
The internal circulation system in Internal Circulation (IC) reactor can increase the volumetric loading rate and biogas production obviously [1]. It provides the high superficial liquid velocity to make the granular sludge fluidized and promotes the mixing between sludge and sewage. The vertical riser supplies the driving force for the internal circulation, which is closely related to its structure and flow patterns. The detailed study of flow behaviors inside vertical riser could provide the guidelines for the design of IC reactor.
Since 1990s, the design of IC reactor relies on the empirical correlations for gas holdup and superficial liquid velocity. Pereboom [2] proposed a mathematical equation for the calculation of the superficial liquid velocity of vertical riser. Based on this equation, the internal circulation flow rate can be obtained. Hu [3] improved Pereboom's equation considering the actual operational condition and a deepened mechanical analysis. However, these equations were derived based on pilot tests or mathematical inference. In addition, some parameters are set as constant. Therefore, the empirical correlations may not be suitable for all conditions and result in large deviations for actual operation.
The calculation of frictional pressure drop is important for the force balance analysis in the design of IC anaerobic reactor. Many researchers have focused on the calculation of frictional pressure drop, and different calculation formulas were developed based on experimental results [4][5][6]. Up to now, where is the mean density of liquid and gas, is the length of the vertical riser, α is the gas holdup, is the density of gas, is the density of liquid, is the mass flux of liquid, is the mass flux of gas. Martinelli models [30] are widely used to predict the gas holdup α of gas-liquid two-phase flow. = . · . · . (8) where is the Martinelli parameter, is the mass flow rate of liquid, and is the mass flow rate of gas (kg/s).
Martinelli gave the relationship between the parameter defined by the above formulas and the gas holdup α in the gas-liquid two-phase flow as a graph, while Chisholm [31] expressed it as the following formula: The frictional pressure drop of gas-liquid two-phase flow in the pipeline can be expressed as follows: where ϕ is the two-phase friction multiplier, (∆ ) is the frictional pressure drop when there is only liquid in the pipe. The two-phase gravitational and accelerational pressure drop are determined as follows [29]: where ρ m is the mean density of liquid and gas, L is the length of the vertical riser, α is the gas holdup, ρ g is the density of gas, ρ l is the density of liquid, G l is the mass flux of liquid, G g is the mass flux of gas. Martinelli models [30] are widely used to predict the gas holdup α of gas-liquid two-phase flow.
where X tt is the Martinelli parameter, W l is the mass flow rate of liquid, and W g is the mass flow rate of gas (kg/s). Martinelli gave the relationship between the parameter X tt defined by the above formulas and the gas holdup α in the gas-liquid two-phase flow as a graph, while Chisholm [31] expressed it as the following formula: The frictional pressure drop of gas-liquid two-phase flow in the pipeline can be expressed as follows: where φ 2 is the two-phase friction multiplier, (∆P F ) l the frictional pressure drop when there is only liquid in the pipe.
φ 2 can be calculated as follows: The internal circulation flow rate can be calculated as flow: According to Equation (11), if the gas flow rate in the riser is known, the internal circulation flow rate can be calculated by iteration. However, for a lack of better prediction method for gas holdup and frictional pressure drop, the internal circulation flow rate and the design of IC anaerobic reactor are still difficult to achieve. In the present study, the gas holdup and frictional pressure drop can be obtained by CFD method.

Multiple Flow Regimes (MFR) Model
The gas-liquid two-phase flow in the vertical riser was simulated by MFR model. This model could describe the dispersed and segregated phase flows within a single framework. The topological structure and physical characteristics of phase interface could be redefined. The calculation process and equations in this section for two-phase flow can be summarized as below [23]:

1.
For the interaction of multiphase flow, a primary phase and a secondary phase are defined. For example, water is designated as the first phase and air as the second phase.

2.
According to the definition of interphase action, three topological structures are defined: the main phase flow pattern, air as the dispersed phase (bubble); the interface flow pattern, water and air are not the dispersed phase; the secondary phase flow pattern, water as the dispersed phase (droplet).

3.
The mesh with large interface is detected, the general topological structure in each mesh is determined, and the weight functions of each topological structure are calculated, which are M f r , M ir , and M sr , respectively. 4.
Calculation of the interphase forces of each sub-topology in all meshes: F f r , F ir , F sr . The subscripts fr, ir, and sr represent first regime, interface regime, and second regime, respectively. The relationships between weight function of three regimes and volume fraction of second phase can be found in Figure 2.

5.
The total interaction force of each sub-topology is calculated by weight function and interaction force: In the MFR model, the large scale interface detection model (LIM) [28] is used to detect the grid in a large interface. A detailed description can be found in STARCCM+ User Guide 12.02 [23].
For a pair of primary and secondary phases, the flow topology is defined based on the volume fraction of the secondary phase (ε s )as: 0.0 < ε s < ε f r : First Regime ε f r < ε s < ε sr : Interface Regime ε sr < ε s < 1.0 : Second Regime When the first regime terminus (ε f r ) is the value of ε s for which the first regime transits to the interface regime. The default value is taken as 0.3 [23]. When the second regime onset (ε sr ) is the value of ε s for which the interface regime transits to the second regime. The default value is taken as 0.3. For a pair of primary and secondary phases, the flow topology is defined based on the volume fraction of the secondary phase ( )as: 0.0 < < : First Regime < < : Interface Regime < < 1.0: Second Regime When the first regime terminus ( ) is the value of for which the first regime transits to the interface regime. The default value is taken as 0.3 [23]. When the second regime onset ( ) is the value of for which the interface regime transits to the second regime. The default value is taken as 0.3.

Geometry Model and Mesh Generation
The geometric model of internal circulation system is shown in Figure 3, which includes a vertical riser, gas collector, and the second reaction chamber. The bottom surface is set as the velocity inlet of gas. The inlet diameter is 1.75 m and the diameter of vertical riser is 0.2 m. The inlet gas velocity is between 0.04486 ~ 0.314 m/s, i.e., 0.108 ~ 0.755 kg/s. The top outlet is the pressure outlet (101325 Pa), and the liquid level of second reaction chamber is the pressure inlet (101325 Pa). All walls are set as no-slip condition (fluid velocity is zero). Under different inlet velocities, the apparent gas velocity(vg) in the vertical riser is between 3.115 ~ 21.81 m/s, and the immersion ratio is between 0.5 ~ 0.9.

Geometry Model and Mesh Generation
The geometric model of internal circulation system is shown in Figure 3, which includes a vertical riser, gas collector, and the second reaction chamber. The bottom surface is set as the velocity inlet of gas. The inlet diameter is 1.75 m and the diameter of vertical riser is 0.2 m. The inlet gas velocity is between 0.04486~0.314 m/s, i.e., 0.108~0.755 kg/s. The top outlet is the pressure outlet (101,325 Pa), and the liquid level of second reaction chamber is the pressure inlet (101,325 Pa). All walls are set as no-slip condition (fluid velocity is zero). Under different inlet velocities, the apparent gas velocity(v g ) in the vertical riser is between 3.115~21.81 m/s, and the immersion ratio is between 0.5~0.9. For a pair of primary and secondary phases, the flow topology is defined based on the volume fraction of the secondary phase ( )as: 0.0 < < : First Regime < < : Interface Regime < < 1.0: Second Regime When the first regime terminus ( ) is the value of for which the first regime transits to the interface regime. The default value is taken as 0.3 [23]. When the second regime onset ( ) is the value of for which the interface regime transits to the second regime. The default value is taken as 0.3.

Geometry Model and Mesh Generation
The geometric model of internal circulation system is shown in Figure 3, which includes a vertical riser, gas collector, and the second reaction chamber. The bottom surface is set as the velocity inlet of gas. The inlet diameter is 1.75 m and the diameter of vertical riser is 0.2 m. The inlet gas velocity is between 0.04486 ~ 0.314 m/s, i.e., 0.108 ~ 0.755 kg/s. The top outlet is the pressure outlet (101325 Pa), and the liquid level of second reaction chamber is the pressure inlet (101325 Pa). All walls are set as no-slip condition (fluid velocity is zero). Under different inlet velocities, the apparent gas velocity(vg) in the vertical riser is between 3.115 ~ 21.81 m/s, and the immersion ratio is between 0.5 ~ 0.9.  The geometric model of internal circulation system was discretized by using polyhedral mesh. The mesh size was determined through the mesh independence test. As shown in Figure 4, the liquid flow rate varied with the base mesh size. It can be seen that when the base mesh size is less than 0.24 m, the liquid flow rate has no obvious change. Therefore, the base mesh size of 0.24 m was adopted for all The geometric model of internal circulation system was discretized by using polyhedral mesh. The mesh size was determined through the mesh independence test. As shown in Figure 4, the liquid flow rate varied with the base mesh size. It can be seen that when the base mesh size is less than 0.24 m, the liquid flow rate has no obvious change. Therefore, the base mesh size of 0.24 m was adopted for all cases in this work. The detailed parameters used for the surface and volume mesh generation can be found in Eppinger's work [32]. The mesh topology of computational domain was shown in Figure 5.   Figure 6 shows the distributions of gas holdup of the riser at different superficial gas velocities. It can be seen that as the superficial gas velocity increases, the flow patterns gradually change from a bubble flow, via a slug flow to a churn flow. As shown in Figure 7, under the same superficial gas velocity of riser (vg), the driving force for the liquid rise increases as the immersion ratio increase. The flow pattern changes from a slug flow to a bubble flow. The flow patterns are determined by the relative magnitudes of superficial gas velocity and liquid velocity. For higher superficial gas velocity, the gas bubbles tend to be aggregated and larger local gas holdups can be found. The geometric model of internal circulation system was discretized by using polyhedral mesh. The mesh size was determined through the mesh independence test. As shown in Figure 4, the liquid flow rate varied with the base mesh size. It can be seen that when the base mesh size is less than 0.24 m, the liquid flow rate has no obvious change. Therefore, the base mesh size of 0.24 m was adopted for all cases in this work. The detailed parameters used for the surface and volume mesh generation can be found in Eppinger's work [32]. The mesh topology of computational domain was shown in Figure 5.   Figure 6 shows the distributions of gas holdup of the riser at different superficial gas velocities. It can be seen that as the superficial gas velocity increases, the flow patterns gradually change from a bubble flow, via a slug flow to a churn flow. As shown in Figure 7, under the same superficial gas velocity of riser (vg), the driving force for the liquid rise increases as the immersion ratio increase. The flow pattern changes from a slug flow to a bubble flow. The flow patterns are determined by the relative magnitudes of superficial gas velocity and liquid velocity. For higher superficial gas velocity, the gas bubbles tend to be aggregated and larger local gas holdups can be found.  Figure 6 shows the distributions of gas holdup of the riser at different superficial gas velocities. It can be seen that as the superficial gas velocity increases, the flow patterns gradually change from a bubble flow, via a slug flow to a churn flow. As shown in Figure 7, under the same superficial gas velocity of riser (v g ), the driving force for the liquid rise increases as the immersion ratio increase. The flow pattern changes from a slug flow to a bubble flow. The flow patterns are determined by the relative magnitudes of superficial gas velocity and liquid velocity. For higher superficial gas velocity, the gas bubbles tend to be aggregated and larger local gas holdups can be found.   For the bubble flow, the interaction between bubbles is smaller and its distribution is more uniform. Therefore, the liquid flow rate is more stable under the bubble flow. As can be seen in Figure 8B, when t = 15.23 s, a large gas slug appears at the exit. The outlet of the vertical riser is occupied by gas, resulting in a decrease of liquid flow rate. When the gas   For the bubble flow, the interaction between bubbles is smaller and its distribution is more uniform. Therefore, the liquid flow rate is more stable under the bubble flow. As can be seen in Figure 8B, when t = 15.23 s, a large gas slug appears at the exit. The outlet of the vertical riser is occupied by gas, resulting in a decrease of liquid flow rate. When the gas  the bubble flow. As can be seen in Figure 8B, when t = 15.23 s, a large gas slug appears at the exit. The outlet of the vertical riser is occupied by gas, resulting in a decrease of liquid flow rate. When the gas slug passes through the riser outlet, the liquid flow rate changed with time into a relatively stable stage, in which the riser outlet mainly presents a bubble flow. When t = 16.4 s, the exit of the riser presents the slug flow pattern again. As shown in Figure 8C, due to the blocky gases, the liquid flow rate is extremely unstable. In conclusion, the flow pattern in the riser should be controlled within the bubble flow as far as possible for the design of internal circulation system.  Figure 9 shows the liquid flow rate under different superficial gas velocities and immersion ratios. It can be seen that with the increase of immersion ratio, the internal circulation flow rate increases gradually. The larger immersion ratios supply a larger driving force of internal circulation. With the increase of the superficial gas velocity, the internal circulation flow rate does not always increase. When the superficial gas velocity is less than 4.36 m/s, the increase of internal circulation flow rate is obvious. For the bubble flow (smaller superficial gas velocity), the interaction between gas and liquid is stronger than that of slug flow. However, when the superficial gas velocity is greater than 7.27 m/s, the internal circulation flow rate decreases with the increase of the superficial gas velocity. Under this circumstance, the interaction between gas and liquid is weak, and the gas slug  Figure 9 shows the liquid flow rate under different superficial gas velocities and immersion ratios. It can be seen that with the increase of immersion ratio, the internal circulation flow rate increases gradually. The larger immersion ratios supply a larger driving force of internal circulation. With the increase of the superficial gas velocity, the internal circulation flow rate does not always increase. When the superficial gas velocity is less than 4.36 m/s, the increase of internal circulation flow rate is obvious. For the bubble flow (smaller superficial gas velocity), the interaction between gas and liquid is stronger than that of slug flow. However, when the superficial gas velocity is greater than 7.27 m/s, the internal circulation flow rate decreases with the increase of the superficial gas velocity. Under this circumstance, the interaction between gas and liquid is weak, and the gas slug or churn is easy to get out of the riser. The maximum flow rate of the internal circulation appears when the immersion ratio is 0.9 and the superficial gas velocity is between 3.16 m/s~4.36 m/s. The high superficial gas velocity will lead to the center position of the riser being occupied by gas, resulting in a decrease of liquid flow rate. When the immersion ratio is 0.8 and the superficial gas velocity is 7.27 m/s, the liquid flow rate is 110.71 kg/s. When the immersion ratio is 0.9 and the superficial gas velocity is 0.778 m/s, the liquid flow rate is 108.35 kg/s. There is little difference in liquid flow rate between these two conditions, but the difference in superficial gas velocity is about 9.5 times. This result shows that it is not an effective solution to increase the internal circulation flow rate by increasing the superficial gas velocity (reducing the diameter of the riser). The measures to be taken is increasing the height from the bottom of the second reaction chamber to the outlet of the wastewater as much as possible, and then adjust the superficial gas velocity between 3.16 m/s and 4.36 m/s to achieve a higher internal circulation flow rate.

Flow Patterns and Internal Circulation Flow Rate
Processes 2019, 7, x FOR PEER REVIEW 9 of 14 or churn is easy to get out of the riser. The maximum flow rate of the internal circulation appears when the immersion ratio is 0.9 and the superficial gas velocity is between 3.16 m/s ~ 4.36 m/s. The high superficial gas velocity will lead to the center position of the riser being occupied by gas, resulting in a decrease of liquid flow rate. When the immersion ratio is 0.8 and the superficial gas velocity is 7.27 m/s, the liquid flow rate is 110.71 kg/s. When the immersion ratio is 0.9 and the superficial gas velocity is 0.778 m/s, the liquid flow rate is 108.35 kg/s. There is little difference in liquid flow rate between these two conditions, but the difference in superficial gas velocity is about 9.5 times. This result shows that it is not an effective solution to increase the internal circulation flow rate by increasing the superficial gas velocity (reducing the diameter of the riser). The measures to be taken is increasing the height from the bottom of the second reaction chamber to the outlet of the wastewater as much as possible, and then adjust the superficial gas velocity between 3.16 m/s and 4.36 m/s to achieve a higher internal circulation flow rate. The flow patterns and stripping ratios (Vl/Vg) under different immersion ratios and superficial gas velocities are shown in Figures 10 and 11. Vl and Vg represent the volume flow rate of liquid and gas phase inside a vertical riser, which was obtained by CFD simulations. It can be seen that with the increase of immersion ratio and the decrease of superficial gas velocity, the stripping ratio increases and the flow pattern tends to be bubble flow.
For the design of IC anaerobic reactor, it is necessary to keep the internal circulation flow rate stable. From Figure 8, it can be seen that the most favorable flow pattern is bubble flow. As shown in Figures 9 and 10, the internal circulation flow rate decreases with the increase of the superficial gas velocity when the flow pattern is churn flow. Therefore, the bubble flow pattern can be realized by adjusting the superficial gas velocity and immersion ratio according to Figure 9. Gas stripping ratios can be obtained from Figure 11.  The flow patterns and stripping ratios (V l /V g ) under different immersion ratios and superficial gas velocities are shown in Figures 10 and 11. V l and V g represent the volume flow rate of liquid and gas phase inside a vertical riser, which was obtained by CFD simulations. It can be seen that with the increase of immersion ratio and the decrease of superficial gas velocity, the stripping ratio increases and the flow pattern tends to be bubble flow.
Processes 2019, 7, x FOR PEER REVIEW 9 of 14 or churn is easy to get out of the riser. The maximum flow rate of the internal circulation appears when the immersion ratio is 0.9 and the superficial gas velocity is between 3.16 m/s ~ 4.36 m/s. The high superficial gas velocity will lead to the center position of the riser being occupied by gas, resulting in a decrease of liquid flow rate. When the immersion ratio is 0.8 and the superficial gas velocity is 7.27 m/s, the liquid flow rate is 110.71 kg/s. When the immersion ratio is 0.9 and the superficial gas velocity is 0.778 m/s, the liquid flow rate is 108.35 kg/s. There is little difference in liquid flow rate between these two conditions, but the difference in superficial gas velocity is about 9.5 times. This result shows that it is not an effective solution to increase the internal circulation flow rate by increasing the superficial gas velocity (reducing the diameter of the riser). The measures to be taken is increasing the height from the bottom of the second reaction chamber to the outlet of the wastewater as much as possible, and then adjust the superficial gas velocity between 3.16 m/s and 4.36 m/s to achieve a higher internal circulation flow rate. The flow patterns and stripping ratios (Vl/Vg) under different immersion ratios and superficial gas velocities are shown in Figures 10 and 11. Vl and Vg represent the volume flow rate of liquid and gas phase inside a vertical riser, which was obtained by CFD simulations. It can be seen that with the increase of immersion ratio and the decrease of superficial gas velocity, the stripping ratio increases and the flow pattern tends to be bubble flow.
For the design of IC anaerobic reactor, it is necessary to keep the internal circulation flow rate stable. From Figure 8, it can be seen that the most favorable flow pattern is bubble flow. As shown in Figures 9 and 10, the internal circulation flow rate decreases with the increase of the superficial gas velocity when the flow pattern is churn flow. Therefore, the bubble flow pattern can be realized by adjusting the superficial gas velocity and immersion ratio according to Figure 9. Gas stripping ratios can be obtained from Figure 11.

Gas Holdup of Vertical Riser
As shown in Figure 12, with the increase of superficial liquid velocity (vl), the gas holdup decreases continuously. There exists approximate linear relationship between the gas holdup and superficial liquid velocity.
The constants (a and b) are set as 1 and 21, respectively [30]. Based on the CFD simulation results for gas holdup, the Martinelli parameter ( ) is fitted and the improved form is given below: Figure 11. Gas stripping ratios V l /V g at different immersion ratios.
For the design of IC anaerobic reactor, it is necessary to keep the internal circulation flow rate stable. From Figure 8, it can be seen that the most favorable flow pattern is bubble flow. As shown in Figures 9 and 10, the internal circulation flow rate decreases with the increase of the superficial gas velocity when the flow pattern is churn flow. Therefore, the bubble flow pattern can be realized by adjusting the superficial gas velocity and immersion ratio according to Figure 9. Gas stripping ratios can be obtained from Figure 11.

Gas Holdup of Vertical Riser
As shown in Figure 12, with the increase of superficial liquid velocity (v l ), the gas holdup decreases continuously. There exists approximate linear relationship between the gas holdup and superficial liquid velocity.
Processes 2019, 7, x FOR PEER REVIEW 10 of 14 Figure 11. Gas stripping ratios / at different immersion ratios.

Gas Holdup of Vertical Riser
As shown in Figure 12, with the increase of superficial liquid velocity (vl), the gas holdup decreases continuously. There exists approximate linear relationship between the gas holdup and superficial liquid velocity.
The constants (a and b) are set as 1 and 21, respectively [30]. Based on the CFD simulation results for gas holdup, the Martinelli parameter ( ) is fitted and the improved form is given below:  Figure 13 shows the relationship between gas holdup obtained by CFD method and Martinelli parameter (X tt ). The Martinelli parameter (X tt ) is given as [24]: Gas hold up X tt Figure 13. Gas holdup as a function of Martinelli parameter ( ).
The gas holdup obtained by Hu's work [2,33] was compared with that by Martinelli and improved gas holdup model (Equation (14)). The calculated deviation for the gas holdup by Martinelli model and the experimental result is 50.2%, while the deviation is 13.7% when calculated by Equation (14). The improved gas holdup model can effectively improve the accuracy of gas holdup prediction in gas-liquid two-phase flow of vertical riser.

Frictional Pressure Drop of Vertical Riser
The most important contribution of this paper is the improvement of the calculation of the frictional pressure drop. The relationship between gas hold up and the Martinelli parameter ( ) can be given as: According to the simulation results in Figure 14, the two-phase friction multiplier (Equation (8)) can be calculated. The relationship between (Equation (8)) and Martinelli parameter ( ) was plotted as shown in Figure 15. The parameters in Equation (18) of a and b was fitted and the improved two-phase friction multiplier is given as below: Therefore, the improved frictional pressure drop model consists of Equations (7)-(9) and Equation (19).  The constants (a and b) are set as 1 and 21, respectively [30]. Based on the CFD simulation results for gas holdup, the Martinelli parameter (X tt ) is fitted and the improved form is given below: The gas holdup obtained by Hu's work [2,33] was compared with that by Martinelli and improved gas holdup model (Equation (14)). The calculated deviation for the gas holdup by Martinelli model and the experimental result is 50.2%, while the deviation is 13.7% when calculated by Equation (14). The improved gas holdup model can effectively improve the accuracy of gas holdup prediction in gas-liquid two-phase flow of vertical riser.

Frictional Pressure Drop of Vertical Riser
The most important contribution of this paper is the improvement of the calculation of the frictional pressure drop. The relationship between gas hold up and the Martinelli parameter (X tt ) can be given as: According to the simulation results in Figure 14, the two-phase friction multiplier φ 2 (Equation (8)) can be calculated. The relationship between φ 2 (Equation (8)) and Martinelli parameter (X tt ) was plotted as shown in Figure 15. The parameters in Equation (15) of a and b was fitted and the improved two-phase friction multiplier φ 2 is given as below: Therefore, the improved frictional pressure drop model consists of Equations (7)-(9) and Equation (16).
The frictional pressure drops calculated by Martinelli model and improved model in this study were compared with the experimental results of Zhou et al. [34]. The deviation between the frictional pressure drop obtained by Martinelli model and the experimental value is 71.7%, while the deviation between the improved model and the experimental value is 13.2%. The improved frictional pressure drop model can significantly improve the accuracy of the prediction of pressure drop in vertical riser.
two-phase friction multiplier is given as below: Therefore, the improved frictional pressure drop model consists of Equations (7)-(9) and Equation (19).  The frictional pressure drops calculated by Martinelli model and improved model in this study were compared with the experimental results of Zhou et al. [34]. The deviation between the frictional pressure drop obtained by Martinelli model and the experimental value is 71.7%, while the deviation between the improved model and the experimental value is 13.2%. The improved frictional pressure drop model can significantly improve the accuracy of the prediction of pressure drop in vertical riser.

Conclusions
In this study, the hydrodynamic behaviors of riser in IC anaerobic reactor were studied by multiple flow regimes model. The results show that the bubble flow has a smaller fluctuation range and shorter fluctuation period compared with the slug flow and the bulk flow. When the immersion ratio is constant, there is a peak value of liquid flow rate with the increase of superficial gas velocity. According to the analysis of the flow pattern and the liquid flow rate, bubble flow is considered as the best flow pattern in riser. Based on the simulation data, the parameters of Martinelli gas holdup model and the two-phase friction multiplier model have been corrected. Compared with the experimental data, the calculated deviations of the gas holdup by the corrected model are reduced from 44% and 50% to 13.7% and 14.4%, respectively. The deviation between the calculated friction pressure drop and the experimental data is reduced from 71.7% to 13.2%, which effectively improves the accuracy of prediction of the gas holdup and frictional pressure drop of vertical riser. The corrected model can be used to design the internal circulation system and IC reactor.

Conclusions
In this study, the hydrodynamic behaviors of riser in IC anaerobic reactor were studied by multiple flow regimes model. The results show that the bubble flow has a smaller fluctuation range and shorter fluctuation period compared with the slug flow and the bulk flow. When the immersion ratio is constant, there is a peak value of liquid flow rate with the increase of superficial gas velocity. According to the analysis of the flow pattern and the liquid flow rate, bubble flow is considered as the best flow pattern in riser. Based on the simulation data, the parameters of Martinelli gas holdup model and the two-phase friction multiplier model have been corrected. Compared with the experimental data, the calculated deviations of the gas holdup by the corrected model are reduced from 44% and 50% to 13.7% and 14.4%, respectively. The deviation between the calculated friction pressure drop and the experimental data is reduced from 71.7% to 13.2%, which effectively improves the accuracy of prediction of the gas holdup and frictional pressure drop of vertical riser. The corrected model can be used to design the internal circulation system and IC reactor.