Research Article | Open Access

Volume 2013 |Article ID 135140 | https://doi.org/10.1155/2013/135140

Fazle Mabood, Nopparat Pochai, "Asymptotic Solution for a Water Quality Model in a Uniform Stream", International Journal of Engineering Mathematics, vol. 2013, Article ID 135140, 4 pages, 2013. https://doi.org/10.1155/2013/135140

# Asymptotic Solution for a Water Quality Model in a Uniform Stream

Accepted02 Oct 2013
Published28 Nov 2013

#### Abstract

We employ approximate analytical method, namely, Optimal Homotopy Asymptotic Method (OHAM), to investigate a one-dimensional steady advection-diffusion-reaction equation with variable inputs arises in the mathematical modeling 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.

#### 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 . Several numerical and semianalytical methods have been developed for solving high order boundary value problems .

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  have provided a mathematical model of water pollution with the help of numerical method. Furthermore, Pochai and coworkers  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 .

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.

#### 2. Dispersion in a Stream

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

The boundary conditions are

#### 3. Basic Principles of OHAM

We review the basic principles of OHAM as illustrated in  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 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/m3. 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/m3 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 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 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 Runge-Kutta-Fehlberg-fourth fifth order method is shown in Figure 1; an excellent agreement can be observed.

 RKF45 OHAM 0.0 1.2500 1.2500 0.1 1.2031 1.2031 0.2 1.1580 1.1580 0.3 1.1147 1.1146 0.4 1.0731 1.0731 0.5 1.0333 1.0332 0.6 0.9952 0.9951 0.7 0.9588 0.9588 0.8 0.9242 0.9242 0.9 0.8914 0.8913 1.0 0.8604 0.8602 1.1 0.8313 0.8310 1.2 0.8042 0.8037 1.3 0.7792 0.7784 1.4 0.7564 0.7553 1.5 0.7362 0.7347 1.6 0.7187 0.7167 1.7 0.7042 0.7016 1.8 0.6832 0.6899 1.9 0.6862 0.6821

#### 5. Concluding Remarks

In this paper, we have presented the solution of the one-dimensional 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.

1. 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
2. S. Awang Kechil and I. Hashim, “Approximate analytical solution for MHD stagnation-point 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
3. 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
4. B. Bradie, A Friendly Introduction To Numerical Analysis, Pearson Education, Upper Saddle River, NJ, USA, 2006.
5. 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
6. 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
7. 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
8. S. J. Liao, The proposed homotopy analysis technique for the solution of nonlinear problem [Ph.D. thesis], Shanghai Jiao Tong University, 1992.
9. 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
10. N. Pochai and S. Tangmanee, “A mathematical model of water pollution using finite element method,” Contributions in Mathematics and Applications, East-West J. Math. Spec., 143–154, 2007. View at: Google Scholar
11. 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
12. 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
13. N. Pochai, “A numerical computation of the non-dimensional form of a non-linear hydrodynamic model in a uniform reservoir,” Nonlinear Analysis: Hybrid Systems, vol. 3, no. 4, pp. 463–466, 2009. View at: Publisher Site | Google Scholar
14. N. Pochai, “A numerical computation of a non-dimensional form of stream water quality model with hydrodynamic advection-dispersion-reaction equations,” Nonlinear Analysis: Hybrid Systems, vol. 3, no. 4, pp. 666–673, 2009. View at: Publisher Site | Google Scholar
15. N. Herişanu and V. Marinca, “Accurate analytical solutions to oscillators with discontinuities and fractional-power 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
16. 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
17. 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
18. 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 Transfer-Asian Research, 2013. View at: Publisher Site | Google Scholar
19. 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

#### More related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.