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.     

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.     

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     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

Construction of Lyapunov functionals for delay differential equations in virology and epidemiology

In the present paper, we present a method for constructing a Lyapunov functional for some delay differential equations in virology and epidemiology. Here some delays are incorporated to the original ordinary differential equations, for which a Lyapunov function is already obtained. We present simple and clear explanation of our method using some models whose Lyapunov functionals are already obtained. Moreover, we present several new results for constructing Lyapunov functionals using our method.     

Global existence and blow-up to a degenerate reaction–diffusion system with nonlinear memory

In this paper, we consider a degenerate reaction–diffusion system coupled by nonlinear memory. Under appropriate hypotheses, we prove that the solution either exists globally or blows up in finite time. Furthermore, the blow-up rates are obtained.     

Nonclassical diffusion with memory

We consider the nonclassical diffusion equation with hereditary memory on a bounded three‐dimensional domain. Setting the problem in the past history framework, and with very general assumptions on the memory kernel κ, we prove that the related solution semigroup possesses a global attractor of optimal regularity. Copyright © 2014 John Wiley & Sons, Ltd.     

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.     

Pullback attractors for the norm-to-weak continuous process and application to the nonautonomous reaction–diffusion equations

In this paper, we introduce the concept of norm-to-weak continuous process in a Banach space, and obtain the existence of pullback attractors for this kind of process. Then we give a new method for proving the existence of the pullback attractors. As an application, we obtain the existence of pullback attractors for nonautonomous reaction–diffusion equation in with exponential growth of the external force.     

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     

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     

Amplitude equation for the stochastic reaction‐diffusion equations with random Neumann boundary conditions

In this paper, we consider a quite general class of reaction‐diffusion equations with cubic nonlinearities and with random Neumann boundary conditions. We derive rigorously amplitude equations, using the natural separation of time‐scales near a change of stability and investigate whether additive degenerate noise and random boundary conditions can lead to stabilization of the solution of the stochastic partial differential equation or not. The nonlinear heat equation (Ginzburg–Landau equation) is used to illustrate our result. Copyright © 2015 John Wiley & Sons, Ltd.