Design of a large-scale surface-aerated bioreactor for biomass production using a VOC substrate

The design of a large-scale bioreactor for the production of bacterial biomass adapted to the biodegradation of volatile organic compounds was carried out. The bioreactor model used integrated the microbial kinetics and fluid dynamics described by the compartment model approach. The process conditions and kinetic parameters were adopted from the laboratory experimental study of (León, E., Seignez, C., Adler, N., Péringer, P., 1999. Growth inhibition of biomass adapted to the degradation of toluene and xylenes in mixture in a batch reactor with substrates supplied by pulses. Biodegradation 10, 245-250). The performance of the pulsed-batch stirred bioreactor under surface aeration conditions was simulated for different mixing configurations and conditions such as the impeller diameter, number of impellers, stirring speed, and oxygen pressure. The simulations were used for the cost analysis which resulted in the optimal design of the bioreactor.     

CFD predicted pH gradients in lactic acid bacteria cultivations

The formation of pH gradients in a 700 L batch fermentation of Streptococcus thermophilus was studied using multi-position pH measurements and computational fluid dynamics (CFD) modeling. To this end, a dynamic, kinetic model of S. thermophilus and a pH correlation were integrated into a validated one-phase CFD model, and a dynamic CFD simulation was performed. First, the fluid dynamics of the CFD model were validated with NaOH tracer pulse mixing experiments. Mixing experiments and simulations were performed whereas multiple pH sensors, which were placed vertically at different locations in the bioreactor, captured the response. A mixing time of about 46 s to reach 95% homogeneity was measured and predicted at an impeller speed of 242 rpm. The CFD simulation of the S. thermophilus fermentation captured the experimentally observed pH gradients between a pH of 5.9 and 6.3, which occurred during the exponential growth phase. A pH higher than 7 was predicted in the vicinity of the base solution inlet. Biomass growth, lactic acid production, and substrate consumption matched the experimental observations. Moreover, the biokinetic results obtained from the CFD simulation were similar to a single-compartment simulation, for which a homogeneous distribution of the pH was assumed. This indicates no influence of pH gradients on growth in the studied bioreactor. This study verified that the pH gradients during a fermentation in the pilot-scale bioreactor could be accurately predicted using a coupled simulation of a biokinetic and a CFD model. To support the understanding and optimization of industrial-scale processes, future biokinetic CFD studies need to assess multiple types of environmental gradients, like pH, substrate, and dissolved oxygen, especially at industrial scale.     

Scale‐up analysis for a CHO cell culture process in large‐scale bioreactors

Bioprocess scale‐up is a fundamental component of process development in the biotechnology industry. When scaling up a mammalian cell culture process, it is important to consider factors such as mixing time, oxygen transfer, and carbon dioxide removal. In this study, cell‐free mixing studies were performed in production scale 5,000‐L bioreactors to evaluate scale‐up issues. Using the current bioreactor configuration, the 5,000‐L bioreactor had a lower oxygen transfer coefficient, longer mixing time, and lower carbon dioxide removal rate than that was observed in bench scale 5‐ and 20‐L bioreactors. The oxygen transfer threshold analysis indicates that the current 5,000‐L configuration can only support a maximum viable cell density of 7 × 10⁶ cells mL^?1. Moreover, experiments using a dual probe technique demonstrated that pH and dissolved oxygen gradients may exist in 5,000‐L bioreactors using the current configuration. Empirical equations were developed to predict mixing time, oxygen transfer coefficient, and carbon dioxide removal rate under different mixing‐related engineering parameters in the 5,000‐L bioreactors. These equations indicate that increasing bottom air sparging rate is more efficient than increasing power input in improving oxygen transfer and carbon dioxide removal. Furthermore, as the liquid volume increases in a production bioreactor operated in fed‐batch mode, bulk mixing becomes a challenge. The mixing studies suggest that the engineering parameters related to bulk mixing and carbon dioxide removal in the 5,000‐L bioreactors may need optimizing to mitigate the risk of different performance upon process scale‐up. Biotechnol. Bioeng. 2009;103: 733–746. © 2009 Wiley Periodicals, Inc.     

Estimation of axial dispersion in horizontal rotating tubular bioreactor by means of a structured model

