Lattice Boltzmann model for the bimolecular autocatalytic reaction–diffusion equation

A lattice Boltzmann model for the bimolecular autocatalytic reaction–diffusion equation is proposed. By using multi-scale technique and the Chapman–Enskog expansion on complex lattice Boltzmann equation, we obtain a series of complex partial differential equations, complex equilibrium distribution function and its complex moments. Then, the complex reaction–diffusion equation is recovered with higher-order accuracy of the truncation error. This equation can be used to describe the bimolecular autocatalytic reaction–diffusion systems, in which a rich variety of behaviors have been observed. Based on this model, the Fitzhugh–Nagumo model and the Gray–Scott model are simulated. The comparisons between the LBM results and the Alternative Direction Implicit results are given in detail. The numerical examples show that assumptions of source term can be used to raise the accuracy of the truncation error of the lattice Boltzmann scheme for the complex reaction–diffusion equation.     

Turing pattern formation for reaction–convection–diffusion systems in fixed domains submitted to toroidal velocity fields

We have studied the effect of advection on reaction–diffusion equations by using toroidal velocity fields. Turing patterns formation in diffusion–advection–reaction problems was studied specifically, considering the Schnackenberg and glycolysis reaction kinetics models. Four cases were analyzed and solved numerically using finite elements. For glycolysis models, the advective effect modified the form of Turing patterns obtained with diffusion–reaction; whereas for Schnackenberg problems, the original patterns distorted themselves slightly, making them rotate in direction of the velocity field. We have also determined that the advective effect surpassed the diffusive one for high values of velocity and instability driven by diffusion was eliminated. On the other hand the advective effect is not considerable for very low values in the velocity field, and there was no modification in the original Turing pattern.     

Computational examples of reaction–convection–diffusion equations solution under the influence of fluid flow: First example

This paper presents several numerical tests on reaction–diffusion equations in the Turing space, affected by convective fields present in incompressible flows under the Schnakenberg reaction mechanism. The tests are performed in 2D on square unit, to which we impose an advective field from the solution of the problem of the flow in a cavity. The model developed consists of a decoupled system of equations of reaction–advection–diffusion, along with the Navier–Stokes equations of incompressible flow, which is solved simultaneously using the finite element method. The results show that the pattern generated by the concentrations of the reacting system varies both in time and space due to the effect exerted by the advective field.     

Reaction–diffusion in viscoelastic materials

In this paper we study initial boundary value problems that describe reaction–diffusion phenomena in viscoelastic materials. The mathematical model, represented by an integro-differential equation coupled with an ordinary differential equation, is analyzed from theoretical and numerical viewpoints.     

On the derivation of fractional diffusion equation with an absorbent term and a linear external force

Recently, the generalized fractional reaction–diffusion equation subject to an external linear force field has been proposed to describe the transport processes in disordered systems. The solution of this generalized model can be formally expressed in closed form through the Fox function. For the sack of completeness, we dedicate this work to construct a neatly derivation of the generalized fractional reaction–diffusion equation. Remarkably, such derivation could in general offer some novel and inspiring inspection to the phenomena of anomalous transport. For instance, there is a strong evidence that the fractional calculus offers some physical insight into the origin of fractional dynamics for a systems which exhibit multiple trapping.     

On the uniqueness of semi-wavefronts for non-local delayed reaction–diffusion equations

We establish the uniqueness of semi-wavefront solution for a non-local delayed reaction–diffusion equation. This result is obtained by using a generalization of the Diekmann–Kaper theory for a nonlinear convolution equation. Several applications to the systems of non-local reaction–diffusion equations with distributed time delay are also considered.     

IMEX Runge–Kutta schemes for reaction–diffusion equations

A fundamental research is carried out into convergence and stability properties of IMEX (implicit–explicit) Runge–Kutta schemes applied to reaction–diffusion equations. It is shown that a fully discrete scheme converges if it satisfies certain conditions using a technique of the B-convergence analysis, developed by Burrage, Hundsdorfer and Verwer in 1986. Stability of the schemes is also examined on the basis of a scalar test equation, proposed by Frank, Hundsdorfer and Verwer in 1997.     

