Next Article in Journal
Temporal Evolution of Cooling by Natural Convection in an Enclosed Magma Chamber
Next Article in Special Issue
Modern Sensor Tools and Techniques for Monitoring, Controlling, and Improving Cell Culture Processes
Previous Article in Journal
Research on Segmented Belt Acceleration Curve Based on Automated Mechanical Transmission
Previous Article in Special Issue
On the Use of Surface Plasmon Resonance-Based Biosensors for Advanced Bioprocess Monitoring
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Insights from Locally Resolved Hydrodynamics in Stirred Cell Culture Reactors

1
Institute of Bioprocess and Biosystems Engineering, Hamburg University of Technology, Denickestr. 15, 21073 Hamburg, Germany
2
Institute of Multiphase Flows, Hamburg University of Technology, Eißendorfer Str. 38, 21073 Hamburg, Germany
*
Author to whom correspondence should be addressed.
Processes 2022, 10(1), 107; https://doi.org/10.3390/pr10010107
Submission received: 11 November 2021 / Revised: 11 December 2021 / Accepted: 22 December 2021 / Published: 5 January 2022

Abstract

:
The link between hydrodynamics and biological process behavior of antibody-producing mammalian cell cultures is still not fully understood. Common methods to describe dependencies refer mostly to averaged hydrodynamic parameters obtained for individual cultivation systems. In this study, cellular effects and locally resolved hydrodynamics were investigated for impellers with different spatial hydrodynamics. Therefore, the hydrodynamics, mainly flow velocity, shear rate and power input, in a single- and a three-impeller bioreactor setup were analyzed by means of CFD simulations, and cultivation experiments with antibody-producing Chinese hamster ovary (CHO) cells were performed at various agitation rates in both reactor setups. Within the three-impeller bioreactor setup, cells could be cultivated successfully at much higher agitation rates as in the single-impeller bioreactor, probably due to a more uniform flow pattern. It could be shown that this different behavior cannot be linked to parameters commonly used to describe shear effects on cells such as the mean energy dissipation rate or the Kolmogorov length scale, even if this concept is extended by locally resolved hydrodynamic parameters. Alternatively, the hydrodynamic heterogeneity was statistically quantified by means of variance coefficients of the hydrodynamic parameters fluid velocity, shear rate, and energy dissipation rate. The calculated variance coefficients of all hydrodynamic parameters were higher in the setup with three impellers than in the single impeller setup, which might explain the rather stable process behavior in multiple impeller systems due to the reduced hydrodynamic heterogeneity. Such comprehensive insights lead to a deeper understanding of the bioprocess.

1. Introduction

Mammalian cell culture processes are state-of-the-art for the production of therapeutic antibodies. However, the influence of the bioreactor hydrodynamics on the cell culture process are still not fully understood [1,2,3,4,5,6]. Therefore, bioreactor design and scale-up in today’s biopharma industry rely mostly on empirical correlations, experience, and engineering heuristics. Common methods for scaling-up mammalian cell-based production processes aim to keep the reactor geometry and certain process parameters such as the average volumetric power input constant [7,8,9,10,11].
These “rules of thumb” methods are limited as they can hardly be used to describe or estimate an appropriate operation range with respect to shear effects on cellular behavior. Platas et al. [12] observed an approx. constant cell specific growth rate, µmax, over a broader stirrer operation range in reactor systems with multiple impellers compared to reactor systems with only one impeller. The cell growth rate in reactors with multiple impellers decreased much slower with increasing agitation rates and corresponding mean power inputs than in reactors with less impellers. Therefore, shear effects seem to be less pronounced in the case of multiple impellers.
An approach frequently discussed in the literature to predict critical conditions with respect to shear effects to be expected in stirred tank bioreactors is the Kolmogorov eddy length scale [12,13], for which cell harm is predicted if the eddies are in the order of magnitude as the cell diameter. This approach should predict an acceptable power input—and thus an adequate agitation rate—to be estimated without extensive experiments. Conventionally, the system averaged energy dissipation rate is used to calculate the Kolmogorov length scale, which does not consider local gradients, even if the local power input can differ by several orders of magnitude within the reactor. Furthermore, this length scale is usually compared considering an average cell diameter, even though the cell diameter is widely distributed. However, while this hypothesis seems to be proven for microcarrier cultures, it is still under discussion for suspension cells [14,15,16,17,18,19]. A broad variety of studies on lethal and sub-lethal responses of mammalian cell cultures to shear stress has been published where a detailed overview can be found in Chalmers et al. The investigated magnitude of the average volumetric power input ranges from 101 to almost 109 W m−3 [3]. However, the vast majority of these studies only considered the estimated or averaged hydrodynamics.
Taken together, more insights into locally resolved hydrodynamics are needed. Since local maxima can reach much higher values than the averaged, locally resolved hydrodynamics could provide more detailed information on the hydrodynamics. These can be identified and quantified by the use of computational fluid dynamics (CFD) methods, which have become an established tool for cell culture process development [5,10,20,21,22]. CFD can be used to predict the locally resolved hydrodynamic behavior of cell cultivation systems and has already been successfully applied for the design of stem cell expansion processes in stirred tank reactors and wave bags with microcarriers [23,24,25].
In this study, we investigated to what extent the cellular and hydrodynamic effects change with different spatial hydrodynamics of different stirrers. Furthermore, it was discussed whether locally resolved hydrodynamics could help to explain the cellular effects. Therefore, the hydrodynamics, mainly flow velocity, shear rate, and power input, in a single- and a three-impeller bioreactor setup were analyzed by means of CFD simulations. Then, cultivation experiments with antibody-producing Chinese hamster ovary (CHO) cells were performed at various agitation rates in both reactor setups. Finally, conventional process design parameters such as the average volumetric power input and the Kolmogorov length scale were evaluated by means of the obtained data. As the model system, a reactor setup was chosen that could be equipped with various numbers of impellers and operated without baffles at constant culture volumes. Because surface aeration is sufficient for a wide range of operation conditions, the impact of bubble aeration was neglected in these investigations.

2. Materials and Methods

2.1. Cultivation Procedures

The cell line used in all cultivation experiments was CHO DP-12, which was kindly provided by Prof. Dr. Thomas Noll (University of Bielefeld, Bielefeld, Germany). All cultivation experiments were performed in TC-42 medium (Xell AG, Bielefeld, Germany), supplemented with 8 mmol L−1 l-glutamine and 200 nmol L−1 methotrexate.

2.1.1. Pre-Culture

Pre-cultivations were carried out in 125 mL shake flasks with 40 mL culture volume for the first culture after thawing and 250 mL shake flasks with 80 mL culture volume for expanding the cells [26]. Culture conditions were adjusted to 37 °C, 5% CO2, 85% relative humidity, and 200 rpm. From pre-cultures, 15 × 107 cells in total were harvested, centrifuged for 10 min at 300× g, and resuspended in 10 mL fresh medium for inoculating the stirred tank bioreactor with an initial cell concentration of 1 × 106 cells mL−1.

2.1.2. Main Culture

