Estimation of Bubble Properties in Bubbling Fluidized Bed Using

May 30, 2018 - In bubbling fluidized beds, important bubble properties can be estimated using electrical capacitance volume tomography (ECVT)...
0 downloads 0 Views 3MB Size
Subscriber access provided by Kaohsiung Medical University

Thermodynamics, Transport, and Fluid Mechanics

Estimation of bubble properties in bubbling fluidized bed using ECVT measurements Vaibhav Agrawal, Yogesh H Shinde, Milinkumar T Shah, Ranjeet P. Utikar, Vishnu Pareek, and Jyeshtharaj B Joshi Ind. Eng. Chem. Res., Just Accepted Manuscript • DOI: 10.1021/acs.iecr.8b00349 • Publication Date (Web): 30 May 2018 Downloaded from http://pubs.acs.org on May 31, 2018

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 38 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

Estimation of bubble properties in bubbling fluidized bed using ECVT measurements Vaibhav Agrawala, Yogesh H Shindea, Milinkumar T Shaha*, Ranjeet P Utikara, Vishnu K Pareeka, Jyeshtharaj B Joshib a

Western Australia School of Mines: Minerals, Energy and Chemical Engineering, Curtin University, Perth, WA 6102, Australia b

Homi Bhabha National Institute, Mumbai 400 094, India

*

Corresponding author: [email protected]

Highlights 1. Bubbling fluidized bed experiments were conducted using ECVT. 2. A method to calculate the optimum threshold to post-process each 2D frame was developed. 3. Size, velocity, and frequency of bubbles were estimated and analyzed. Abstract In bubbling fluidized beds, important bubble properties can be estimated using electrical capacitance volume tomography (ECVT), resolution of which depends on sensor configuration and image reconstruction algorithm. At low resolution, the estimation of bubble size critically relies on threshold solids volume fraction. The majority of previous studies used a constant threshold, however, present study determined that the use of a constant threshold could fail to identify small bubbles and lead to an erroneous estimation of bubble size. A series of calibration experiments were conducted, and the data from which were used to develop an iterative method to determine the optimum threshold. The data from subsequent

Page 1 of 38 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 38

dynamic experiments were post-processed using LBP algorithm and the optimum threshold method. The use of optimum threshold not only resolved the smaller bubbles but it also estimated the increase in the bubble size more accurately as it moved up through the fluidized bed. Keywords: Bubbling fluidized bed, experiments, ECVT, bubble properties Abstract graphic

1

Introduction

Bubbling fluidized beds (BFBs) are widely used as a gas-solid reactor in catalytic cracking, reforming, polymerization, pyrolysis, combustion, and numerous

other

chemical

processes.

In

such

reactors,

gas-solid

hydrodynamics depend on the formation, growth, and eruption of bubbles. Thus, it is vital to understand the inherent bubble hydrodynamics for design and scale-up. The present study investigates size, velocity, and frequency of

Page 2 of 38 ACS Paragon Plus Environment

Page 3 of 38 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

bubbles in a lab-scale BFB using measurements obtained using electrical capacitance volume tomography (ECVT). BFB experiments have been conducted using different invasive and noninvasive measurement techniques. The non-invasive techniques are more desirable as they neither interfere with flow nor require transparent equipment. Typically, non-invasive techniques such as -ray tomography,

-ray tomography, electrical capacitance tomography (ECT), magnetic resonance imaging (MRI), radioactive particle tracking (RPT), and positron emission particle tracking (PEPT) are used. Each of these techniques has advantages and disadvantages in terms of operational feasibility and resolution of data.1 ECT has advantages such as fast imaging speed and low operational hazards.1 Consequently, ECT (data obtained for a single 2D plane) and ECVT (data acquired for a 3D volume) were used in several BFB studies as summarized in Table 1.2–7 Despite its extensive use, capacitance tomography also has several limitations including variation in capacitance at one location causing the signal reconstruction to be sensitive to errors and noise.1 Further, the resolution of the reconstructed images is significantly influenced by the number of electrodes and applied image reconstruction algorithm.2,8 Several algorithms such as linear back propagation (LBP)9, Landweber iterative (LI)10, Tikhonov regularization principle (TRP)11,12, algebraic reconstruction technique

(ART)13,

simultaneous

iterative

reconstruction

technique

14

(SIRT) , neural network multi-criteria optimization image reconstruction technique (NN-MOIRT)15,16, combination of the TRP and ART17, modified LI18,19, total variation iterative soft thresholding (TV-IST)2,20 and others have been proposed in the literature. Of these algorithms, LBP, NNMOIRT, and TV-IST were used in previous BFB studies (Table 1). Cui et al.21 highlighted that the performance of the algorithms varied significantly

Page 3 of 38 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 4 of 38

with imaging tasks and hence, the selection of an algorithm must be taskspecific. LBP algorithm is the simplest and fastest4 and hence, it has been most commonly used in the previous studies (Table 1). However, it results in relatively low-resolution images that demonstrate a surface enclosing bubbles with a gradual transition of volume fraction.22,23 Owing to the blurring of the boundaries, the estimation of the bubble size critically relies on a cut-off value of solids volume fraction (also referred as a threshold) that distinguishes a bubble from the dense phase. Halow and Nicoletti22 discussed the limitation concerning the selection of the threshold value. The previous studies used a threshold value between 0.3 and 0.35 to post-process their captured data. For example, Hulme and Kantzas24 used a threshold of 0.3, whereas Asegehegn et al.25 used a threshold of 0.33 to calculate the bubble diameter in a 2D fluidized bed. Holland et al.26 determined that measurements from ECVT and MRI were consistent with each other at a threshold of 0.315. Weber and Mei7 also used 0.315 to calculate the bubble diameter; however, they also highlighted the unavailability of a method to determine the threshold. White27,28 performed experiments by maintaining a hollow sphere inside a static bed and repeated these experiments using spheres of different size. He also used LBP algorithm and constructed contours of solids volume fraction at three thresholds (0.2, 0.5, and 0.8). He could not estimate the size of a 6.5 cm diameter sphere with any of these values. He attributed this deficiency to poor spatial resolution of the ECT images and concluded that a variable threshold must be used to estimate the size of different spheres. The present study proposes an iterative method to determine an optimum threshold value for post-processing of each 2D image obtained by ECVT. For this, static fluidized bed experiments with a known size beaker

Page 4 of 38 ACS Paragon Plus Environment

Page 5 of 38 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

