Formation and growth of heterotrophic aerobic biofilms on small suspended particles in airlift reactors

In this article, the conditions for aerobic biofilm formation on suspended particles, the dynamics of biofilm formation, and the biomass production during the start-up of a Biofilm Airlift Suspension reactor (BAS reactor) have been studied. The dynamics of biofilm formation during start up in the biofilm airlift suspension reactor follows three consecutive stages: bare carrier, microcolonies or patchy biofilms on the carrier, and biofilms completely covering the carrier. The effect of hydraulic retention time and of substrate loading rate on the formation of biofilms were investigated. To obtain in a BAS reactor a high biomass concentration and predominantly continuous biofilms, which completely surround the carrier, the hydraulic retention time must be shorter than the inverse of the maximum growth rate of the suspended bacteria. At longer hydraulic retention times, a low amount of attached biomass can be present on the carrier material as patchy biofilms. During the start-up at short hydraulic retention times the bare carrier concentration decreases, the amount of biomass per biofilm particle remains constant, and biomass increase in the reactor is due to increasing numbers of biofilm particles. The substrate surface loading rate has effect only on the amount of biomass on the biofilm particle. A higher surface load leads to a thicker biofilm.A strong nonlinear increase of the concentration of attached biomass in time was observed. This can be explained by a decreased abrasion of the biofilm particles due to the decreasing concentration of bare carriers. The detachment rate per biofilm area during the start-up is independent of the substrate loading rate, but depends strongly upon the bare carrier concentration.The Pirt-maintenance concept is applicable to BAS reactors. Surplus biomass production is diminished at high biomass concentrations. The average maximal yield of biomass on substrate during the experiments presented in this article was 0.44 +/- 0.08 C-mol/C-mol, the maintenance value 0.019 +/- 0.012 C-mol/(C-mol h). The lowest actual biomass yield measured in this study was 0.15 C-mol/C-mol. (c) 1994 John Wiley & Sons, Inc.     

Influence of biomass production and detachment forces on biofilm structures in a biofilm airlift suspension reactor

The influence of process conditions (substrate loading rate and detachment force) on the structure of biofilms grown on basalt particles in a Biofilm Airlift Suspension (BAS) reactor was studied. The structure of the biofilms (density, surface shape, and thickness) and microbial characteristics (biomass yield) were investigated at substrate loading rates of 5, 10, 15, and 20 kg COD/m3. day with basalt concentrations of 60 g/L, 150 g/L, and 250 g/L. The basalt concentration determines the number of biofilm particles in steady state, which is the main determining factor for the biofilm detachment in these systems. In total, 12 experimental runs were performed. A high biofilm density (up to 67 g/L) and a high biomass concentration was observed at high detachment forces. The higher biomass content is associated with a lower biomass substrate loading rate and therefore with a lower biomass yield (from 0.4 down to 0.12 gbiomass/gacetate). Contrary to general beliefs, the observed biomass detachment decreased with increasing detachment force. In addition, smoother (fewer protuberances), denser and thinner compact biofilms were obtained when the biomass surface production rate decreased and/or the detachment force increased. These observations confirmed a hypothesis, postulated earlier by Van Loosdrecht et al. (1995b), that the balance between biofilm substrate surface loading (proportional to biomass surface production rate, when biomass yield is constant) and detachment force determines the biofilm structure. When detachment forces are relatively high only a patchy biofilm will develop, whereas at low detachment forces, the biofilm becomes highly heterogeneous with many pores and protuberances. With the right balance, smooth, dense and stable biofilms can be obtained. Copyright 1998 John Wiley & Sons, Inc.     

