Search

The scientific community has recognized that almost 99% of the microbial life on earth is represented by biofilms. Considering the impacts of their sessile lifestyle on both natural and human activities, extensive experimental activity has been carried out to understand how biofilms grow and interact with the environment. Many mathematical models have also been developed to simulate and elucidate the main processes characterizing the biofilm growth. Two main mathematical approaches for biomass representation can be distinguished: continuum and discrete. This review is aimed at exploring the main characteristics of each approach. Continuum models can simulate the biofilm processes in a quantitative and deterministic way. However, they require a multidimensional formulation to take into account the biofilm spatial heterogeneity, which makes the models quite complicated, requiring significant computational effort. Discrete models are more recent and can represent the typical multidimensional structural heterogeneity of biofilm reflecting the experimental expectations, but they generate computational results including elements of randomness and introduce stochastic effects into the solutions. [ABSTRACT FROM AUTHOR]

A multiscale mathematical model is presented to describe de novo granulation, and the evolution of multispecies granular biofilms, in a continuously fed bioreactor. The granule is modelled as a spherical free boundary domain with radial symmetry. The equation governing the free boundary is derived from global mass balance considerations and takes into account the growth of sessile biomass as well as exchange fluxes with the bulk liquid. Starting from a vanishing initial value, the expansion of the free boundary is initiated by the attachment process, which depends on the microbial species concentrations within the bulk liquid and their specific attachment velocity. Nonlinear hyperbolic PDEs model the growth of the sessile microbial species, while quasi-linear parabolic PDEs govern the dynamics of substrates and invading species within the granular biofilm. Nonlinear ODEs govern the evolution of soluble substrates and planktonic biomass within the bulk liquid. The model is applied to an anaerobic, granular-based bioreactor system, and solved numerically to test its qualitative behaviour and explore the main aspects of de novo anaerobic granulation: ecology, biomass distribution, relative abundance, dimensional evolution of the granules and soluble substrates, and planktonic biomass dynamics within the bioreactor. The numerical results confirm that the model accurately describes the ecology and the concentrically layered structure of anaerobic granules observed experimentally, and that it can predict the effects on the process of significant factors, such as influent wastewater composition; granulation properties of planktonic biomass; biomass density; detachment intensity; and number of granules. [ABSTRACT FROM AUTHOR]

The dynamics of biofilm lifecycle are deeply influenced by the surrounding environment and the interactions between sessile and planktonic phenotypes. Bacterial biofilms typically develop in three distinct stages: attachment of cells to a surface, growth of cells into colonies, and detachment of cells from the colony into the surrounding medium. The attachment of planktonic cells from the surrounding environment plays a prominent role in the initial phase of biofilm lifecycle as it initiates the colony formation. During the maturation stage, biofilms harbor numerous microenvironments which lead to metabolic heterogeneity. Such microniches provide conditions suitable for the growth of new species, which are present in the bulk liquid as planktonic cells and can penetrate the porous biofilm matrix. We present a 1D continuum model on the interaction of sessile and planktonic phenotypes in biofilm lifestyle. Such a model is able to reproduce the key role of planktonic cells in the formation and development of biofilms by considering the initial attachment and colonization phenomena. The model is formulated as a hyperbolic–elliptic free boundary value problem with vanishing initial value which considers the concentrations of planktonic and sessile cells as state variables. Hyperbolic equations reproduce the transport and growth of sessile species, while elliptic equations model the diffusion and conversion of planktonic cells and dissolved substrates. The attachment is modeled as a continuous, deterministic process which depends on the concentrations of the attaching species. The growth of new species is modeled through a reaction term in the hyperbolic equations which depends on the concentration of planktonic species within the biofilm. Existence and uniqueness of solutions are discussed and proved for the attachment regime. Finally, some numerical examples show that the proposed model correctly reproduces the growth of new species within the biofilm and overcomes the ecological restrictions characterizing the Wanner–Gujer-type models. [ABSTRACT FROM AUTHOR]