maintained at different positions of the bed are conducted. These experiments are repeated for different sizes of beakers and the captured ECVT data are used to develop an iterative method to calculate an optimum threshold value. Then, experiments with three fluidization velocities are conducted and the captured data are post-processed using the LBP algorithm with the developed method for optimum thresholds to estimate bubble sizes. Further, rise velocities, frequency, and evolution of bubbles at varying fluidization velocity are also estimated and analyzed.

2 2.1

Experiments Setup

Figure 1 displays a schematic of the experimental setup consisting of a cylindrical fluidized bed, 20 cm diameter (D) and 120 cm height (H). The fluidized bed column was attached to a 15 cm long plenum chamber, which was filled with 2 cm glass marbles for uniform gas distribution. Between the column and plenum chamber, a perforated distributor with 302 uniformly spaced holes of 1 mm was placed. Further, a nylon mesh of 250 µm was placed on top of the distributor to prevent the solids from sinking. Sand was used as the solid, which had a mean diameter of 320 µm with a size distribution of 100–600 µm measured by a Malvern Mastersizer (Figure 2(a)) and can be classified as Geldart B type. The average sphericity of the sand particles was 0.7 as determined by image analysis. The minimum fluidization velocity ( ) of sand was calculated by measuring the pressure drop across the bed at different air velocities. From the pressure drop profile (Figure 2(b)),  was determined to be 5.8 cm/s. Consequently, the experiments were performed with an initial bed height of 35 cm using three different fluidizing air velocities, 1.25, 1.85, and 2.5  . 2.2

Measurement technique

Page 5 of 38 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 6 of 38

The experimental setup with an ECVT sensor embedded 5 cm above the distributor is displayed in Figure 1(b). The sensor, developed and supplied by Tech4Imaging29, had a 20 cm diameter and 35 cm length. It consisted of four planes of electrodes with six electrodes on each plane (total 24 electrodes) arranged along the wall. The capacitance between an electrode and another electrode surrounding the bed was captured by the data acquisition system. Such measurements were performed by 24 electrodes. The captured capacitance data were translated to the distribution of solids volume fraction using an image reconstruction algorithm, which provided an image of the flow domain covered by the sensor at a resolution of a 1 cm × 1.75 cm × 1 cm (x × y × z) with 8000 voxels. For calibration of the ECVT sensor, capacitance was measured for an empty test section and then the section entirely filled with sand. These capacitance data were used to normalize the measurements between 0 and 1 with 0 representing an empty column and 1 representing a completely filled column having the solids volume fraction of 0.63. Before conducting further experiments, tests were conducted by measuring solids volume fraction in an empty, completely filled, and half filled column. As recommended by Makkawi and Wright30, solids volume fraction data were recorded for 60 s at 50 Hz frequency for each experiment. The instantaneous data was time-averaged to report steady-state values. All experiments were conducted three times and the time-averaged values with error bars representing maximum and minimum values were reported. 2.3

Selection of image reconstruction algorithm

The translation of capacitance data to solids volume fraction is a critical post-processing step that is performed using an image reconstruction algorithm. Algorithms such as LBP, NN-MOIRT, and TV-IST were used in previous studies (Table 1). In the present study, the effect of the image reconstruction algorithm was investigated by comparing the images from

Page 6 of 38 ACS Paragon Plus Environment

Page 7 of 38 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

the LBP, LI, and NN-MOIRT algorithms. For this, the bed was completely filled with sand to the upper end of the sensor and an empty glass beaker of known size was maintained at the center of bed as displayed in Figure 3(a). Then, data were recorded and images were constructed. As indicated in Figure 3(b), LBP resulted in clearer images of the object (beaker) in both horizontal and vertical planes compared to those from the other two algorithms. NN-MOIRT resulted in blurry images, whereas LI resulted in a sharp image in the horizontal plane with a distorted image in the vertical plane. In the LBP images, the transition from dense to lean phase was blurry; whereas an oval shape of the object was visible for the cylindrical beaker. It is clear from Figure 3(b) that a sharp transition between two phases was not captured by any of the algorithms. Of the three algorithms, LBP yielded sharper images that could be used with an optimized threshold to reconstruct the object. Hence, LBP was selected for post-processing of the data. 2.4

Calculation of bubble diameter

White27 conducted static fluidized bed experiments with spherical balls inside the bed and determined that a constant threshold could not capture all the spheres. He recommended a variable threshold for different size balls. To determine a relation between the threshold and size of the bubble, extensive calibration experiments were performed as follows. ECVT data were recorded for the static bed with an empty glass beaker as indicated in Figure 3(a). These experiments were repeated for four different glass beakers (3.5, 5.8, 8.8, and 11.6 cm diameter). It is noteworthy that the selected range of the beaker diameters represents a range of possible bubble diameters in the experimented BFB. Horizontal plane images (Figure 3(b)) at the center of the measurement volume were extracted for all the glass beakers; they were further post-processed using MATLAB to obtain the diameter of the beakers. The contour plot (Figure 3(b)) indicates a gradual

Page 7 of 38 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 8 of 38

transition of solids volume fractions at the bubble surface. MATLAB code was used to calculate the number of pixels having a solids fraction less than a given threshold in a given horizontal image captured at a particular axial location. The calculated number of pixels was multiplied by the pixel size of 1 × 1 cm2 and then the calculated value was converted to the beaker diameter. Table 2 presents the calculated beaker diameters using different thresholds. The use of a constant threshold resulted in a considerable error. For example, threshold of 0.189, 0.252, or 0.315 could not detect the presence of a 3.5 cm beaker, whereas threshold above 0.441 significantly overestimated the beaker diameter. Consequently, optimum threshold that minimize the error between the predicted and actual beaker diameters were determined by trial and error (Table 2). It is known that ECVT is a “soft-field” technique, where measurements depend on the location of an object in the sensor domain. Thus, the optimum threshold found for the beakers maintained at the center of the bed may not be applicable to bubbles at different locations. Consequently, additional calibration experiments were conducted by keeping the beakers at various axial and radial locations. For these experiments, the axial and radial locations were selected in such a way that the base of beakers was located at 10 and 20 cm from the base of the sensor and the sidewall of beakers were positioned at 5 and 10 cm from the wall of the sensor. Figure 4(a) shows positions of beakers. It should be noted that the selected locations were not near the wall or upper end or bottom part of the sensor. This is because the images of the beakers kept near the sensor ends were highly distorted and least accurate in terms of resolution owing to electric field fringing effect (Tech4imaging, private communication). Figure 4(b) shows optimum threshold for different beaker sizes kept at various locations. Depending on the size of the beaker, the variation in optimum threshold with location was within 10%.