For bioreactor cultivations, the stirred tank bioreactor Vario 1000 (MDX Biotech GmbH, Nörten-Hardenberg, Germany) was used with 150 mL culture volume and headspace aeration by gas mixtures of air and CO2. Only at the end of the cultivations at cell densities of around 10 × 106 cells mL−1, the reactor was additionally bubble-aerated with pure oxygen. Thus, effects of gas bubbles and additional turbulences, caused by baffles could be neglected and observed effects could be referred to the impellers. The dissolved oxygen tension (DO) was controlled at 40% of air saturation in all cultivations. pH was controlled at 7.1 by the addition of CO2 via headspace or by adding 0.5 M Na2CO3. The temperature was set to 37 °C. To investigate the effect of hydrodynamics caused by stirring on the cell culture process, cultivation experiments were performed at different agitation rates for a reactor setup with one pitched blade impeller (MDX Biotech GmbH, Nörten-Hardenberg, Germany) (single impeller system, SIS) and a setup with two additional six blade impellers (self-made at Hamburg University of Technology, Nörten-Hardenberg, Germany) (triple impeller system, TIS) (see Appendix A Figure A1). The agitation rates chosen for the SIS were 200 to 1400 rpm in 200 rpm steps, while TIS cultivations at agitation rates of 770, 930, 1080, 1200 and 1400 rpm were run. In all experiments, the pitched blade impeller was operated in the down flow direction. Agitation rates for the TIS were chosen while following a discarded working hypothesis.

2.2. Analytical Methods

Cell densities were determined with the Z2 Particle Counter (Beckman Coulter, Brea, CA, USA). All particle distributions were recorded in triplicate. The cell viability was measured using the flow cytometer CytoFLEX (Beckman Coulter, Brea, CA, USA) after staining the cells with 1 µg mL−1 DAPI. Samples above 2 × 106 cells mL−1 were diluted ten-fold before staining. Antibody concentrations were determined with Protein-A binding sensors in the OCTET according to the manufacturer’s protocol (Pall Corporation, Port Washington, NY, USA). For a comparison of the processes, in addition to growth curves, the space time yield (STY) related to the antibody concentration was calculated as the quotient of the maximum antibody concentration and the cultivation time until the maximum concentration was reached.

2.3. CFD Simulations

All CFD models were developed in COMSOL Multiphysics 5.3 to 5.5 (COMSOL AB, Stockholm, Sweden). The reactor geometry was fully implemented in COMSOL, except for the six blade impeller, which was imported from a CAD file. Reactor vessel and impellers were modeled as accurately as possible and necessary. Additionally, all reactor inserts, namely the impeller shaft, temperature, pH and DO probes as well as sampling tubes and aeration tubes were implemented. As the physics module, the mixture module was chosen with turbulent flow conditions. The module contains the Reynolds Averaged Navier–Stokes equation and continuity equation as governing equations. Standard parameters were used to model the turbulence with the k-ε-model. The geometry was meshed to 1.69 million mesh elements for the SIS and 1.78 million mesh elements for the TIS. To solve the stationary studies, the PARDISO solver was set up as a direct block structured solver. For further, detailed information about the CFD models, see the model reports attached as Supplementary Materials (Files S1 and S2). Due to numerical reasons, agitation rates of up to 800 rpm for the reactor setup with a single impeller and 600 rpm for the reactor setup with three impellers were simulated. Velocities and shear rates above the given agitation rates were extrapolated linearly while the energy dissipation rates were extrapolated quadratically, according to the simplified dependency given by Hu et al. [1]. The computed values from the CFD simulations and extrapolations can be found in Appendix A (Table A1 and Table A2). A mesh refinement study was performed for all shown models (not shown). No dependencies on the degree of meshing were observed.
To highlight local distributions of hydrodynamic parameters, the computed velocities, shear rates, and energy dissipation rates were averaged over horizontal cut planes along the reactor height in steps of one millimeter. The obtained mean values were then plotted in comparison to the reactor height.

2.4. Calculation of the Critical Energy Dissipation Rate

The Kolmogorov length scale represents the size of the smallest turbulence eddies in a fluid flow and can be calculated from the energy dissipation rate ε, the kinematic viscosity ν, and the density ρ of the fluid.
λ = ν 3 ρ ε 1 4
According to the hypothesis, significant cell damage occurs when the smallest turbulence eddies are within the same order of magnitude as the cell diameter dc. Thus, the equation for the Kolmogorov length was solved for the critical energy dissipation rate εkrit where the cell diameter dc is the cell diameter, which is in the same order as the Kolmogorov length λ, ν is the kinematic viscosity, and ρ is the density of the fluid.
ε krit = ν 3 ρ d c 4
Distributions of the cell diameter dc were determined alongside the cell number in the Z2 Particle Counter. More details on the relation of turbulence eddies and cell damage according to the Kolmogorov length scale hypothesis can be found in Section 3.3.

3. Results

In the following, the results from the CFD study as well as from the cell cultivation experiments are shown for the single-impeller setup (SIS) and the three-impeller setup (TIS). Cultivation results are presented for both reactor setups to demonstrate the biological responses to different hydrodynamic conditions in the different setups. Finally, it was examined to what extent the observed effects can be explained on the basis of the hydrodynamic parameters determined via CFD simulations, where mainly the fluid velocity, power input, or energy dissipation rate were considered. It was investigated whether locally resolved parameters and parameter distributions had an advantage in this respect compared to the averaged values. This is further discussed in Section 4.

3.1. Characterization of Hydrodynamics of Single- and Multiple-Impeller Setups with CFD Methods

