Particle-laden bubble size and its distribution in microstructured

Jan 21, 2019 - Bubble size distribution (BSD) in three different experimental ... ratio, and interfacial area in the three-phase system by considering...
0 downloads 0 Views 3MB Size
Subscriber access provided by EKU Libraries

General Research

Particle-laden bubble size and its distribution in microstructured bubbling bed in presence and absence of a surface active agent Ritesh Prakash, Subrata Majumder, and Anugrah Singh Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.8b05625 • Publication Date (Web): 21 Jan 2019 Downloaded from http://pubs.acs.org on January 24, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Particle-laden bubble size and its distribution in microstructured bubbling bed in presence and absence of a surface active agent Ritesh Prakash*, Subrata Kumar Majumder* and Anugrah Singh Department of Chemical Engineering, Indian Institute of Technology Guwahati, Guwahati781039, Assam, India *Corresponding author Prof. S. K. Majumder/ Mr. Ritesh Prakash Email: [email protected] or [email protected] Abstract This work reports experimental measurements of the bubble size and its distribution in the threephase (air-water-coal, air-water-coal-surfactant) system. A thorough investigation of the effect of particle concentration, particle size, viscosity, and axial positions on bubble size and its distribution in the presence and absence of the surfactant is reported. Bubble size distribution (BSD) in three different experimental conditions is studied which follows the beta, Weibull, and log-logistic nature of the distribution function. Correlations are developed for estimation of parameters of each distribution function. A generalized correlation model is also developed for the estimation of the Sauter mean bubble diameter, aspect ratio, and interfacial area in the threephase system by considering operating, geometric variables, and physical properties of the system. The current work will be useful for process intensification of chemical and biochemical processes based on the interfacial phenomena. Keywords: bubble size; bubble distribution function; particle concentration; particle size, surfactant

1 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 59

1. Introduction Slurry columns are widely used three-phase contactors due to their simple construction, better transport efficiency, and low operating cost

1, 2.

The bubble size and its distribution in a non-

Newtonian systems are of great importance in many areas of chemical 3, mineral 4, and biological industries 5. Typical applications of the slurry columns are flotation, petroleum, wastewater treatment, chemical, and metallurgical industries

6.

Over the decades, the

hydrodynamic characteristics of slurry systems gathered considerable research and attracted the attention of multiphase flow reactor investigation 4, 7. The flow behaviour can significantly affect the performance of the slurry column. Accurate prediction of operating variables that govern the flow characteristics in multiphase flow is a necessary step for efficiency and design evaluation. The presence of particles substantially impacts the hydrodynamics of the multiphase system, but their effects are still not completely understood 8. Hydrodynamic characteristics of a slurry column are highly dependent on the operating variables (superficial gas velocity, superficial slurry velocity), physical properties (slurry density, slurry viscosity, slurry surface tension) and geometric variable (column length, diameter). Moreover, the efficiency of the slurry column is crucially dependent on the bubble size, rising velocity of the bubble and bubble wake phenomena 9. The comprehensive understanding of bubble generation technique, bubble size & its distribution is of fundamental significance, and they essentially affect the transport efficiency of the multiphase system. Though, in comparison with the understanding of bubbles size and its distribution in the Newtonian system, essential information about the effect of particle concentration, particle size, viscosity, axial height, and superficial gas velocity in the presence and absence of surfactant still needs to be understood in three-phase systems. As a result of the complicated flow behaviour of bubbles and the effect of the physical properties of the system, a thorough investigation of the bubble size and its distribution is a major challenge. Bubble size and its distribution significantly affect the hydrodynamics

10

and mass transfer characteristics of a slurry column. For scale-up,

modelling, and accurate design of the column it is essentially important to understand the flow and transport characteristics between different phases by studying the bubble size and its distribution

9, 11.

Numerous techniques for the measurement of bubble size and its distribution

have been used by many researchers. Some of the techniques as Mie scattering technique 2 ACS Paragon Plus Environment

12,

Page 3 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

pressure fluctuations 13, Bayesian magnetic resonance technique 14, resistivity probes 15, Coulter count method or pore electrical resistance method method

18,

photographic technique

19, 20,

16,

optical probe

and x-ray tomography

21

17,

two-point conductivity

are used for estimation of

bubble size. Intrusive techniques (e.g., optical probe) disturb the flow hydrodynamics, therefore non-intrusive methods (e.g., photographic technique) are favoured over intrusive technique

22.

The simplest technique to obtain the bubble size from an image is to manually draw lines on bubble interface 23. Interfacial area available for a process is a crucial variable for the design and scale-up purpose. The lower the bubble size is, the higher the interfacial area is, which improves the process performance. The interfacial area depends on the bubble size, gas holdup, sparger design, property of phases, and column geometry. Bubble size in the three-phase system depends on the physico-chemical property of the system and the type of gas distributor used 24. Passos et al.,

25

reported the measurement of the bubble size by an optical technique in a non-Newtonian solution. Addition of non-ionic surfactant (Triton X-100) results in a reduction in bubble size as well as the transition of homogenous to the heterogeneous regime is shifted to a higher gas flow rate. When smaller bubbles coalesce to form a bigger bubble, surface area for bubble-particle interaction decreases, thus process efficiency reduces

26.

Grau and Heiskanen

27

studied the

bubble size distribution and observed that reagents control the bubble size by decreasing the coalescence behaviour. The mean bubble size was found to increase with quartz particle addition and an increase in air flow rate. Yianatos et al.,

28

observed increased slurry density by an

increase in the particle concentration, which results in an increase in the bubble rise velocity 29. The rising velocity of the larger bubble is higher than the smaller bubble

30, 31,

hence the

residence time of the bubble decreases which impact the system performance 32. Since the past few years, several efforts have been carried out to investigate the bubble size and its distributions in the slurry system. The influence of slurry concentration and particle diameter on the bubble size has been studied by many investigators

30, 31, 33.

Numerous researchers

deduced that the addition of particles causes generation of bigger bubbles 11, 30-34. The increase in bubble size was because of the rise in the viscosity and mixture density of the slurry after particle addition

31.

The increase in liquid viscosity generally increases the bubble coalescence rate,

consequently enlarge the gas bubble size

34-36.

In contrast, Götz et al., 3

ACS Paragon Plus Environment

8

reported that the small

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 59

