Mathematical modelling of a biofilm : The Adomian decomposition method

A mathematical modelling by a biofilm under steady state conditions is discussed. The nonlinear differential Equations in biofilm reaction is solved using the Adomian decomposition method. Approximate analytical expressions for substrate concentration have been derived for all values of parameters δ and SL. These analytical results are compared with the available numerical results and are found to be in good agreement.


INTRODUCTION
Microorganisms biofilms adhere to the interfaces between gas and liquid phases, liquid and solid phases, or two liquid phases [1].Their activity can have an adverse effect, e.g., biofilms damage materials [2], and water purification technology [3].The situation is largely similar to the case of heterogeneous reaction in a porous layer [4,5]; however, for the biofilm kinetics, there are a number of specific features.The dependence of the biochemical reaction rate on a substrate concentration is described by a saturating curve and can not be characterized by a power function [3,6].Dueck et al. [7-9] described the mathematical modeling of evaluation of biofilm with allowance for its erosion.Recently, Min kov et al. [10] obtained the analytical expression for the substrate flux in to the biofilm as the function of volumetric flow rate of the substrate solution.To my knowledge no rigorous analytical solution of concentration of substrate and flux into the biofilm for a square law of microbial death rate with steady-state conditions has been reported.The purpose of this communication is to derive approximate analytical expressions for the steady-state concentration of substrate and flux into the biofilm using the Adomian decomposition method for all values of biofilm thickness and substrate concentration outside the biofilm.

MATHEMATICAL FORMULATION OF THE PROBLEM
It is assumed that the substrate consumption is described by the Michaelis-Menten kinetics.The Equation of which are derived on the basis of the theory of enzymatic reactions [3,5] on the particle surface of biofilm [10] is of the following form: The boundary conditions are The biomass balance [10].
From Eq.4, the concentration of active biomass can be expressed through the substrate concentration.Now the Eq.1 can be written in the form where S f is the substrate concentration in the biofilm, K is the Michaelis-Menten constant, z is the co-ordinate, L f is the biofilm thickness, D f is the diffusion coefficient within the biofilm, b is the Microbial death constant, q is the substrate consumption rate constant, S 1 is the substrate concentration outside the bioflim and Y is the biomass yield per unit amount of substrate consumed respectively.The non-linear ODE (Eq.5) is made dimensionless by defining the following parameters: The above Eq.5 reduces to the following dimensionless form: The dimensionless concentration flux into the biofilm is given by

SOLUTION OF BOUNDARY VALUE PROBLEM BY THE ADOMIAN DECOMPOSITION METHOD
The Adomian's decomposition method has been successfully applied to linear and nonlinear problems.One of its advantages is that it provides a rapid convergent series solution.However, in this method, some modifications are proposed by several authors [11][12][13][14][15].By applying the Adomian's decomposition method, a new iterative method to compute nonlinear Equations are developed.The Adomian decomposition method is an extremely simple method [11][12][13][14][15] to solve the non-linear differential Equations.First iteration is enough.Furthermore, the obtained result is of high accuracy.Using this Adomian decomposition method (see Appendix A and B), the solution of Eq.7 becomes: The solution of concentration flux into the biofilm is obtained as

NUMERICAL SIMULATION
The non-linear differentials Eq.7 is also solved by numerical methods.The function bvp4c in Matlab software which is a function of solving two-point boundary value problems (BVPs) for ordinary differential equations is used to solve this equation.The Matlab program is also given in Appendix C. Its numerical solution is compared with Adomian decomposition method in Tables 1 and 2 and Figures 1-3 for various value of parameters.

RESULTS AND DISCUSSION
An approximate analytical expression of concentrations S is given in the Eq.11.The concentration         log .From this figure, it is inferred that, the value of the flux is high when S L is large and then decreases slowly and reaches   

CONCLUSION
This paper reports a mathematical treatment for analyzing biofilm for a square law of microbial death rate.In this paper, we have evaluated a theoretical model for an investigation of the dynamic behavior of substrate consumption by a biofilm.The approximate analytical expressions for the steady state substrate concentrations for all values of biochemical parameters (δ and S L ) were obtained using Adomian decomposition method.Further- more, an analytical expression corresponding to the steady state flux response is also presented.A satisfactory agreement with the existing results is noted.This theoretical result is useful to further develop the model involving the balance of production of active biomass and bioilm erosion.

