Search PubMed⌕ Search

Biomedical subjects

Ellák Somfai

Publications and source records attributed to Ellák Somfai.

9 recordsLinked to original sources

Scale invariance and universality of force networks in static granular matter.

Force networks form the skeleton of static granular matter. They are the key factor that determines mechanical properties such as stability, elasticity and sound transmission, which are important for civil engineering and industrial processing. Previous studies have focused on investigations of the global structure of external forces (the boundary condition) and on the probability distribution of individual contact forces. So far, however, precise knowledge of the disordered spatial structure of the force network has remained elusive. Here we report that molecular dynamics simulations of realistic granular packings reveal scale invariance of clusters of particles interacting by means of relatively strong forces. Despite visual variation, force networks for various values of the confining pressure and other parameters have identical scaling exponents and scaling function, thereby determining a universality class. Unexpectedly, the flat ensemble of force configurations (a simple generalization of equilibrium statistical mechanics) belongs to this universality class, whereas some widely studied simplified models do not. This implies that the elasticity of the grains and their geometrical disorder do not affect the universal mechanical properties.

Journal Article↗

Anisotropic diffusion limited aggregation in three dimensions: universality and nonuniversality.

We explore the macroscopic consequences of lattice anisotropy for diffusion limited aggregation (DLA) in three dimensions. Simple cubic and bcc lattice growths are shown to approach universal asymptotic states in a coherent fashion, and the approach is accelerated by the use of noise reduction. These states are strikingly anisotropic dendrites with a rich hierarchy of structure. For growth on an fcc lattice, our data suggest at least two stable fixed points of anisotropy, one matching the bcc case. Hexagonal growths, favoring six planar and two polar directions, appear to approach a line of asymptotic states with continuously tunable polar anisotropy. The more planar of these growths visually resembles real snowflake morphologies. Our simulations use a new and dimension-independent implementation of the DLA model. The algorithm maintains a hierarchy of sphere coverings of the growth, supporting efficient random walks onto the growth by spherical moves. Anisotropy was introduced by restricting growth to certain preferred directions.

Journal Article↗

Elastic wave propagation in confined granular systems.

We present numerical simulations of acoustic wave propagation in confined granular systems consisting of particles interacting with the three-dimensional Hertz-Mindlin force law. The response to a short mechanical excitation on one side of the system is found to be a propagating coherent wave front followed by random oscillations made of multiply scattered waves. We find that the coherent wave front is insensitive to details of the packing: force chains do not play an important role in determining this wave front. The coherent wave propagates linearly in time, and its amplitude and width depend as a power law on distance, while its velocity is roughly compatible with the predictions of macroscopic elasticity. As there is at present no theory for the broadening and decay of the coherent wave, we numerically and analytically study pulse propagation in a one-dimensional chain of identical elastic balls. The results for the broadening and decay exponents of this system differ significantly from those of the random packings. In all our simulations, the speed of the coherent wave front scales with pressure as p1/6; we compare this result with experimental data on various granular systems where deviations from the p1/6 behavior are seen. We briefly discuss the eigenmodes of the system and effects of damping are investigated as well.

Journal Article↗

Random walks, diffusion limited aggregation in a wedge, and average conformal maps.

We investigate diffusion-limited aggregation (DLA) in a wedge geometry. Arneodo and collaborators have suggested that the ensemble average of DLA cluster density should be close to the noise-free selected Saffman-Taylor finger. We show that a different, but related, ensemble average, that of the conformal maps associated with random clusters, yields a nontrivial shape which is also not far from the Saffman-Taylor finger. However, we have previously demonstrated that the same average of DLA in a channel geometry is not the Saffman-Taylor finger. This casts doubt on the idea that the average of noisy diffusion-limited growth is governed by a simple transcription of noise-free results.

Cluster Analysis↗

Growth by random walker sampling and scaling of the dielectric breakdown model.

Random walkers absorbing on a boundary sample the harmonic measure linearly and independently: we discuss how the recurrence times between impacts enable nonlinear moments of the measure to be estimated. From this we derive a technique to simulate dielectric breakdown model growth, which is governed nonlinearly by the harmonic measure. For diffusion-limited aggregation, recurrence times are shown to be accurate and effective in probing the multifractal growth measure in its active region. For the dielectric breakdown model our technique grows large clusters efficiently and we are led to significantly revise earlier exponent estimates. Previous results by two conformal mapping techniques were less converged than expected, and in particular a recent theoretical suggestion of superuniversality is firmly refuted.

Journal Article↗

Packing geometry and statistics of force networks in granular media.

The relation between packing geometry and force network statistics is studied for granular media. Based on simulations of two-dimensional packings of Hertzian spheres, we develop a geometrical framework relating the distribution of interparticle forces P(f) to the weight distribution P(w), which is measured in experiments. We apply this framework to reinterpret recent experimental data on strongly deformed packings and suggest that the observed changes of P(w) are dominated by changes in contact network while P(f) remains relatively unaltered. We furthermore investigate the role of packing disorder in the context of the q model and address the question of how force fluctuations build up as a function of the distance beneath the top surface.

Journal Article↗

Diffusion-limited aggregation in channel geometry.

We performed extensive numerical simulation of diffusion-limited aggregation in two-dimensional channel geometry. Contrary to earlier claims, the measured fractal dimension D=1.712+/-0.002 and its leading correction to scaling are the same as in the radial case. The average cluster, defined as the average conformal map, is similar but not identical to Saffman-Taylor fingers.

Journal Article↗

Force and weight distributions in granular media: effects of contact geometry.

The influence of the local contact network on interparticle forces and effective particle weights is studied in simulations of two-dimensional packings of frictionless, Hertzian spheres. The weight distribution P(w) changes qualitatively when approaching a boundary and differs for regular and irregular packings, while the interparticle force distribution P(f) is robust. We provide examples where P(w) at the boundary, which is the quantity probed experimentally, deviates substantially from P(f) in the bulk. Discrepancies between the P(w)'s predicted by the q model and measured in experiments are due to differences in the contact geometry.

Journal Article↗

Off-lattice noise reduction and the ultimate scaling of diffusion-limited aggregation in two dimensions.

Off-lattice diffusion-limited aggregation (DLA) clusters grown with different levels of noise reduction are found to be consistent with a simple fractal fixed point. Cluster shapes and their ensemble variation exhibit a dominant slowest correction to scaling, and this also accounts for the apparent "multiscaling" in the DLA mass distribution. We interpret the correction to scaling in terms of renormalized noise. The limiting value of this variable is strikingly small and is dominated by fluctuations in cluster shape. Earlier claims of anomalous scaling in DLA were misled by the slow approach to this small fixed point value.

Journal Article↗