Hopf bifurcation of a diffusive Gause‐type predator‐prey model induced by time fractional‐order derivatives

Since population behaviors possess the characteristic of history memory, we, in this paper, introduce time fractional‐order derivatives into a diffusive Gause‐type predator‐prey model, which is time fractional‐order reaction‐diffusion equations and a generalized form of its corresponding first‐derivative model. For this kind of model, we prove the existence and uniqueness of a global positive solution by using the theory of evolution equations and the comparison principle of time fractional‐order partial differential equations. Besides, we obtain the stability and Hopf bifurcation of the Gause‐type predator‐prey model in the forms of the time fractional‐order ordinary equations and of the time fractional‐order reaction‐diffusion equations, respectively. Our results show that the stable region of the parameters in these 2 models can be enlarged by the time fractional‐order derivatives. Some numerical simulations are made to verify our results.     

Wave simulations of Gray‐Scott reaction‐diffusion system

In this work, we study the numerical simulation of the one‐dimensional reaction‐diffusion system known as the Gray‐Scott model. This model is responsible for the spatial pattern formation, which we often meet in nature as the result of some chemical reactions. We have used the trigonometric quartic B‐spline (T4B) functions for space discretization with the Crank‐Nicolson method for time integration to integrate the nonlinear reaction‐diffusion equation into a system of algebraic equations. The solutions of the Gray‐Scott model are presented with different wave simulations. Test problems are chosen from the literature to illustrate the stationary waves, pulse‐splitting waves, and self‐replicating waves.     

High‐order characteristic‐tracking strategy for simulation of a nonlinear advection–diffusion equation

In this study, new high‐order backward semi‐Lagrangian methods are developed to solve nonlinear advection–diffusion type problems, which are realized using high‐order characteristic‐tracking strategies. The proposed characteristic‐tracking strategies are second‐order L‐stable and third‐order L(α)‐stable methods, which are based on a classical implicit multistep method combined with a error‐correction method. We also use backward differentiation formulas and the fourth‐order finite‐difference scheme for diffusion problem discretization in the temporal and spatial domains, respectively. To demonstrate the adaptability and efficiency of these time‐discretization strategies, we apply these methods to nonlinear advection–diffusion type problems such as the viscous Burgers' equation. Through simulations, not only the temporal and spatial accuracies are numerically evaluated but also the proposed methods are shown to be superior to the compared existing characteristic‐tracking methods under the same rates of convergence in terms of accuracy and efficiency. Finally, we have shown that the proposed method well preserves the energy and mass when the viscosity coefficient becomes zero.     

Existence and non‐existence of steady states to a cross‐diffusion system arising in a Leslie predator–prey model

This paper is concerned with a cross‐diffusion system arising in a Leslie predator–prey population model in a bounded domain with no flux boundary condition. We investigate sufficient condition for the existence and the non‐existence of non‐constant positive solution. We obtain that if natural diffusion coefficient of predator is large enough and cross‐diffusion coefficients are fixed, then under some conditions there exists non‐constant positive solution. Furthermore, we show that if natural diffusion coefficients of predator and prey are both large enough, and cross‐diffusion coefficients are small enough, then there exists no non‐constant positive solution. Copyright © 2012 John Wiley & Sons, Ltd.     

Uniform superconvergence of a Galerkin finite element method on Shishkin‐type meshes

We consider a Galerkin finite element method that uses piecewise bilinears on a class of Shishkin‐type meshes for a model singularly perturbed convection‐diffusion problem on the unit square. The method is shown to be convergent, uniformly in the diffusion parameter ϵ, of almost second order in a discrete weighted energy norm. As a corollary, we derive global L₂‐norm error estimates and local L_∞‐norm estimates. Numerical experiments support our theoretical results. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16:426–440, 2000     

A high‐order Godunov method for one‐dimensional convection‐diffusion‐reaction problems

In this article we develop a high‐order Godunov method for one‐dimensional convection‐diffusion‐reaction problems where convection dominates diffusion. The heart of this method comes from incorporating the diffusion term via the slope of the linear representation (recovery) of the solution on each grid cell. The method is conservative and explicit. Therefore, it is efficient in computing time. For constant coefficient linear convection, diffusion, and Lipschitz‐type reaction, the properties of the total variation stability and monotonicity preservation are proved. An error estimation is derived. Computational examples are presented and compared with the exact solutions. © 2000 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 16: 495–512, 2000     

Nonexistence of coexisting steady‐state solutions of Dirichlet problem for a cross‐diffusion model

We study a special case of Shigesada–Kawasaki–Teramoto (SKT) model for two competing species with the Dirichlet boundary conditions. In our case, one of the species is not influenced by self‐diffusion or cross‐diffusion. We specify the explicit range of parameters by contradiction such that there are no coexisting steady‐state solutions to the model.     

