Biofilm development and the dynamics of preferential flow paths in porous media

A two-dimensional pore-scale numerical model was developed to evaluate the dynamics of preferential flow paths in porous media caused by bioclogging. The liquid flow and solute transport through the pore network were coupled with a biofilm model including biomass attachment, growth, decay, lysis, and detachment. Blocking of all but one flow path was obtained under constant liquid inlet flow rate and biomass detachment caused by shear forces only. The stable flow path formed when biofilm detachment balances growth, even with biomass weakened by decay. However, shear forces combined with biomass lysis upon starvation could produce an intermittently shifting location of flow channels. Dynamic flow pathways may also occur when combined liquid shear and pressure forces act on the biofilm. In spite of repeated clogging and unclogging of interconnected pore spaces, the average permeability reached a quasi-constant value. Oscillations in the medium permeability were more pronounced for weaker biofilms.     

The effect of biofilm permeability on bio-clogging of porous media

A 3D Biofilm model, appropriate for complex porous media support structures, is successfully modified such that non‐zero permeability of biofilms structures is enabled. A systematic study is then conducted into the influence of biofilm permeability on overall biomass growth rate. This reveals a significant influence at large biofilm concentrations; even when the permeability of the biomass is 1.25% of that of the free pore space, biomass accumulation increased by a factor of ～3 over 40 h. The effect is shown to be retained when allowing for biomass detachment or erosion as a consequence of adjacent velocity shear. We conclude that biofilm permeability should be included in biofilm models and that further experimental work is required to better describe the link between biofilm permeability and local microstructure. Biotechnol. Bioeng. 2012; 109:1031–1042. © 2011 Wiley Periodicals, Inc.     

Biomass evolution in porous media and its effects on permeability under starvation conditions

The purpose of this study was to understand bacteria profile modification and its applications in subsurface biological operations such as biobarrier formation, in situ bioremediation, and microbial-enhanced oil recovery. Biomass accumulation and evolution in porous media were investigated both experimentally and theoretically. To study both nutrient-rich and carbon-source-depleted conditions, Leuconostoc mesenteroides was chosen because of its rapid growth rate and exopolymer production rate. Porous micromodels were used to study the effects of biomass evolution on the permeability of a porous medium. Bacterial starvation was initiated by switching the feed from a nutrient solution to a buffer solution in order to examine biofilm stability under nutrient-poor conditions. Four different evolution patterns were identified during the nutrient-rich and nutrient-depleted conditions used in the micromodel experiments. In phase I, the permeability of the porous micromodel decreased as a result of biomass accumulation in pore bodies and pore throats. In phase II, starvation conditions were initiated. The depletion of nutrient in the phase II resulted in slower growth of the biofilm causing the permeability to reach a minimum as all the remaining nutrients were consumed. In phase III, permeability began to increase due to biofilm sloughing caused by shear stress. In phase IV, shear stress remained below the critical shear stress for sloughing and the biofilm remained stable for long periods of time during starvation. The critical shear stress for biofilm sloughing provided an indication of biofilm strength. Shear removal of biofilms occurred when shear stress exceeded critical shear stress. A network model was used to describe the biofilm formation phenomenon and the existence of a critical shear stress. Simulations were in qualitative agreement with the experimental results, and demonstrate the existence of a critical shear stress.     

Towards optimum permeability reduction in porous media using biofilm growth simulations

While biological clogging of porous systems can be problematic in numerous processes (e.g., microbial enhanced oil recovery—MEOR), it is targeted during bio‐barrier formation to control sub‐surface pollution plumes in ground water. In this simulation study, constant pressure drop (CPD) and constant volumetric flow rate (CVF) operational modes for nutrient provision for biofilm growth in a porous system are considered with respect to optimum (minimum energy requirement for nutrient provision) permeability reduction for bio‐barrier applications. Biofilm growth is simulated using a Lattice‐Boltzmann (LB) simulation platform complemented with an individual‐based biofilm model (IbM). A biomass detachment technique has been included using a fast marching level set (FMLS) method that models the propagation of the biofilm–liquid interface with a speed proportional to the adjacent velocity shear field. The porous medium permeability reduction is simulated for both operational modes using a range of biofilm strengths. For stronger biofilms, less biomass deposition and energy input are required to reduce the system permeability during CPD operation, whereas CVF is more efficient at reducing the permeability of systems containing weaker biofilms. Biotechnol. Bioeng. 2009;103: 767–779. © 2009 Wiley Periodicals, Inc.     

A new approach to model the spatiotemporal development of biofilm phase in porous media

