Numerical solution of linear one dimensional Fisher’s reaction diffusion equations by finite difference methods

2021 ◽  
Author(s):  
Nur Aishah Jamaludin ◽  
Amirah Azmi
Keyword(s):  
Finite Difference ◽  
Numerical Solution ◽  
Finite Difference Methods ◽  
Reaction Diffusion ◽  
Reaction Diffusion Equations ◽  
Diffusion Equations ◽  
One Dimensional ◽  
Difference Methods

A conservative, piecewise–analytical, transversal method of lines for reaction–diffusion equations

2019 ◽  
Vol 29 (11) ◽  
pp. 4093-4129 ◽  
Author(s):  
J.I. Ramos
Keyword(s):  
Finite Difference ◽  
Boundary Layers ◽  
Diffusion Coefficients ◽  
Method Of Lines ◽  
Reaction Diffusion ◽  
Reaction Rates ◽  
Reaction Diffusion Equations ◽  
Diffusion Equations ◽  
Content Type ◽  
One Dimensional

Purpose The purpose of this paper is to develop a new transversal method of lines for one-dimensional reaction–diffusion equations that is conservative and provides piecewise–analytical solutions in space, analyze its truncation errors and linear stability, compare it with other finite-difference discretizations and assess the effects of the nonlinear diffusion coefficients, reaction rate terms and initial conditions on wave propagation and merging. Design/methodology/approach A conservative, transversal method of lines based on the discretization of time and piecewise analytical integration of the resulting two-point boundary-value problems subject to the continuity of the dependent variables and their fluxes at the control-volume boundaries, is presented. The method provides three-point finite difference expressions for the nodal values and continuous solutions in space, and its accuracy has been determined first analytically and then assessed in numerical experiments of reaction-diffusion problems, which exhibit interior and/or boundary layers. Findings The transversal method of lines presented here results in three-point finite difference equations for the nodal values, treats the diffusion terms implicitly and is unconditionally stable if the reaction terms are treated implicitly. The method is very accurate for problems with the interior and/or boundary layers. For a system of two nonlinearly-coupled, one-dimensional reaction–diffusion equations, the formation, propagation and merging of reactive fronts have been found to be strong function of the diffusion coefficients and reaction rates. For asymmetric ignition, it has been found that, after front merging, the temperature and concentration profiles are almost independent of the ignition conditions. Originality/value A new, conservative, transversal method of lines that treats the diffusion terms implicitly and provides piecewise exponential solutions in space without the need for interpolation is presented and applied to someone.

Sign in / Sign up

Export Citation Format

Share Document