A horizontal rotating tubular bioreactor (HRTB) is designed as the combination of a "thin layer bioreactor" and a "biodisc" reactor. The investigation of mixing in HRTB was done by the temperature step method in a wide range of process conditions [residence time (t_z=360036000 s) and bioreactor rotation speed (n=0.0830.917 s^у)]. In all experiments heat losses were detected. A mathematical model based on "tank in series" concept was developed to describe the mixing in HRTB - a "spiral flow" model (SFM) which has incorporated heat losses. However, the simulations of SFM could be used for calculation of temperature response curves for the case when there is no heat losses. These corrected curves were used then to estimate Bodenstein number as a parameter of standard dispersion model (SDM). The obtained Bodenstein numbers were in the range 10-17. The simulations showed that SFM was more capable to describe the mixing in HRTB giving better fitting with experimental measurements than SDM, indicating that mixing pattern in HRTB is too complex to be described with this relatively simple, one-parameter model.     

Dynamic modeling and optimal fed-batch feeding strategies for a two-phase partitioning bioreactor

A dynamic model for the degradation of phenol in a two-phase partitioning bioreactor has been developed based on mechanistic balances around the bioreactor. The key process characteristics including substrate transfer between the organic and aqueous phases, substrate inhibition, oxygen limitation, and cell entrainment were incorporated into the model. The model predictions were validated against existing experimental data obtained for a 2-L bioreactor, and good correlation was observed for the time frames of the simulations, as well as for trends in cell and substrate concentrations. Optimal fed-batch, phenol feeding strategies were then developed based on two approaches: (1) maximization of phenol consumption in a fixed time interval and (2) consumption of a fixed amount of phenol in minimal time. The optimal feeding policies, determined using the Iterative Dynamic Programming algorithm, provided substantial improvements in the amount of phenol consumed when compared to a typical experimental heuristic approach. For example, 45.73 g of phenol was predicted to be consumed in 50 h (not including lag phase) using the optimal feeding profile compared to 10.26 g of phenol consumed in the simulated experimental approach. Oxygen limitation was predicted to be a recurring operational challenge in the partitioning bioreactor, and had a strong impact on the optimization results.     

Studies on mixing in horizontal rotating tubular bioreactor

A horizontal rotating tubular bioreactor (HRTB) is a combination of a “thin-layer bioreactor” and a “biodisc” reactor. Its interior is divided by O-ring shaped partition walls. Mixing properties of this new type of the bioreactor were investigated by using a temperature step method. The mixing simulations were done by Runge-Kutta-Fehlberg numerical integration. Adjustable parameters of the “spiral flow” model were optimised by Monte-Carlo method. In this investigation, the structured “spiral flow” model (containing four adjustable parameters) was tested in a wide range of experimental conditions. The results show that the structured “spiral flow” model is capable to describe the mixing in HRTB in the whole range of both bioreactor operational parameters (n and D).     

CFD of mixing of multi‐phase flow in a bioreactor using population balance model

Mixing in bioreactors is known to be crucial for achieving efficient mass and heat transfer, both of which thereby impact not only growth of cells but also product quality. In a typical bioreactor, the rate of transport of oxygen from air is the limiting factor. While higher impeller speeds can enhance mixing, they can also cause severe cell damage. Hence, it is crucial to understand the hydrodynamics in a bioreactor to achieve optimal performance. This article presents a novel approach involving use of computational fluid dynamics (CFD) to model the hydrodynamics of an aerated stirred bioreactor for production of a monoclonal antibody therapeutic via mammalian cell culture. This is achieved by estimating the volume averaged mass transfer coefficient (k_La) under varying conditions of the process parameters. The process parameters that have been examined include the impeller rotational speed and the flow rate of the incoming gas through the sparger inlet. To undermine the two‐phase flow and turbulence, an Eulerian‐Eulerian multiphase model and k‐ε turbulence model have been used, respectively. These have further been coupled with population balance model to incorporate the various interphase interactions that lead to coalescence and breakage of bubbles. We have successfully demonstrated the utility of CFD as a tool to predict size distribution of bubbles as a function of process parameters and an efficient approach for obtaining optimized mixing conditions in the reactor. The proposed approach is significantly time and resource efficient when compared to the hit and trial, all experimental approach that is presently used. © 2016 American Institute of Chemical Engineers Biotechnol. Prog., 32:613–628, 2016     

Studies on mixing in horizontal rotating tubular bioreactor