Page 8 of 38 ACS Paragon Plus Environment

Page 9 of 38 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

Calibration was performed using average optimal threshold for a given beaker, whereas, the minimum and maximum of optimum threshold were used to establish error bounds on the calculated bubble diameters. An error of ±20% was observed for smaller bubbles (diameter ~ 3.5 cm). The error gradually decreased with increasing bubble diameter. For 11.6 cm bubble, the error was less than ±5%. Consequently, a polynomial function correlating the averaged optimum threshold and bubble diameter was derived.

 = 0.0054  − 0.0755 + 0.5829

(1)

Where,  is the averaged optimum threshold and  is the bubble diameter. This function was used in the iterative algorithm (Figure 5) to calculate bubble diameters. The convergence criteria of the algorithm had two conditions. Firstly, the difference between the bubble diameter calculated by image analysis (, ) applying a given threshold and that from the polynomial equation (, ) should be less than one pixel area. Secondly, the , must be less than or equal to the maximum possible diameter (, ) which is defined as the diameter of bubble formed by available excess gas in that particular horizontal plane ((єg - 0.37) × Ac/s). 2.5

Calculation of bubble velocity

Makkawi and Wright31 proposed a simple method to calculate bubble rise velocity from instantaneous solids volume fraction data. Following this method, time-series of cross-sectional averaged solids volume fractions at two axial heights (h1 and h2) were plotted (Figure 6(a)). Both these timeseries were observed oscillating. Assuming a single bubble evolution in the bed (discussed in Section 3.2), each large dip in the time series implied the detection of a bubble. When the two time-series were compared, a particular dip at a time t1 in the time-series at height h1 could be mapped to a dip at

Page 9 of 38 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 10 of 38

time t2 in the time-series at height h2 (with t2 > t1). For all mapped dips, this implied that the bubble detected at height h1 was again detected at height h2, however, after some time ∆t (= t2 − t1). Both (h2 − h1) and ∆t were used to calculate the rise velocity of the bubble. However, the time series also included fluctuations within a small number of larger dips as displayed in Figure 6(b). For these dips, the identification of a single minima and mapping of the dips was rather difficult. To overcome this challenge, the time series were further post-processed as described below. For each time series, the consecutive maximum of different peaks (as indicated in colored dots in Figure 6(b)) was determined. These maxima are written as

!"# $%  = &, ' , &, ' , … … &, ' , … , &, ' (

)

,

+

(2)

Where  is the solids volume fraction, &, ' is the maximum solids +

volume fraction of a peak at time - , the subscript . denotes a time interval number, and -/ to -0 are the time when , occurs. These maximum values were used to normalize the time series as follows:

, = 1 2 − &, '

- ≤ - < -5/

+

(3)

Where , is the normalized solids volume fraction and 1 2 is the solids volume fraction at any time -. Eq (3) implies that the signal is normalized by using successive maxima. Assume &, ' and &, ' (

)

are two consecutive maxima. Then, values between times -/ and - are normalized by using &, ' . The resulting normalized time series at h1 (

and h2 are displayed in Figure 6(c) with all values less than zero. The normalized time series included dips of different magnitudes with each having distinct minima, and therefore, the mapping of the dips was easier.

Page 10 of 38 ACS Paragon Plus Environment

Page 11 of 38 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

Rise velocity was calculated for all mapped dips in the two normalized time-series and the average of the calculated velocities was reported as bubble rise velocity at a given height (h1 + h2)/2. For data collected over 60 s, 90–100 mapped dips were obtained and it was determined that 50 mapped dips were sufficient to obtain statistically invariant data.

3 3.1

Results and discussion Distribution of solids volume fraction

Time-averaged radial distributions of solids volume fraction at three different axial locations (10.25, 20.75, and 31.25 cm above the distributor) are displayed in Figures 7(a), (b), and (c), respectively. At 10.25 and 20.75 cm, profiles at 1.25  were flat with values close to maximum packing of 0.63. This was mainly due to low excess gas available for bubble formation at such a low fluidization velocity. At higher velocities (1.85 and 2.5  ), higher excess gas caused bubble formation leading to lower solids volume fraction and a core-annulus structure with high solids volume fraction near the wall and low values at the center of the bed. Escudero and Heindel32 also observed a core-annulus structure at higher fluidization velocity at a height above H/D of one. They observed that rising gas bubbles coalesced and migrated toward the center of the bed, causing high gas holdup in the center region. They also observed that larger bubbles erupt at the top of the bed and throw solids towards the wall. Both the coalescence and eruption of bubbles causes a core-annulus profile in fluidized beds. Solids volume fraction decreased with an increase in gas velocity at all three velocities (Figure 7). However, the decrease near the wall was minor compare to that at the center of the bed. This implies that higher fluidization velocity does not promote mixing near the wall; rather it causes channeling of gas in the center region. The variation of solids volume fraction along the height is indicated in Figure 7(d). It remained virtually constant; however, it

Page 11 of 38 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 12 of 38

decreased with an increase in fluidization velocity. After a height of approximately 27–30 cm, solids volume fraction declined, with the decline being rather steep for 1.25  compared to the other two velocities. This decline can be attributed to the eruption of bubbles and expansion of the bed near the top. At a low velocity of 1.25  , bed expansion was less than the sensor height and hence, the solids volume fraction began declining at approximately 27 cm height and approached zero at 40 cm height. At 2.5  , bed expansion covered the entire sensor height (up to 40 cm) and thus, the decline in solids volume fraction began at approximately 32 cm height; however, the bed had more than 40% solids at 40 cm height. The axial profile of solids volume fraction at 1.25  was also compared with the data of Escudero and Heinden32, who conducted BFB experiments of Geldart B particles using X-ray. A good qualitative agreement was found for the solids volume fraction profile (Figure 7(e)). The discrepancy in values can be attributed to the differences in the column diameter, initial bed height, and particle properties. 3.2

Bubble evolution