concentration of fine particles ( 0. The shape parameter and scale 36 ACS Paragon Plus Environment

Page 37 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

parameters are obtained from the experimental bubble size. The log-logistic distributions provide the closest match to the size distribution of the bubbles in (Air-water + SDS + particle) which can be represented as

 F d b   

 db     

 1

  d   1   b         

2

(13)

where  is the bubble shape parameter and  is the scale parameter. The shift parameter and scale parameter are obtained from the experiment. The experimental BSD in (Air-water + Tween 20 + particle) are well represented by a beta distribution function in agreement with the findings of many other researchers

F (d b ) 

1 (d b  a ) p 1 (b  d b ) q 1 B ( p, q ) (b  a ) p  q 1

(14)

where a is the minimum value of db, b is the maximum value of db, p is the lower shape parameter > 0, q is the upper shape parameter > 0, and B (p, q) is the beta function. The parameters b, p, and q are obtained from the experimental conditions. Generalized correlations based on the experimental conditions are developed for all the parameters of the log-logistic, beta and Weibull distribution function. Developed correlations for distribution function parameters are reported in Table 5. Table 5 Developed correlations for distribution parameters for different system. System

Correlation for distributions parameters

Air-water+ SDS+ particle (log-logistics distribution)

c    775.11 s   m 

0.477

 u sg d p  m       eff  

1.363

 u sg2     gd   p

0.092

 u sg2 d p  m       l  

Region of measure ment A

0.887

,

R2= 0.95

c    24.49 s   m 

0.176

 u sg d p  m       eff  

0.593

 u sg2     gd   p

0.419

0.98

37 ACS Paragon Plus Environment

 u sg2 d p  m       l  

0.260

, R2= (15a)

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

c    1.192 s   m 

0.473

1.367

 u sg d p  m       eff  

 u sg2     gd  p  

Page 38 of 59

 u sg2 d p  m       l  

1.085

 u sg2 d p  m       l  

0.0002

0.410

B R2=

,

0.94

c    0.32 s   m 

0.079

 u sg d p  m       eff  

0.067

 u sg2     gd  p  

0.354

, R2=

(15b)

0.98

c    1.192 s   m 

0.187

 u sg d p  m       eff  

0.373

 u sg2     gd   p

0.313

 u sg2 d p  m       l  

C

0.158

, R2=

0.94

c    0.32 s   m 

0.343

 u sg d p  m       eff  

0.907

 u sg2     gd   p

0.183

 u sg2 d p  m       l  

0.599

, R2= (15c)

0.92 For Airwater+Tween20 +particle (beta distribution)

c  b  1.03  10  s   m 

0.569

4

 u sg d p  m       eff  

2.179

1.066

 u sg2     gd   p

 u sg2 d p  m       l  

A

1.963

,

R2 = 0.98

c  p  2.12  10  s   m 

0.1074

3

1.063

 u sg d p  m       eff  

 u sg2     gd   p

0.426

 u sg2 d p  m       l  

0.916

,

R2 = 0.94

c  q  1.99 s   m 

0.10

 u sg d p  m       eff  

c  b  2.78  103  s   m 

0.534

0.051

 u sg2     gd   p

 u sg d p  m       eff  

0.238

1.407

 u sg2 d p  m       l  

u     gd  p   2 sg

0.133

0.319

, R2 = 0.98 (16a) B

1.204

 u d p m       l   2 sg

,

R2= 0.97

c  p  19.23 s   m 

0.138

 u sg d p  m       eff  

0.530

 u sg2     gd  p  

0.123

 u sg2 d p  m       l  

0.142

, R2=

0.96

c  q  4.18 s   m 

0.006

 u sg d p  m       eff  

0.055

 u sg2     gd  p  

0.268

 u sg2 d p  m       l  

0.169

, R2= (16b)

0.99

c  b  7.73  10  s   m  4

0.80

 u sg d p  m       eff  

2.260

 u sg2     gd  p  

R2= 0.95

38 ACS Paragon Plus Environment

0.469

C

1.844

 u sg2 d p  m       l  

,

Page 39 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

c  p  36.87 s   m 

0.272

 u sg d p  m       eff  

0.672

 u sg2     gd  p  

0.192

 u sg2 d p  m       l  

0.434

, R2=

0.94

c  q  1.21  10  s   m 

0.562

4

 u sg d p  m       eff  

1.872

 u sg2     gd  p  

0.471

1.405

 u sg2 d p  m       l  

, (16c)

R2= 0.97 Air-water + particle (Weibull distribution)

c    408.844 s   m 

0.406

 u sg d p  m       eff  

1.106

0.981

 u sg2 d p  m       l  

 u sg2     gd   p

A

0.418

,

R2= 0.91

c    2.37  10  s   m 

1.895

9

 u sg d p  m       eff  

5.283

 u sg2 d p  m       l  

3.929

 u sg2     gd   p

0.936

(17a)

, R2= 0.99

c    8.65  103  s   m 

0.604

 u sg d p  m       eff  

1.776

1.484

 u sg2 d p  m       l  

 u sg2     gd   p

B

0.518

,

R2= 0.89

c    2.54  105  s   m 

1.038

 u sg d p  m       eff  

2.854

 u sg d p  m       eff  

3.050

 u sg2 d p  m       l  

2.198

 u sg2 d p  m       l  

2.342

 u sg2     gd   p

0.405

 u sg2     gd  p  

0.797

(17b)

, R2= 0.96

c    6.90  10  s   m 

1.072

5

C

, R2= 0.93

c    1403.12 s   m 

0.399

 u sg d p  m       eff  

1.361

1.071

 u sg2 d p  m       l  

R2= 0.95

 u sg2     gd   p

0.046

, (17c)

Typical results of a comparison of experimental and predicted cumulative distribution function (CDF) are reported in Figure 17. Experimental values were calculated from the Eq. 12, Eq. 13, and Eq. 14 and predicted values were obtained from the developed correlations provided in the Table 5. The Figure 17a, reports the log-logistic distribution function where comparison between experimentally measured quantities (Eq. 13) and predicted values (Eq. 15a) are shown in system 1 (Table 4) with the AARE of 5.60% (SDS, 1 ppm, dp: 242.72 m) and 22.16% (SDS, 6 ppm, dp: 408.31 m). In Figure 17b, experimental (Eq. 14) and predicted quantities (Eq. 16a and Eq. 16b)

39 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

are obtained for system 3, in which AARE is found to be 25.28% (Tween 20, 1 ppm, dp: 408.31 m) and 9.46% (Tween 20, 6 ppm, dp: 242.72 m). Similarly, in Figure 17c, Weibull distribution function for the slurry system is reported in which experimental quantities are calculated using Eq. 12 whereas predicted values are obtained from Eq. 17a and Eq. 17c. In region A, AARE was found to be 6.60%, whereas in region C it was 14.92%. It can be concluded that the developed empirical correlations in Table 4 are suitable to predict the bubble size distribution function well within AARE of 30%. 1.05 SDS (1 ppm) dp: 242.72 m

Loglogistic distribution function

Beta distribution function

Experimentally obtained from Eq. 14 Predicted from Eq. 16a Experimentally obtained from Eq. 14 Predicted from Eq. 16b

1.0

Qg: 1 l/min

Region A AARE: 5.60% 0.75 ws: 0.2 wt.%

0.8

SDS (6 ppm) dp: 408.31 m

0.60

Qg: 5 l/min

CDF (-)

0.90

CDF (-)

Region A AARE: 22.16% ws: 0.8 wt.%

0.45

0.6

Tween 20 (6 ppm) dp: 242.72 m

0.4

Ql: 4 l/min

0.35

0.70

1.05

1.40

1.75

2.10

2.45

2.80

3.15

Ql: 2 l/min Region A AARE: 25.28%

0.2

Experimentally obtained from Eq. 13 Predicted from (Eq. 15a) Experimentally obtained from Eq. 13 Predicted from (Eq. 15a)

0.00 0.00

ws: 0.2 wt.%

Region B AARE: 9.46%

0.15

(a)

Tween 20 (1 ppm) dp: 408.31 m

ws: 0.8 wt.%

0.30

0.0

(b) 0.35

3.50

0.70

1.05

1.40

1.75

2.10

2.45

2.80

3.15

3.50

db (mm)

db (mm)

1.0

0.8

Weibull distribution Air-water-coal ws: 0.8 wt.%

Region A AARE: 6.60%

Ql: 3 l/min

0.6

CDF (-)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 40 of 59

Region C AARE: 14.92%

0.4

Experimentally obtained from Eq. 12 Predicted from Eq. 17a Experimentallly obtained from Eq. 12 Predicted from Eq. 17c

0.2

0.0

(c)

0.69

1.38

2.07

2.76

3.45

4.14

4.83

5.52

6.21

6.90

db (mm)

Figure 17. Comparison between experimental and predicted quantities (a) log logistic, (b) beta distribution, and (c) Weibull distribution. 6.4. Approach for prediction of Sauter mean bubble diameter, aspect ratio, and interfacial area 40 ACS Paragon Plus Environment

Page 41 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Bubble formation in the three-phase system is a very complex phenomena and there is no availability of reliable correlations or model which can accurately predict the bubble diameter. Hence, a novel bubble size measurement correlation in three-phase is required and can be highly useful in the slurry system. Predicting the mean bubble size is essential for increasing the understanding of the formation of bubbles in the three-phase system. In the current investigation, an effort has been made to estimate the Sauter mean bubble diameter by developing a generalized empirical correlation. Empirical correlations are developed for the three-phase system. In the three-phase system, it was observed that Sauter mean bubble diameter depends on the superficial gas velocity, particle diameter, mixture density, slurry viscosity, slurry surface tension, axial height, particle concentration, and acceleration due to gravity. The general functional form considering the operating variables and physical properties of the system is expressed as

d 32 f u sg , d p ,  m ,  sl ,  sl , z, cs , g 

(18)

Using Buckingham’s pi theorem of dimensional analysis, the following relation is obtained

u d   d 32  A  sg p m  Z   sl 

b

 u sg2 d p  m       sl  

c

 u sg2     gd   p

d

 z    d   p

e

 cs      m

f

(19)

The values of A, b, c, d, e, and f were obtained by multiple linear regression analysis. An empirical correlation for the three-phase system can be expressed as

u d   d 32  1.79  10 11  sg p m  Z   sl 

4.033

 u sg2 d p  m       sl  

3.408

1.786

 u sg2     gd   p

 z    d   p

0.936

 cs      m

0.006

(20)

The Eq. (20) is valid for a definite range of the following parameters:

2.74  (u sg d p  m /  sl )  18.20 ;

141.18  ( z / d p )  1100 ;

2.49  10 4  (u sg2 d p  m /  sl )  9.02  10 3 ;

1.50  10 2  (u sg2 / gd p )  5.34  10 1 ;

1.90  10 3  cs /  m   1.13  10 2 .

41 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

The developed correlation (Eq. 20) fits the experimental data well with the correlation coefficient (R2 = 0.99) and standard error (0.070) and suitable for prediction of Sauter mean bubble diameter in three-phase with AARE of 5.59%. Developed correlation is suitable to predict the experimentally measured Sauter mean bubble diameter within the error range of 16.93%. The experimental values of Sauter mean bubble diameter were also compared with the correlations of Pohorecki et al., 97, Wilkinson et al., 98, Zhang et al., 99, and Lemoine et al.,

100

to evaluate their

range of validity. The comparison of the present correlation (Eq. 20) with those correlations is shown in Figure 18. It is noticed that correlations suggested by Pohorecki et al., 97 and Wilkinson et al., 98 overpredicts the Sauter mean bubble diameter. As per correlation given by Zhang et al., 99,

some values are in the range of present experimental data while the other data points are

underpredcited. The overprediction and underprediction of the values are due to variation in experimental conditions. The important relevant correlations suggested by different investigators, which were compared with our present developed correlation for Sauter mean bubble diameter is given in the Table S2 as supporting information.

5

34 6. +1

%

er

r ro

4

d32- predicted (mm)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 42 of 59

%

.92 -16

3

e rr

or

Air-water (dp = 242.72 m) Air-water (dp = 408.31 m) SDS (dp = 242.72 m) SDS (dp = 408.31 m)

2

Tween 20 (dp = 242.72 m) Tween 20 (dp = 408.31 m) 97

Pohorecki et al. 98 Wilkinson et al. 99 Zhanag et al. 100 Lemoine et al.

1

AARE: 5.59% (Eq. 20)

0 0.0

0.5

1.0

1.5

2.0

2.5

3.0

3.5

d32- experimental (mm)

42 ACS Paragon Plus Environment

4.0

4.5

5.0

5.5

Page 43 of 59

Figure 18. Parity of Sauter mean bubble diameter in the three-phase system. The aspect ratio of a gas bubble can be expressed as the ratio of its minor and major axis. Estimation of aspect ratio is a matter of concern for assessment of bubble size and its distribution, flow regime transition101, and deformation of bubbles during the rising process102. Many researchers have attempted to make a correlation that relates the aspect ratio () with a dimensionless number. Some researchers have established a correlation based on Eotvos number (Eo) 101-104, Weber number (We)104-108, whereas some of the other researchers used dimensionless group which is the function of Reynolds number (Re) and Morton number (Mo)109 and Eotvos number110. A typical plot of the aspect ratio of gas bubbles for a discrete bubbly flow regime at different experimental condition is shown in Figure 19a  19d. It is observed that aspect ratio values of bubble shape are distributed widely. Based on the aspect ratio, the size discretizes to different equivalent bubble size classes. Each bubble size class represents the average aspect ratio of the gas bubbles belonging to that class. AR = 6 Qg = 3 l/min

1.0

AR: 7.8 Qg: 3 l/min

1.2

0.9 1.0

 (-)

0.8

 (-)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

0.7

0.6

0.6

0.5

0.4 (a)

0.8

0.4 (b)

1.26

1.68

2.10

2.52

2.94

3.36

3.78

4.20

4.62

1.44

1.80

2.16

2.52

2.88

deq (mm)

deq (mm)

43 ACS Paragon Plus Environment

3.24

3.60

3.96

4.32

Industrial & Engineering Chemistry Research

1.1

1.1

AR: 6 Qg: 5 l/min

1.0

AR: 7.8 Qg: 5 l/min

1.0

0.9

0.9

0.8

0.8

0.7

 (-)

 (-)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 44 of 59

0.7

0.6

0.6 0.5

0.5 0.4

0.4 (c)

0.3 1.56

1.95

2.34

2.73

3.12

3.51

3.90

4.29

4.68

1.98

2.31

2.64

2.97

deq (mm)

3.30

3.63

3.96

4.29

4.62

deq (mm)

Figure 19. Distribution of aspect ratio of bubbles at different experimental conditions. As per report of Besagni and Inzoli

101

the aspect ratio can be well predicted statistically by

single dimensionless number only instead of considering various dimensionless numbers simultaneously. Therefore, based on the experimental data and numerical fitting, a correlation for prediction of aspect ratio is proposed in terms of Weber number, which is expressed as 

1

(21)

1.0012

1.2  0.0099We

Experimentally calculated aspect ratio is also compared with the correlations of Moore Taylor and Acrivos Maeda

109,

Li et al.,

106,Wellek 102,

et al.,

104,Fan

and Tsuchiya

and Besagni and Inzoli

101

correlations except the correlation of Li et al.,

103,Okawa

et al.,

111,

105,

Tadaki and

as shown in Figure 20. It is observed that all 102

either over predict or underpredict the

experimental gas bubbles aspect ratio. The correlation of Li et al.,

102

predicts the aspect ratio

with the AARE of 6.79%, whereas the present correlation (Eq. 21) predicts within the AARE of 4.48%. The correlations for the bubble aspect ratio against which the present developed correlation is compared is provided in the Supporting Information as a Table S3.

44 ACS Paragon Plus Environment

Page 45 of 59

1.1 1.0 0.9 0.8 0.7 0.6

 (-)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

Present experiemntal data Present correlation (Eq. 21) 105 Moore et al. 106 Taylor and Acrivos 104 Wellek et al. 103 Fan and Tsuchiya 111 Okawa et al. 109 Tadaki and Maeda 102 Li et al. 101 Besagni and Inzoli

0.5 0.4 0.3 0.2 0.1 0.0 1.5

2.0

2.5

3.0

3.5

4.0

4.5

5.0

deq (mm) Figure 20. Comparison of the experimental aspect ratio with predicted and others available aspect ratio correlation. Similarly, a generalized correlation for prediction of interfacial area of gas bubble is also developed by applying Buckingham’s Pi theorem. The generalized correlation for the interfacial area is represented as

u d   ai  1.61  1015  sg p m    sl 

4.783

 u sg2 d p  m       sl  

3.762

 u sg2     gd   p

1.126

 cs     m 

0.146

 z    d   p

0.080

(22)

The correlation coefficient and standard error for developed correlation (Eq. 22) are 0.92 and 0.174 respectively. The Eq. (22) is suitable to estimate the interfacial area within the error range of  29.90% and AARE of 13.90%. The developed correlation (Eq. 22) is also compared with the correlations present in the literature

47, 112, 113

to analyse their range of validity as shown in

Figure 21. All the literature correlations for the prediction of bubble interfacial area against which the current developed correlation is compared is provided in the Supporting Information 45 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research

(Table S4). The Eq. (22) is valid for a definite range of the following parameters:

2.90  (u sg d p  m /  sl )  18.21 ;

141.18  ( z / d p )  1100 ;

2.50  10 4  (u sg2 d p  m /  sl )  9.02  10 3 ;

1.50  10 2  (u sg2 / gd p )  5.34  10 1 ;

1.99  10 3  cs /  m   1.17  10 2 . 900 Air-water (dp = 242.72 m) Air-water (dp = 408.31 m)

800

er ro r

SDS (dp = 242.72 m)

700

%

SDS (dp = 408.31 m)

.7 2

Tween 20 (dp = 242.72 m)

29

Tween 20 (dp = 408.31 m) 47

2

+

Besagni and Inzoli 112 Akita and Yoshida 113 Gestrich and Krauss

600

3

ai- predicted (m /m )

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 46 of 59

500

9 -2

400 300

.90

%

or err

AARE: 13.90% (Eq. 22)

200 100 0 0

100

200

300

400

500

600 2

700

800

900

3

ai- experimental (m /m ) Figure 21. Comparison of experimental and predicted interfacial area in the three-phase system.

Conclusion A photographic technique was used to measure the bubble size and its distribution in three axial regions of the microstructured slurry column. A detailed study was conducted to enunciate the effect of particle concentration, particle size, surfactant type, and surfactant doses on bubble size, bubble size distribution, and the type of the distribution function in a discrete bubbly flow 46 ACS Paragon Plus Environment

Page 47 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

regime. A linear relationship was observed between superficial gas velocity and bubble size. Addition of surfactant inhibits the rate of coalescence of bubbles, which result in smaller bubble size and narrower bubble size distribution compared to three-phase system. A slight difference in bubble size distribution observed when solid loading increases from 0.2 to 0.8 wt.% in the presence of the surfactant. The bubble diameter increased with increasing distance from the base of the column due to the coalescence of bubbles. The size distributions of bubbles in the presence of solid particles were wider than the two-phase system. It was also observed that as axial height increases the bubble size distribution curve shifted from smaller bubble diameter to the larger bubble diameter. In all systems, the obtained bubble size distributions were not symmetrical as normal distributions. The obtained bubble size distribution observed to follow distribution function, namely, beta, log-logistic, and Weibull in different experimental conditions. Correlations were developed for the estimation of each distribution parameters of all the distribution functions for the three axial regions of the column. A generalized correlation model is developed to predict the Sauter mean bubble diameter, bubble aspect ratio, and bubble interfacial area in the three-phase system and also compared with the published correlation in the literature. All the developed correlations are suitable to predict the experimental data well within the error range of  29.90% and AARE of 13.90%. The present experiment will be helpful for process intensification of a multiphase system. Supporting Information: 

Table S1 shows discretization of bubble classes in the axial region B in the presence of Tween 20 surfactant at 0.008 m/s superficial gas velocity is shown in Table S1



Table S2 shows the different correlations proposed by various investigators for Sauter mean bubble diameter against which the present proposed correlation is compared



Table S3 shows the different correlations proposed by various investigators for aspect ratio against which the present proposed correlation is compared



Table S4 shows the different correlations for the bubble interfacial area against which the present proposed correlation is compared

Nomenclature ai

Interfacial area (m2/m3) 47 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

a

Minimum value of bubble size in Eq. (14) (-)

AC

Cross section area of the column (m2)

AARE

Absolute average relative error (%)

AR

Aspect ratio (-)

b

Maximum value of bubble size in Eq. (14) (-)

B (p, q)

Beta Function in Eq. (14) (-)

BSD

Bubble size distribution (-)

CDF

Cumulative distribution function (-)

cs

Particle concentration (kg/m3)

cf

Surfactant concentration (ppm)

cv

Volumetric particle concentration (vol.%)

dp

Particle diameter (m)

db

Bubble diameter (m)

dc

Column diameter (m)

d32

Sauter mean bubble diameter (m)

deq

Equivalent bubble diameter (m)

d0

Sparger diameter (m)

F(db)

Bubble distribution function (-)

g

Acceleration due to gravity (m/s2)

hm

Gas-liquid-solid mixture height (m)

L

Column length (m) 48 ACS Paragon Plus Environment

Page 48 of 59

Page 49 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

lmajor

Maximum axial length of the bubble (m)

lminor

Minimum axial length of the bubble (m)

ms

Mass of solid (kg)

Mw-gas

Molecular weight of gas (kg/kmol)

n

Flow behaviour index (-)

ni

Number of bubbles (-)

N

Total number of variables (-)

p

Lower shape parameter in Eq. (14) (-)

q

Upper shape parameter in Eq. (14) (-)

R2

Correlation coefficient (-)

RF

Relative frequency (-)

RPM

Rotation per minute (-)

STDEV

Standard deviation (-)

T

Temperature (k)

usg

Superficial gas velocity (m/s)

U

Standard uncertainty (-)

Ur

Relative uncertainty (-)

ws

Particle concentration (wt.%)

Xw

Concentration of purest component in the liquid mixture (wt.%)

xi

ith component of the variables (-) 49 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

x

Average of data points (-)

z

Axial height (m)

Non-dimensional number Eo

2 Eotvos number g (  l   g )d eq /  l (-)

Mo

4 2 3 Morton number g (  l   g )  l /  l  l (-)

Re

Reynolds number  sl u sg d eq / l ) (-)

