Research Papers

Numerical Solution of Nonlinear Reaction–Advection–Diffusion Equation

[+] Author and Article Information
Anup Singh

Department of Mathematical Sciences,
Indian Institute of Technology (BHU),
Varanasi 221005, India

S. Das

Department of Mathematical Sciences,
Indian Institute of Technology (BHU),
Varanasi 221005, India
e-mail: sdas.apm@iitbhu.ac.in

S. H. Ong

Institute of Mathematical Sciences,
University of Malaya,
Kuala Lumpur 50603, Malaysia;
Department of Actuarial Science
and Applied Statistics,
UCSI University,
Kuala Lumpur 56000, Malaysia

H. Jafari

Department of Mathematical Sciences,
University of South Africa (UNISA),
Pretoria 0003, South Africa

1Corresponding author.

Contributed by the Design Engineering Division of ASME for publication in the JOURNAL OF COMPUTATIONAL AND NONLINEAR DYNAMICS. Manuscript received April 17, 2018; final manuscript received January 16, 2019; published online February 15, 2019. Assoc. Editor: Dumitru Baleanu.

J. Comput. Nonlinear Dynam 14(4), 041003 (Feb 15, 2019) (6 pages) Paper No: CND-18-1173; doi: 10.1115/1.4042687 History: Received April 17, 2018; Revised January 16, 2019

In the present article, the advection–diffusion equation (ADE) having a nonlinear type source/sink term with initial and boundary conditions is solved using finite difference method (FDM). The solution of solute concentration is calculated numerically and also presented graphically for conservative and nonconservative cases. The emphasis is given for the stability analysis, which is an important aspect of the proposed mathematical model. The accuracy and efficiency of the proposed method are validated by comparing the results obtained with existing analytical solutions for a conservative system. The novelty of the article is to show the damping nature of the solution profile due to the presence of the nonlinear reaction term for different particular cases in less computational time by using the reliable and efficient finite difference method.

Copyright © 2019 by ASME
Your Session has timed out. Please sign back in to continue.


Thangarajan, M. , 2007, Groundwater Models and Their Role in Assessment and Management of Groundwater Resources and Pollution, Springer, Dordrecht, The Netherlands.
Fried, J. J. , 1981, “Groundwater Pollution Mathematical Modelling: Improvement or Stagnation?,” Stud. Environ. Sci., 17, pp. 807–822. [CrossRef]
Bear, J. , and Verruijt, A. , 1987, Modeling Groundwater Flow and Pollution, Springer, Dordrecht, The Netherlands.
Gómez-Aguilar, J. F. , Miranda-Hernández, M. , López-López, M. G. , Alvarado-Martínez, V. M. , and Baleanu, D. , 2016, “Modeling and Simulation of the Fractional Space-Time Diffusion Equation,” Commun. Nonlinear Sci. Numer. Simul., 30(1–3), pp. 115–127. [CrossRef]
Jaiswal, S. , Chopra, M. , and Das, S. , 2018, “Numerical Solution of Two-Dimensional Solute Transport System Using Operational Matrices,” Transp. Porous, Media, 122(1), pp. 1–23. [CrossRef]
Yaseen, M. , Abbas, M. , Nazir, T. , and Baleanu, D. , 2017, “A Finite Difference Scheme Based on Cubic Trigonometric B-Splines for a Time Fractional Diffusion-Wave Equation,” Adv. Differ. Equations, 2017(1), p. 274. [CrossRef]
Wexler, E. J. , 1992, Analytical Solutions for One-, Two-, and Three-Dimensional Solute Transport in Ground-Water Systems With Uniform Flow, U.S. Government Printing Office, Denver, CO.
De Smedt, F. , 1984, “Analytical Solutions of the One-Dimensional Convective-Dispersive Solute Transport Equation,” Agric. Water Manage., 9(1), pp. 79–80. [CrossRef]
De La Hoz, F. , and Vadillo, F. , 2013, “The Solution of Two-Dimensional Advection-Diffusion Equations Via Operational Matrices,” Appl. Numer. Math., 72, pp. 172–187. [CrossRef]
Dehghan, M. , 2004, “Numerical Solution of the Three-Dimensional Advection-Diffusion Equation,” Appl. Math. Comput., 150(1), pp. 5–19.
San Jose Martinez, F. , Pachepsky, Y. A. , and Rawls, W. J. , 2010, “Modelling Solute Transport in Soil Columns Using Advective-Dispersive Equations With Fractional Spatial Derivatives,” Adv. Eng. Software, 41(1), pp. 4–8. [CrossRef]
Savović, S. , and Djordjevich, A. , 2012, “Finite Difference Solution of the One-Dimensional Advection-Diffusion Equation With Variable Coefficients in Semi-Infinite Media,” Int. J. Heat Mass Transfer, 55(15–16), pp. 4291–4294. [CrossRef]
Djordjevich, A. , and Savović, S. , 2013, “Solute Transport With Longitudinal and Transverse Diffusion in Temporally and Spatially Dependent Flow From a Pulse Type Source,” Int. J. Heat Mass Transfer, 65, pp. 321–326. [CrossRef]
Jaiswal, S. , Chopra, M. , Ong, S. H. , and Das, S. , 2017, “Numerical Solution of One-Dimensional Finite Solute Transport System With First Type Source Boundary Condition,” Int. J. Appl. Comput. Math., 3(4), pp. 3035–3045. [CrossRef]
Causon, D. M. , and Mingham, C. G. , 2010, Introductory Finite Difference Methods for PDEs, BookBoon, London.
LeVeque, R. J. , 2007, Finite Difference Methods for Ordinary and Partial Differential Equations, Society for Industrial and Applied Mathematics, Philadelphia, PA.
Siddiqi, S. S. , and Arshed, S. , 2014, “Quintic b-Spline for the Numerical Solution of the Good Boussinesq Equation,” J. Egypt. Math. Soc., 22(2), pp. 209–213. [CrossRef]
Seydaoğlu, M. , 2018, “An Accurate Approximation Algorithm for Burgers' Equation in the Presence of Small Viscosity,” J. Comput. Appl. Math., 344, pp. 473–481. [CrossRef]


Grahic Jump Location
Fig. 1

A partition of the (x,t)-plane into uniform cells of size h×k

Grahic Jump Location
Fig. 2

Comparison between analytical and numerical results at t =2.5, 10 when λ=0

Grahic Jump Location
Fig. 3

(a) Plots of u(x,t)/u0 versus x at t=2.5, 5, 10, 15, 20, D =0.6 and V =0.6 when λ=0 and (b) plots of u(x,t)/u0 versus x at t=2.5, 5, 10, 15, 20, D =0.6 and V =0.6 when λ≠0

Grahic Jump Location
Fig. 4

(a) Plots of u(x,t)/u0 versus x at t=1, 2.5, 5, 10,D =0.6 and V =0 when λ=0 and (b) plots of u(x,t)/u0 versus x at t=1,  2.5, 5, 10,D =0.6 and V =0 when λ≠0

Grahic Jump Location
Fig. 5

(a) Plots of u(x,t)/u0 versus x at t=2.5,  5,  10,  15,  20, D =1 and V =1 when λ=0 and (b) plots of u(x,t)/u0 versus x at t=2.5,  5,  10,  15,  20,D =1 and V =1 when λ≠0



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In