**Natural Science** Vol.5 No.4(2013), Article ID:29934,7 pages DOI:10.4236/ns.2013.54059

Mathematical modelling of a biofilm: The Adomian decomposition method

^{ }^{ } ^{}^{}

^{1}District Institute of Education and Training, Pudukkottai, India

^{2}Department of Mathematics, K. L. N. College of Engineering, Sivagangai, India

^{3}Department of Mathematics, The Madura College, Madurai, India; ^{*}Corresponding Author: raj_sms@rediffmail.com

Copyright © 2013 Subramanian Muthukaruppan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received 5 December 2012; revised 10 January 2013; accepted 20 February 2013

**Keywords:** Mathematical Modeling; Simulation; Adomian Decomposition Method; Initial Boundary Value Problems

ABSTRACT

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 S_{L}. These analytical results are compared with the available numerical results and are found to be in good agreement.

1. 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.

2. 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:

(1)

The boundary conditions are

(2)

(3)

The biomass balance [10].

(4)

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

(5)

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:

(6)

The above Eq.5 reduces to the following dimensionless form:

(7)

(8)

(9)

The dimensionless concentration flux into the biofilm is given by

(10)

3. 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-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-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:

(11)

The solution of concentration flux into the biofilm is obtained as

(12)

4. 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.

5. RESULTS AND DISCUSSION

An approximate analytical expression of concentrations S is given in the Eq.11. The concentration is plotted in Figures 1-3 for various values of δ and. 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 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, the

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.

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.

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. 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 4. 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 δ = 10, (—) denotes Eq.11 and (…) denotes the numerical simulation.

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. 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

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. 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.

6. 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. Furthermore, 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.

7. ACKNOWLEDGEMENTS

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.

REFERENCES

- Christensen, B.E. and Characklis, W.G. (1990) Physical and chemical properties of biofilms. In: Characklis, W.G. and Marshall, K.C., Eds., Biofilm, Wiley, New York, 93- 130.
- Skorokhodov, V.D. and Shestakova, S.I. (2004) Protection of nonmetalic materials from biocorrosion. Vysshaya Shkola, Moscow.
- Rittmann, B.E. and McCarty, P. L. (2000) Environmental biotechnology: Principles and applications. McGraw-Hill, Boston.
- Frank-Kamenetskii, D.A. (1987) Diffusion and heat transfer in chemical kinetics. Nauka, Moscow.
- Poltorak, O.M., Pryakhin, A.N. and Shhaitan, K.V. (1975) A general approach to solving kinetic problems. Vestnik Moskovskogo Universiteta, Seriya 2: Khimiya, 5, 536-543.
- Reznichenko, G.Yu. (2002) Lectures on mathematical models in biology. NITs Regulyarnaya i Khaoticheskaya Dinamika Part 1, Izhevsk.
- Dueck, J.H., Pyl’nik, S.V. and Min’kov, L.L. (2002) Modeling of the evolution of a biofilm with allowance for its erosion. In: Modeling of Processes in Synergetic Systems, Tomskogo Gosudarstvennogo Universiteta, Tomsk, 245-248.
- Dueck, J., Pylnik, S. and Minkov, L. (2004) Mathematical modelling of biofilm dynamics. In: Vestrate, W., Ed., Proceedings of European Symposium on Environmental Biotechnology ESEB, Balkema, Belgium, Leiden, 545-547.
- 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.
- Minkov, L.L., Pylnik, S.V. and Dueck, J.H. (2006) Steady state problem of substrate consumption in a biofilm for a square law of microbial death rate. Theoretical Foundations of Chemical Engineering, 40, 496-502. doi:10.1134/S004057950605006X
- Jaradat, O.K. (2008) Adomian decomposition method for solving Abelian differential equations. Journal of Applied Sciences, 8, 1962-1966. doi:10.3923/jas.2008.1962.1966
- Majid Wazwaz, A. and Gorguis, A. (2004) Solution of wave equation by Adomian decomposition method and the restrictions of the method. Applied Mathematics and Computation, 149, 807-814. doi:10.1016/S0096-3003(03)00186-3
- Makinde, O.D. (2007) Adomian decomposition approach to a SIR epidemic model with constant vaccination strategy. Applied Mathematics and Computation, 184, 842- 848. doi:10.1016/j.amc.2006.06.074
- Siddiquia, A.M., Hameedb, M., Siddiquic, B.M. and Ghoric, Q.K. (2010) Use of Adomian decomposition method in the study of parallel plate flow of a third grade fluid. Communications in Nonlinear Science and Numerical Simulation, 15, 2388-2399. doi:10.1016/j.cnsns.2009.05.073
- Mohamed, M.A. (2010) Comparison differential transformation technique with Adomian decomposition method for dispersive long-wave equations in (2+1)-dimensions. Applications & Applied Mathematics, 5, 148-166.