In previous investigations on mixing in a horizontal rotating tubular bioreactor (HRTB) the structured “spiral flow” model was developed which contained four adjustable parameters [1, 2]. In order to incorporate the mixing model in a semifundamental scale-up procedure it was necessary to make a relation between the adjustable model parameters and process parameters of the bioreactor expressed as dimensionless numbers. Mathematical equations which relate adjustable model parameters with dimensionless numbers were developed by non-linear and surface regression methods. These equations were applied to develop the prediction systems for adjustable model parameters. In total, nine systems of equations for the prediction of the adjustable model parameters were established and examined by simulation. Three of them (SC-2, SC-6 and SC-9) were selected as adequate to describe the mixing performance of HRTB in a wide range of process conditions.     

Dual-chambered membrane bioreactor for coculture of stratified cell populations

We have developed a dual-chambered bioreactor (DCB) that incorporates a membrane to study stratified 3D cell populations for skin tissue engineering. The DCB provides adjacent flow lines within a common chamber; the inclusion of the membrane regulates flow layering or mixing, which can be exploited to produce layers or gradients of cell populations in the scaffolds. Computational modeling and experimental assays were used to study the transport phenomena within the bioreactor. Molecular transport across the membrane was defined by a balance of convection and diffusion; the symmetry of the system was proven by its bulk convection stability, while the movement of molecules from one flow line to the other is governed by coupled convection-diffusion. This balance allowed the perfusion of two different fluids, with the membrane defining the mixing degree between the two. The bioreactor sustained two adjacent cell populations for 28 days, and was used to induce indirect adipogenic differentiation of mesenchymal stem cells due to molecular cross-talk between the populations. We successfully developed a platform that can study the dermis–hypodermis complex to address limitations in skin tissue engineering. Furthermore, the DCB can be used for other multilayered tissues or the study of communication pathways between cell populations.     

Mixing characteristics and liquid circulation in a new multi-environment bioreactor

The theoretical and experimental aspects of the hydrodynamics and mixing in a new multi-environment bioreactor that uses the air-lift design were investigated. This study focused on the mixing characteristics, residence time distribution, liquid circulation between three zones of aerobic, microaerophilic and anoxic, and liquid displacement in the bioreactor at influent flow rates of 720–1,450 L/day and air flow rates of 15–45 L/min. The theoretical analysis of liquid displacement led to the estimation of the specific rate of liquid discharge from the bioreactor at any given influent flow rate, and the number of liquid circulations between various bioreactor zones before the discharge of a given quantity of wastewater. The ratio of mean residence time to the overall hydraulic retention time (t _m/HRT) decreased with the increase of air flow rate at any given influent flow rate, and approached unity at higher air flow rates. Mixing was characterized in terms of the axial dispersion coefficient and Bodenstein number, demonstrating a linear relationship with the superficial gas velocity. A correlation was developed between the Bodenstein number and the Froude number. The study of liquid circulation between the zones showed that less than 1.5 % of the circulating liquid escapes circulation at each cycle and flows towards the outer clarifier, while the percentage of escaped liquid decreases with increasing air flow rate at a given influent flow rate. The specific rate of liquid discharge from the bioreactor increased from 0.19 to 0.69 h^?1 with the increase of air and influent flow rates from 15 to 45 L/min and 500 to 1,450 L/day, respectively. Under the examined operating conditions, mixed liquor circulates between 364 and 1,698 times between the aerobic, microaerophilic and anoxic zones before 99 % of its original volume is replaced by the influent wastewater.     

Numerical simulation of global hydro-dynamics in a pulsatile bioreactor for cardiovascular tissue engineering

Previous numerical simulations of the hydro-dynamic response in the various bioreactor designs were mostly concentrated on the local flow field analysis using computational fluid dynamics, which cannot provide the global hydro-dynamics information to assist the bioreactor design. In this research, a mathematical model is developed to simulate the global hydro-dynamic changes in a pulsatile bioreactor design by considering the flow resistance, the elasticity of the vessel and the inertial effect of the media fluid in different parts of the system. The developed model is used to study the system dynamic response in a typical pulsatile bioreactor design for the culturing of cardiovascular tissues. Simulation results reveal the detailed pressure and flow-rate changes in the different positions of the bioreactor, which are very useful for the evaluation of hydro-dynamic performance in the bioreactor designed. Typical pressure and flow-rate changes simulated agree well with the published experimental data, thus validates the mathematical model developed. The proposed mathematical model can be used for design optimization of other pulsatile bioreactors that work under different experimental conditions and have different system configurations.     

