Characterization of Turing diffusion-driven instability on evolving domains

In this paper we establish a general theoretical framework for Turing diffusion-driven instability for reaction-diffusion systems on time-dependent evolving domains. The main result is that Turing diffusion-driven instability for reaction-diffusion systems on evolving domains is characterised by Lyapunov exponents of the evolution family associated with the linearised system (obtained by linearising the original system along a spatially independent solution). This framework allows for the inclusion of the analysis of the long-time behavior of the solutions of reaction-diffusion systems. Applications to two special types of evolving domains are considered: (i) time-dependent domains which evolve to a final limiting fixed domain and (ii) time-dependent domains which are eventually time periodic. Reaction-diffusion systems have been widely proposed as plausible mechanisms for pattern formation in morphogenesis

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

Evolving reaction systems

Algorithms and the Foundations of Software technolog

Exhibiting cross-diffusion-induced patterns for reaction-diffusion systems on evolving domains and surfaces

The aim of this manuscript is to present for the first time the application of the finite element method for solving reaction-diffusion systems with cross-diffusion on continuously evolving domains and surfaces. Furthermore we present pattern formation generated by the reaction-diffusion systemwith cross-diffusion on evolving domains and surfaces. A two-component reaction-diffusion system with linear cross-diffusion in both u and v is presented. The finite element method is based on the approximation of the domain or surface by a triangulated domain or surface consisting of a union of triangles. For surfaces, the vertices of the triangulation lie on the continuous surface. A finite element space of functions is then defined by taking the continuous functions which are linear affine on each simplex of the triangulated domain or surface. To demonstrate the role of cross-diffusion to the theory of pattern formation, we compute patterns with model kinetic parameter values that belong only to the cross-diffusion parameter space; these do not belong to the standard parameter space for classical reaction-diffusion systems. Numerical results exhibited show the robustness, flexibility, versatility, and generality of our methodology; the methodology can deal with complicated evolution laws of the domain and surface, and these include uniform isotropic and anisotropic growth profiles as well as those profiles driven by chemical concentrations residing in the domain or on the surface

Reaction-diffusion systems on evolving domains with applications to the theory of biological pattern formation

In this thesis we investigate a model for biological pattern formation during growth development. The pattern formation phenomenon is described by a reaction-diffusion system on a time-dependent domain. We prove the global existence of solutions to reaction-diffusion systems on time-dependent domains.We extend global existence results for a class of reaction-diffusion systems on fixed domains to the same systems posed on spatially linear isotropically evolving domains. We demonstrate that the analysis is applicable to many systems that commonly arise in the theory of pattern formation. Our results give a mathematical justification to the widespread use of computer simulations of reaction-diffusion systems on evolving domains. We propose a finite element method to approximate the solutions to reaction-diffusion systems on time-dependent domains. We prove optimal convergence rates for the error in the method and we derive a computable error estimator that provides an upper bound for the error in the semidiscrete (space) scheme. We have implemented the method in the C programming language and we verify our theoretical results with benchmark computations. The method is a robust tool for the study of biological pattern formation, as it is applicable to domains with irregular geometries and nonuniform evolution. This versatility is illustrated with extensive computer simulations of reaction-diffusion systems on evolving domains. We observe varied pattern transitions induced by domain evolution, such as stripe to spot transitions, spotsplitting, spot-merging and spot-annihilation. We also illustrate the striking effects of spatially nonuniform domain evolution on the position, orientation and symmetry of patterns generated by reaction-diffusion systems. To improve the efficiency of the method, we have implemented a space-time adaptive algorithm where spatial adaptivity is driven by an error estimator and temporal adaptivity is driven by an error indicator.We illustrate with numerical simulations the dramatic improvements in accuracy and efficiency that are achieved via adaptivity. To demonstrate the applicability and generality of our methodology, we examine the process of parr mark pattern formation during the early development of the Amago trout. By assuming the existence of chemical concentrations residing on the surface of the Amago fish which react and diffuse during surface evolution, we model the pattern formation process with reactiondiffusion systems posed on evolving surfaces. An important generalisation of our study is the experimentally driven modelling of the fish’s developing body surface. Our results add weight to the feasibility of reaction-diffusion system models of fish skin patterning, by illustrating that a reaction-diffusion system posed on an evolving surface generates transient patterns consistent with those experimentally observed on the developing Amago trout. Furthermore, we conclude that the surface evolution profile, the surface geometry and the curvature are key factors which play a pivotal role in pattern formation via reaction-diffusion systems