Dynamics of virus infection models with density‐dependent diffusion and Beddington–DeAngelis functional response

In this work, we integrate both density‐dependent diffusion process and Beddington–DeAngelis functional response into virus infection models to consider their combined effects on viral infection and its control. We perform global analysis by constructing Lyapunov functions and prove that the system is well posed. We investigated the viral dynamics for scenarios of single‐strain and multi‐strain viruses and find that, for the multi‐strain model, if the basic reproduction number for all viral strains is greater than 1, then each strain persists in the host. Our investigation indicates that treating a patient using only a single type of therapy may cause competitive exclusion, which is disadvantageous to the patient's health. For patients infected with several viral strains, the combination of several therapies is a better choice. Copyright © 2017 John Wiley & Sons, Ltd.     

Taylor‐Galerkin B‐spline finite element method for the one‐dimensional advection‐diffusion equation

The advection‐diffusion equation has a long history as a benchmark for numerical methods. Taylor‐Galerkin methods are used together with the type of splines known as B‐splines to construct the approximation functions over the finite elements for the solution of time‐dependent advection‐diffusion problems. If advection dominates over diffusion, the numerical solution is difficult especially if boundary layers are to be resolved. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show the behavior of the method with emphasis on treatment of boundary conditions. Taylor‐Galerkin methods have been constructed by using both linear and quadratic B‐spline shape functions. Results shown by the method are found to be in good agreement with the exact solution. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

Finite‐difference schemes for transport‐dominated equations using special collocation nodes

We introduce finite‐difference schemes based on a special upwind‐type collocation grid, in order to obtain approximations of the solution of linear transport‐dominated advection‐diffusion problems. The method is well suited when the diffusion parameter is very small compared to the discretization parameter. A theory is developed and many numerical experiments are shown. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

An Eulerian‐Lagrangian substructuring domain decomposition method for unsteady‐state advection‐diffusion equations

We develop an Eulerian‐Lagrangian substructuring domain decomposition method for the solution of unsteady‐state advection‐diffusion transport equations. This method reduces to an Eulerian‐Lagrangian scheme within each subdomain and to a type of Dirichlet‐Neumann algorithm at subdomain interfaces. The method generates accurate and stable solutions that are free of artifacts even if large time‐steps are used in the simulation. Numerical experiments are presented to show the strong potential of the method. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17:565–583, 2001     

Travelling Wavefronts in Reaction‐Diffusion Equations with Convection Effects and Non‐Regular Terms

This paper deals with the appearance of monotone bounded travelling wave solutions for a parabolic reaction‐diffusion equation which frequently meets both in chemical and biological systems. In particular, we prove the existence of monotone front type solutions for any wave speed c ≥ c* and give an estimate for the threshold value c*. Our model takes into account both of a density dependent diffusion term and of a non‐linear convection effect. Moreover, we do not require the main non‐linearity g to be a regular C¹ function; in particular we are able to treat both the case when g′(0) = 0, giving rise to a degenerate equilibrium point in the phase plane, and the singular case when g′(0) = +∞. Our results generalize previous ones due to Aronson and Weinberger [Adv. Math. 30 (1978), pp. 33–76 ], Gibbs and Murray (see Murray [Mathematical Biology, Springer‐Verlag, Berlin, 1993 ]) and McCabe , Leach and Needham [SIAM J. Appl. Math. 59 (1998), pp. 870–899 ]. Finally, we obtain our conclusions by means of a comparison‐type technique which was introduced and developed in this framework in a recent paper by the same authors.     

An upwind finite volume scheme and its maximum‐principle‐preserving ADI splitting for unsteady‐state advection‐diffusion equations

We develop an upwind finite volume (UFV) scheme for unsteady‐state advection‐diffusion partial differential equations (PDEs) in multiple space dimensions. We apply an alternating direction implicit (ADI) splitting technique to accelerate the solution process of the numerical scheme. We investigate and analyze the reason why the conventional ADI splitting does not satisfy maximum principle in the context of advection‐diffusion PDEs. Based on the analysis, we propose a new ADI splitting of the upwind finite volume scheme, the alternating‐direction implicit, upwind finite volume (ADFV) scheme. We prove that both UFV and ADFV schemes satisfy maximum principle and are unconditionally stable. We also derive their error estimates. Numerical results are presented to observe the performance of these schemes. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 211–226, 2003     

Numerical treatment of two‐parameter singularly perturbed parabolic convection diffusion problems with non‐smooth data

In the present work, we consider a parabolic convection‐diffusion‐reaction problem where the diffusion and convection terms are multiplied by two small parameters, respectively. In addition, we assume that the convection coefficient and the source term of the partial differential equation have a jump discontinuity. The presence of perturbation parameters leads to the boundary and interior layers phenomena whose appropriate numerical approximation is the main goal of this paper. We have developed a uniform numerical method, which converges almost linearly in space and time on a piecewise uniform space adaptive Shishkin‐type mesh and uniform mesh in time. Error tables based on several examples show the convergence of the numerical solutions. In addition, several numerical simulations are presented to show the effectiveness of resolving layer behavior and their locations.     