Fr

2 Froude number (u sg / d eq g ) (-)

We

2 Weber number (  l u sg d eq /  l ) (-)

Greek Letter 

Aspect ratio of bubble (-)

g

Gas holdup (-)



Bubble shape parameter in Eq. (12) (-)

l

Liquid holdup (-)

s

Solid holdup (-)



Energy dissipation rate per unit mass (m2/s3)

β

Scale parameter of bubble in Eq. (12) (-)

l

Surface tension of liquid (kg/m.s)

sl

Surface tension of slurry (kg/m.s)

m

Mixture (gas-liquid-solid) density (kg/m3)

50 ACS Paragon Plus Environment

Page 50 of 59

Page 51 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

l

Density of liquid (kg/m3)

s

Density of particle (kg/m3)

sl

Density of slurry (kg/m3)

g

Density of gas (kg/m3)

eff

Effective viscosity (Pa.s)

l

Viscosity of liquid (Pa.s)

References (1) Kanaris, A. G.; Pavlidis, T. I.; Chatzidafni, A. P.; Mouza, A. A., The Effects of the Properties of Gases on the Design of Bubble Columns Equipped with a Fine Pore Sparger. ChemEngineering 2018, 2 (1), 11. (2) Sarhan, A. R.; Naser, J.; Brooks, G., Effects of particle size and concentration on bubble coalescence and froth formation in a slurry bubble column. Particuology 2018, 36, 82-95. (3) Deckwer, W. D., Bubble Column Reactors. Wiley: Chichester: UK, 1992. (4) Sarhan, A. R.; Naser, J.; Brooks, G., CFD simulation on influence of suspended solid particles on bubbles' coalescence rate in flotation cell. Int. J. Miner. Process. 2016, 146, 54-64. (5) Finch, J. A.; Dobby, G. S., Column Flotation. Pergamon Press: Oxford: UK, 1990. (6) Nedeltchev, S., Theoretical prediction of mass transfer coefficients in both gas–liquid and slurry bubble columns. Chem. Eng. Sci. 2017, 157, 169-181. (7) Ojima, S.; Hayashi, K.; Tomiyama, A., Effects of hydrophilic particles on bubbly flow in slurry bubble column. Int. J. Multiphase Flow 2014, 58, 154-167. (8) Götz, M.; Lefebvre, J.; Mörs, F.; Reimert, R.; Graf, F.; Kolb, T., Hydrodynamics of organic and ionic liquids in a slurry bubble column reactor operated at elevated temperatures. Chem. Eng. J. 2016, 286, 348-360. (9) Nedeltchev, S.; Nigam, K. D. P.; Schumpe, A., Prediction of mass transfer coefficients in a slurry bubble column based on the geometrical characteristics of bubbles. Chem. Eng. Sci 2014, 106, 119-125. (10) Moshtari, B.; Babakhani, E. G.; Moghaddas, J. S., Experimental study of gas hold-up and bubble behavior in gas-liquid bubble column. Pet. Coal 2009, 51 (1), 27-32. 51 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(11) Li, H.; Prakash, A., Influence of slurry concentrations on bubble population and their rise velocities in a three-phase slurry bubble column. Powder Technol. 2000, 113 (1-2), 158-167. (12) Hansen, G. M., Mie scattering as a technique for the sizing of air bubbles. Appl. Opt. 1985, 24 (19), 3214-3220. (13) Chilekar, V. P.; Warnier, M. J. F.; Van Der Schaaf, J.; Kuster, B. F. M.; Schouten, J. C.; Van Ommen, J. R., Bubble size estimation in slurry bubble columns from pressure fluctuations. AIChE J. 2005, 51 (7), 1924-1937. (14) Holland, D. J.; Blake, A.; Tayler, A. B.; Sederman, A. J.; Gladden, L. F., Bubble size measurement using Bayesian magnetic resonance. Chem. Eng. Sci. 2012, 84, 735-745. (15) Zhang, K.; Qi, N.; Jin, J.; Lu, C.; Zhang, H., Gas holdup and bubble dynamics in a threephase internal loop reactor with external slurry circulation. Fuel 2010, 89 (7), 1361-1369. (16) Fujiwara, H., Spectroscopic ellipsometry: principles and applications. John Wiley & Sons: 2007. (17) Xue, J.; Al‐Dahhan, M.; Dudukovic, M. P.; Mudde, R. F., Bubble velocity, size, and interfacial area measurements in a bubble column by four‐point optical probe. AIChE J. 2008, 54 (2), 350-363. (18) Buchholz, R.; Zakrzewski, W.; Schugerl, K., Techniques for determining the properties of bubbles in a bubble columns. Chem. Ing. Tech. 1979, 51 (6), 568-575. (19) Lau, R.; Mo, R.; Sim, W. S. B., Bubble characteristics in shallow bubble column reactors. Chem. Eng. Res. Des. 2010, 88 (2), 197-203. (20) Prakash, R.; Majumder, S. K.; Singh, A., Flotation technique: Its mechanisms and design parameters. Chem. Eng. Process. 2018, 127, 249-270. (21) Rabha, S.; Schubert, M.; Wagner, M.; Lucas, D.; Hampel, U., Bubble size and radial gas hold‐up distributions in a slurry bubble column using ultrafast electron beam X‐ray tomography. AIChE J. 2013, 59 (5), 1709-1722. (22) Guet, S.; Luther, S.; Ooms, G., Bubble shape and orientation determination with a fourpoint optical fibre probe. Exp. Therm. Fluid Sci. 2005, 29 (7), 803-812. (23) Winkel, E. S.; Ceccio, S. L.; Dowling, D. R.; Perlin, M., Bubble-size distributions produced by wall injection of air into flowing freshwater, saltwater and surfactant solutions. Exp. Fluids 2004, 37 (6), 802-810. (24) Camarasa, E.; Vial, C. H.; Poncin, S.; Wild, G.; Midoux, N.; Bouillard, J., Influence of coalescence behaviour of the liquid and of gas sparging on hydrodynamics and bubble characteristics in a bubble column. Chem. Eng. Process. 1999, 38 (4-6), 329-344.