For both bioreactor setups (SIS and TIS, respectively), the hydrodynamic parameters fluid velocity u, shear rate γ, and energy dissipation rate ε were evaluated in the CFD study. Locally resolved values of these parameters were determined along the reactor height.
For example, the hydrodynamic parameters u and ε for 400 rpm are shown in Figure 1 and Figure 2a (SIS) and Figure 2b (TIS), plotted in longitudinal section (for shear rate γ, see Figure A2). To gain a deeper numerical insight into locally resolved values, the CFD computed values for velocity, shear rate, and energy dissipation rate were averaged for horizontal slices along the reactor height. The mean values were then plotted over these slices to quantitatively visualize locally resolved values for all simulated agitation rates (slice plots in Figure 1 and Figure 2c,d).
From the plots of the fluid velocity u (see Figure 1a,b), the flow pattern in the respective reactor setups can be derived. Both setups were simulated and operated in down flow mode. In the SIS (Figure 1a), the fluid is pressed downward and sideways by the impeller at an approximately 2.5 cm reactor height. Then, it rises to the surface at the vessel walls and comes down again close to the stirrer shaft, yielding a typical axial flow pattern. The fluid velocity u reaches maximal values near the impeller blades up to 0.55 m s−1. With increasing distance from the impeller, the fluid velocity u decreases to approx. 0.05 m s−1. In the TIS (Figure 1b), the fluid in the upper parts of the reactor is pushed to the side by the six blade impellers, known as radial flow pattern at approx. 5.5 cm and 7 cm reactor height. At the vessel walls, it flows upward and back downward to the impeller close to the stirrer shaft, resulting in a typical radial flow pattern. The system averaged fluid velocity uav was about twice as high as in the SIS with averaged 0.11 m s−1.
The plots in Figure 1c,d show general high local values as well as high maxima in narrow areas close to the impeller at about a 2.5 cm reactor height, while large parts of the reactor stayed at values up to five-fold lower. In the TIS, additional peaks could be observed at 5.5 cm and 7 cm reactor height due to the added impellers. It can be seen that the velocity in both setups reached the highest values in regions close to the impellers. Nevertheless, due to the additional impellers, the velocity profile in the TIS was more homogenized than in the SIS.
The local energy dissipation rates are shown in Figure 2. In the SIS, the energy dissipation rate εmax reached values up to 139 kW m−3 directly at the impeller (see Figure 2a). The slice plots (see Figure 2c,d) revealed that locally resolved values for the energy dissipation rate showed quite sharp and narrow peaks (also compare Figure A2 for the shear rate γ). At 800 rpm, the peaks of the energy dissipation rate close to the impeller in 2 cm to 3 cm reactor height are up to 280 times higher compared to the rest of the cultivation system. In this area, the mean energy dissipation rate ε across the diameter showed a maximum εmax,SP, with extent values of over 14 kW m−3, despite it staying between 50 and 100 W m−3 in the largest parts of the reactor. The system average volumetric energy dissipation rate εav, which corresponds to the averaged volumetric power input P/V, did not exceed 110 W m−3 for 800 rpm. In the TIS, additional peaks could again be found in the height of the six blade impellers at 5.5 cm and 7.0 cm reactor height. At 600 rpm, the slice plot maximum εmax,SP of the pitched blade impeller at 3 cm reactor height was over 6 kW m−3, and the maximum of the six blade impellers was only slightly over 3 kW m−3.
The slice plots (averaged values across the diameter) were obtained for different agitation rates to evaluate the spatial distribution of the investigated hydrodynamic parameters within both reactor systems (see Figure 1 and Figure 2c,d). They all showed peaks for all parameters in the area of the impellers, which were magnitudes higher than in the remaining parts of the system. From the plots, it can be concluded that the characteristic profiles of the six blade impellers were added to the profile of the pitched blade impeller. Pre-investigations to this study showed that these profiles were moved but not altered with a change in the impeller position. In the TIS, maxima from slice plots εmax,SP for the pitched blade impeller tended to have the same magnitudes as in the SIS for corresponding agitation rates. However, the characteristic peaks of the six blade impellers were added to the peak of the pitched blade impeller, which led to a higher energy dissipation rate εav averaged for the whole system, but lower peaks for the same system averaged energy dissipation rate εav. Hence, setups with multiple impellers yielded flatter energy dissipation rate ε distributions for the same averaged power inputs P/V in the same reactor vessels. Further values for agitation rates, averaged dissipation rates εav, and maximum energy dissipation rates εmax for both systems can be found in Table A1 and Table A2 in Appendix A.

3.2. Cultivation Results

Cultivation experiments with CHO DP-12 cells were run in both reactor systems at different agitation rates. For all samples of the cultivation experiments, cell densities and antibody concentrations were quantified. Mean exemplary data containing viable cell densities and viabilities for cultivations at 400, 800, and 1200 rpm in the SIS are shown in Figure 3.
Figure 4 and Figure 5 show the results of the cultivation experiments performed in the TIS. Surprisingly, the effect of the agitation rate on maximum cell density, viability, or growth rate seemed to be less pronounced as that for the SIS, even at 1400 rpm. Obviously, the stable operational range was broader when using multiple impellers.
From Figure 3 and Figure 5 (The maximum cell densities are summarized in Figure 5), it can be concluded that the maximal cell density and the cell specific growth rate (not shown) were more or less unaffected by the agitation rate up to 800 rpm. For 1000 rpm and higher agitation rates, a decrease in maximum cell densities and cell specific growth rate was observed. Furthermore, the viability had already started to decrease in processes at high agitation rates above 1000 rpm before reaching the maximum cell density.
With respect to the antibody production behavior of CHO DP-12 cells, the space time yield (STY) was calculated for all performed processes (Figure 6). For the SIS experiments, the highest STYs with up to 2.9 mg L−1 h−1 are observed between agitation rates of 400 rpm and 800 rpm; for higher agitation rates, the STY declined and at 1400 rpm, the SIS of the antibody concentration in all samples was even below the detection limit of 0.15 g L−1. For the TIS, the STYs were approx. constant at around 1.5 mg L−1 h−1 for agitation rates between 770 rpm and 1400 rpm. These results also support the idea of a broader stable operational range when using multiple impellers.

3.3. Comparison of Cellular Behavior and Hydrodynamic Parameters