The time series of solids volume fraction at 22.5 cm height for all three fluidization velocities are displayed in Figure 8(a). The time series indicated fluctuations around a certain average value that decreased with an increase in the velocity. The scale of fluctuations increased with an increase in the air velocity. This suggests that the bed became increasingly chaotic with a rise in air velocity. For 2.5  , fluctuations in the time series with a small number of peaks between 30–32 s are indicated in Figure 8(b). Each dip represents a bubble. Corresponding to the first dip in Figure 8(b), contours of solids volume fraction in a 2D slice at 22.5 cm height are displayed in Figure 8(c). These contours demonstrate evolution of a single bubble passing through the 22.5 cm bed height. The bubble that formed between

Page 12 of 38 ACS Paragon Plus Environment

Page 13 of 38 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

30.2 and 30.32 s was attached to the wall and then increased in size between 30.32 and 30.44 s. It then detached from the wall between 30.48 and 30.50 s, and began growing opposite to its previous position. The detachment of the bubble may have caused the minor dip beginning at approximately 30.5 s in the time series. After the detachment, the bubble continued growing between 30.50 and 30.60 s. It then gradually disappeared from the 2D frame during last few fractions of a second. The contours clearly suggest a single bubble passing through a given height with both bubble size and position varying with time. The evolution of 3D bubbles corresponding to 2D contours of Figure 8 is shown in Figure 9, where 2D contours in a vertical plane and corresponding 3D bubbles obtained by applying a constant (0.315) and optimum thresholds are presented for a period of 0.30 s, starting from 30.2 s. It can be seen that both the 2D contours and reconstructed bubbles captured the complex flow behavior within the fluidized bed. However, the images reconstructed using the constant threshold missed a large number of bubbles as highlighted in Figure 9(a), (e), (f), (g), and (h). The reconstructed bubbles using the constant threshold were also much smaller compared to voids seen in the 2D contours. On the other hand, the reconstructed bubbles using optimal threshold were consistent with the 2D contours. Using optimum threshold, several additional bubbles were identified. However, optimal threshold was not able to identify bubbles that were close to each other. For example, in Figures 9(g), one can visually distinguish three bubbles but the optimum threshold method resulted in a single elongated bubble. Such discrepancy was also noticed when two bubbles were very close to each other at the same axial plane. When two bubbles were sufficiently apart (such as in Figure 9(i)), the optimal threshold method could identify two distinct bubbles. It is possible to further increase the number of bubbles identified by using high-resolution ECVT sensor and further tuning of the algorithm.

Page 13 of 38 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

3.3

Page 14 of 38

Bubble frequency

Makkawi and Wright30 used time series of average solids volume fraction at a particular height to calculate the bubble cycle frequency, which is defined as half of the total number of times a signal passes over its mean (6 ) in a given time interval. Using this method, bubble frequency 1% 2 is equal to 6 ⁄ 27 , where 7 represents the total sampling period in seconds. The variation in bubble frequency with distance from the distributor is displayed in Figure 10(a). Bubble frequency decreased with bed height and increased with superficial air velocity. Many researchers30,33–35 analyzed bubble frequency by plotting the power spectra density (PSD) of solids volume fraction time series. In PSD analysis, the dominant frequency is considered as the bubble frequency in the case of BFB and amplitude can be linked to bubble size. For a single bubble of uniform size, a sharp single peak appears in the PSD plots.35 In the case of BFB, a wide band of dominant frequencies is observed; this signifies that different sizes of bubbles appear either at the same locations or at different times. Figures 10(b), (c), and (d) display the PSDs of solids volume fraction time series at 22.5 cm height for the three air velocities. The amplitude in PSDs increases with an increase in air velocity. This was consistent with higher bubble size resulting from higher air velocity. The PSDs indicated no clear peak or single dominant frequency; however, they had several peaks appearing over the band of 1 to 4 Hz. This implies that different sizes of bubbles pass through a given bed height as indicated in Figure 8(c). 3.4

Bubble diameter

Bubble diameters at axial locations between 12 – 27 cm bed heights, where the calibration experiments were conducted, were calculated for superficial air velocities of 1.85 and 2.5Umf by applying the iterative method described in Section 2.4. Figure 11(a) displays a variation of bubble diameter with

Page 14 of 38 ACS Paragon Plus Environment

Page 15 of 38 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

axial distance from the distributor. At both 1.85 and 2.5Umf, bubble size increased between ~12 to 20 cm bed heights and then decreased. The increase in bubble size can be attributed to coalescence of smaller bubbles formed at the bottom section, whereas the decrease in bubble size can be attributed to splitting of bubbles. A comparison of the calculated bubble diameters from the ECVT data with those calculated from different empirical correlations36–41 is displayed in Figure 11(b). The growth of the bubble diameter between 12 and 20 cm bed height reasonably agreed with that calculated by the correlation of Kato and Wen40. The estimated bubble diameters were closer to the correlations of Choi et al.37, Horio and Nonaka39, and Mori and Wen41, whereas the correlation of Geldart38 significantly underestimated the estimated values. In 20–27 cm bed height, a declining trend of bubble diameter from the ECVT data did not match with monotonically increasing profile of the empirical correlations. The calculated bubble diameters of the present study were also compared with that reported by Weber and Mei7 (Figure 11(c)). Up to the bed height of 1.25 h/D, both the data sets showed increasing bubble diameter with values comparable to each other. After h/D of 1.25, the current data showed decrease in bubble diameter whereas those of Weber and Mei7 showed a steady increase. The declining trend of bubble diameter after the bed height of 20 cm can be attributed to splitting of bubbles as captured in Figure 9. Shen et al.42 also observed a sharp decline in bubble diameter at a top section of the bed. He defined a critical height h* after which bubbles became unstable due to its size and then split into more than one. However, it should be noted that the current data suggest the bubble splitting at the middle height of the bed (~20 cm) whereas Shen et al.42 reported such phenomena at extremely top part of the bed. Werther36 reported a stable bubble diameter after a certain height,

Page 15 of 38 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 16 of 38

where the equilibrium between coalescence of smaller bubbles and splitting of larger bubbles established. However, the current data did not show a stable bubble diameter. Besides splitting of bubble at the middle of the bed, another possible reason for this discrepancy might be the wide particle size distribution used in the present study. Beestra et al.43, Rautenbach et al.44,45 and Brouwer et al.46 have reported decrease in bubble diameter with presence of fines in the bed. Furthermore, bi-dispersed fluidized bed experiments47–49 suggested that fine particles segregated at the top of the bed at low fluidization velocity. Such segregation can also affect bubble properties in the different region of the bed. In the present study, particles size varied from 100 to 600 µm with 6% particles less than 100 µm. 3.5

