Controlling the hydraulic resistance of membrane biofilms by engineering biofilm physical structure

The application of membrane technology for water treatment and reuse is hampered by the development of a microbial biofilm. Biofilm growth in micro-and ultrafiltration (MF/UF) membrane modules, on both the membrane surface and feed spacer, can form a secondary membrane and exert resistance to permeation and crossflow, increasing energy demand and decreasing permeate quantity and quality. In recent years, exhaustive efforts were made to understand the chemical, structural and hydraulic characteristics of membrane biofilms. In this review, we critically assess which specific structural features of membrane biofilms exert resistance to forced water passage in MF/UF membranes systems applied to water and wastewater treatment, and how biofilm physical structure can be engineered by process operation to impose less hydraulic resistance ("below-the-pain threshold"). Counter-intuitively, biofilms with greater thickness do not always cause a higher hydraulic resistance than thinner biofilms. Dense biofilms, however, had consistently higher hydraulic resistances compared to less dense biofilms. The mechanism by which density exerts hydraulic resistance is reported in the literature to be dependant on the biofilms' internal packing structure and EPS chemical composition (e.g., porosity, polymer concentration). Current reports of internal porosity in membrane biofilms are not supported by adequate experimental evidence or by a reliable methodology, limiting a unified understanding of biofilm internal structure. Identifying the dependency of hydraulic resistance on biofilm density invites efforts to control the hydraulic resistance of membrane biofilms by engineering internal biofilm structure. Regulation of biofilm internal structure is possible by alteration of key determinants such as feed water nutrient composition/concentration, hydraulic shear stress and resistance and can engineer biofilm structural development to decrease density and therein hydraulic resistance. Future efforts should seek to determine the extent to which the concept of "biofilm engineering" can be extended to other biofilm parameters such as mechanical stability and the implication for biofilm control/removal in engineered water systems (e.g., pipelines and/or, cooling towers) susceptible to biofouling.

Lattice Boltzmann Method in Modeling Biofilm Formation, Growth and Detachment

Biofilms are a complex and heterogeneous aggregation of multiple populations of microorganisms linked together by their excretion of extracellular polymer substances (EPS). Biofilms can cause many serious problems, such as chronic infections, food contamination and equipment corrosion, although they can be useful for constructive purposes, such as in wastewater treatment, heavy metal removal from hazardous waste sites, biofuel production, power generation through microbial fuel cells and microbially enhanced oil recovery; however, biofilm formation and growth are complex due to interactions among physicochemical and biological processes under operational and environmental conditions. Advanced numerical modeling techniques using the lattice Boltzmann method (LBM) are enabling the prediction of biofilm formation and growth and microbial community structures. This study is the first attempt to perform a general review on major contributions to LBM-based biofilm models, ranging from pioneering efforts to more recent progress. We present our understanding of the modeling of biofilm formation, growth and detachment using LBM-based models and present the fundamental aspects of various LBM-based biofilm models. We describe how the LBM couples with cellular automata (CA) and individual-based model (IbM) approaches and discuss their applications in assessing the spatiotemporal distribution of biofilms and their associated parameters and evaluating bioconversion efficiency. Finally, we discuss the main features and drawbacks of LBM-based biofilm models from ecological and biotechnological perspectives and identify current knowledge gaps and future research priorities.

Non-Invasive Measurement, Mathematical Simulation and In Situ Detection of Biofilm Evolution in Porous Media: A Review

The development of biofilms and the related changes in porous media in the subsurface cannot be directly observed and evaluated. The primary reason that the mechanism of biofilm clogging in porous media cannot be clearly demonstrated is due to the opacity and structural complexity of three-dimensional pore space. Interest in exploring methods to overcome this limitation has been increasing. In the first part of this review, we introduce the underlying characteristics of biofilm in porous media. Then, we summarize two approaches, non-invasive measurement methods and mathematical simulation strategies, for studying fluid–biofilm–porous medium interaction with spatiotemporal resolution. We also discuss the advantages and limitations of these approaches. Lastly, we provide a perspective on opportunities for in situ monitoring at the field site.

Modeling the impact of evolving biofilms on flow in porous media inside a microfluidic channel

This study integrates microfluidic experiments and mathematical modeling to study the impacts of biofilms on flow in porous media and to explore approaches to simplify modeling permeability with complicated biofilm geometries. E. coli biofilms were grown in a microfluidic channel packed with a single layer of glass beads to reach three biofilm levels: low, intermediate, and high, with biofilm ratios (β_r) of 2.7%, 17.6%, and 55.2%, respectively. Two-dimensional biofilm structures and distributions in the porous medium were modeled by digitizing confocal images and considering broad ranges of biofilm permeability (k_b) (from 10^-15 m² to 10^-7 m²) and biofilm porosity (ε_b) (from 0.2 to 0.8). The overall permeability of the porous medium (k), the flow pathways and the overall/local pressure gradients were found to be highly dependent on β_r and k_b but were moderately impacted by ε_b when the biofilm levels were high and intermediate with k_b>10^-11 m². When biofilm structures are well developed, simplified biofilm geometries, such as uniform coating and symmetric contact filling, can provide reasonable approximations of k.