52 ACS Paragon Plus Environment

Page 52 of 59

Page 53 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

(25) Passos, A. D.; Voulgaropoulos, V. P.; Paras, S. V.; Mouza, A. A., The effect of surfactant addition on the performance of a bubble column containing a non-Newtonian liquid. Chem. Eng. Res. Des. 2015, 95, 93-104. (26) Ata, S.; Ahmed, N.; Jameson, G. J., A study of bubble coalescence in flotation froths. Int. J. Miner. Process. 2003, 72 (1-4), 255-266. (27) Grau, R. A.; Heiskanen, K., Bubble size distribution in laboratory scale flotation cells. Miner. Eng. 2005, 18 (12), 1164-1172. (28) Yianatos, J. B.; Finch, J. A.; Dobby, G. S.; Xu, M., Bubble size estimation in a bubble swarm. J. Colloid Interface Sci. 1988, 126 (1), 37-44. (29) Bhunia, K.; Kundu, G.; Mukherjee, D., Gas holdup characteristics in a flotation column with different solids. Sep. Sci. Technol. 2017, 52 (7), 1298-1309. (30) Moshtari, B.; Moghaddas, J. S.; Gangi, E., A hydrodynamic experimental study of slurry bubble column. Stud. Surf. Sci. Catal. Elsevier: 2007; Vol. 167, pp 67-72. (31) Manjrekar, O. N.; Dudukovic, M. P., Application of a 4-point optical probe to a Slurry Bubble Column Reactor. Chem. Eng. Sci. 2015, 131, 313-322. (32) Luo, X.; Lee, D. J.; Lau, R.; Yang, G.; Fan, L. S., Maximum stable bubble size and gas holdup in high‐pressure slurry bubble columns. AIChE J. 1999, 45 (4), 665-680. (33) Tyagi, P.; Buwa, V. V., Dense gas–liquid–solid flow in a slurry bubble column: Measurements of dynamic characteristics, gas volume fraction and bubble size distribution. Chem. Eng. Sci. 2017, 173, 346-362. (34) Li, H.; Prakash, A., Heat transfer and hydrodynamics in a three-phase slurry bubble column. Ind. Eng. Chem. Res. 1997, 36 (11), 4688-4694. (35) Sehabiague, L.; Morsi, B. I., Hydrodynamic and mass transfer characteristics in a largescale slurry bubble column reactor for gas mixtures in actual Fischer–Tropsch cuts. Int. J. Chem. Reactor Eng. 2013, 11 (1), 83-102. (36) Kato, Y.; Nishiwaki, A.; Fukuda, T.; Tanaka, S., The behavior of suspended solid particles and liquid in bubble columns. J. Chem. Eng. Jpn. 1972, 5 (2), 112-118. (37) Rabha, S.; Schubert, M.; Hampel, U., Intrinsic flow behavior in a slurry bubble column: a study on the effect of particle size. Chem. Eng. Sci. 2013, 93, 401-411. (38) Banisi, S.; Finch, J. A.; Laplante, A. R.; Weber, M. E., Effect of solid particles on gas holdup in flotation columns—I. Measurement. Chem. Eng. Sci. 1995a, 50 (14), 2329-2334. (39) Banisi, S.; Finch, J. A.; Laplante, A. R.; Weber, M. E., Effect of solid particles on gas holdup in flotation columns—II. Investigation of mechanisms of gas holdup reduction in presence of solids. Chem. Eng. Sci. 1995b, 50 (14), 2335-2342. 53 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(40) Vazirizadeh, A.; Bouchard, J.; Chen, Y., Effect of particles on bubble size distribution and gas hold-up in column flotation. Int. J. Miner. Process. 2016, 157, 163-173. (41) Wang, T.; Wang, J.; Yang, W.; Jin, Y., Bubble behavior in gas–liquid–solid three-phase circulating fluidized beds. Chem. Eng. J. 2001, 84 (3), 397-404. (42) Zhang, L.-j.; Li, T.; Ying, W.-y.; Fang, D.-y., Rising and descending bubble size distributions in gas–liquid and gas–liquid–solid slurry bubble column reactor. Chem. Eng. Res. Des. 2008, 86 (10), 1143-1154. (43) Ojima, S.; Sasaki, S.; Hayashi, K.; Tomiyama, A., Effects of particle diameter on bubble coalescence in a slurry bubble column. J. Chem. Eng. Jpn. 2015, 48 (3), 181-189. (44) Shukla, S. C.; Kukade, S.; Mandal, S. K.; Kundu, G., Coal–oil–water multiphase fuel: Rheological behavior and prediction of optimum particle size. Fuel 2008, 87 (15-16), 34283432. (45) Orvalho, S.; Hashida, M.; Zednikova, M.; Stanovsky, P.; Ruzicka, M. C.; Sasaki, S.; Tomiyama, A., Flow Regimes in Slurry Bubble Column: Effect of Column Height and Particle Concentration. Chem. Eng. J. 2018. (46) Kim, J. Y.; Kim, B.; Nho, N.-S.; Go, K.-S.; Kim, W.; Bae, J. W.; Jeong, S. W.; Epstein, N.; Lee, D. H., Gas holdup and hydrodynamic flow regime transition in bubble columns. J. Ind. Eng. Chem. 2017, 56, 450-462. (47) Besagni, G.; Inzoli, F., The effect of liquid phase properties on bubble column fluid dynamics: Gas holdup, flow regime transition, bubble size distributions and shapes, interfacial areas and foaming phenomena. Chem. Eng. Sci. 2017, 170, 270-296. (48) Mota, A.; Vicente, A. A.; Teixeira, J., Effect of spent grains on flow regime transition in bubble column. Chem. Eng. Sci. 2011, 66 (14), 3350-3357. (49) Hooshyar, N.; Hamersma, P. J.; Mudde, R. F.; van Ommen, J. R., Intensified operation of slurry bubble columns using structured gas injection. Can. J. Chem. Eng. 2010, 88 (4), 533-542. (50) Shaikh, A.; Al-Dahhan, M. H., A review on flow regime transition in bubble columns. Int. J. Chem. Reactor Eng. 2007, 5 (1). (51) Angeles, M.; Leyva, C.; Ancheyta, J.; Ramírez, S., A review of experimental procedures for heavy oil hydrocracking with dispersed catalyst. Catal. Today 2014, 220, 274-294. (52) Nedeltchev, S.; Shaikh, A.; Al‐Dahhan, M., Flow regime identification in a bubble column based on both statistical and chaotic parameters applied to computed tomography data. Chem. Eng. Technol. 2006, 29 (9), 1054-1060. (53) Shiea, M.; Mostoufi, N.; Sotudeh-Gharebagh, R., Comprehensive study of regime transitions throughout a bubble column using resistivity probe. Chem. Eng. Sci. 2013, 100, 1522. 54 ACS Paragon Plus Environment

