Mathematical for model microbial fuel cells with anodic biofilms
and anaerobic digestion
Cristian Picioreanu*^{1}, Krishna P.
Katuri^{2}, Ian M. Head^{2},
Mark C.M. van Loosdrecht^{1}, Keith Scott^{2}
^{1} Department of Biotechnology, Delft University of
Technology
Julianalaan 67, 2628 BC Delft, The Netherlands.
Tel: +31(0)152781551, Fax: +31(0)152782355
^{2} University of Newcastle upon Tyne, Newcastle, UK
*Corresponding author. Email: 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):965971.
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
AE
4.2  3d simulations in Figure 3
FG
Supplementary Material
1.
Spatial model scales
Figure S1
Compartments and subdomains 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 subdomains 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:
nonelectroactive reactant, electroactive reactant, nonelectroactive
product, and electroactive product.
2.
Model solution algorithms
Figure
S2. Model solution
algorithm and software implementation. The threedimensional (3d) model
construction allows a general description of the timedependent 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 twodimensional and onedimensional
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)
Particlebased multidimensional multispecies biofilm model. Appl Environ Microb 70(5), 30243040.
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),
10851103.
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.
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 AE)

Download
files
MFC2d_AD11_Fig3.avi 6.4 Mb, Microsoft RLE MFC2d_AD11_Fig3.gif 2.5 Mb, animated gif MFC2d_AD11_Fig3.mov 9.0 Mb, Quicktime (higher quality  recommended) 
Download
files
