Reaction-diffusion Equations (reaction-diffusion + equation)

Distribution by Scientific Domains


Selected Abstracts


THE DIFFUSIVE SPREAD OF ALLELES IN HETEROGENEOUS POPULATIONS

EVOLUTION, Issue 3 2004
Garrick T. Skalski
Abstract The spread of genes and individuals through space in populations is relevant in many biological contexts. I study, via systems of reaction-diffusion equations, the spatial spread of advantageous alleles through structured populations. The results show that the temporally asymptotic rate of spread of an advantageous allele, a kind of invasion speed, can be approximated for a class of linear partial differential equations via a relatively simple formula, c= 2,rD, that is reminiscent of a classic formula attributed to R. A. Fisher. The parameters r and D, represent an asymptotic growth rate and an average diffusion rate, respectively, and can be interpreted in terms of eigenvalues and eigenvectors that depend on the population's demographic structure. The results can be applied, under certain conditions, to a wide class of nonlinear partial differential equations that are relevant to a variety of ecological and evolutionary scenarios in population biology. I illustrate the approach for computing invasion speed with three examples that allow for heterogeneous dispersal rates among different classes of individuals within model populations. [source]


A fourth-order accurate, Numerov-type, three-point finite-difference discretization of electrochemical reaction-diffusion equations on nonuniform (exponentially expanding) spatial grids in one-dimensional space geometry

JOURNAL OF COMPUTATIONAL CHEMISTRY, Issue 12 2004
aw K. Bieniasz
Abstract The validity for finite-difference electrochemical kinetic simulations, of the extension of the Numerov discretization designed by Chawla and Katti [J Comput Appl Math 1980, 6, 189,196] for the solution of two-point boundary value problems in ordinary differential equations, is examined. The discretization is adapted to systems of time-dependent reaction-diffusion partial differential equations in one-dimensional space geometry, on nonuniform space grids resulting from coordinate transformations. The equations must not involve first spatial derivatives of the unknowns. Relevant discrete formulae are outlined and tested in calculations on two example kinetic models. The models describe potential step chronoamperometry under limiting current conditions, for the catalytic EC, and Reinert-Berg CE reaction mechanisms. Exponentially expanding space grid is used. The discretization considered proves the most accurate and efficient, out of all the three-point finite-difference discretizations on such grids, that have been used thus far in electrochemical kinetics. Therefore, it can be recommended as a method of choice. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1515,1521, 2004 [source]


A fourth-order compact algorithm for nonlinear reaction-diffusion equations with Neumann boundary conditions

NUMERICAL METHODS FOR PARTIAL DIFFERENTIAL EQUATIONS, Issue 3 2006
Wenyuan Liao
Abstract In this article, we discuss a scheme for dealing with Neumann and mixed boundary conditions using a compact stencil. The resulting compact algorithm for solving systems of nonlinear reaction-diffusion equations is fourth-order accurate in both the temporal and spatial dimensions. We also prove that the standard second-order approximation to zero Neumann boundary conditions provides fourth-order accuracy when the nonlinear reaction term is independent of the spatial variables. Numerical examples, including an application of this algorithm to a mathematical model describing frontal polymerization process, are presented in the article to demonstrate the accuracy and efficiency of the scheme. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005 [source]


An efficient high-order algorithm for solving systems of reaction-diffusion equations

NUMERICAL METHODS FOR PARTIAL DIFFERENTIAL EQUATIONS, Issue 3 2002
Wenyuan Liao
Abstract An efficient higher-order finite difference algorithm is presented in this article for solving systems of two-dimensional reaction-diffusion equations with nonlinear reaction terms. The method is fourth-order accurate in both the temporal and spatial dimensions. It requires only a regular five-point difference stencil similar to that used in the standard second-order algorithm, such as the Crank-Nicolson algorithm. The Padé approximation and Richardson extrapolation are used to achieve high-order accuracy in the spatial and temporal dimensions, respectively. Numerical examples are presented to demonstrate the efficiency and accuracy of the new algorithm. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 340,354, 2002; Published online in Wiley InterScience (www.interscience.wiley.com). DOI 10.1002/num.10012 [source]


Linear vs. nonlinear selection for the propagation speed of the solutions of scalar reaction-diffusion equations invading an unstable equilibrium

COMMUNICATIONS ON PURE & APPLIED MATHEMATICS, Issue 5 2004
Marcello Lucia
We revisit the classical problem of speed selection for the propagation of disturbances in scalar reaction-diffusion equations with one linearly stable and one linearly unstable equilibrium. For a wide class of initial data this problem reduces to finding the minimal speed of the monotone traveling wave solutions connecting these two equilibria in one space dimension. We introduce a variational characterization of these traveling wave solutions and give a necessary and sufficient condition for linear versus nonlinear selection mechanism. We obtain sufficient conditions for the linear and nonlinear selection mechanisms that are easily verifiable. Our method also allows us to obtain efficient lower and upper bounds for the propagation speed. © 2004 Wiley Periodicals, Inc. [source]