Page 54 of 59

Page 55 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

(54) Mena, P. C.; Ruzicka, M. C.; Rocha, F. A.; Teixeira, J. A.; Drahoš, J., Effect of solids on homogeneous–heterogeneous flow regime transition in bubble columns. Chem. Eng. Sci. 2005, 60 (22), 6013-6026. (55) Gan, Z., Holdup and velocity profiles of monosized spherical solids in a three-phase bubble column. Chem. Eng. Sci. 2013, 94, 291-301. (56) Li, W.; Zhong, W.; Jin, B.; Lu, Y.; He, T., Flow patterns and transitions in a rectangular three-phase bubble column. Powder Technol. 2014, 260, 27-35. (57) Hooshyar, N.; Hamersma, P. J.; Mudde, R. F.; van Ommen, J. R., Gas fraction and bubble dynamics in structured slurry bubble columns. Ind. Eng. Chem. Res. 2010, 49 (21), 1068910697. (58) Ruthiya, K. C.; Chilekar, V. P.; Warnier, M. J.; van der Schaaf, J.; Kuster, B. F.; Schouten, J. C.; van Ommen, J. R., Detecting regime transitions in slurry bubble columns using pressure time series. AIChE J. 2005, 51 (7), 1951-1965. (59) Nygren, A., Simulation of bubbly flow in a flat bubble column. 2014. (60) Sivaiah, M.; Majumder, S. K., Gas holdup and frictional pressure drop in a modified slurry bubble column. Int. J. Chem. Reactor Eng. 2012, 10 (1). (61) Brian, B. W.; Chen, J. C., Surface tension of solid‐liquid slurries. AIChE J. 1987, 33 (2), 316-318. (62) Prakash, R.; Majumder, S. K.; Singh, A., Gas holdup and frictional pressure drop contributions in microstructured two-and three-phase bubbling bed with Newtonian and nonNewtonian liquids: Effect of coarse and fine particles with surface active agent. Chem. Eng. Process. 2018, 133, 40-57. (63) Couvert, A.; Roustan, M.; Chatellier, P., Two-phase hydrodynamic study of a rectangular air-lift loop reactor with an internal baffle. Chem. Eng. Sci. 1999, 54 (21), 5245-5252. (64) Haberman, W. L.; Morton, R. K. An experimental investigation of the drag and shape of air bubbles rising in various liquids; David Taylor Model Basin Washington DC: 1953. (65) Besagni, G.; Inzoli, F.; Ziegenhein, T., Two-Phase Bubble Columns: A Comprehensive Review. ChemEngineering 2018, 2 (2), 13. (66) Chao, B. T., Motion of spherical gas bubbles in a viscous liquid at large Reynolds numbers. Phys. Fluids 1962, 5 (1), 69-79. (67) Hyndman, C. L.; Larachi, F.; Guy, C., Understanding gas-phase hydrodynamics in bubble columns: a convective model based on kinetic theory. Chem Eng. Sci. 1997, 52 (1), 63-77. (68) Kantarci, N.; Borak, F.; Ulgen, K. O., Bubble column reactors. Process Biochem. 2005, 40 (7), 2263-2283. 55 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(69) Azad, M. A. K.; Syeda, S. R., A numerical model for bubble size distribution in turbulent gas-liquid dispersion. J. Chem. Eng. 2006, 24, 25-34. (70) Geary, N. W.; Rice, R. G., Bubble size prediction for rigid and flexible spargers. AIChE J. 1991, 37 (2), 161-168. (71) Bothe, M.; Schlüter, M., Modellierungsparameter für Blasenströmungen. Chem. Ing. Tech. 2013, 7 (85), 1023-1035. (72) Besagni, G.; Di Pasquali, A.; Gallazzini, L.; Gottardi, E.; Colombo, L. P. M.; Inzoli, F., The effect of aspect ratio in counter-current gas-liquid bubble columns: Experimental results and gas holdup correlations. Int. J. Multiphase Flow 2017, 94, 53-78. (73) Wilkinson, P. M.; Spek, A. P.; van Dierendonck, L. L., Design parameters estimation for scale‐up of high‐pressure bubble columns. AIChE J. 1992, 38 (4), 544-554. (74) Shaikh, A.; Al-Dahhan, M., Scale-up of bubble column reactors: a review of current stateof-the-art. Ind. Eng. Chem. Res. 2013, 52 (24), 8091-8108. (75) Vadlakonda, B.; Mangadoddy, N., Hydrodynamic study of three-phase flow in column flotation using electrical resistance tomography coupled with pressure transducers. Sep. Purif. Technol. 2018, 203, 274-288. (76) Sada, E.; Kumazawa, H.; Lee, C. H., Influences of suspended fine particles on gas holdup and mass transfer characteristics in a slurry bubble column. AIChE J. 1986, 32 (5), 853-856. (77) Luo, X.; Yang, G.; Lee, D. J.; Fan, L.-S., Single bubble formation in high pressure liquid— solid suspensions. Powder Technol. 1998, 100 (2-3), 103-112. (78) Akita, K.; Yoshida, F., Gas holdup and volumetric mass transfer coefficient in bubble columns. Effects of liquid properties. Ind. Eng. Chem. Process Des. Dev. 1973, 12 (1), 76-80. (79) Kara, S.; Kelkar, B. G.; Shah, Y. T.; Carr, N. L., Hydrodynamics and axial mixing in a three-phase bubble column. Ind. Eng. Chem. Process Des. Dev. 1982, 21 (4), 584-594. (80) Tavera, F. J.; Escudero, R., Effect of solids on gas dispersion characteristics: Addition of hydrophobic and hydrophilic solids. J. Mex. Chem. Soc. 2012, 56 (2), 217-221. (81) Xu, X.; Zhang, J.; Liu, F.; Wang, X.; Wei, W.; Liu, Z., Rising behavior of single bubble in infinite stagnant non-Newtonian liquids. Int. J. Multiphase Flow 2017, 95, 84-90. (82) Yoshida, Y.; Katsumoto, T.; Taniguchi, S.; Shimosaka, A.; Shirakawa, Y.; Hidaka, J., Prediction of viscosity of slurry suspended fine particles using coupled DEM-DNS simulation. Chem. Eng. Trans. 2013, 32, 2089-2094. (83) Kawatra, S. K.; Eisele, T. C., Rheological effects in grinding circuits. Int. J. Miner. Process. 1988, 22, 251-259.

