Uniformly convergent novel finite difference methods for singularly perturbed reaction–diffusion equations

In this paper, we construct a kind of novel finite difference (NFD) method for solving singularly perturbed reaction–diffusion problems. Different from directly truncating the high‐order derivative terms of the Taylor's series in the traditional finite difference method, we rearrange the Taylor's expansion in a more elaborate way based on the original equation to develop the NFD scheme for 1D problems. It is proved that this approach not only can highly improve the calculation accuracy but also is uniformly convergent. Then, applying alternating direction implicit technique, the newly deduced schemes are extended to 2D equations, and the uniform error estimation based on Shishkin mesh is derived, too. Finally, numerical experiments are presented to verify the high computational accuracy and theoretical prediction.     

Nonstandard finite difference schemes for Michaelis–Menten type reaction‐diffusion equations

We compare and investigate the performance of the exact scheme of the Michaelis–Menten (M–M) ordinary differential equation with several new nonstandard finite difference (NSFD) schemes that we construct using Mickens' rules. Furthermore, the exact scheme of the M–M equation is used to design several dynamically consistent NSFD schemes for related reaction‐diffusion equations, advection‐reaction equations, and advection‐reaction‐diffusion equations. Numerical simulations that support the theory and demonstrate computationally the power of NSFD schemes are presented. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Parallel characteristic finite difference method for convection–diffusion equations

Based on the overlapping domain decomposition, an efficient parallel characteristic finite difference scheme is proposed for solving convection‐diffusion equations numerically. We give the optimal convergence order in error estimate analysis, which shows that we just need to iterate once or twice at each time level to reach the optimal convergence order. Numerical experiments also confirm the theoretical analysis. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 854–866, 2011     

Nonstandard finite difference schemes for reaction‐diffusion equations

We construct finite difference schemes for a particular class of one‐space dimension, nonlinear reaction‐diffusion PDEs. The use of nonstandard finite difference methods and the imposition of a positivity condition constrain the schemes to be explicit and allow the determination of functional relations between the space and time step‐sizes. The general procedure is illustrated by applying it to several important model systems of PDEs © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 201–214, 1999     

Lyapunov functionals for reaction–diffusion equations with memory

We consider a reaction‐diffusion equation in which the usual diffusion term also depends on the past history of the diffusion itself. This equation has been analysed by several authors, with an emphasis on the longtime behaviour of the solutions. In this respect, the first results have been obtained by using the past history approach. They show that the equation, subject to a suitable boundary condition, defines a dissipative dynamical system which possesses a global attractor. A similar theorem has been recently proved by Chepyzhov and Miranville, using a different method based on the notion of trajectory attractors. In addition, those authors provide sufficient conditions that ensure the existence of a Lyapunov functional. Here we show that a similar result can be demonstrated within the past history approach, with less restrictive conditions. Copyright © 2005 John Wiley & Sons, Ltd.     

A study on a second order finite difference scheme for fractional advection–diffusion equations

Second order finite difference schemes for fractional advection–diffusion equations are considered in this paper. We note that, when studying these schemes, advection terms with coefficients having the same sign as those of diffusion terms need additional estimates. In this paper, by comparing generating functions of the corresponding discretization matrices, we find that sufficiently strong diffusion can dominate the effects of advection. As a result, convergence and stability of schemes are obtained in this situation.     

Convergence of an implicit Voronoi finite volume method for reaction–diffusion problems

We investigate the convergence of an implicit Voronoi finite volume method for reaction–diffusion problems including nonlinear diffusion in two space dimensions. The model allows to handle heterogeneous materials and uses the chemical activities of the involved species as primary variables. The numerical scheme works with boundary conforming Delaunay meshes and preserves positivity and the dissipative property of the continuous system. Starting from a result on the global stability of the scheme (uniform, mesh‐independent global upper, and lower bounds), we prove strong convergence of the chemical activities and their gradients to a weak solution of the continuous problem. To illustrate the preservation of qualitative properties by the numerical scheme, we present a long‐term simulation of the Michaelis–Menten–Henri system. Especially, we investigate the decay properties of the relative free energy over several magnitudes of time, and obtain experimental orders of convergence for this quantity. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 141–174, 2016     

Numerical modelling of linear and nonlinear diffusion equations by compact finite difference method

In this work, accurate solutions to linear and nonlinear diffusion equations were introduced. A combination of a sixth-order compact finite difference scheme in space and a low-storage third-order total variation diminishing Runge-Kutta scheme in time have been used for treatment of these equations. The computed results with the use of this technique have been compared with the exact solution to show the accuracy of it. Here, the approximate solution to the diffusion equations has been obtained easily and elegantly with neither transforming nor linearizing the equation. The present method is seen to be a very good alternative method to some existing techniques for realistic problems.     

Superconvergence of finite difference approximations for convection–diffusion problems

An Erratum has been published for this article in Numerical Linear Algebra with Applications 8 (4) 2001, iii–iv. We are concerned with numerical solutions of convection–diffusion equations. The convergence behaviour of numerical solutions is considered by using the finite difference approximation with respect to spatial variables and implicit method with respect to time variable. It is shown that superconvergence occurs near a part of the boundary which has Dirichlet's data. Copyright © 2001 John Wiley & Sons, Ltd.     

A second‐order linearized difference scheme for a strongly coupled reaction‐diffusion system

This article deals with the numerical solution to some models described by the system of strongly coupled reaction–diffusion equations with the Neumann boundary value conditions. A linearized three‐level scheme is derived by the method of reduction of order. The uniquely solvability and second‐order convergence in L₂‐norm are proved by the energy method. A numerical example is presented to demonstrate the accuracy and efficiency of the proposed method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Numerical and theoretical study of weak Galerkin finite element solutions of Turing patterns in reaction–diffusion systems