A convergence result for the periodic unfolding method related to fast diffusion on manifolds

Based on the periodic unfolding method in periodic homogenization, we deduce a convergence result for gradients of functions defined on connected, smooth, and periodic manifolds. Under the assumption of certain a-priori estimates of the gradient, which are typical for fast diffusion, the sum of a term involving a gradient with respect to the slow variable and one with respect to the fast variable is obtained in the homogenization limit. In addition, we show in a brief example how to apply this result and find for a reaction–diffusion equation defined on a periodic manifold that the homogenized equation contains a term describing macroscopic diffusion.     

Existence of periodic travelling wave solutions of non-autonomous reaction–diffusion equations with lambda–omega type

Periodic travelling wave solutions of reaction–diffusion equations were studied by many authors. The λ–ω$\lambda \text{–}\omega$ type reaction–diffusion system is a notable special model that admits explicit periodic travelling wave solutions and was introduced by Kopell and Howard in 1973. There are now similar systems which are investigated by means of autonomous dynamics. In contrast, there are few papers which are concerned with non-autonomous cases. For this reason, we apply Mawhin’s continuation theorem to derive the existence of periodic travelling wave solutions for non-autonomous λ–ω$\lambda \text{–}\omega$ systems, and we describe the ‘disappearance’ of periodic travelling wave solutions under special situations. Our main result is also illustrated by examples.     

Stability analysis of impulsive Cohen–Grossberg neural networks with distributed delays and reaction–diffusion terms

In this paper, we investigate a class of impulsive Cohen–Grossberg neural networks with distributed delays and reaction–diffusion terms. By establishing an integro-differential inequality with impulsive initial conditions and applying M-matrix theory, we find some sufficient conditions ensuring the existence, uniqueness, global exponential stability and global robust exponential stability of equilibrium point for impulsive Cohen–Grossberg neural networks with distributed delays and reaction–diffusion terms. An example is given to illustrate the results obtained here.     

Wick-stochastic finite element solution of reaction–diffusion problems

Real life reaction–diffusion problems are characterized by their inherent or externally induced uncertainties in the design parameters. This paper presents a finite element solution of reaction–diffusion equations of Wick type. Using the Wick-product properties and the Wiener–Itô chaos expansion, the stochastic variational problem is reformulated to a set of deterministic variational problems. To obtain the chaos coefficients in the corresponding deterministic reaction–diffusion, we implement the usual Galerkin finite element method using standard techniques. Once this representation is computed, the statistics of the numerical solution can be easily evaluated. Computational results are shown for one- and two-dimensional test examples.     

Existence and uniqueness of solutions to nonlinear evolution equations with locally monotone operators

In this paper we establish the existence and uniqueness of solutions for nonlinear evolution equations on a Banach space with locally monotone operators, which is a generalization of the classical result for monotone operators. In particular, we show that local monotonicity implies pseudo-monotonicity. The main results are applied to PDE of various types such as porous medium equations, reaction–diffusion equations, the generalized Burgers equation, the Navier–Stokes equation, the 3D Leray-α model and the p-Laplace equation with non-monotone perturbations.     

Large time behavior for some nonlinear degenerate parabolic equations

We study the asymptotic behavior of Lipschitz continuous solutions of nonlinear degenerate parabolic equations in the periodic setting. Our results apply to a large class of Hamilton–Jacobi–Bellman equations. Defining Σ as the set where the diffusion vanishes, i.e., where the equation is totally degenerate, we obtain the convergence when the equation is uniformly parabolic outside Σ and, on Σ, the Hamiltonian is either strictly convex or satisfies an assumption similar of the one introduced by Barles–Souganidis (2000) for first-order Hamilton–Jacobi equations. This latter assumption allows to deal with equations with nonconvex Hamiltonians. We can also release the uniform parabolic requirement outside Σ. As a consequence, we prove the convergence of some everywhere degenerate second-order equations.     

On the existence of non-monotone non-oscillating wavefronts

