Quasisolutions and dynamics of time-periodic nonquasimonotone reaction-diffusion systems

In this paper, using the upper-lower solution method, we investigate the properties of periodic quasisolutions for time-periodic nonquasimonotone reaction-diffusion systems, and present some existence, asymptotic behavior results for two ecological models.     

An almost third order finite difference scheme for singularly perturbed reaction-diffusion systems

This paper addresses the numerical approximation of solutions to coupled systems of singularly perturbed reaction-diffusion problems. In particular a hybrid finite difference scheme of HODIE type is constructed on a piecewise uniform Shishkin mesh. It is proved that the numerical scheme satisfies a discrete maximum principle and also that it is third order (except for a logarithmic factor) uniformly convergent, even for the case in which the diffusion parameter associated with each equation of the system has a different order of magnitude. Numerical examples supporting the theory are given.     

非线性时滞反应扩散有限差分方程组的高阶单调迭代方法

对一类非线性时滞反应扩散方程的有限差分方程组建立了一类高阶单调迭代方法.这类方法给出了一个有效的线性迭代算法.迭代序列单调收敛于方程组的唯一解,并且序列的单调性使得每一步迭代都给出了解的改进的上下界.迭代收敛率具有p+2阶,这里p≥1是一个正整数,它依赖于迭代方法的构造.数值结果显示了方法的有效性.     

A new fitted operator finite difference method to solve systems of evolutionary reaction-diffusion equations

Abstract

In recent years, fitted operator finite difference methods (FOFDMs) have been developed for numerous types of singularly perturbed ordinary differential equations. The construction of most of these methods differed though the final outcome remained similar. The most crucial aspect was how the difference operator was designed to approximate the differential operator in question. Very often the approaches for constructing these operators had limited scope in the sense that it was difficult to extend them to solve even simple one-dimensional singularly perturbed partial differential equations. However, in some of our most recent work, we have successfully designed a class of FOFDMs and extended them to solve singularly perturbed time-dependent partial differential equations. In this paper, we design and analyze a robust FOFDM to solve a system of coupled singularly perturbed parabolic reaction-diffusion equations. We use the backward Euler method for the semi-discretization in time. An FOFDM is then developed to solve the resulting set of boundary value problems. The proposed method is analyzed for convergence. Our method is uniformly convergent with order one and two, respectively, in time and space, with respect to the perturbation parameters. Some numerical experiments supporting the theoretical investigations are also presented.     

Monotone iterative methods for finite difference system of reaction-diffusion equations

Summary This paper presents an existence-comparison theorem and an iterative method for a nonlinear finite difference system which corresponds to a class of semilinear parabolic and elliptic boundary-value problems. The basic idea of the iterative method for the computation of numerical solutions is the monotone approach which involves the notion of upper and lower solutions and the construction of monotone sequences from a suitable linear discrete system. Using upper and lower solutions as two distinct initial iterations, two monotone sequences from a suitable linear system are constructed. It is shown that these two sequences converge monotonically from above and below, respectively, to a unique solution of the nonlinear discrete equations. This formulation leads to a well-posed problem for the nonlinear discrete system. Applications are given to several models arising from physical, chemical and biological systems. Numerical results are given to some of these models including a discussion on the rate of convergence of the monotone sequences.     

Finite difference reaction-diffusion systems with coupled boundary conditions and time delays

This paper is concerned with finite difference solutions of a coupled system of reaction-diffusion equations with nonlinear boundary conditions and time delays. The system is coupled through the reaction functions as well as the boundary conditions, and the time delays may appear in both the reaction functions and the boundary functions. The reaction-diffusion system is discretized by the finite difference method, and the investigation is devoted to the finite difference equations for both the time-dependent problem and its corresponding steady-state problem. This investigation includes the existence and uniqueness of a finite difference solution for nonquasimonotone functions, monotone convergence of the time-dependent solution to a maximal or a minimal steady-state solution for quasimonotone functions, and local and global attractors of the time-dependent system, including the convergence of the time-dependent solution to a unique steady-state solution. Also discussed are some computational algorithms for numerical solutions of the steady-state problem when the reaction function and the boundary function are quasimonotone. All the results for the coupled reaction-diffusion equations are directly applicable to systems of parabolic-ordinary equations and to reaction-diffusion systems without time delays.     