ACKNOWLEDGEMENTS
Then substituting Eqs.A.5 and A.6 in Eq.A. 4 gives Then equating the terms in the linear system of Eq.A.8 gives the recurrent relation However, in practice all the terms of series in Eq.A.7 cannot be determined, and the solution is approximated by the truncated series 0 .
n n y    This method has been proven to be very efficient in solving various types of non-linear boundary and initial value problems.

APPENDIX B Analytical Solutions of Concentrations of Substrate Using ADM
In this appendix, we derive the general solution of nonlinear Eq.11 by using Adomian decomposition method.We write the Eq.11 in the operator form,   where A and B are the constants of integration.We let, We identify the zeroth component as and the remaining components as the recurrence relation where A n are the Adomian polynomials of We can find the first few A n as follows: t=linspace(0,1000); 2  (B.10) sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t); u1 = sol(:,:,1); figure plot(x,u1(end,:)) title('u1(x,t


Figures 1-3 for various values of δ and L S .From these figures, it is evident that the value of concentration gradually increases as the dimensionless biofilm thickness δ decreases.Figures 4 to 5 represent the concentration   S x for various values of S L .From these this figures it is observed that, the value of the concentration increases when S L increases.When 1

Figure 1 .
Figure 1.Normalized concentration profile S(x) as a function of dimensionless distance x.The concentrations were computed using Eq.11 for various values of the δ when S L = 0.05, (-) denotes Eq.11 and (…) denotes the numerical simulation.

Figure 2 .
Figure 2. Normalized concentration profile S(x) as a function of dimensionless distance x.The concentrations were computed using Eq.11 for various values of the δ when S L = 0.5, (-) denotes Eq.11 and (…) denotes the numerical simulation.

Figure 3 .Figure 4 .
Figure 3. Normalized concentration profile S(x) as a function of dimensionless distance x.The concentrations were computed using Eq.11 for various values of the δ when S L = 5, (-) denotes Eq.11 and (…) denotes the numerical simulation.

Figure 5 .
Figure 5. Normalized concentration profile S(x) as a function of dimensionless distance x.The concentrations were computed using Eq.11 for various values of the S L when δ = 1, (-) denotes Eq.11 and (…) denotes the numerical simulation.concentration is uniform and the uniform value depends upon S L .It is clear that as dimensionless substrate concentration outside the biofilm S L increases when the value of dimensionless concentration increases.Eq.12 represents the normalized concentration flux into the biofilm.Figure 6 represents flux versus S L (dimensionless substrate concentration outside the biofilm).From this figure, it is inferred that the value of concentration flux decreases when the thickness of biofilm increases.Figure 7 represents flux versus

Figure 6
represents flux versus S L (dimensionless substrate concentration outside the biofilm).From this figure, it is inferred that the value of concentration flux decreases when the thickness of biofilm increases.Figure 7 represents flux versus   S x

Figure 6 .
Figure 6.Normalized concentration flux into the biofilm ψ as a function of dimensionless substrate concentration outside the biofilm S L .The concentrations were computed using Eq.12 for various values of the δ (-) denotes Eq.12 and (…) denotes the numerical solution.

Figure 7 .
Figure 7. Normalized concentration flux into the biofilm ψ as a function of dimensionless biofilm thickness δ.The concentrations were computed using Eq12 for various alues of the S L (-) denotes Eq.12 and (…) denotes the numerical simulation.the minimum value when .2 log 10

[ 9 ]
Dueck, I.G., Pyl'nik, S.V. and Min'kov, L.L. (2005) Modeling the evolution of a water remediation biofilm with account for its erosion.Biophysics, 50, 445-453.This work was supported by the Council of Scientific and Industrial Research (CSIR No. 01(2442)/10/EMR-II), Government of India.The authors are thankful to Dr. R. Murali, the Principal, the Mudura College, Madurai and Mr. S. Natanagobal, the Secretary, Madura College Board, Madurai for their encouragement.The author S. Muthukaruppan is very thankful to the Manonmaniam Sundaranar University, Tirunelveli for allowing to do the research work.

1 L
 on both sides of Eq.B.1 yields

Table 1 .
Comparison of normalized steady-state concentration S(x) with simulation results for various values of x and for some fixed values of S L = 0.5.

Table 2 .
Comparison of normalized steady-state concentration S(x) with simulation results for various values of x and for some fixed values of S L = 5.