56 ACS Paragon Plus Environment

Page 56 of 59

Page 57 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

(84) Zhou, Z. A.; Plitt, L. R.; Egiebor, N. O., The effects of solids and reagents on the characteristics of coal flotation in columns. Miner. Eng. 1993, 6 (3), 291-306. (85) Senapati, P. K.; Panda, D.; Parida, A., Predicting viscosity of limestone–water slurry. J. Miner. Mater. Char. Eng. 2009, 8 (03), 203. (86) Smith, J. S.; Burns, L. F.; Valsaraj, K. T.; Thibodeaux, L. J., Bubble column reactors for wastewater treatment. 2. The effect of sparger design on sublation column hydrodynamics in the homogeneous flow regime. Ind. Eng. Chem. Res. 1996, 35 (5), 1700-1710. (87) Lau, Y. M.; Sujatha, K. T.; Gaeini, M.; Deen, N. G.; Kuipers, J. A. M., Experimental study of the bubble size distribution in a pseudo-2D bubble column. Chem. Eng. Sci. 2013, 98, 203211. (88) Wongsuchoto, P.; Charinpanitkul, T.; Pavasant, P., Bubble size distribution and gas–liquid mass transfer in airlift contactors. Chem. Eng. J. 2003, 92 (1-3), 81-90. (89) Parthasarathy, R.; Ahmed, N., Bubble size distribution in a gas sparged vessel agitated by a Rushton turbine. Ind. Eng. Chem. Res. 1994, 33 (3), 703-711. (90) Honkanen, M.; Saarenrinne, P.; Stoor, T.; Niinimäki, J., Recognition of highly overlapping ellipse-like bubble images. Meas. Sci. Technol. 2005, 16 (9), 1760. (91) Lage, P.; Espósito, R., Experimental determination of bubble size distributions in bubble columns: prediction of mean bubble diameter and gas hold up. Powder Technol. 1999, 101 (2), 142-150. (92) Al-Oufi, F. M. An investigation of gas void fraction and transition conditions for two-phase flow in an annular gap bubble column. © Fahd Aloufi, 2011. (93) Fan, L. S.; Yang, G. Q.; Lee, D. J.; Tsuchiya, K.; Luo, X., Some aspects of high-pressure phenomena of bubbles in liquids and liquid–solid suspensions. Chem. Eng. Sci. 1999, 54 (21), 4681-4709. (94) Sarkar, S. K.; Machniewski, P. M.; Evans, G. M., Modelling and measurement of bubble formation and growth in electroflotation processes. Chem. Process Eng. 2013, 34 (3), 327-336. (95) Wilkinson, P. M.; Dierendonck, L. L. v., Pressure and gas density effects on bubble breakup and gas hold-up in bubble columns. Chem. Eng. Sci. 1990, 45 (8), 2309-2315. (96) Cui, Z. Hydrodynamics in a bubble column at elevated pressures and turbulence energy distribution in bubbling gas-liquid and gas-liquid-solid flow systems. The Ohio State University, 2005. (97) Pohorecki, R.; Moniuk, W.; Bielski, P.; Sobieszuk, P.; Dąbrowiecki, G., Bubble diameter correlation via numerical experiment. Chem. Eng. J. 2005, 113 (1), 35-39.

