Cellular automaton model for molecular traffic jams

We consider the time evolution of an exactly solvable cellular automaton with random initial conditions both in the large-scale hydrodynamic limit and on the microscopic level. This model is a version of the totally asymmetric simple exclusion process with sublattice parallel update and thus may serve as a model for studying traffic jams in systems of self-driven particles. We study the emergence of shocks from the microscopic dynamics of the model. In particular, we introduce shock measures whose time evolution we can compute explicitly, both in the thermodynamic limit and for open boundaries where a boundary-induced phase transition driven by the motion of a shock occurs. The motion of the shock, which results from the collective dynamics of the exclusion particles, is a random walk with an internal degree of freedom that determines the jump direction. This type of hopping dynamics is reminiscent of some transport phenomena in biological systems.

Modeling the Potential Spread of the Recently Identified Non-Native Panther Grouper (Chromileptes altivelis)in the Atlantic Using a Cellular Automaton Approach

The Indo-pacific panther grouper (Chromileptes altiveli) is a predatory fish species and popular imported aquarium fish in the United States which has been recently documented residing in western Atlantic waters. To date, the most successful marine invasive species in the Atlantic is the lionfish (Pterois volitans/miles), which, as for the panther grouper, is assumed to have been introduced to the wild through aquarium releases. However, unlike lionfish, the panther grouper is not yet thought to have an established breeding population in the Atlantic. Using a proven modeling technique developed to track the lionfish invasion, presented is the first known estimation of the potential spread of panther grouper in the Atlantic. The employed cellular automaton-based computer model examines the life history of the subject species including fecundity, mortality, and reproductive potential and combines this with habitat preferences and physical oceanic parameters to forecast the distribution and periodicity of spread of this potential new invasive species. Simulations were examined for origination points within one degree of capture locations of panther grouper from the United States Geological Survey Nonindigenous Aquatic Species Database to eliminate introduction location bias, and two detailed case studies were scrutinized. The model indicates three primary locations where settlement is likely given the inputs and limits of the model; Jupiter Florida/Vero Beach, the Cape Hatteras Tropical Limit/Myrtle Beach South Carolina, and Florida Keys/Ten Thousand Islands locations. Of these locations, Jupiter Florida/Vero Beach has the highest settlement rate in the model and is indicated as the area in which the panther grouper is most likely to become established. This insight is valuable if attempts are to be made to halt this potential marine invasive species

An extinction-survival-type phase transition in the probabilistic cellular automaton p182-q200