Accelerated monotone iterative methods for finite difference equations of reaction-diffusion

This paper is concerned with numerical methods for a finite difference system of reaction-diffusion-convection equation under nonlinear boundary condition. Various monotone iterative methods are presented, and each of these methods leads to an existence-comparison theorem as well as a computational algorithm for numerical solutions. The monotone property of the iterations gives improved upper and lower bounds of the solution in each iteration, and the rate of convergence of the iterations is either quadratic or nearly quadratic depending on the property of the nonlinear function. Application is given to a model problem from chemical engineering, and some numerical results, including a test problem with known analytical solution, are presented to illustrate the various rates of convergence of the iterations. Received November 2, 1995 / Revised version received February 10, 1997     

Nonlinear stability in reaction-diffusion systems via optimal Lyapunov functions

We define optimal Lyapunov functions to study nonlinear stability of constant solutions to reaction-diffusion systems. A computable and finite radius of attraction for the initial data is obtained. Applications are given to the well-known Brusselator model and a three-species model for the spatial spread of rabies among foxes.     

On stability of interconnected finite difference systems

We give a method of construction of Lyapunov functions in the form of a linear form with respect to moduli of variables, for which there exist Krasovskii constants in the case of asymptotic stability, for linear systems with constant coefficients and some types of nonlinear systems of finite-difference equations. An application of the above functions as components of a vector Lyapunov function allowed us to obtain conditions on asymptotic stability for interrelated finite-difference systems.Translated from Dinamicheskie Sistemy, No. 8, pp. 68–71, 1989.     

A high order HODIE finite difference scheme for 1D parabolic singularly perturbed reaction-diffusion problems

This paper deals with the numerical approximation of the solution of 1D parabolic singularly perturbed problems of reaction-diffusion type. The numerical method combines the standard implicit Euler method on a uniform mesh to discretize in time and a HODIE compact fourth order finite difference scheme to discretize in space, which is defined on a priori special meshes condensing the grid points in the boundary layer regions. The method is uniformly convergent having first order in time and almost fourth order in space. The analysis of the uniform convergence is made in two steps, splitting the contribution to the error from the time and the space discretization. Although this idea has been previously used to prove the uniform convergence for parabolic singularly perturbed problems, here the proof is based on a new study of the asymptotic behavior of the exact solution of the semidiscrete problems obtained after the time discretization by using the Euler method. Some numerical results are given corroborating in practice the theoretical results.     

Efficient finite difference derivatives for algebraic modelling systems

A procedure is described which determines Jacobian incidence structure and the constant/nonconstant nature of each Jacobian element via examination of the text of function expression strings. This procedure may be used to minimize the effort required to evaluate by finite differences the Jacobian of a set of functions. Target applications involve algebraic modelling systems and other systems with interpreted functions which require evaluation of first derivatives. Computational experience is presented and discussed.     

A finite difference method on layer-adapted meshes for an elliptic reaction-diffusion system in two dimensions

An elliptic system of singularly perturbed linear reaction-diffusion equations, coupled through their zero-order terms, is considered on the unit square. This system does not in general satisfy a maximum principle. It is solved numerically using a standard difference scheme on tensor-product Bakhvalov and Shishkin meshes. An error analysis for these numerical methods shows that one obtains nodal convergence on the Bakhvalov mesh and convergence on the Shishkin mesh, where mesh intervals are used in each coordinate direction and the convergence is uniform in the singular perturbation parameter. The analysis is much simpler than previous analyses of similar problems, even in the case of a single reaction-diffusion equation, as it does not require the construction of an elaborate decomposition of the solution. Numerical results are presented to confirm our theoretical error estimates.

     

On the uniform convergence of a finite difference scheme for time dependent singularly perturbed reaction-diffusion problems

In this work we are interested in the numerical approximation of 1D parabolic singularly perturbed problems of reaction-diffusion type. To approximate the multiscale solution of this problem we use a numerical scheme combining the classical backward Euler method and central differencing. The scheme is defined on some special meshes which are the tensor product of a uniform mesh in time and a special mesh in space, condensing the mesh points in the boundary layer regions. In this paper three different meshes of Shishkin, Bahkvalov and Vulanovic type are used, proving the uniform convergence with respect to the diffusion parameter. The analysis of the uniform convergence is based on a new study of the asymptotic behavior of the solution of the semidiscrete problems, which are obtained after the time discretization by the Euler method. Some numerical results are showed corroborating in practice the theoretical results on the uniform convergence and the order of the method.     

