Inter-compartment interaction in multi-impeller mixing. Part II. Experiments, sliding mesh and large Eddy simulations

An Experimental and Computational Approach to Development of Mixing Processes for Miscible Liquids in Unbaffled Tanks

PURPOSE

Mixing of liquids is a critical unit operation in the biopharmaceutical drug product manufacturing. It commonly consists of mixing miscible liquids to dilute bulk drug substance (DS) or pool multiple lots of drug substance. In the past, at-scale mixing studies have been conducted to determine the mixing parameters, namely mixing speed, and mixing time. At-scale studies have historically been utilized to overcome the challenges associated with geometric dissimilarity of mixing systems found when scaling up. In addition, such studies are quite costly, as they often use actual DS to overcome a lack of representativeness associated with simple salt trace models often employed. As a result, there is a significant need for alternative cost-effective methods that can predict mixing parameters with close agreement to actual experiments and operations.

METHOD

At-scale mixing experiments were conducted using full-sized tanks and surrogate solutions. Several computational fluid dynamic (CFD) simulation methods were conducted and compared with the experiments to determine the most reliable computational techniques.

RESULTS

The experiments demonstrate that surrogate solutions can be used reliably to determine mixing parameters in at-scale studies instead of the valuable drug products. Studying different CFD methods also showed that transient simulations that use a large eddy simulation (LES) viscous model and a sliding mesh can correctly predict the mixing parameters.

CONCLUSION

Results of this study establish a practical and reliable methodology to perform mixing studies for miscible liquids with different kinematic viscosities. The methods discussed herein greatly reduce the routine mixing study costs in the biopharmaceutical industry and increase efficiency and accuracy of the results, allowing proactive scale-up of mixing operations.

An analysis of organism lifelines in an industrial bioreactor using Lattice-Boltzmann CFD

Euler-Lagrange CFD simulations, where the biotic phase is represented by computational particles (parcels), provide information on environmental gradients inside bioreactors from the microbial perspective. Such information is highly relevant for reactor scale-down and process optimization. One of the major challenges is the computational intensity of CFD simulations, especially when resolution of dynamics in the flowfield is required. Lattice-Boltzmann large-eddy simulations (LB-LES) form a very promising approach for simulating accurate, dynamic flowfields in stirred reactors, at strongly reduced computation times compared to finite volume approaches. In this work, the performance of LB-LES in resolving substrate gradients in large-scale bioreactors is explored, combined with the inclusion of a Lagrangian biotic phase to provide the microbial perspective. In addition, the hydrodynamic performance of the simulations is confirmed by verification of hydrodynamic characteristics (radial velocity, turbulent kinetic energy, energy dissipation) in the impeller discharge stream of a 29 cm diameter stirred tank. The results are compared with prior finite volume simulation results, both in terms of hydrodynamic and biokinetic observations, and time requirements.

Investigating a Stirred Bioreactor: Impact of Evolving Fermentation Broth Pseudoplastic Rheology on Mixing Mechanisms

The culture medium in many fermentations is a non-Newtonian fluid. In bacterial alginate batch production, the broth becomes more pseudoplastic as the alginate concentration increases, which impairs the mixing process. This work characterizes the effect of the interaction between changing broth rheology and impeller mixing on a bioreactor fluid dynamics. Experimentally, a fermentation with evolving broth pseudoplastic rheology is reproduced. Three fermentation stages are mimicked using appropriate solutions of water and xanthan gum. Impeller torque measurements are reported. The weakening of the impellers’ interaction over the fermentation process is identified. To overcome the experimental limitations, CFD is applied to study the evolution of the fermentation fluid flow patterns, velocity field, dead zones, and vortical structures. Precessional vortex macro-instabilities are identified as being responsible for the unstable flow patterns identified at the earlier stages of the fermentation. A stable parallel flow pattern accounts for the weakest impellers’ interaction at the final stage. Overall, this work contributes with a complete workflow to adapt CFD models for characterization and aided design of stirred tanks with changing broth pseudoplastic rheology as well as an evolving flow regime.

Validation of Novel Lattice Boltzmann Large Eddy Simulations (LB LES) for Equipment Characterization in Biopharma

