Bi-SOC-states in one-dimensional random cellular automaton.

Two statistically stationary states with power-law scaling of avalanches are found in a simple 1 D cellular automaton. Features of the fixed points, the spiral saddle and the saddle with index 1, are investigated. The migration of states of the automaton between these two self-organized criticality states is demonstrated during evolution of the system in computer simulations. The automaton, being a slowly driven system, can be applied as a toy model of earthquake supercycles.

Effect of mobility ratio on interaction between the fingers in unstable growth processes.

We investigate interactions between thin fingers formed as a result of an instability of an advancing front in growth processes. We show that the fingers can both attract and repel each other, depending on their lengths and the mobility ratio between the invading and displaced phase. To understand the origin of these interactions we introduce a simple resistor model of the fingers. The predictions of the model are then compared to the numerical simulations of two unstable growth processes: dissolution of partially cemented rock fracture and viscous fingering in a regular network of channels.

Network models of dissolution of porous media.

We investigate the chemical dissolution of porous media using a 2D network model in which the system is represented as a series of interconnected pipes with the diameter of each segment increasing in proportion to the local reactant consumption. Moreover, the topology of the network is allowed to change dynamically during the simulation: As the diameters of the eroding pores become comparable with the interpore distances, the pores are joined together, thus changing the interconnections within the network. With this model, we investigate different growth regimes in an evolving porous medium, identifying the mechanisms responsible for the emergence of specific patterns. We consider both the random and regular network and study the effect of the network geometry on the patterns. Finally, we consider practically important problem of finding an optimum flow rate that gives a maximum increase in permeability for a given amount of reactant.