Received 3 February 2016; accepted 25 April 2016; published 28 April 2016
Microbial fuel cells (MFC) can be defined as a microbial catalyzed electrochemical system which can facilitate the direct conversion of substrate to electricity through a cascade of redox reactions, especially in the absence of oxygen  . The applications of MFC are widespread in different fields including waste water remediation, toxic pollutants/xenobiotics removal, recovery of commercially viable products, sequestration of CO2 harvesting the energy stored in marine sediments, desalination, etc.  . Microbial electrochemical cells are recognized as a modern technology to directly utilize bioenergy stored in organic substances, especially in wastewater  . Several experiments have been conducted to evaluate MXCs main performance as a current or hydrogen generator fed with different organic matters  and  . A simple mediator based model with suspended cells was investigated  . A simple model with rapid implementation and computations is used to describe the effect of some operational conditions such as temperature and substrate concentration on in the MFC performance  . Pinto developed a time-dependent mathematical model with the uniform distribution of bacteria in the anode chamber. Although a number of MFC mathematical models have been developed and discussed. To the best of the knowledge, only one MEC model has been proposed  . Yahya modified this model for a fed-batch reactor. It is a multi population mediator-based model developed based on the Bernard’s anaerobic digestion kinetics model  . Alavijeh  used a variety of approaches to develop the first generalized conduction-based model for MXCs including both MFCs and MECs. It is a one-dimensional spatial distribution and time-dependent model using Bernard’s anaerobic digestion kinetics model and both biofilm and liquid bulk simulation. The purpose of this communication is to derive the analytical expression for acetate and methane concentration using the Adomain decomposition method. We also provide the tabular complication of concentration of acetate with limiting case results (first order and zero order kinetics).
2. Mathematical Formulation of the Problem
Anolyte contains fermentative microorganisms and acetoclastic methanogens. Biofilm contains acetoclastic methanogens and anode respiring bacteria (electrogens). Acetate is produced during fermentation process and then diffuses to the biofilm where electrogens consume it and conduct electrons to the anode surface  . The systematic diagram of the model is represented in Figure 1.
The acetate and methane mass transfer equations through the biofilm are described as follows  :
Figure 1. Schematic representation of the model  .
where and the diffusion coefficient of acetate and methane in the bioflim. and are the concentration of acetate on the biofilm and methane concentration on the biofilm., and are the density of biomass yield coefficient. The acetate consumption rate by electromagnetic microorganism in the biofuel is represented by the Nenst-Monoid equation 
where is the maximum uptake, are the volume fraction, and Half saturated constant. are the Faraday constant, universal gas constant, temperature and local electrical potential of the biofilm respectively. The acetate consumption rate by acetoclastic methangen microorganism in the bioflim is
where and are the volume fraction and Half saturated constant. At the anode surface, there is no substrate flux and at the surface of the biofilm, there is an interface transfer. The boundary conditions for the above equations are given by 
where represents the space coordinate in the biofilm (m), thickness of the biofilm (m) and boundary layer thickness (m) respectively. and are the diffusion coefficient of acetate and methane in the liquid. is the acetate concentration in the liquid bulk and in the biofilm interface
and is the methane concentration in the liquid bulk and in the biofilm interface. We introduce the following set of dimensionless variables:
Using the above dimensionless variables the non-linear reaction-diffusion Equations ((1) and (2)) are expressed in the following dimensionless form:
The boundary conditions can be written as follows:
3. Approximate Analytical Expression of Concentration of Acetate and Methane
Recently, many authors have been applied the Adomain decomposition method (ADM) to various problems and demonstrated the efficiency of the ADM for handling non-linear problem in physics and engineering sciences  -  . The modified Adomain decomposition method  is used to give the approximate solutions of the non-linear Equations ((7) and (8)). Many researchers find that the ADM requires less computational work than traditional approaches  -  . Other advantages include the ability to solve nonlinear problems without linearization, the wide applicability to several types of problems and scientific fields, and the development of a reliable, analytic solution. Many researchers find that the ADM requires less computational work than traditional approaches  -  . Other advantages include the ability to solve nonlinear problems without linearization, the wide applicability to several types of problems and scientific fields, and the development of a reliable, analytic solution. Using this method (refer Appendix A), we can obtain the concentrations acetate and methane as follows:
where the constants
Equations ((13) and (14)) are valid provided and. This is the only limitations in this method.
4. Limiting Case
4.1. Unsaturated (First Order) Catalysis
We initially consider the situation where the concentration of acetate and methane is less than the half saturation constants and. In this case Hence, Equations ((7) and (8)) reduces to
Hence, the non-linear Equations ((7) and (8)) have been reduces to linear equations. Now, the concentration of acetate and methane for corresponding boundary conditions (9a) to (9c) becomes as follows:
where and. Equations ((18) and (19)) are the exact solution of Equations ((16) and (17)).
4.2. Saturated (Zero-Order) Catalysis
We now consider that the second major limiting situation found in practice, when the concentration of acetate and methane is very much greater than the half saturation constants and. In this case,. Hence, the non-linear Equations ((7) and (8)) have been reduces to
The above Equations ((17) and (18)) are linear reaction-diffusion equations which are exactly solvable. By solving the above Equations ((17) and (18)), we can obtain the concentration of Acetate (16), and Methane (17).
Equations ((22) and (23)) are the exact solution of Equations ((20) and (21)).
4.3. Saturated Electrogenic Microorganism and Acetoclastic Methanogens Are Equal (α = β)
In this case, Equations ((7) and (8)) become as follows:
In this case, the above non-linear equation can be solved using Adomain decomposition method. Now, the concentrations become
where the constants to and k are given in Equation (12), when replacing. Equations ((26) and (27)) are the approximate analytical expression of concentration of acetate and methane.
Equations ((7) and (8)) represent the general closed-form of analytical expression for the concentrations of acetate and methane for non steady state condition and for various system parameters (potential, saturation parameter of electrogenic microorganism and acetoelastic methanogenes, the diffusion coefficient of acetate, ratio of the thickness of the biofilm and boundary layer). It is of interest to compare the influence of each parameter on the concentration of acetate and methane for various realistic experimental parameters.
Influence of Potential on the Concentration of Acetate. The influence of dimensionless potential on the concentration of the acetate for some experimental values of parameters is shown in Figure 2(a). The microbial activity is strongly dependent on the redox potential of the anode. From this figure it is observed that the concentration of acetate decreases when increases or potential decreases.
Figure 2. Plot of dimensionless concentration versus dimensionless distance x using Equation (10) for various values of the parameter (a), (b), (c) and, (d) and for some fixed experimental values of other parameters.
acetate increases when saturation parameter of electrogenic microrganism or bulk concentration of the acetate increases and saturation parameter of acetoelastic methanogenes decreases. From Figure 2(c) it is also observed that the concentration of acetate is inversely proportional to diffusion coefficient of acetate.
Influnce of the Ratio of Thickness of the Biofilm and the Boundary Layer. Figure 2(d) represents the concentration verses distance from the anode surface for various values of or ratio of biofilm thickness and boundary layer thickness. From this figure, it is inferred that the concentration of acetate increases when the ratio of thickness increases.
Influence of Other Parameters of the Concentration of Methane and Acetate. The concentration of methane versus dimensionless distance x for various experimental values of parameters is plotted in Figure 3. From these figure, it is interfered that the concentration of methane increases when increases or decreases.
Figure 4 represents the dimensionless concentrations of acetate versus potential for various values of and. From these figures, it is observed that the concentration of acetate increases when thickness and decreases. From Figure 4(d), it is observed that the concentration of acetate does not differ significantly about the parameter.
Figure 3. Plot of dimensionless concentration of versus dimensionless distance x using Equation (11). (a) For various values of the parameter (a), (b), (c) and, (d) and for some fixed experimental values of other parameters.
Figure 4. Plot of dimensionless concentration versus dimensionless potential using Equation (10) for various values of the parameter (a), (b), (c), (d) and for some fixed experimental values of other parameters.
6. Comparison with Numerical Data and Limiting Case Results
The non-linear differential Equations ((9) and (10)) for the given initial-boundary conditions are being solved numerically. The function pdex, in Matlab software which is a function of solving the initial-boundary value problems for non-linear ordinary differential equations is used to solve this equation. Its numerical solution is compared with analytical results in Table 1. The maximum relative error between our analytical result and the numerical result is 0.32%. The Matlab program is also given in Appendix C. The concentration of acetate and methane are also obtained for the following limiting cases, that is zero order kinetics, first order kinetics and saturated electrogenic microorganism and saturated acetoclastic methanogens are equal. Also our analytical results are compared with limiting case results in Figure 5 and it gives a satisfactory agreement.
7. Determination of Kinetic Parameters KAC,E, qAC,AM,max, φE,a, α and β
The acetate consumption rate by electromagnetic microorganism in the microbial fuel cell (Equation (3)) can be written as follows:
Table 1. Comparison of the acetate concentration calculated using Equation (10) with the numerical simulation for various experimental values of when.
Figure 5. Plot of two dimensional comparative diagram of general solution of dimensionless Acetate concentration and methane concentration with the limiting cases (Zero order and in (a) and (b) respectively) for the experimental values Comparison of dimensionless concentration of and using Equations ((10), (11)), Equations ((22), (23)), Equations ((26), (27)).
As shown in Figure 6, is plotted against to obtain the straight line with the slope and intercept. The slope and intercept yields the value of the parameters and.
From Equation (10), we can obtain the concentration of acetate at bioflim and anode interface as
Now the plot of versus gives the slope and intercept. From these plot, we can obtain the kinetic constant and (Table 2).
Figure 6. Estimation of kinetic parameters and from Equation (28).
Table 2. Nomenclature.
A theoretical model describing the bio energy production using microbial electrochemical cell via Nernst-Mo- noid kinetics is analyzed. The time independent non-linear partial differential equations have been solved analytically using the Adomain decomposition method. The primary result of this work is the approximate analytical expression of concentration of acetate and methane for all values of parameters. The influence of potential, ratio of thickness of biofilm and boundary layer, etc. on the concentration of acetate and methane is discussed. Our results are in excellent agreement with stimulation and limiting case results. Also two graphical procedures are suggested for estimating the kinetic parameters.
This work was supported by the DST SB/SI/PC-50/2012, New Delhi, India. The authors are thankful to Mr. S. Mohamed Jaleel, The Chairman, Dr. A. Senthilkumar, The Principal, Dr. P. G. Jansi Rani, Head of the Department of Mathematics, SethuInistitute of Technology, Kariapatti-626115, Tamilnadu, India for their encouragement.
Appendix A: Basic Concept of the Adomain Decomposition Method
This is given in the supplementary material of the manuscript.
Appendix B: Approximate Analytical Solution of Non Linear Equation (7) Using ADM
In this appendix, we indicate how Equation (6) in this paper is derived. Furthermore, an ADM is constructed to determine the solution of Equation (4) in the operator form,
where. Applying the inverse operator on both sides of Equation (B.1) yields
where A and B are the constants of integration. We let,
From Equations ((B.3) to (B.5)), Equation (B.2) becomes
We identify the zeroth component as
and the remaining components as the recurrence relation
where are the Adomain polynomials of.We can find the first two as follows:
Adding (B.9) and (B.10), we can obtain the concentration of acetate as described in Equation (10) in the text. By substituting the values of in Equation (8), we get the concentration of the methane (Equation (11)) in the text.
Appendix C: Scilab/Matlab Program for the Numerical Solution of Equation
This is given in the supplementary material of the manuscript.
Supplementary Material of the Manuscript
Appendix A: Basic Concept of the Modified Adomain Decomposition Method
Consider the singular boundary value problem of order nonlinear differential equation in the form
where N is a non-linear differential operator of order less than n, is given function and are given constants. We propose the new differential operator, as below
where, so, the problem can be written as
The inverse operator is therefore considered a fold integral operator, as below  
By applying on (A.3), we have
The Adomian decomposition method introduce the solution and the nonlinear function by infinite series
where the components of the solution will be determined recurrently. Specific algorithms were seen in   to formulate Adomian polynomials. The following algorithm:
can be used constant Adomian polynomials, when is a nonlinear function. By substituting (A.6) and (A.7) in to (A.5)
Through using modified Adomian decomposition method, the components can be determined as
From (A.8) and (A.11), we can determine the components, and hence the series solution of in (A.6) can be immediately obtained. For numerical purposes, the n-term approximate
can be used to approximate the exact solution. The approach presented above can be validated by testing it on a variety of several linear and nonlinear initial value problems.
Scilab/MatlabProgram for the Numerical Solution of Equation (4)
m = 0;
x = linspace (0, 1);
t = linspace (0, 100000);
sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t);
u1 = sol(:,:,1);
function [c,f,s] = pdex4pde(x,t,u,DuDx)
f =1.* DuDx;
function u0 = pdex4ic(x);
u0 = ;
function [pl,ql,pr,qr] = pdex4bc(xl,ul,xr,ur,t)
pl = ;
ql = ;
pr = [-j*(1-ur(1))];
qr = ;
 Mohan, S.V., Srikanth, S., Velvizhi, G. and Lenin Babu, M. (2013) Microbial Fuel Cells Forsustainable Bioenergy Generati-On: Principles and Perspective Applications. Biofuels Technologies, 14, 335-368.
 Alavijeh, M.K., Mardanpour, M.M. and Yaghmaei. S. (2015) A Generalized Model for Complex Waste Water Treatment with Simultaneous Bioenergy Production Using the Microbial Electrochemical Cell. Electrochimica Acta, 167, 84-96.
 Pant, D., Van Bogaert, G., Diels, L. and Vanbroekhoven, K. (2010) A Review of the Substrates Used in Microbial Fuel Cells (MFCs) for Sustainable Energy Production. Bioresource Technology, 101, 1533-1543. http://dx.doi.org/10.1016/j.biortech.2009.10.017
 Kadier, A., Simayi, Y. Kalil, M.S. Abdeshahian, P. and Hamid, A.A. (2014) Review of the Substrates Used in Microbialele-Ctrolysis Cells (MECs) for Producing Sustainable and Clean Hydrogen Gas. Renewable Energy, 71, 466-472. http://dx.doi.org/10.1016/j.renene.2014.05.052
 Pinto, R.P., Srinivasan, B., Manuel, M.F. and Tartakovsky, B. (2010) A Two-Population Bioelectrochemical Model of a Microbial Fuel Cell. Bioresource Technology, 101, 5256-5265. http://dx.doi.org/10.1016/j.biortech.2010.01.122
 Yahya, A.M., Hussain, M.A. and Abdul Wahab, A.K. (2015) Modeling Optimization, and Control of Microbial Electrolysis Cells in a Fed-Batch Reactor for Production of Renewable Bio Hydrogen Gas. International Journal of Energy Research, 39, 557-572. http://dx.doi.org/10.1002/er.3273
 Saravanakumar, K. and Rajendran, L. (2011) Analytical Solution of the Concentration of Substrate and Effectiveness Factor for Acetophenone in Packed Bed Reactor. International Journal of Mathematical Archive, 2, 2347-2357.
 Yahya, Q.H. and Liu, M.Z. (2009) Solving Singular Boundary Value Problems of Higher Order Ordinary Differential Equation by Modifiedadomain Decomposition Method. Communications in Nonlinear Science and Numerical Simulation, 14, 2592-2596. http://dx.doi.org/10.1016/j.cnsns.2008.09.027
 Yahya, Q.H. and Liu, M.Z. (2008) Modified Adomian Decomposition Method for Singular initial Value Problems in the Second-Order Ordinary Differential Equations. Surveys in Mathematics and Its Applications, 3, 183-193.