Bubble rise velocity

Bubble rise velocities at three axial locations (10.25, 17.25, and 24.25 cm from the distributor) were calculated using the method described in Section 2.5. Figure 12(a) indicates that bubble rise velocity increased with an increase in both air velocity and bed height. This can be attributed to the growth in bubble size between 10–25 cm height (Figure 11), resulting in a higher buoyancy of the bubbles. At 1.25  , the increase in the rise velocity was observed between 10.25 and 17.25 cm heights; the increase was muted between 17.25 and 24.25 cm. For both 1.85 and 2.5  , the bubble rise velocity increased steadily from 10.25 to 24.25 cm height. Different empirical correlations are available in the literature to predict bubble rise velocity in gas-solid fluidized beds.50–52 The bubble rise velocities estimated from the ECVT data were compared with those calculated using empirical correlations (Figure 12(b)). The trend of the estimated values agreed with that calculated using the correlations.36,50–53 Quantitatively, the estimated values were closer to those from the correlation of Werther36; the correlation of Dry et al.51 resulted in significantly higher values. The discrepancy between the experimental

Page 16 of 38 ACS Paragon Plus Environment

Page 17 of 38 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

values and those from the correlation of Dry et al.51 was also observed at the other two air velocities; whereas experimental values were consistently in agreement with those from the correlation of Werther36. Previous studies4,54 also used the cross-correlation (CC) method to calculate the bubble rise velocity. In this method, the CC function (91:22 (eq. 4) from the solids volume fraction time series captured at two heights (h1 and h2) was calculated as /

9 1:2 =  ∑?/ # 1.∆2 × >1.∆ + :∆2 : = 1,2,3, … , !

(4)

Where, # and > are the solids volume fraction time series at h1 and h2 respectively, Δ is the sampling interval, B is the number of sample intervals,

! is the number of samples in the CC function and : is the time delay. The CC function will give the time interval 1C2 when the similarity index between the two time series is maximum. The time lag (D = C∆2 is then used to calculate the bubble rise velocity (ℎ − ℎ/ ⁄D 2. In the current study, rise velocity calculated from CC and signal analysis method described in Section 2.5 were compared in Figure 12. It should be noted that the CC function values for 1.25Umf is not reported in Figure 12, as no time lag could be determined for two time series. For 1.85Umf and 2.5Umf, the CC method resulted in significantly higher values than those calculated by signal analysis. Further, the values from the CC method did not show an increasing bubble velocity as resulted by the empirical correlations. In the CC method, the lag was found at the highest similarity between two time series, whereas the signal analysis method uses the average lag calculated for all mapped dips in two time series. The CC method is influenced by noise and frequency of measurement, whereas the signal analysis method filters out the noise before mapping of dips. Thus, the values from the signal

Page 17 of 38 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 18 of 38

analysis method is more reliable and resulted in values that were qualitatively consistent with the empirical correlations.

4

Conclusion

BFB experiments (sand-air, D = 20 cm, and H/D = 1.75) were performed at three fluidization velocities (1.25, 1.85, and 2.5  ). ECVT was used to measure the distribution of solids volume fraction, which was further analyzed to estimate diameter, rise velocity, and frequency of bubbles. The previous studies used a constant threshold value between 0.3 and 0.35 to post-process their data and estimate the size of bubbles. The present study conducted extensive calibration experiments using a static fluidized bed with a known size beaker. It was determined that a low threshold value between 0.189 and 0.315 could not detect the presence of small bubbles (< 3.5 cm). The data from the calibration experiments were used to develop an iterative method to determine the optimum threshold for each 2D image. The profiles of the bubble diameter along the bed height indicated bubble growth up to the middle height (~20 cm) of the bed, after which bubble diameter decreased. The growth of the bubbles was consistent with that calculated by the empirical correlation of Kato and Wen40. The growth of bubbles also caused an increase in the bubble rise velocity. The estimated bubble rise velocities reasonably agreed with the correlation of Werther36. The estimated bubble frequency was in a band of 2–4 Hz at three fluidization velocities. The contours of solids volume fraction and 3D bubbles calculated by optimum threshold suggested that the size and shape of the bubbles continuously varied owing to the coalescence, elongation, attachment and detachment with the wall, breakup, and eruption of the bubbles.

5

Acknowledgement

Page 18 of 38 ACS Paragon Plus Environment

Page 19 of 38 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

The authors would like to acknowledge the Department of Industry, Innovation and Science – Government of Australia for funding the research under Australia-India Strategic Research Fund (AISRF) Grant No. GCF020010.

6

Notations

F  D fb h H j k m n Nc R t T  x, y

Area (m2) Bubble diameter (m) Column diameter (m) Bubble frequency (s-1) Axial height above distributor (m) Column height (m) Time interval (-) Time interval at which cross-correlation is maximum (-) Number of samples in cross-correlation function (-) Number of sample interval (-) Number of times a signal passes over its mean (-) Cross – correlation function (-) Time (s) Total sampling period (s) Fluidization velocity (m/s) Solids volume fraction time series at height h1 and h2 (-)

Greek letters Δ  D

Sampling interval (s) Volume fraction (-) Time lag (s)

Subscript i b calc calibration c/s g max mf norm s t

Time interval number Bubble Calculated from image analysis Calculated from equation (1) Cross section Gas Maximum Minimum fluidization condition Normalized Solids Time

Page 19 of 38 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

Abbreviation 2D 3D ART BFB CC ECT ECVT LBP LI MRI NN-MOIRT

PEPT PSD RPT SIRT TRP TV-IST

7

(1)

Page 20 of 38

Two dimensional Three dimensional Algebraic reconstruction technique Bubbling fluidized bed Cross- correlation Electrical capacitance tomography Electrical capacitance volume tomography Linear back projection Landweber iteration Magnetic resonance imaging Neural network multi-criteria optimization reconstruction technique Positron emission particle tracking Power spectral density Radioactive particle tracking Simultaneous iterative reconstruction technique Tikhonov regularization principle Total variation iterative soft thresholding

image

References

Sun, J.; Yan, Y. Non-Intrusive Measurement and Hydrodynamics Characterization of Gas–solid Fluidized Beds: A Review. Meas. Sci. Technol. 2016, 27 (11), 112001.

(2)

Chandrasekera, T. C.; Li, Y.; Moody, D.; Schnellmann, M. A.; Dennis, J. S.; Holland, D. J. Measurement of Bubble Sizes in Fluidised Beds Using Electrical Capacitance Tomography. Chem. Eng. Sci. 2015, 126, 679.

(3)

Du, B.; Warsito, W.; Fan, L.-S. ECT Studies of Gas−Solid Fluidized Beds of Different Diameters. Ind. Eng. Chem. Res. 2005, 44 (14), 5020.

(4)

Makkawi, Y. T.; Wright, P. C. Fluidization Regimes in a Conventional Fluidized Bed Characterized by Means of Electrical Capacitance Tomography. Chem. Eng. Sci. 2002,

57 (13), 2411. (5)

McKeen, T.; Pugsley, T. Simulation and Experimental Validation of a Freely Bubbling Bed of FCC Catalyst. Powder Technol. 2003, 129 (1–3), 139.

Page 20 of 38 ACS Paragon Plus Environment

Page 21 of 38 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

(6)

Rautenbach, C.; Melaaen, M. C.; Halvorsen, B. M. Statistical Diagnosis of a Gas–solid Fluidized Bed Using Electrical Capacitance Tomography. Int. J. Multiph. Flow 2013,

49, 70. (7)

Weber, J. M.; Mei, J. S. Bubbling Fluidized Bed Characterization Using Electrical Capacitance Volume Tomography (ECVT). Powder Technol. 2013, 242, 40.

(8)

Guo, Q.; Meng, S.; Wang, D.; Zhao, Y.; Ye, M.; Yang, W.; Liu, Z. Investigation of Gas-Solid Bubbling Fluidized Beds Using ECT with a Modified Tikhonov Regularization Technique. AIChE J. 2017.

(9)

Xie, C. G.; Huang, S. M.; Beck, M. S.; Hoyle, B. S.; Thorn, R.; Lenn, C.; Snowden, D. Electrical Capacitance Tomography for Flow Imaging: System Model for Development of Image Reconstruction Algorithms and Design of Primary Sensors.

IEE Proc. G Circuits, Devices Syst. 1992, 139 (1), 89. (10)

Yang, W. Q.; Spink, D. M.; York, T. A.; McCann, H. An Image-Reconstruction Algorithm Based on Landweber’s Iteration Method for Electrical-Capacitance Tomography. Meas. Sci. Technol. 1999, 10 (11), 1065.

(11)

Tikhonov, A. N.; Goncharsky, A. V.; Stepanov, V. V.; Yagola, A. G. Numerical

Methods for the Solution of Ill-Posed Problems; Springer Netherlands: Dordrecht, 1995. (12)

Peng, L.; Merkus, H.; Scarlett, B. Using Regularization Methods for Image Reconstruction of Electrical Capacitance Tomography. Part. Part. Syst. Charact. 2000, 17 (3), 96.

(13)

Yang, W. Q.; Peng, L. Image Reconstruction Algorithms for Electrical Capacitance Tomography. Meas. Sci. Technol. 2003, 14 (1), R1.

(14)

Bangliang, S.; Yiheng, Z.; Lihui, P.; Danya, Y.; Baofen, Z. The Use of Simultaneous Iterative Reconstruction Technique for Electrical Capacitance Tomography. Chem.

Eng. J. 2000, 77 (1–2), 37. (15)

Marashdeh, Q.; Warsito, W.; Fan, L.-S.; Teixeira, F. L. A Nonlinear Image Reconstruction Technique for ECT Using a Combined Neural Network Approach. Page 21 of 38 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 22 of 38

Meas. Sci. Technol. 2006, 17 (8), 2097. (16)

Warsito, W.; Fan, L.-S. Neural Network Based Multi-Criterion Optimization Image Reconstruction Technique for Imaging Two- and Three-Phase Flow Systems Using Electrical Capacitance Tomography. Meas. Sci. Technol. 2001, 12 (12), 2198.

(17)

Huang, Z.; Wang, B.; Li, H. Dynamic Voidage Measurements in a Gas Solid Fluidized Bed by Electrical Capacitance Tomography. Chem. Eng. Commun. 2003, 190 (10), 1395.

(18)

Yan, H.; Wang, Y. F.; Zhou, Y. G.; Sun, Y. H. 3D ECT Reconstruction by an Improved Landweber Iteration Algorithm. Flow Meas. Instrum. 2014, 37, 92.

(19)

Ye, J.; Wang, H.; Yang, W. A Sparsity Reconstruction Algorithm for Electrical Capacitance Tomography Based on Modified Landweber Iteration. Meas. Sci.

Technol. 2014, 25 (11), 115402. (20)

Chandrasekera, T. C.; Li, Y.; Dennis, J. S.; Holland, D. J. Total Variation Image Reconstruction for Electrical Capacitance Tomography. In 2012 IEEE International

Conference on Imaging Systems and Techniques Proceedings; IEEE, 2012; Vol. 50, pp 584–589. (21)

Cui, Z.; Wang, Q.; Xue, Q.; Fan, W.; Zhang, L.; Cao, Z.; Sun, B.; Wang, H.; Yang, W. A Review on Image Reconstruction Algorithms for Electrical Capacitance/resistance Tomography. Sens. Rev. 2016, 36 (4), 429.

(22)

Halow, J. S.; Nicoletti, P. Observations of Fluidized Bed Coalescence Using Capacitance Imaging. Powder Technol. 1992, 69 (3), 255.

(23)

Halow, J. S.; Fasching, G. E.; Nicoletti, P.; Spenik, J. L. Observations of a Fluidized Bed Using Capacitance Imaging. Chem. Eng. Sci. 1993, 48 (4), 643.

(24)

Hulme, I.; Kantzas, A. Determination of Bubble Diameter and Axial Velocity for a Polyethylene Fluidized Bed Using X-Ray Fluoroscopy. Powder Technol. 2004, 147 (1–3), 20.

(25)

Asegehegn, T. W.; Schreiber, M.; Krautz, H. J. Investigation of Bubble Behavior in Fluidized Beds with and without Immersed Horizontal Tubes Using a Digital Image Page 22 of 38 ACS Paragon Plus Environment

Page 23 of 38 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

Analysis Technique. Powder Technol. 2011, 210 (3), 248. (26)

Holland, D. J.; Marashdeh, Q.; Müller, C. R.; Wang, F.; Dennis, J. S.; Fan, L.-S.; Gladden, L. F. Comparison of ECVT and MR Measurements of Voidage in a GasFluidized Bed. Ind. Eng. Chem. Res. 2009, 48 (1), 172.

(27)

White, R. B. Using Electrical Capacitance Tomography to Monitor Gas Voids in a Packed Bed of Solids. Meas. Sci. Technol. 2002, 13 (12), 1842.

(28)

White, R. B. Using Electrical Capacitance Tomography to Investigate Gas Solid Contacting. Can. J. Chem. Eng. 2005, 83 (1), 64.

(29)

Warsito, W.; Marashdeh, Q.; Fan, L. S. 3D and Real Time Electrical Capacitance Volume-Tomography Sensor Design and Image Reconstruction. US 8614707 B2, 2013.

(30)

Makkawi, Y. T.; Wright, P. C. Optimization of Experiment Span and Data Acquisition Rate for Reliable Electrical Capacitance Tomography Measurement in Fluidization Studies a Case Study. Meas. Sci. Technol. 2002, 13 (12), 1831.

(31)

Makkawi, Y. T.; Wright, P. C. Electrical Capacitance Tomography for Conventional Fluidized Bed Measurements—remarks on the Measuring Technique. Powder

Technol. 2004, 148 (2–3), 142. (32)

Escudero, D.; Heindel, T. J. Bed Height and Material Density Effects on Fluidized Bed Hydrodynamics. Chem. Eng. Sci. 2011, 66 (16), 3648.

(33)

Acosta-Iborra, A.; Sobrino, C.; Hernández-Jiménez, F.; de Vega, M. Experimental and Computational Study on the Bubble Behavior in a 3-D Fluidized Bed. Chem. Eng. Sci. 2011, 66 (15), 3499.

(34)

Johnsson, F.; Zijerveld, R. .; Schouten, J. .; van den Bleek, C. .; Leckner, B. Characterization of Fluidization Regimes by Time-Series Analysis of Pressure Fluctuations. Int. J. Multiph. Flow 2000, 26 (4), 663.

(35)

Svensson, A.; Johnsson, F.; Leckner, B. Fluidization Regimes in Non-Slugging Fluidized Beds: The Influence of Pressure Drop across the Air Distributor. Powder

Technol. 1996, 86 (3), 299. Page 23 of 38 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

(36)

Page 24 of 38

Werther, J. Effect of Gas Distributor on the Hydrodynamics of Gas Fluidized Beds.

Ger. Chem. Eng. 1978, 1, 166. (37)

Choi, J.; Son, J.; Kim, S. Generalized Model for Bubble Size and Frequency in GasFluidized Beds. Ind. Eng. Chem. Res. 1998, 37 (6), 2559.

(38)

Geldart, D. The Effect of Particle Size and Size Distribution on the Behaviour of GasFluidised Beds. Powder Technol. 1972, 6 (4), 201.

(39)

Horio, M.; Nonaka, A. A Generalized Bubble Diameter Correlation for Gas-Solid Fluidized Beds. AIChE J. 1987, 33 (11), 1865.

(40)

Kato, K.; Wen, C. Y. Bubble Assemblage Model for Fluidized Bed Catalytic Reactors.

Chem. Eng. Sci. 1969, 24 (8), 1351. (41)

Mori, S.; Wen, C. Y. Estimation of Bubble Diameter in Gaseous Fluidized Beds.

AIChE J. 1975, 21 (1), 109. (42)

Shen, L.; Johnsson, F.; Leckner, B. Digital Image Analysis of Hydrodynamics TwoDimensional Bubbling Fluidized Beds. Chem. Eng. Sci. 2004, 59 (13), 2607.

(43)

Beetstra, R.; Nijenhuis, J.; Ellis, N.; van Ommen, J. R. The Influence of the Particle Size Distribution on Fluidized Bed Hydrodynamics Using High-Throughput Experimentation. AIChE J. 2009, 55 (8), 2013.

(44)

Rautenbach, C.; Mudde, R. F.; Yang, X.; Melaaen, M. C.; Halvorsen, B. M. A Comparative Study between Electrical Capacitance Tomography and Time-Resolved X-Ray Tomography. Flow Meas. Instrum. 2013, 30, 34.

(45)

Rautenbach, C.; Melaaen, M. C.; Halvorsen, B. M. Investigating the Influence of Fines in Fluidized Bed Reactors Using 3D ECT Images. In WIT Transactions on

Engineering Sciences,; 2011; Vol. 70, pp 141–151. (46)

Brouwer, G. C.; Wagner, E. C.; van Ommen, J. R.; Mudde, R. F. Effects of Pressure and Fines Content on Bubble Diameter in a Fluidized Bed Studied Using Fast X-Ray Tomography. Chem. Eng. J. 2012, 207–208, 711.

(47)

Huilin, L.; Yurong, H.; Gidaspow, D.; Lidan, Y.; Yukun, Q. Size Segregation of Page 24 of 38 ACS Paragon Plus Environment

Page 25 of 38 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

Binary Mixture of Solids in Bubbling Fluidized Beds. Powder Technol. 2003, 134 (1– 2), 86. (48)

Marzocchella, A.; Salatino, P.; Di Pastena, V.; Lirer, L. Transient Fluidization and Segregation of Binary Mixtures of Particles. AIChE J. 2000, 46 (11), 2175.

(49)

Formisani, B.; Cristofaro, G. D.; Girimonte, R. A Fundamental Approach to the Phenomenology of Fluidization of Size Segregating Binary Mixtures of Solids. Chem.

Eng. Sci. 2001, 56 (1), 109. (50)

Davidson, J. F.; Harrison, D.; Jackson, R. Fluidized Particles: Cambridge University Press, 1963. 155 Pp. 35s. Pergamon 1964.

(51)

Dry, R. J.; Judd, M. R.; Shingles, T. Bubble Velocities in Fluidized Beds of Fine, Dense Powders. Powder Technol. 1984, 39 (1), 69.

(52)

Hillegardt, K.; Werther, J. Local Bubble Gas Hold-up and Expansion of Gas/solid Fluidized Beds. Ger. Chem. Eng. 1986, 9 (4), 215.

(53)

Kunii, D.; Levenspiel, O. Fluidization Engineering; Elsevier, 2013.

(54)

Liu, S.; Chen, Q.; Wang, H. G.; Jiang, F.; Ismail, I.; Yang, W. Q. Electrical Capacitance Tomography for Gas–solids Flow Measurement for Circulating Fluidized Beds. Flow Meas. Instrum. 2005, 16 (2–3), 135.

(55)

Yang, W. Q.; Liu, S. Role of Tomography in Gas/solids Flow Measurement. Flow

Meas. Instrum. 2000, 11 (3), 237. (56)

Chandrasekera, T. C.; Wang, A.; Holland, D. J.; Marashdeh, Q.; Pore, M.; Wang, F.; Sederman, A. J.; Fan, L. S.; Gladden, L. F.; Dennis, J. S. A Comparison of Magnetic Resonance Imaging and Electrical Capacitance Tomography: An Air Jet through a Bed of Particles. Powder Technol. 2012, 227, 86.

(57)

Ye, J.; Wang, H.; Yang, W. Characterization of a Multi-Plane Electrical Capacitance Tomography Sensor with Different Numbers of Electrodes. Meas. Sci. Technol. 2016,

27 (3), 35103.

8

Figures Page 25 of 38 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 26 of 38

(a)

(b) Figure 1: (a) Schematic diagram of the experimental setup and (b) fluidized bed with ECVT sensor and data acquisition system

Page 26 of 38 ACS Paragon Plus Environment

Page 27 of 38 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

(a)

(b)

Figure 2: (a) Particle size distribution and (b) variation of pressured drop across bed with fluidization velocity

LBP

NN-MOIRT

Landweber

(b) Horizontal plane

(a)

LBP

NN-MOIRT

Landweber

(c) Vertical plane

Figure 3: (a) Schematic diagram of stationary bed with empty beaker (Diameter – 5.8 cm), and (b) and (c) comparison of solids volume fraction contours from three image reconstruction algorithms

Page 27 of 38 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)

Page 28 of 38

(b)

Figure 4: (a) Positions of beakers in the sensor domain, and (b) optimum thresholds for different beakers kept at variation locations

Page 28 of 38 ACS Paragon Plus Environment

Page 29 of 38 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

Figure 5: Algorithm to calculate bubble diameter

Page 29 of 38 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 30 of 38

(a)

(b)

(c)

Figure 6: (a) Variation of solids volume fraction at two heights, (b) fluctuations within a dip, and (c) normalized time series

Page 30 of 38 ACS Paragon Plus Environment

Page 31 of 38 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

(a)

(b)

(c)

(d)

(e) Figure 7: Radial profiles of solids volume fraction at (a) 10.25 cm, (b) 20.75, (c) 31.25 cm, (d) axial profile of solids volume fraction (e) comparison of axial profile of solids volume fraction with literature data Page 31 of 38 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)

Page 32 of 38

(b)

30.2 s

30.24 s

30.28 s

30.32 s

30.36 s

30.40 s

30.44 s

30.48 s

30.50 s

30.52 s

30.56 s

30.60 s

30.64 s

30.68 s

30.72 s

30.76 s

(c) Figure 8: (a) Time series of solids volume fraction at 22.5 cm height, (b) fluctuations in the time series between 30 and 32 s, and (c) evolution of a bubble in a cross section at 22.5 cm height.

Page 32 of 38 ACS Paragon Plus Environment

Page 33 of 38 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

(a) 30.2 s

(b) 30.24 s

(c) 30.28 s

(d) 30.32 s

(e) 30.36 s

(f) 30.40 s

(g) 30.44 s

(h) 30.48 s Figure 9: Contour of solid volume fraction and calculated 3D bubbles from constant (0.315) and optimum threshold

(i) 30.50 s

Page 33 of 38 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 34 of 38

(a)

(b)

(c)

(d)

Figure 10: (a) Variation of bubble frequency along bed height; (b), (c), and (d) power spectra density plots of solids volume fraction time series at 1.25, 1.85, and 2.5Umf, respectively

Page 34 of 38 ACS Paragon Plus Environment

Page 35 of 38 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

(a)

(b)

(c) Figure 11: (a) Variation of bubble diameter along bed height (b) comparison between calculated bubble diameters from ECVT data and those calculated from empirical correlations, and (c) comparison of axial profile of bubble diameter with literature data.

Page 35 of 38 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)