Pore-Scale Hydrodynamics in a Progressively Bioclogged Three-Dimensional Porous Medium: 3-D Particle Tracking Experiments and Stochastic Transport Modeling

Biofilms are ubiquitous bacterial communities that grow in various porous media including soils, trickling, and sand filters. In these environments, they play a central role in services ranging from degradation of pollutants to water purification. Biofilms dynamically change the pore structure of the medium through selective clogging of pores, a process known as bioclogging. This affects how solutes are transported and spread through the porous matrix, but the temporal changes to transport behavior during bioclogging are not well understood. To address this uncertainty, we experimentally study the hydrodynamic changes of a transparent 3-D porous medium as it experiences progressive bioclogging. Statistical analyses of the system's hydrodynamics at four time points of bioclogging (0, 24, 36, and 48 h in the exponential growth phase) reveal exponential increases in both average and variance of the flow velocity, as well as its correlation length. Measurements for spreading, as mean-squared displacements, are found to be non-Fickian and more intensely superdiffusive with progressive bioclogging, indicating the formation of preferential flow pathways and stagnation zones. A gamma distribution describes well the Lagrangian velocity distributions and provides parameters that quantify changes to the flow, which evolves from a parallel pore arrangement under unclogged conditions, toward a more serial arrangement with increasing clogging. Exponentially evolving hydrodynamic metrics agree with an exponential bacterial growth phase and are used to parameterize a correlated continuous time random walk model with a stochastic velocity relaxation. The model accurately reproduces transport observations and can be used to resolve transport behavior at intermediate time points within the exponential growth phase considered.

Biofilm imaging in porous media by laboratory X-Ray tomography: Combining a non-destructive contrast agent with propagation-based phase-contrast imaging tools

X-ray tomography is a powerful tool giving access to the morphology of biofilms, in 3D porous media, at the mesoscale. Due to the high water content of biofilms, the attenuation coefficient of biofilms and water are very close, hindering the distinction between biofilms and water without the use of contrast agents. Until now, the use of contrast agents such as barium sulfate, silver-coated micro-particles or 1-chloronaphtalene added to the liquid phase allowed imaging the biofilm 3D morphology. However, these contrast agents are not passive and potentially interact with the biofilm when injected into the sample. Here, we use a natural inorganic compound, namely iron sulfate, as a contrast agent progressively bounded in dilute or colloidal form into the EPS matrix during biofilm growth. By combining a very long source-to-detector distance on a X-ray laboratory source with a Lorentzian filter implemented prior to tomographic reconstruction, we substantially increase the contrast between the biofilm and the surrounding liquid, which allows revealing the 3D biofilm morphology. A comparison of this new method with the method proposed by Davit et al (Davit et al., 2011), which uses barium sulfate as a contrast agent to mark the liquid phase was performed. Quantitative evaluations between the methods revealed substantial differences for the volumetric fractions obtained from both methods. Namely, contrast agent—biofilm interactions (e.g. biofilm detachment) occurring during barium sulfate injection caused a reduction of the biofilm volumetric fraction of more than 50% and displacement of biofilm patches elsewhere in the column. Two key advantages of the newly proposed method are that passive addition of iron sulfate maintains the integrity of the biofilm prior to imaging, and that the biofilm itself is marked by the contrast agent, rather than the liquid phase as in other available methods. The iron sulfate method presented can be applied to understand biofilm development and bioclogging mechanisms in porous materials and the obtained biofilm morphology could be an ideal basis for 3D numerical calculations of hydrodynamic conditions to investigate biofilm-flow coupling.

Microbial competition in porous environments can select against rapid biofilm growth

Microbes often live in dense communities called biofilms, where competition between strains and species is fundamental to both evolution and community function. Although biofilms are commonly found in soil-like porous environments, the study of microbial interactions has largely focused on biofilms growing on flat, planar surfaces. Here, we use microfluidic experiments, mechanistic models, and game theory to study how porous media hydrodynamics can mediate competition between bacterial genotypes. Our experiments reveal a fundamental challenge faced by microbial strains that live in porous environments: cells that rapidly form biofilms tend to block their access to fluid flow and redirect resources to competitors. To understand how these dynamics influence the evolution of bacterial growth rates, we couple a model of flow-biofilm interaction with a game theory analysis. This investigation revealed that hydrodynamic interactions between competing genotypes give rise to an evolutionarily stable growth rate that stands in stark contrast with that observed in typical laboratory experiments: cells within a biofilm can outcompete other genotypes by growing more slowly. Our work reveals that hydrodynamics can profoundly affect how bacteria compete and evolve in porous environments, the habitat where most bacteria live.