In the following, the presented experimental findings from cell cultures and hydrodynamic parameters are compared. In particular, it must be noted that the cells in the multi-stage reactors can be cultivated stably over a much broader operation range. This comparison will be made on the basis of various definitions for the energy dissipation rate including averaged and locally resolved ones. On one hand, parameters frequently used for process design, average power input P/V, and Kolmogorov length scale λ based on average cell diameters dc,av are discussed. On the other hand, locally resolved parameters such as energy dissipation rate ε distributions from CFD simulations, critical energy dissipation rate εkrit distributions calculated from cell size dc distributions, maximum energy dissipation rates εmax, maximum energy dissipation rates from slice plots εmax,SP and variance coefficients of hydrodynamic parameters u, γ, and ε will be considered.
The Kolmogorov length scale λ [12,13] is a commonly used tool to evaluate the impact of hydrodynamics on mammalian cell cultures. With respect to flow induced cell damage, it is assumed that significant cell damage has to be expected when the size of the smallest turbulence eddy length is in the order of magnitude of the cell size. Therefore, for a known cell diameter dc, the respective critical energy dissipation rate εkrit and hence, an appropriate agitation rate, can be estimated. Usually, the system averaged energy dissipation rate εav is used to calculate the Kolmogorov length scale λ and compared to a mean cell diameter dc. The concept seems to work for microcarrier cultures, but is under discussion for suspension cells [14,15,16,17,18,19]
As an example, the average cell size dc of CHO DP-12 cells with about 10 µm would result in a critical energy dissipation rate εkrit of about 32 kW m−3, which would be far above the calculated average energy dissipation rates εav even at 1400 rpm, where a strong decrease in maximum viable cell density and antibody productivity were observed. Vice versa, an agitation rate calculated for process design with this approach would be way too high and might lead to undesired process behavior. Therefore, in order to evaluate a possible relationship between the cell size dc distribution and the Kolmogorov length scale λ, critical energy dissipation rate εkrit distributions (see Figure 7 and Figure 8) were calculated from actual cell size dc distributions. Exemplary cell size dc distributions obtained from cultivation data are shown in the Appendix (Figure A3). It was observed that cell size distributions neither expressively differed with changing agitation rates nor were dependent on the cultivation system.
The distribution of the normalized frequency of the energy dissipation rate ε from CFD data for different agitation rates is shown in Figure 7 and Figure 8. Furthermore, the averaged and maximum energy dissipation rates εav and εmax from CFD simulations as well as the maxima of the slice plots are indicated in the diagrams.
In the cultivation experiments (see Section 3.2) with the SIS, no decrease in cell growth or antibody productivity was observed up to 800 rpm. The averaged energy dissipation rate ε stayed a magnitude below the critical energy dissipation rate εkrit distribution curve, calculated from the cell size distributions. The maximum energy dissipation rates εmax from the CFD simulations were already above the main cell peak of the critical energy dissipation rate εkrit distribution at the lowest agitation rate. At 800 rpm, the maximum value was right below the peak that represents the main cell population. With increasing agitation rates, the maximum energy dissipation rate εmax also increases and, according to the Kolmogorov length scale hypothesis, for half of the cell population, shear related effects would have to be expected.
With CFD tools, it is possible to obtain normalized volumetric distributions of hydrodynamic parameters for the simulated systems [21,27]. Since the so far compared parameters might not provide satisfying information on the impact of shear forces, the actual normalized energy dissipation rate ε distribution in the bioreactor was compared to the critical energy dissipation rate εkrit distribution, calculated from the cell size distributions (see Figure 7 and Figure 8).
It was found that the energy dissipation rate in the largest part of the vessel volume stayed below the critical energy dissipation rate distribution curve. Only in a very small fraction did the energy dissipation rate in the reactor reach the area of the critical energy dissipation rate distribution. Although agitation rates above 800 rpm cannot be simulated, it can be concluded that the energy dissipation rate distribution curve of the reactor will move further into the critical energy dissipation rate distribution curve, according to extrapolations of the average and maximum values.
Additionally, for the TIS, the critical energy dissipation rate distributions were calculated and compared to the actual energy dissipation rate distribution in the cultivation system as well as the maximum energy dissipation rate from the slice plots. The results are shown in Figure 8.
For comparison of the different reactor systems, the plots for 600 rpm and 800 rpm in the SIS and 450 rpm in the TIS are shown in Figure 9. Their average volumetric power inputs were within the same order of magnitude (47 W m−3 to 107 W m−3) and clear differences are recognized.
In all cultivations in the TIS at 770 rpm to 1400 rpm, no noteworthy changes in cell growth or antibody productivity were present. From extrapolations of the slice plot maxima, which resulted in 64 kW m−3 at 1400 rpm, which was just at the end of the main cell population peak (modal value at 10 kW m−3 to 20 kW m−3), it can be derived that the maximum from slice plots might not be a suitable parameter to make a connection between hydrodynamics and biology. Additionally, for this setup, the energy dissipation rate stayed below the critical energy dissipation rate distribution curve in the vast majority of the reactor volume. To have an optical impression of critical volume fractions, the areas in the cultivation system with energy dissipation rates above 1000 W m−3 were identified. Detailed pictures can be found in the Appendix (Figure A4).
For further analysis of both reactor systems, the hydrodynamic homogeneity was evaluated. Therefore, variance coefficients were calculated as the standard deviation from slice plot data divided by their mean value, which are shown in Figure 10.
It can be seen that the variance coefficient of all three hydrodynamic parameters (fluid velocity u, shear rate γ, and energy dissipation rate ε) was higher in the SIS. This means that the SIS is a more heterogeneous system than the TIS. This might possibly be an explanation for the more stable process behavior of the TIS with respect to a changing agitation rate.

4. Discussion

The previously shown results from this study are discussed in the following under different aspects. The process stability that depends on hydrodynamics is debated and the application of the Kolmogorov length scale hypothesis with modern CFD methods is highlighted. Finally, consequences for the scale-up and scale-down of the bioreactors are drawn.

4.1. Hydrodynamic Differences for Single and Multi-Stage Impellers

Certain differences in flow pattern, velocity, and other hydrodynamic parameters have been observed for multiple impellers in the same stirred reactor system. It was found from CFD studies that the system averaged fluid velocity, shear rate, and energy dissipation increased with the addition of multiple impellers. The maximum shear rate and energy dissipation rate also increased in the reactor setup with multiple impellers. Only the maximum fluid velocity remained almost constant as it depends on the agitation rate solely. Overall, it can be concluded that the use of multiple impellers leads to a more uniform fluid flow and therefore to a lower variance coefficient of the hydrodynamic parameters.

4.2. Stability of Cell Culture Processes Related to Locally Resolved Hydrodynamics

Cultivation results for both reactor setups showed clear differences in the process behavior with respect to changing agitation rates. While the maximum cell density decreased up to 60% in the SIS, and antibody productivity decreased up to 100% above 1000 rpm (averaged 203 W m−3), no changes were recorded for agitation rates up to 1400 rpm (averaged 4742 W m−3) in the TIS. Platas et al. [28] reported a rather stable process behavior at higher agitation rates when using a bioreactor with multiple impellers instead of a single impeller. We identified that local hydrodynamic maxima and hydrodynamic heterogeneities, represented by variant coefficients were less pronounced in the TIS at the same averaged energy dissipation rate, resulting in a more homogeneous fluid flow. From the cultivation data, it could be concluded that this leads to a more robust process behavior.

4.3. Evaluation of Concepts for Estimation of Shear Related Parameters

From the results of the cultivation experiments, it can be concluded that the average volumetric power input or energy dissipation rate is clearly not suitable as a universal shear-related parameter for bioprocess design. The two examined reactor systems showed a different process behavior for the same averaged energy dissipation rates. Even if the average dissipation rate increases way faster in the TIS with increasing agitation rate, the process is still more stable with respect to the agitation rate in the TIS. Maximum energy dissipation rates are also not suitable.
None of the parameters calculated from CFD data corresponded to the cultivation process behavior with changing agitation rate for both cultivation systems. Hence, it can be assumed that these parameters are not sufficient as indicators for the impact of hydrodynamics on the cells. However, solely, the maximum energy dissipation rates from the slice plots εmax,SP could possibly show a connection to the cellular effects.
The Kolmogorov length scale is commonly used to describe the size of smallest eddies that might lead to cell damage [14,15,16,17,18,19]. In this study, it was calculated in the form of distributions instead of one fixed length for the whole bioreactor due to strong hydrodynamic heterogeneities in the cultivation systems. In addition, cell size distributions instead of a mean cell size were considered. The resulting critical energy dissipation rate distribution, calculated from cell size distributions, was compared to the averaged and maximum energy dissipation rate as well as to the maxima from slice plots and volumetric distributions of the energy dissipation rate, all from CFD simulations.
It was found that the Kolmogorov length scale is not suitable to describe a link between hydrodynamics and cell damage, if the averaged or maximum energy dissipation rate is used since they do not correlate with any biological observations (Figure 3, Figure 4, Figure 5 and Figure 6). Furthermore, critical volume fractions of the cultivation vessels do not deliver the desired insights, despite the fact that the energy dissipation rate distribution curve moves into the critical region with increasing agitation rate. Nevertheless, extrapolations of the maximum energy dissipation rates from slice plots showed a slightly faster increase for the TIS. Since cells grew with roughly the same growth rate up to the highest agitation rate in this setup, this parameter does not seem to be suitable to describe a link between cell growth and hydrodynamics, but still seems to be a closer estimation than the other parameters.
The hydrodynamic heterogeneity of a bioreactor, which was quantified in this study with variance coefficients of the hydrodynamic parameters fluid velocity, shear rate, and energy dissipation rate, might be a useful parameter to estimate the suitability of a cultivation system. The calculated variance coefficients of all hydrodynamic parameters were higher in the TIS than in the SIS, which might explain the rather stable process behavior in multiple impeller systems due to the improved hydrodynamic homogeneity.

