Scaling up stirred tanks is a significant challenge because of the research gaps between laboratory and industrial-scale setups. It is necessary to understand the effects of scale-up on the mass transfer in stirred tanks, and this requires meticulous experimental analysis. The present study investigates the effects of tank size and aspect ratio ( ) on the volumetric mass transfer coefficients of shear-thinning fluids. The experiments were conducted in three stirred tanks of different sizes (laboratory and pilot scale) and geometries (standard and nonstandard). was 1 for the standard tanks and 3.5 for the nonstandard stirred tanks. Three sizes of stirred tanks were used: 11 L with of 1, 40 L with of 3.5, and 47 L with of 1. Impeller stirring speeds and gas flow rates were in the range of 800–900 rev min−1 and 8–10 L min−1, respectively. The volumetric mass transfer coefficient was estimated based on the dissolved oxygen concentration in the fluids, and the effects of rheology and operating conditions on the volumetric mass transfer coefficient were observed. The volumetric mass transfer coefficient decreased as tank size increased and increased with an increase in operating conditions, but these effects were also clearly influenced by fluid rheology. The impacts of scale-up and operating conditions on the volumetric mass transfer coefficient decreased as liquid viscosity increased.
Scaling up stirred tanks in the presence of non-Newtonian fluids is highly challenging because of difficulties in translating the mass transfer process from the laboratory to the industrial scale. Fluid rheology, geometric parameters, and gas-distribution system are among the factors that control gas–liquid mass transfer rate in a stirred tank. Volumetric mass transfer coefficient (k L a) is a standard parameter to characterize mass transport by correlating the mass transfer rate to concentration change (García-Ochoa and Gómez 1998; Garcia-Ochoa and Gomez 2004). k L a must be known to design and scale-up stirred tanks. Dynamic (Bandyopadhyay, Humphrey, and Taguchi 1967) and steady-state (Imai, Takei, and Matsumura 1987) methods are generally used to measure k L a in stirred tanks. Dynamic methods (Bandyopadhyay, Humphrey, and Taguchi 1967) begin with deoxygenating the given liquid by using nitrogen. The cleaned air is then passed through the liquid, and a change in the dissolved oxygen (DO) concentration is recorded. k L a is estimated by fitting the transient curve of the DO concentration to an unsteady-state mass balance model. Dynamic methods (Bandyopadhyay, Humphrey, and Taguchi 1967) assume that there are well-mixed dispersed gas and liquid phases in the stirred tank, which is valid for small-scale reactors but not for large vessels, where achieving well-mixed gas and liquid phases is not always possible. Thus, large deviations in k L a measurements are expected when using dynamic methods (Bandyopadhyay, Humphrey, and Taguchi 1967) in large-scale reactors (Dunn and Einsele 2007). Additionally, the time response ( ) of O2 sensors poses a problem with dynamic methods (Bandyopadhyay, Humphrey, and Taguchi 1967) when the is higher than the characteristic time of mass transfer (1/k L a). Another concern about applying the dynamic method to viscous non-Newtonian fluids is the appearance of tiny gas bubbles suspended in such liquids. The contribution of these microbubbles to k L a is significantly smaller than that of large bubbles owing to low driving force (Muller and Davidson 1992; Shetty, Kantak, and Kelkar 1992). Thus, microbubbles are considered perfectly mixed, and their effect on k L a is mostly neglected (Gabelle et al. 2011; Heijnen, Riet, and Wolthuis 1980). Steady-state methods (Imai, Takei, and Matsumura 1987) require a chemical or biochemical reaction occurring in the liquid to reduce the DO concentration below the saturation level. k L a is then estimated by interpreting the recorded DO concentration with a steady-state mass balance model. The use of a chemical in steady-state methods can strongly affect k L a. Moreover, steady-state methods are prone to errors, and a large error of nearly ±100% can occur if measurements are not conducted correctly. Therefore, steady-state methods are not recommended for large-scale reactors and viscous fluids because of related issues to their cost, accuracy, and validation (Gogate and Pandit 1999; Linek, Beneš, and Vacek 1989).
Geometric parameters (e.g., size and aspect ratio of tank) play a decisive role in scaling up stirred tanks because of their direct influence on k L a. Oxygen transfer due to surface aeration contributes nearly 50% of k L a, especially in small tanks. Surface aeration decreases with increase in tank size that consequently reduces k L a (Humphrey 1998). The aspect ratio of the stirred tank, defined as the ratio of liquid height to tank diameter ( ), can be varied according to the application area. For instance, stirred tanks with the standard geometry ( ) are extensively used in the chemical industry. Standard stirred tanks have low construction costs and usually operate with a single impeller. Conversely, stirred tanks with a non-standard geometry ( ) are preferred for bioprocess applications, such as aerobic fermentation. of non-standard stirred tanks used in bioprocesses is suggested to range within 2–5 (Doran 2013). A higher provides larger residence time to gas bubbles, which increases k L a by improving oxygen solubility. Non-standard stirred tanks are equipped with multiple impellers to provide decent gas dispersion. Experimental studies on standard stirred tanks are extensive, but those on non-standard stirred tanks are limited, especially with non-Newtonian fluids. The authors also recently estimated k L a in standard stirred tanks with working volumes of 11 (Ali and Solsvik 2020) and 50 L (Ali and Solsvik 2021). A detailed literature review of standard stirred tanks has been carried out by Nere, Patwardhan, and Joshi (2003), Grenville and Nienow (2003), and Jaszczur and Młynarczykowska (2020). However, this work focused on a non-standard stirred tank in the presence of non-Newtonian fluids. Studies on k L a measurement in non-standard stirred tanks with non-Newtonian fluids are reviewed and summarized in Table 1.
|Author||V (L)||(−)||N s (rev min−1)||Q (L min−1)||Impeller number and type||Non-Newtonian polymer type and wt%||k L a method|
|Linek, Kordač, and Moucha (2005)||6.0||1.2||330–1130||3–9||Single: RT||CMC: 0.2–0.6%||Steady-state|
|Arjunwadkar et al. (1998)||5.0||1.2||400–750||1.5–5.0||Double: RT + PBT||CMC: 0.7%||Dynamic|
|Puthli, Rathod, and Pandit (2005)||2.0||1.7||300–600||0.5–2.0||Single: RT||CMC: 0.25–0.5%||Dynamic|
|Double: RT + PBT|
|Triple: RT + PBT + PBT|
|Liu et al. (2019)||56||1.4||60–150||17–33||Single: LDB, FZ, MB||XG: 0.5%||Dynamic|
|Liu et al. (2020)||56||1.4||60–150||0–50||Single: LDB, FZ, MB||XG: 0.25–1.0%||Dynamic|
|Jamshidzadeh, Ein-Mozaffari, and Lohi (2020)||62||1.25||100–450||6.0–12.5||Triple: SC + PBT + Anchor||CMC: 0.5–1.5%||Dynamic|
|Gabelle et al. (2011)||21‒42, 170–340||1, 2||100–1000||8–330||Single: RT, PBT, MP||CMC, XG: 0.25, 0.5%||Dynamic|
|Double: RT + PBT, RT + MP, RT + RT, MP + MP|
|This work||11, 47, 40||1, 3.5||800–1000||8–10||Single: RT||CMC: 0.2–0.6%||Dynamic|
A limited number of studies on k L a have been conducted in laboratory-scale stirred tanks with a working volume (V) of less than 10 L (Arjunwadkar et al. 1998; Linek, Kordač, and Moucha 2005; Puthli, Rathod, and Pandit 2005) (Table 1). Linek, Kordač, and Moucha (2005) used the steady-state method to measure k L a in a 6 L non-standard stirred tank ( ). Aqueous solutions of carboxymethylcellulose sodium salt (CMC) with different concentrations were used as non-Newtonian fluids. The rheological data of the non-Newtonian fluids are not reported. The physical properties (i.e., density, viscosity, and molecular weight) of CMC polymer used to prepare non-Newtonian fluids are also not provided, thus preventing the identification of the polymer type (Nishikawa et al. 1981). Polymers for preparing non-Newtonian fluids are generally graded on different viscosities (or molecular weights). A small change in polymer concentration can drastically influence the rheological behavior of non-Newtonian fluids. Therefore, providing the physical properties of the polymer enables us to compare the reported results with those of future work. Puthli, Rathod, and Pandit (2005) and Arjunwadkar et al. (1998) conducted experiments in a 2 L ( ) and 5 L ( ) non-standard stirred tank, respectively. k L a was measured in aqueous CMC solutions of different concentrations by using the dynamic method. However, the rheological data of non-Newtonian fluids and physical properties of CMC polymer are not provided (Arjunwadkar et al. 1998; Puthli, Rathod, and Pandit 2005).
Very few studies have also focused pilot-scale stirred tanks with working volumes ranging from 56 to 340 L (Table 1). Liu et al. (2019, 2020 measured k L a in a 56 L non-standard stirred tank ( ). The dynamic method is used to measure k L a in aqueous xanthan gum (XG) solutions of different concentrations. The flow behavior of XG is described by the power-law model. Similarly, Jamshidzadeh, Ein-Mozaffari, and Lohi (2020) conducted experiments on a 62 L stirred tank with a of 1.25. The shear-thinning behavior of CMC solutions with different concentrations is described by the power-law model. However, the physical properties of XG and CMC polymers are not reported by Liu et al. (2019, 2020 and Jamshidzadeh, Ein-Mozaffari, and Lohi (2020), respectively. Furthermore, Gabelle et al. (2011) studied the effects of tank size on k L a by using two stirred tanks with standard ( ) and non-standard ( ) geometries. k L a of aqueous CMC and XG solutions is obtained using a dynamic method. The maximum and minimum working volumes of the tanks are 42 and 340 L, respectively. The rheological behavior of shear-thinning fluids is expressed by the power-law model. The physical properties of CMC and XG polymers are also provided (Gabelle et al. 2011).
Studies focusing on the effects of geometric parameters on k L a with non-Newtonian fluids are few (Table 1), particularly those using non-standard pilot-scale stirred tanks limited to of 1–2. Therefore, in view of the bioprocess industry, the effect of geometric parameters on k L a of non-Newtonian fluids must be investigated using pilot-scale stirred tanks with higher than 2. The results can help determine the behavior of k L a during scale-up, as well as the optimal tank size and . To fill this research gap, the present study aimed to investigate the effects of geometric parameters (tank size and ) on k L a of shear-thinning fluids by operating standard and non-standard stirred tanks. Three stirred tanks with different sizes (laboratory and pilot-scale) and geometry (standard and non-standard) were used. The values of were 1 for standard and 3.5 for nonstandard stirred tanks. The volumes of the stirred tanks were 11 L ( ), 40 L ( ), and 47 L ( ). k L a was estimated using a dynamic method based on recorded DO concentrations in liquids. Additionally, k L a was evaluated with the effects of fluid rheology and operating conditions. The effects of different types of polymers (for preparing shear-thinning fluids) on k L a were examined. The obtained results of k L a were compared with existing experimental data in literature and those obtained by empirical correlation.
2 Materials and methods
2.1 Experimental setup
The experiments were conducted in three flat-bottomed stirred tanks made of plexiglass with different geometrical configurations (Figure 1). Tank A was a laboratory-scale model with a standard geometry ( ), and working volume of 11 L. Tank B was a pilot-scale model, which was scaled up from Tank A. Tank B had a standard geometry ( ) with a working volume of 47 L. Tank C was also a pilot-scale model, but with nonstandard geometry ( ) and a working volume of 40 L. The detailed geometric configurations of the stirred tanks are presented in Table 2. Recirculation zones were reduced by equipping the stirred tanks with four equispaced (90° apart) rectangular baffles. Mixing in the tanks was achieved with a six-blade Rushton turbine mounted to IKA EUROSTAR 200 control (Tank A) and IKA RW 47 digital (Tank B and Tank C) overhead stirrers. Tank A and Tank B both have one impeller, whereas Tank C had four, as shown in Figure 1(c). The gas was introduced to Tank A and Tank C with a ring sparger. By contrast, a C-type ring sparger was used in Tank B. Impeller speeds and sparger gas flow rates ranged from 800 rev min−1 to 900 rev min−1 and from 8 L min−1 to 10 L min−1, respectively. The selected operating conditions facilitated the supply of gas from each sparger hole and the dispersion of the gas bubbles across the entire tank. The bubbly flow regime was observed within the range at which experiments were conducted (Figure 1). The measurement positions (location of O2 probe) of DO concentrations inside the stirred tanks are shown in Figure 1. All the experiments were repeated three times for each selected impeller speed and gas flow rate. The repeatability of the presented results was checked through statistical analysis. Standard error bars were plotted for all the experimental results, but the error bars may not be visible for a specific result because of the low standard errors for the repeatability of measurements. The experiments were conducted at room temperature conditions (22 °C ± 1 °C).
|Geometry||Volume (L)||(−)||T (cm)||H (cm)||H L (cm)||C (cm)||H p (cm)||D (cm)||d (cm)|
The tested fluids were deionized water and non-Newtonian solution made of CMC. Four types of CMC polymers (Sigma-Aldrich, Norway) based on different polymer viscosities were used in preparing non-Newtonian fluids: low (viscosity = 0.05–0.2 Pa s), medium (viscosity = 0.4–1.0 Pa s), high (viscosity = 1.5–3.0 Pa s), and ultra-high (viscosity = 1.5–4.5 Pa s). Differences in the viscosities of the CMC polymers influenced the apparent viscosities and transparency of the non-Newtonian fluids. Low, medium, and high-viscous CMCs are suspending agents, i.e., the particles stay suspended in the non-Newtonian solutions that makes them opaque. However, ultra-high-viscous CMC is highly purified and less suspended; thus, the non-Newtonian solutions made from ultra-high-viscous CMC are transparent. Moreover, four mass concentrations (wt%) of CMC were used in examining shear-thinning behavior on mass transfer rate: 0.2%, 0.3%, 0.4%, and 0.6%. The non-Newtonian fluid samples of CMC were prepared with deionized water. The apparent viscosities of the samples were measured using a rheometer (Physcia MCR 301, Anton Paar GmbH, Austria) with a concentric cylindrical geometry. The samples were given a rest time of 10 min before each measurement for them to reach thermal equilibrium (temperature = 22 °C) and for the removal of the smallest trapped bubbles in the fluids. Each rheology measurement was repeated twice. The Metzner–Otto correlation was used in relating average shear rate in the stirred tank to stirring speed:
where is the applied shear rate (s−1), k s is the Metzner constant (−), and N s is the stirring speed (s−1). The k s values for six-blade Rushton turbine were assumed 11.5 (Metzner and Otto 1957). Given this k s value, a shear rate in the range of 0.001–1000 s−1 was used in investigating the shear-thinning behaviors of the non-Newtonian samples. Rheological data (shear stress and apparent viscosity) for the CMC solutions are plotted in Figure 2. Non-Newtonian solutions made from high-viscous and ultra-high-viscous CMC polymers showed shear-thinning behavior, as shown in Figures 2(c) and 2(d). A noise from the rheometer was recorded for the liquids whose apparent viscosities are comparable to water [Figure 2(a) and (b)]. Thus, the solutions prepared from high-viscous and ultra-high-viscous CMC polymers were selected for the measurement of k L a. The apparent viscosities of the CMC solutions are represented using the Carreau model:
where represents apparent fluid viscosity (Pa s), and are fluid viscosities at minimum and maximum shear rates (Pa s), is the time at which the behavior of a fluid changes from Newtonian to power-law (s), and n is the power-law index (−). The physical properties of the experimental fluids are presented in Table 3. Reynolds number based on impeller diameter was in the range of 400–75,000 for Tanks A and C, whereas 1000–182,000 for Tank B.
|Fluid||CMC polymer||μ a × 10−3 (Pa-s)||Carreau model|
2.2.2 Volumetric mass transfer coefficient
The k L a values in the stirred tanks were predicted using a dynamic gassing-in method (Garcia-Ochoa and Gomez 2009). Each tank was first deoxygenated with sparging nitrogen gas until DO concentration was less than 2%. Cleaned compressed air was then supplied until DO concentration reached 98%. The opening and closing of the nitrogen and air valves were controlled with the LabVIEW software. DO concentration was recorded every 5 s with an optical O2 probe (InPro® 6870i, Mettler Toledo, Switzerland). Recorded data were monitored and saved using a multiparameter transmitter (M800, Mettler Toledo, Switzerland). The positioning of the O2 probe inside each stirred tank is illustrated in Figure 1. The O2 probe was calibrated with a two-point calibration method before each experiment. Two-point calibration involved calibration in air (DO concentration = 100%) followed by calibration in nitrogen (DO concentration = 0%). This calibration procedure provided the most accurate calibration curve over the entire measuring range. of the O2 probe was significantly less than 1/k L a. Therefore, the effects of the O2 probe were neglected in this work (de Jesus, Moreira Neto, and Maciel Filho 2017; Michelin et al. 2013). For small k L a values, such as those in this work, the dynamic gassing-in method can be employed provided that O2 only is considered during air absorption (Linek et al. 1982; Scargiali et al. 2007, 2010). Thus, the rate of change in oxygen concentration in a liquid can be described using the following equation (Badino, Facciotti, and Schmidell 2001; Nauha et al. 2015):
where C t is the measured local DO concentration in a liquid (mol L−1) and is the saturation concentration of oxygen (mol L−1). Henry’s law was used in determining .
where K H represents Henry’s Law constant (mol L−1 atm), is the partial pressure of oxygen (atm), is oxygen molar fraction (%), P air is the partial pressure of air (atm), and is the partial pressure of water vapors for 100% relative humidity (atm). The value of varied for each experiment depending on the fluid rheology and operating conditions. The variation in was negligible with liquid height (Ali and Solsvik 2020). Equation (3) was integrated as a linear equation.
where t is the time (s), C t = 0 is the DO concentration that is less than 5% at t = 0 s (mol L−1), and k L a is estimated by plotting against t with a linear regression. MATLAB (R2019b, MathWorks, USA) was used (de Jesus, Moreira Neto, and Maciel Filho 2017; Michelin et al. 2013).
2.2.3 Power input
The power input of each stirred tank was estimated using the correlation proposed by Prasher and Wills (1973)
where is the power input per unit liquid volume (W m−3), P o is the power number (−), N s is the impeller speed (s−1), D is the diameter of a stirrer (m), T is tank inner diameter (m), H is liquid height (m), and Q is gas flow rate (m3 s−1). P o of the Rushton turbine stirrer was estimated from the characteristic curve of the Power number–Reynolds number for each experimental fluid (Bates, Fondy, and Corpstein 1963). Taghavi et al. (2011) and Devi and Kumar (2017) showed that the correlation (Eq. 6) can be used in estimating the total power input of aerated stirred tanks with different scales (Devi and Kumar 2017; Taghavi et al. 2011).
3 Results and discussion
The effects of different types of fluids and CMC polymers on the volumetric mass transfer coefficient (k L a) are plotted in Figure 3. The bubbles produced in water were large compared to that of CMC, whereas the bubble count was significantly higher in CMC than in water (Ali and Solsvik 2020, 2021). However, the higher viscosity of CMC solutions decreased the fluid mixing (turbulence), which consequently reduced the rate of bubbles dissolution (DO concentration) in the liquid. Thus, the k L a value in water was approximately two and three times greater than the values in high-viscous CMC 0.3% and ultra-high-viscous CMC 0.3%, respectively (Figure 3a). Under the same operating conditions, the decrease in k L a in the scaled-up tanks (Tanks B and C) was attributed to high bubble coalescence caused by low stirrer-induced turbulence (large eddies). Large bubbles cannot dissolve easily in liquids, thereby reducing DO concentration and k L a. However, the higher number of impellers in Tank C prompted bubble breakage by increasing the production of small-scale eddies, thereby improving k L a. The k L a of water in Tank A was 13 and 1% higher than the values in Tank B and Tank C, respectively. The k L a value in Tank C was 9.8% higher than that in Tank B. Similarly, the k L a of the high-viscous CMC 0.3% in Tank A was 11 and 3.7% higher than the values in Tanks B and C, respectively, and the value in Tank C was 10% higher than that in Tank B. However, the k L a of ultra-high-viscous CMC 0.3% in Tank A was 11% lower and 3.3% higher than the values in Tanks B and C, respectively. The k L a value in Tank C was 14% lower than that in Tank B. The type and concentration of CMC polymer negatively affected k L a, as shown in Figure 3b. The k L a values in high-viscous CMC solutions were two times higher than the values in ultra-high-viscous CMC. In addition, k L a decreased as CMC polymer concentration increased. An increase in CMC concentration (in both types of CMC polymers) from 0.2% to 0.3% significantly impacted k L a. However, k L a was only marginally impacted by an increase in CMC concentration from 0.4% to 0.6%, especially in ultra-high-viscous CMC. This suggests that the system reached a limit, and that additional increases in CMC concentration did not have significant effects on k L a due to the reduced driving force. The k L a values showed somewhat similar trends; the highest values were found in Tank A and the lowest in Tank C. However, from a physical point of view, the deviation of k L a from the general trend was due to experimental error.
3.2 Operating conditions
The effects of impeller speed (N s) on k L a in different CMC polymers and concentrations are presented in Figure 4. With a fixed Q of 8 L min−1, an increase in N s enhanced the turbulence which in return promoted the bubble breakage. The size of the bubbles decreased, thereby increasing their dissolution into the liquid. Thus, increasing the magnitude of N s improved DO concentration and subsequently increased k L a. The effect of increasing N s on k L a was considerably stronger in water; this effect was least pronounced in the ultra-high-viscous CMC. An increase in N s at low CMC concentrations (0.2 and 0.3%) yielded significant differences in k L a. However, at high CMC concentrations (0.4 and 0.6%), especially in ultra-high-viscous CMC, the impact of N s on k L a was small due to substantial increase in liquid viscosity. The k L a values were plotted as a function of Q for different CMC polymers and concentrations (Figure 5). The measurement of bubble size was out of the scope of the present study and thus was not performed. The understanding of the relation of k L a with bubble size was inferred from our previous work (Ali and Solsvik 2020). N s was maintained constant at 800 rev min−1, thereby increasing Q supported bubble coalescence, and accordingly, the bubble size increased (Ali and Solsvik 2020). Large bubbles are difficult to dissolve in the liquid, thus improving DO concentrations marginally. Moreover, an increase in Q also increased the size of gas cavities behind the impeller blades (Cappello et al. 2020). These effects explain the slight increase in the magnitude of k L a by increasing Q. Increasing Q had a more pronounced effect on k L a in water; this effect was much smaller in viscous fluids. The effect of Q on k L a was considerable at low CMC concentrations (0.2 and 0.3%) and nominal at high CMC concentrations (0.4 and 0.6%). These findings imply that the impacts of operating conditions (N s and Q) onk L a depend on rheology (liquid viscosity). Therefore, increasing the operating conditions in high CMC concentrations will not significantly improve k L a.
3.3 Experimental comparison
k L a measured in the present study was compared with the results of our laboratory-scale (V = 11 L) (Ali and Solsvik 2020) and pilot-scale (V = 50 L) (Ali and Solsvik 2021) standard geometry stirred tanks and with results for standard and non-standard geometry stirred tanks found in the literature (Figure 6). Student’s t-test (SigmaPlot, 14.0, Systat Software Inc., USA) was used to evaluate the standardized mean difference (SMD) between the results. The comparison with k L a presented in Figure 6(a) uses experimental data obtained in the standard stirred tanks. For Newtonian fluids, the values of k L a reported in this work were close to those of our laboratory-scale stirred tank (Ali and Solsvik 2020) (SMD = 2%). The SMD between the values of k L a obtained in this work and those from our pilot-scale stirred tank was approximately 3.1% (Ali and Solsvik 2021). The deviation between the results was due to differences in impeller bottom clearance (C) and operating conditions.
The SMDs between the results of this work and those of Petříček et al. (2018) and Cappello et al. (2020) were approximately 10.4 and 15.4%, respectively. Comparing k L a in non-Newtonian fluids, the SMD between the results of the present study and our laboratory-scale stirred tank was 1.1% (Ali and Solsvik 2020). Similarly, the SMD between the k L a reported in this work and that of the pilot-scale stirred tank was 1.3% (Ali and Solsvik 2021), but a significant difference was observed between the results of the present study and those of Tecante and Choplin (1993) (SMD = 4.72%). Comparing the values of k L a obtained in this work with those of Cappello et al. (2020) and García-Ochoa and Castro (2001) produced SMDs of 3.51 and 3.03%, respectively.
The k L a results obtained in the non-standard stirred tank (Tank C) were also compared with experimental data reported in the literature, as shown in Figure 6(b). Regarding Newtonian fluids, the measured k L a closely matched the results of both Arjunwadkar et al. (1998) (SMD = 1.1%) and Bouaifi and Roustan (1998) (SMD = 1.2%). The differences between the values of k L a reported in the present study and those of Badino, Facciotti, and Schmidell (2001) and Linek et al. (2004) were 2 and 3.6%, respectively. The SMD when comparing the results with those of Kazemzadeh et al. was high (7.1%) (Kazemzadeh et al. 2020), but the mean difference was largest (16.3%) between the measured k L a in non-Newtonian fluids and the experimental data reported by Linek, Kordač, and Moucha (2005). The results of the present study were in close agreement with Jamshidzadeh, Ein-Mozaffari, and Lohi (2020) and Liu et al. (2020), with SMDs of less than 1%, and the mean difference between the results reported in the present work and those of Gabelle et al. (2011) was about 2.5%. Different operating conditions, reactor geometry, water quality, and polymer type were among the reasons that the results of this work varied from those reported in the literature.
3.4 Empirical correlation
The experimental k L a values reported in this work were compared to those estimated from the standard empirical correlation (Figure 7). The standard correlation for reliably predicting k L a in non-Newtonian fluids was defined as follows (Gabelle et al. 2011; Jamshidzadeh, Ein-Mozaffari, and Lohi 2020):
where C 1 is the empirical coefficient; a, b, and c are the empirical exponents; is the power input per unit liquid volume (W m−3); U g is the superficial gas velocity (m s−1); and μ a is apparent fluid viscosity (Pa s). The values of C 1, a, b, and c were calculated using nonlinear regression analysis and are presented in Table 4. The results predicted by the empirical correlation correlated to a reasonable degree with the experimental data reported here. When the empirical correlation was used to predict k L a values in water, the standard errors were 15.2%, 10.2%, and 12.4% for Tanks A, B, and C, respectively. The empirical correlation predicted our findings quite well in the case of high-viscous CMC; the standard errors were 7.8%, 4.6%, and 7.2% for Tanks A, B, and C, respectively. The standard errors decreased to 4.5% (Tank A), 4.2% (Tank B), and 3.8% (Tank C) for ultra-high-viscous CMC. These results indicate that the empirical correlation depends on tank size (scale). The empirical correlation cannot predict k L a for all experimental data (Tanks A, B, and C) with the same values of coefficients and exponents. Thus, the correlation must be modified to include the effects of tank size or scale on mass transfer.
|Fluids||CMC polymer||C 1||a||b||c||SE|
|Water||–||7 × 10−3||0.120||0.004||−0.079||0.027|
|CMC 0.2%||High-viscous||7.2 × 10−6||0.122||−0.115||−1.829||0.039|
|CMC 0.3%||High-viscous||1.2 × 10−5||0.039||0.117||−2.438||0.016|
|CMC 0.4%||High-viscous||2.0 × 10−7||−0.145||0.051||−4.437||0.015|
|CMC 0.6%||High-viscous||2.1 × 10−4||−0.020||−0.105||−1.536||0.017|
|CMC 0.2%||Ultra-high-viscous||7.4 × 10−8||0.049||0.0515||−3.429||0.023|
|CMC 0.3%||Ultra-high-viscous||5 × 10−9||0.0132||0.0350||−4.829||0.027|
|CMC 0.4%||Ultra-high-viscous||9.4 × 10−9||0.0062||0.0466||−5.446||0.024|
|CMC 0.6%||Ultra-high-viscous||6.4 × 10−6||0.1761||−0.246||−2.381||0.037|
The presents study investigated the effects of geometric parameters on k L a using three stirred tanks with different sizes and geometries. The k L a values in shear-thinning fluids were estimated from the logged data of DO concentration. The effects of rheology and operating conditions on k L a were measured. The results of k L a obtained in this work were then compared to previously published experimental data and empirical correlation calculations.
Rheology significantly affected k L a, since a significant decrease in k L a was observed in viscous shear-thinning fluids. CMC concentration and polymer type also influenced k L a. An increase in CMC concentration decreased k L a considerably. High-viscous CMC resulted in a k L a that was two times higher than that in ultra-high-viscous CMC. However, the k L a values changed only slightly at higher CMC concentrations (0.4 and 0.6%) because of the reduced driving force once the system reached a limit. A significant decrease in k L a was observed as tank size increased. The value of k L a in Tank A was higher than those in Tanks B and C. The higher number of impellers in Tank C substantially improved k L a by increasing turbulence. Increasing the operating conditions (N s and Q) increased k L a, and N s impacted k L a more significantly than Q. High operating conditions increased bubble breakage, which consequently enhanced k L a. The effect of operating conditions on k L a was high at low CMC concentrations but small at high CMC concentrations. The k L a values obtained in this work align closely with those found with our laboratory (Ali and Solsvik 2020) and pilot-scale (Ali and Solsvik 2021) stirred tanks. However, they differ significantly to the findings of other published experimental data. Operating conditions, tank size and geometry, water quality, and polymer type (for non-Newtonian fluids) are decisive parameters that must be carefully considered when predicting k L a in stirred tanks. Empirical correlation calculations aligned with our experimental findings when different coefficient and exponent values were used for each stirred tank. Thus, empirical correlation calculations must be modified to fit all experimental data, and the same values of coefficient and exponents should be used. The results of the present study can be used to validate studies of computational fluid dynamics (CFD), such as the CFD-based population model (Niño et al. 2020). More detailed investigations of the effects of on k L a are needed to reliably describe the relationship between scale-up and mass transfer.
empirical exponent (−)
empirical exponent (−)
empirical exponent (−)
impeller clearance (m)
- C t
oxygen concentration in the liquid (mol L−1)
- C 1
empirical coefficient (−)
saturation concentration of oxygen (mol L−1)
- C t=0
oxygen concentration at t = 0 (mol L−1)
sparger diameter (m)
stirrer diameter (m)
tank height (m)
- H L
liquid height (m)
aspect ratio (−)
- H p
probe height (m)
- k s
Metzner constant (−)
- k L a
volumetric mass transfer coefficient (s−1)
- 1/k L a
characteristic time (s)
- K H
Henry’s Law constant (mol L−1 atm)
- N s
stirring speed (s−1)
power-law index (−)
- P g
power input (W)
- P o
power number (−)
- P air
partial pressure of air (atm)
partial pressure of water vapors (atm)
partial pressure of oxygen (atm)
gas flow rate (m3 s−1)
tank diameter (m)
- U g
superficial gas velocity (m s−1)
liquid volume (m−3)
oxygen molar fraction (%)
- Greek symbols
shear rate (s−1)
time response of O2 probe (s)
time constant of Carreau model (s)
- μ a
apparent fluid viscosity (Pa s)
- μ 0
fluid viscosity at minimum shear rate (Pa s)
- μ ∞
fluid viscosity at maximum shear rate (Pa s)
Funding source: Research Council of Norway
Award Identifier / Grant number: 274398
Author contributions: All the authors have accepted responsibility for the entire content of this submitted manuscript and approved submission.
Research funding: The work was supported by the Research Council of Norway (Grant no. 274398).
Conflict of interest statement: The authors declare no conflicts of interest regarding this article.
Ali, H., and J. Solsvik. 2020. “Axial Distributions of Bubble-Liquid Mass Transfer Coefficient in Laboratory-Scale Stirred Tank with Viscous Newtonian and Non-Newtonian Fluids.” Physics of Fluids 32 (12): 123308, https://doi.org/10.1063/5.0031574. Search in Google Scholar
Ali, H., and J. Solsvik. 2021. “Bubble Hydrodynamics and Mass Transfer in Stirred Tank with Non-Newtonian Fluids: Scale-Up from Laboratory to Pilot-Scale.” Physics of Fluids 33 (3): 033319, https://doi.org/10.1063/5.0045425. Search in Google Scholar
Arjunwadkar, S. J., K. Sarvanan, P. R. Kulkarni, and A. B. Pandit. 1998. “Gas–Liquid Mass Transfer in Dual Impeller Bioreactor.” Biochemical Engineering Journal 1 (2): 99–106, https://doi.org/10.1016/s1385-8947(97)00083-1. Search in Google Scholar
Badino, A. C., M. C. R. Facciotti, and W. Schmidell. 2001. “Volumetric Oxygen Transfer Coefficients (kLa) in Batch Cultivations Involving Non-Newtonian Broths.” Biochemical Engineering Journal 8 (2): 111–9, https://doi.org/10.1016/s1369-703x(01)00092-4. Search in Google Scholar
Bandyopadhyay, B., A. E. Humphrey, and H. Taguchi. 1967. “Dynamic Measurement of the Volumetric Oxygen Transfer Coefficient in Fermentation Systems.” Biotechnology and Bioengineering 9 (4): 533–44, https://doi.org/10.1002/bit.260090408. Search in Google Scholar
Bates, R. L., P. L. Fondy, and R. R. Corpstein. 1963. “Examination of Some Geometric Parameters of Impeller Power.” Industrial & Engineering Chemistry Process Design and Development 2 (4): 310–4, https://doi.org/10.1021/i260008a011. Search in Google Scholar
Bouaifi, M., and M. Roustan. 1998. “Bubble Size and Mass Transfer Coefficients in Dual-Impeller Agitated Reactors.” The Canadian Journal of Chemical Engineering 76 (3): 390–7, https://doi.org/10.1002/cjce.5450760307. Search in Google Scholar
Cappello, V., C. Plais, C. Vial, and F. Augier. 2020. “Bubble Size and Liquid-Side Mass Transfer Coefficient Measurements in Aerated Stirred Tank Reactors with Non-Newtonian Liquids.” Chemical Engineering Science 211: 115280, https://doi.org/10.1016/j.ces.2019.115280. Search in Google Scholar
Devi, T. T., and B. Kumar. 2017. “Mass Transfer and Power Characteristics of Stirred Tank with Rushton and Curved Blade Impeller.” Engineering Science and Technology, an International Journal 20 (2): 730–7, https://doi.org/10.1016/j.jestch.2016.11.005. Search in Google Scholar
Doran, P. M. 2013. “Chapter 8 – Mixing.” In Bioprocess Engineering Principles, edited by P. M. Doran, 255–332. Oxford: Academic Press. Search in Google Scholar
Dunn, I. J., and A. Einsele. 2007. “Oxygen Transfer Coefficients by the Dynamic Method.” Journal of Applied Chemistry and Biotechnology 25 (9): 707–20, https://doi.org/10.1002/jctb.5020250910. Search in Google Scholar
Gabelle, J. C., F. Augier, A. Carvalho, R. Rousset, and J. Morchain. 2011. “Effect of Tank Size on kLa and Mixing Time in Aerated Stirred Reactors with Non-Newtonian Fluids.” Canadian Journal of Chemical Engineering 89 (5): 1139–53, https://doi.org/10.1002/cjce.20571. Search in Google Scholar
García-Ochoa, F., and E. G. Castro. 2001. “Estimation of Oxygen Mass Transfer Coefficient in Stirred Tank Reactors Using Artificial Neural Networks.” Enzyme and Microbial Technology 28 (6): 560–9. Search in Google Scholar
Garcia-Ochoa, F., and E. Gomez. 2009. “Bioreactor Scale-Up and Oxygen Transfer Rate in Microbial Processes: An Overview.” Biotechnology Advances 27 (2): 153–76, https://doi.org/10.1016/j.biotechadv.2008.10.006. Search in Google Scholar
Garcia-Ochoa, F., and E. Gomez. 2004. “Theoretical Prediction of Gas–Liquid Mass Transfer Coefficient, Specific Area and Hold-Up in Sparged Stirred Tanks.” Chemical Engineering Science 59 (12): 2489–501, https://doi.org/10.1016/j.ces.2004.02.009. Search in Google Scholar
García-Ochoa, F., and E. Gómez. 1998. “Mass Transfer Coefficient in Stirred Tank Reactors for Xanthan Gum Solutions.” Biochemical Engineering Journal 1 (1): 1–10, https://doi.org/10.1016/s1369-703x(97)00002-8. Search in Google Scholar
Gogate, P. R., and A. B. Pandit. 1999. “Survey of Measurement Techniques for Gas–Liquid Mass Transfer Coefficient in Bioreactors.” Biochemical Engineering Journal 4 (1): 7–15, https://doi.org/10.1016/s1369-703x(99)00033-9. Search in Google Scholar
Grenville, R. K., and A. W. Nienow. 2003. “Blending of Miscible Liquids.” In Handbook of Industrial Mixing, 507–42. Oxford: Academic Press. Search in Google Scholar
Heijnen, J. J., K. V. Riet, and A. J. Wolthuis. 1980. “Influence of Very Small Bubbles on the dynamic kLA Measurement in Viscous Gas–Liquid Systems.” Biotechnology and Bioengineering 22 (9): 1945–56, https://doi.org/10.1002/bit.260220912. Search in Google Scholar
Imai, Y., H. Takei, and M. Matsumura. 1987. “A Simple Na2SO3 Feeding Method for KLa Measurement in Large-Scale Fermentors.” Biotechnology and Bioengineering 29 (8): 982–93, https://doi.org/10.1002/bit.260290810. Search in Google Scholar
Jamshidzadeh, M., F. Ein-Mozaffari, and A. Lohi. 2020. “Experimental Analysis of the Mass Transfer Coefficient and Interfacial Area in an Aerated Coaxial Mixing System Comprising a Non-Newtonian Solution.” Industrial & Engineering Chemistry Research 59 (49): 21530–47, https://doi.org/10.1021/acs.iecr.0c03641. Search in Google Scholar
Jaszczur, M., and A. Młynarczykowska. 2020. “A General Review of the Current Development of Mechanically Agitated Vessels.” Processes 8 (8): 982, https://doi.org/10.3390/pr8080982. Search in Google Scholar
de Jesus, S. S., J. Moreira Neto, and R. Maciel Filho. 2017. “Hydrodynamics and Mass Transfer in Bubble Column, Conventional Airlift, Stirred Airlift and Stirred Tank Bioreactors, Using Viscous Fluid: A Comparative Study.” Biochemical Engineering Journal 118: 70–81, https://doi.org/10.1016/j.bej.2016.11.019. Search in Google Scholar
Kazemzadeh, A., C. Elias, M. Tamer, A. Lohi, and F. Ein-Mozaffari. 2020. “Mass Transfer in a Single-Use Angled-Shaft Aerated Stirred Bioreactor Applicable for Animal Cell Culture.” Chemical Engineering Science 219: 115606, https://doi.org/10.1016/j.ces.2020.115606. Search in Google Scholar
Linek, V., P. Beneš, and V. Vacek. 1989. “Dynamic Pressure Method for kLa Measurement in Large‐Scale Bioreactors.” Biotechnology and Bioengineering 33 (11): 1406–12, https://doi.org/10.1002/bit.260331107. Search in Google Scholar
Linek, V., P. Beneš, V. Vacek, and F. Hovorka. 1982. “Analysis of Differences in kLa Values Determined by Steady-State and Dynamic Methods in Stirred Tanks.” The Chemical Engineering Journal 25 (1): 77–88, https://doi.org/10.1016/0300-9467(82)85024-7. Search in Google Scholar
Linek, V., M. Kordač, M. Fujasová, and T. Moucha. 2004. “Gas–Liquid Mass Transfer Coefficient in Stirred Tanks Interpreted Through Models of Idealized Eddy Structure of Turbulence in the Bubble Vicinity.” Chemical Engineering and Processing: Process Intensification 43 (12): 1511–7, https://doi.org/10.1016/j.cep.2004.02.009. Search in Google Scholar
Linek, V., M. Kordač, and T. Moucha. 2005. “Mechanism of Mass Transfer from Bubbles in Dispersions Part II: Mass Transfer Coefficients in Stirred Gas–Liquid Reactor and Bubble Column.” Chemical Engineering and Processing: Process Intensification 44 (1): 121–30, https://doi.org/10.1016/j.cep.2004.05.009. Search in Google Scholar
Liu, B., Q. Xiao, P. Gao, B. Sunden, and F. Fan. 2020. “Investigation of Gas–Liquid Dispersion and Mass Transfer Performance of Wide-Viscosity-Range Impellers in Water Solutions of Xanthan Gum.” Chemical Engineering Research and Design 154: 60–9, https://doi.org/10.1016/j.cherd.2019.12.005. Search in Google Scholar
Liu, B., Q. Xiao, N. Sun, P. Gao, F. Fan, and B. Sunden. 2019. “Effect of Gas Distributor on Gas–Liquid Dispersion and Mass Transfer Characteristics in Stirred Tank.” Chemical Engineering Research and Design 145: 314–22, https://doi.org/10.1016/j.cherd.2019.03.035. Search in Google Scholar
Michelin, M., A. M. de Oliveira Mota, M. L. T. M. Polizeli, D. P. da Silva, A. A. Vicente, and J. A. Teixeira. 2013. “Influence of Volumetric Oxygen Transfer Coefficient (kLa) on Xylanases Batch Production by Aspergillus niger Van Tieghem in Stirred Tank and Internal-Loop Airlift Bioreactors.” Biochemical Engineering Journal 80: 19–26, https://doi.org/10.1016/j.bej.2013.09.002. Search in Google Scholar
Muller, F. L., and J. F. Davidson. 1992. “On the Contribution of Small Bubbles to Mass Transfer in Bubble Columns Containing Highly Viscous Liquids.” Chemical Engineering Science 47 (13–14): 3525–32, https://doi.org/10.1016/0009-2509(92)85066-k. Search in Google Scholar
Nauha, E. K., O. Visuri, R. Vermasvuori, and V. Alopaeus. 2015. “A New Simple Approach for the Scale-Up of Aerated Stirred Tanks.” Chemical Engineering Research and Design 95: 150–61, https://doi.org/10.1016/j.cherd.2014.10.015. Search in Google Scholar
Nere, N. K., A. W. Patwardhan, and J. B. Joshi. 2003. “Liquid-phase Mixing in Stirred Vessels: Turbulent Flow Regime.” Industrial and Engineering Chemistry Research 42 (12): 2661–98, https://doi.org/10.1021/ie0206397. Search in Google Scholar
Niño, L., R. Gelves, H. Ali, J. Solsvik, and H. Jakobsen. 2020. “Applicability of a Modified Breakage and Coalescence Model Based on the Complete Turbulence Spectrum Concept for CFD Simulation of Gas–Liquid Mass Transfer in a Stirred Tank Reactor.” Chemical Engineering Science 211: 115272, https://doi.org/10.1016/j.ces.2019.115272. Search in Google Scholar
Nishikawa, M., M. Nakamura, K. Hashimoto, and M. Nishikawa. 1981. “Gas Absorption in Aerated Mixing Vessels with Non-Newtonian Liquid.” Journal of Chemical Engineering of Japan 14 (3): 227–32, https://doi.org/10.1252/jcej.14.227. Search in Google Scholar
Petříček, R., T. Moucha, F. J. Rejl, L. Valenz, J. Haidl, and T. Čmelíková. 2018. “Volumetric Mass Transfer Coefficient, Power Input and Gas Hold-Up in Viscous Liquid in Mechanically Agitated Fermenters. Measurements and Scale-Up.” International Journal of Heat and Mass Transfer 124: 1117–35. Search in Google Scholar
Prasher, B. D., and G. B. Wills. 1973. “Mass Transfer in an Agitated Vessel.” Industrial & Engineering Chemistry Process Design and Development 12 (3): 351–4, https://doi.org/10.1021/i260047a023. Search in Google Scholar
Puthli, M. S., V. K. Rathod, and A. B. Pandit. 2005. “Gas–Liquid Mass Transfer Studies with Triple Impeller System on a Laboratory Scale Bioreactor.” Biochemical Engineering Journal 23 (1): 25–30, https://doi.org/10.1016/j.bej.2004.10.006. Search in Google Scholar
Scargiali, F., A. Busciglio, F. Grisafi, and A. Brucato. 2010. “Simplified Dynamic Pressure Method for kLa Measurement in Aerated Bioreactors.” Biochemical Engineering Journal 49 (2): 165–72, https://doi.org/10.1016/j.bej.2009.12.008. Search in Google Scholar
Scargiali, F., R. Russo, F. Grisafi, and A. Brucato. 2007. “Mass Transfer and Hydrodynamic Characteristics of a High Aspect Ratio Self-Ingesting Reactor for Gas–Liquid Operations.” Chemical Engineering Science 62 (5): 1376–87, https://doi.org/10.1016/j.ces.2006.11.040. Search in Google Scholar
Taghavi, M., R. Zadghaffari, J. Moghaddas, and Y. Moghaddas. 2011. “Experimental and CFD Investigation of Power Consumption in a Dual Rushton Turbine Stirred Tank.” Chemical Engineering Research and Design 89 (3): 280–90, https://doi.org/10.1016/j.cherd.2010.07.006. Search in Google Scholar
Tecante, A., and L. Choplin. 1993. “Gas–Liquid Mass Transfer in Non-Newtonian Fluids in a Tank Stirred with a Helical Ribbon Screw Impeller.” The Canadian Journal of Chemical Engineering 71 (6): 859–65, https://doi.org/10.1002/cjce.5450710606. Search in Google Scholar
© 2021 Haider Ali et al., published by De Gruyter, Berlin/Boston
This work is licensed under the Creative Commons Attribution 4.0 International License.