CFD simulation of mixing in tall gas–liquid stirred vessel: Role of local flow patterns

Validation of a CFD model for cell culture bioreactors at large scale and its application in scale-up

Among all the operating parameters that control the cell culture environment inside bioreactors, appropriate mixing and aeration are crucial to ensure sufficient oxygen supply, homogeneous mixing, and CO2 stripping. A model-based manufacturing facility fit approach was applied to define agitation and bottom air flow rates during the process scale-up from laboratory to manufacturing, of which computational fluid dynamics (CFD) was the core modeling tool. The realizable k-ε turbulent dispersed Eulerian gas-liquid flow model was established and validated using experimental values for the volumetric oxygen transfer coefficient (kLa). Model validation defined the process operating parameter ranges for application of the model, identified mixing issues (e.g., impeller flooding, dissolved oxygen gradients, etc.) and the impact of antifoam on kLa. Using the CFD simulation results as inputs to the models for oxygen demand, gas entrance velocity, and CO2 stripping aided in the design of the agitation and bottom air flow rates needed to meet cellular oxygen demand, control CO2 levels, mitigate risks for cell damage due to shear, foaming, as well as fire hazards due to high O2 levels in the bioreactor gas outlet. The recommended operating conditions led to the completion of five manufacturing runs with a 100% success rate. This model-based approach achieved a seamless scale-up and reduced the required number of at-scale development batches, resulting in cost and time savings of a cell culture commercialization process.

CFD evaluation of hydrophobic feedstock bench-scale fermenters for efficient high agitation volumetric mass transfer

A new biomanufacturing platform combining intracellular metabolic engineering of the oleaginous yeast Yarrowia lipolytica and extracellular bioreaction engineering provides efficient bioconversion of plant oils/animal fats into high-value products. However, predicting the hydrodynamics and mass transfer parameters is difficult due to the high agitation and sparging required to create dispersed oil droplets in an aqueous medium for efficient yeast fermentation. In the current study, commercial computational fluid dynamic (CFD) solver Ansys CFX coupled with the MUSIG model first predicts two-phase system (oil/water and air/water) mixing dynamics and their particle size distributions. Then, a three-phase model (oil, air, and water) utilizing dispersed air bubbles and a polydispersed oil phase was implemented to explore fermenter mixing, gas dispersion efficiency, and volumetric mass transfer coefficient estimations (kL a). The study analyzed the effect of the impeller type, agitation speed, and power input on the tank's flow field and revealed that upward-pumping pitched blade impellers (PBI) in the top two positions (compared to Rushton-type) provided advantageous oil phase homogeneity and similar estimated kL a values with reduced power. These results show good agreement with the experimental mixing and kL a data.

Gas Dispersion in Non-Newtonian Fluids with Mechanically Agitated Systems: A Review

Gas dispersion in non-Newtonian fluids is encountered in a broad range of chemical, biochemical, and food industries. Mechanically agitated vessels are commonly employed in these processes because they promote high degree of contact between the phases. However, mixing non-Newtonian fluids is a challenging task that requires comprehensive knowledge of the mixing flow to accurately design stirred vessels. Therefore, this review presents the developments accomplished by researchers in this field. The present work describes mixing and mass transfer variables, namely volumetric mass transfer coefficient, power consumption, gas holdup, bubble diameter, and cavern size. It presents empirical correlations for the mixing variables and discusses the effects of operating and design parameters on the mixing and mass transfer process. Furthermore, this paper demonstrates the advantages of employing computational fluid dynamics tools to shed light on the hydrodynamics of this complex flow. The literature review shows that knowledge gaps remain for gas dispersion in yield stress fluids and non-Newtonian fluids with viscoelastic effects. In addition, comprehensive studies accounting for the scale-up of these mixing processes still need to be accomplished. Hence, further investigation of the flow patterns under different process and design conditions are valuable to have an appropriate insight into this complex system.

CFD simulations of stirred-tank reactors for gas-liquid and gas-liquid-solid systems using OpenFOAM®

An open-source CFD software OpenFOAM® is used to simulate two multiphase stirred-tank reactors relevant to industrial processes such as slurry polymerization and fuel production. Gas-liquid simulations are first performed in a single-impeller stirred-tank reactor, studied experimentally by Ford, J. J., T. J. Heindel, T. C. Jensen, and J. B. Drake. 2008. “X-Ray Computed Tomography of a Gas-Sparged Stirred-Tank Reactor.” Chemical Engineering Science 63: 2075–85. Three impeller rotation speeds (200, 350 and 700 rpm) with three different bubble diameters (0.5, 1.5 and 2.5 mm) are investigated. Flow patterns compared qualitatively to those from experiments. Compared to the experimental data, the simulations are in relatively good agreement for gas holdup in the reactor. The second multiphase system is a multi-impeller stirred-tank reactor, studied experimentally by Shewale, S. D., and A. B. Pandit. 2006. “Studies in Multiple Impeller Agitated Gas-Liquid Contractors.” Chemical Engineering Science 61: 486–504. Gas-liquid simulations are performed at two impeller rotation speeds (3.75 and 5.08 RPS). The simulated flow patterns agree with published pictures from the experiments. Gas-liquid-solid simulations of the multi-impeller stirred-tank reactor are also carried out at impeller rotation speed 5.08 RPS. The addition of solid particles with a volume fraction characteristic of slurry reactors changes the flow pattern significantly. The bottom Rushton turbine becomes flooded, while the upper pitched-blade downflow turbines present a radial-pumping flow pattern instead of down-pumping. Nonetheless, the solid phase has a similar flow pattern to the liquid phase, indicating that the particles modify the effective density of the fluid.