Bacteria can exist within biofilms that are attached to the solid matrix of a porous medium. Under certain conditions, the biomass can fully occupy the pore space leading to reduced hydraulic conductivity and mass transport. Here, by treating biofilm as a growing, high-viscosity phase, a novel macroscopic approach to model biofilm spatial expansion and its corresponding effects on porous medium hydraulic properties is presented. The separate yet coupled flow of the water and biofilm phases is handled by using relative permeability curves that allow for biofilm movement within the porous medium and bioclogging effects. Fluid flow is governed by Darcy's law and component transport is set by the convection-diffusion equation reaction terms for each component. Here, the system of governing equations is solved by using a commercial multiphase flow reservoir simulator, which is used to validate the model against published laboratory experiments. A comparison of the model and experimental observations reveal that the model provides a reasonable means to predict biomass development in the porous medium. The results reveal that coupled flow of water and movement of biofilm, as described by relative permeability curves, is complex and has a large impact on the development of biomass and consequent bioclogging in the porous medium.     

Modeling of an aerobic biofilm reactor with double-limiting substrate kinetics: bifurcational and dynamical analysis

A mathematical model of an aerobic biofilm reactor is presented to investigate the bifurcational patterns and the dynamical behavior of the reactor as a function of different key operating parameters. Suspended cells and biofilm are assumed to grow according to double limiting kinetics with phenol inhibition (carbon source) and oxygen limitation. The model presented by Russo et al. is extended to embody key features of the phenomenology of the granular‐supported biofilm: biofilm growth and detachment, gas–liquid oxygen transport, phenol, and oxygen uptake by both suspended and immobilized cells, and substrate diffusion into the biofilm. Steady‐state conditions and stability, and local dynamic behavior have been characterized. The multiplicity of steady states and their stability depend on key operating parameter values (dilution rate, gas–liquid mass transfer coefficient, biofilm detachment rate, and inlet substrate concentration). Small changes in the operating conditions may be coupled with a drastic change of the steady‐state scenario with transcritical and saddle‐node bifurcations. The relevance of concentration profiles establishing within the biofilm is also addressed. When the oxygen level in the liquid phase is <10% of the saturation level, the biofilm undergoes oxygen starvation and the active biofilm fraction becomes independent of the dilution rate. © 2011 American Institute of Chemical Engineers Biotechnol. Prog., 2011     

A novel approach to investigate biofilm accumulation and bacterial transport in porous matrices

Knowledge of bacterial transport through, and biofilm growth in, porous media is vitally important in numerous natural and engineered environments. Despite this, porous media systems are generally oversimplified and the local complexity of cell transport, biofilm formation and the effect of biofilm accumulation on flow patterns is lost. In this study, cells of the sulphate-reducing bacterium, Desulfovibrio sp. EX265, accumulated primarily on the leading faces of obstructions and developed into biofilm, which grew to narrow and block pore throats (at a rate of 12 micro m h(-1) in one instance). This pore blocking corresponded to a decrease in permeability from 9.9 to 4.9 Darcy. Biofilm processes were observed in detail and quantitative data were used to describe the rate of biofilm accumulation temporally and spatially. Accumulation in the inlet zone of the micromodel was 10% higher than in the outlet zone and a mean biofilm height of 28.4 micro m was measured in a micromodel with an average pore height of 34.9 microm. Backflow (flow reversal) of fluid was implemented on micromodels blocked with biofilm growth. Although biofilm surface area cover did immediately decrease (approximately 5%), the biofilm quickly re-established and permeability was not significantly affected (9.4 Darcy). These results demonstrate that the glass micromodel used here is an effective tool for in situ analysis and quantification of bacteria in porous media.     

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.     

Flowing biofilms as a transport mechanism for biomass through porous media under laminar and turbulent conditions in a laboratory reactor system

Abstract

Fluid flow has been shown to be important in influencing biofilm morphology and causing biofilms to flow over surfaces in flow cell experiments. However, it is not known whether similar effects may occur in porous media. Generally, it is assumed that the primary transport mechanism for biomass in porous media is through convection, as suspended particulates (cells and flocs) carried by fluid flowing through the interstices. However, the flow of biofilms over the surfaces of soils and sediment particles, may represent an important flux of biomass, and subsequently affect both biological activity and permeability. Mixed species bacterial biofilms were grown in glass flow cells packed with 1 mm diameter glass beads, under laminar or turbulent flow (porous media Reynolds number = 20 and 200 respectively). The morphology and dynamic behavior reflected those of biofilms grown in the open flow cells. The laminar biofilm was relatively uniform and after 23 d had inundated the majority of the pore spaces. Under turbulent flow the biofilm accumulated primarily in protected regions at contact points between the beads and formed streamers that trailed from the leeward face. Both biofilms caused a 2 to 3-fold increase in friction factor and in both cases there were sudden reductions in friction factor followed by rapid recovery, suggesting periodic sloughing and regrowth events. Time-lapse microscopy revealed that under both laminar and turbulent conditions biofilms flowed over the surface of the porous media. In some instances ripple structures formed. The velocity of biofilm flow was on the order of 10 μm h^?1 in the turbulent flow cell and 1.0 μm h^?1 in the laminar flow cell.     

Determination of pollutant diffusion coefficients in naturally formed biofilms using a single tube extractive membrane bioreactor

A novel technique has been used to determine the effective diffusion coefficients for 1,1,2-trichloroethane (TCE), a nonreacting tracer, in biofilms growing on the external surface of a silicone rubber membrane tube during degradation of 1,2-dichloroethane (DCE) by Xanthobacter autotrophicus GJ10 and monochlorobenzene (MCB) by Pseudomonas JS150. Experiments were carried out in a single tube extractive membrane bioreactor (STEMB), whose configuration makes it possible to measure the transmembrane flux of substrates. A video imaging technique (VIT) was employed for in situ biofilm thickness measurement and recording. Diffusion coefficients of TCE in the biofilms and TCE mass transfer coefficients in the liquid films adjacent to the biofilms were determined simultaneously using a resistances-in-series diffusion model. It was found that the flux and overall mass transfer coefficient of TCE decrease with increasing biofilm thickness, showing the importance of biofilm diffusion on the mass transfer process. Similar fluxes were observed for the nonreacting tracer (TCE) and the reactive substrates (MCB or DCE), suggesting that membrane-attached biofilm systems can be rate controlled primarily by substrate diffusion. The TCE diffusion coefficient in the JS150 biofilm appeared to be dependent on biofilm thickness, decreasing markedly for biofilm thicknesses of >1 mm. The values of the TCE diffusion coefficients in the JS150 biofilms <1-mm thick are approximately twice those in water and fall to around 30% of the water value for biofilms >1-mm thick. The TCE diffusion coefficients in the GJ10 biofilms were apparently constant at about the water value. The change in the diffusion coefficient for the JS150 biofilms is attributed to the influence of eddy diffusion and convective flow on transport in the thinner (<1-mm thick) biofilms.     

Flowing biofilms as a transport mechanism for biomass through porous media under laminar and turbulent conditions in a laboratory reactor system

Fluid flow has been shown to be important in influencing biofilm morphology and causing biofilms to flow over surfaces in flow cell experiments. However, it is not known whether similar effects may occur in porous media. Generally, it is assumed that the primary transport mechanism for biomass in porous media is through convection, as suspended particulates (cells and flocs) carried by fluid flowing through the interstices. However, the flow of biofilms over the surfaces of soils and sediment particles, may represent an important flux of biomass, and subsequently affect both biological activity and permeability. Mixed species bacterial biofilms were grown in glass flow cells packed with 1 mm diameter glass beads, under laminar or turbulent flow (porous media Reynolds number = 20 and 200 respectively). The morphology and dynamic behavior reflected those of biofilms grown in the open flow cells. The laminar biofilm was relatively uniform and after 23 d had inundated the majority of the pore spaces. Under turbulent flow the biofilm accumulated primarily in protected regions at contact points between the beads and formed streamers that trailed from the leeward face. Both biofilms caused a 2 to 3-fold increase in friction factor and in both cases there were sudden reductions in friction factor followed by rapid recovery, suggesting periodic sloughing and regrowth events. Time-lapse microscopy revealed that under both laminar and turbulent conditions biofilms flowed over the surface of the porous media. In some instances ripple structures formed. The velocity of biofilm flow was on the order of 10 mum h(-1) in the turbulent flow cell and 1.0 mum h(-1) in the laminar flow cell.     

Magnetic resonance microscopy analysis of advective transport in a biofilm reactor

In this article we present magnetic resonance microscopy (MRM) characterization of the advective transport in a biofilm capillary reactor. The biofilm generates non-axial flows that are up to 20% of the maximum axial velocity. The presence of secondary velocities of this magnitude alters the mass transport in the bioreactor relative to non-biofilm fouled reactors and questions the applicability of empirical mass transfer coefficient approaches. The data are discussed in the context of simulations and models of biofilm transport and conceptual aspects of transport modeling in complex flows are also discussed. The variation in the residence time distribution due to biofilm growth is calculated from the measured propagator of the motion. Dynamical systems methods applied to model fluid mixing in complex flows are indicated as a template for extending mass transport theory to quantitatively incorporate microscale data on the advection field into macroscale mass transfer models.     

Modelling biofilm‐induced formation damage and biocide treatment in subsurface geosystems

Biofilm growth in subsurface porous media, and its treatment with biocides (antimicrobial agents), involves a complex interaction of biogeochemical processes which provide non‐trivial mathematical modelling challenges. Although there are literature reports of mathematical models to evaluate biofilm tolerance to biocides, none of these models have investigated biocide treatment of biofilms growing in interconnected porous media with flow. In this paper, we present a numerical investigation using a pore network model of biofilm growth, formation damage and biocide treatment. The model includes three phases (aqueous, adsorbed biofilm, and solid matrix), a single growth‐limiting nutrient and a single biocide dissolved in the water. Biofilm is assumed to contain a single species of microbe, in which each cell can be a viable persister, a viable non‐persister, or non‐viable (dead). Persisters describe small subpopulation of cells which are tolerant to biocide treatment. Biofilm tolerance to biocide treatment is regulated by persister cells and includes ‘innate’ and ‘biocide‐induced’ factors. Simulations demonstrate that biofilm tolerance to biocides can increase with biofilm maturity, and that biocide treatment alone does not reverse biofilm‐induced formation damage. Also, a successful application of biological permeability conformance treatment involving geologic layers with flow communication is more complicated than simply engineering the attachment of biofilm‐forming cells at desired sites.     

Experimental and theoretical evidence for convective nutrient transport in an immobilized cell support.

Even though immobilized-cell reactors possess several engineering advantages over free-cell reactors, their full potential has not been realized because mass transfer often limits the rate of nutrient supply and product removal from immobilized cell supports. We studied the interaction between mass transfer and reaction kinetics in the anaerobic conversion of glucose to CO2 and ethanol by yeast immobilized in a porous rotating disk on the agitator shaft of a conventional CSTR. A Sherwood number correlation was used to show that external mass-transfer resistances were negligible under typical operating conditions. The modulus of Weisz based on observable reaction parameters was used to gauge the importance of pore diffusion limitations. Under conditions for which significant pore diffusion effects and hence low effectiveness factors (eta = ca. 0.1) would be predicted, the observed reaction rates were much higher than expected (eta = ca. 1), suggesting that pore diffusion limitations were at least partially relieved by convective transport of glucose into the support. Two possible mechanisms of convective transport are discussed. We hypothesize that gas evolution was responsible for the convective enhancement of glucose supply.     

Biofilm development in a membrane-aerated biofilm reactor: effect of flow velocity on performance

The effect of liquid flow velocity on biofilm development in a membrane-aerated biofilm reactor was investigated both by mathematical modeling and by experiment, using Vibrio natriegens as a test organism and acetate as carbon substrate. It was shown that velocity influenced mass transfer in the diffusion boundary layer, the biomass detachment rate from the biofilm, and the maximum biofilm thickness attained. Values of the overall mass transfer coefficient of a tracer through the diffusion boundary layer, the biofilm, and the membrane were shown to be identical during different experiments at the maximum biofilm thickness. Comparison of the results with published values of this parameter in membrane attached biofilms showed a similar trend. Therefore, it was postulated that this result might indicate the mechanism that determines the maximum biofilm thickness in membrane attached biofilms. In a series of experiments, where conditions were set so that the active layer of the membrane attached biofilm was located close to the membrane biofilm interface, it was shown that the most critical effect on process performance was the effect of velocity on biofilm structure. Biofilm thickness and effective diffusivity influenced reaction and diffusion in a complex manner such that the yield of biomass on acetate was highly variable. Consideration of endogenous respiration in the mathematical model was validated by direct experimental measurements of yield coefficients. Good agreement between experimental measurements of acetate and oxygen uptake rates and their prediction by the mathematical model was achieved.     

Microelectrode measurements of local mass transport rates in heterogeneous biofilms

