Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCellular automaton model of reaction-transport processes
122
0
0.0
(
0
)
Authors
T. Karapiperis
Ask ChatGPT about the research
No Arabic abstract
The transport and chemical reactions of solutes are modelled as a cellular automaton in which molecules of different species perform a random walk on a regular lattice and react according to a local probabilistic rule. The model describes advection and diffusion in a simple way, and as no restriction is placed on the number of particles at a lattice site, it is also able to describe a wide variety of chemical reactions. Assuming molecular chaos and a smooth density function, we obtain the standard reaction-transport equations in the continuum limit. Simulations on one- and two-dimensional lattices show that the discrete model can be used to approximate the solutions of the continuum equations. We discuss discrepancies which arise from correlations between molecules and how these disappear as the continuum limit is approached. Of particular interest are simulations displaying long-time behaviour which depends on long-wavelength statistical fluctuations not accounted for by the standard equations. The model is applied to the reactions $a + b rightleftharpoons c$ and $a + b rightarrow c$ with homogeneous and inhomogeneous initial conditions as well as to systems subject to autocatalytic reactions and displaying spontaneous formation of spatial concentration patterns.
rate research
Read More
Precipitation/dissolution reactions coupled with solute transport are modelled as a cellular automaton in which solute molecules perform a random walk on a regular lattice and react according to a local probabilistic rule. Stationary solid particles dissolve with a certain probability and, provided solid is already present or the solution is saturated, solute particles have a probability to precipitate. In our simulation of the dissolution of a solid block inside uniformly flowing water we obtain solid precipitation downstream from the original solid edge, in contrast to the standard reaction-transport equations. The observed effect is the result of fluctuations in solute density and diminishes when we average over a larger ensemble. The additional precipitation of solid is accompanied by a substantial reduction in the relatively small solute concentration. The model is appropriate for the study of the r^ole of intrinsic fluctuations in the presence of reaction thresholds and can be employed to investigate porosity changes associated with the carbonation of cement.
The state of structural balance (termed also `Heider balance) of a social network is often discussed in social psychology and sociophysics. In this state, actors at network nodes classify other individuals as enemies or friends. Hence, the network contains two kinds of links: positive and negative. Here a new cellular automaton is designed and investigated, which mimics the time evolution towards the structural balance. The automaton is deterministic and synchronous. The medium is the triangular lattice with some fraction $f$ of links removed. We analyse the number of unbalanced triads (parameterized as `energy), the frequencies of balanced triads and correlations between them. The time evolution enhances the local correlations of balanced triads. Local configurations of unbalanced triads are found which are blinking with period of two time steps.
A cellular automaton model of pulsar glitches is described, based on the superfluid vortex unpinning paradigm. Recent analyses of pulsar glitch data suggest that glitches result from scale-invariant avalanches citep{Melatos07a}, which are consistent with a self-organized critical system (SOCS). A cellular automaton provides a computationally efficient means of modelling the collective behaviour of up to $10^{16}$ vortices in the pulsar interior, whilst ensuring that the dominant aspects of the microphysics are not lost. The automaton generates avalanche distributions that are qualitatively consistent with a SOCS and with glitch data. The probability density functions of glitch sizes and durations are power laws, and the probability density function of waiting times between successive glitches is Poissonian, consistent with statistically independent events. The output of the model depends on the physical and computational paramters used. The fitted power law exponents for the glitch sizes ($a$) and durations ($b$) decreases as the strength of the vortex pinning increases. Similarly the exponents increase as the fraction of vortices that are pinned decreases. For the physical and computational parameters considered, one finds $-4.3leq a leq -2.0$ and $-5.5leq bleq -2.2$, and mean glitching rates in the range $0.0023leqlambdaleq0.13$ in units of inverse time.
Gliders in one-dimensional cellular automata are compact groups of non-quiescent and non-ether patterns (ether represents a periodic background) translating along automaton lattice. They are cellular-automaton analogous of localizations or quasi-local collective excitations travelling in a spatially extended non-linear medium. They can be considered as binary strings or symbols travelling along a one-dimensional ring, interacting with each other and changing their states, or symbolic values, as a result of interactions. We analyse what types of interaction occur between gliders travelling on a cellular automaton `cyclotron and build a catalog of the most common reactions. We demonstrate that collisions between gliders emulate the basic types of interaction that occur between localizations in non-linear media: fusion, elastic collision, and soliton-like collision. Computational outcomes of a swarm of gliders circling on a one-dimensional torus are analysed via implementation of cyclic tag systems.
We present theoretical arguments and simulation data indicating that the scaling of earthquake events in models of faults with long-range stress transfer is composed of at least three distinct regions. These regions correspond to three classes of earthquakes with different underlying physical mechanisms. In addition to the events that exhibit scaling, there are larger ``breakout events that are not on the scaling plot. We discuss the interpretation of these events as fluctuations in the vicinity of a spinodal critical point.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
