Biological system for elimination of volatile organics have been explored on experimental studies [1-3]. However, researches into the theoretical studies regarding biofilter models is rather limited. The pioneering contribution of Ottengraf and co-workers [1,2] of this model is based on some rather simplistic assumptions. The closed analytical expression have been used in validating-scale experimental data and actual design of pilot-scale biofilter units. Recently Zarook et al. [3], extended the work of Ottengraf et al. [1] and presented a detailed steady-state biofilteration model for single volume. Allen and Phatak [4] have extended their model to describe the biofilteration of VOC mixtures under steady-state conditions. Dehusses and Dunn [5] reports a transient biofilteration model which is based on the assumption that oxygen is in excess and the kinetics are of the Michaelis-Menten or Monod type. The steady-state model of Zarook et al. [3] was extended to describe the transient performance [6] of the biofilters. All the steady-state and transient biofilteration models [1-6] are based on the assumption that substrates are transported into the biofilm through diffusion. Biological systems for elimination of VOCs have been explored both on the experimental and mathematical modeling levels primarily in the Netherlands by Ottengraf et al. [6-8] followed by many researches even though land area requirements and lack of process control still restrict the industrial use of these systems. Several researchers [1,3,6,9-13] developed models to predict biodegradability of organic compounds in biofilters. The three general plans for biological treatment systems are biofilters, biotrickling filters and bioscrubbers. In biofilters, the porous medium is kept damp by maintaining the humidity of the incoming air and by occasional sprinkling. The reliability of biological processes and in particular of biofilteration for the treatment of waste gas streams containing VOC has been demonstrated by a very large number of experimental studies.

Recently Zarook et al. [14] obtained the concentration of VOC and oxygen only for the limiting cases (zeroorder kinetics and first-order kinetics) for monoid kinetics. However, to the best our knowledge, no analytical expressions pertaining to the steady state concentrations of VOC, oxygen and effectiveness factor have been reported. The purpose of this paper is to derive the analytical expression of concentration of VOC and oxygen for all values of parameters and all reaction mechanisms, using the Adomian decomposition method.

A steady-state biofilteration model (Figure 1) constitutes a set of mass balances within the biofilm. The mass balance equations in the biofilm are [14]:

with the boundary conditions

where and are the concentration of VOC and oxygen at a position x in the biofilm, and are the effective diffusion coefficient of VOC and oxygen in the biofilm. denotes biofilm density, and is the amount of biomass produced per amount of VOC consumed and amount of biomass produced per amount of oxygen consumed. For biological systems, The growth rate, for various reaction kinetics are given as follows:

Monod kinetics:

Andrews kinetics:

When oxygen limits the biodegradation rate, the growth rate is given by interactive model. The above Equations (5) and (6) are written as follows:

Interactive model from Monod kinetics:

Interactive model from Andrews kinetics:

In order to obtain numerical solution of model these equations are brought in dimensionless form through the dimensionless variables and groups. We make the above non-linear partial differential Equations (1) and (2) in dimensionless form by defining the following dimensionless parameters:

where and denotes the dimensionless concentration of VOC and oxygen, X is the dimensionless position in the biolayer. represents the Thiele modulus, , M, L and N are dimensionless constants. By substituting the Equation (9) in Equations (1) and (2), we can obtain the following dimensionless non-linear equation for Monod kinetics:

Using Equation (9), in the non-linear Equations (1) and (2), we can obtain the following dimensionless nonlinear equation for Andrews kinetics:

Using Equation (9), the dimensionless non-linear equation for Interactive model of Monod kinetics becomes

The dimensionless non-linear equation for Interactive model of Andrews kinetics of the Equations (1) and (2) is

Now the boundary condition in dimensionless form may be represented as follows:

For all the above cases, we can obtain the relation between and as follows:

Integrating the above equation twice and using the boundary condition (18) and (19) we get

When, the concentration of VOC and oxygen becomes equal. When, the concentration of oxygen. When, the concentration of VOC becomes one.