Detailed process and equipment knowledge is crucial for the successful production of biopharmaceuticals. An essential part is the characterization of equipment for which Computational Fluid Dynamics (CFD) is an important tool. While the steady, Reynolds-averaged Navier–Stokes (RANS) k − ε approach has been extensively reviewed in the literature and may be used for fast equipment characterization in terms of power number determination, transient schemes have to be further investigated and validated to gain more detailed insights into flow patterns because they are the method of choice for mixing time simulations. Due to the availability of commercial solvers, such as M-Star CFD, Lattice Boltzmann simulations have recently become popular in the industry, as they are easy to set up and require relatively low computing power. However, extensive validation studies for transient Lattice Boltzmann Large Eddy Simulations (LB LES) are still missing. In this study, transient LB LES were applied to simulate a 3 L bioreactor system. The results were compared to novel 4D particle tracking (4D PTV) experiments, which resolve the motion of thousands of passive tracer particles on their journey through the bioreactor. Steady simulations for the determination of the power number followed a structured workflow, including grid studies and rotating reference frame volume studies, resulting in high prediction accuracy with less than 11% deviation, compared to experimental data. Likewise, deviations for the transient simulations were less than 10% after computational demand was reduced as a result of prior grid studies. The time averaged flow fields from LB LES were in good accordance with the novel 4D PTV data. Moreover, 4D PTV data enabled the validation of transient flow structures by analyzing Lagrangian particle trajectories. This enables a more detailed determination of mixing times and mass transfer as well as local exposure times of local velocity and shear stress peaks. For the purpose of standardization of common industry CFD models, steady RANS simulations for the 3 L vessel were included in this study as well.

Euler-Lagrangian Simulations: A Proper Tool for Predicting Cellular Performance in Industrial Scale Bioreactors

Eulerian-Lagrangian approach to investigate cellular responses in a bioreactor has become the center of attention in recent years. It was introduced to biotechnological processes about two decades ago, but within the last few years, it proved itself as a powerful tool to address scale-up and -down topics of bioprocesses. It can capture the history of a cell and reveal invaluable information for, not only, bioprocess control and design but also strain engineering. This way it will be possible to shed light on the actual environment that cell experiences throughout its lifespan. Lifelines of a microorganism in a bioreactor can serve as the missing link that encompasses the biological timescales and the physical timescales. For this purpose digitalization of bioreactors provides us with new insights that are not achievable in industrial reactors easily if at all, namely, substrate and product gradients; high-shear regions are among the most interesting factors that can be reproduced adequately with help of a digital twin. In this chapter basic principles of this method will be introduced, and later on some practical aspects of particle tracking technique will be illustrated. In the final section, some of the advantages and challenges associated with this method will be discussed.

Coupled metabolic-hydrodynamic modeling enabling rational scale-up of industrial bioprocesses

Metabolomics aims to address what and how regulatory mechanisms are coordinated to achieve flux optimality, different metabolic objectives as well as appropriate adaptations to dynamic nutrient availability. Recent decades have witnessed that the integration of metabolomics and fluxomics within the goal of synthetic biology has arrived at generating the desired bioproducts with improved bioconversion efficiency. Absolute metabolite quantification by isotope dilution mass spectrometry represents a functional readout of cellular biochemistry and contributes to the establishment of metabolic (structured) models required in systems metabolic engineering. In industrial practices, population heterogeneity arising from fluctuating nutrient availability frequently leads to performance losses, that is reduced commercial metrics (titer, rate, and yield). Hence, the development of more stable producers and more predictable bioprocesses can benefit from a quantitative understanding of spatial and temporal cell-to-cell heterogeneity within industrial bioprocesses. Quantitative metabolomics analysis and metabolic modeling applied in computational fluid dynamics (CFD)-assisted scale-down simulators that mimic industrial heterogeneity such as fluctuations in nutrients, dissolved gases, and other stresses can procure informative clues for coping with issues during bioprocessing scale-up. In previous studies, only limited insights into the hydrodynamic conditions inside the industrial-scale bioreactor have been obtained, which makes case-by-case scale-up far from straightforward. Tracking the flow paths of cells circulating in large-scale bioreactors is a highly valuable tool for evaluating cellular performance in production tanks. The "lifelines" or "trajectories" of cells in industrial-scale bioreactors can be captured using Euler-Lagrange CFD simulation. This novel methodology can be further coupled with metabolic (structured) models to provide not only a statistical analysis of cell lifelines triggered by the environmental fluctuations but also a global assessment of the metabolic response to heterogeneity inside an industrial bioreactor. For the future, the industrial design should be dependent on the computational framework, and this integration work will allow bioprocess scale-up to the industrial scale with an end in mind.