4.4. Learning for Reactor Scale-Up and Scale-Down

When performing scale-up of bioreactors, common criteria that are kept constant are the impeller tip speed or the averaged power input [7,8,9,10]. In this study, the averaged power input increased much faster with increasing agitation rate in the TIS than in the SIS. However, no prominent differences in process behavior were found in cultivations at different agitation rates in the TIS, while in the SIS, the maximum cell densities and antibody productivities were reduced at high agitation rates. Consequently, the cultivation data showed that the biological process behavior does not solely depend on single process parameters such as the averaged power input, which are commonly used in process design, but rather on their homogeneity. Since maxima of hydrodynamic parameters also increased faster in the TIS and no relevant changes in process behavior were recorded, it is likely that the whole reactor system and locally resolved hydrodynamics need to be considered. Therefore, locally resolved hydrodynamic parameters were carved out of the CFD data as slice-plots and as volume fraction specific distributions. It was observed that strong hydrodynamic heterogeneities were present in stirred tank reactor systems, which became steeper with increasing agitation rate. A multi-parametrical approach considering locally resolved hydrodynamics resulting from geometrical characteristics could be the consequence for process scale-up. In turn, the investigation of local hydrodynamic effects has been the focus of late. Local distributions of the shear stress and the Kolmogorov length scale have been evaluated for spinner flasks, but have not been compared to cell size distributions thus far [27]. Furthermore, the calculation of cellular residence times in differently mixed reactor areas could provide deeper insights on the actual impact of hydrodynamic phenomena on cellular behavior [22,29]. Another recent study showed that it is worth considering the local distributions of conventional process design parameters such as kLa [10]. In addition, the local determination of the power input shows promising results for the design of cultivation processes for human mesenchymal stem cells. These approaches could be combined with novel uncertainty-based model evaluation methods [6,30]. It is noteworthy that the six blade impellers, specifically designed for the TIS in this study, worked immediately for the cultivation of CHO DP-12 cells.

5. Conclusions

In this study, parameters obtained from CFD simulations were linked to cell culture cultivation experiments to investigate the influence of hydrodynamic indifferent reactor setups on cell growth and antibody productivity. It was shown that hydrodynamically more uniform conditions and the resulting flatter hydrodynamic profiles might be the reason for the broader stable operational space regarding the agitation rate of a stirred tank with multiple impellers, which was previously observed by Platas et al. [28]. The evaluation of the CFD results clearly showed stronger pronounced hydrodynamic heterogeneities for the same power input in single impeller setups. Hence, the use of conventional, mostly averaged process design parameters, needs to be questioned and rather, local gradients should be considered, especially for scale-up [22,29]. Furthermore, it was found that the Kolmogorov length scale hypothesis might not be appropriate to describe the influence of hydrodynamics on mammalian cell culture processes without a hydrodynamic characterization of the full cultivation system including all hydrodynamic gradients. Overall, a deeper insight into local hydrodynamic gradients in cell culture reactors was obtained, but the reliable prediction of design parameters prior to an experimental evaluation still remains difficult.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pr10010107/s1, File S1: Supplementary File SIS Model Report.pdf, File S2: Supplementary File TIS Model Report.pdf.

Author Contributions

