Efficient numerical treatment of nonlinearities in the advection-diffusion-reaction equations

Erdogan U., Sari M., Kocak H.

INTERNATIONAL JOURNAL OF NUMERICAL METHODS FOR HEAT & FLUID FLOW, cilt.29, sa.1, ss.132-145, 2019 (SCI-Expanded) identifier identifier

  • Yayın Türü: Makale / Tam Makale
  • Cilt numarası: 29 Sayı: 1
  • Basım Tarihi: 2019
  • Doi Numarası: 10.1108/hff-05-2017-0198
  • Dergi Adı: INTERNATIONAL JOURNAL OF NUMERICAL METHODS FOR HEAT & FLUID FLOW
  • Derginin Tarandığı İndeksler: Science Citation Index Expanded (SCI-EXPANDED), Scopus
  • Sayfa Sayıları: ss.132-145
  • Anahtar Kelimeler: Linearization, Advection-diffusion-reaction equations, Frechet derivative, Hybrid spline difference method, Space discretization, Time discretization, GENERALIZED BURGERS-FISHER, DIFFERENCE METHOD, HUXLEY, COLLOCATION
  • İstanbul Teknik Üniversitesi Adresli: Hayır

Özet

Purpose The purpose of this study is to propose a non-classical method to obtain efficient and accurate numerical solutions of the advection-diffusion-reaction equations. Design/methodology/approach Unlike conventional numerical methods, this study proposes a numerical scheme using outer Newton iteration applied to a time-dependent PDE. The linearized time dependent PDE is discretized by trapezoidal rule, which is second order in time, and by spline-based finite difference method of fourth order in space. Findings Using the proposed technique, even when relatively large time step sizes are used in computations, the efficiency of the proposed procedure is very clear for the numerical examples in comparison with the existing classical methods. Originality/value This study, unlike these classical methods, proposes an alternative approach based on linearizing the nonlinear problem at first, and then discretizing it by an appropriate scheme. This technique helps to avoid considering the convergence issues of Newton iteration applied to nonlinear algebraic system containing many unknowns at each time step if an implicit method is used in time discretization. The linearized PDE can be solved by implicit time integrator, which enables the use of large time step size.