Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn efficient, nonlinear stability analysis for detecting pattern formation in reaction diffusion systems
231
0
0.0
(
0
)
Authors
William R. Holmes
Ask ChatGPT about the research
No Arabic abstract
Reaction diffusion systems are often used to study pattern formation in biological systems. However, most methods for understanding their behavior are challenging and can rarely be applied to complex systems common in biological applications. I present a relatively simple and efficient, non-linear stability technique that greatly aids such analysis when rates of diffusion are substantially different. This technique reduces a system of reaction diffusion equations to a system of ordinary differential equations tracking the evolution of a large amplitude, spatially localized perturbation of a homogeneous steady state. Stability properties of this system, determined using standard bifurcation techniques and software, describe both linear and non-linear patterning regimes of the reaction diffusion system. I describe the class of systems this method can be applied to and demonstrate its application. Analysis of Schnakenberg and substrate inhibition models is performed to demonstrate the methods capabilities in simplified settings and show that even these simple models have non-linear patterning regimes not previously detected. Analysis of a protein regulatory network related to chemotaxis shows its application in a more complex setting where other non-linear methods become intractable. Predictions of this method are verified against results of numerical simulation, linear stability, and full PDE bifurcation analyses.
rate research
Read More
A nonlinear PDE featuring flux limitation effects together with those of the porous media equation (nonlinear Fokker-Planck) is presented in this paper. We analyze the balance of such diverse effects through the study of the existence and qualitative behavior of some admissible patterns, namely traveling wave solutions, to this singular reaction-difusion equation. We show the existence and qualitative behavior of different types of traveling waves: classical profiles for wave speeds high enough, and discontinuous waves that are reminiscent of hyperbolic shock waves when the wave speed lowers below a certain threshold. Some of these solutions are of particular relevance as they provide models by which the whole solution (and not just the bulk of it, as it is the case with classical traveling waves) spreads through the medium with finite speed.
We establish global-in-time existence results for thermodynamically consistent reaction-(cross-)diffusion systems coupled to an equation describing heat transfer. Our main interest is to model species-dependent diffusivities, while at the same time ensuring thermodynamic consistency. A key difficulty of the non-isothermal case lies in the intrinsic presence of cross-diffusion type phenomena like the Soret and the Dufour effect: due to the temperature/energy dependence of the thermodynamic equilibria, a nonvanishing temperature gradient may drive a concentration flux even in a situation with constant concentrations; likewise, a nonvanishing concentration gradient may drive a heat flux even in a case of spatially constant temperature. We use time discretisation and regularisation techniques and derive a priori estimates based on a suitable entropy and the associated entropy production. Renormalised solutions are used in cases where non-integrable diffusion fluxes or reaction terms appear.
Using pointwise semigroup techniques, we establish sharp rates of decay in space and time of a perturbed reaction diffusion front to its time-asymptotic limit. This recovers results of Sattinger, Henry and others of time-exponential convergence in weighted $L^p$ and Sobolev norms, while capturing the new feature of spatial diffusion at Gaussian rate. Novel features of the argument are a point-wise Green function decomposition reconciling spectral decomposition and short-time Nash-Aronson estimates and an instantaneous tracking scheme similar to that used in the study of stability of viscous shock waves.
The convergence to equilibrium for renormalised solutions to nonlinear reaction-diffusion systems is studied. The considered reaction-diffusion systems arise from chemical reaction networks with mass action kinetics and satisfy the complex balanced condition. By applying the so-called entropy method, we show that if the system does not have boundary equilibria, then any renormalised solution converges exponentially to the complex balanced equilibrium with a rate, which can be computed explicitly up to a finite dimensional inequality. This inequality is proven via a contradiction argument and thus not explicitly. An explicit method of proof, however, is provided for a specific application modelling a reversible enzyme reaction by exploiting the specific structure of the conservation laws. Our approach is also useful to study the trend to equilibrium for systems possessing boundary equilibria. More precisely, to show the convergence to equilibrium for systems with boundary equilibria, we establish a sufficient condition in terms of a modified finite dimensional inequality along trajectories of the system. By assuming this condition, which roughly means that the system produces too much entropy to stay close to a boundary equilibrium for infinite time, the entropy method shows exponential convergence to equilibrium for renormalised solutions to complex balanced systems with boundary equilibria.
A self-consistent equation to derive a discreteness-induced stochastic steady state is presented for reaction-diffusion systems. For this formalism, we use the so-called Kuramoto length, a typical distance over which a molecule diffuses in its lifetime, as was originally introduced to determine if local fluctuations influence globally the whole system. We show that this Kuramoto length is also relevant to determine whether the discreteness of molecules is significant or not. If the number of molecules of a certain species within the Kuramoto length is small and discrete, localization of some other chemicals is brought about, which can accelerate certain reactions. When this acceleration influences the concentration of the original molecule species, it is shown that a novel, stochastic steady state is induced that does not appear in the continuum limit. A theory to obtain and characterize this state is introduced, based on the self-consistent equation for chemical concentrations. This stochastic steady state is confirmed by numerical simulations on a certain reaction model, which agrees well with the theoretical estimation. Formation and coexistence of domains with different stochastic states are also reported, which is maintained by the discreteness. Relevance of our result to intracellular reactions is briefly discussed.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
