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.     

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.     

Global boundedness of solutions to a reaction–diffusion system

We obtain in this paper the global boundedness of solutions to a Fujita‐type reaction–diffusion system. This global boundedness results from diffusion effect, homogeneous Dirichlet boundary value conditions and appropriate reactions. Copyright © 1999 John Wiley & Sons, Ltd.     

Positive solutions and dynamics of a finite difference reaction–diffusion system

We investigate the dynamics and methods of computation for some nonlinear finite difference systems that are the discretized equations of a time-dependent and a steady-state reaction–diffusion problem. The formulation of the discrete equations for the time-dependent problem is based on the implicit method for parabolic equations, and the computational algorithm is based on the method of monotone iterations using upper and lower solutions as the initial iterations. The monotone iterative method yields improved upper and lower bounds of the solution in each iteration, and the sequence of iterations converges monotonically to a solution for both the time-dependent and the steady-state problems. An important consequence of this method is that it leads to a bifurcation point that determines the dynamic behavior of the time-dependent problem in relation to the corresponding steady-state problem. This bifurcation point also determines whether the steady-state problem has one or two non-negative solutions, and is explicitly given in terms of the physical parameters of the system and the type of boundary conditions. Numerical results are presented for both the time-dependent and the steady-state problems under various boundary conditions, including a test problem with known analytical solution. These numerical results exhibit the predicted dynamic behavior of the time-dependent solution given by the theoretical analysis. Also discussed are the numerical stability of the computational algorithm and the convergence of the finite difference solution to the corresponding continuous solution of the reaction–diffusion problem. © 1993 John Wiley & Sons, Inc.     

Global asymptotic stability of BAM neural networks with distributed delays and reaction–diffusion terms

The global asymptotic stability of bi-directional associative memory neural networks with distributed delays and reaction–diffusion terms are studied by using the analysis technique and Lyapunov functional. A sufficient condition is proposed. Two numerical examples are given to show the correctness of our analysis.     

Algebraic multigrid preconditioners for the bidomain reaction–diffusion system

The so-called bidomain system is possibly the most complete model for the cardiac bioelectric activity. It consists of a reaction–diffusion system, modeling the intra, extracellular and transmembrane potentials, coupled through a nonlinear reaction term with a stiff system of ordinary differential equations describing the ionic currents through the cellular membrane. In this paper we address the problem of efficiently solving the large linear system arising in the finite element discretization of the bidomain model, when a semiimplicit method in time is employed. We analyze the use of structured algebraic multigrid preconditioners on two major formulations of the model, and report on our numerical experience under different discretization parameters and various discontinuity properties of the conductivity tensors. Our numerical results show that the less exercised formulation provides the best overall performance on a typical simulation of the myocardium excitation process.     

Blow‐up estimates for a quasi‐linear reaction–diffusion system

In this paper, some sufficient conditions under which the quasilinear elliptic system ‐div(∣?_u∣^p‐2?_u) = uv, ‐div(∣?_u∣^q‐2?_u) = uv in ?^N(N≥3) has no radially symmetric positive solution is derived. Then by using this non‐existence result, blow‐up estimates for a class of quasilinear reaction–diffusion systems u_t = div (∣?_u∣^p‐2?_u)+uv,vt = div(∣?_v∣^q‐2?_v) +uv with the homogeneous Dirichlet boundary value conditions are obtained. Copyright © 2003 John Wiley & Sons, Ltd.     

A reaction–diffusion problem with a reaction front

We state a 1D model with quasi-stationary gas flows approximation for a carbon reactivity test in the production of silicon. The mathematical problem we formulate is a non-linear boundary value problem for a third-order ordinary differential equation with non-linear boundary conditions, which are non-local in time. We prove existence and uniqueness of a classical solution and provide a numerical example. © 1998 B. G. Teubner Stuttgart–John Wiley & Sons Ltd.     

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.     

A parameter-uniform Schwarz method for a coupled system of reaction–diffusion equations

We consider an arbitrarily sized coupled system of one-dimensional reaction–diffusion problems that are singularly perturbed in nature. We describe an algorithm that uses a discrete Schwarz method on three overlapping subdomains, extending the method in [H. MacMullen, J.J.H. Miller, E. O’Riordan, G.I. Shishkin, A second-order parameter-uniform overlapping Schwarz method for reaction-diffusion problems with boundary layers, J. Comput. Appl. Math. 130 (2001) 231–244] to a coupled system. On each subdomain we use a standard finite difference operator on a uniform mesh. We prove that when appropriate subdomains are used the method produces ε-uniform results. Furthermore we improve upon the analysis of the above-mentioned reference to show that, for small ε, just one iteration is required to achieve the expected accuracy.     