APPENDIX A

Basic Concept of the Adomian Decomposition Method (ADM)

Adomian decomposition method [9-13] depends on the non-linear differential Equation

(A.1)

into the two components

(A.2)

where L and N are the linear and non-linear parts of F respectively. The operator L is assumed to be an invertible operator. Solving for leads to

(A.3)

Applying the inverse operator L on both sides of Eq. A.3 yields

(A.4)

where is the constant of integration which satisfies the condition. Now assuming that the solution y can be represented as infinite series of the form

(A.5)

Furthermore, suppose that the non-linear term can be written as infinite series in terms of the Adomian polynomials of the form

(A.6)

where the Adomian polynomials of are evaluated using the formula:

(A.7)

Then substituting Eqs.A.5 and A.6 in Eq.A.4 gives

(A.8)

Then equating the terms in the linear system of Eq. A.8 gives the recurrent relation

(A.9)

However, in practice all the terms of series in Eq.A.7 cannot be determined, and the solution is approximated by the truncated series 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,

(B.1)

where and. Applying the inverse operator on both sides of Eq.B.1 yields

(B.2)

where A and B are the constants of integration. We let,

(B.3)

(B.4)

In view of Eqs.B.2-B.4, gives

(B.5)

We identify the zeroth component as

(B.6)

and the remaining components as the recurrence relation

(B.7)

where A_{n} are the Adomian polynomials of We can find the first few A_{n} as follows:

(B.8)

(B.9)

(B.10)

The remaining polynomials can be generated easily, and so,

(B.11)

(B.12)

(B.13)

(B.14)

Adding (B.11) to (B.14) we get Eq.11 in the text.

APPENDIX C

Scilab/Matlab Program to Find the Numerical Solution of Eqs.7-9

function pdex4 m = 0;

x = linspace(0,1);

t=linspace(0,1000);

sol = pdepe(m,@pdex4pde,@pdex4ic,@pdex4bc,x,t);

u1 = sol(:,:,1);

figure plot(x,u1(end,:))

title('u1(x,t)')

xlabel('Distance x')

ylabel('u1(x,2)')

function [c,f,s] = pdex4pde(x,t,u,DuDx)

c = 1;

f = 1 .* DuDx;

y = (u(1)/(1+u(1)))^2;

d=0.5;

s=-d*y;

function u0 = pdex4ic(x);

u0 =1;

APPENDIX D

Nomenclature

Symbols

: Microbial death constant, cm^{3}/(mg day)

: Diffusion coefficient within the biofilm, cm^{2}/day

: Substrate flux into the biofilm, (mg cm^{2})/day

: Michaelis-Menten constant, mg/cm^{3}

: Biofilm thickness,cm

: Substrate consumption rate constant, day^{−}^{1}

: Dimensionless substrate in the biofilm

: Substrate concentration in the biofilm, mg/cm^{3}

: Dimensionless substrate concentration outside the bioflim S_{1}: Substrate concentration outside the biofilm, mg/cm^{3}

T: Time, days

: Dimensionless co-ordinates

: Biomass yield per unit amount of substrate consumed, mg/mg

: Co-ordinate, cm

: Concentration of physiologically active microorganisms, mg/cm^{3}

: Dimensionless biofilm thickness

: Concentration flux into the biofilm