On nonlinear reaction-diffusion systems

This paper presents a qualitative analysis for a coupled system of two reaction-diffusion equations under various boundary conditions which arises from a number of physical problems. The nonlinear reaction functions are classified into three basic types according to their relative quasi-monotone property. For each type of reaction functions, an existence-comparison theorem, in terms of upper and lower solutions, is established for the time-dependent system as well as some boundary value problems. Three concrete physical systems arising from epidemics, biochemistry and engineering are taken as representatives of the basic types of reacting problems. Through suitable construction of upper and lower solutions, various qualitative properties of the solution for each system are obtained. These include the existence and bounds of time-dependent solutions, asymptotic behavior of the solution, stability and instability of nontrivial steady-state solutions, estimates of stability regions, and finally the blowing-up property of the solution. Special attention is given to the homogeneous Neumann boundary condition.     

Resilience in reaction-diffusion systems

Reaction-diffusion systems with zero-flux Neumann boundariesare widely used to model various kinds of interaction in, forexample, the scientific fields of ecology, biology, chemistry,medicine and industry. The physical systems within these fieldsare often known to be (conditionally or unconditionally) resilientwith respect to shocks, disturbances or catastrophies in theimmediate environment. In order to be good mathematical modelsof such situations the reaction-diffusion systems must havethe same resilient or asymptotic behaviour as that of the physicalsituation. Three fundamentally different kinds of reaction termsare usually distinguished according to the entry signs of thereaction Jacobian: mutualism, mixed (predator-prey) interactionand competition. The asymptotic stability (in the Poincarésense) of mutualistic systems has already been studied extensively,but the results cannot be generalized (globally) to the othertwo fundamental types, which are not order-preserving. A partial(local) generalization is, however given here for these twotypes, involving simple Jacobian inequalities and knowledge(often prompted by the underlying physical situation) of invariantsets in solution space. The return time of resilient systemsand the approach rate of asymptotically stable solutions arealso estimated.     

Functions with nonzero finite difference

Vilnius Civil Engineering Institute. Translated from Litovskii Matematicheskii Sbornik (Lietuvos Matematikos Rinkinys), Vol. 30, No. 2, pp. 275–287, April–June, 1990.     

Conservative finite difference scheme for a singularly perturbed elliptic reaction-diffusion equation: Approximation of solutions and derivatives

A boundary value problem for a singularly perturbed elliptic reaction-diffusion equation in a vertical strip is considered. The derivatives are written in divergent form. The derivatives in the differential equation are multiplied by a perturbation parameter ɛ², where ɛ takes arbitrary values in the interval (0, 1]. As ɛ → 0, a boundary layer appears in the solution of this problem. Using the integrointerpolational method and the condensing grid technique, conservative finite difference schemes on flux grids are constructed that converge ɛ-uniformly at a rate of O(N ₁⁻²ln² N ₁ + N ₂⁻²), where N ₁ + 1 and N ₂ + 1 are the number of mesh points on the x ₁-axis and the minimal number of mesh points on a unit interval of the x ₂-axis respectively. The normalized difference derivatives ɛ ^k (∂ ^k /∂x ₁ ^k )u(x) (k = 1, 2), which are ɛ-uniformly bounded and approximate the normalized derivatives in the direction across the boundary layer, and the derivatives along the boundary layer (∂ ^k /∂ x ₂ ^k )u(x) (k = 1, 2) converge ɛ-uniformly at the same rate.     

Baxter’s difference systems and orthogonal rational functions

Rational functions orthogonal on the unit circle are considered beginning with their recurrence formulas. Various summability conditions are imposed on the recurrence coefficients and the asymptotics of the solutions are studied and the orthogonality measure is recovered. The techniques developed by Baxter and Benzaid and Lutz are used.     

Neural networks and reaction-diffusion systems with prescribed dynamics

Reaction-diffusion systems and neural networks are considered. We prove that they can produce any structurally stable inertial dynamics.