Finite element and finite volume‐element simulation of pseudo‐ECGs and cardiac alternans

In this paper, we are interested in the spatio‐temporal dynamics of the transmembrane potential in paced isotropic and anisotropic cardiac tissues. In particular, we observe a specific precursor of cardiac arrhythmias that is the presence of alternans in the action potential duration. The underlying mathematical model consists of a reaction–diffusion system describing the propagation of the electric potential and the nonlinear interaction with ionic gating variables. Either conforming piecewise continuous finite elements or a finite volume‐element scheme are employed for the spatial discretization of all fields, whereas operator splitting strategies of first and second order are used for the time integration. We also describe an efficient mechanism to compute pseudo‐ECG signals, and we analyze restitution curves and alternans patterns for physiological and pathological cardiac rhythms. Copyright © 2014 John Wiley & Sons, Ltd.     

A non‐standard finite difference scheme for an advection‐diffusion‐reaction equation

In this note, a non‐standard finite difference (NSFD) scheme is proposed for an advection‐diffusion‐reaction equation with nonlinear reaction term. We first study the diffusion‐free case of this equation, that is, an advection‐reaction equation. Two exact finite difference schemes are constructed for the advection‐reaction equation by the method of characteristics. As these exact schemes are complicated and are not convenient to use, an NSFD scheme is derived from the exact scheme. Then, the NSFD scheme for the advection‐reaction equation is combined with a finite difference space‐approximation of the diffusion term to provide a NSFD scheme for the advection‐diffusion‐reaction equation. This new scheme could preserve the fixed points, the positivity, and the boundedness of the solution of the original equation. Numerical experiments verify the validity of our analytical results. Copyright © 2014 JohnWiley & Sons, Ltd.     

Analysis of algebraic systems arising from fourth‐order compact discretizations of convection‐diffusion equations

We study the properties of coefficient matrices arising from high‐order compact discretizations of convection‐diffusion problems. Asymptotic convergence factors of the convex hull of the spectrum and the field of values of the coefficient matrix for a one‐dimensional problem are derived, and the convergence factor of the convex hull of the spectrum is shown to be inadequate for predicting the convergence rate of GMRES. For a two‐dimensional constant‐coefficient problem, we derive the eigenvalues of the nine‐point matrix, and we show that the matrix is positive definite for all values of the cell‐Reynolds number. Using a recent technique for deriving analytic expressions for discrete solutions produced by the fourth‐order scheme, we show by analyzing the terms in the discrete solutions that they are oscillation‐free for all values of the cell Reynolds number. Our theoretical results support observations made through numerical experiments by other researchers on the non‐oscillatory nature of the discrete solution produced by fourth‐order compact approximations to the convection‐diffusion equation. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 155–178, 2002; DOI 10.1002/num.1041     

Graded Galerkin methods for the high‐order convection‐diffusion problem

We develop a Galerkin method using the Hermite spline on an admissible graded mesh for solving the high‐order singular perturbation problem of the convection‐diffusion type. We identify a special function class to which the solution of the convection‐diffusion problem belongs and characterize the approximation order of the Hermite spline for such a function class. The approximation order is then used to establish the optimal order of uniform convergence for the Galerkin method. Numerical results are presented to confirm the theoretical estimate.© 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Bifurcations and spatiotemporal patterns in a ratio‐dependent diffusive Holling‐Tanner system with time delay

In this paper, we concentrate on the spatiotemporal patterns of a delayed reaction‐diffusion Holling‐Tanner model with Neumann boundary conditions. In particular, the time delay that is incorporated in the negative feedback of the predator density is considered as one of the principal factors to affect the dynamic behavior. Firstly, a global Turing bifurcation theorem for τ = 0 and a local Turing bifurcation theorem for τ > 0 are given. Then, further considering the degenerated situation, we derive the existence of Bogdanov‐Takens bifurcation and Turing‐Hopf bifurcation. The normal form method is used to study the explicit dynamics near the Turing‐Hopf singularity. It is shown that a pair of stable nonconstant steady states (stripe patterns) and a pair of stable spatially inhomogeneous periodic solutions (spot patterns) could be bifurcated from a positive equilibrium. Moreover, the Turing‐Turing‐Hopf–type spatiotemporal patterns, that is, a subharmonic phenomenon with two spatial wave numbers and one temporal frequency, are also found and explained theoretically. Our results imply that the interaction of Turing and Hopf instabilities can be considered as the simplest mechanism for the appearance of complex spatiotemporal dynamics.