Global a posteriori error estimates for convection–reaction–diffusion problems

In this note we propose a nonstandard technique for constructing global a posteriori error estimates for the stationary convection–reaction–diffusion equation. In order to estimate the approximation error in appropriate weighted energy norms, which measures the overall quality of the approximations, the underlying bilinear form is decomposed into several terms which can be directly computed or easily estimated from above using elementary tools of functional analysis. Several auxiliary parameters are introduced to construct such a splitting and tune the resulting upper error bound. It is demonstrated how these parameters can be chosen in some natural and convenient way for computations so that the weighted energy norm of the error is almost recovered, which shows that the estimates proposed are, in fact, quasi-sharp. The presented methodology is completely independent of numerical techniques used to compute approximate solutions. In particular, it is applicable to approximations which fail to satisfy the Galerkin orthogonality, e.g., due to an inconsistent stabilization, flux limiting, low-order quadrature rules, round-off and iteration errors etc. Moreover, the only constant that appears in the proposed error estimates is of global nature and comes from the Friedrichs–Poincaré inequality.     

Numerical exploration of a system of reaction–diffusion equations with internal and transient layers

A reaction pathway for a classical two-species reaction is considered with one reaction that is several orders of magnitudes faster than the other. To sustain the fast reaction, the transport and reaction effects must balance in such a way as to give an internal layer in space. For the steady-state problem, existing singular perturbation analysis rigorously proves the correct scaling of the internal layer. This work reports the results of exploratory numerical simulations that are designed to provide guidance for the analysis to be performed for the transient problem. The full model is comprised of a system of time-dependent reaction–diffusion equations coupled through the non-linear reaction terms with mixed Dirichlet and Neumann boundary conditions. In addition to internal layers in space, the time-dependent problem possesses an initial transient layer in time. To resolve both types of layers as accurately as possible, we design a finite element method with analytic evaluation of all integrals. This avoids all errors associated with the evaluation of the non-linearities and allows us to provide an analytic Jacobian matrix to the implicit time stepping method. The numerical results show that the method resolves the localized sharp gradients accurately and can predict the scaling of the internal layers for the time-dependent problem.     

Transient and asymptotic dynamics of a prey–predator system with diffusion

In this paper, we study a prey–predator system associated with the classical Lotka–Volterra nonlinearity. We show that the dynamics of the system are controlled by the ODE part. First, we show that the solution becomes spatially homogeneous and is subject to the ODE part as t → ∞ . Next, we take the shadow system to approximate the original system as D → ∞ . The asymptotics of the shadow system are also controlled by those of the ODE. The transient dynamics of the original system approaches to the dynamics of its ODE part with the initial mean as D → ∞ . Although the asymptotic dynamics of the original system are also controlled by the ODE, the time periods of these ODE solutions may be different. Concerning this property, we have that any δ > 0 admits D₀ > 0 such that if , the time period of the ODE, satisfies , then the solution to the original system with D ≥ D₀ cannot approach the stationary state. Copyright © 2012 John Wiley & Sons, Ltd.     

Pulsating traveling waves in the singular limit of a reaction–diffusion system in solid combustion

We consider a coupled system of parabolic/ODE equations describing solid combustion. For a given rescaling of the reaction term (the high activation energy limit), we show that the limit solution solves a free boundary problem which is to our knowledge new.In the time-increasing case, the limit coincides with the Stefan problem with spatially inhomogeneous coefficients. In general it is a parabolic equation with a memory term.In the first part of our paper we give a characterization of the limit problem in one space dimension. In the second part of the paper, we construct a family of pulsating traveling waves for the limit one phase Stefan problem with periodic coefficients. This corresponds to the assumption of periodic initial concentration of reactant.     

A family of numerical methods for diffusion and reaction–diffusion equations

A family of methods is developed for the numerical solution of second-order parabolic partial differential equations in one space dimension. The methods are second-, third-, or fourth-order accurate in time; five of them are seen to be L₀-stable in the sense of Gourlay and Morris, while the sixth is seen to be A₀-stable, The methods are tested on four problems from the literature, three diffusion problems and one reaction–diffusion problem.     

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.     

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.     

Multiple bifurcations generated by mode interactions in a reaction–diffusion problem

We study multiple bifurcations in a system of reaction–diffusion equations defined on a unit square with Robin boundary conditions. First we investigate linear stabilities of the system at the uniform steady state solution. Then we discuss how multiple bifurcations can be generated by mode interactions of the system, and how these multiple bifurcations can be preserved in the associated discrete system. A continuation-unsymmetric Lanczos algorithm is described to trace discrete solution curves. Numerical experiments on the Brusselator equations are reported.     

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.