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 weak solution approach to a reaction–diffusion system modeling pattern formation on seashells

We investigate a reaction–diffusion system proposed by H. Meinhardt as a model for pattern formation on seashells. We give a new proof for the existence of a local weak solution for general initial conditions and parameters upon using an iterative approach. Furthermore, the solution is shown to exist globally for suitable initial data. The behavior of the solution in time and space is illustrated through numerical simulations. Copyright © 2009 John Wiley & Sons, Ltd.     

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     

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.     

Global existence to a reaction–diffusion equation through a self‐similar‐like upper solution

A demonstration method is presented, which will ensure the existence of positive global solutions in time to the reaction–diffusion equation ?u_t+Δu+u^p=0 in ?ⁿ×[0, ∞), for exponents p?3 and space dimensions n?3. This method does not require the initial value to have a specific uniform smallness condition, but rather to satisfy a bell‐like form. The method is based on a specific upper solution, which models the diffusion process of the heat equation. The upper solution is not self‐similar, but does have a self‐similar‐like form. After transforming the reaction–diffusion problem into an equivalent one, whose initial value is uniformly very small, a local solution is obtained in the time interval [0, 1] by the use of this upper solution. This local solution is then extended to [0, ∞) through an infinite sequence of extensions. At each step, an appropriate change of variables will transform the extension into a problem nearly identical to the local problem in [0, 1]. These transformations exploit the diffusive and self‐similar‐like nature of the upper solution. Copyright © 2009 John Wiley & Sons, Ltd.     

Turing instability and Hopf bifurcation in a diffusive Leslie–Gower predator–prey model

In this paper, a reaction‐diffusion predator–prey system that incorporates the Holling‐type II and a modified Leslie‐Gower functional responses is considered. For ODE, the local stability of the positive equilibrium is investigated and the specific conditions are obtained. For partial differential equation, we consider the dissipation and persistence of solutions, the Turing instability of the equilibrium solutions, and the Hopf bifurcation. By calculating the normal form, we derive the formulae, which can determine the direction and the stability of Hopf bifurcation according to the original parameters of the system. We also use some numerical simulations to illustrate our theoretical results. Copyright © 2016 John Wiley & Sons, Ltd.     

A Lyapunov functional for a triangular reaction–diffusion system with nonlinearities of exponential growth

The aim of this work is to study the global existence of solutions to a triangular system of reaction–diffusion equations, which describes epidemiological or chemical situations. On the basis of the construction of a suitable Lyapunov functional, we show that for any initial data, classical global solutions exist even when the nonlinearities are of exponential growth. Copyright © 2012 John Wiley & Sons, Ltd.     

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.     

Stability and Hopf bifurcation of a delayed reaction–diffusion neural network

In this paper, a delayed reaction–diffusion neural network with Neumann boundary conditions is investigated. By analyzing the corresponding characteristic equations, the local stability of the trivial uniform steady state is discussed. The existence of Hopf bifurcation at the trivial steady state is established. Using the normal form theory and the center manifold reduction of partial function differential equations, explicit formulae are derived to determine the direction and stability of bifurcating periodic solutions. Numerical simulations are carried out to illustrate the main results. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Well‐posedness,positivity, and time asymptotics properties for a reaction–diffusion model of plankton communities,involving a rational nonlinearity with singularity

In this work, we consider a reaction–diffusion system, modeling the interaction between nutrients, phytoplankton, and zooplankton. Using a semigroup approach in , we prove global existence, uniqueness, and positivity of the solutions. The nonlinearity is handled by providing estimates in , allowing to deal with most of the functional responses that describe predator/prey interactions (Holling I, II, III, Ivlev) in ecology. The paper finally exhibits some time asymptotic properties of the solutions.     

Global existence of solutions to reaction–diffusion systems without conditions on the nonlinearities growth

It is well known that, for reaction–diffusion systems, if the nonlinearities grow faster than a polynomial, nothing seems to be known for instance. The purpose of this paper is to give sufficient conditions guaranteeing global existence, uniqueness and uniform boundedness of solutions for coupled reaction–diffusion equations without condition growth on the reactions terms f and g in case f + g ≠ 0. These systems possess many and various applications in physics as the diffusion of the Phosphorus in the Silicone or some models describing some nuclear reactions; there have also been other applications in chemistry and biology. Our techniques are based on the Lyapunov functional methods. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