We investigate the critical behaviour of a probabilistic mixture of cellular automata (CA) rules 182 and 200 (in Wolfram`s enumeration scheme) by mean-field analysis and Monte Carlo simulations. We found that as we switch off one CA and switch on the other by the variation of the single parameter of the model, the probabilistic CA (PCA) goes through an extinction-survival-type phase transition, and the numerical data indicate that it belongs to the directed percolation universality class of critical behaviour. The PCA displays a characteristic stationary density profile and a slow, diffusive dynamics close to the pure CA 200 point that we discuss briefly. Remarks on an interesting related stochastic lattice gas are addressed in the conclusions.

A Cellular Automaton Approach to Spatial Electric Load Forecasting

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Recurrence interval statistics of cellular automaton seismicity models

The recurrence interval statistics for regional seismicity follows a universal distribution function, independent of the tectonic setting or average rate of activity (Corral, 2004). The universal function is a modified gamma distribution with power-law scaling of recurrence intervals shorter than the average rate of activity and exponential decay for larger intervals. We employ the method of Corral (2004) to examine the recurrence statistics of a range of cellular automaton earthquake models. The majority of models has an exponential distribution of recurrence intervals, the same as that of a Poisson process. One model, the Olami-Feder-Christensen automaton, has recurrence statistics consistent with regional seismicity for a certain range of the conservation parameter of that model. For conservation parameters in this range, the event size statistics are also consistent with regional seismicity. Models whose dynamics are dominated by characteristic earthquakes do not appear to display universality of recurrence statistics.

A cellular automata model to investigate immune cell–tumor cell interactions in growing tumors in two spatial dimensions

We develop a hybrid cellular automata model to describe the effect of the immune system and chemokines on a growing tumor. The hybrid cellular automata model consists of partial differential equations to model chemokine concentrations, and discrete cellular automata to model cell–cell interactions and changes. The computational implementation overlays these two components on the same spatial region. We present representative simulations of the model and show that increasing the number of immature dendritic cells (DCs) in the domain causes a decrease in the number of tumor cells. This result strongly supports the hypothesis that DCs can be used as a cancer treatment. Furthermore, we also use the hybrid cellular automata model to investigate the growth of a tumor in a number of computational “cancer patients.” Using these virtual patients, the model can explain that increasing the number of DCs in the domain causes longer “survival.” Not surprisingly, the model also reflects the fact that the parameter related to tumor division rate plays an important role in tumor metastasis.

A hybrid cellular automata model of multicellular tumour spheroid growth in hypoxic microenvironment

A three-dimensional hybrid cellular automata (CA) model is developed to study the dynamic process of multicellular tumour spheroid (MTS) growth by introducing hypoxia as an important microenvironment factor which influences cell migration and cell phenotype expression. The model enables us to examine the effects of different hypoxic environments on the growth history of the MTS and to study the dynamic interactions between MTS growth and chemical environments. The results include the spatial distribution of different phenotypes of tumour cells and associated oxygen concentration distributions under hypoxic conditions. The discussion of the model system responses to the varied hypoxic conditions reveals that the improvement of the resistance of tumour cells to a hypoxic environment may be important in the tumour normalization therapy.

Analysis of the compressive behaviour of the three-dimensional printed porous titanium for dental implants using a modified cellular solid model

A set of cylindrical porous titanium test samples were produced using the three-dimensional printing and sintering method with samples sintered at 900 °C, 1000 °C, 1100 °C, 1200 °C or 1300 °C. Following compression testing, it was apparent that the stress-strain curves were similar in shape to the curves that represent cellular solids. This is despite a relative density twice as high as what is considered the threshold for defining a cellular solid. As final sintering temperature increased, the compressive behaviour developed from being elastic-brittle to elastic-plastic and while Young's modulus remained fairly constant in the region of 1.5 GPa, there was a corresponding increase in 0.2% proof stress of approximately 40-80 MPa. The cellular solid model consists of two equations that predict Young's modulus and yield or proof stress. By fitting to experimental data and consideration of porous morphology, appropriate changes to the geometry constants allow modification of the current models to predict with better accuracy the behaviour of porous materials with higher relative densities (lower porosity).

Coupling biomechanics to a cellular level model: An approach to patient-specific image driven multi-scale and multi-physics tumor simulation

Modeling of tumor growth has been performed according to various approaches addressing different biocomplexity levels and spatiotemporal scales. Mathematical treatments range from partial differential equation based diffusion models to rule-based cellular level simulators, aiming at both improving our quantitative understanding of the underlying biological processes and, in the mid- and long term, constructing reliable multi-scale predictive platforms to support patient-individualized treatment planning and optimization. The aim of this paper is to establish a multi-scale and multi-physics approach to tumor modeling taking into account both the cellular and the macroscopic mechanical level. Therefore, an already developed biomodel of clinical tumor growth and response to treatment is self-consistently coupled with a biomechanical model. Results are presented for the free growth case of the imageable component of an initially point-like glioblastoma multiforme tumor. The composite model leads to significant tumor shape corrections that are achieved through the utilization of environmental pressure information and the application of biomechanical principles. Using the ratio of smallest to largest moment of inertia of the tumor material to quantify the effect of our coupled approach, we have found a tumor shape correction of 20\% by coupling biomechanics to the cellular simulator as compared to a cellular simulation without preferred growth directions. We conclude that the integration of the two models provides additional morphological insight into realistic tumor growth behavior. Therefore, it might be used for the development of an advanced oncosimulator focusing on tumor types for which morphology plays an important role in surgical and/or radio-therapeutic treatment planning.

The influence of connectivity on the firing rate in a neuronal network with electrical and chemical synapses

We study the firing rate properties of a cellular automaton model for a neuronal network with chemical synapses. We propose a simple mechanism in which the nonlocal connections are included, through electrical and chemical synapses. In the latter case, we introduce a time delay which produces self-sustained activity. Nonlocal connections, or shortcuts, are randomly introduced according to a specified connection probability. There is a range of connection probabilities for which neuron firing occurs, as well as a critical probability for which the firing ceases in the absence of time delay. The critical probability for nonlocal shortcuts depends on the network size according to a power-law. We also compute the firing rate amplification factor by varying both the connection probability and the time delay for different network sizes. (C) 2011 Elsevier B.V. All rights reserved.

Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

Stochastic lattice gas model describing the dynamics of the SIRS epidemic process

We study a stochastic process describing the onset of spreading dynamics of an epidemic in a population composed of individuals of three classes: susceptible (S), infected (I), and recovered (R). The stochastic process is defined by local rules and involves the following cyclic process: S -> I -> R -> S (SIRS). The open process S -> I -> R (SIR) is studied as a particular case of the SIRS process. The epidemic process is analyzed at different levels of description: by a stochastic lattice gas model and by a birth and death process. By means of Monte Carlo simulations and dynamical mean-field approximations we show that the SIRS stochastic lattice gas model exhibit a line of critical points separating the two phases: an absorbing phase where the lattice is completely full of S individuals and an active phase where S, I and R individuals coexist, which may or may not present population cycles. The critical line, that corresponds to the onset of epidemic spreading, is shown to belong in the directed percolation universality class. By considering the birth and death process we analyze the role of noise in stabilizing the oscillations. (C) 2009 Elsevier B.V. All rights reserved.

Modelling fungus dispersal scenarios using cellular automata

This study focused on representing spatio-temporal patterns of fungal dispersal using cellular automata. Square lattices were used, with each site representing a host for a hypothetical fungus population. Four possible host states were allowed: resistant, permissive, latent or infectious. In this model, the probability of infection for each of the healthy states (permissive or resistant) in a time step was determined as a function of the host's susceptibility, seasonality, and the number of infectious sites and the distance between them. It was also assumed that infected sites become infectious after a pre-specified latency period, and that recovery is not possible. Several scenarios were simulated to understand the contribution of the model's parameters and the spatial structure on the dynamic behaviour of the modelling system. The model showed good capability for representing the spatio-temporal pattern of fungus dispersal over planar surfaces. With a specific problem in mind, the model can be easily modified and used to describe field behaviour, which can contribute to the conservation and development of management strategies for both natural and agricultural systems. © 2012 Elsevier B.V.