Confined sandpile in two dimensions: Percolation and singular diffusion.

We investigate the properties of a two-state sandpile model subjected to a confining potential in two dimensions. From the microdynamical description, we derive a diffusion equation, and find a stationary solution for the case of a parabolic confining potential. By studying the systems at different confining conditions, we observe two scale-invariant regimes. At a given confining potential strength, the cluster size distribution takes the form of a power law. This regime corresponds to the situation in which the density at the center of the system approaches the critical percolation threshold. The analysis of the fractal dimension of the largest cluster frontier provides evidence that this regime is reminiscent of gradient percolation. By increasing further the confining potential, most of the particles coalesce in a giant cluster, and we observe a regime where the jump size distribution takes the form of a power law. The onset of this second regime is signaled by a maximum in the fluctuation of energy.

Fragmentation processes in two-phase materials.

We investigate the fragmentation process of solid materials with crystalline and amorphous phases using the the discrete element method. Damage initiates inside spherical samples above the contact zone in a region where the circumferential stress field is tensile. Cracks initiated in this region grow to form meridional planes. If the collision energy exceeds a critical value which depends on the material's internal structure, cracks reach the sample surface resulting in fragmentation. We show that this primary fragmentation mechanism is very robust with respect to the internal structure of the material. For all configurations, a sharp transition from the damage to the fragmentation regime is observed, with smaller critical collision energies for crystalline samples. The mass distribution of the fragments follows a power law for small fragments with an exponent that is characteristic for the branching merging process of unstable cracks. Moreover this exponent depends only on the dimensionally of the system and not on the microstructure.

Fragmentation processes in impact of spheres.

We study the brittle fragmentation of spheres by using a three-dimensional discrete element model. Large scale computer simulations are performed with a model that consists of agglomerates of many particles, interconnected by beam-truss elements. We focus on the detailed development of the fragmentation process and study several fragmentation mechanisms. The evolution of meridional cracks is studied in detail. These cracks are found to initiate in the inside of the specimen with quasiperiodic angular distribution. The fragments that are formed when these cracks penetrate the specimen surface give a broad peak in the fragment mass distribution for large fragments that can be fitted by a two-parameter Weibull distribution. This mechanism can only be observed in three-dimensional models or experiments. The results prove to be independent of the degree of disorder in the model. Our results significantly improve the understanding of the fragmentation process for impact fracture since besides reproducing the experimental observations of fragment shapes, impact energy dependence, and mass distribution, we also have full access to the failure conditions and evolution.

Computer simulation of fatigue under diametrical compression.

We study the fatigue fracture of disordered materials by means of computer simulations of a discrete element model. We extend a two-dimensional fracture model to capture the microscopic mechanisms relevant for fatigue and we simulate the diametric compression of a disc shape specimen under a constant external force. The model allows us to follow the development of the fracture process on the macrolevel and microlevel varying the relative influence of the mechanisms of damage accumulation over the load history and healing of microcracks. As a specific example we consider recent experimental results on the fatigue fracture of asphalt. Our numerical simulations show that for intermediate applied loads the lifetime of the specimen presents a power law behavior. Under the effect of healing, more prominent for small loads compared to the tensile strength of the material, the lifetime of the sample increases and a fatigue limit emerges below which no macroscopic failure occurs. The numerical results are in a good qualitative agreement with the experimental findings.

Counterion correlations and attraction between like-charged macromolecules.

A simple model is presented for the appearance of attraction between two like-charged polyions inside a polyelectrolyte solution. The polyions are modeled as rigid cylinders in a continuum dielectric solvent. The strong electrostatic interaction between the polyions and the counterions results in counterion condensation. If the two polyions are sufficiently close to each other their layers of condensed counterions can become correlated resulting in attraction between the macromolecules. To explore the counterion induced attraction we calculate the correlation functions for the condensed counterions. It is found that the correlations are of very short range. For the parameters specific to the double stranded DNA, the correlations and the attraction appear only when the surface-to-surface separation is less than 7 A.