We present a monostable delayed reaction–diffusion equation with the unimodal birth function which admits only non-monotone wavefronts. Moreover, these fronts are either eventually monotone (in particular, such is the minimal wave) or slowly oscillating. Hence, for the Mackey–Glass type diffusive equations, we answer affirmatively the question about the existence of non-monotone non-oscillating wavefronts. As it was recently established by Hasik et al. and Ducrot et al., the same question has a negative answer for the KPP-Fisher equation with a single delay.     

An analytical solution for diffusion and nonlinear uptake of oxygen in a spherical cell

We develop a new analytical solution for a reactive transport model that describes the steady-state distribution of oxygen subject to diffusive transport and nonlinear uptake in a sphere. This model was originally reported by Lin [S.H. Lin, Oxygen diffusion in a spherical cell with nonlinear oxygen uptake kinetics, J. Theor. Biol. 60 (1976) 449–457] to represent the distribution of oxygen inside a cell and has since been studied extensively by both the numerical analysis and formal analysis communities. Here we extend these previous studies by deriving an analytical solution to a generalized reaction–diffusion equation that encompasses Lin’s model as a particular case. We evaluate the solution for the parameter combinations presented by Lin and show that the new solutions are identical to a grid-independent numerical approximation.     

Application of He's variational iteration method for solving the Cauchy reaction–diffusion problem

In this paper, the solution of Cauchy reaction–diffusion problem is presented by means of variational iteration method. Reaction–diffusion equations have special importance in engineering and sciences and constitute a good model for many systems in various fields. Application of variational iteration technique to this problem shows the rapid convergence of the sequence constructed by this method to the exact solution. Moreover, this technique does not require any discretization, linearization or small perturbations and therefore it reduces significantly the numerical computations.     

Rise and Fall of Periodic Patterns for a Generalized Klausmeier–Gray–Scott Model

In this paper we introduce a conceptual model for vegetation patterns that generalizes the Klausmeier model for semi-arid ecosystems on a sloped terrain (Klausmeier in Science 284:1826–1828, 1999). Our model not only incorporates downhill flow, but also linear or nonlinear diffusion for the water component. To relate the model to observations and simulations in ecology, we first consider the onset of pattern formation through a Turing or a Turing–Hopf bifurcation. We perform a Ginzburg–Landau analysis to study the weakly nonlinear evolution of small amplitude patterns and we show that the Turing/Turing–Hopf bifurcation is supercritical under realistic circumstances. In the second part we numerically construct Busse balloons to further follow the family of stable spatially periodic (vegetation) patterns. We find that destabilization (and thus desertification) can be caused by three different mechanisms: fold, Hopf and sideband instability, and show that the Hopf instability can no longer occur when the gradient of the domain is above a certain threshold. We encounter a number of intriguing phenomena, such as a ‘Hopf dance’ and a fine structure of sideband instabilities. Finally, we conclude that there exists no decisive qualitative difference between the Busse balloons for the model with standard diffusion and the Busse balloons for the model with nonlinear diffusion.     

Viscosity solutions of a class of degenerate quasilinear parabolic equations with Dirichlet boundary condition

This paper is concerned with viscosity solutions for a class of degenerate quasilinear parabolic equations in a bounded domain with homogeneous Dirichlet boundary condition. The equation under consideration arises from a number of practical model problems including reaction–diffusion processes in a porous medium. The degeneracy of the problem appears on the boundary and possibly in the interior of the domain. The goal of this paper is to establish some comparison properties between viscosity upper and lower solutions and to show the existence of a continuous viscosity solution between them. An application of the above results is given to a porous-medium type of reaction–diffusion model which demonstrates some distinctive properties of the solution when compared with the corresponding semilinear problem.     

Log-Harnack inequality for mild solutions of SPDEs with multiplicative noise

Due to technical reasons, existing results concerning Harnack type inequalities for SPDEs with multiplicative noise apply only to the case where the coefficient in the noise term is a Hilbert–Schmidt perturbation of a constant bounded operator. In this paper we obtained gradient estimates, log-Harnack inequality for mild solutions of general SPDEs with multiplicative noise whose coefficient is even allowed to be unbounded which cannot be Hilbert–Schmidt. Applications to stochastic reaction–diffusion equations driven by space–time white noise are presented.