Nonlinear phenomena play a crucial role in physical chemistry and biology (heat and mass transfer, filtration of liquids, diffusion in chemical reactions, etc.). Constructing a particular, exact solution for these equations remains an important problem. Finding an exact solution that has a physicochemical or biological interpretation is of fundamental importance. This model is based on a non-stationary system of diffusion equations containing a nonlinear reaction term. It is not possible to solve these equations using standard analytical techniques. The investigation of an exact solution of nonlinear equations is interesting and important. In the past several decades, many authors mainly paid attention to studying the solution of nonlinear equations by using various methods, such as the Backlund and the Darboux transformation [15,16], the inverse scattering method [17], the bilinear method [18], the tanh method [19], the variational iteration method [20], the HPM [21-25], ADM [26-30]. The ADM was successfully applied to autonomous ordinary differential equations for nonlinear polycrystalline solids and to other fields. This method has been proved by many authors to be a powerful mathematical tool for various kinds of nonlinear problems. It is a promising and evolving method. The ADM is unique in its applicability, accuracy and efficiency. In this method, the solution procedure is very simple and only few iterations lead to highly accurate solutions that are valid for the whole solution domain. Using this method (see Appendix A), the concentration of VOC and oxygen for all the four cases can be obtained. The concentration of VOC for Monod kinetics in the biofilm is,

By solving the Equation (12), we can obtain the concentration of VOC for Andrews kinetics as follows,

Solving the dimensionless form of Interactive model from Monod kinetics Equation (14), we get the concentration as,

Solving the dimensionless form of Interactive model from Andrews kinetics Equation (16), we get the concentration as,

where

,

,

, and

.

The effectiveness factor is defined as the ratio of actual rate of reaction to the rate of reaction that would result if the entire biofilm was exposed to the concentration at the gas/biofilm interface. The effectiveness factor of various kinetics are as follows:

Monod kinetics

Andrews kinetics

Interactive model from Monod kinetics

Interactive model from Andrews kinetics

The dimensionless form of Equations (10)-(17) corresponding to the boundary conditions (18) and (19) were solved by numerical methods. We have used pdex4 to solve these equations (Pdex4 in MATLAB is a function to solve the initial-boundary value problems of differential equations. Matlab program to find the numerical solution of Equations (10) and (11) is given in the Appendix B. The numerical solution is compared with our analytical results and is shown in Figures 2-9. A satisfactory agreement is noticed for various values of the Thiele modulus and possible small values of reaction/ diffusion parameters.

Equations (21)-(24) represents the new analytical expressions of the concentration of VOC for Monoid, Andrews, Interactive Monoid and Interactive Andrews kinetics for all values of the parameter. Using the relation (Equation (20)), we can also obtain the concentration of oxygen for all the kinetics. Zarook et al. [17] obtained the

analytical expressions of concentration of VOC and oxygen only for the limiting cases (Zero kinetics and First-order kinetics). Concentration of VOC and oxygen depends upon the value of parameters, , L, N and.

The non-linear differential equations in biofilter models have been solved analytically for various kinetics using the Adomian decomposition method. Analytical expression of concentration of VOC and oxygen and corresponding effectiveness factor have been obtained for Monoid, Andrews, Interactive Monoid and Andrews kinetics and for all values of parameters. These analytical reactions very much useful for designing or scaling-up of biofilters.

This work was supported by the University Grants Commission (F. No. 39-58/2010(SR)), New Delhi, India and Council of Scientific and Industrial Research (CSIR No.: 01(2442)/10/EMR-II), New Delhi, India. The authors are thankful to Dr. R. Murali, The Principal, The Madura College, Madurai and The Secretary, Madura College Board, Madurai for their encouragement.

Analytical Solution of Non-Linear (Equation (16)) Using The Adomian Decomposition Method In the operator form, Equation (16) becomes

where

and

Applying to both sides of (B1) yields

where a and b are constants of integration. To solve (B3) by the Adomian method, we get

In view of the Equations (B4)-(B5), Equation (B3) gives

we identify the zeroth component as

and the remaining components as the recurrence relation,

where are the Adomian polynomials that represent the non-linear term in (B8).

Using (B9) we can find the first few as follows:

The remaining polynomials can be generated easily. The corresponding boundary condition becomes

Substitution of (B10) and (B11) in (B8) and operating with in conjunction with the boundary conditions (B12) in each case separately, we obtain

Substituting the Equations (B13)-(B15) in

we can obtain the Equation (24) in the text. Similarly, applying the above same procedure, we obtain the Equations (21)-(23).

Matlab/Scilab program to find the numerical solution of the Equations (10)-(11).