Influence of detachment,substrate loading and reactor scale on the formation of biofilms in airlift reactors

 For a stable and reliable operation of the biofilm airlift suspension reactor (BAS reactor) means to control biomass concentration, biofilm thickness and biofilm morphology are required. For this reason, the influence of applied detachment forces and surface substrate loading on the formation of heterotrophic biofilms in laboratory-scale BAS reactors was studied. Detachment forces were altered by variation of the initial bare carrier concentration or the superficial air velocity. In addition, the dynamics of biofilm formation during start-up of a full scale BAS reactor (300 m³) was monitored and compared with the laboratory-scale start-up (3 l). This study shows that the biofilm morphology and strength were influenced to a large extent by the surface substrate loading and applied detachment forces. A moderate surface substrate loading and a high detachment force yielded smooth and strong biofilms. The combination of a high surface substrate loading and low detachment forces did lead to rough biofilms, but did not lead to the expected high amount of biomass on the carrier, apparently because of the formation of weaker biofilms. The strength of the bio-films appeared to be related to the detachment forces applied during biofilm formation, in combination with the surface substrate loading. The biofilm morphology and biomass on carrier in the BAS reactor can be controlled using the carrier concentration, substrate loading rate and the superficial air velocity as parameters. The dynamics of biofilm formation during the start-up of a full-scale BAS reactor proved to be similar to heterotrophic biofilm formation in laboratory-scale reactors. This indicates that a model system on the laboratory scale can successfully be applied to predict dynamic phenomena in the full-scale reactor. Received: 31 March 1995/Received revision: 11 August 1995/Accepted: 22 August 1995     

A new method to solve a non-steady-state multispecies biofilm model

A transient multispecies model for quantifying microbial space competition in biofilm is derived from existing models, introducing a new approach to biomass detachment modelling. This model includes inert biomass, substrate diffusion and utilization rate within the biofilm and diffusional layers. It predicts the evolution of biofilm thickness, bulk substrate concentration, species distribution and substrate concentration within the biofilm. A zero-dimensional transient model is described. Its steady-state solution is used to set up initial conditions of the one-dimensional model and case computation towards steady-state solution. Some numerical tools have been developed, enabling fast computation on microcomputers. Simulations show the validity of a zero-dimensional model and perturbated systems are also simulated. Simulations with experimental data give acceptable results.     

A general description of detachment for multidimensional modelling of biofilms

A general method for describing biomass detachment in multidimensional biofilm modelling is introduced. Biomass losses from processes acting on the entire surface of the biofilm, such as erosion, are modelled using a continuous detachment speed function F(det). Discrete detachment events, i.e. sloughing, are implicitly derived from simulations. The method is flexible to allow F(det) to take several forms, including expressions dependent on any state variables such as the local biofilm density. This methodology for biomass detachment was integrated with multidimensional (2D and 3D) particle-based multispecies biofilm models by using a novel application of the level set method. Application of the method is illustrated by trends in the dynamics of biofilms structure and activity derived from simulations performed on a simple model considering uniform biomass (case study I) and a model discriminating biomass composition in heterotrophic active mass, extracellular polymeric substances (EPS) and inert mass (case study II). Results from case study I demonstrate the effect of applied detachment forces as a fundamental factor influencing steady-state biofilm activity and structure. Trends from experimental observations reported in literature were correctly described. For example, simulation results indicated that biomass sloughing is reduced when erosion forces are increased. Case study II illustrates the application of the detachment methodology to systems with non-uniform biomass composition. Simulations carried out at different bulk concentrations of substrate show changes in biofilm structure (in terms of shape, density and spatial distribution of biomass components) and activity (in terms of oxygen and substrate consumption) as a consequence of either oxygen-limited or substrate-limited growth.     

Effect of diffusive and convective substrate transport on biofilm structure formation: a two-dimensional modeling study

