Volume 26, Issue 1
A Novel Method for Solving Time-Dependent 2D Advection-Diffusion-Reaction Equations to Model Transfer in Nonlinear Anisotropic Media

Ji Lin ,  Sergiy Reutskiy ,  C. S. Chen and Jun Lu


Commun. Comput. Phys., 26 (2019), pp. 233-264.

Preview Full PDF BiBTex 345 994
  • Abstract

This paper presents a new numerical technique for solving initial and boundary value problems with unsteady strongly nonlinear advection diffusion reaction (ADR) equations. The method is based on the use of the radial basis functions (RBF) for the approximation space of the solution. The Crank-Nicolson scheme is used for approximation in time. This results in a sequence of stationary nonlinear ADR equations. The equations are solved sequentially at each time step using the proposed semianalytical technique based on the RBFs. The approximate solution is sought in the form of the analytical expansion over basis functions and contains free parameters. The basis functions are constructed in such a way that the expansion satisfies the boundary conditions of the problem for any choice of the free parameters. The free parameters are determined by substitution of the expansion in the equation and collocation in the solution domain. In the case of a nonlinear equation, we use the well-known procedure of quasilinearization. This transforms the original equation into a sequence of the linear ones on each time layer. The numerical examples confirm the high accuracy and robustness of the proposed numerical scheme.

  • History

Published online: 2019-02

  • AMS Subject Headings

65N35, 65N40, 65Y20

  • Cited by