Electronic Journal of Biotechnology
versión On-line ISSN 0717-3458
Electron. J. Biotechnol. v.11 n.1 Valparaíso ene. 2008
Oxygen mass transfer studies on batch cultivation of P. aeruginosa in a biocalorimeter
Financial support: Council for Scientific and Industrial Research (CSIR) Fellowship.
Keywords: biocalorimeter, metabolic heat, P. aeruginosa, scale-up, volumetric mass transfer coefficient.
In the present work volumetric mass transfer coefficient (kLa) was investigated during batch cultivations of Pseudomonas aeruginosa on a nutrient media. The effects of process variables (viz. impeller speed, oxygen flow and geometry of impeller) on the volumetric mass transfer coefficient of oxygen, kLa, in a biocalorimeter (Bio-RC1) was investigated and reported in this research work. The experimental data have been analyzed employing MATLAB to obtain the influences of the process parameters on kLa. An attempt was made to correlate volumetric mass transfer coefficient with metabolic heat production rate at optimized process conditions. The correlation reported in this work would be useful to control and scale up of bioprocesses.
With the recent activities for developing biological processes on industrial scale, there is a need for evolving robust design and scale-up methods for bioreactors. Bioreactor design is a complex engineering task (Garcia-Ochoa and Gomez, 2000; Garcia-Ochoa and Gomez, 2005). Under optimum conditions the microorganisms or cells will reproduce at an astounding rate. The vessel's geometry and operating conditions like gas (i.e., air, oxygen, nitrogen, carbon dioxide) flow rates, temperature, pH and dissolved oxygen (DO) levels, and agitation speed have to be closely monitored and controlled. One of the key issues in aerobic bioreactor design is oxygen transfer (Galaction et al. 2004). Oxygen transfer in aerobic bioprocesses is essential and shortage of it would drastically affect the process performance. Bioprocesses are usually carried out in aqueous media, where the solubility of oxygen is very low owing to the presence of ionic salts and nutrients. At the same time the rate of oxygen utilization by the microorganisms is rather high. Oxygen transfer through microbial cells controls the most of aerated biological systems. The amount of DO into the broths is limited by its solubility and mass transfer rate, as well as by its consumption rate on cells metabolic pathways (Calik et al. 1997).Oxygen transfer is usually helped by agitation, and is also for mixing nutrients and to keep the system homogeneous. There are, however, limits to the speed of agitation, due to high power consumption (which is proportional to the cube of the speed of the electric motor) and the damage to organisms due to excessive tip speed. Gas-liquid mass transfer is commonly the rate-limiting step in industrial-scale biochemical reactions. For this reason, bioreactor design and scale-up focuses on providing ample gas mass transfer. The oxygen mass transfer can be described and analyzed by means of volumetric mass transfer coefficient, kLa. It represents the most important parameter implied on the design and operation of mixing-sparging equipment of the bioreactors. Maintaining a constant volumetric mass-transfer coefficient is a common scale-up strategy (Ozbek and Gayik, 2001). The oxygen mass transfer rate in a bioreactor depends on several factors, such as the geometry of the bioreactor (type of bioreactor, distributor or stirrer design, etc.), the liquid properties (viscosity, superficial tension, etc.), and the dissipated energy in the fluid, which depends on the airflow rate, the stirrer speed, etc. Therefore, composition and rheological properties of the medium changing with time have an important effect on mass transfer rate (Montes et al. 1999; Chisti et al. 2002). Thus the coefficient (kLa) plays an important role in design, scale-up and economy of the process.
Most of the volumetric mass transfer coefficient correlations in the literature belong to two categories. In the first kind kLa is correlated to stirrer speed (or power input per unit volume) and superficial gas velocity, and in the case of the non-Newtonian fluids the broth viscosity must be also taken into account (Garcia-Ochoa and Gomez, 1998). The second category includes empirical correlations based on dimensionless numbers (Nishikawa et al. 1981), i.e., expressing the dependence of the Sherwood number as a function of the Reynolds, Weber and Aeration numbers. Biocalorimeter (RC1) is a bench scale bioreactor widely applied to investigate several biochemical reactions using the metabolic heat generated due to metabolic activity (Marison et al. 1998). Though lot of studies has been dealt on thermodynamics of microbial growth using biocalorimeter, knowledge of mass transfer occurring simultaneously with metabolic heat generation is feeble (which is more important when the calorimetric results are used for bioreactor scale-up work). The present study aims to determine the influence of process variables such as impeller speed and oxygen flow rate on the gas-liquid mass transfer of oxygen in a biocalorimeter. Correlations for kLa are needed to be developed for biocalorimeter, which could immensely help for efficient scale up of calorimetric results for large scale level. Using a bench-scale biocalorimeter, with a single impeller and a capacity of 1 l (total volume), the effects of process variables on the kLa were investigated. Later, kLa values were analyzed to establish correlations with impeller speed and oxygen flow rate. From the established correlations, the numerical values for the exponents (α, β) were determined. A correlation was also developed between kLa and metabolic heat (q) at optimized process conditions.
The organism used for oxygen transfer studies was P. aeruginosa (BC1) isolated from tannery soak liquor. Isolate BC1 was confirmed as P. aeruginosa by its positive oxidase reaction (Kovacs, 1956), oxidative utilization of glucose (Hugh and Leifson, 1953) and production of pyocyanin on Pseudomonas Isolation Agar (PIA; Himedia, Mumbai, India). The strain was further confirmed with biochemical analysis performed at Institute of Biomedical Sciences (IBMS, Chennai, India). The isolated strain was halotolerant (up to 4% NaCl by wt) in optimized growth conditions and its high proteolytic activity even at saline conditions made it a suitable organism for soak liquor treatment.
P. aeruginosa was maintained in glycerol stock under frozen conditions (-20ºC). A loop full of colonies was streaked on nutrient agar plates and incubated for 24 hrs at 37ºC. Well grown individual pure colonies were transferred aseptically to NB media and incubated in a shaker till it reached its exponential phase, which would serve as inoculum for all the experiments.
Nutrient broth (a complex growth media) was used for cultivation of P. aeruginosa in all oxygen uptake experiments. NB medium consists of Peptone 10 g l-1; Yeast extract 3 g l-1; NaCl 5 g l-1; (pH 7.0). pH adjusted to 7.0 by adding NaOH and HCl solutions respectively using a pH meter. Nutrient agar medium (Himedia, MO12) was used for preservation of pure strain. All mediums were sterilized in an autoclave at 120ºC and 15 psi pressure respectively.
All the experiments have been carried out in a biological reaction calorimeter (Bio-RC1, Mettler-Toledo, AG, Switzerland) with 1 l of working volume (Internal dia 8.2 cm & Height of reaction broth on reactor 15 cm), equipped with DO (Schott, Oxy12) and temperature probes (Figure 1). Agitation was provided with a 4-blade turbine impeller (0-1000 rpm) type with a ratio of stirrer/tank diameters (D/T) of 0.44. Oxygenation was performed by supplying iolar grade pure oxygen from an oxygen cylinder and the flow was controlled by a rotameter. The inlet oxygen was sterilized through a membrane filter (0.2 µm) and sparged in to the bottom of the reactor. In this study, 2 different stirrers viz. Turbine-impeller type (4-bladed and 4 cm dia) and a Rushton-turbine type (4-bladed and 4 cm dia) of 1 mm of thickness were used to investigate the oxygen mass transfer ‘L-bend’ type gas sparger (4 cm long) with holes of 0.5 mm diameter was employed for oxygen supply. Oxygen uptake studies were performed at isothermal mode and the reactant temperature (Tr) was maintained constant by controlling the jacket temperature (Tj). This was done by circulating at high rate (2 ls-1) low-viscosity silicone oil through the reactor jacket. The jacket temperature was carefully controlled by blending oils from ‘hot’ and ‘cold’ circuits, via an electronically controlled metering valve. Heat generated from the biological reactions was monitored by WINRC software (Marison et al. 1998), which acquired the temperatures of both reactor (Tr) and jacket (Tj) and evaluated the value for heat production rate.
The measured heat production rate (q in W) is given by
q = UA (Tr - Tj) 
Where U is the global heat transfer coefficient in (W K-1 m-2), A is the heat transfer surface in m2, Tr is the reactor temperature in K and Tj is the jacket temperature in ºC.
All the inlet and outlet streams were completely insulated to avoid heat loss due to heat bridge phenomenon. Baseline heat signal (qb) value is sum of heat effects due to aeration and stirring power apart from the measured heat of reaction (qr). After reaching a stable baseline heat signal (qb), the reaction was initiated by addition of inoculum and similar kind of stable baseline signal was observed at the end of the experiment (endogenous phase). This baseline heat value (qb) was carefully eliminated from the heat evolved due to biological reaction. The value of metabolic heat evolution during the bioprocess (qr) could simply be determined by deducting the measured q (from Equation 1) from that of the base-line signal. Hence the metabolic heat production rate (qr) is given as follows:
q - qb = UA (ΔT - ΔTb) = qr 
The UA factor was determined for each run at regular time intervals by means of in situ calibration electrical heater provided in calorimeter (Hariri et al. 1996).
The experimental work has been carried out in a complex growth (NB) media, maintaining optimum operational conditions (pH 7.0, 37ºC and 2% seed culture conc.). The oxygen flow rate was varied in the range of 0.5-1.5 l min-1 and the stirrer speed was changed from 50 rpm to 250 rpm. A minimum critical DO concentration up to 10% saturation value (2 ppm) was maintained to ensure pure aerobic conditions i.e. the minimum requirement of DO for obligatory aerobic organism such as P. aeruginosa.
Biomass concentration was determined by centrifuging the culture sample taken from the calorimeter in an ultracentrifuge (Sigma, Model no. 3K30) at 10000 rpm for 10 min. The harvested biomass was dried in an oven for 24 hrs at 105ºC. The cell dry biomass was estimated gravimetrically. DO was monitored online by means of a polarographic electrode obtained from M/s. Schott, Germany (OXY12 Model).
The following mass balance equation for the DO in batch reaction can be established:
where the first term on the right hand side of Equation (3) is the oxygen transfer rate (OTR) and the second term is the volumetric OUR of the culture. The measurement of OUR and OTR can be made using the dynamic technique, in two stages (Figure 2). In the first stage, the inlet of airflow to the broth is interrupted, and a decrease of DO concentration due to cellular respiration is observed, which was recorded by a polarographic oxygen probe. The volumetric OUR is determined by noting the changes in the DO concentration when the oxygen flow was stopped.
Equation (3) can be simplified to:
The DO concentration should be maintained higher than 10% of saturation value to ensure that the microorganisms are not damaged due to lack of oxygen. The volumetric OUR is obtained from the slope of the linear regression of the change in DO concentration vs. time. This first stage is employed to obtain the value of the specific OUR (qO2) by the microorganism, by dividing the volumetric OUR by the biomass concentration. In the second stage the inlet of oxygen flow to the broth is restarted (stirrer speed and airflow rate are fixed at predetermined values) and this lead to increase in DO concentration. Oxygen flow was maintained till the DO concentration () reached the steady state value. Under these conditions Equation (3) can be integrated, taking into account the moment at which the oxygen supply of the culture began () and another time (), yielding the equation:
Equation (5) was used to determine the volumetric mass transfer coefficient. Equation 5 was solved by plotting against steady state DO concentration value () for each data set of experimental values of DO concentration and time, to obtain volumetric mass transfer coefficient. Experiments were carried out for measuring both OUR and OTR at different stages of growth: lag, exponential and stationary stages. Volumetric mass transfer coefficient values, kLa, have been determined in a wide range of operational conditions viz. impeller speed (50-250 rpm), oxygen flow (0.5-1.5 l min-1) and impeller type (Rushton-Turbine & Turbine). The measurements have been carried out changing the superficial gas velocity, Vs, as 0.00165, 0.0003304 and 0.004976 ms-1, which corresponding to oxygen flow rates between 0.5, 1 and 1.5 l min-1. The stirrer speed, N, was varied between 50 and 250 rpm and simultaneously OUR, OTR and kLa values were determined.
The oxygen mass transfer rate in the biocalorimeter with mechanical agitation of broth is a function of many variables, such as the physical properties of the liquid (viscosity, surface tension, etc.), the geometry of the vessel and stirrer, the type of sparger and the operational conditions. In the literature, both dimensional and dimensionless equations for the volumetric mass transfer coefficient have been proposed. In this work, the effect of impeller speed (N), superficial gas velocity (Vs, oxygen flow) and geometry of stirrer on oxygen mass transfer were extensively studied. Predictions about the rate of absorption of a gaseous species in a stirred tank is usually based on correlations of overall volumetric mass transfer coefficient (kLa) with mechanical agitation power per unit volume (P/VL) and gas sparging rate expressed as the superficial velocity (Vs). The power input per unit volume (P/VL) and, superficial gas velocity Vs are the major correlation coefficients for kLa. Therefore, equations of the following type are frequently found in the literature (Vant Riet, 1979).
where, P, is the mechanical agitation power in gas liquid dispersion (W); VL, liquid volume (m3); Vs, gas superficial velocity (m s-1); A, constant; α and β are exponents. When P is omitted an equation like the one following can usually be adapted.
kLa = ANDα1,α2 [Vs]β 
Where, D is the impeller diameter (m) and N is the impeller speed. However, since the accurate measurement of the actual power input into the system may be quite cumbersome and many investigators (Wu, 1995) have chosen to use the Power factor term instead, i.e. N3D2 to represent the power input per volume. In this present work, kLa values obtained from the experimental data were plotted against the operating variables (stirrer speed and superficial gas velocity) and the correlation was obtained for power number (N3D2) by performing multiple regression analysis in MATLAB (Version 7.0).
Attempts were made to correlate the heat production rate (q) with volumetric mass transfer coefficient (kLa) and a linear correlation was obtained at optimized process conditions.
Figure 3 shows the plot of kLa values for these 2 types of stirrer (Turbine and Rushton-turbine type). It was clear from the Figure 3 that turbine-impeller type stirrer shows a maximum kLa value at all the growth phases of the culture and was fixed as a suitable impeller for oxygen mass transfer experiments. Oxygen mass transfer in Rushton-turbine type impeller was found to decrease with the biomass growth. This may be attributed to the increase in viscosity of broth media due to biomass growth and the inefficiency of Rushton-turbine type to break the bubbles effectively for enhanced mass transfer. Further, the high shear force caused a negative effect on growth of culture and made it unsuitable for further studies. Excessive foaming was also observed when Rushton-Turbine type impeller was used.
The effects of impeller speed (turbine) on the kLa was investigated using the impeller speeds in the range of 50-200 rpm by maintaining constant working volume of 0.75 l, oxygen flow rate of 1 l min-1, temperature of 37ºC and pH 7.0. Experiments were performed at respective impeller speeds starting from lag-phase until the culture reached its stationary phase. The values of kLa for different impeller speeds are plotted in Figure 4. The volumetric mass transfer coefficient, kLa, increased from 0.1212 to 0.582 min-1 when the impeller speed was increased from 50 to 200 rpm. It can be also seen from the Figure 5 that further increase in impeller speed to 250 rpm showed no significant rise in kLa value (results not shown). Hence an optimum stirrer speed of 200 rpm was fixed for efficient oxygen mass transfer for the cultivation of P. aeruginosa. Analysis of biomass growth of P. aeruginosa at different impeller speeds also proved that 200 rpm is an optimum value (Figure 5). Further increase in impeller speed to 250 rpm caused a significant decrease in growth, perhaps due to inefficient oxygen transfer under high turbulence conditions. Further we also observed high foaming formation at 250 rpm. In order to prove the decrease in biomass growth rates at high stirrer speeds, the specific growth rates (μ) were estimated for all the evaluated conditions and shown in Figure 6. It can be seen that at 250 rpm the specific growth rate decreases and supports our decision for fixing the stirrer speed to 200 rpm.
The effects of oxygen flow rates on kLa were investigated and the values of kLa were plotted in Figure 7. Volumetric mass transfer coefficient, kLa, increased from 0.181 to 0.41 min-1 when the oxygen flow rate was increased from 0.5 to 1 l min-1. Further increase in oxygen flow to 1.5 l min-1 show no improvement in mass transfer (results not shown). Cell dry weight was estimated at regular time intervals and specific OUR (qO2, min-1) was calculated. Figure 8 shows variation inqO2profiles with change in oxygen flow rates. From Figure 8, it was found thatqO2results also showing similar trend like kLa plot, where 1 l min-1 was found to be optimum value.
Metabolic heat generated due to the growth of P. aeruginosa in NB media at varying impeller speeds and oxygen flow rates were simultaneously measured in biocalorimeter and plotted against time. Figure 9 shows a typical power-time curve obtained for the growth of P. aeruginosa in NB media at optimized conditions (1 l min-1 oxygen flow rate and 200 rpm impeller speed). From the Figure 9, it was observed that the biomass growth profile was following closely the power-time curve and this clearly shows that growth of P. aeruginosa is directly proportional to metabolic heat generation. BothqO2and heat values seem to be influenced by growth phases of P. aeruginosa.qO2, start to increase from lag phase and dramatically increases in exponential phase to a maximum value and then starts falling in endogenous phase (death phase). The same trend was observed in heat values during the growth of P. aeruginosa in NB media at each different operating condition. From Figure 9, the maximum value ofqO2was found to be 2.888 min-1 at the middle of the exponential phase and the max. heat generated also occurred at the same phase (3.01 W). Usually both the respirometric and calorimetric data provides similar information (Voisard et al. 2002) and was in good agreement with the experimental results. By analyzing the oxygen mass transfer and calorimetric results, the optimum impeller speed and oxygen flow rate were determined for the growth of P. aeruginosa.
Total metabolic heat generated (cumulative) till the exponential growth phase at varying impeller speed conditions were evaluated using WINRC software (Mettler-Toledo, AG, Switzerland) and compared in the plot shown in Figure 10. Heat added due to stirring power was systematically eliminated from the cumulative metabolic heat generated by performing baseline experiments (results not given). Figure 10 clearly depicts that the metabolic heat production was increased approximately 4-fold when the impeller speed was raised from 50 to 200 rpm. It can be also noted from the Figure 10, that further increase in impeller speed (to 250 rpm) caused a considerable fall in heat production. In the similar way, heat productions at varying oxygen flow rates were analyzed and results were shown in Figure 11. From the Figure 11, it can be seen that increasing the oxygen flow from 0.5 l min-1 to 1 l min-1 caused a two fold increase in heat generation. Increasing the oxygen flow to 1.5 l min-1 resulted an observable fall in heat generation. Hence 1 l min-1 oxygen flow rate was found to be optimum for the better growth of P. aeruginosa in NB media. Thus from the calorimetric studies it can be stated that the optimum conditions for growth of P. aeruginosa are at 200 rpm (impeller speed) and 1 l min-1 (oxygen flow) respectively. The kLa value obtained under the above conditions was found to be maximum (0.559 min-1).
It can be found from experimental data that kLa depends on impeller speed and attains a maximum value at 200 rpm. 1 l min-1 oxygen flow showed a maximum kLa value when compared to other flow rates.
In order to determine the influence of operational conditions (Vs, N) on volumetric mass transfer coefficient values have been considered and kLa values have been correlated according to the relationship given below Equation 8,
kLaα (N3D2)α (Vs)β 
In the above equation, the influence and the relative importance of the considered variables are determined by the coefficients α and β values. The values of these coefficients are specific for each microorganism type or morphology. The exponential relationship between kLa and impeller speed (N) for this study is given by Equation (9),
kLaα (N3D2)0.24 
The statistical values of R2 and standard error (σ) obtained by using the data were 0.9997 and 0.00016 respectively. In the absence of power measurement equipment, the exponents of N3D2 were obtained as representation of the power number per unit volume . The exponent of N3D2 obtained from this study (α = 0.24), is within the range reported in the literature (0.16-0.68). The wide variation of α is related to the geometry of the vessel (the proportionality of impeller diameter to vessel diameter, tank bottom geometry, tank liquid height, etc.) employed by different researchers.
The exponential relationship between kLa and superficial oxygen velocity (Vs) for this study is given by Equation (10),
kLaα (Vs)0.41 
The statistical values of R2 and standard error (σ) obtained by using the data were 0.9995 and 0.00029 respectively. The exponent (β = 0.41) representing the dependence of kLa on the superficial gas velocity (Vs) was compared with the literature values. Previous investigators (Hyman and Bogaerde, 1960) obtained the values of the exponent on the superficial gas velocity in the range 0.17 to 0.67, depending on the agitator speed and the geometry of the equipment. The exponent obtained from this study again well corroborates to the reported values. A satisfactory linear correlation (Figure 12) was observed between volumetric mass transfer coefficient (kLa) and metabolic heat (q) at optimized process conditions (200 rpm & 1 l min-1) as follows
kLa = 0.0268 (q) + 0.3176 
Evaluation of the experimental data showed that kLa values are affected by process variables such as impeller speed, oxygen flow rate and geometry of the stirrer. Increase in impeller speed resulted in four-fold increase in kLa value and the increase in oxygen flow rate caused a two-fold rise in kLa values. Turbine-type impeller showed maximum oxygen transfer compared to Rushton-turbine type. The kLa values obtained from the data of oxygen absorption were used to establish correlations as a function of the process variables. The exponents of power number (N3D2) and superficial gas velocity (Vs) found for NB growth media were 0.24 and 0.41 respectively, are in good agreement with the exponents given in the literature. Calorimetric experiments revealed that the growth of P. aeruginosa in NB media was found to follow the heat generated due to metabolic activity. BothqO2and metabolic heat are found to be influenced greatly by growth phases of P. aeruginosa and attained their maximum value at the middle of the exponential phase. Heat data (q) were found to correlate well with volumetric mass transfer coefficient (kLa) and a linear correlation was obtained at optimized process conditions. The correlations obtained in this study clearly indicates that both respirometric and calorimetric results from biocalorimetric experiments could be used for effective scale-up and design of biological reactions to large scale.
One of the authors Mr. S. Senthilkumar is grateful to CSIR, New Delhi for the SRF fellowship. The author thanks the Director, CLRI for his kind permission to publish this work. The author also wishes to express his gratitude to Prof. K.S. Gandhi for his useful suggestion and Prof. N.R. Rajagopal for his constant encouragement.
CALIK, Güzide; VURAL, Halil and OZDAMAR, Tunçer. Bioprocess parameters and oxygen transfer effects in the growth of Pseudomonas dacunhae for L-alanine production. Chemical Engineering Journal, May 1997, vol. 65, no. 2, p. 109-116. [ [ Links ]CrossRef]
CHISTI, Yusuf and JAUREGUI-HAZA, Ulises J. Oxygen transfer and mixing in mechanically agitated airlift bioreactors. Biochemical Engineering Journal, March 2002, vol. 10, no. 2, p. 143-153. [ [ Links ]CrossRef]
GALACTION, A.-I.; CASCAVAL, D.; ONISCU, C. and TURNEA, M. Prediction of oxygen mass transfer coefficients in stirred bioreactors for bacteria, yeasts and fungus broths. Biochemical Engineering Journal, August 2004, vol. 20, no. 1, p. 85-94. [ [ Links ]CrossRef]
GARCIA-OCHOA, Félix and GOMEZ-CASTRO, Emilio. Mass transfer coefficient in stirrer tank reactors for xanthan gum solutions. Biochemical Engineering Journal, January 1998, vol. 1, no. 1, p. 1-10. [ [ Links ]CrossRef]
GARCIA-OCHOA, F.; GOMEZ-CASTRO, E. and SANTOS, V.E. Oxygen transfer and uptake rates during xanthan gum production. Enzyme and Microbial Technology, November 2000, vol. 27, no. 9, p. 680-690. [ [ Links ]CrossRef]
GARCIA-OCHOA, Felix and GOMEZ-CASTRO, Emilio. Prediction of gas-liquid mass transfer coefficient in sparged stirred tank bioreactors. Biotechnology and Bioengineering, December 2005, vol. 92, no. 6, p. 761-772. [ [ Links ]CrossRef]
HARIRI, M. Hossein; LEWIS, Robert A.; LONIS, Timothy R.; PARKER, Bruce E. and STEPHENSON, Kay A. The effect of operating conditions in a reaction calorimeter. Thermochimica Acta, December 1996, vol. 289, no. 2, p. 343-349. [ [ Links ]CrossRef]
HUGH, Rudolph and LEIFSON, Einar. The taxonomic significance of fermentative versus oxidative metabolism of carbohydrates by various gram-negative bacteria. Journal of Bacteriology, January 1953, vol. 66, no. 1, p. 24-26. [ Links ]
MARISON, I.W.; LIU, J.S.; AMPUERO, S.; VON STOCKAR, U. and SCHENKER, U. Biological reaction calorimetry: Development of high sensitivity bio-calorimeters. Thermochimica Acta, January 1998, vol. 309, no. 1, p. 157-173. [ [ Links ]CrossRef]
NISHIKAWA, Masabumi; NAKAMURA, Masao and HASHIMOTO, Kenji. Gas absorption in aerated mixing vessels with non-Newtonian liquid. Journal of Chemical Engineering of Japan, June 1981, vol. 14, no. 3, p. 227-232. [ Links ]
VANT RIET, Klass. Review of measuring methods and results in nonviscous gas-liquid mass transfer in stirred vessels. Industrial and Engineering Chemistry Process Design and Development, 1979, vol. 18, no. 3, p. 357-364. [ [ Links ]CrossRef]
VOISARD, D.; VON STOCKAR, U. and MARISON, I.W. Quantitative calorimetric investigation of fed-batch cultures of Bacillus sphaericus 1593M. Thermochimica Acta, October 2002, vol. 394, no. 1-2, p. 99-111. [ [ Links ]CrossRef]
Note: Electronic Journal of Biotechnology is not responsible if on-line references cited on manuscripts are not available any more after the date of publication.