A two-dimensional model for quantitative evaluation of the effect of convective and diffusive substrate transport on biofilm heterogeneity was developed. The model includes flow computation around the irregular biofilm surface, substrate mass transfer by convection and diffusion, biomass growth, and biomass spreading. It was found that in the absence of detachment, biofilm heterogeneity is mainly determined by internal mass transfer rate of substrates and by the initial percentage of carrier-surface colonization. Model predictions show that biofilm structures with highly irregular surface develop in the mass transfer-limited regime. As the nutrient availability increases, there is a gradual shift toward compact and smooth biofilms. A smaller fraction of colonized carrier surface leads to a patchy biofilm. Biofilm surface irregularity and deep vertical channels are, in this case, caused by the inability of the colonies to spread over the whole substratum surface. The maximum substrate flux to the biofilm was greatly influenced by both internal and external mass transfer rates, but not affected by the inoculation density. In general, results of the present model were similar to those obtained by a simple diffusion-reaction-growth model.     

Effects of stratification in a fluidized bed bioreactor during treatment of metalworking wastewater

During wastewater treatment, biofilm-coated sand particles stratified in a fluidized bed bioreactor (FBB); particles coated by thicker biofilm segregated toward the top of the bed. Stratification was so well developed that at least two co-existing regions of significantly different mean biofilm thickness were visually distinct within the operating FBB. The observed stratification is attributed to differences in forces of drag, buoyancy, shear, and collisional impact, as well as differences of collision rate within the different regions. Particles with thick biofilm (thickness >100 μm) near the top of the bed consumed substrate at significantly lower rates per unit biomass than particles with thin biofilm (10-20 μm) near the bottom of the bed, thereby suggesting that substrate mass-transfer resistance through biofilm may limit biodegradation rates in the upper portion of the FBB. Large agglomerates of biomass floc and sand, which formed at the top of the fluidized bed, and sand particles with thick biofilm were susceptible to washout from the FBB, causing operational and treatment instability. Radial injection of supplemental liquid feed near the top of the bed increased shear and mixing, thereby preventing formation and washout of agglomerates and thickly coated sand particles. Supplemental liquid injection caused the mean specific biomass loading on the sand to increase and also increased the total biomass inventory in the FBB. Rates of biodegradation in the FBB appeared to be limited by penetration of substrates into the biofilm and absorption of oxygen from air into the wastewater. Copyright 1999 John Wiley & Sons, Inc.     

Toluene degradation kinetics for planktonic and biofilm-grown cells of Pseudomonas putida 54G

Toluene degradation kinetics by biofilm and planktonic cells of Pseudomonas putida 54G were compared in this study. Batch degradation of (14)C toluene was used to evaluate kinetic parameters for planktonic cells. The kinetic parameters determined for toluene degradation were: specific growth rate, mu(max) = 10.08 +/- 1.2/day; half-saturation constant, K(S) = 3.98 +/- 1.28 mg/L; substrate inhibition constant, K(I) = 42.78 +/- 3.87 mg/L. Biofilm cells, grown on ceramic rings in vapor phase bioreactors, were removed and suspended in batch cultures to calculate (14)C toluene degradation rates. Specific activities measured for planktonic and biofilm cells were similar based on toluene degrading cells and total biomass. Long-term toluene exposure reduced specific activities that were based on total biomass for both biofilm and planktonic cells. These results suggest that long-term toluene exposure caused a large portion of the biomass to become inactive, even though the biofilm was not substrate limited. Conversely, specific activities based on numbers of toluene-culturable cells were comparable for both biofilm and planktonically grown cultures. Planktonic cell kinetics are often used in bioreactor models to model substrate degradation and growth of bacteria in biofilms, a procedure we found to be appropriate for this organism. For superior bioreactor design, however, changes in cellular activity that occur during biofilm development should be investigated under conditions relevant to reactor operation before predictive models for bioreactor systems are developed. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 53: 535-546, 1997.     

Development and experimental evaluation of a steady-state, multispecies biofilm model