Flow dynamics within a bioreactor for tissue engineering by residence time distribution analysis combined with fluorescence and magnetic resonance imaging to investigate forced permeability and apparent diffusion coefficient in a perfusion cell culture chamber

This study reveals that residence time distribution (RTD) analysis with pH monitoring after acid bolus injection can be used to globally study the flow dynamics of a perfusion bioreactor, while fluorescence microscopy and magnetic resonance imaging (MRI) were used to locally investigate mass transport within a hydrogel scaffold seeded or not with cells. The bioreactor used in this study is a close‐loop tubular reactor. A dispersion model in one dimension has been used to describe the non‐ideal behavior of the reactor. From open‐loop experiments (single‐cycle analysis), the presence of stagnant zones and back mixing were observed. The impact of the flow rate, the compliance chamber volume and mixing were investigated. Intermediate flows (30, 45, 60, and 90 mL min⁻¹) had no effect over RTD function expressed in reduced time (θ). Lower flow rates (5 and 15 mL min⁻¹) were associated to smaller extent of dispersion. The compliance chamber volume greatly affected the dynamics of the RTD function, while the effects of mixing and flow were small to non‐significant. An empirical equation has been proposed to localize minima of the RTD function and to predict Pe_r. Finally, cells seeded in a gelatin gel at a density of 800,000 cells mL⁻¹ had no effect over the permeability and the apparent diffusion coefficient, as revealed by fluorescent microscopy and MRI experiments. Biotechnol. Bioeng. 2011;108: 2488–2498. © 2011 Wiley Periodicals, Inc.     

Substrate gradient formation in the large-scale bioreactor lowers cell yield and increases by-product formation

A heterogeneous micro-environment was identified in a 12 m³ bioreactor with a height-to-diameter ratio of 2.5. The reactor was aerated by a ring sparger and stirred by three Rushton turbines. E. coli cells were cultivated in minimal medium to a cell density in the order of 30?g/l. Samples of glucose, the growth limiting component fed to the process, were taken at three levels in the bioreactor (top/middle/bottom). These showed that glucose concentration declined away from the feedpoint. The gradients depended on the mixing characteristics of the feedpoint, and concentrations of up to 400 times the mean value were found when feed was added to a relatively stagnant mixing zone. This resulted in up to 20% lower biomass yield compared to the bench scale. Gradients also affected the by-product formation, resulting in acetate formation in the large-scale bioreactor. IPTG induction of a recombinant protein was shown to influence important cell parameters and considerably increased the yield of carbon dioxide per glucose added, indicating an increased maintenance. The product formation rate was, however, not notably affected by the scale-up.     

A noninvasive technique for the measurement of the energetic state of free‐suspension mammalian cells

A perfusion small‐scale bioreactor allowing on‐line monitoring of the cell energetic state was developed for free‐suspension mammalian cells. The bioreactor was designed to perform in vivo nuclear magnetic resonance (NMR) spectroscopy, which is a noninvasive and nondestructive method that permits the monitoring of intracellular nutrient concentrations, metabolic precursors and intermediates, as well as metabolites and energy shuttles, such as ATP, ADP, and NADPH. The bioreactor was made of a 10‐mm NMR tube following a fluidized bed design. Perfusion flow rate allowing for adequate oxygen supply was found to be above 0.79 mL min^?1 for high‐density cell suspensions (10⁸ cells). Chinese hamster ovary (CHO) cells were studied here as model system. Hydrodynamic studies using coloration/decoloration and residence time distribution measurements were realized to perfect bioreactor design as well as to determine operating conditions bestowing adequate homogeneous mixing and cell retention in the NMR reading zone. In vivo ³¹P NMR was performed and demonstrated the small‐scale bioreactor platform ability to monitor the cell physiological behavior for 30‐min experiments. © 2009 American Institute of Chemical Engineers Biotechnol. Prog., 2010     

Scale down of recombinant protein production: a comparative study of scaling performance

