An invariant-region-preserving scheme for a convection-reaction-Cahn-Hilliard multiphase model of biofilm growth in slow sand filters

A multidimensional model of biofilm growth present in the supernatant water of a Slow Sand Filter is derived. The multiphase model, consisting of solid and liquid phases, is written as a convection-reaction system with a Cahn-Hilliard-type equation with degenerate mobility coupled to a Stokes-flow equation for the mixture velocity. An upwind discontinuous Galerkin approach is used to approximate the convection-reaction equations, whereas an $H^1$-conforming primal formulation is proposed for the Stokes system. By means of a splitting procedure due to the reaction terms, an invariant-region principle is shown for the concentration unknowns, namely non-negativity for all phases and an upper bound for the total concentration of the solid phases. Numerical examples with reduced biofilm reactions are presented to illustrate the performance of the model and numerical scheme.