2,748 research outputs found
A Semi-Automatic Numerical Algorithm for Turing Patterns Formation in a Reaction-Diffusion Model
The Turing pattern formation is modeled by reaction - diffusion (RD) type partial differential equations , and it plays a crucial role in ecological studies. Big data analytics and suitable frameworks to manage and predict structures and configurations are mandatory. The processing and resolution procedures of mathematical models relies upon numerical schemes, and concurrently upon the related automated algorithms. Starting from a RD model for vegetation patterns, we propose a semi-automatic algorithm based on a smart numerical criterion for observing ecological reliable results. Numerical experiments are carried out in the case of spot's formations
Extended patchy ecosystems may increase their total biomass through self-replication
Patches of vegetation consist of dense clusters of shrubs, grass, or trees,
often found to be circular characteristic size, defined by the properties of
the vegetation and terrain. Therefore, vegetation patches can be interpreted as
localized structures. Previous findings have shown that such localized
structures can self-replicate in a binary fashion, where a single vegetation
patch elongates and divides into two new patches. Here, we extend these
previous results by considering the more general case, where the plants
interact non-locally, this extension adds an extra level of complexity and
shrinks the gap between the model and real ecosystems, where it is known that
the plant-to-plant competition through roots and above-ground facilitating
interactions have non-local effects, i.e. they extend further away than the
nearest neighbor distance. Through numerical simulations, we show that for a
moderate level of aridity, a transition from a single patch to periodic pattern
occurs. Moreover, for large values of the hydric stress, we predict an opposing
route to the formation of periodic patterns, where a homogeneous cover of
vegetation may decay to spot-like patterns. The evolution of the biomass of
vegetation patches can be used as an indicator of the state of an ecosystem,
this allows to distinguish if a system is in a self-replicating or decaying
dynamics. In an attempt to relate the theoretical predictions to real
ecosystems, we analyze landscapes in Zambia and Mozambique, where vegetation
forms patches of tens of meters in diameter. We show that the properties of the
patches together with their spatial distributions are consistent with the
self-organization hypothesis. We argue that the characteristics of the observed
landscapes may be a consequence of patch self-replication, however, detailed
field and temporal data is fundamental to assess the real state of the
ecosystems.Comment: 38 pages, 12 figures, 1 tabl
Parameter estimation for a morphochemical reaction-diffusion model of electrochemical pattern formation
The process of electrodeposition can be described in terms of a reaction-diffusion PDE system that models the dynamics of the morphology profile and the chemical composition. Here we fit such a model to the different patterns present in a range of electrodeposited and electrochemically modified alloys using PDE constrained optimization. Experiments with simulated data show how the parameter space of the model can be divided into zones corresponding to the different physical patterns by examining the structure of an appropriate cost function. We then use real data to demonstrate how numerical optimization of the cost function can allow the model to fit the rich variety of patterns arising in experiments. The computational technique developed provides a potential tool for tuning experimental parameters to produce desired patterns
Analysis of pattern dynamics for a nonlinear model of the human cortex via bifurcation theories
This thesis examines the bifurcations, i.e., the emergent behaviours, for the Waikato cortical model under the influence of the gap-junction inhibitory diffusion Dâ (identified as the Turing bifurcation parameter) and the time-to-peak for hyperpolarising GABA response Îłi (i.e., inhibitory rate-constant, identified as the Hopf bifurcation parameter). The cortical model simplifies the entire cortex to a cylindrical macrocolumn (⌠1 mmÂł) containing ⌠10â” neurons (85% excitatory, 15% inhibitory) communicating via both chemical and electrical (gap-junction) synapses. The linear stability analysis of the model equations predict the emergence of a Turing instability (in which separated areas of the cortex become activated) when gap-junction diffusivity is increased above a critical level. In addition, a Hopf bifurcation (oscillation) occurs when the inhibitory rate-constant is sufficiently small. Nonlinear interaction between these instabilities leads to spontaneous cortical patterns of neuronal activities evolving in space and time. Such model dynamics of delicately balanced interplay between Turing and Hopf instabilities may be of direct relevance to clinically observed brain dynamics such as epileptic seizure EEG spikes, deep-sleep slow-wave oscillations and cognitive gamma-waves.
The relationship between the modelled brain patterns and model equations can normally be inferred from the eigenvalue dispersion curve, i.e., linear stability analysis. Sometimes we experienced mismatches between the linear stability analysis and the formed cortical patterns, which hampers us in identifying the type of instability corresponding to the emergent patterns. In this thesis, I investigate the pattern-forming mechanism of the Waikato cortical model to better understand the model nonlinearities. I first study the pattern dynamics via analysis of a simple pattern-forming system, the Brusselator model, which has a similar model structure and bifurcation phenomena as the cortical model. I apply both linear and nonlinear perturbation methods to analyse the near-bifurcation behaviour of the Brusselator in order to precisely capture the dominant mode that contributes the most to the final formed-patterns. My nonlinear analysis of the Brusselator model yields Ginzburg-Landau type amplitude equations that describe the dynamics of the most unstable mode, i.e., the dominant mode, in the vicinity of a bifurcation point. The amplitude equations at a Turing point unfold three characteristic spatial structures: honeycomb HÏ, stripes, and reentrant honeycomb Hâ. A codimension-2 TuringâHopf point (CTHP) predicts three mixed instabilities: stable TuringâHopf (TH), chaotic TH, and bistable TH. The amplitude equations precisely determine the bifurcation conditions for these instabilities and explain the pattern-competition mechanism once the bifurcation parameters cross the thresholds, whilst driving the system into a nonlinear region where the linear stability analysis may not be applicable.
Then, I apply the bifurcation theories to the cortical model for its pattern predictions. Analogous to the Brusselator model, I find cortical Turing pattens in HÏ, stripes and Hâ spatial structures. Moreover, I develop the amplitude equations for the cortical model, with which I derive the envelope frequency for the beating-waves of a stable TH mode; and propose ideas regarding emergence of the cortical chaotic mode. Apart from these pattern dynamics that the cortical model shares with the Brusselator system, the cortical model also exhibits âeye-blinkingâ TH patterns latticed in hexagons with localised oscillations. Although we have not found biological significance of these model pattens, the developed bifurcation theories and investigated pattern-forming mechanism may enrich our modelling strategies and help us to further improve model performance.
In the last chapter of this thesis, I introduce a TuringâHopf mechanism for the anaesthetic slow-waves, and predict a coherence drop of such slow-waves with the induction of propofol anaesthesia. To test this hypothesis, I developed an EEG coherence analysing algorithm, EEG coherence, to automatically examine the clinical EEG recordings across multiple subjects. The result shows significantly decreased coherence along the fronto-occipital axis, and increased coherence along the left- and right-temporal axis. As the Waikato cortical model is spatially homogenous, i.e., there are no explicit front-to-back or right-to-left directions, it is unable to produce different coherence changes for different regions. It appears that the Waikato cortical model best represents the cortical dynamics in the frontal region. The theory of pattern dynamics suggests that a mode transition from waveâTuringâwave to TuringâwaveâTuring introduces pattern coherence changes in both positive and negative directions. Thus, a further modelling improvement may be the introduction of a cortical bistable mode where Turing and wave coexist
Turing Patterning in Stratified Domains
Reaction-diffusion processes across layered media arise in several scientific
domains such as pattern-forming E. coli on agar substrates,
epidermal-mesenchymal coupling in development, and symmetry-breaking in cell
polarisation. We develop a modelling framework for bi-layer reaction-diffusion
systems and relate it to a range of existing models. We derive conditions for
diffusion-driven instability of a spatially homogeneous equilibrium analogous
to the classical conditions for a Turing instability in the simplest nontrivial
setting where one domain has a standard reaction-diffusion system, and the
other permits only diffusion. Due to the transverse coupling between these two
regions, standard techniques for computing eigenfunctions of the Laplacian
cannot be applied, and so we propose an alternative method to compute the
dispersion relation directly. We compare instability conditions with full
numerical simulations to demonstrate impacts of the geometry and coupling
parameters on patterning, and explore various experimentally-relevant
asymptotic regimes. In the regime where the first domain is suitably thin, we
recover a simple modulation of the standard Turing conditions, and find that
often the broad impact of the diffusion-only domain is to reduce the ability of
the system to form patterns. We also demonstrate complex impacts of this
coupling on pattern formation. For instance, we exhibit non-monotonicity of
pattern-forming instabilities with respect to geometric and coupling
parameters, and highlight an instability from a nontrivial interaction between
kinetics in one domain and diffusion in the other. These results are valuable
for informing design choices in applications such as synthetic engineering of
Turing patterns, but also for understanding the role of stratified media in
modulating pattern-forming processes in developmental biology and beyond.Comment: 25 pages, 7 figure
- âŠ