Conceptualization, F.F., J.M., R.P., and M.S.; Methodology, F.F.; Software, F.F. and J.B.; Validation, F.F. and J.M.; Formal analysis, F.F. and E.A.; Investigation, F.F., E.A. and J.B.; Resources, R.P.; Data curation, F.F., E.A. and J.B.; Writing—original draft preparation, F.F., J.M. and R.P.; Writing—review and editing, M.S.; Project administration, J.M. and R.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Pitched blade impeller (above) and six blade impeller (below) with corresponding CAD models. For the pitched blade impeller, a Newton number Ne of 0.35 was determined with an empiric correlation [7], while for the six blade impeller, 3.95 was considered, following another empiric correlation [31].
Figure A1. Pitched blade impeller (above) and six blade impeller (below) with corresponding CAD models. For the pitched blade impeller, a Newton number Ne of 0.35 was determined with an empiric correlation [7], while for the six blade impeller, 3.95 was considered, following another empiric correlation [31].
Processes 10 00107 g0a1
Figure A2. Shear rate γ in the SIS and TIS: (a) shear rate in the SIS; (b) shear rate in the TIS. The shown data correspond to an agitation rate of 400 rpm. (d) Mean shear rate along the reactor height (slice plot) in the SIS for different agitation rates; (d) mean shear rate along the reactor height (slice plot) in the TIS for different agitation rates. In (c,d), the shear rate was averaged across the diameter in steps of one millimeter along the reactor height.
Figure A2. Shear rate γ in the SIS and TIS: (a) shear rate in the SIS; (b) shear rate in the TIS. The shown data correspond to an agitation rate of 400 rpm. (d) Mean shear rate along the reactor height (slice plot) in the SIS for different agitation rates; (d) mean shear rate along the reactor height (slice plot) in the TIS for different agitation rates. In (c,d), the shear rate was averaged across the diameter in steps of one millimeter along the reactor height.
Processes 10 00107 g0a2
In both cases, for most parts of the reactor systems, a quite uniform distribution of the shear rate could be observed. The mean shear rate did not exceed 100 s−1, but reached almost 600 s−1 at the height of the stirrer in the SIS, and the shear rate showed values of up to 9167 s−1 at 400 rpm. Only small parts of the reactor reached values above 1000 s−1, namely at the stirrer. In the largest part along the reactor height, the mean shear rate did not exceed 100 s−1, but reached almost 600 s−1 at the height of the stirrer. The maximum values in the TIS were similar for equal agitation rates with 9383 s−1 at 400 rpm, but again, two additional peaks could be found in the slice plots, representing the added impellers. However, the values between the stirrers were still comparably low, below 100 s−1.
Figure A3. Exemplary cell size distributions from cultivation data in the SIS (800 rpm and 1200 rpm) and in the TIS (770 rpm and 1400 rpm). Distributions for all cultivations were taken from samples at 72 h cultivation time.
Figure A3. Exemplary cell size distributions from cultivation data in the SIS (800 rpm and 1200 rpm) and in the TIS (770 rpm and 1400 rpm). Distributions for all cultivations were taken from samples at 72 h cultivation time.
Processes 10 00107 g0a3
Figure A4. Areas in the TIS with energy dissipation rates above 1000 W m−3.
Figure A4. Areas in the TIS with energy dissipation rates above 1000 W m−3.
Processes 10 00107 g0a4
The regions above 1000 W m−3 are all located in thin layers on the pitched blade impellers, small areas around the six blade impellers, and very small parts at the edges of probes and other inserts.
With a defined threshold for the critical energy dissipation rate, the critical volume fractions of the cultivation system can be calculated from CFD data. They are linearly dependent on the power input or energy dissipation rate. Due to the respective number of impellers, the different reactor setups result in different critical reactor volume fractions for the same power input. The critical areas were smaller for equal averaged power inputs in the TIS due to flatter gradients. Cultivation experiments in the SIS showed a decrease in cell growth and antibody productivity above 1000 rpm. In the TIS, no changes were recorded. Since the critical reactor volume for 1400 rpm was about five-fold higher in the TIS (0.45% at 4.7 kW m−3) than in the SIS (0.07% at 465 W m−3), it can be concluded that the critical volume fraction might not be a suitable parameter to also evaluate the influence of hydrodynamics on the process.
Figure A5. Critical reactor volume fractions with energy dissipation rates above 20 kW m−3 plotted over the corresponding averaged power input. Black circles represent values for the SIS, grey squares those of the TIS.
Figure A5. Critical reactor volume fractions with energy dissipation rates above 20 kW m−3 plotted over the corresponding averaged power input. Black circles represent values for the SIS, grey squares those of the TIS.
Processes 10 00107 g0a5
Table A1. Hydrodynamic parameters in the SIS based on CFD simulations. All values above 800 rpm were extrapolated linear (velocity and shear rate) or quadratic, respectively (energy dissipation rate).
Table A1. Hydrodynamic parameters in the SIS based on CFD simulations. All values above 800 rpm were extrapolated linear (velocity and shear rate) or quadratic, respectively (energy dissipation rate).
N (rpm)uav (m s−1)γav (s−1)εav (W m−3)εmax (W m−3)εmax,SP (W m−3)
2000.03210.062.3922,169214
4000.06821.2315.35146,5101575
6000.10332.2247.35491,7674959
8000.13843.15107.121,159,62811,001
10000.17153.56203.562,265,06620,239
12000.20564.27345.053,932,42533,194
14000.24074.99540.506,286,05050,385
Table A2. Hydrodynamic parameters in the TIS based on CFD simulations. All values above 600 rpm were extrapolated linearly (velocity and shear rate) or quadratically, respectively (energy dissipation rate).
Table A2. Hydrodynamic parameters in the TIS based on CFD simulations. All values above 600 rpm were extrapolated linearly (velocity and shear rate) or quadratically, respectively (energy dissipation rate).
N (rpm)uav (m s−1)γav (s−1)εav (W m−3)εmax (W m−3)εmax,SP (W m−3)
1500.03913.243.6345,738221
3000.08128.6223.86194,0001118
4500.12543.9372.90479,0003014
6000.15955.53176.931,010,0006247
7700.20872.98454.352,007,73212,103
9300.25288.15964.253,401,19120,258
10800.292102.361744.795,172,10430,686
12000.325113.742633.296,952,68441,252
14000.379132.694742.4910,725,75263,867

