Direct discontinuous Galerkin method for nonlinear reaction–diffusion systems in pattern formation

Nonlinear reaction–diffusion systems are often employed in mathematical modeling for pattern formation. Most of the work to date has been concerned within one-dimensional or rectangular domains. However, it is recognised that in most applications multidimensional complex geometrical domains are typically more important. In this paper we solve reaction–diffusion systems by combining direct discontinuous Galerkin (DDG) finite element methods with implicit integration factor (IIF) time integration method, on triangular meshes. This allows us solve the nonlinear algebraic systems on an element-by-element bases with significant gains in computational time. Numerical solutions of two reaction–diffusion systems, the well-studied Schnakenberg model and chloride–iodide–malonic acid (CIMA) reactive model, are presented to demonstrate effects of various domain geometries on the resulting biological patterns. Our numerical results are in good agreement with other numerical and analytical results, and with experimental results.     

Turing pattern formation in a predator–prey system with cross diffusion

The paper explores the impacts of cross-diffusion on the formation of spatial patterns in a ratio-dependent predator–prey system with zero-flux boundary conditions. Our results show that under certain conditions, cross-diffusion can trigger the emergence of spatial patterns which is however impossible under the same conditions when cross-diffusion is absent. We give a rigorous proof that the model has at least one spatially heterogenous steady state by means of the Leray–Schauder degree theory. In addition, numerical simulations are performed to visualize the complex spatial patterns.     

A waiting time phenomenon for modulations of pattern in reaction–diffusion systems

In order to describe slow modulations in time and space of stable or slightly unstable spatially periodic stationary solutions of pattern forming reaction–diffusion systems, so-called phase diffusion equations and Cahn–Hilliard equations can be derived via multiple scaling analysis as formal approximation equations. In the case that these equations degenerate, waiting time phenomena are well known to occur. In this paper, we prove that such waiting time phenomena can also occur approximately in the original reaction–diffusion systems by proving estimates between the formal approximations and the exact solutions of the original systems.     

A self-organized mesh generator using pattern formation in a reaction–diffusion system

A new type of mesh generator is developed by using a self-organized pattern in a reaction–diffusion system. The system is the Gray–Scott model, which creates a spot pattern in a specific parameter region. The spots correspond to nodes of a mesh. The mesh generator has several advantages: the algorithm is simple and processes to improve the mesh, such as smoothing, (locally) addition, and removal of nodes, are automatically performed by the system.     

Turing pattern formation in a predator–prey–mutualist system

In this paper, we develop a theoretical framework about spatial patterns in a three-species predator–prey–mutualist system with cross-diffusion. We concentrate on three aspects of Turing pattern formation: (1) what conditions enable the occurrence of Turing patterns? (2) what are the underlying mechanisms? (3) what are the corresponding configurations? For the first two questions, by use of the stability analysis for the positive uniform solution and the Leray–Schauder degree theory, we prove that under some conditions, the system admits at least a nonhomogeneous stationary solution. For the third question, we carry out numerical simulations for a Turing pattern, and we show that the configurations of Turing pattern are stable spotted patterns, which resemble a real ecosystem.     

Stabilization to a positive equilibrium for some reaction–diffusion systems

The aim of this paper is to study the asymptotic behavior of solutions for some reaction–diffusion systems in biology. First, we establish a Liouville type theorem for entire solutions of these reaction–diffusion systems. Based on this theorem, we derive the stabilization of the solutions of the reaction–diffusion system to the unique positive constant state, under the condition that this positive constant state is globally stable in the corresponding kinetic systems. Two specific examples about spreading phenomena from ecology and epidemiology are given to illustrate the application of this theory.     

Viability for nonlinear multi-valued reaction–diffusion systems

In this paper we consider a nonlinear evolution reaction–diffusion system governed by multi-valued perturbations of m-dissipative operators, generators of nonlinear semigroups of contractions. Let X and Y be real Banach spaces, ${\mathcal{K}}In this paper we consider a nonlinear evolution reaction–diffusion system governed by multi-valued perturbations of m-dissipative operators, generators of nonlinear semigroups of contractions. Let X and Y be real Banach spaces, K{\mathcal{K}} be a nonempty and locally closed subset in \mathbbR ×X×Y, A:D(A) í X\rightsquigarrow X, B:D(B) í Y\rightsquigarrow Y{\mathbb{R} \times X\times Y,\, A:D(A)\subseteq X\rightsquigarrow X, B:D(B)\subseteq Y\rightsquigarrow Y} two m-dissipative operators, F:K ? X{F:\mathcal{K} \rightarrow X} a continuous function and G:K \rightsquigarrow Y{G:\mathcal{K} \rightsquigarrow Y} a nonempty, convex and closed valued, strongly-weakly upper semi-continuous (u.s.c.) multi-function. We prove a necessary and a sufficient condition in order that for each (t,x,h) ? K{(\tau,\xi,\eta)\in \mathcal{K}}, the next system

