Simulation of the Burridge-Knopoff model of earthquakes with variable range stress transfer

Simple models of earthquake faults are important for understanding the mechanisms for their observed behavior, such as Gutenberg-Richter scaling and the relation between large and small events, which is the basis for various forecasting methods. Although cellular automaton models have been studied extensively in the long-range stress transfer limit, this limit has not been studied for the Burridge-Knopoff model, which includes more realistic friction forces and inertia. We find that the latter model with long-range stress transfer exhibits qualitatively different behavior than both the long-range cellular automaton models and the usual Burridge-Knopoff model with nearest-neighbor springs, depending on the nature of the velocity-weakening friction force. These results have important implications for our understanding of earthquakes and other driven dissipative systems.     

Model for the distribution of aftershock interoccurrence times

In this work the distribution of interoccurrence times between earthquakes in aftershock sequences is analyzed and a model based on a nonhomogeneous Poisson (NHP) process is proposed to quantify the observed scaling. In this model the generalized Omori's law for the decay of aftershocks is used as a time-dependent rate in the NHP process. The analytically derived distribution of interoccurrence times is applied to several major aftershock sequences in California to confirm the validity of the proposed hypothesis.     

Ergodic dynamics in a natural threshold system

Numerical simulations suggest that certain driven, dissipative mean-field threshold systems, including earthquake models, can be characterized by statistical properties often associated with ergodic dynamics, in the same sense as stochastic Brownian motion. We applied a fluctuation metric proposed by Thirumalai and Mountain [Phys. Rev. E 47, 479 (1993)]] for statistically stationary systems and find that the natural earthquake fault system in California demonstrates similar ergodic dynamics.     

Nonlinear network dynamics on earthquake fault systems

Earthquake faults occur in interacting networks having emergent space-time modes of behavior not displayed by isolated faults. Using simulations of the major faults in southern California, we find that the physics depends on the elastic interactions among the faults defined by network topology, as well as on the nonlinear physics of stress dissipation arising from friction on the faults. Our results have broad applications to other leaky threshold systems such as integrate-and-fire neural networks.     

Scaling and critical phenomena in a cellular automaton slider-block model for earthquakes

The dynamics of a general class of two-dimensional cellular automaton slider-block models of earthquake faults is studied as a function of the failure rules that determine slip and the nature of the failure threshold. Scaling properties of clusters of failed sites imply the existence of a mean-field spinodal line in systems with spatially random failure thresholds, whereas spatially uniform failure thresholds produce behavior reminiscent of self-organized critical behavior. This model can describe several classes of faults, ranging from those that only exhibit creep to those that produce large events.     

Preconditioned iterative methods for coupled discretizations of fluid flow problems

Computational fluid dynamics, where simulations require largecomputation times, is one of the areas of application of highperformance computing. Schemes such as the SIMPLE (semi-implicitmethod for pressure-linked equations) algorithm are often usedto solve the discrete Navier-Stokes equations. Generally theseschemes take a short time per iteration but require a largenumber of iterations. For simple geometries (or coarser grids)the overall CPU time is small. However, for finer grids or morecomplex geometries the increase in the number of iterationsmay be a drawback and the decoupling of the differential equationsinvolved implies a slow convergence of rotationally dominatedproblems that can be very time consuming for realistic applications.So we analyze here another approach, DIRECTO, that solves theequations in a coupled way. With recent advances in hardwaretechnology and software design, it became possible to solvecoupled Navier-Stokes systems, which are more robust but implyincreasing computational requirements (both in terms of memoryand CPU time). Two approaches are described here (band blockLU factorization and preconditioned GMRES) for the linear solverrequired by the DIRECTO algorithm that solves the fluid flowequations as a coupled system. Comparisons of the effectivenessof incomplete factorization preconditioners applied to the GMRES(generalized minimum residual) method are shown. Some numericalresults are presented showing that it is possible to minimizeconsiderably the CPU time of the coupled approach so that itcan be faster than the decoupled one.     

Preparation, structural analysis and dielectric properties of Bi x La1−x FeO3 perovskite

Rare earth element substituted bismuth ferrites (BiFeO₃) are of enormous importance as magnetoelectric materials. The polycrystalline samples of Bi _x La_1−x FeO₃ (x=0, 0.2, 0.4, 0.6, 0.8) were prepared by solid-state reaction using standard ceramic method. The single-phase formation of these compounds was confirmed by X-ray diffraction (XRD) studies. The samples with x=0, 0.2, 0.4, 0.6 are found to be orthorhombic while the sample with x=0.8 is triclinic. The dielectric constant (ε′) and dissipation factor (tan δ) were measured in the frequency range 100 Hz to 1 MHz at room temperature and as a function of temperature at certain fixed frequencies (1 kHz, 10 kHz, 100 kHz, 1 MHz). All the samples showed dielectric dispersion. The dielectric constant with temperature shows a broad peak; the peak temperature shifts with frequency which reflects the relaxor-type behavior. The peak above 600 K in the measured temperature range corresponds to antiferromagnetic ordering temperature (Néel temperature). The broadness of the peak changes with composition. The ac conductivity as well as ε′ are found to be maximum for the sample x=0.2 at room temperature.