References

  1. Hu, W.; Berdugo, C.; Chalmers, J.J. The Potential of Hydrodynamic Damage to Animal Cells of Industrial Relevance: Current Understanding. Cytotechnology 2011, 63, 445–460. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Godoy-Silva, R.; Chalmers, J.J.; Casnocha, S.A.; Bass, L.A.; Ma, N. Physiological responses of CHO cells to repetitive hydrodynamic stress. Biotechnol. Bioeng. 2009, 103, 1103–1117. [Google Scholar] [CrossRef] [PubMed]
  3. Chalmers, J.J.; Ma, N. Hydrodynamic Damage to Animal Cells. In Animal Cell Culture; Al-Rubeai, M., Ed.; Springer International Publishing: Basel, Switzerland, 2015; pp. 169–183. [Google Scholar]
  4. Nienow, A.W.; Scott, W.H.; Hewitt, C.J.; Thomas, C.R.; Lewis, G.; Amanullah, A.; Kiss, R.; Meier, S.J. Scale-down studies for assessing the impact of different stress parameters on growth and product quality during animal cell culture. Chem. Eng. Res. Des. 2013, 91, 2265–2274. [Google Scholar] [CrossRef] [Green Version]
  5. Jossen, V.; Eibl, D.; Eibl, R. Numerical Methods for the Design and Description of In Vitro Expansion Processes of Human Mesenchymal Stem Cells. Adv. Biochem. Eng. Biotechnol. 2020, 177, 185–228. [Google Scholar] [CrossRef]
  6. Arndt, L.; Wiegmann, V.; Kuchemüller, K.B.; Baganz, F.; Pörtner, R.; Möller, J. Model-based workflow for scale-up of process strategies developed in miniaturized bioreactor systems. Biotechnol. Prog. 2021, 37, e3122. [Google Scholar] [CrossRef]
  7. Doran, P.M. Bioprocess Engineering Principles, 2nd ed.; Academic Press: Cambridge, MA, USA, 2013. [Google Scholar]
  8. Storhas, W. Bioverfahrensentwicklung; John Wiley & Sons: Hoboken, NJ, USA, 2003; ISBN 3-527-28866-X. [Google Scholar]
  9. Xu, S.; Hoshan, L.; Jiang, R.; Gupta, B.; Brodean, E.; O’Neill, K.; Seamans, T.C.; Bowers, J.; Chen, H. A practical approach in bioreactor scale-up and process transfer using a combination of constant P/V and vvm as the criterion. Biotechnol. Prog. 2017, 33, 1146–1159. [Google Scholar] [CrossRef]
  10. Scully, J.; Considine, L.B.; Smith, M.T.; McAlea, E.; Jones, N.; O’Connell, E.; Madsen, E.; Power, M.; Mellors, P.; Crowley, J.; et al. Beyond heuristics: CFD-based novel multiparameter scale-up for geometrically disparate bioreactors demonstrated at industrial 2kL–10kL scales. Biotechnol. Bioeng. 2020, 117, 1710–1723. [Google Scholar] [CrossRef]
  11. Rosseburg, A.; Fitschen, J.; Wutz, J.; Wucherpfennig, T.; Schlüter, M. Hydrodynamic inhomogeneities in large scale stirred tanks—Influence on mixing time. Chem. Eng. Sci. 2018, 188, 208–220. [Google Scholar] [CrossRef]
  12. Kolmogorov, A. Dissipation of Energy in the Locally Isotropic Turbulence. Dokl. Akad. Nauk SSSR 1941, 32, 538–541. [Google Scholar]
  13. Kolmogorov, A. The Local Structure of Turbulence in Incompressible Viscous Fluid for Very Large Reynolds Numbers. Dokl. Akad. Nauk SSSR 1941, 30, 301–305. [Google Scholar]
  14. Cherry, R.S.; Papoutsakis, E.T. Hydrodynamic Effects on Cells in Agitated Tissue Culture Reactors. Bioprocess Eng. 1986, 1, 29–41. [Google Scholar] [CrossRef]
  15. Cherry, R.S.; Papoutsakis, E.T. Physical mechanisms of cell damage in microcarrier cell culture bioreactors. Biotechnol. Bioeng. 1988, 32, 1001–1014. [Google Scholar] [CrossRef]
  16. Cherry, R.S.; Papoutsakis, E.T. Growth and death rates of bovine embryonic kidney cells in turbulent microcarrier bioreactors. Bioprocess Eng. 1989, 4, 81–89. [Google Scholar] [CrossRef]
  17. Hua, J.; Erickson, L.E.; Yiin, T.Y.; Glasgow, L.A. A review of the effects of shear and interfacial phenomena on cell viability. Crit. Rev. Biotechnol. 1993, 13, 305–328. [Google Scholar] [CrossRef]
  18. Cherry, R.S. Animal cells in turbulent fluids: Details of the physical stimulus and the biological response. Biotechnol. Adv. 1993, 11, 279–299. [Google Scholar] [CrossRef]
  19. Croughan, M.S.; Hamel, J.F.; Wang, D.I. Hydrodynamic Effects on Animal Cells Grown in Microcarrier Cultures. Biotechnol. Bioeng. 1987, 29, 130–141. [Google Scholar] [CrossRef]
  20. Wang, G.; Haringa, C.; Tang, W.; Noorman, H.; Chu, J.; Zhuang, Y.; Zhang, S. Coupled metabolic-hydrodynamic modeling enabling rational scale-up of industrial bioprocesses. Biotechnol. Bioeng. 2020, 117, 844–867. [Google Scholar] [CrossRef]
  21. Ueki, M.; Tansho, N.; Sato, M.; Kanamori, H.; Ito, Y.; Kato, Y. Improved cultivation of CHO cells in bioreactor with reciprocal mixing. Authorea 2020. [Google Scholar] [CrossRef]
  22. Kuschel, M.; Fitschen, J.; Hoffmann, M.; von Kameke, A.; Schlüter, M.; Wucherpfennig, T. Validation of Novel Lattice Boltzmann Large Eddy Simulations (LB LES) for Equipment Characterization in Biopharma. Processes 2021, 9, 950. [Google Scholar] [CrossRef]
  23. Jossen, V.; Schirmer, C.; Mostafa Sindi, D.; Eibl, R.; Kraume, M.; Pörtner, R.; Eibl, D. Theoretical and practical issues that are relevant when scaling up hMSC microcarrier production porocesses. Stem Cells Int. 2016, 2016, 4760414. [Google Scholar] [CrossRef] [Green Version]
  24. Isu, G.; Morbiducci, U.; Nisco, G.d.; Kropp, C.; Marsano, A.; Deriu, M.A.; Zweigerdt, R.; Audenino, A.; Massai, D. Modeling methodology for defining a priori the hydrodynamics of a dynamic suspension bioreactor. Application to human induced pluripotent stem cell culture. J. Biomech. 2019, 94, 99–106. [Google Scholar] [CrossRef]
  25. Kaiser, S.; Jossen, V.; Schirmaier, C.; Eibl, D.; Brill, S.; van den Bos, C.; Eibl, R. Fluid Flow and Cell Proliferation of Mesenchymal Adipose-Derived Stem Cells in Small-Scale, Stirred, Single-Use Bioreactors. Chem. Ing. Tech. 2013, 85, 95–102. [Google Scholar] [CrossRef]
  26. Möller, J.; Kuchemüller, K.B.; Steinmetz, T.; Koopmann, K.S.; Pörtner, R. Model-assisted Design of Experiments as a concept for knowledge-based bioprocess development. Bioprocess Eng. 2019, 42, 867–882. [Google Scholar] [CrossRef]
  27. Ghasemian, M.; Layton, C.; Nampe, D.; Zur Nieden, N.I.; Tsutsui, H.; Princevac, M. Hydrodynamic characterization within a spinner flask and a rotary wall vessel for stem cell culture. Biochem. Eng. J. 2020, 157, 107533. [Google Scholar] [CrossRef]
  28. Platas Barradas, O.; Jandt, U.; Minh Phan, L.D.; Villanueva, M.E.; Schaletzky, M.; Rath, A.; Freund, S.; Reichl, U.; Skerhutt, E.; Scholz, S.; et al. Evaluation of Criteria for Bioreactor Comparison and Operation Standardization for Mammalian Cell Culture. Eng. Life Sci. 2012, 12, 518–528. [Google Scholar] [CrossRef]
  29. Fitschen, J.; Hofmann, S.; Wutz, J.; Hoffmann, M.; Wucherpfennig, T.; Schlüter, M. Novel evaluation method to determine the local mixing time distribution in stirred tank reactors. Chem. Eng. Sci. X 2021, 10, 100098. [Google Scholar] [CrossRef]
  30. Möller, J.; Hernández Rodríguez, T.; Müller, J.; Arndt, L.; Kuchemüller, K.B.; Frahm, B.; Eibl, R.; Eibl, D.; Pörtner, R. Model uncertainty-based evaluation of process strategies during scale-up of biopharmaceutical processes. Comput. Chem. Eng. 2020, 134, 106693. [Google Scholar] [CrossRef]
  31. Bates, R.L.; Fondy, P.L.; Corpstein, R.R. Examination of Some Geometric Parameters of Impeller Power. Ind. Eng. Chem. Proc. Des. Dev. 1963, 2, 310–314. [Google Scholar] [CrossRef]
