Numerical schemes for mixture theory models with filling constraint: application to biofilm ecosystems

¹ BIOCORE, Inria, Sorbonne Université, Nice Sophia-Antipolis
² IRMA, Université de Strasbourg, CNRS UMR-7501, Strasbourg, France
³ Université Côte d’Azur, Inria, CNRS, LJAD
⁴ IDP, Université d’Orléans, CNRS, CNRS UMR-7013, Orléans, France
⁵ Université Côte d’Azur, Inria, CNRS, LJAD
⁶ IMB, Université de Bordeaux, UMR CNRS 5251, Bordeaux
⁷ LmB, Université de Franche-Comté, CNRS UMR-6623, Besançon, France

Abstract

The mixture theory framework is a powerful way to describe multi-phasic systems at an intermediary scale between microscopic and macroscopic scales. In particular, mixture theory reveals a powerful approach to represent microbial biofilms where a consortium of cells is embedded in a polymeric structure. To simulate a model of microalgal biofilm, we propose an upgraded numerical scheme, consolidating the one proposed by Berthelin et al. (2016) to enforce the volume-filling constraint in mixture models including mass exchanges. The strategy consists in deducing the discrete version of the incompressibility constraint from the discretized mass balance equations. Numerical simulations show that this method constrains the total volume filling constraint, even at the discrete level. Moreover, we add viscous terms in the biofilm model to properly represent biofilms interactions with its fluidic environment. It turns out that a well-balanced numerical scheme becomes of outmost importance to capture the biofilm dynamic when including the viscosity. This modelling upgrade also involves recalibrating model parameters. In particular, the elastic tensors to recover realistic front features. With the new parameters, the numerical set-up becomes more demanding to reach convergence.