{ lc u￠(t) ? Au(t)+F(t,u(t),v(t))    t 3 tv￠(t) ? Bv(t)+G(t,u(t),v(t))    t 3 tu(t)=x,    v(t)=h, \left\{ \begin{array}{lc} u'(t)\in Au(t)+F(t,u(t),v(t))\quad t\geq\tau \\ v'(t)\in Bv(t)+G(t,u(t),v(t))\quad t\geq\tau \\ u(\tau)=\xi,\quad v(\tau)=\eta, \end{array} \right.      8. Blow-up rates for semi-linear reaction–diffusion systems      Henghui Zou 《Journal of Differential Equations》2014       9. On the Dirichlet problem for reaction–diffusion equations in non-smooth domains      《Nonlinear Analysis: Theory, Methods & Applications》2001,47(2):765-776       10. Bifurcation to locked fronts in two component reaction–diffusion systems      Grégory Faye  Matt Holzer 《Annales de l'Institut Henri Poincaré (C) Analyse Non Linéaire》2019,36(2):545-584 We study invasion fronts and spreading speeds in two component reaction–diffusion systems. Using a variation of Lin's method, we construct traveling front solutions and show the existence of a bifurcation to locked fronts where both components invade at the same speed. Expansions of the wave speed as a function of the diffusion constant of one species are obtained. The bifurcation can be sub or super-critical depending on whether the locked fronts exist for parameter values above or below the bifurcation value. Interestingly, in the sub-critical case numerical simulations reveal that the spreading speed of the PDE system does not depend continuously on the coefficient of diffusion.      11. Dead cores in nonlinear reaction–diffusion systems      《Nonlinear Analysis: Theory, Methods & Applications》1987,11(10):1219-1227       12. Global existence and blowing-up of solutions to a class of coupled reaction–convection–diffusion systems      《Applied Mathematics Letters》2014 This paper concerns the global existence and blowing-up of solutions to the homogeneous Neumann problem of a coupled reaction–convection–diffusion system. The critical Fujita curve is determined and blowing-up theorem of Fujita type is established. An interesting phenomenon is that the critical Fujita curve even could be the infinite due to the convection.      13. Simultaneous identification of convection velocity and source strength in a convection–diffusion equation      Ting Cheng  Junjun Hu 《Applicable analysis》2020,99(12):2170-2189 In this work we are concerned with the analysis on a simultaneous finite element reconstruction of the convection velocity and source strength in a time-dependent convection–diffusion equation. The ill-posed problem is formulated into an output least-squares nonlinear minimization by an appropriately selected Tikhonov regularization. The regularizing effect and mathematical properties of the regularized system are justified and demonstrated. The nonlinear optimization problem is approximated by a fully discrete finite element method, whose convergence is rigorously established.      14. Global existence for degenerate quadratic reaction–diffusion systems      M. Pierre  R. Texier-Picard 《Annales de l'Institut Henri Poincaré (C) Analyse Non Linéaire》2009 We consider a class of degenerate reaction–diffusion systems with quadratic nonlinearity and diffusion only in the vertical direction. Such systems can appear in the modeling of photochemical generation and atmospheric dispersion of pollutants. The diffusion coefficients are different for all equations. We study global existence of solutions.      15. Instability induced by cross-diffusion in reaction–diffusion systems      Canrong Tian  Zhigui Lin  Michael Pedersen 《Nonlinear Analysis: Real World Applications》2010,11(2):1036-1045 In this paper the instability of the uniform equilibrium of a general strongly coupled reaction–diffusion is discussed. In unbounded domain and bounded domain the sufficient conditions for the instability are obtained respectively. The conclusion is applied to the ecosystem, it is shown that cross-diffusion can induce the instability of an equilibrium which is stable for the kinetic system and for the self-diffusion–reaction system.      16. Conservation laws for a variable coefficient nonlinear diffusion–convection–reaction equation      Zhijie Cao  Yiping Lin  Jianwei Shen 《Journal of Mathematical Analysis and Applications》2014 We find the Lie point symmetries of a class of second-order nonlinear diffusion–convection–reaction equations containing an unspecified coefficient function of the independent variable t and determine the subclasses of these equations which are nonlinearly self-adjoint. By using a general theorem on conservation laws proved recently by N.H. Ibragimov we establish conservation laws corresponding to the aforementioned Lie point symmetries, one by one, for the simultaneous system of the original equation together with its adjoint equation through a formal Lagrangian. Particularly, for the nonlinearly self-adjoint subclasses, we construct conservation laws for the corresponding equations themselves.      17. Generalized travelling waves for perturbed monotone reaction–diffusion systems      《Nonlinear Analysis: Theory, Methods & Applications》2001,46(6):757-776       18. Global dynamics of delayed reaction–diffusion equations in unbounded domains      Taishan Yi  Yuming Chen  Jianhong Wu 《Zeitschrift für Angewandte Mathematik und Physik (ZAMP)》2012,63(5):793-812 We consider a nonlocal delayed reaction–diffusion equation in an unbounded domain that includes some special cases arising from population dynamics. Due to the non-compactness of the spatial domain, the solution semiflow is not compact. We first show that, with respect to the compact open topology for the natural phase space, the solutions induce a compact and continuous semiflow ${\Phi}$ on a bounded and positively invariant set Y in C +?=?C([?1, 0], X +) that attracts every solution of the equation, where X + is the set of all bounded and uniformly continuous functions from ${\mathbb{R}}$ to [0, ∞). Then, to overcome the difficulty in describing the global dynamics, we establish a priori estimate for nontrivial solutions after describing the delicate asymptotic properties of the nonlocal delayed effect and the diffusion operator. The estimate enables us to show the permanence of the equation with respect to the compact open topology. With the help of the permanence, we can employ standard dynamical system theoretical arguments to establish the global attractivity of the nontrivial equilibrium. The main results are illustrated with the diffusive Nicholson’s blowfly equation and the diffusive Mackey–Glass equation.      19. A gradient structure for systems coupling reaction–diffusion effects in bulk and interfaces      Annegret Glitzky  Alexander Mielke 《Zeitschrift für Angewandte Mathematik und Physik (ZAMP)》2013,64(1):29-52 We derive gradient-flow formulations for systems describing drift-diffusion processes of a finite number of species which undergo mass-action type reversible reactions. Our investigations cover heterostructures, where material parameter may depend in a nonsmooth way on the space variable. The main results concern a gradient-flow formulation for electro-reaction–diffusion systems with active interfaces permitting drift-diffusion processes and reactions of species living on the interface and transfer mechanisms allowing bulk species to jump into an interface or to pass through interfaces. The gradient flows are formulated in terms of two functionals: the free energy and the dissipation potential. Both functionals consist of a bulk and an interface integral. The interface integrals determine the interface dynamics as well as the self-consistent coupling to the model in the bulk. The advantage of the gradient structure is that it automatically generates thermodynamically consistent models.      20. Time periodic traveling wave solutions for periodic advection–reaction–diffusion systems      Guangyu Zhao  Shigui Ruan 《Journal of Differential Equations》2014 We study the existence, uniqueness, and asymptotic stability of time periodic traveling wave solutions to a class of periodic advection–reaction–diffusion systems. Under certain conditions, we prove that there exists a maximal wave speed c?$c^{?}$ such that for each wave speed c≤c?$c\le c^{?}$, there is a time periodic traveling wave connecting two periodic solutions of the corresponding kinetic system. It is shown that such a traveling wave is unique modulo translation and is monotone with respect to its co-moving frame coordinate. We also show that the traveling wave solutions with wave speed c≤c?$c\le c^{?}$ are asymptotically stable in certain sense. In addition, we establish the nonexistence of time periodic traveling waves with speed c>c?$c>c^{?}$.

Bifurcation to locked fronts in two component reaction–diffusion systems

We study invasion fronts and spreading speeds in two component reaction–diffusion systems. Using a variation of Lin's method, we construct traveling front solutions and show the existence of a bifurcation to locked fronts where both components invade at the same speed. Expansions of the wave speed as a function of the diffusion constant of one species are obtained. The bifurcation can be sub or super-critical depending on whether the locked fronts exist for parameter values above or below the bifurcation value. Interestingly, in the sub-critical case numerical simulations reveal that the spreading speed of the PDE system does not depend continuously on the coefficient of diffusion.     

Global existence and blowing-up of solutions to a class of coupled reaction–convection–diffusion systems

This paper concerns the global existence and blowing-up of solutions to the homogeneous Neumann problem of a coupled reaction–convection–diffusion system. The critical Fujita curve is determined and blowing-up theorem of Fujita type is established. An interesting phenomenon is that the critical Fujita curve even could be the infinite due to the convection.     

Simultaneous identification of convection velocity and source strength in a convection–diffusion equation

In this work we are concerned with the analysis on a simultaneous finite element reconstruction of the convection velocity and source strength in a time-dependent convection–diffusion equation. The ill-posed problem is formulated into an output least-squares nonlinear minimization by an appropriately selected Tikhonov regularization. The regularizing effect and mathematical properties of the regularized system are justified and demonstrated. The nonlinear optimization problem is approximated by a fully discrete finite element method, whose convergence is rigorously established.     

Global existence for degenerate quadratic reaction–diffusion systems

We consider a class of degenerate reaction–diffusion systems with quadratic nonlinearity and diffusion only in the vertical direction. Such systems can appear in the modeling of photochemical generation and atmospheric dispersion of pollutants. The diffusion coefficients are different for all equations. We study global existence of solutions.     

Instability induced by cross-diffusion in reaction–diffusion systems

In this paper the instability of the uniform equilibrium of a general strongly coupled reaction–diffusion is discussed. In unbounded domain and bounded domain the sufficient conditions for the instability are obtained respectively. The conclusion is applied to the ecosystem, it is shown that cross-diffusion can induce the instability of an equilibrium which is stable for the kinetic system and for the self-diffusion–reaction system.     

Conservation laws for a variable coefficient nonlinear diffusion–convection–reaction equation

We find the Lie point symmetries of a class of second-order nonlinear diffusion–convection–reaction equations containing an unspecified coefficient function of the independent variable t and determine the subclasses of these equations which are nonlinearly self-adjoint. By using a general theorem on conservation laws proved recently by N.H. Ibragimov we establish conservation laws corresponding to the aforementioned Lie point symmetries, one by one, for the simultaneous system of the original equation together with its adjoint equation through a formal Lagrangian. Particularly, for the nonlinearly self-adjoint subclasses, we construct conservation laws for the corresponding equations themselves.     

Global dynamics of delayed reaction–diffusion equations in unbounded domains

We consider a nonlocal delayed reaction–diffusion equation in an unbounded domain that includes some special cases arising from population dynamics. Due to the non-compactness of the spatial domain, the solution semiflow is not compact. We first show that, with respect to the compact open topology for the natural phase space, the solutions induce a compact and continuous semiflow ${\Phi}$ on a bounded and positively invariant set Y in C ₊?=?C([?1, 0], X ₊) that attracts every solution of the equation, where X ₊ is the set of all bounded and uniformly continuous functions from ${\mathbb{R}}$ to [0, ∞). Then, to overcome the difficulty in describing the global dynamics, we establish a priori estimate for nontrivial solutions after describing the delicate asymptotic properties of the nonlocal delayed effect and the diffusion operator. The estimate enables us to show the permanence of the equation with respect to the compact open topology. With the help of the permanence, we can employ standard dynamical system theoretical arguments to establish the global attractivity of the nontrivial equilibrium. The main results are illustrated with the diffusive Nicholson’s blowfly equation and the diffusive Mackey–Glass equation.     

A gradient structure for systems coupling reaction–diffusion effects in bulk and interfaces

We derive gradient-flow formulations for systems describing drift-diffusion processes of a finite number of species which undergo mass-action type reversible reactions. Our investigations cover heterostructures, where material parameter may depend in a nonsmooth way on the space variable. The main results concern a gradient-flow formulation for electro-reaction–diffusion systems with active interfaces permitting drift-diffusion processes and reactions of species living on the interface and transfer mechanisms allowing bulk species to jump into an interface or to pass through interfaces. The gradient flows are formulated in terms of two functionals: the free energy and the dissipation potential. Both functionals consist of a bulk and an interface integral. The interface integrals determine the interface dynamics as well as the self-consistent coupling to the model in the bulk. The advantage of the gradient structure is that it automatically generates thermodynamically consistent models.     

Time periodic traveling wave solutions for periodic advection–reaction–diffusion systems

We study the existence, uniqueness, and asymptotic stability of time periodic traveling wave solutions to a class of periodic advection–reaction–diffusion systems. Under certain conditions, we prove that there exists a maximal wave speed c^?$c^{?}$ such that for each wave speed c≤c^?$c\le c^{?}$, there is a time periodic traveling wave connecting two periodic solutions of the corresponding kinetic system. It is shown that such a traveling wave is unique modulo translation and is monotone with respect to its co-moving frame coordinate. We also show that the traveling wave solutions with wave speed c≤c^?$c\le c^{?}$ are asymptotically stable in certain sense. In addition, we establish the nonexistence of time periodic traveling waves with speed c>c^?$c>c^{?}$.