In this paper, we introduce numerical schemes and their analysis based on weak Galerkin finite element framework for solving 2‐D reaction–diffusion systems. Weak Galerkin finite element method (WGFEM) for partial differential equations relies on the concept of weak functions and weak gradients, in which differential operators are approximated by weak forms through the Green's theorem. This method allows the use of totally discontinuous functions in the approximation space. In the current work, the WGFEM solves reaction–diffusion systems to find unknown concentrations (u, v) in element interiors and boundaries in the weak Galerkin finite element space WG(P₀, P₀, RT₀) . The WGFEM is used to approximate the spatial variables and the time discretization is made by the backward Euler method. For reaction–diffusion systems, stability analysis and error bounds for semi‐discrete and fully discrete schemes are proved. Accuracy and efficiency of the proposed method successfully tested on several numerical examples and obtained results satisfy the well‐known result that for small values of diffusion coefficient, the steady state solution converges to equilibrium point. Acquired numerical results asserted the efficiency of the proposed scheme.     

Nonstandard finite difference schemes for reaction‐‐diffusion equations having linear advection

We extend previous work on nonstandard finite difference schemes for one‐space dimension, nonlinear reaction–diffusion PDEs to the case where linear advection is included. The use of a positivity condition allows the determination of a functional relation between the time and space step‐sizes, and provides schemes that are explicit. The Fisher equation is used to illustrate the method. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 361–364, 2000     

Asymptotic stability of traveling wave fronts in nonlocal reaction–diffusion equations with delay

This paper is concerned with the asymptotic stability of traveling wave fronts of a class of nonlocal reaction–diffusion equations with delay. Under monostable assumption, we prove that the traveling wave front is exponentially stable by means of the (technical) weighted energy method, when the initial perturbation around the wave is suitable small in a weighted norm. The exponential convergent rate is also obtained. Finally, we apply our results to some population models and obtain some new results, which recover, complement and/or improve a number of existing ones.     

Solving the reaction–diffusion equations with nonlocal boundary conditions based on reproducing kernel space

The reaction–diffusion equations with initial condition and nonlocal boundary conditions are discussed in this article. A reproducing kernel space is constructed, in which an arbitrary function satisfies the initial condition and nonlocal boundary conditions of the reaction‐diffusion equations. Based on the reproducing kernel space, a new algorithm for solving the reaction–diffusion equations with initial condition and nonlocal boundary conditions is presented. Some examples are displayed to demonstrate the validity and applicability of the proposed method. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Normal deviations from the averaged motion for some reaction–diffusion equations with fast oscillating perturbation

We study the normalized difference between the solution u of a reaction–diffusion equation in a bounded interval [0,L], perturbed by a fast oscillating term arising as the solution of a stochastic reaction–diffusion equation with a strong mixing behavior, and the solution of the corresponding averaged equation. We assume the smoothness of the reaction coefficient and we prove that a central limit type theorem holds. Namely, we show that the normalized difference converges weakly in C([0,T];L²(0,L)) to the solution of the linearized equation, where an extra Gaussian term appears. Such a term is explicitly given.     

Traveling wave solutions of n‐dimensional delayed reaction–diffusion systems and application to four‐dimensional predator–prey systems

This paper deals with the existence of traveling wave solutions for n‐dimensional delayed reaction–diffusion systems. By using Schauder's fixed point theorem, we establish the existence result of a traveling wave solution connecting two steady states by constructing a pair of upper–lower solutions that are easy to construct. As an application, we apply our main results to a four‐dimensional delayed predator–prey system and obtain the existence of traveling wave solutions. Copyright © 2016 John Wiley & Sons, Ltd.     

Adequacy of finite difference schemes for convection‐diffusion equations

Finite difference schemes for the numerical solution of singularly perturbed convection problems on uniform grids are studied in the limit case where the viscosity and the meshsize approach zero at the same time. The present error estimates are given in terms of order of magnitude in the above limit process and are useful in a priori choosing adequate schemes and meshsizes for boundary‐layer problems and problems with closed characteristics. Published 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 280–295, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10007     

Numerical simulation of reaction–diffusion equations on spherical domains

The two-variable reaction diffusion equations on the spherical domain is considered and simulated, using the semi-implicit Euler finite difference method. It is shown that the method keeps the kinetics from overshooting the stable branches when a large time step is used in the simulation.     

Spike patterns in a reaction–diffusion ODE model with Turing instability

We explore a mechanism of pattern formation arising in processes described by a system of a single reaction–diffusion equation coupled with ordinary differential equations. Such systems of equations arise from the modeling of interactions between cellular processes and diffusing growth factors. We focus on the model of early carcinogenesis proposed by Marciniak‐Czochra and Kimmel, which is an example of a wider class of pattern formation models with an autocatalytic non‐diffusing component. We present a numerical study showing emergence of periodic and irregular spike patterns because of diffusion‐driven instability. To control the accuracy of simulations, we develop a numerical code on the basis of the finite‐element method and adaptive mesh grid. Simulations, supplemented by numerical analysis, indicate a novel pattern formation phenomenon on the basis of the emergence of nonstationary structures tending asymptotically to a sum of Dirac deltas. Copyright © 2013 John Wiley & Sons, Ltd.     

The decomposition method applied to systems of partial differential equations and to the reaction–diffusion Brusselator model

In this work, systems of linear and nonlinear partial differential equations and the reaction–diffusion Brusselator model are handled by applying the decomposition method. The advantage of this work is twofold. Firstly, the decomposition method reduces the computational work. Secondly, in comparison with existing techniques, the decomposition method is an improvement with regard to its accuracy and rapid convergence. The decomposition method has the advantage of being more concise for analytical and numerical purposes.