A large bioreactor is heterogeneous with respect to concentration gradients of substrates fed to the reactor such as oxygen and growth limiting carbon source. Gradient formation will highly depend on the fluid dynamics and mass transfer capacity of the reactor, especially in the area in which the substrate is added. In this study, some production-scale (12 m³ bioreactor) conditions of a recombinant Escherichia coli process were imitated on a laboratory scale. From the large-scale cultivations, it was shown that locally high concentration of the limiting substrate fed to the process, in this case glucose, existed at the level of the feedpoint. The large-scale process was scaled down from: (i) mixing time experiments performed in the large-scale bioreactor in order to identify and describe the oscillating environment and (ii) identification of two distinct glucose concentration zones in the reactor. An important parameter obtained from mixing time experiments was the residence time in the feed zone of about 10 seconds. The size of the feed zone was estimated to 10%. Based on these observations the scale-down reactor with two compartments was designed. It was composed of one stirred tank reactor and an aerated plug flow reactor, in which the effect of oscillating glucose concentration on biomass yield and acetate formation was studied. Results from these experiments indicated that the lower biomass yield and higher acetate formation obtained on a large scale compared to homogeneous small-scale cultivations were not directly caused by the cell response to the glucose oscillation. This was concluded since no acetate was accumulated during scale-down experiments. An explanation for the differences in results between the two reactor scales may be a secondary effect of high glucose concentration resulting in an increased glucose metabolism causing an oxygen consumption rate locally exceeding the transfer rate. The results from pulse response experiments and glucose concentration measurements, at different locations in the reactor, showed a great consistency for the two feeding/pulse positions used in the large-scale bioreactor. Furthermore, measured periodicity from mixing data agrees well with expected circulation times for each impeller volume. Conclusions are drawn concerning the design of the scale-down reactor.     

Influence of oxygen concentration on the metabolism of Penicillium chrysogenum

In large-scale bioreactors, there is often insufficient mixing and as a consequence, cells experience uneven substrate and oxygen levels that influence product formation. In this study, the influence of dissolved oxygen (DO) gradients on the primary and secondary metabolism of a high producing industrial strain of Penicillium chrysogenum was investigated. Within a wide range of DO concentrations, obtained under chemostat conditions, we observed different responses from P. chrysogenum: (i) no influence on growth or penicillin production (>0.025 mmol L⁻¹); (ii) reduced penicillin production, but no growth limitation (0.013–0.025 mmol L⁻¹); and (iii) growth and penicillin production limitations (<0.013 mmol L⁻¹). In addition, scale down experiments were performed by oscillating the DO concentration in the bioreactor. We found that during DO oscillation, the penicillin production rate decreased below the value observed when a constant DO equal to the average oscillating DO value was used. To understand and predict the influence of oxygen levels on primary metabolism and penicillin production, we developed a black box model that was linked to a detailed kinetic model of the penicillin pathway. The model simulations represented the experimental data during the step experiments; however, during the oscillation experiments the predictions deviated, indicating the involvement of the central metabolism in penicillin production.     

Incomplete mixing in large bioreactors — a study of its role in the fermentative production of streptokinase

The production of streptokinase in a batch fermentation has been analysed for the role of incomplete macromixing of the broth. The analysis is based on a kinetic model exhibiting inhibition by the substrate and a primary metabolite (lactic acid), and a mixing model comprising two continuous flow reactors (CFRs) with closed-loop recycle. The inoculum is introduced into one region (one CFR) and the mixing process determines its distribution, growth and reactivity. By varying the dilution rates of the CFRs, any degree of macromixing can be simulated. For dilution rates larger than 1.0 h^?1 almost complete macromixing is achieved, for which an analogy has been drawn with micromixing. Increasing the volume of the inoculated region relative to the noninoculated region improves the maximum attainable activity of streptokinase and shortens the time for this. In such a situation an imperfectly mixed bioreactor is superior to a perfectly mixed one, implying that good productivity requires a large inoculated region and incomplete macromixing. These inferences are supported by earlier studies of fluid mixing and relaxation times in bioreactors.     

Fast dynamic response of the fermentative metabolism of Escherichia coli to aerobic and anaerobic glucose pulses

