Discrete-continuum continuum multiscale model for transport, biofilm development and solid restructuring in porous media

Upscaling transport in porous media including both biomass development and simultaneous structural changes in the solid matrix is extremely challenging. This is because both affect the medium’s porosity as well as mass transport parameters and flow paths. We address this challenge by means of a multiscale model. At the pore scale, the local discontinuous Galerkin (LDG) method is used to solve differential equations describing particularly the bacteria’s and the nutrient’s development. Likewise, a sticky agent tightening together solid or bio cells is considered. This is combined with a cellular automaton method (CAM) capturing structural changes of the underlying computational domain stemming from biomass development and solid restructuring. Findings from standard homogenization theory are applied to determine the medium’s characteristic timeand space-dependent properties. Investigating these results enhances our understanding of the strong interplay between a medium’s functional properties and its geometric structure. Finally, integrating such properties as model parameters into models defined on a larger scale enables reflecting the impact of pore scale processes on the larger scale.

[1]  A. Bensoussan,et al.  Asymptotic analysis for periodic structures , 1979 .

[2]  G. Allaire Homogenization and two-scale convergence , 1992 .

[3]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[4]  Michel Quintard,et al.  Biofilms in porous media: Development of macroscopic transport equations via volume averaging with closure for local mass equilibrium conditions , 2009 .

[5]  Michel Quintard,et al.  Upscaling of transport processes in porous media with biofilms in non-equilibrium conditions , 2010 .

[6]  Tianyu Zhang,et al.  Review of mathematical models for biofilms , 2010 .

[7]  Rainer Helmig,et al.  An upscaled model for biofilm growth in a thin strip , 2010 .

[8]  L. Deacon,et al.  Microbial diversity affects self-organization of the soil–microbe system with consequences for function , 2012, Journal of The Royal Society Interface.

[9]  A. Valocchi,et al.  An improved cellular automaton method to model multispecies biofilms. , 2013, Water research.

[10]  A. Valocchi,et al.  An improved pore‐scale biofilm model and comparison with a microfluidic flow cell experiment , 2013 .

[11]  Albert J. Valocchi,et al.  A hybrid pore‐scale and continuum‐scale model for solute diffusion, reaction, and biofilm development in porous media , 2015 .

[12]  Balthasar Reuter,et al.  FESTUNG: A MATLAB/GNU Octave toolbox for the discontinuous Galerkin method, Part I: Diffusion operator , 2014, Comput. Math. Appl..

[13]  D. Wildenschild,et al.  Biofilm growth in porous media: Experiments, computational modeling at the porescale, and upscaling , 2015, 1510.01218.

[14]  Ali Ebrahimi,et al.  Microbial community dynamics in soil aggregates shape biogeochemical gas fluxes from soil profiles – upscaling an aggregate biophysical model , 2016, Global change biology.

[15]  N. Vasilyeva,et al.  Nonlinear Dynamical Model of Microorganism Growth in Soil , 2016 .

[16]  Peter Knabner,et al.  Derivation and analysis of an effective model for biofilm growth in evolving porous media , 2017 .