57 ACS Paragon Plus Environment

Industrial & Engineering Chemistry Research 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(98) Wilkinson, P. M.; Haringa, H.; Van Dierendonck, L. L., Mass transfer and bubble size in a bubble column under pressure. Chem. Eng. Sci. 1994, 49 (9), 1417-1427. (99) Zhang, X.; Dong, H.; Huang, Y.; Li, C.; Zhang, X., Experimental study on gas holdup and bubble behavior in carbon capture systems with ionic liquid. Chem. Eng. J. 2012, 209, 607-615. (100) Lemoine, R.; Behkish, A.; Sehabiague, L.; Heintz, Y. J.; Oukaci, R.; Morsi, B. I., An algorithm for predicting the hydrodynamic and mass transfer parameters in bubble column and slurry bubble column reactors. Fuel Process. Technol. 2008, 89 (4), 322-343. (101) Besagni, G.; Inzoli, F., Bubble size distributions and shapes in annular gap bubble column. Exp. Therm. Fluid Sci. 2016, 74, 27-48. (102) Li, H.; Liu, Z.; Chen, J.; Sun, B.; Guo, Y.; He, H., Correlation of aspect ratio and drag coefficient for hydrate-film-covered methane bubbles in water. Exp. Therm. Fluid Sci. 2017, 88, 554-565. (103) Fan, L.; Tsuchiya, K., Bubble Wake Dynamics in Liquid and Liquid-Solid Suspensions,(1990). Butter-worth-Heinemann Series in Chemical Engineering. (104) Wellek, R.; Agrawal, A.; Skelland, A., Shape of liquid drops moving in liquid media. AIChE J. 1966, 12 (5), 854-862. (105) Moore, D., The rise of a gas bubble in a viscous liquid. J. Fluid Mech. 1959, 6 (1), 113130. (106) Taylor, T.; Acrivos, A., On the deformation and drag of a falling viscous drop at low Reynolds number. J. Fluid Mech. 1964, 18 (3), 466-476. (107) Raymond, F.; Rosant, J.-M., A numerical and experimental study of the terminal velocity and shape of bubbles in viscous liquids. Chem. Eng. Sci. 2000, 55 (5), 943-955. (108) Kelbaliyev, G.; Ceylan, K., Development of new empirical equations for estimation of drag coefficient, shape deformation, and rising velocity of gas bubbles or liquid drops. Chem. Eng. Commun. 2007, 194 (12), 1623-1637. (109) Tadaki, T.; Maeda, M., On the shape and velocity of single air bubbles rising in various liquids. Kagaku Kogaku 1961, 25, 254-264. (110) Gumerov, N. A.; Chahine, G. L. In Dynamics of bubbles in conditions of gas hydrate formation, The Eighth International Offshore and Polar Engineering Conference, 1998; International Society of Offshore and Polar Engineers: 1998. (111) Okawa, T.; Tanaka, T.; Kataoka, I.; Mori, M., Temperature effect on single bubble rise characteristics in stagnant distilled water. Int. J. Heat Mass Transfer 2003, 46 (5), 903-913. (112) Akita, K.; Yoshida, F., Bubble size, interfacial area, and liquid-phase mass transfer coefficient in bubble columns. Ind. Eng. Chem. Process Des. Dev. 1974, 13 (1), 84-91. 58 ACS Paragon Plus Environment

Page 58 of 59

Page 59 of 59 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Industrial & Engineering Chemistry Research

(113) Gestrich, W.; Krauss, W., Die spezifische Phasengrenzfläche in Blasenschichten. Chem. Ing. Tech. 1975, 47 (9), 360-367.

For Table of Contents Only

59 ACS Paragon Plus Environment