Figure 1. Fluid velocity u in the SIS and TIS. (a) Fluid velocity in the SIS; (b) Fluid velocity in the TIS. The shown data correspond to an agitation rate of 400 rpm. (c) Mean fluid velocity along the reactor height (slice plot) in the SIS for different agitation rates; (d) mean fluid velocity along the reactor height (slice plot) in the TIS for different agitation rates. In (c,d), the flow velocity was averaged across the diameter in steps of one millimeter along the reactor height.
Figure 1. Fluid velocity u in the SIS and TIS. (a) Fluid velocity in the SIS; (b) Fluid velocity in the TIS. The shown data correspond to an agitation rate of 400 rpm. (c) Mean fluid velocity along the reactor height (slice plot) in the SIS for different agitation rates; (d) mean fluid velocity along the reactor height (slice plot) in the TIS for different agitation rates. In (c,d), the flow velocity was averaged across the diameter in steps of one millimeter along the reactor height.
Processes 10 00107 g001
Figure 2. Energy dissipation rate ε in the SIS and TIS: (a) energy dissipation rate in the SIS; (b) energy dissipation rate in the TIS. The shown data correspond to an agitation rate of 400 rpm. (d) Mean energy dissipation rate along the reactor height (slice plot) in the SIS for different agitation rates; (d) mean energy dissipation rate along the reactor height (slice plot) in the TIS for different agitation rates. In (c,d), the energy dissipation rate was averaged across the diameter in steps of one millimeter along the reactor height.
Figure 2. Energy dissipation rate ε in the SIS and TIS: (a) energy dissipation rate in the SIS; (b) energy dissipation rate in the TIS. The shown data correspond to an agitation rate of 400 rpm. (d) Mean energy dissipation rate along the reactor height (slice plot) in the SIS for different agitation rates; (d) mean energy dissipation rate along the reactor height (slice plot) in the TIS for different agitation rates. In (c,d), the energy dissipation rate was averaged across the diameter in steps of one millimeter along the reactor height.
Processes 10 00107 g002
Figure 3. Mean growth curves (viable cell density X and viability) for the cultivations in the SIS at different agitation rates and respective averaged power inputs: 400 rpm (16 W m−3); 800 rpm (110 W m−3); 1200 rpm (356 W m−3). Data were averaged from three cultivations each.
Figure 3. Mean growth curves (viable cell density X and viability) for the cultivations in the SIS at different agitation rates and respective averaged power inputs: 400 rpm (16 W m−3); 800 rpm (110 W m−3); 1200 rpm (356 W m−3). Data were averaged from three cultivations each.
Processes 10 00107 g003
Figure 4. Growth curves (viable cell density X and viability) for the cultivations in the TIS at different agitation rates and respective averaged power inputs. The agitation rate settings range from 770 rpm or 454 W m−3 to 1400 rpm or 4742 W m−3. One experiment was performed at each agitation rate.
Figure 4. Growth curves (viable cell density X and viability) for the cultivations in the TIS at different agitation rates and respective averaged power inputs. The agitation rate settings range from 770 rpm or 454 W m−3 to 1400 rpm or 4742 W m−3. One experiment was performed at each agitation rate.
Processes 10 00107 g004
Figure 5. Maximum viable cell densities (VCD) compared for experiments in both reactor setups: cultivation results from the SIS (grey bars); cultivation results from the TIS (blue bars). Data for 400, 800, and 1200 rpm in the SIS were averaged from three cultivations each. Dashed arrows indicate qualitative trends for each reactor setup.
Figure 5. Maximum viable cell densities (VCD) compared for experiments in both reactor setups: cultivation results from the SIS (grey bars); cultivation results from the TIS (blue bars). Data for 400, 800, and 1200 rpm in the SIS were averaged from three cultivations each. Dashed arrows indicate qualitative trends for each reactor setup.
Processes 10 00107 g005
Figure 6. Antibody space time yields (STY) compared for experiments in both reactor setups: cultivation results from the SIS (grey bars); cultivation results from the TIS (blue bars). Data for 400, 800, and 1200 rpm in the SIS were averaged from three cultivations each. Dashed arrows indicate expected trends for each reactor setup.
Figure 6. Antibody space time yields (STY) compared for experiments in both reactor setups: cultivation results from the SIS (grey bars); cultivation results from the TIS (blue bars). Data for 400, 800, and 1200 rpm in the SIS were averaged from three cultivations each. Dashed arrows indicate expected trends for each reactor setup.
Processes 10 00107 g006
Figure 7. Critical energy dissipation rate εkrit distributions (black circles) calculated from the cell size distributions in the SIS compared to the energy dissipation rate distribution in the reactor system (grey squares) obtained from CFD data. The averaged energy dissipation rates from CFD simulations are indicated in blue, their maxima in red, maximum energy dissipation rates calculated from slice plots are indicated in purple.
Figure 7. Critical energy dissipation rate εkrit distributions (black circles) calculated from the cell size distributions in the SIS compared to the energy dissipation rate distribution in the reactor system (grey squares) obtained from CFD data. The averaged energy dissipation rates from CFD simulations are indicated in blue, their maxima in red, maximum energy dissipation rates calculated from slice plots are indicated in purple.
Processes 10 00107 g007
Figure 8. Critical energy dissipation rate distributions (black circles) in the TIS compared to the energy dissipation rate distribution in the reactor system (grey squares). The averaged energy dissipation rates from CFD simulations are indicated in blue, their maxima in red, maximum energy dissipation rates calculated from slice plots are indicated in purple.
Figure 8. Critical energy dissipation rate distributions (black circles) in the TIS compared to the energy dissipation rate distribution in the reactor system (grey squares). The averaged energy dissipation rates from CFD simulations are indicated in blue, their maxima in red, maximum energy dissipation rates calculated from slice plots are indicated in purple.
Processes 10 00107 g008
Figure 9. Critical energy dissipation rate distributions (black circles) compared to the energy dissipation rate distribution in the SIS (600 rpm and 800 rpm) and in the TIS (450 rpm) (grey squares). The averaged energy dissipation rates from CFD simulations are indicated in blue, their maxima in red. Maximum energy dissipation rates calculated from slice plots are indicated in purple.
Figure 9. Critical energy dissipation rate distributions (black circles) compared to the energy dissipation rate distribution in the SIS (600 rpm and 800 rpm) and in the TIS (450 rpm) (grey squares). The averaged energy dissipation rates from CFD simulations are indicated in blue, their maxima in red. Maximum energy dissipation rates calculated from slice plots are indicated in purple.
Processes 10 00107 g009
Figure 10. Variance coefficients for hydrodynamic parameters: (a) fluid velocity u, (b) shear rate γ, and (c) energy dissipation rate ε in both cultivation systems. The coefficients were calculated from standard deviations from slice plot data divided by their mean value. They represent a measure of hydrodynamic homogeneity of the respective reactor system.
Figure 10. Variance coefficients for hydrodynamic parameters: (a) fluid velocity u, (b) shear rate γ, and (c) energy dissipation rate ε in both cultivation systems. The coefficients were calculated from standard deviations from slice plot data divided by their mean value. They represent a measure of hydrodynamic homogeneity of the respective reactor system.
Processes 10 00107 g010
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Freiberger, F.; Budde, J.; Ateş, E.; Schlüter, M.; Pörtner, R.; Möller, J. New Insights from Locally Resolved Hydrodynamics in Stirred Cell Culture Reactors. Processes 2022, 10, 107. https://doi.org/10.3390/pr10010107

AMA Style

Freiberger F, Budde J, Ateş E, Schlüter M, Pörtner R, Möller J. New Insights from Locally Resolved Hydrodynamics in Stirred Cell Culture Reactors. Processes. 2022; 10(1):107. https://doi.org/10.3390/pr10010107

Chicago/Turabian Style

Freiberger, Fabian, Jens Budde, Eda Ateş, Michael Schlüter, Ralf Pörtner, and Johannes Möller. 2022. "New Insights from Locally Resolved Hydrodynamics in Stirred Cell Culture Reactors" Processes 10, no. 1: 107. https://doi.org/10.3390/pr10010107

APA Style

Freiberger, F., Budde, J., Ateş, E., Schlüter, M., Pörtner, R., & Möller, J. (2022). New Insights from Locally Resolved Hydrodynamics in Stirred Cell Culture Reactors. Processes, 10(1), 107. https://doi.org/10.3390/pr10010107

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop