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.     

Pore‐network modeling of biofilm evolution in porous media

The influence of bacterial biomass on hydraulic properties of porous media (bioclogging) has been explored as a viable means for optimizing subsurface bioremediation and microbial enhanced oil recovery. In this study, we present a pore network simulator for modeling biofilm evolution in porous media including hydrodynamics and nutrient transport based on coupling of advection transport with Fickian diffusion and a reaction term to account for nutrient consumption. Biofilm has non‐zero permeability permitting liquid flow and transport through the biofilm itself. To handle simultaneous mass transfer in both liquid and biofilm in a pore element, a dual‐diffusion mass transfer model is introduced. The influence of nutrient limitation on predicted results is explored. Nutrient concentration in the network is affected by diffusion coefficient for nutrient transfer across biofilm (compared to water/water diffusion coefficient) under advection dominated transport, represented by mass transport Péclet number >1. The model correctly predicts a dependence of rate of biomass accumulation on inlet concentration. Poor network connectivity shows a significantly large reduction of permeability, for a small biomass pore volume. Biotechnol. Bioeng. 2011;108: 2413–2423. © 2011 Wiley Periodicals, Inc.     

Comparison of bioclogging effects in saturated porous media within one- and two-dimensional flow systems

The accumulation of microbial biomass in a porous medium can lead to a reduction of pore space and an associated decrease in the hydraulic conductivity of the medium – an effect called bioclogging. This phenomenon may occur in several natural and engineered subsurface systems and can be relevant in fields ranging from contaminant hydrology to civil and environmental engineering, as well as for enhanced oil recovery. During the last decades bioclogging has been studied in various laboratory and theoretical studies. Most of these studies considered only one-dimensional flow fields inside a porous medium. Although these studies provided valuable information on bioclogging and factors controlling it, recent studies showed that an extrapolation of these results to multi-dimensional flow fields is not straight forward. This paper reviews the experimental results obtained for one- and two-dimensional flow fields and compares the modeling results obtained using different conceptualizations of the pore space.     

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.     

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.     

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.     

Real time monitoring of biofilm development under flow conditions in porous media

Biofilm growth can impact the effectiveness of industrial processes that involve porous media. To better understand and characterize how biofilms develop and affect hydraulic properties in porous media, both spatial and temporal development of biofilms under flow conditions was investigated in a translucent porous medium by using Pseudomonas fluorescens HK44, a bacterial strain genetically engineered to luminesce in the presence of an induction agent. Real-time visualization of luminescent biofilm growth patterns under constant pressure conditions was captured using a CCD camera. Images obtained over 8 days revealed that variations in bioluminescence intensity could be correlated to biofilm cell density and hydraulic conductivity. These results were used to develop a real-time imaging method to study the dynamic behavior of biofilm evolution in a porous medium, thereby providing a new tool to investigate the impact of biological fouling in porous media under flow conditions.     

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.     

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.     

Real time monitoring of biofilm development under flow conditions in porous media

Biofilm growth can impact the effectiveness of industrial processes that involve porous media. To better understand and characterize how biofilms develop and affect hydraulic properties in porous media, both spatial and temporal development of biofilms under flow conditions was investigated in a translucent porous medium by using Pseudomonas fluorescens HK44, a bacterial strain genetically engineered to luminesce in the presence of an induction agent. Real-time visualization of luminescent biofilm growth patterns under constant pressure conditions was captured using a CCD camera. Images obtained over 8 days revealed that variations in bioluminescence intensity could be correlated to biofilm cell density and hydraulic conductivity. These results were used to develop a real-time imaging method to study the dynamic behavior of biofilm evolution in a porous medium, thereby providing a new tool to investigate the impact of biological fouling in porous media under flow conditions.     

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.     

Permeability of a growing biofilm in a porous media fluid flow analyzed by magnetic resonance displacement‐relaxation correlations

Biofilm growth in porous media is difficult to study non‐invasively due to the opaqueness and heterogeneity of the systems. Magnetic resonance is utilized to non‐invasively study water dynamics within porous media. Displacement‐relaxation correlation experiments were performed on fluid flow during biofilm growth in a model porous media of mono‐dispersed polystyrene beads. The spin–spin T₂ magnetic relaxation distinguishes between the biofilm phase and bulk fluid phase due to water–biopolymer interactions present in the biofilm, and the flow dynamics are measured using PGSE NMR experiments. By correlating these two measurements, the effects of biofilm growth on the fluid dynamics can be separated into a detailed analysis of both the biofilm phase and the fluid phase simultaneously within the same experiment. Within the displacement resolution of these experiments, no convective flow was measured through the biomass. An increased amount of longitudinal hydrodynamic dispersion indicates increased hydrodynamic mixing due to fluid channeling caused by biofilm growth. The effect of different biofilm growth conditions was measured by varying the strength of the bacterial growth medium. Biotechnol. Bioeng. 2013; 110: 1366–1375. © 2012 Wiley Periodicals, Inc.     

A biphasic model for micro-indentation of a hydrogel-based contact lens

The stiffness and hydraulic permeability of soft contact lenses may influence its clinical performance, e.g., on-eye movement, fitting, and wettability, and may be related to the occurrence of complications; e.g., lesions. It is therefore important to determine these properties in the design of comfortable contact lenses. Micro-indentation provides a nondestructive means of measuring mechanical properties of soft, hydrated contact lenses. However, certain geometrical and material considerations must be taken into account when analyzing output force-displacement (F-D) data. Rather than solely having a solid response, mechanical behavior of hydrogel contact lenses can be described as the coupled interaction between fluid transport through pores and solid matrix deformation. In addition, indentation of thin membranes ( approximately 100 microm) requires special consideration of boundary conditions at lens surfaces and at the indenter contact region. In this study, a biphasic finite element model was developed to simulate the micro-indentation of a hydrogel contact lens. The model accounts for a curved, thin hydrogel membrane supported on an impermeable mold. A time-varying boundary condition was implemented to model the contact interface between the impermeable spherical indenter and the lens. Parametric studies varying the indentation velocities and hydraulic permeability show F-D curves have a sensitive region outside of which the force response reaches asymptotic limits governed by either the solid matrix (slow indentation velocity, large permeability) or the fluid transport (high indentation velocity, low permeability). Using these results, biphasic properties (Young's modulus and hydraulic permeability) were estimated by fitting model results to F-D curves obtained at multiple indentation velocities (1.2 and 20 microm/s). Fitting to micro-indentation tests of Etafilcon A resulted in an estimated permeability range of 1.0 x 10(-15) to 5.0 x 10(-15) m(4)N s and Young's modulus range of 130 to 170 kPa.     

Reduction of porous media permeability from in situ Leuconostoc mesenteroides growth and dextran production

In situ growth of bacteria in a porous medium can alter the permeability of that media. This article reveals that the rate of permeability alteration can be controlled by the inoculation strategy, nutrient concentrations, and injection rates. Based on experimental observations a phenomenological model has been developed to describe the inoculation of the porous medium, the in situ growth of bacteria, and the permeability decline of the porous medium. This model consists of two phases that describe the bacteria in the porous medium: (1) the nongrowth phase in which cell transport and retention are occurring; and (2) the growth phase in which the retained cells grow and plug the porous media. Transition from the transport phase to the growth phase is governed by the growth lag time of the cells within the porous medium. The importance of the inoculum injection strategy and the nutrient injection strategy is illustrated by the model. (c) 1996 John Wiley & Sons, Inc.     

A porous elastic model for bacterial biofilms: application to the simulation of deformation of bacterial biofilms under microfluidic jet impingement

The presence of bacterial biofilms is detrimental in a wide range of healthcare situations especially wound healing. Physical debridement of biofilms is a method widely used to remove them. This study evaluates the use of microfluidic jet impingement to debride biofilms. In this case, a biofilm is treated as a saturated porous medium also having linear elastic properties. A numerical modeling approach is used to calculate the von Mises stress distribution within a porous medium under fluid-structure interaction (FSI) loading to determine the initial rupture of the biofilm structure. The segregated model first simulates the flow field to obtain the FSI interface loading along the fluid-solid interface and body force loading within the porous medium. A stress-strain model is consequently used to calculate the von Mises stress distribution to obtain the biofilm deformation. Under a vertical jet, 60% of the deformation of the porous medium can be accounted for by treating the medium as if it was an impermeable solid. However, the maximum deformation in the porous medium corresponds to the point of maximum shear stress which is a different position in the porous medium than that of the maximum normal stress in an impermeable solid. The study shows that a jet nozzle of 500 μm internal diameter (ID) with flow of Reynolds number (Re) of 200 can remove the majority of biofilm species.     

Modeling water flow through arterial tissue

A simple model of, water flow through deformable porous media has been developed with emphasis on application to arterial walls. The model incorporates a strain-dependent permeability function into Darcy's Law which is coupled, to the force balance for the bulk material. A simple analytical expression relating water flux (volume flux) to pressure differential is developed which shows how strain-dependent permeability can lead to a reduction in hydraulic conductivity with increasing differential pressure as observed in experiments with arteries. The variation of permeability with position in the wall, which may influence the convective diffusion of macromolecules, is determined for both cylindrical and planar segments and a marked influence of geometry is noted.     

Numerical simulation of blood and interstitial flow through a solid tumor

A theoretical framework is presented for describing blood flow through the irregular vasculature of a solid tumor. The tumor capillary bed is modeled as a capillary tree of bifurcating segments whose geometrical construction involves deterministic and random parameters. Blood flow along the individual capillaries accounts for plasma leakage through the capillary walls due to the transmural pressure according to Sterling’s law. The extravasation flow into the interstitium is described by Darcy’s law for a biological porous medium. The pressure field developing in the interstitium is computed by solving Laplace’s equation subject to derived boundary conditions at the capillary vessel walls. Given the arterial, venous, and tumor surface pressures, the problem is formulated as a coupled system of integral and differential equations arising from the interstitium and capillary flow transport equations. Numerical discretization yields a system of linear algebraic equations for the interstitial and capillary segment pressures whose solution is found by iterative methods. Results of numerical computations document the effect of the interstitial hydraulic and vascular permeability on the fractional plasma leakage. Given the material properties, the fractional leakage reaches a maximum at a particular grade of the bifurcating vascular tree.