A General Review of the Current Development of Mechanically Agitated Vessels

The mixing process in a mechanically agitated vessel is a widespread phenomenon which plays an important role among industrial processes. In that process, one of the crucial parameters, the mixing efficiency, depends on a large number of geometrical factors, as well as process parameters and complex interactions between the phases which are still not well understood. In the last decade, large progress has been made in optimisation, construction and numerical and experimental analysis of mechanically agitated vessels. In this review, the current state in this field has been presented. It shows that advanced computational fluid dynamic techniques for multiphase flow analysis with reactions and modern experimental techniques can be used with success to analyse in detail mixing features in liquid-liquid, gas-liquid, solid-liquid and in more than two-phase flows. The objective is to show the most important research recently carried out.

A Comparative CFD Study on Gas-Liquid Dispersion in A Stirred Tank with Low and High Gas Loadings

The computational fluid dynamics (CFD) combined with a population balance model (PBM) was applied to simulate gas-liquid dispersion in a stirred tank with low and high gas loadings. The model predictions were validated by using the data in the literature. The simulation results show that the flow patterns and gas dispersion characteristics are very different in the stirred tank for low and high gas loadings. A typical two-loop flow pattern forms as that in single-phase stirred tank for low gas loadings, while a triple-loop flow pattern, with two recirculation loops above and one below the impeller is found in the tank for high gas loadings. Shaft power input of impeller agitation plays a major role for gas dispersion with low gas loadings. For high gas loadings, the potential energy due to gas sparging has significant effect on gas dispersion and can not be neglected. Compared to low gas loading, high gas loading causes average gas holdup increased in the stirred tank, while relative local gas holdup in the lower circulation-loop region and near-wall region reduced. The ability of impeller agitation for gas dispersion reduces with high gas loadings, and mean bubble size becomes larger and the volume-averaged bubble size distribution is wider.

Hydrodynamic performance of a single-use aerated stirred bioreactor in animal cell culture: applications of tomography, dynamic gas disengagement (DGD), and CFD

The hydrodynamics of gas-liquid two-phase flow in a single-use bioreactor were investigated in detail both experimentally and numerically. Electrical resistance tomography (ERT) and dynamic gas disengagement (DGD) combined with computational fluid dynamics (CFD) were employed to assess the effect of the volumetric gas flow rate and impeller speed on the gas-liquid flow field, local and global gas holdup values, and Sauter mean bubble diameter. From the results obtained from DGD coupled with ERT, the bubble sizes were determined. The experimental data indicated that the total gas holdup values increased with increasing both the rotational speed of impeller and volumetric gas flow rate. Moreover, the analysis of the flow field generated inside the aerated stirred bioreactor was conducted using CFD results. Overall, a more uniform distribution of the gas holdup was obtained at impeller speeds ≥ 100 rpm for volumetric gas flow rates ≥ 1.6 × 10^-5 m³/s.

CFD Simulation of Gas Dispersion in a Stirred Tank of Dual Rushton Turbines

3D Eulerian-Eulerian model was applied to simulate the gas-liquid two-phase flow in a stirred tank of dual Rushton turbines using computational fluid dynamics (CFD). The effects of two different bubble treatment methods (constant bubble sizevs. population balance model, PBM) and two different coalescence models (Luo modelvs. Zaichik model) on the prediction of liquid flow field, local gas holdup or bubble size distribution were studied. The results indicate that there is less difference between the predictions of liquid flow field and gas holdup using the above models, and the use of PBM did not show any advantage over the constant bubble size model under lower gas holdup. However, bubble treatment methods have great influence on the local gas holdup under larger gas flow rate. All the models could reasonably predict the gas holdup distribution in the tank operated at a low aeration rate except the region far from the shaft. Different coalescence models have great influence on the prediction of bubble size distribution (BSD). Both the Luo model and Zaichik model could qualitatively estimate the BSD, showing the turning points near the impellers along the height, but the quantitative agreement with experiments is not achieved. The former over-predicts the BSD and the latter under-predicts, showing that the existing PBM models need to be further developed to incorporate more physics.

SIMULATION OF SOLID–LIQUID SUSPENSION AND SCALE-UP OF AGAROSE GEL ACTIVATION REACTOR

Agarose gel activation reaction, which is of great importance in preparing the carrier of the column packing material for blood purification, would be significantly influenced by the configuration and parameter of reactor. In order to optimize the structure design of the reactor and operating parameters of the process, the characteristics of suspension system composed of agarose gel, NaOH, water, 3-allyl bromide, and activated alumina were simulated numerically utilizing an Eulerian multiphase flow model and multi-reference frame (MRF) approach. The effect of impeller configuration was studied with three typical impellers, including Rushton disk turbine (DT), pitched blade downflow turbine (PBTD45), and pitched blade upflow turbine (PBTU45). The results showed that the optimum solid suspension was obtained using PBTD45 impeller with a diameter of 100[Formula: see text]mm and critical suspension speed of 570[Formula: see text]rpm in a 20[Formula: see text]L stirred reactor. In addition, the critical suspension speeds using the three impellers were calculated and the errors were all within the engineering allowance. Finally, the feasibility of scale-up design for agarose gel activation reactor was preliminarily discussed. The results would be useful to the optimization and scale-up of relevant reactor design.