The response of Escherichia coli cells to transient exposure (step increase) in substrate concentration and anaerobiosis leading to mixed‐acid fermentation metabolism was studied in a two‐compartment bioreactor system consisting of a stirred tank reactor (STR) connected to a mini‐plug‐flow reactor (PFR: BioScope, 3.5 mL volume). Such a system can mimic the situation often encountered in large‐scale, fed‐batch bioreactors. The STR represented the zones of a large‐scale bioreactor that are far from the point of substrate addition and that can be considered as glucose limited, whereas the PFR simulated the region close to the point of substrate addition, where glucose concentration is much higher than in the rest of the bioreactor. In addition, oxygen‐poor and glucose‐rich regions can occur in large‐scale bioreactors. The response of E. coli to these large‐scale conditions was simulated by continuously pumping E. coli cells from a well stirred, glucose limited, aerated chemostat (D = 0.1 h^?1) into the mini‐PFR. A glucose pulse was added at the entrance of the PFR. In the PFR, a total of 11 samples were taken in a time frame of 92 s. In one case aerobicity in the PFR was maintained in order to evaluate the effects of glucose overflow independently of oxygen limitation. Accumulation of acetate and formate was detected after E. coli cells had been exposed for only 2 s to the glucose‐rich (aerobic) region in the PFR. In the other case, the glucose pulse was also combined with anaerobiosis in the PFR. Glucose overflow combined with anaerobiosis caused the accumulation of formate, acetate, lactate, ethanol, and succinate, which were also detected as soon as 2 s after of exposure of E. coli cells to the glucose and O₂ gradients. This approach (STR‐mini‐PFR) is useful for a better understanding of the fast dynamic phenomena occurring in large‐scale bioreactors and for the design of modified strains with an improved behavior under large‐scale conditions. Biotechnol. Bioeng. 2009; 104: 1153–1161. © 2009 Wiley Periodicals, Inc.     

A novel planar flow cell for studies of biofilm heterogeneity and flow-biofilm interactions

Biofilms are microbial communities growing on surfaces, and are ubiquitous in nature, in bioreactors, and in human infection. Coupling between physical, chemical, and biological processes is known to regulate the development of biofilms; however, current experimental systems do not provide sufficient control of environmental conditions to enable detailed investigations of these complex interactions. We developed a novel planar flow cell that supports biofilm growth under complex two-dimensional fluid flow conditions. This device provides precise control of flow conditions and can be used to create well-defined physical and chemical gradients that significantly affect biofilm heterogeneity. Moreover, the top and bottom of the flow chamber are transparent, so biofilm growth and flow conditions are fully observable using non-invasive confocal microscopy and high-resolution video imaging. To demonstrate the capability of the device, we observed the growth of Pseudomonas aeruginosa biofilms under imposed flow gradients. We found a positive relationship between patterns of fluid velocity and biofilm biomass due to faster microbial growth under conditions of greater local nutrient influx, but this relationship eventually reversed because high hydrodynamic shear leads to the detachment of cells from the surface. These results reveal that flow gradients play a critical role in the development of biofilm communities. By providing new capability for observing biofilm growth, solute and particle transport, and net chemical transformations under user-specified environmental gradients, this new planar flow cell system has broad utility for studies of environmental biotechnology and basic biofilm microbiology, as well as applications in bioreactor design, environmental engineering, biogeochemistry, geomicrobiology, and biomedical research.     

New miniaturized hollow-fiber bioreactor for in vivo like cell culture, cell expansion, and production of cell-derived products

We have developed a miniaturized hollow-fiber bioreactor system for mammalian cell culture with a volume of 1 mL. Cell and medium compartments of the bioreactor are separated by a semipermeable membrane, and oxygenation of the cell compartment is accomplished using an oxygenation membrane. As a result of the geometry of the transparent housing, cells can be observed by microscopy during culture. The leukemic cell lines CCRF-CEM, HL-60, and REH were cultivated up to densities of 3.5 x 10(7)/mL without medium change or manipulation of the cells. As shown using CCRF-CEM cells, growth in the bioreactor was strongly influenced and could be controlled by the medium flow rate. As a consequence, consumption of glucose and generation of lactate varied with flow rate. Depending on the molecular size cutoff of the membranes used, added growth factors such as GM-CSF, as well as factors secreted from the cells, are retained in the cell compartment for up to 1 week. This new miniaturized hollow-fiber bioreactor offers advantages in tissue engineering by continuous nutrient supply for cells in high density, retention of added or autocrine produced factors, and undisturbed long-term culture in a closed system.