Microelectrodes were used to measure oxygen profiles and local mass transfer coefficient profiles in biofilm clusters and interstitial voids. Both profiles were measured at the same location in the biofilm. From the oxygen profile, the effective diffusive boundary layer thickness (DBL) was determined. The local mass transfer coefficient profiles provided information about the nature of mass transport near and within the biofilm. All profiles were measured at three different average flow velocities, 0.62, 1.53, and 2.60 cm sec-1, to determine the influence of flow velocity on mass transport. Convective mass transport was active near the biofilm/liquid interface and in the upper layers of the biofilm, independent of biofilm thickness and flow velocity. The DBL varied strongly between locations for the same flow velocities. Oxygen and local mass transfer coefficient profiles collected through a 70 micrometer thick cluster revealed that a cluster of that thickness did not present any significant mass transport resistance. In a 350 micrometer thick biofilm cluster, however, the local mass transfer coefficient decreased gradually to very low values near the substratum. This was hypothetically attributed to the decreasing effective diffusivity in deeper layers of biofilms. Interstitial voids between clusters did not seem to influence the local mass transfer coefficients significantly for flow velocities of 1.53 and 2.60 cm sec-1. At a flow velocity of 0.62 cm sec-1, interstitial voids visibly decreased the local mass transfer coefficient near the bottom.     

Diffusion and reaction within porous packing media: a phenomenological model

A phenomenological model has been developed to describe biomass distribution and substrate depletion in porous diatomaceous earth (DE) pellets colonized by Pseudomonas aeruginosa. The essential features of the model are diffusion, attachment and detachment to/from pore walls of the biomass, diffusion of substrate within the pellet, and external mass transfer of both substrate and biomass in the bulk fluid of a packed bed containing the pellets. A bench-scale reactor filled with DE pellets was inoculated with P. aeruginosa and operated in plug flow without recycle using a feed containing glucose as the limiting nutrient. Steady-state effluent glucose concentrations were measured at various residence times, and biomass distribution within the pellet was measured at the lowest residence time. In the model, microorganism/substrate kinetics and mass transfer characteristics were predicted from the literature. Only the attachment and detachment parameters were treated as unknowns, and were determined by fitting biomass distribution data within the pellets to the mathematical model. The rate-limiting step in substrate conversion was determined to be internal mass transfer resistance; external mass transfer resistance and microbial kinetic limitations were found to be nearly negligible. Only the outer 5% of the pellets contributed to substrate conversion. (c) 1993 Wiley & Sons, Inc.     

Xanthan degragation by biofilm in porous media

Biological activity in oil reservoirs can cause significant problems such as souring and plugging. This study focuses on the problem of polymer degradation and permeability reduction due to biofilm formation during polymer injection for improved oil recovery. Polymers are included in injection fluids to increase their viscosity. Results relating biological processes and polymer degradation to fluid‐dynamic conditions in a laboratory model porous medium are presented.

A transparent flow cell with an etched two‐dimensional network of pores served as a model porous medium. A sterile xanthan polymer and natural sea water solution were continuously injected into the porous medium. A bacterial culture capable of xanthan degradation was introduced into the cell by a single injection. Some of the cells from this culture attached to the pore walls forming an immobile bacterial culture, termed biofilm. The development of this biofilm, its xanthan degradation and its effect on permeability were measured.

The effects of injection rate and rate transitions were analyzed. Injection fluid viscosity was reduced by 30% after 5 min flow through the porous medium at the maximum steady state degradation rate observed. Permeability was significantly reduced by the xanthan degrading biofilm, causing an increase in pressure drop through the porous medium of up to 80%. Polymer injection in oil reservoirs may, therefore, have negative effects on oil recovery, unless efficient biofouling control is applied. The methodology presented may serve as a tool in the development of biofouling control measures in porous media.     

Dynamic volume-averaged model of heat and mass transport within a compost biofilter: I. Model development

Successful, long-term operation of a biofilter system depends on maintaining a suitable biofilm environment within a porous medium reactor. In this article a mathematical study was conducted of the spatial and temporal changes of biofilter performance due to interphase heat and mass transport. The method of volume averaging was used to spatially smooth the three-phase (solid, liquid, and gas) conservation equations over the biofilter domain. The packing medium was assumed to be inert, removing the solid phase mass continuity equation from the system. The finite element method was used to integrate the resulting nonlinear-coupled partial differential equations, tracking eight state variables: temperature, water vapor, dry air, liquid water, biofilm, gas and liquid phase organic pollutant, and nutrient densities, through time and space. A multiphase, gas and liquid flow model was adapted to the biofilter model from previous studies of unsaturated groundwater flow. Newton's method accelerated by an LU direct solver was used to iterate the model for solutions. Effects of packing media on performance were investigated to illustrate the utility of the model. The moisture dynamics and nutrient cycling are presented in Part II of this article.