Review of pore network modelling of porous media: Experimental characterisations, network constructions and applications to reactive transport

Pore network models have been applied widely for simulating a variety of different physical and chemical processes, including phase exchange, non-Newtonian displacement, non-Darcy flow, reactive transport and thermodynamically consistent oil layers. The realism of such modelling, i.e. the credibility of their predictions, depends to a large extent on the quality of the correspondence between the pore space of a given medium and the pore network constructed as its representation. The main experimental techniques for pore space characterisation, including direct imaging, mercury intrusion porosimetry and gas adsorption, are firstly summarised. A review of the main pore network construction techniques is then presented. Particular focus is given on how such constructions are adapted to the data from experimentally characterised pore systems. Current applications of pore network models are considered, with special emphasis on the effects of adsorption, dissolution and precipitation, as well as biomass growth, on transport coefficients. Pore network models are found to be a valuable tool for understanding and predicting meso-scale phenomena, linking single pore processes, where other techniques are more accurate, and the homogenised continuum porous media, used by engineering community.

Dissolved Organic Carbon Mobilisation in a Groundwater System Stressed by Pumping

The concentration and flux of organic carbon in aquifers is influenced by recharge and abstraction, and surface and subsurface processing. In this study groundwater was abstracted from a shallow fractured rock aquifer and dissolved organic carbon (DOC) was measured in observation bores at different distances from the abstraction bore. Groundwater abstraction at rates exceeding the aquifers yield resulted in increased DOC concentration up to 3,500 percent of initial concentrations. Potential sources of this increased DOC were determined using optical fluorescence and absorbance analysis. Groundwater fluorescent dissolved organic material (FDOM) were found to be a combination of terrestrial-derived humic material and microbial or protein sourced material. Relative molecular weight of FDOM within four metres of the abstraction well increased during the experiment, while the relative molecular weight of FDOM between four and ten metres from the abstraction well decreased. When the aquifer is not being pumped, DOC mobilisation in the aquifer is low. We hypothesise that the physical shear stress on aquifer materials caused by intense abstraction significantly increases the temporary release of DOC from sloughing of biofilms and release of otherwise bound colloidal and sedimentary organic carbon (SOC).

A mathematical model and computational framework for three-dimensional chondrocyte cell growth in a porous tissue scaffold placed inside a bi-directional flow perfusion bioreactor

The in vitro chondrocyte cell culture for cartilage tissue regeneration in a perfusion bioreactor is a complex process. Mathematical modeling and computational simulation can provide important insights into the culture process, which would be helpful for selecting culture conditions to improve the quality of the developed tissue constructs. However, simulation of the cell culture process is a challenging task due to the complicated interaction between the cells and local fluid flow and nutrient transport inside the complex porous scaffolds. In this study, a mathematical model and computational framework has been developed to simulate the three-dimensional (3D) cell growth in a porous scaffold placed inside a bi-directional flow perfusion bioreactor. The model was developed by taking into account the two-way coupling between the cell growth and local flow field and associated glucose concentration, and then used to perform a resolved-scale simulation based on the lattice Boltzmann method (LBM). The simulation predicts the local shear stress, glucose concentration, and 3D cell growth inside the porous scaffold for a period of 30 days of cell culture. The predicted cell growth rate was in good overall agreement with the experimental results available in the literature. This study demonstrates that the bi-directional flow perfusion culture system can enhance the homogeneity of the cell growth inside the scaffold. The model and computational framework developed is capable of providing significant insight into the culture process, thus providing a powerful tool for the design and optimization of the cell culture process.

Effect of permeable biofilm on micro- and macro-scale flow and transport in bioclogged pores

Simulations of coupled flow around and inside biofilms in pores were conducted to study the effect of porous biofilm on micro- and macro-scale flow and transport. The simulations solved the Navier-Stokes equations coupled with the Brinkman equation representing flow in the pore space and biofilm, respectively, and the advection-diffusion equation. Biofilm structure and distribution were obtained from confocal microscope images. The bulk permeability (k) of bioclogged porous media depends on biofilm permeability (kbr) following a sigmoidal curve on a log-log scale. The upper and lower limits of the curve are the k of biofilm-free media and of bioclogged media with impermeable biofilms, respectively. On the basis of this, a model is developed that predicts k based solely on kbr and biofilm volume ratio. The simulations show that kbr has a significant impact on the shear stress distribution, and thus potentially affects biofilm erosion and detachment. The sensitivity of flow fields to kbr directly translated to effects on the transport fields by affecting the relative distribution of where advection and diffusion dominated. Both kbr and biofilm volume ratio affect the shape of breakthrough curves.

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.

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.