Asymptotic Solution for a Water Quality Model in a Uniform Stream

We employ approximate analytical method, namely, Optimal Homotopy Asymptotic Method (OHAM), to investigate a onedimensional steady advection-diffusion-reaction equationwith variable inputs arises in themathematicalmodeling of dispersion of pollutants in water is proposed. Numerical values are obtained via Runge-Kutta-Fehlberg fourth-fifth order method for comparison purpose. It was found that OHAM solution agrees well with the numerical solution. An example is included to demonstrate the efficiency, accuracy, and simplicity of the proposed method.


Introduction
Differential equations have been the focus of many studies due to their frequent appearance in various applications in physics, fluid mechanics, biology, and engineering.Consequently, considerable attention has been given to the solutions of higher order ordinary differential equations, integral equations, and fractional order partial differential equations of physical interest.Number of literatures concerning the application of higher order differential equations in nonlinear dynamics has grown rapidly in the recent years [1][2][3][4][5].Several numerical and semianalytical methods have been developed for solving high order boundary value problems [6][7][8][9].
A mathematical model for the dispersion of pollutants in a river is presented.The optimal homotopy asymptotic method for assessment of the chemical oxygen demand (COD) concentration in a river is considered.Pochai and Tangmanee [10] have provided a mathematical model of water pollution with the help of numerical method.Furthermore, Pochai and coworkers [11][12][13][14] have used numerical methods for the solution of hydrodynamic model with constant coefficients in the uniform reservoir and stream.
The optimal homotopy asymptotic method is an approximate analytical tool that is simple and straightforward and does not require the existence of any small or large parameter as does traditional perturbation method.Optimal Homotopy Asymptotic Method (OHAM) has been successfully applied to a number of nonlinear problems arising in fluid mechanics and heat transfer by various researchers [15][16][17][18][19].
This paper is organized as follows.First in Section 2, advection-diffusion-reaction equation is presented.In Section 3 we described the basic principles of OHAM.The OHAM solution of the problem is given in Section 4. Section 5 is devoted for the concluding remarks.

Dispersion in a Stream
The dispersion of chemical oxygen demand (COD) is described by the advection-diffusion-reaction equation (ADRE) [11] in the domain [, ]: where () is the concentration of COD at the point  ∈ [, ] (kg/m 3 ),  = () is the flow velocity in the  direction (m/s),   = () is the diffusion coefficient (m 2 /s),  = () is the substance decay rate (s −1 ), and  = () is the rate of change of substance concentration due to a source (kg/m 3 s).

International Journal of Engineering Mathematics
The boundary conditions are (2)

Basic Principles of OHAM
We review the basic principles of OHAM as illustrated in [3] and other works.
(i) Consider the following differential equation: where where ,  are linear and nonlinear operators, V() is an unknown function, and () is a known function.
(ii) Construct an optimal homotopy equation as where 0 ≤  ≤ 1 is an embedding parameter and () = ∑  =1     is auxiliary function on which the convergence of the solution greatly dependents.The auxiliary function () also adjusts the convergence domain and controls the convergence region.
For nonlinear problems, generally, this will not be the case.For determining   ( = 1, 2, . . ., ), collocation method, Ritz method, or the method of least squares can be used.
(v) Finally, substitute these constants in (6) and one can get the approximate solution.

Application of OHAM
Consider the advection-dispersion-reaction equation (1) in the form We assume that there is a plant which discharges waste water into the channel at the starting point 0.0 km and that the COD concentrations of the waste water are 1.2500 kg/m 3 .Let the physical parameter values be diffusion coefficient 2, flow velocity () = 5 −  m/s, where  ∈ [0, 2], substance decay rate 3 s −1 , and rate of change of substance concentration due to the source 1 kg/m 3 s; we can obtain variable coefficients of convection-diffusion equation (8) as Equation ( 8) becomes subject to the boundary conditions: where primes denote differentiation with respect to .
According to OHAM, we have Zeroth order problem is with boundary conditions: The solution of ( 13) with boundary conditions ( 14) is First order problem is with boundary conditions: The solution of ( 16) with boundary conditions (17) is 343 (6 The terms of second order problem and its solution are too large to be written above; therefore the final three-term solution via OHAM for  = 1 is We use the method of least squares to obtain  1 and  2 , the unknown convergent constants in c.The values of the convergent constants are  1 = 0.001881235287,  2 = −1.103025681. By substituting the values of  1 and  2 in (19) and after simplification, we obtain the second order approximate solution via OHAM.To check the accuracy of the OHAM solution, a comparison between the solutions obtained by OHAM and numerical method was made and is presented in Table 1.Graphical representation of the solution using OHAM and Runge-Kutta-Fehlberg-fourth fifth order method is shown in Figure 1; an excellent agreement can be observed.

Concluding Remarks
In this paper, we have presented the solution of the onedimensional steady advection-diffusion-reaction equation with variable inputs using homotopy approach and Runge-Kutta-Fehlberg fourth-fifth order method.Both approximate analytical and numerical results are obtained for the given problem.The validity of the proposed procedure, called the Optimal Homotopy Asymptotic Method (OHAM), was demonstrated on an example, and very good agreement was found between the approximate analytic results and numerical simulation results.The proposed scheme provides us with a simple and accurate way to optimally control and adjust the convergence of a solution and can give very good approximations in a few terms.

Table 1 :
Comparison of () via OHAM and numerical method.