A steady-state model for quantifying the space competition in multispecies biofilms is developed. The model includes multiple active species, inert biomass, substrate utilization and diffusion within the biofilm, external mass transport, and detachment phenomena. It predicts the steady-state values of biofilm thickness, species distribution, and substrate fluxes. An experimental evaluation is carried out in completely mixed biofilm reactors in which slow-growing nitrifying bacteria compete with acetate-utilizing heterotrophs. The experimental results show that the model successfully describes the space competition. In particular, increasing acetate concentrations causes NH(4) (+)-N fluxes to decrease, because nitrifiers are forced deeper into the biofilm, where they experience greater mass-transport resistance.     

Biodegradation Kinetic Studies on Phenol in Internal Draft Tube (Inverse Fluidized Bed) Biofilm Reactor Using Pseudomonas fluorescens: Performance Evaluation of Biofilm and Biomass Characteristics

The bioremediation potential of Pseudomonas fluorescens was studied in an internal draft tube (inverse fluidized bed) biofilm reactor (IDTBR) under batch recirculation conditions using synthetic phenol of various concentrations (400, 600, 800, 1000, and 1200 mg/L). The performance of IDTBR was investigated and the characteristics of biomass and biofilm were determined by evaluating biofilm dry density and thickness, bioparticle density, suspended and attached biomass concentration, chemical oxygen demand, and phenol removal efficiency. Biodegradation kinetics had been studied for the suspended biomass culture and biofilm systems. Suspended biomass followed substrate inhibition kinetics, and the experimental data fitted well with the Haldane model. The correlation coefficient, R ², and root-mean-square error (RMSE) obtained for the Haldane model with respect to specific growth rate were .9389 and .00729, respectively, and with respect to specific phenol consumption rate were .9259 and .00972, respectively. It was also observed experimentally that biofilm overcame substrate inhibition effect and fitted the same to the Monod model (R ² = .9831, RMSE = .00884 for specific growth rate and R ² = .9686, RMSE = .00912 for specific phenol consumption rate).     

A kinetic model for suspended and attached growth of a defined mixed culture

Kinetic experiments were carried out in a semicontinuous wastewater treatment process called self-cycling fermentation (SCF) using a defined mixed culture and various concentrations of synthetic brewery wastewater. The same consortium, which had been previously identified as Acinetobacter sp., Enterobacter sp., and Candida sp., were used in these experiments. The overall rate of substrate removal was attributable to both suspended microbes and the biofilm that formed during the treatment process. A rate expression was developed for the SCF system for a range of synthetic wastewaters containing glucose and various initial concentrations of ethanol and maltose. The data indicated that substrate removal by the suspended cells was directly related to the biomass concentration. However, substrate removal by the biofilm was apparently not affected by the biofilm thickness and was a function of substrate concentration only.     

A theoretical study on the effect of surface roughness on mass transport and transformation in biofilms

This modeling study evaluates the influence of biofilm geometrical characteristics on substrate mass transfer and conversion rates. A spatially two-dimensional model was used to compute laminar fluid flow, substrate mass transport, and conversion in irregularly shaped biofilms. The flow velocity above the biofilm surface was varied over 3 orders of magnitude. Numerical results show that increased biofilm roughness does not necessarily lead to an enhancement of either conversion rates or external mass transfer. The average mass transfer coefficient and Sherwood numbers were found to decrease almost linearly with biofilm area enlargement in the flow regime tested. The influence of flow, biofilm geometry and biofilm activity on external mass transfer could be quantified by Sh-Re correlations. The effect of biofilm surface roughness was incorporated in this correlation via area enlargement. Conversion rates could be best correlated to biofilm compactness. The more compact the biofilm, the higher the global conversion rate of substrate. Although an increase of bulk fluid velocity showed a large effect on mass transfer coefficients, the global substrate conversion rate per carrier area was less affected. If only diffusion occurs in pores and channels, then rough biofilms behave as if they were compact but having less biomass activity. In spite of the fact that the real biofilm area is increased due to roughness, the effective mass transfer area is actually decreased because only biofilm peaks receive substrate. This can be explained by the fact that in the absence of normal convection in the biofilm valleys, the substrate gradients are still largely perpendicular to the carrier. Even in the cases where convective transport dominates the external mass transfer process, roughness could lead to decreased conversion rates. The results of this study clearly indicate that only evaluation of overall conversion rates or mass fluxes can describe the correct biofilm conversion, whereas interpretation of local concentration or flow measurements as such might easily lead to erroneous conclusions.     

