An ADI extrapolated Crank-Nicolson orthogonal spline collocation method for nonlinear reaction-diffusion systems: a computational study

An alternating direction implicit (ADI) orthogonal spline collocation (OSC) method is described for the approximate solution of a class of nonlinear reaction-diffusion systems. Its efficacy is demonstrated on the solution of well-known examples of such systems, specifically the Brusselator, Gray-Scott, Gierer-Meinhardt and Schnakenberg models, and comparisons are made with other numerical techniques considered in the literature. The new ADI method is based on an extrapolated Crank-Nicolson OSC method and is algebraically linear. It is efficient, requiring at each time level only $O({\cal N})$ operations where ${\cal N}$ is the number of unknowns. Moreover,it is shown to produce approximations which are of optimal global accuracy in various norms, and to possess superconvergence properties

Global existence for semilinear reaction-diffusion systems on evolving domains

We present global existence results for solutions of reaction-diffusion systems on evolving domains. Global existence results for a class of reaction-diffusion systems on fixed domains are extended to the same systems posed on spatially linear isotropically evolving domains. The results hold without any assumptions on the sign of the growth rate. The analysis is valid for many systems that commonly arise in the theory of pattern formation. We present numerical results illustrating our theoretical findings.Comment: 24 pages, 3 figure

Adaptive multiresolution computations applied to detonations

A space-time adaptive method is presented for the reactive Euler equations describing chemically reacting gas flow where a two species model is used for the chemistry. The governing equations are discretized with a finite volume method and dynamic space adaptivity is introduced using multiresolution analysis. A time splitting method of Strang is applied to be able to consider stiff problems while keeping the method explicit. For time adaptivity an improved Runge--Kutta--Fehlberg scheme is used. Applications deal with detonation problems in one and two space dimensions. A comparison of the adaptive scheme with reference computations on a regular grid allow to assess the accuracy and the computational efficiency, in terms of CPU time and memory requirements.Comment: Zeitschrift f\"ur Physicalische Chemie, accepte

Integrating actin and myosin II in a viscous model for cell migration

This article presents a mathematical and computational model for cell migration that couples a system of reaction-advection-diffusion equations describing the bio-molecular interactions between F-actin and myosin II to a force balance equation describing the structural mechanics of the actin-myosin network. In eukaryotic cells, cell migration is largely powered by a system of actin and myosin dynamics. We formulate the model equations on a two-dimensional cellular migrating evolving domain to take into account the convective and dilution terms for the biochemical reaction-diffusion equations, with hypothetically proposed reaction-kinetics. We employ the evolving finite element method to compute approximate numerical solutions of the coupled biomechanical model in two dimensions. Numerical experiments exhibit cell polarization through symmetry breaking which are driven by the F-actin and myosin II. This conceptual hypothetical proof-of-concept framework set premises for studying experimentally-driven actin-myosin reaction-kinetic network interactions with generalizations to multi-dimensions

Stability analysis of reaction-diffusion models on evolving domains: the effects of cross-diffusion

This article presents stability analytical results of a two component reaction-diffusion system with linear cross-diffusion posed on continuously evolving domains. First the model system is mapped from a continuously evolving domain to a reference stationary frame resulting in a system of partial differential equations with time-dependent coefficients. Second, by employing appropriately asymptotic theory, we derive and prove cross-diffusion-driven instability conditions for the model system for the case of slow, isotropic domain growth. Our analytical results reveal that unlike the restrictive diffusion-driven instability conditions on stationary domains, in the presence of cross-diffusion coupled with domain evolution, it is no longer necessary to enforce cross nor pure kinetic conditions. The restriction to activator-inhibitor kinetics to induce pattern formation on a growing biological system is no longer a requirement. Reaction-cross-diffusion models with equal diffusion coefficients in the principal components as well as those of the short-range inhibition, long-range activation and activator-activator form can generate patterns only in the presence of cross-diffusion coupled with domain evolution. To confirm our theoretical findings, detailed parameter spaces are exhibited for the special cases of isotropic exponential, linear and logistic growth profiles. In support of our theoretical predictions, we present evolving or moving finite element solutions exhibiting patterns generated by a short-range inhibition, long-range activation reaction-diffusion model with linear cross-diffusion in the presence of domain evolution