Page 36 of 38

(b)

Figure 12: (a) Variation of bubble rise velocity along bed height and (b) comparison between calculated bubble velocity from ECVT data and those calculated from empirical correlations

Page 36 of 38 ACS Paragon Plus Environment

Page 37 of 38 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

9

Tables

Table-1: Previous BFB studies using either ECT or ECVT

Halow et al.23

Sensor details Dia. (cm)

Freq. (Hz)

Algorithm

Particle Dia. (µm)

15.24

62.5

-

704, 300, 4475, 3075, 3175

1.14

LBP

--

--

Yang and Liu55

Experiment details Solid density (gm/cm3)

Estimated properties Geldart classification B and D

Contribution

Umf (cm/s)

Gas velocity (U/Umf)

19.0, 8.0, 82.2, 80.1, 84.1

1-4

(b), (c)

Empirical correlations of bubble size and velocity

--

--

(e)

Fluidization regimes were identified and a method to calculate velocity of solids was proposed.

Makkawi and Wright4

13.8

100

LBP

530

2.6

B

37

0.54-6.4

(e)

Fluidization regimes were identified based on the calculated transition velocity.

Makkawi and Wright30

13.8

20-100

LBP

530

2.6

B

27

3.33

(a), (c), (d)

The effect of measurement frequency was investigated. The authors recommended at least 4000 data points to achieve statistically invariant estimations of bubble velocity and frequency.

White27

20

LBP

113, 153

3.333, 4.645

A, B

(b)

Need of a variable threshold to capture various sphere diameters was emphasized.

McKeen and Pugsley5

14

100

LBP

79

1.4

A

U = 5, 10, 15 and 20

(b), (c)

Data from ECT experiments were compared with CFD simulation results.

Makkawi and Wright31

13.8

100

LBP

90, 350, 530, 1850

2.5

B/A, B and D

-

(U-Ub)~0.2 to 0.6

(c), (d)

The effect of measurement span was studied and minimum 60 s measurement was recommended to get reproducible results.

Du et al.3

5,10, 30

100

NNMOIRT

60

1.4

A

-

(a), (b)

Data from ECVT were compared with those from an optical probe. Flow behavior in void phase and emulsion phase was characterized based on the two-region model.

Holland et al.26

5

80

NNMOIRT

58

0.6

A

0.014

2.8-7.0

(a)

Images of solids volume fraction distribution obtained from ECVT were compared with those from MRI.

Chandrasekera et al.56

5

80

NNMOIRT

1200

0.96

D

30