Biofilm development in a membrane-aerated biofilm reactor: effect of intra-membrane oxygen pressure on performance

The effect on intra-membrane oxygen pressure at a constant carbon substrate loading rate on the development of biofilms of Vibrio natrigens in a membrane aerated biofilm reactor (MABR) was investigated experimentally and by mathematical modelling. A recently reported technique (Zhang et al., 1998. Biotechnol. Bioeng. 59: 80-89) for the in situ measurement of the substrate diffusion coefficients in a growing biofilm and the mass transfer coefficients in the boundary layer at the biofilm liquid interface was used. This aided the study of the effect of the heterogeneous biofilm structure and also improved the reliability of the model predictions. The different intra-membrane oxygen pressures used, 12.5, 25 and 50 kPa, with acetate as the carbon substrate, showed a marked effect on the initial biofilm growth rate, on acetate removal rate, particularly in thick biofilms and on biofilm structure. The model predicted the substrate limitation regimes, the location of the active biomass layer within the biofilms and the trends in oxygen uptake rate through the membrane into the biofilms. During the development of the biofilms, the biofilm thickness and the intra-membrane oxygen pressure were found to be the most important parameters influencing the MABR performance while the effect of biofilm structure was less marked.     

Pressure drop increase by biofilm accumulation in spiral wound RO and NF membrane systems: role of substrate concentration,flow velocity,substrate load and flow direction

In an earlier study, it was shown that biofouling predominantly is a feed spacer channel problem. In this article, pressure drop development and biofilm accumulation in membrane fouling simulators have been studied without permeate production as a function of the process parameters substrate concentration, linear flow velocity, substrate load and flow direction. At the applied substrate concentration range, 100–400 μg l^?1 as acetate carbon, a higher concentration caused a faster and greater pressure drop increase and a greater accumulation of biomass. Within the range of linear flow velocities as applied in practice, a higher linear flow velocity resulted in a higher initial pressure drop in addition to a more rapid and greater pressure drop increase and biomass accumulation. Reduction of the linear flow velocity resulted in an instantaneous reduction of the pressure drop caused by the accumulated biomass, without changing the biofilm concentration. A higher substrate load (product of substrate concentration and flow velocity) was related to biomass accumulation. The effect of the same amount of accumulated biomass on the pressure drop increase was related to the linear flow velocity. A decrease of substrate load caused a gradual decline in time of both biomass concentration and pressure drop increase. It was concluded that the pressure drop increase over spiral wound reverse osmosis (RO) and nanofiltration (NF) membrane systems can be reduced by lowering both substrate load and linear flow velocity. There is a need for RO and NF systems with a low pressure drop increase irrespective of the biomass formation. Current efforts to control biofouling of spiral wound membranes focus in addition to pretreatment on membrane improvement. According to these authors, adaptation of the hydrodynamics, spacers and pressure vessel configuration offer promising alternatives. Additional approaches may be replacing heavily biofouled elements and flow direction reversal.     

Modeling bisubstrate removal by biofilms

A bisubstrate secondary utilization model is based on the concept that an individual substrate can be utilized not only by the biomass by its utilization but also by the biomass made from the utilization of the other substrate. When substrate concentrations are low, a key factor is having sufficient substrate to initiate biofilm growth. Modeling results for three characteristic cases demonstrate that satisfying a total S(min) concentration for a bisubstrate system is the necessary condition for initiating biofilm growth and simultaneous utilization of both substrates. Because having more than one substrate supporting biofilm growth enhances the removal of each compound, the utilization rate of a specific compound can be increased by the concentration of other compounds, and the total S(min) concentration can be less than the weighted average of individual S(min) values.     

