Uniform grid approximation of nonsmooth solutions to the mixed boundary value problem for a singularly perturbed reaction-diffusion equation in a rectangle

A mixed boundary value problem for a singularly perturbed reaction-diffusion equation in a square is considered. A Neumann condition is specified on one side of the square, and a Dirichlet condition is set on the other three. It is assumed that the coefficient of the equation, its right-hand side, and the boundary values of the desired solution or its normal derivative on the sides of the square are smooth enough to ensure the required smoothness of the solution in a closed domain outside the neighborhoods of the corner points. No compatibility conditions are assumed to hold at the corner points. Under these assumptions, the desired solution in the entire closed domain is of limited smoothness: it belongs only to the Hölder class C ^μ, where μ ∈ (0, 1) is arbitrary. In the domain, a nonuniform rectangular mesh is introduced that is refined in the boundary domain and depends on a small parameter. The numerical solution to the problem is based on the classical five-point approximation of the equation and a four-point approximation of the Neumann boundary condition. A mesh refinement rule is described under which the approximate solution converges to the exact one uniformly with respect to the small parameter in the L _∞ ^h norm. The convergence rate is O(N ^?2ln² N), where N is the number of mesh nodes in each coordinate direction. The parameter-uniform convergence of difference schemes for mixed problems without compatibility conditions at corner points was not previously analyzed.     

On the Structure of Spectra of Modulated Travelling Waves

Modulated travelling waves are solutions to reaction-diffusion equations that are time-periodic in an appropriate moving coordinate frame. They may arise through Hopf bifurcations or essential instabilities from pulses o fronts. In this article, a framework for the stability analysis of such solutions is presented: the reaction-diffusion equation is cast as an ill-posed elliptic dynamical system in the spatial variable acting upon time-periodic functions. Using this formulation, points in the esolvent set, the point spectrum, and the essential spectrum of the linearization about a modulated travelling wave are related to the existence of exponential dichotomies on appopriate intervals for the associated spatial elliptic eigenvalue problem. Fredholm properties of the linearized operator are characterized by a relative Morse-Floe index of the elliptic equation. These results are proved without assumptions on the asymptotic shape of the wave. Analogous results are true for the spectra of travelling waves to parabolic equations on unbounded cylinders. As an application, we study the existence and stability of modulated spatially-periodic patterns with long-wavelength that accompany modulated pulses.     

Synchronization of the Glycolysis Reaction-Diffusion Model via Linear Control Law

The Selkov system, which is typically employed to model glycolysis phenomena, unveils some rich dynamics and some other complex formations in biochemical reactions. In the present work, the synchronization problem of the glycolysis reaction-diffusion model is handled and examined. In addition, a novel convenient control law is designed in a linear form and, on the other hand, the stability of the associated error system is demonstrated through utilizing a suitable Lyapunov function. To illustrate the applicability of the proposed schemes, several numerical simulations are performed in one- and two-spatial dimensions.     

The Richardson scheme for the singularly perturbed parabolic reaction-diffusion equation in the case of a discontinuous initial condition

The Dirichlet problem for a singularly perturbed parabolic reaction-diffusion equation with a piecewise continuous initial condition in a rectangular domain is considered. The higher order derivative in the equation is multiplied by a parameter ?², where ? ∈ (0, 1]. When ? is small, a boundary and an interior layer (with the characteristic width ?) appear, respectively, in a neighborhood of the lateral part of the boundary and in a neighborhood of the characteristic of the reduced equation passing through the discontinuity point of the initial function; for fixed ?, these layers have limited smoothness. Using the method of additive splitting of singularities (induced by the discontinuities of the initial function and its low-order derivatives) and the condensing grid method (piecewise uniform grids that condense in a neighborhood of the boundary layers), a finite difference scheme is constructed that converges ?-uniformly at a rate of O(N ^?2ln² N + n ₀ ^?1 ), where N + 1 and N ₀ + 1 are the numbers of the mesh points in x and t, respectively. Based on the Richardson technique, a scheme that converges ?-uniformly at a rate of O(N ^?3 + N ₀ ^?2 ) is constructed. It is proved that the Richardson technique cannot construct a scheme that converges in ?-uniformly in x with an order greater than three.     

非线性反应扩散方程的广义条件对称及其精确解

利用广义条件对称,考虑非线性反应扩散方程的精确解,对应于不同的参数讨论,得到相应的方程及其允许的广义条件对称,进而得到方程的精确解.     

Effect of a Local Source or Sink of Inhibitor on Turing Patterns

We investigate a reaction-diffusion model in which a Turing pattern develops and reproduces the formation of periodic segments behind a propagating chemical wave front. The chemical scheme involves two species known as activator and inhibitor. The model can be used to mimic the formation of prevertebrae during the early development of vertebrate embryo. Deterministic and stochastic analyses of the reaction-diffusion processes are performed for two typical sets of parameter values, far from and close to the Turing bifurcation. The effects of a local source or sink of inhibitor on the growing structure are studied and successfully compared with experiments performed on chick embryos. We show that fluctuations may lead to the formation of additional prevertebra.     

Computational analysis of the roles of biochemical reactions in anomalous diffusion dynamics

Most biochemical processes in cells are usually modeled by reaction–diffusion(RD) equations. In these RD models,the diffusive process is assumed to be Gaussian. However, a growing number of studies have noted that intracellular diffusion is anomalous at some or all times, which may result from a crowded environment and chemical kinetics. This work aims to computationally study the effects of chemical reactions on the diffusive dynamics of RD systems by using both stochastic and deterministic algorithms. Numerical method to estimate the mean-square displacement(MSD) from a deterministic algorithm is also investigated. Our computational results show that anomalous diffusion can be solely due to chemical reactions. The chemical reactions alone can cause anomalous sub-diffusion in the RD system at some or all times.The time-dependent anomalous diffusion exponent is found to depend on many parameters, including chemical reaction rates, reaction orders, and chemical concentrations.     

Symmetry Analysis and Conservation Laws to the(2+1)-Dimensional Coupled Nonlinear Extension of the Reaction-Diffusion Equation

In this paper, a detailed Lie symmetry analysis of the(2+1)-dimensional coupled nonlinear extension of the reaction-diffusion equation is presented. The general finite transformation group is derived via a simple direct method,which is equivalent to Lie point symmetry group actually. Similarity reduction and some exact solutions of the original equation are obtained based on the optimal system of one-dimensional subalgebras. In addition, conservation laws are constructed by employing the new conservation theorem.     

Asymptotic solution of nonlocal nonlinear reaction-diffusion Robin problems with two parameters

In this paper, the nonlocal nonlinear reaction-diffusion singularly perturbed problems with two parameters are studied. Using a singular perturbation method, the structure of the solutions to the problem is discussed in relation to two small parameters. The asymptotic solutions of the problem are given.     

Bifurcation in infinite dimensional spaces and applications in spatiotemporal biological and chemical models

Recent advances in abstract local and global bifurcation theory is briefly reviewed. Several applications are included to illustrate the applications of abstract theory, and it includes Turing instability of chemical reactions, pattern formation in water limited ecosystems, and diffusive predator-prey models.