Research Article  Open Access
Asymptotic Solution for a Water Quality Model in a Uniform Stream
Abstract
We employ approximate analytical method, namely, Optimal Homotopy Asymptotic Method (OHAM), to investigate a onedimensional steady advectiondiffusionreaction equation with variable inputs arises in the mathematical modeling of dispersion of pollutants in water is proposed. Numerical values are obtained via RungeKuttaFehlberg fourthfifth 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.
1. 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–5]. Several numerical and semianalytical methods have been developed for solving high order boundary value problems [6–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–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–19].
This paper is organized as follows. First in Section 2, advectiondiffusionreaction 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.
2. Dispersion in a Stream
The dispersion of chemical oxygen demand (COD) is described by the advectiondiffusionreaction 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).
The boundary conditions are
3. 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 is problem domain, , where , are linear and nonlinear operators, is an unknown function, and is a known function.(ii)Construct an optimal homotopy equation as where is an embedding parameter and 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.(iii)Expand in Taylor’s series about ; one has an approximate solution: Many researchers have observed that the convergence of the series in (5) depends upon ; if it is convergent then, we obtain (iv)Substitute (6) in (3); we have the following residual: If , then will be the exact solution. For nonlinear problems, generally, this will not be the case. For determining , 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.
4. Application of OHAM
Consider the advectiondispersionreaction 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 m/s, where , 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 convectiondiffusion 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 The terms of second order problem and its solution are too large to be written above; therefore the final threeterm solution via OHAM for is
We use the method of least squares to obtain and , the unknown convergent constants in . The values of the convergent constants are , .
By substituting the values of and 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 RungeKuttaFehlbergfourth fifth order method is shown in Figure 1; an excellent agreement can be observed.

5. Concluding Remarks
In this paper, we have presented the solution of the onedimensional steady advectiondiffusionreaction equation with variable inputs using homotopy approach and RungeKuttaFehlberg fourthfifth 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.
References
 A.M. Wazwaz, “The modified decomposition method and Padé approximants for a boundary layer equation in unbounded domain,” Applied Mathematics and Computation, vol. 177, no. 2, pp. 737–744, 2006. View at: Publisher Site  Google Scholar
 S. Awang Kechil and I. Hashim, “Approximate analytical solution for MHD stagnationpoint flow in porous media,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1346–1354, 2009. View at: Publisher Site  Google Scholar
 V. Marinca and N. Herisanu, “Optimal homotopy perturbation method for strongly nonlinear differential equations,” Nonlinear Science Letters A, vol. 1, no. 3, pp. 273–280, 2010. View at: Google Scholar
 B. Bradie, A Friendly Introduction To Numerical Analysis, Pearson Education, Upper Saddle River, NJ, USA, 2006.
 F. Mabood, A. I. M. Ismail, and I. Hashim, “The application of optimal homotopy asymptotic method for the approximate solution of Riccati equation,” Sains Malaysiana, vol. 42, no. 6, pp. 863–867, 2013. View at: Google Scholar
 G. Adomian, “A new approach to nonlinear partial differential equations,” Journal of Mathematical Analysis and Applications, vol. 102, no. 2, pp. 420–434, 1984. View at: Google Scholar
 J.H. He, “Asymptotology by homotopy perturbation method,” Applied Mathematics and Computation, vol. 156, no. 3, pp. 591–596, 2004. View at: Publisher Site  Google Scholar
 S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problem [Ph.D. thesis], Shanghai Jiao Tong University, 1992.
 V. Marinca and N. Herişanu, “Application of optimal homotopy asymptotic method for solving nonlinear equations arising in heat transfer,” International Communications in Heat and Mass Transfer, vol. 35, no. 6, pp. 710–715, 2008. View at: Publisher Site  Google Scholar
 N. Pochai and S. Tangmanee, “A mathematical model of water pollution using finite element method,” Contributions in Mathematics and Applications, EastWest J. Math. Spec., 143–154, 2007. View at: Google Scholar
 N. Pochai and R. Depana, “An optimal control of water pollution in a stream using a finite difference method,” World Academy of Science, Engineering and Technology, vol. 6, no. 56, pp. 1186–1188, 2011. View at: Google Scholar
 N. Pochai, S. Tangmanee, L. J. Crane, and J. J. H. Miller, “A Mathematical model of water pollution using finite element method,” Proceedings in Applied Mathematics and Mechanics, vol. 6, pp. 755–756, 2006. View at: Google Scholar
 N. Pochai, “A numerical computation of the nondimensional form of a nonlinear hydrodynamic model in a uniform reservoir,” Nonlinear Analysis: Hybrid Systems, vol. 3, no. 4, pp. 463–466, 2009. View at: Publisher Site  Google Scholar
 N. Pochai, “A numerical computation of a nondimensional form of stream water quality model with hydrodynamic advectiondispersionreaction equations,” Nonlinear Analysis: Hybrid Systems, vol. 3, no. 4, pp. 666–673, 2009. View at: Publisher Site  Google Scholar
 N. Herişanu and V. Marinca, “Accurate analytical solutions to oscillators with discontinuities and fractionalpower restoring force by means of the optimal homotopy asymptotic method,” Computers and Mathematics with Applications, vol. 60, no. 6, pp. 1607–1615, 2010. View at: Publisher Site  Google Scholar
 S. Islam, F. Mabood, G. Zaman, X.Z. Li, and I. H. Jung, “Optimal homotopy asymptotic method solution to convection heat transfer flow,” International Journal of Physical Sciences, vol. 6, no. 23, pp. 5511–5519, 2011. View at: Google Scholar
 S. Islam, R. Ali Shah, I. Ali, and N. M. Allah, “Optimal homotopy asymptotic solutions of couette and poiseuille flows of a third grade fluid with heat transfer analysis,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 11, no. 6, pp. 389–400, 2010. View at: Google Scholar
 F. Mabood, W. A. Khan, and A. I. M. Ismail, “Optimal homotopy asymptotic method for heat transfer in hollow sphere with robin boundary conditions,” Heat TransferAsian Research, 2013. View at: Publisher Site  Google Scholar
 F. Mabood, W. A. Khan, and A. I. M. Ismail, “Analytical solution for radiation effects on heat transfer in blasius flow,” International Journal of Engineering Science, vol. 2, no. 2, pp. 63–72, 2013. View at: Google Scholar
Copyright
Copyright © 2013 Fazle Mabood and Nopparat Pochai. 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.