Specific ATP and specific oxygen uptake rate in immobilized cell aggregates: Experimental results and theoretical analysis using a structured model of immobilized cell growth

In this work the quality and activity of immobilized Beneckea natriegens have been measured using Specific ATP (SATP) [mg (ATP) g(-1) (dry biomass)] and Specific Oxygen Uptake Rate (SOUR) [mg O(2) g(-1) (dry biomass) h(-1)]. The cells were grown in a 3 L Three Phase Air Lift Bioreactor (TPALB) and were immobilized on diatomaceous earth (silica) support particles; sterile conditions were employed during the experiment, with n-propanol as the sole carbon source. Two sets of experiments were performed, one with varying dilution rate and the other with varying inlet substrate concentration. The average SATP and SOUR of the immobilized biomass was 3-4 times lower than the values obtained for suspended Beneckea natriegens growing at its maximum growth rate. The suspended biomass present in the TPALB was generated primarily through attrition from the outer layer of the biofilm, and maintained higher levels of SATP and SOUR than the immobilized biomass. This result indicates that the immobilized biomass quality/activity is higher at the external layer of the biofilm. A structured model in which biomass is divided into two compartments, active and inert, was used to describe the experimental results. This model predicts the biomass quality/activity and substrate concentration distributions in the biofilm. These distributions were integrated to give overall values of SATP and SOUR for the immobilized biomass which compared favorably with experimental data. The primary implication of the results is that the location of immobilized biomass in the biofilm affects its biocatalytic activity, and should be taken into account when modeling immobilized biomass bioreactors. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 55: 660-673, 1997.     

A framework for multidimensional modelling of activity and structure of multispecies biofilms

Concepts from previous biofilm models were integrated to create a framework for the implementation of multidimensional (2D and 3D) multispecies biofilm models. The framework is here described at three levels: (i) mathematical representation of the processes involved in biofilm formation, (ii) numerical implementation into a computer program (freely available from our website http://www.biofilms.bt.tudelft.nl/frameworkMaterial) and (iii) using the program for the creation of biofilm models with multiple bacterial and solute species. An improved version of the individual-based modelling (IbM) that allows structured biomass was used. In this approach biomass composition may be discriminated into any number of particulate species, including extracellular polymeric substances (EPS) for which specific functionality was included. Detachment is also included, described as occurring at the biofilm surface with variable local rates derived from functions of state variables. The application of this modelling framework to a multispecies system with structured biomass is illustrated in a case study where the competition between an organism capable of accumulating polyhydroxybutyrate (PHB, an internal storage compound) and an EPS-producing organism in a two-species biofilm is analysed. Results illustrate that biofilms enriched in PHB-producing organisms may be obtained by supplying substrate intermittently in feast/famine cycles.     

Biofilm model calibration and microbial diversity study using Monte Carlo simulations

Mathematical models are useful tools for studying and exploring biological conversion processes as well as microbial competition in biological treatment processes. A single‐species biofilm model was used to describe biofilm reactor operation at three different hydraulic retention times (HRT). The single‐species biofilm model was calibrated with sparse experimental data using the Monte Carlo filtering method. This calibrated single‐species biofilm model was then extended to a multi‐species model considering 10 different heterotrophic bacteria. The aim was to study microbial diversity in bulk phase biomass and biofilm, as well as the competition between suspended and attached biomass. At steady state and independently of the HRT, Monte Carlo simulations resulted only in one unique dominating bacterial species for suspended and attached biomass. The dominating bacterial species was determined by the highest specific substrate affinity (ratio of µ/K_S). At a short HRT of 20 min, the structure of the microbial community in the bulk liquid was determined by biomass detached from the biofilm. At a long HRT of 8 h, both biofilm detachment and microbial growth in the bulk liquid influenced the microbial community distribution. Biotechnol. Bioeng. 2013; 110: 1323–1332. © 2012 Wiley Periodicals, Inc.