Self-and cross-diffusion are important nonlinear spatial derivative terms that are included into biological models of predator-prey interactions. Self-diffusion models overcrowding effects, while cross-diffusion incorporates the response of one species in light of the concentration of another. In this paper, a novel nonlinear operator splitting method is presented that directly incorporates both self-and cross-diffusion into a computational efficient design. The numerical analysis guarantees the accuracy and demonstrates appropriate criteria for stability. Numerical experiments display its efficiency and accuracy.
KEYWORDS
cross-diffusion, nonlinear operator splitting, reactiondiffusion equations, self-diffusionHere Ω = [0, L] × [0, L], Δ is the two-dimensional Laplacian operator, and the reactive functions f and g are in C 1 (R 2 ). We define x to be the spatial coordinate vector in two dimensions.Numer Methods Partial Differential Eq. 2019;35:597-614. wileyonlinelibrary.com/journal/num