TU Delft

Mathematical for model microbial fuel cells with anodic biofilms and anaerobic digestion

Cristian Picioreanu*1, Krishna P. Katuri2, Ian M. Head2, Mark C.M. van Loosdrecht1, Keith Scott2

1 Department of Biotechnology, Delft University of Technology
Julianalaan 67, 2628 BC Delft, The Netherlands.
Tel: +31-(0)15-2781551, Fax: +31-(0)15-2782355

2 University of Newcastle upon Tyne, Newcastle, UK

*Corresponding author. E-mail: c.picioreanu@tudelft.nl

This web page contains supplementary material for the paper "Mathematical model for microbial fuel cells with anodic biofilms and anaerobic digestion" submitted to the Anaerobic Digestion 11 conference in Brisbane, September 2007.
This paper appeared in Water Science and Technology, 57(7):965-971.

Go to:
1. Description of spatial model scales
2. Model solution algorithms
3. List of parameters used in the simulations
    3.1 - 2d/3d simulations in Figure 3
    3.2 - 2d simulations in Figure 5
4. Animation of simulation results
    4.1 - 2d simulations in Figure 3 A-E
    4.2 - 3d simulations in Figure 3 F-G

Supplementary Material

1. Spatial model scales

Figure S1    Compartments and sub-domains of the computational model. (A) The large scale comprises an ideally mixed bulk liquid and a biofilm compartment attached to the electrode surface. (B) The small scale biofilm model contains three sub-domains characterized by the presence of different processes: (1) a completely mixed zone B which is connected to the bulk liquid, (2) a mass transfer boundary layer L and, (3) the biofilm matrix F developing on the electrode support. Typical solute concentration profiles averaged along the x direction are represented as a function of distance z from the electrode surface for: non-electroactive reactant, electroactive reactant, non-electroactive product, and electroactive product.


2. Model solution algorithms

Figure S2.  Model solution algorithm and software implementation. The three-dimensional (3d) model construction allows a general description of the time-dependent spatial distributions of dissolved and biomass components, and of currents on the anode surface. The model was implemented in this full 3d setup, but also in computationally faster two-dimensional and one-dimensional reductions. The software implementation is written on a modular base in C/C++, and it runs on ordinarily available personal computers. The program code is based on the multidimensional biofilm modelling software described in Picioreanu et al. (2004) and Xavier et al. (2005), with major additions for the specific MFC processes. Being a rather flexible framework in which specific model setups can be built, the software allows an easy choice of any number of components and processes. Therefore, it is relatively easy to construct different scenarios and easily test various hypotheses, within the general model assumptions described.

Picioreanu, C., Kreft, J.-U. and van Loosdrecht, M.C.M. (2004) Particle-based multidimensional multispecies biofilm model. Appl Environ Microb 70(5), 3024-3040.

Xavier, J.B., Picioreanu, C. and van Loosdrecht, M.C.M. (2005) A framework for multidimensional modelling of activity and structure of multispecies biofilms. Environ Microbiol 7(8), 1085-1103.

Press, W.H., Teukolsky, S.A., Vetterling, W.T. and Flannery, B.P. (1997) Numerical recipes in C: The art of scientific computing. Cambridge University Press, NY.


3. Parameters used in the simulations

3.1    Figure 3.      Parameter file used as model input for the 2d simulations in Figure 3.

3.2    Figure 5.     

Parameter file used as model input for the 2d simulations in Figure 5A (Low R, 100 Ohm).

Parameter file used as model input for the 2d simulations in Figure 5B (High R, 1000 Ohm).

4. Animations of simulation results

4.1  Case #1,   2d
(Figure 3 A-E)


Download files

6.4 Mb,  Microsoft RLE

2.5 Mb,  animated gif

9.0 Mb,  Quicktime
(higher quality - recommended)
4.2  Case #1,   3d
(Figure 3 F-G)

Download files

3.6 Mb,  Microsoft AVI

2007 - Updated on April 10, 2007, by Cristian Picioreanu