Large‐time behavior of a two‐scale semilinear reaction–diffusion system for concrete sulfatation

We study the large‐time behavior of (weak) solutions to a two‐scale reaction–diffusion system coupled with a nonlinear ordinary differential equations modeling the partly dissipative corrosion of concrete (or cement)‐based materials with sulfates. We prove that as t → ∞ , the solution to the original two‐scale system converges to the corresponding two‐scale stationary system. To obtain the main result, we make use essentially of the theory of evolution equations governed by subdifferential operators of time‐dependent convex functions developed combined with a series of two‐scale energy‐like time‐independent estimates. Copyright © 2014 John Wiley & Sons, Ltd.     

Higher-order compact finite difference method for systems of reaction–diffusion equations

This paper is concerned with a compact finite difference method for solving systems of two-dimensional reaction–diffusion equations. This method has the accuracy of fourth-order in both space and time. The existence and uniqueness of the finite difference solution are investigated by the method of upper and lower solutions, without any monotone requirement on the nonlinear term. Three monotone iterative algorithms are provided for solving the resulting discrete system efficiently, and the sequences of iterations converge monotonically to a unique solution of the system. A theoretical comparison result for the various monotone sequences is given. The convergence of the finite difference solution to the continuous solution is proved, and Richardson extrapolation is used to achieve fourth-order accuracy in time. An application is given to an enzyme–substrate reaction–diffusion problem, and some numerical results are presented to demonstrate the high efficiency and advantages of this new approach.     

A variable nonlinear splitting algorithm for reaction diffusion systems with self‐ and cross‐ diffusion

Self‐ and cross‐diffusion are important nonlinear spatial derivative terms that are included into biological models of predator–prey interactions. Self‐diffusion models overcrowding effects, while cross‐diffusion incorporates the response of one species in light of the concentration of another. In this paper, a novel nonlinear operator splitting method is presented that directly incorporates both self‐ and cross‐diffusion into a computational efficient design. The numerical analysis guarantees the accuracy and demonstrates appropriate criteria for stability. Numerical experiments display its efficiency and accuracy.     

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     

Transcendental Bernstein series for solving reaction–diffusion equations with nonlocal boundary conditions through the optimization technique

In this paper, we apply transcendental Bernstein series (TBS) for solving reaction–diffusion equations with nonlocal boundary conditions which is the novel approximation tool. To carry out the method, we firstly expand the solution of the system in the term of TBS through the operational matrix scheme. To determine the unknown free coefficients and control parameters appeared in TBS expansion, we define an optimization problem which combines the reaction–diffusion equation with its nonlocal boundary conditions. Then we use the Lagrange multipliers technique for converting the problem under study into a system of algebraic equations. High accuracy and simplicity in reducing the integral boundary conditions are some privileges of the proposed scheme. We emphasize that Bernstein polynomials is the particular case of transcendental Bernstein series. Theoretical discussion about convergence confirms the reliability of the proposed method. Some test problems are chosen to investigate the applicability and computational efficiency. The experimental results confirm that the obtained results are in good agreement with the exact solutions with high rate of convergence.     

Reaction–diffusion systems for the microscopic cellular model of the cardiac electric field

The paper deals with a mathematical model for the electric activity of the heart at microscopic level. The membrane model used to describe the ionic currents is a generalization of the phase‐I Luo–Rudy, a model widely used in 2‐D and 3‐D simulations of the action potential propagation. From the mathematical viewpoint the model is made up of a parabolic reaction diffusion system coupled with an ODE system. We derive existence and some regularity results. Copyright © 2006 John Wiley & Sons, Ltd.     

Quenching for a reaction–diffusion system with logarithmic singularity

In this paper, we study the quenching phenomenon for a reaction–diffusion system with singular logarithmic source terms and positive Dirichlet boundary conditions. Some sufficient conditions for quenching of the solutions in finite time are obtained, and the blow-up of time-derivatives at the quenching point is verified. Furthermore, under appropriate hypotheses, the non-simultaneous quenching of the system is proved, and the estimates of quenching rate is given.