Domain-growth-induced patterning for reaction-diffusion systems with linear cross-diffusion

In this article we present, for the first time, domain-growth induced pat- tern formation for reaction-diffusion systems with linear cross-diffusion on evolving domains and surfaces. Our major contribution is that by selecting parameter values from spaces induced by domain and surface evolution, patterns emerge only when domain growth is present. Such patterns do not exist in the absence of domain and surface evolution. In order to compute these domain-induced parameter spaces, linear stability theory is employed to establish the necessary conditions for domain- growth induced cross-diffusion-driven instability for reaction-diffusion systems with linear cross-diffusion. Model reaction-kinetic parameter values are then identified from parameter spaces induced by domain-growth only; these exist outside the classical standard Turing space on stationary domains and surfaces. To exhibit these patterns we employ the finite element method for solving reaction-diffusion systems with cross-diffusion on continuously evolving domains and surfaces

Transient behavior in Single-File Systems

We have used Monte-Carlo methods and analytical techniques to investigate the influence of the characteristics, such as pipe length, diffusion, adsorption, desorption and reaction rates on the transient properties of Single-File Systems. The transient or the relaxation regime is the period in which the system is evolving to equilibrium. We have studied the system when all the sites are reactive and when only some of them are reactive. Comparisons between Mean-Field predictions, Cluster Approximation predictions, and Monte Carlo simulations for the relaxation time of the system are shown. We outline the cases where Mean-Field analysis gives good results compared to Dynamic Monte-Carlo results. For some specific cases we can analytically derive the relaxation time. Occupancy profiles for different distribution of the sites both for Mean-Field and simulations are compared. Different results for slow and fast reaction systems and different distribution of reactive sites are discussed.Comment: 18 pages, 19 figure

Non-equilibrium transitions in multiscale systems with a bifurcating slow manifold

Noise-induced transitions between metastable fixed points in systems evolving on multiple time scales are analyzed in situations where the time scale separation gives rise to a slow manifold with bifurcation. This analysis is performed within the realm of large deviation theory. It is shown that these non-equilibrium transitions make use of a reaction channel created by the bifurcation structure of the slow manifold, leading to vastly increased transition rates. Several examples are used to illustrate these findings, including an insect outbreak model, a system modeling phase separation in the presence of evaporation, and a system modeling transitions in active matter self-assembly. The last example involves a spatially extended system modeled by a stochastic partial differential equation

Post-Newtonian gravitational radiation and equations of motion via direct integration of the relaxed Einstein equations. III. Radiation reaction for binary systems with spinning bodies

Using post-Newtonian equations of motion for fluid bodies that include radiation-reaction terms at 2.5 and 3.5 post-Newtonian (PN) order (O[(v/c)^5] and O[(v/c)^7] beyond Newtonian order), we derive the equations of motion for binary systems with spinning bodies. In particular we determine the effects of radiation-reaction coupled to spin-orbit effects on the two-body equations of motion, and on the evolution of the spins. For a suitable definition of spin, we reproduce the standard equations of motion and spin-precession at the first post-Newtonian order. At 3.5PN order, we determine the spin-orbit induced reaction effects on the orbital motion, but we find that radiation damping has no effect on either the magnitude or the direction of the spins. Using the equations of motion, we find that the loss of total energy and total angular momentum induced by spin-orbit effects precisely balances the radiative flux of those quantities calculated by Kidder et al. The equations of motion may be useful for evolving inspiraling orbits of compact spinning binaries.Comment: 19 pages, small corrections, equivalent to published versio