Catalog footprint

What is connected

103works
25topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

103 published item(s)

preprint2022arXiv

Curvature-induced quantum spin-Hall effect on a Möbius strip

The quantum Hall effect has been predicted and discovered in various condensed-matter systems. A promising quantum material for such topological effects is graphene. We report the numerical observation of a curvature-induced spin-Hall effect in a monolayer graphene Möbius strip. The solution of the Dirac equation on the nontrivial and non-Euclidean manifold reveals that despite the absence of a Hall current, a spin-Hall current is a natural consequence for such a topology, as predicted from symmetry arguments.

preprint2022arXiv

On the scaling of avalanche shape and activity power spectrum in neuronal networks

Many systems in Nature exhibit avalanche dynamics with scale-free features. A general scaling theory has been proposed for critical avalanche profiles in crackling noise, predicting the collapse onto a universal avalanche shape, as well as the scaling behaviour of the activity power spectrum as Brown noise. Recently, much attention has been given to the profile of neuronal avalanches, measured in neuronal systems in vitro and in vivo. Although a universal profile was evidenced, confirming the validity of the general scaling theory, the parallel study of the power spectrum scaling under the same conditions was not performed. The puzzling observation is that in the majority of healthy neuronal systems the power spectrum exhibits a behaviour close to $1/f$, rather than Brown, noise. Here we perform a numerical study of the scaling behaviour of avalanche shape and power spectrum for a model of integrate and fire neurons with a short-term plasticity parameter able to tune the system to criticality. We confirm that, at criticality, the average avalanche size and the avalanche profile fulfill the general avalanche scaling theory. However, the power spectrum consistently exhibits Brown noise behaviour, for both fully excitatory networks and systems with 30\% inhibitory networks. Conversely, a behaviour closer to $1/f$ noise is observed in systems slightly off-criticality. Results suggest that the power spectrum is a good indicator to determine how close neuronal activity is to criticality.

preprint2021arXiv

Phase transitions in a conservative Game of Life

We investigate the dynamics of a conservative version of Conway's Game of Life, in which a pair consisting of a dead and a living cell can switch their states following Conway's rules but only by swapping their positions, irrespective of their mutual distance. Our study is based on square-lattice simulations as well as a mean-field calculation. As the density of dead cells is increased, we identify a discontinuous phase transition between an inactive phase, in which the dynamics freezes after a finite time, and an active phase, in which the dynamics persists indefinitely in the thermodynamic limit. Further increasing the density of dead cells leads the system back to an inactive phase via a second transition, which is continuous on the square lattice but discontinuous in the mean-field limit.

preprint2021arXiv

Self-similar but not conformally invariant traces obtained by modified Loewner forces

The two-dimensional Loewner exploration process is generalized to the case where the random force is self-similar with positively correlated increments. We model this random force by a fractional Brownian motion with Hurst exponent $H\geq \frac{1}{2}\equiv H_{\text{BM}}$, where $H_{\text{BM}}$ stands for the one-dimensional Brownian motion. By manipulating the deterministic force, we design a scale-invariant equation describing self-similar traces which lack conformal invariance. The model is investigated in terms of the "input diffusivity parameter" $κ$, which coincides with the one of the ordinary Schramm-Loewner evolution (SLE) at $H=H_{\text{BM}}$. In our numerical investigation, we focus on the scaling properties of the traces generated for $κ=2,3$, $κ=4$ and $κ=6,8$ as the representatives, respectively, of the dilute phase, the transition point and the dense phase of the ordinary SLE. The resulting traces are shown to be scale-invariant. Using two equivalent schemes, we extract the fractal dimension, $D_f(H)$, of the traces which decrease monotonically with increasing $H$, reaching $D_f=1$ at $H=1$ for all $κ$ values. The left passage probability (LPP) test demonstrates that, for $H$ values not far from the uncorrelated case (small $ε_H\equiv \frac{H-H_{\text{BM}}}{H_{\text{BM}}}$) the prediction of the ordinary SLE is applicable with an effective diffusivity parameter $κ_{\text{eff}}$. Not surprisingly, the $κ_{\text{eff}}$'s do not fulfill the prediction of SLE for the relation between $D_f(H)$ and the diffusivity parameter.

preprint2020arXiv

Simulating isothermal aging of snow

A Monte Carlo algorithm to simulate the isothermal recrystallization process of snow is presented. The snow metamorphism is approximated by two mass redistribution processes, surface diffusion and sublimation-deposition. The algorithm is justified and its parametrization is determined. The simulation results are compared to experimental data, in particular, the temporal evolution of the specific surface area and the ice thickness. We find that the two effects of surface diffusion and sublimation-deposition can accurately model many aspects of the isothermal metamorphism of snow. Furthermore, it is shown that sublimation-deposition is the dominant contribution for temperatures close to the melting point, whereas surface diffusion dominates at temperatures far below the melting point. A simple approximation of gravitational compaction is implemented to simulate density change.

preprint2020arXiv

Spatio-temporal characteristics of dengue outbreaks

After their re-emergence in the last decades, dengue fever and other vector-borne diseases are a potential threat to the lives of millions of people. Based on a data set of dengue cases in the Brazilian city of Fortaleza, collected from 2011 to 2016, we study the spatio-temporal characteristics of dengue outbreaks to characterize epidemic and non-epidemic years. First, we identify regions that show a high prevalence of dengue cases and mosquito larvae in different years and also analyze their corresponding correlations. Our results show that the characteristic correlation length of the epidemic is of the order of the system size, suggesting that factors such as citizen mobility may play a major role as a drive for spatial spreading of vector-borne diseases. Inspired by this observation, we perform a mean-field estimation of the basic reproduction number and find that our estimated values agree well with the values reported for other regions, pointing towards similar underlying spreading mechanisms. These findings provide insights into the spreading characteristics of dengue in densely populated areas and should be of relevance for the design of improved disease containment strategies.

preprint2019arXiv

Landau levels in wrinkled and rippled graphene sheets

We study the discrete energy spectrum of curved graphene sheets in the presence of a magnetic field. The shifting of the Landau levels is determined for complex and realistic geometries of curved graphene sheets. The energy levels follow a similar square root dependence on the energy quantum number as for rippled and flat graphene sheets. The Landau levels are shifted towards lower energies proportionally to the average deformation and the effect is larger compared to a simple uni-axially rippled geometry. Furthermore, the resistivity of wrinkled graphene sheets is calculated for different average space curvatures and shown to obey a linear relation. The study is carried out with a quantum lattice Boltzmann method, solving the Dirac equation on curved manifolds.

preprint2016arXiv

A mesoscopic approach to subcritical fatigue crack growth

We investigate a model for fatigue crack growth in which damage accumulation is assumed to follow a power law of the local stress amplitude, a form which can be generically justified on the grounds of the approximately self-similar aspect of microcrack distributions. Our aim is to determine the relation between model ingredients and the Paris exponent governing subcritical crack-growth dynamics at the macroscopic scale, starting from a single small notch propagating along a fixed line. By a series of analytical and numerical calculations, we show that, in the absence of disorder, there is a critical damage-accumulation exponent $γ$, namely $γ_c=2$, separating two distinct regimes of behavior for the Paris exponent $m$. For $γ>γ_c$, the Paris exponent is shown to assume the value $m=γ$, a result which proves robust against the separate introduction of various modifying ingredients. Explicitly, we deal here with (i) the requirement of a minimum stress for damage to occur; (ii) the presence of disorder in local damage thresholds; (iii) the possibility of crack healing. On the other hand, in the regime $γ<γ_c$ the Paris exponent is seen to be sensitive to the different ingredients added to the model, with rapid healing or a high minimum stress for damage leading to $m=2$ for all $γ<γ_c$, in contrast with the linear dependence $m=6-2γ$ observed for very long characteristic healing times in the absence of a minimum stress for damage. Upon the introduction of disorder on the local fatigue thresholds, which leads to the possible appearance of multiple cracks along the propagation line, the Paris exponent tends to $m\approx 4$ for $γ\lesssim 2$, while retaining the behavior $m=γ$ for $γ\gtrsim 4$.

preprint2016arXiv

Evolution of Strength and Failure of SCC during Early Hydration

The early strength evolution of self-consolidating concrete (SCC) is studied by a set of non-standard mechanical tests for compressive, tensile, shear and bending failure. The results are applicable in an industrial environment for process control, e.g. of slip casting with adaptive molds in robotic fabrication. A procedure for collapsing data to a master evolution curve is presented that allows to distinguish two regimes in the evolution. In the first, the material is capable of undergoing large localized plastic deformation, as expected from thixotropic yield stress fluids. This is followed by a transition to cohesive frictional material behavior dominated by crack growth. The typical differences in tensile and compressive strength of hardened concrete are observed to originate at the transition. Finally, the evolution of a limit surface in principal stress space is constructed and discussed.

preprint2016arXiv

Itinerant conductance in fuse-antifuse networks

We report on a novel dynamic phase in electrical networks, in which current channels perpetually change in time. This occurs when the individual elements of the network are fuse-antifuse devices, namely, become insulators within a certain finite interval of local applied voltages. As a consequence, the macroscopic current exhibits temporal fluctuations which increase with system size. We determine the conditions under which this exotic situation appears by establishing a phase diagram as a function of the applied field and the size of the insulating window. Besides its obvious application as a versatile electronic device, due to its rich variety of behaviors, this network model provides a possible description for particle-laden flow through porous media leading to dynamical clogging and reopening of the local channels in the pore space.

preprint2016arXiv

Micro-mechanical Failure Analysis of Wet Granular Matter

We employ a novel fluid-particle model to study the shearing behavior of granular soils under different saturation levels, ranging from the dry material via the capillary bridge regime to higher saturation levels with percolating clusters. The full complexity of possible liquid morphologies is taken into account, implying the formation of isolated arbitrary-sized liquid clusters with individual Laplace pressures that evolve by liquid exchange via films on the grain surface. Liquid clusters can grow in size, shrink, merge and split, depending on local conditions, changes of accessible liquid and the pore space morphology determined by the granular phase. This phase is represented by a discrete particle model based on Contact Dynamics, where capillary forces exerted from a liquid phase add to the motion of spherical particles. We study the macroscopic response of the system due to an external compression force at various liquid contents with the help of triaxial shear tests. Additionally, the change in liquid cluster distributions during the compression due to the deformation of the pore space is evaluated close to the critical load.

preprint2015arXiv

A percolation model with continuously varying exponents

This work analyzes a percolation model on the diamond hierarchical lattice (DHL), where the percolation transition is retarded by the inclusion of a probability of erasing specific connected structures. It has been inspired by the recent interest on the existence of other universality classes of percolation models. The exact scale invariance and renormalization properties of DHL leads to recurrence maps, from which analytical expressions for the critical exponents and precise numerical results in the limit of very large lattices can be derived. The critical exponents $ν$ and $β$ of the investigated model vary continuously as the erasing probability changes. An adequate choice of the erasing probability leads to the result $ν=\infty$, like in some phase transitions involving vortex formation. The percolation transition is continuous, with $β>0$, but $β$ can be as small as desired. The modified percolation model turns out to be equivalent to the $Q\rightarrow1$ limit of a Potts model with specific long range interactions on the same lattice.

preprint2015arXiv

Continuous Wire Reinforcement for Jammed Granular Architecture

The mechanical behavior of continuous fiber reinforced granular columns is simulated by means of a Discrete Element Model. Spherical particles are randomly deposited simultaneously with a wire, that is deployed following different patterns inside of a flexible cylinder for triaxial compression testing. We quantify the effect of three different fiber deployment patterns on the failure envelope, represented by Mohr-Coulomb cones, and derive suggestions for improved deployment strategies.

preprint2015arXiv

Discrete Element Method simulations of the saturation of aeolian sand transport

The saturation length of aeolian sand transport ($L_s$), characterizing the distance needed by wind-blown sand to adapt to changes in the wind shear, is essential for accurate modeling of the morphodynamics of Earth's sandy landscapes and for explaining the formation and shape of sand dunes. In the last decade, it has become a widely-accepted hypothesis that $L_s$ is proportional to the characteristic distance needed by transported particles to reach the wind speed (the ``drag length''). Here we challenge this hypothesis. From extensive numerical Discrete Element Method simulations, we find that, for medium and strong winds, $L_s\propto V_s^2/g$, where $V_s$ is the saturated value of the average speed of sand particles traveling above the surface and $g$ the gravitational constant. We show that this proportionality is consistent with a recent analytical model, in which the drag length is just one of four similarly important length scales relevant for sand transport saturation.

preprint2015arXiv

Disease-induced resource constraints can trigger explosive epidemics

Advances in mathematical epidemiology have led to a better understanding of the risks posed by epidemic spreading and informed strategies to contain disease spread. However, a challenge that has been overlooked is that, as a disease becomes more prevalent, it can limit the availability of the capital needed to effectively treat those who have fallen ill. Here we use a simple mathematical model to gain insight into the dynamics of an epidemic when the recovery of sick individuals depends on the availability of healing resources that are generated by the healthy population. We find that epidemics spiral out of control into "explosive" spread if the cost of recovery is above a critical cost. This can occur even when the disease would die out without the resource constraint. The onset of explosive epidemics is very sudden, exhibiting a discontinuous transition under very general assumptions. We find analytical expressions for the critical cost and the size of the explosive jump in infection levels in terms of the parameters that characterize the spreading process. Our model and results apply beyond epidemics to contagion dynamics that self-induce constraints on recovery, thereby amplifying the spreading process.

preprint2015arXiv

Fractal frontiers of bursts and cracks in a fiber bundle model of creep rupture

We investigate the geometrical structure of breaking bursts generated during the creep rupture of heterogeneous materials. Based on a fiber bundle model with localized load sharing we show that bursts are compact geometrical objects, however, their external frontier has a fractal structure which reflects their growth dynamics. The perimeter fractal dimension of bursts proved to have the universal value 1.25 independent of the external load and of the amount of disorder in the system. We conjecture that according to their geometrical features breaking bursts fall in the universality class of loop-erased self-avoiding random walks with perimeter fractal dimension 5/4 similar to the avalanches of Abelian sand pile models. The fractal dimension of the growing crack front along which bursts occur proved to increase from 1 to 1.25 as bursts gradually cover the entire front.

preprint2015arXiv

Generating Robust and Efficient Networks Under Targeted Attacks

Much of our commerce and traveling depend on the efficient operation of large scale networks. Some of those, such as electric power grids, transportation systems, communication networks, and others, must maintain their efficiency even after several failures, or malicious attacks. We outline a procedure that modifies any given network to enhance its robustness, defined as the size of its largest connected component after a succession of attacks, whilst keeping a high efficiency, described in terms of the shortest paths among nodes. We also show that this generated set of networks is very similar to networks optimized for robustness in several aspects such as high assortativity and the presence of an onion-like structure.

preprint2015arXiv

Grain Scale Modeling of Arbitrary Fluid Saturation in Random Packings

We propose a model for increasing liquid saturation in a granular packing which can account for liquid redistribution at saturation levels beyond the well-studied capillary bridge regime. The model is capable of resolving and combining capillary bridges, menisci and fully saturated pores to form local liquid clusters of any shape. They can exchange volume due to the local Laplace pressure gradient via a liquid film on the surfaces of grains. Local instabilities like Haines jumps trigger the discontinuous evolution of the liquid front. The applicability of the model is demonstrated and compared to benchmark experiments on the level of individual liquid structures as well as on larger systems.

preprint2015arXiv

Hierarchical structures for a robustness-oriented capacity design

In this paper, we study the response of 2D framed structures made of rectangular cells, to the sudden removal of columns. We employ a simulation algorithm based on the Discrete Element Method, where the structural elements are represented by elasto-plastic Euler Bernoulli beams with elongation-rotation failure threshold. The effect of structural cell slenderness and of topological hierarchy on the dynamic residual strength after damage $\ROne$ is investigated. Topologically \textit{hierarchical} frames have a primary structure made of few massive elements, while \textit{homogeneous} frames are made of many thin elements. We also show how $\ROne$ depends on the activated collapse mechanisms, which are determined by the mechanical hierarchy between beams and columns, i.e. by their relative strength and stiffness. Finally, principles of robustness-oriented capacity design which seem to be in contrast to the conventional anti-seismic capacity design are addressed.

preprint2015arXiv

Invasion Percolation with a Hardening Interface under Gravity

We propose a modified Invasion Percolation (IP) model to simulate the infiltration of glue into a porous medium under gravity in 2D. Initially, the medium is saturated with air and then invaded by a fluid that has a hardening effect taking place from the interface towards the interior by contact with the air. To take into account that interfacial hardening, we use an IP model where capillary pressures of the growth sites are increased with time. In our model, if a site stays for a certain time at interface, it becomes a hard site and cannot be invaded anymore. That represents the glue interface becoming hard due to exposition with the air. Buoyancy forces are included in this system through the Bond number which represents the competition between the hydrostatic and capillary forces. We then compare our results with results from literature of non-hardening fluids in each regime of Bond number. We see that the invasion patterns change strongly with hardening while the non-hardening behavior remains basically not affected.

preprint2015arXiv

Inverse determination of effective mechanical properties of adhesive bondlines

A new approach for determining effective mechanical bondline properties using a combined experimental-numerical modal analysis technique is proposed. After characterizing clear spruce wood boards, an adhesive layer is applied to the boards surfaces. The shift of the eigenfrequencies resulting from the adhesive layer together with information on the bondline geometry can then be used to inversely determine the mechanical properties of the adhesive layer using Finite Element Models. The calculated values for clear wood as well as for the adhesive layer lie within reasonable ranges, thus demonstrating the methods potential.

preprint2015arXiv

Lattice-Boltzmann Method for Geophysical Plastic Flows

We explore possible applications of the Lattice-Boltzmann Method for the simulation of geophysical flows. This fluid solver, while successful in other fields, is still rarely used for geotechnical applications. We show how the standard method can be modified to represent free-surface realization of mudflows, debris flows, and in general any plastic flow, through the implementation of a Bingham constitutive model. The chapter is completed by an example of a full-scale simulation of a plastic fluid flowing down an inclined channel and depositing on a flat surface. An application is given, where the fluid interacts with a vertical obstacle in the channel.

preprint2015arXiv

Mechanisms in Impact Fragmentation

The brittle fragmentation of spheres is studied numerically by a 3D Discrete Element Model. Large scale computer simulations are performed with models that consist of agglomerates of many spherical particles, interconnected by beam-truss elements. We focus on a detailed description of the fragmentation process and study several fragmentation mechanisms involved. The evolution of meridional cracks is studied in detail. These cracks are found to initiate in the inside of the specimen with quasi-periodic angular distribution and give a broad peak in the fragment mass distribution for large fragments that can be fitted by a two-parameter Weibull distribution. The results prove to be independent of the degree of disorder in the model, but mean fragment sizes scale with velocity. Our results reproduce many experimental observations of fragment shapes, impact energy dependence or mass distribution, and significantly improve the understanding of the fragmentation process for impact fracture since we have full access to the failure conditions and evolution.

preprint2015arXiv

Midair collisions enhance saltation

Here we address the old question in Aeolian particle transport about the role of midair collisions. We find that, surprisingly, these collisions do enhance the overall flux substantially. The effect depends strongly on restitution coefficient and wind speed. We can explain this observation as a consequence of a soft-bed of grains which floats above the ground and reflects the highest flying particles. We make the unexpected observation that the flux is maximized at an intermediate restitution coefficient of about 0.7, which is comparable to values experimentally measured for collisions between sand grains.

preprint2015arXiv

Moisture-induced Damage Evolution in Laminated Beech

Structural elements made of laminated hardwood are increasingly used in timber engineering. In this combined numerical and experimental approach, damage onset and propagation in uni-directional and cross-laminated samples out of European beech due to climatic changes are studied. The inter- and intra-laminar damage evolution is characterized for various configurations adhesively bonded by three structural adhesive systems. Typical situations are simulated by means of a comprehensive moisture-dependent non-linear rheological finite element model for wood with the capability to capture delaminations. The simulations give insight into the role of different strain components such as visco-elastic, mechano-sorptive, plastic, and hygro-elastic deformations under changing moisture content in progressive damage and delamination. We show the stress buildup under cyclic hygric loading resulting in hygro-fatigue and modify an analytical micro-mechanics of damage model, originally developed for cross-ply laminates, to describe the problem of moisture-induced damage in beech lamellae.

preprint2015arXiv

Optimal percentage of inhibitory synapses in multi-task learning

Performing more tasks in parallel is a typical feature of complex brains. These are characterized by the coexistence of excitatory and inhibitory synapses, whose percentage in mammals is measured to have a typical value of 20-30\%. Here we investigate parallel learning of more Boolean rules in neuronal networks. We find that multi-task learning results from the alternation of learning and forgetting of the individual rules. Interestingly, a fraction of 30\% inhibitory synapses optimizes the overall performance, carving a complex backbone supporting information transmission with a minimal shortest path length. We show that 30\% inhibitory synapses is the percentage maximizing the learning performance since it guarantees, at the same time, the network excitability necessary to express the response and the variability required to confine the employment of resources.

preprint2015arXiv

Orthotropic rotation-free thin shell elements

A method to simulate orthotropic behaviour in thin shell finite elements is proposed. The approach is based on the transformation of shape function derivatives, resulting in a new orthogonal basis aligned to a specified preferred direction for all elements. This transformation is carried out solely in the undeformed state leaving minimal additional impact on the computational effort expended to simulate orthotropic materials compared to isotropic, resulting in a straightforward and highly efficient implementation. This method is implemented for rotation-free triangular shells using the finite element framework built on the Kirchhoff--Love theory employing subdivision surfaces. The accuracy of this approach is demonstrated using the deformation of a pinched hemispherical shell (with a 18° hole) standard benchmark. To showcase the efficiency of this implementation, the wrinkling of orthotropic sheets under shear displacement is analyzed. It is found that orthotropic subdivision shells are able to capture the wrinkling behavior of sheets accurately for coarse meshes without the use of an additional wrinkling model.

preprint2015arXiv

Packing of elastic wires in flexible shells

The packing problem of long thin filaments that are injected into confined spaces is of fundamental interest for physicists and biologists alike. How linear threads pack and coil is well known only for the ideal case of rigid containers, though. Here, we force long elastic rods into flexible spatial confinement borne by an elastic shell to examine under which conditions recently acquired knowledge on wire packing in rigid spheres breaks down. We find that unlike in rigid cavities, friction plays a key role by giving rise to the emergence of two distinct packing patterns. At low friction, the wire densely coils into an ordered toroidal bundle with semi-ellipsoidal cross-section, while at high friction, it packs into a highly disordered, hierarchic structure. These two morphologies are shown to be separated by a continuous phase transition. Our findings demonstrate the dramatic impact of friction and confinement elasticity on filamentous packing and might drive future research on such systems in physics, biology and even medical technology toward including these mutually interacting effects.

preprint2015arXiv

Pore space analysis of beech wood: the vessel network

Water transport in wood is vital for the survival of trees. With synchrotron radiation X-ray tomographic microscopy (SRXTM), it becomes possible to characterize and quantify the 3D network formed by vessels that are responsible for longitudinal transport. In the present paper, the spatial size dependence of vessels and the organization inside single growth rings in terms of vessel induced porosity was studied by SRXTM. Network characteristics, such as connectivity, were deduced by digital image analysis from the processed tomographic data and related to known complex network topologies.

preprint2015arXiv

Progressive Collapse Mechanisms of Brittle and Ductile Framed Structures

In this paper, we study the progressive collapse of 3D framed structures made of reinforced concrete after the sudden loss of a column. The structures are represented by elasto-plastic Euler Bernoulli beams with elongation-rotation failure threshold. We performed simulations using the Discrete Element Method considering inelastic collisions between the structural elements. The results show what collapse initiation and impact-driven propagation mechanisms are activated in structures with different geometric and mechanical features. Namely, we investigate the influence of the cross sectional size and reinforcement $α$ and of the plastic capacity $β$ of the structural elements. We also study the final collapse extent and the fragment size distribution and their relation to $α$, $β$ and to the observed collapse mechanisms. Finally, we compare the damage response of structures with symmetric and asymmetric reinforcement in the beams.

preprint2015arXiv

Simulating Thin Sheets: Buckling, Wrinkling, Folding and Growth

Numerical simulations of thin sheets undergoing large deformations are computationally challenging. Depending on the scenario, they may spontaneously buckle, wrinkle, fold, or crumple. Nature's thin tissues often experience significant anisotropic growth, which can act as the driving force for such instabilities. We use a recently developed finite element model to simulate the rich variety of nonlinear responses of Kirchhoff-Love sheets. The model uses subdivision surface shape functions in order to guarantee convergence of the method, and to allow a finite element description of anisotropically growing sheets in the classical Rayleigh-Ritz formalism. We illustrate the great potential in this approach by simulating the inflation of airbags, the buckling of a stretched cylinder, as well as the formation and scaling of wrinkles at free boundaries of growing sheets. Finally, we compare the folding of spatially confined sheets subject to growth and shrinking confinement to find that the two processes are equivalent.

preprint2015arXiv

Study on the fragmentation of shells

Fragmentation can be observed in nature and in everyday life on a wide range of length scales and for all kinds of technical applications. Most studies on dynamic failure focus on the behaviour of bulk systems in one, two and three dimensions under impact and explosive loading, showing universal power law behaviour of fragment size distribution. However, hardly any studies have been devoted to fragmentation of shells. We present a detailed experimental and theoretical study on the fragmentation of closed thin shells of various materials, due to an excess load inside the system and impact with a hard wall. Characteristic fragmentation mechanisms are identified by means of a high speed camera and fragment shapes and mass distributions are evaluated. Theoretical rationalisation is given by means of stochastic break-up models and large-scale discrete element simulations with spherical shell systems under different extreme loading situations. By this we ex-plain fragment shapes and distributions and prove a power law for the fragment mass distribution. Satisfactory agreement between experimental findings and nu-merical predictions of the exponents of the power laws for the fragment shapes is obtained.

preprint2015arXiv

The Role of Contact Angle Hysteresis for Fluid Transport in Wet Granular Matter

The stability of sand castles is determined by the structure of wet granulates. Experimental data about the size distribution of fluid pockets are ambiguous about their origin. We discovered that contact angle hysteresis plays a fundamental role in the equilibrium distribution of bridge volumes, and not geometrical disorder as commonly conjectured, which has substantial consequences on the mechanical properties of wet granular beds, including a history dependent rheology and lowered strength. Our findings are obtained using a novel model where the Laplace pressures, bridge volumes and contact angles are dynamical variables associated to the contact points. While accounting for contact line pinning, we track the temporal evolution of each bridge. We observe a cross-over to a power-law decay of the variance of capillary pressures at late times and a saturation of the variance of bridge volumes to a finite value connected to contact line pinning. Large scale simulations of liquid transport in the bridge network reveal that the equilibration dynamics at early times is well described by a mean field model. The spread of final bridge volumes can be directly related to the magnitude of contact angle hysteresis.

preprint2015arXiv

Why do Particle Clouds Generate Electric Charges?

Grains in desert sandstorms spontaneously generate strong electrical charges; likewise volcanic dust plumes produce spectacular lightning displays. Charged particle clouds also cause devastating explosions in food, drug and coal processing industries. Despite the wide-ranging importance of granular charging in both nature and industry, even the simplest aspects of its causes remain elusive, because it is difficult to understand how inert grains in contact with little more than other inert grains can generate the large charges observed. Here, we present a simple yet predictive explanation for the charging of granular materials in collisional flows. We argue from very basic considerations that charge transfer can be expected in collisions of identical dielectric grains in the presence of an electric field, and we confirm the model's predictions using discrete-element simulations and a tabletop granular experiment.

preprint2014arXiv

Crazy 1d billiards: Behavior of spring-fixated, noisy colliding particles

We study a one-dimensional system of spatially extended particles, which are fixated to regularly spaced locations by means of elastic springs. The particles are assumed to be driven by a Gaussian noise and to have dissipative, energy-conserving or anti-dissipative (flipper-like) interactions, when the particle density exceeds a critical threshold. While each particle in separation shows a wellbehaved behavior characterized by a Gaussian velocity distribution, the interaction of particles at high densities can cause an avalanche-like momentum and energy transfer, which can generate steep power laws without a well-defined variance and mean value. Specifically, the velocity variance increases dramatically towards the free boundaries of the driven-many-particle system. The model might also have some relevance for a better understanding of crowd disasters. Our results suggest that these are most likely caused by passive momentum transfers, and not by active pushing.

preprint2014arXiv

Explosive metallic breakdown

We investigate the metallic breakdown of a substrate on which highly conducting particles are adsorbed and desorbed with a probability that depends on the local electric field. We find that, by tuning the relative strength $q$ of this dependence, the breakdown can change from continuous to explosive. Precisely, in the limit in which the adsorption probability is the same for any finite voltage drop, we can map our model exactly onto the $q$-state Potts model and thus the transition to a jump occurs at $q=4$. In another limit, where the adsorption probability becomes independent of the local field strength, the traditional bond percolation model is recovered. Our model is thus an example of a possible experimental realization exhibiting a truly discontinuous percolation transition.

preprint2014arXiv

Fragmentation processes in two-phase materials

We investigate the fragmentation process of solid materials with crystalline and amorphous phases using 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 micro structure.

preprint2014arXiv

Morphogenesis of filaments growing in flexible confinements

Space-saving design is a requirement that is encountered in biological systems and the development of modern technological devices alike. Many living organisms dynamically pack their polymer chains, filaments or membranes inside of deformable vesicles or soft tissue like cell walls, chorions, and buds. Surprisingly little is known about morphogenesis due to growth in flexible confinements - perhaps owing to the daunting complexity lying in the nonlinear feedback between packed material and expandable cavity. Here we show by experiments and simulations how geometric and material properties lead to a plethora of morphologies when elastic filaments are growing far beyond the equilibrium size of a flexible thin sheet they are confined in. Depending on friction, sheet flexibility and thickness, we identify four distinct morphological phases emerging from bifurcation and present the corresponding phase diagram. Four order parameters quantifying the transitions between these phases are proposed.

preprint2014arXiv

Multi-scale approach to invasion percolation of rock fracture networks

A multi-scale scheme for the invasion percolation of rock fracture networks with heterogeneous fracture aperture fields is proposed. Inside fractures, fluid transport is calculated on the finest scale and found to be localized in channels as a consequence of the aperture field. The channel network is characterized and reduced to a vectorized artificial channel network (ACN). Different realizations of ACNs are used to systematically calculate efficient apertures for fluid transport inside differently sized fractures as well as fracture intersection and entry properties. Typical situations in fracture networks are parameterized by fracture inclination, flow path length along the fracture and intersection lengths in the entrance and outlet zones of fractures. Using these scaling relations obtained from the finer scales, we simulate the invasion process of immiscible fluids into saturated discrete fracture networks, which were studied in previous works.

preprint2014arXiv

Particle-fluid-structure interaction for debris flow impact on flexible barriers

Flexible barriers are increasingly used for the protection from debris flow in mountainous terrain due to their low cost and environmental impact. However, a numerical tool for rational design of such structures is still missing. In this work, a hybrid computational framework is presented, using a total Lagrangian formulation of the Finite Element Method (FEM) to represent a flexible barrier. The actions exerted on the structure by a debris flow are obtained from simultaneous simulations of the flow of a fluid-grain mixture, using two conveniently coupled solvers: the Discrete Element Method (DEM) governs the motion of the grains, while the free-surface non-Newtonian fluid phase is solved using the Lattice-Boltzmann Method (LBM). Simulations on realistic geometries show the dependence of the momentum transfer on the barrier on the composition of the debris flow, challenging typical assumptions made during the design process today. In particular, we demonstrate that both grains and fluid contribute in a non-negligible way to the momentum transfer. Moreover, we show how the flexibility of the barrier reduces its vulnerability to structural collapse, and how the stress is distributed on its fabric, highlighting potential weak points.

preprint2014arXiv

Rheological Model for Wood

Wood as the most important natural and renewable building material plays an important role in the construction sector. Nevertheless, its hygroscopic character basically affects all related mechanical properties leading to degradation of material stiffness and strength over the service life. Accordingly, to attain reliable design of the timber structures, the influence of moisture evolution and the role of time- and moisture-dependent behaviors have to be taken into account. For this purpose, in the current study a 3D orthotropic elasto-plastic, visco-elastic, mechano-sorptive constitutive model for wood, with all material constants being defined as a function of moisture content, is presented. The corresponding numerical integration approach, with additive decomposition of the total strain is developed and implemented within the framework of the finite element method (FEM). Moreover to preserve a quadratic rate of asymptotic convergence the consistent tangent operator for the whole model is derived. Functionality and capability of the presented material model are evaluated by performing several numerical verification simulations of wood components under different combinations of mechanical loading and moisture variation. Additionally, the flexibility and universality of the introduced model to predict the mechanical behavior of different species are demonstrated by the analysis of a hybrid wood element. Furthermore, the proposed numerical approach is validated by comparisons of computational evaluations with experimental results.

preprint2014arXiv

Shock waves on complex networks

Power grids, road maps, and river streams are examples of infrastructural networks which are highly vulnerable to external perturbations. An abrupt local change of load (voltage, traffic density, or water level) might propagate in a cascading way and affect a significant fraction of the network. Almost discontinuous perturbations can be modeled by shock waves which can eventually interfere constructively and endanger the normal functionality of the infrastructure. We study their dynamics by solving the Burgers equation under random perturbations on several real and artificial directed graphs. Even for graphs with a narrow distribution of node properties (e.g., degree or betweenness), a steady state is reached exhibiting a heterogeneous load distribution, having a difference of one order of magnitude between the highest and average loads. Unexpectedly we find for the European power grid and for finite Watts-Strogatz networks a broad pronounced bimodal distribution for the loads. To identify the most vulnerable nodes, we introduce the concept of node-basin size, a purely topological property which we show to be strongly correlated to the average load of a node.

preprint2013arXiv

Efficient algorithm to study interconnected networks

Interconnected networks have been shown to be much more vulnerable to random and targeted failures than isolated ones, raising several interesting questions regarding the identification and mitigation of their risk. The paradigm to address these questions is the percolation model, where the resilience of the system is quantified by the dependence of the size of the largest cluster on the number of failures. Numerically, the major challenge is the identification of this cluster and the calculation of its size. Here, we propose an efficient algorithm to tackle this problem. We show that the algorithm scales as O(N log N), where N is the number of nodes in the network, a significant improvement compared to O(N^2) for a greedy algorithm, what permits studying much larger networks. Our new strategy can be applied to any network topology and distribution of interdependencies, as well as any sequence of failures.

preprint2013arXiv

Finite-size scaling investigation of the liquid-liquid critical point in ST2 water and its stability with respect to crystallization

The liquid-liquid critical point scenario of water hypothesizes the existence of two metastable liquid phases---low-density liquid (LDL) and high-density liquid (HDL)---deep within the supercooled region. The hypothesis originates from computer simulations of the ST2 water model, but the stability of the LDL phase with respect to the crystal is still being debated. We simulate supercooled ST2 water at constant pressure, constant temperature and constant number of molecules N for N<=729 and times up to 1000 ns. We observe clear differences between the two liquids, both structural and dynamical. Using several methods, including finite-size scaling, we confirm the presence of a liquid-liquid phase transition ending in a critical point. We find that the LDL is stable with respect to the crystal in 98% of our runs (we perform 372 runs for LDL or LDL-like states), and in 100% of our runs for the two largest system sizes (N=512 and 729, for which we perform 136 runs for LDL or LDL-like states). In all these runs tiny crystallites grow and then melt within 1000 ns. Only for N<=343 we observe six events (over 236 runs for LDL or LDL-like states) of spontaneous crystallization after crystallites reach an estimated critical size of about 70+/-10 molecules.

preprint2013arXiv

Flux saturation length of sediment transport

Sediment transport along the surface drives geophysical phenomena as diverse as wind erosion and dune formation. The main length-scale controlling the dynamics of sediment erosion and deposition is the saturation length $L_\mathrm{s}$, which characterizes the flux response to a change in transport conditions. Here we derive, for the first time, an expression predicting $L_\mathrm{s}$ as a function of the average sediment velocity under different physical environments. Our expression accounts for both the characteristics of sediment entrainment and the saturation of particle and fluid velocities, and has only two physical parameters which can be estimated directly from independent experiments. We show that our expression is consistent with measurements of $L_\mathrm{s}$ in both aeolian and subaqueous transport regimes over at least five orders of magnitude in the ratio of fluid and particle density, including on Mars.

preprint2013arXiv

Microscopic mechanism for the shear-thickening of non-Brownian suspensions

We propose a simple model, supported by contact-dynamics simulations as well as rheology and friction measurements, that links the transition from continuous to discontinuous shear-thickening in dense granular pastes to distinct lubrication regimes in the particle contacts. We identify a local Sommerfeld number that determines the transition from Newtonian to shear-thickening flows, and then show that the suspension's volume fraction and the boundary lubrication friction coefficient control the nature of the shear-thickening transition, both in simulations and experiments.

preprint2013arXiv

Numerical modeling of the wind flow over a transverse dune

Transverse dunes, which form under unidirectional winds and have fixed profile in the direction perpendicular to the wind, occur on all celestial objects of our solar system where dunes have been detected. Here we perform a numerical study of the average turbulent wind flow over a transverse dune by means of computational fluid dynamics simulations. We find that the length of the zone of recirculating flow at the dune lee --- the {\em{separation bubble}} --- displays a surprisingly strong dependence on the wind shear velocity, $u_{\ast}$: it is nearly independent of $u_{\ast}$ for shear velocities within the range between $0.2\,$m$$s and $0.8\,$m$$s but increases linearly with $u_{\ast}$ for larger shear velocities. Our calculations show that transport in the direction opposite to dune migration within the separation bubble can be sustained if $u_{\ast}$ is larger than approximately $0.39\,$m$$s, whereas a larger value of $u_{\ast}$ (about $0.49\,$m$$s) is required to initiate this reverse transport.

preprint2013arXiv

Origins of barchan dune asymmetry: insights from numerical simulations

Barchan dunes --- crescent-shaped dunes that form in areas of unidirectional winds and low sand availability --- commonly display an asymmetric shape, with one limb extended downwind. Several factors have been identified as potential causes for barchan dune asymmetry on Earth and Mars: asymmetric bimodal wind regime, topography, influx asymmetry and dune collision. However, the dynamics and potential range of barchan morphologies emerging under each specific scenario that leads to dune asymmetry are far from being understood. In the present work, we use dune modeling in order to investigate the formation and evolution of asymmetric barchans. We find that a bimodal wind regime causes limb extension when the divergence angle between primary and secondary winds is larger than $90^{\circ}$, whereas the extended limb evolves into a seif dune if the ratio between secondary and primary transport rates is larger than 25%. Calculations of dune formation on an inclined surface under constant wind direction also lead to barchan asymmetry, however no seif dune is obtained from surface tilting alone. Asymmetric barchans migrating along a tilted surface move laterally, with transverse migration velocity proportional to the slope of the terrain. Limb elongation induced by topography can occur when a barchan crosses a topographic rise. Furthermore, transient asymmetric barchan shapes with extended limb also emerge during collisions between dunes or due to an asymmetric influx. Our findings can be useful for making quantitative inference on local wind regimes or spatial heterogeneities in transport conditions of planetary dune fields hosting asymmetric barchans.

preprint2013arXiv

Subdivision Shell Elements with Anisotropic Growth

A thin shell finite element approach based on Loop's subdivision surfaces is proposed, capable of dealing with large deformations and anisotropic growth. To this end, the Kirchhoff-Love theory of thin shells is derived and extended to allow for arbitrary in-plane growth. The simplicity and computational efficiency of the subdivision thin shell elements is outstanding, which is demonstrated on a few standard loading benchmarks. With this powerful tool at hand, we demonstrate the broad range of possible applications by numerical solution of several growth scenarios, ranging from the uniform growth of a sphere, to boundary instabilities induced by large anisotropic growth. Finally, it is shown that the problem of a slowly and uniformly growing sheet confined in a fixed hollow sphere is equivalent to the inverse process where a sheet of fixed size is slowly crumpled in a shrinking hollow sphere in the frictionless, quasi-static, elastic limit.

preprint2013arXiv

Towards designing robust coupled networks

Natural and technological interdependent systems have been shown to be highly vulnerable due to cascading failures and an abrupt collapse of global connectivity under initial failure. Mitigating the risk by partial disconnection endangers their functionality. Here we propose a systematic strategy of selecting a minimum number of autonomous nodes that guarantee a smooth transition in robustness. Our method which is based on betweenness is tested on various examples including the famous 2003 electrical blackout of Italy. We show that, with this strategy, the necessary number of autonomous nodes can be reduced by a factor of five compared to a random choice. We also find that the transition to abrupt collapse follows tricritical scaling characterized by a set of exponents which is independent on the protection strategy.

preprint2012arXiv

Activity-dependent neuronal model on complex networks

Neuronal avalanches are a novel mode of activity in neuronal networks, experimentally found in vitro and in vivo, and exhibit a robust critical behaviour: These avalanches are characterized by a power law distribution for the size and duration, features found in other problems in the context of the physics of complex systems. We present a recent model inspired in self-organized criticality, which consists of an electrical network with threshold firing, refractory period and activity-dependent synaptic plasticity. The model reproduces the critical behaviour of the distribution of avalanche sizes and durations measured experimentally. Moreover, the power spectra of the electrical signal reproduce very robustly the power law behaviour found in human electroencephalogram (EEG) spectra. We implement this model on a variety of complex networks, i.e. regular, small-world and scale-free and verify the robustness of the critical behaviour.

preprint2012arXiv

Conductivity of Coniglio-Klein clusters

We performed numerical simulations of the $q$-state Potts model to compute the reduced conductivity exponent $t/ ν$ for the critical Coniglio-Klein clusters in two dimensions, for values of $q$ in the range $[1;4]$. At criticality, at least for $q<4$, the conductivity scales as $C(L) \sim L^{-\frac{t}ν}$, where $t$ and $ν$ are, respectively, the conductivity and correlation length exponents. For q=1, 2, 3, and 4, we followed two independent procedures to estimate $t / ν$. First, we computed directly the conductivity at criticality and obtained $t / ν$ from the size dependence. Second, using the relation between conductivity and transport properties, we obtained $t / ν$ from the diffusion of a random walk on the backbone of the cluster. From both methods, we estimated $t / ν$ to be $0.986 \pm 0.012$, $0.877 \pm 0.014$, $0.785 \pm 0.015$, and $0.658 \pm 0.030$, for q=1, 2, 3, and 4, respectively. We also evaluated $t /ν$ for non integer values of $q$ and propose the following conjecture $40gt/ ν=72+20g-3g^2$ for the dependence of the reduced conductivity exponent on $q$, in the range $ 0 \leq q \leq 4$, where $g$ is the Coulomb gas coupling.

preprint2012arXiv

Fatigue in disordered media

We obtain the Paris law of fatigue crack propagation in a disordered solid using a fuse network model where the accumulated damage in each resistor increases with time as a power law of the local current amplitude. When a resistor reaches its fatigue threshold, it burns irreversibly. Over time, this drives cracks to grow until the system is fractured in two parts. We study the relation between the macroscopic exponent of the crack growth rate -- entering the phenomenological Paris law -- and the microscopic damage-accumulation exponent, $γ$, under the influence of disorder. The way the jumps of the growing crack, $Δa$, and the waiting-time between successive breaks, $Δt$, depend on the type of material, via $γ$, are also investigated. We find that the averages of these quantities, $<Δa>$ and $<Δt>/<t_r>$, scale as power laws of the crack length $a$, $<Δa> \propto a^α$ and $<Δt>/<t_r> \propto a^{-β}$, where $<t_r>$ is the average rupture time. Strikingly, our results show, for small values of $γ$, a decrease in the exponent of the Paris law in comparison with the homogeneous case, leading to an increase in the lifetime of breaking materials. For the particular case of $γ=0$, when fatigue is exclusively ruled by disorder, an analytical treatment confirms the results obtained by simulation.

preprint2012arXiv

Finite Element Simulation of Dense Wire Packings

A finite element program is presented to simulate the process of packing and coiling elastic wires in two- and three-dimensional confining cavities. The wire is represented by third order beam elements and embedded into a corotational formulation to capture the geometric nonlinearity resulting from large rotations and deformations. The hyperbolic equations of motion are integrated in time using two different integration methods from the Newmark family: an implicit iterative Newton-Raphson line search solver, and an explicit predictor-corrector scheme, both with adaptive time stepping. These two approaches reveal fundamentally different suitability for the problem of strongly self-interacting bodies found in densely packed cavities. Generalizing the spherical confinement symmetry investigated in recent studies, the packing of a wire in hard ellipsoidal cavities is simulated in the frictionless elastic limit. Evidence is given that packings in oblate spheroids and scalene ellipsoids are energetically preferred to spheres.

preprint2012arXiv

Fluid depletion in shear bands

How does pore liquid reconfigure within shear bands in wet granular media? Conventional wisdom predicts that liquid is drawn into dilating granular media. We, however, find a depletion of liquid in shear bands despite increased porosity due to dilatancy. This apparent paradox is resolved by a microscale model for liquid transport at low liquid contents induced by rupture and reconfiguration of individual liquid bridges. Measured liquid content profiles show macroscopic depletion bands similar to results of numerical simulations. We derive a modified diffusion description for rupture-induced liquid migration.

preprint2012arXiv

How dense can one pack spheres of arbitrary size distribution?

We present the first systematic algorithm to estimate the maximum packing density of spheres when the grain sizes are drawn from an arbitrary size distribution. With an Apollonian filling rule, we implement our technique for disks in 2d and spheres in 3d. As expected, the densest packing is achieved with power-law size distributions. We also test the method on homogeneous and on empirical real distributions, and we propose a scheme to obtain experimentally accessible distributions of grain sizes with low porosity. Our method should be helpful in the development of ultra-strong ceramics and high performance concrete.

preprint2012arXiv

Linear stability analysis of transverse dunes

Sand-moving winds blowing from a constant direction in an area of high sand availability form transverse dunes, which have a fixed profile in the direction orthogonal to the wind. Here we show, by means of a linear stability analysis, that transverse dunes are intrinsically unstable. Any along-axis perturbation on a transverse dune amplify in the course of dune migration due to the combined effect of two main factors, namely: the lateral transport through avalanches along the dune's slip-face, and the scaling of dune migration velocity with the inverse of the dune height. Our calculations provide a quantitative explanation for recent observations from experiments and numerical simulations, which showed that transverse dunes moving on the bedrock cannot exist in a stable form and decay into a chain of crescent-shaped barchans.

preprint2012arXiv

Liquid migration in sheared unsaturated granular media

We show how liquid migrates in sheared unsaturated granular media using a grain scale model for capillary bridges. Liquid is redistributed to neighboring contacts after rupture of individual capillary bridges leading to redistribution of liquid on large scales. The liquid profile evolution coincides with a recently developed continuum description for liquid migration in shear bands. The velocity profiles which are linked to the migration of liquid as well as the density profiles of wet and dry granular media are studied.

preprint2012arXiv

Optimal box-covering algorithm for fractal dimension of complex networks

The self-similarity of complex networks is typically investigated through computational algorithms the primary task of which is to cover the structure with a minimal number of boxes. Here we introduce a box-covering algorithm that not only outperforms previous ones, but also finds optimal solutions. For the two benchmark cases tested, namely, the E. Coli and the WWW networks, our results show that the improvement can be rather substantial, reaching up to 15% in the case of the WWW network.

preprint2012arXiv

The spread of gossip in American schools

We study different mechanisms of gossip propagation on several network topologies and introduce a new network property, the "spread factor", describing the fraction of neighbors that get to know the gossip. We postulate that for scale-free networks the spreading time grows logarithmically with the degree of the victim and prove this statement for the case of the Apollonian network. Applying our concepts to real data from an American school survey, we confirm the logarithmic law and disclose that there exists an ideal number of acquaintances minimizing the fraction attained by the gossip. The similarity between the school survey and scale-free networks remains even for cases when gossip propagation only occurs with some probability $q<1$. The spreading times follow an exponential distribution that can also be calculated analytically for the Apollonian network. When gossip also spreads through strangers the situation changes substantially: the spreading time becomes a constant and there exists no ideal degree of connectivity anymore.

preprint2011arXiv

Density profiles of loose and collapsed cohesive granular structures generated by ballistic deposition

Loose granular structures stabilized against gravity by an effective cohesive force are investigated on a microscopic basis using contact dynamics. We study the influence of the granular Bond number on the density profiles and the generation process of packings, generated by ballistic deposition under gravity. The internal compaction occurs discontinuously in small avalanches and we study their size distribution. We also develop a model explaining the final density profiles based on insight about the collapse of a packing under changes of the Bond number.

preprint2011arXiv

From invasion percolation to flow in rock fracture networks

The main purpose of this work is to simulate two-phase flow in the form of immiscible displacement through anisotropic, three-dimensional (3D) discrete fracture networks (DFN). The considered DFNs are artificially generated, based on a general distribution function or are conditioned on measured data from deep geological investigations. We introduce several modifications to the invasion percolation (MIP) to incorporate fracture inclinations, intersection lines, as well as the hydraulic path length inside the fractures. Additionally a trapping algorithm is implemented that forbids any advance of the invading fluid into a region, where the defending fluid is completely encircled by the invader and has no escape route. We study invasion, saturation, and flow through artificial fracture networks, with varying anisotropy and size and finally compare our findings to well studied, conditioned fracture networks.

preprint2011arXiv

Inverse targeting -- an effective immunization strategy

We propose a new method to immunize populations or computer networks against epidemics which is more efficient than any method considered before. The novelty of our method resides in the way of determining the immunization targets. First we identify those individuals or computers that contribute the least to the disease spreading measured through their contribution to the size of the largest connected cluster in the social or a computer network. The immunization process follows the list of identified individuals or computers in inverse order, immunizing first those which are most relevant for the epidemic spreading. We have applied our immunization strategy to several model networks and two real networks, the Internet and the collaboration network of high energy physicists. We find that our new immunization strategy is in the case of model networks up to 14%, and for real networks up to 33% more efficient than immunizing dynamically the most connected nodes in a network. Our strategy is also numerically efficient and can therefore be applied to large systems.

preprint2011arXiv

Mitigation of Malicious Attacks on Networks

Terrorist attacks on transportation networks have traumatized modern societies. With a single blast, it has become possible to paralyze airline traffic, electric power supply, ground transportation or Internet communication. How and at which cost can one restructure the network such that it will become more robust against a malicious attack? We introduce a unique measure for robustness and use it to devise a method to mitigate economically and efficiently this risk. We demonstrate its efficiency on the European electricity system and on the Internet as well as on complex networks models. We show that with small changes in the network structure (low cost) the robustness of diverse networks can be improved dramatically while their functionality remains unchanged. Our results are useful not only for improving significantly with low cost the robustness of existing infrastructures but also for designing economically robust network systems.

preprint2011arXiv

Model for a dune field with exposed water table

Aeolian transport in coastal areas can be significantly affected by the presence of an exposed water table. In some dune fields, such as in Lençóis Maranhenses, northeastern Brazil, the water table oscillates in response to seasonal changes of rainfall and rates of evapotranspiration, rising above the ground during the wet season and sinking below in the dry period. A quantitative understanding of dune mobility in an environment with varying groundwater level is essential for coastal management as well as for the study of long-term evolution of many dune fields. Here we apply a model for aeolian dunes to study the genesis of coastal dune fields in presence of an oscillating water table. We find that the morphology of the field depends on the time cycle, $T_{\mathrm{w}}$, of the water table and the maximum height, $H_{\mathrm{w}}$, of its oscillation. Our calculations show that long chains of barchanoids alternating with interdune ponds such as found at Lençóis Maranhenses arise when $T_{\mathrm{w}}$ is of the order of the dune turnover time, whereas $H_{\mathrm{w}}$ dictates the growth rate of dune height with distance downwind. We reproduce quantitatively the morphology and size of dunes at Lençóis Maranhenses, as well as the total relative area between dunes.

preprint2011arXiv

Modeling the topology of protein interaction networks

A major issue in biology is the understanding of the interactions between proteins. These interactions can be described by a network, where the proteins are modeled by nodes and the interactions by edges. The origin of these protein networks is not well understood yet. Here we present a two-step model, which generates clusters with the same topological properties as networks for protein-protein interactions, namely, the same degree distribution, cluster size distribution, clustering coefficient and shortest path length. The biological and model networks are not scale free but exhibit small world features. The model allows the fitting of different biological systems by tuning a single parameter.

preprint2011arXiv

Scale free networks by preferential depletion

We show that not only preferential attachment but also preferential depletion leads to scale-free networks. The resulting degree distribution exponents is typically less than two (5/3) as opposed to the case of the growth models studied before where the exponents are larger. Our approach applies in particular to biological networks where in fact we find interesting agreement with experimental measurements. We investigate the most important properties characterizing these networks, as the cluster size distribution, the average shortest path and the clustering coefficient.

preprint2011arXiv

Simulation of particle mixing in turbulent channel flow due to intrinsic fluid velocity fluctuation

We combine a DEM simulation with a stochastic process to model the movement of spherical particles in a turbulent channel flow. With this model we investigate the mixing properties of two species of particles flowing through the channel. We find a linear increase of the mixing zone with the length of the pipe. Flows at different Reynolds number are studied. Below a critical Reynolds number at the Taylor microscale of around $R_{c} \approx 300$ the mixing rate is strongly dependent on the Reynolds number. Above $R_{c}$ the mixing rate stays nearly constant.

preprint2011arXiv

Size distribution and structure of Barchan dune fields

Barchans are isolated mobile dunes often organized in large dune fields. Dune fields seem to present a characteristic dune size and spacing, which suggests a cooperative behavior based on dune interaction. In Duran et al. (2009), we propose that the redistribution of sand by collisions between dunes is a key element for the stability and size selection of barchan dune fields. This approach was based on a mean-field model ignoring the spatial distribution of dune fields. Here, we present a simplified dune field model that includes the spatial evolution of individual dunes as well as their interaction through sand exchange and binary collisions. As a result, the dune field evolves towards a steady state that depends on the boundary conditions. Comparing our results with measurements of Moroccan dune fields, we find that the simulated fields have the same dune size distribution as in real fields but fail to reproduce their homogeneity along the wind direction.

preprint2011arXiv

Suppressing Epidemics with a Limited Amount of Immunization Units

The way diseases spread through schools, epidemics through countries, and viruses through the Internet is crucial in determining their risk. Although each of these threats has its own characteristics, its underlying network determines the spreading. To restrain the spreading, a widely used approach is the fragmentation of these networks through immunization, so that epidemics cannot spread. Here we develop an immunization approach based on optimizing the susceptible size, which outperforms the best known strategy based on immunizing the highest-betweenness links or nodes. We find that the network's vulnerability can be significantly reduced, demonstrating this on three different real networks: the global flight network, a school friendship network, and the internet. In all cases, we find that not only is the average infection probability significantly suppressed, but also for the most relevant case of a small and limited number of immunization units the infection probability can be reduced by up to 55%.

preprint2011arXiv

Transverse instability of dunes

The simplest type of dune is the transverse one, which propagates with invariant profile orthogonally to a fixed wind direction. Here we show numerically and with a linear stability analysis that transverse dunes are unstable with respect to along-axis perturbations in their profile and decay on the bedrock into barchan dunes. Any forcing modulation amplifies exponentially with growth rate determined by the dune turnover time. We estimate the distance covered by a transverse dune before fully decaying into barchans and identify the patterns produced by different types of perturbation.

preprint2011arXiv

Tricritical point in explosive percolation

The suitable interpolation between classical percolation and a special variant of explosive percolation enables the explicit realization of a tricritical percolation point. With high-precision simulations of the order parameter and the second moment of the cluster size distribution a fully consistent tricritical scaling scenario emerges yielding the tricritical crossover exponent $1/ϕ_t=1.8\pm0.1$.

preprint2011arXiv

Watersheds and Explosive percolation

The recent work by Achlioptas, D'Souza, and Spencer opened up the possibility of obtaining a discontinuous (explosive) percolation transition by changing the stochastic rule of bond occupation. Despite the active research on this subject, several questions still remain open about the leading mechanism and the properties of the system. We review the largest cluster and the Gaussian models recently introduced. We show that, to obtain a discontinuous transition it is solely necessary to control the size of the largest cluster, suppressing the growth of a cluster differing significantly, in size, from the average one. As expected for a discontinuous transition, a Gaussian cluster-size distribution and compact clusters are obtained. The surface of the clusters is fractal, with the same fractal dimension of the watershed line.

preprint2010arXiv

A micromechanical model of collapsing quicksand

The discrete element method constitutes a general class of modeling techniques to simulate the microscopic behavior (i.e. at the particle scale) of granular/soil materials. We present a contact dynamics method, accounting for the cohesive nature of fine powders and soils. A modification of the model adjusted to capture the essential physical processes underlying the dynamics of generation and collapse of loose systems is able to simulate "quicksand" behavior of a collapsing soil material, in particular of a specific type, which we call "living quicksand". We investigate the penetration behavior of an object for varying density of the material. We also investigate the dynamics of the penetration process, by measuring the relation between the driving force and the resulting velocity of the intruder, leading to a "power law" behavior with exponent 1/2, i.e. a quadratic velocity dependence of the drag force on the intruder.

preprint2010arXiv

Inter-arrival times of message propagation on directed networks

One of the challenges in fighting cybercrime is to understand the dynamics of message propagation on botnets, networks of infected computers used to send viruses, unsolicited commercial emails (SPAM) or denial of service attacks. We map this problem to the propagation of multiple random walkers on directed networks and we evaluate the inter-arrival time distribution between successive walkers arriving at a target. We show that the temporal organization of this process, which models information propagation on unstructured peer to peer networks, has the same features as SPAM arriving to a single user. We study the behavior of the message inter-arrival time distribution on three different network topologies using two different rules for sending messages. In all networks the propagation is not a pure Poisson process. It shows universal features on Poissonian networks and a more complex behavior on scale free networks. Results open the possibility to indirectly learn about the process of sending messages on networks with unknown topologies, by studying inter-arrival times at any node of the network.

preprint2010arXiv

Learning as a phenomenon occurring in a critical state

Recent physiological measurements have provided clear evidence about scale-free avalanche brain activity and EEG spectra, feeding the classical enigma of how such a chaotic system can ever learn or respond in a controlled and reproducible way. Models for learning, like neural networks or perceptrons, have traditionally avoided strong fluctuations. Conversely, we propose that brain activity having features typical of systems at a critical point, represents a crucial ingredient for learning. We present here a study which provides novel insights toward the understanding of the problem. Our model is able to reproduce quantitatively the experimentally observed critical state of the brain and, at the same time, learns and remembers logical rules including the exclusive OR (XOR), which has posed difficulties to several previous attempts. We implement the model on a network with topological properties close to the functionality network in real brains. Learning occurs via plastic adaptation of synaptic strengths and exhibits universal features. We find that the learning performance and the average time required to learn are controlled by the strength of plastic adaptation, in a way independent of the specific task assigned to the system. Even complex rules can be learned provided that the plastic adaptation is sufficiently slow.

preprint2010arXiv

Protein Interaction Networks are Fragile against Random Attacks and Robust against Malicious Attacks

The capacity to resist attacks from the environment is crucial to the survival of all organisms. We quantitatively analyze the susceptibility of protein interaction networks of numerous organisms to random and malicious attacks. We find for all organisms studied that random rewiring improves protein network robustness, so that actual networks are more fragile than rewired surrogates. This unexpected fragility contrasts with the behavior of networks such as the Internet, whose robustness decreases with random rewiring. We trace this surprising effect to the modular structure of protein networks.

preprint2010arXiv

q-state Potts model on the Apollonian network

The q-state Potts model is studied on the Apollonian network with Monte Carlo simulations and the Transfer Matrix method. The spontaneous magnetization, correlation length, entropy, and specific heat are analyzed as a function of temperature for different number of states, $q$. Different scaling functions in temperature and $q$ are proposed. A quantitative agreement is found between results from both methods. No critical behavior is observed in the thermodynamic limit for any number of states.

preprint2010arXiv

Self-contact and instabilities in the anisotropic growth of elastic membranes

We investigate the morphology of thin discs and rings growing in the circumferential direction. Recent analytical results suggest that this growth produces symmetric excess cones (e-cones). We study the stability of such solutions considering self-contact and bending stress. We show that, contrary to what was assumed in previous analytical solutions, beyond a critical growth factor, no symmetric \textit{e}-cone solution is energetically minimal any more. Instead, we obtain skewed e-cone solutions having lower energy, characterized by a skewness angle and repetitive spiral winding with increasing growth. These results are generalized to discs with varying thickness and rings with holes of different radii.

preprint2010arXiv

Superdiffusion of massive particles induced by multi-scale velocity fields

We study drag-induced diffusion of massive particles in scale-free velocity fields, where superdiffusive behavior emerges due to the scale-free size distribution of the vortices of the underlying velocity field. The results show qualitative resemblance to what is observed in fluid systems, namely the diffusive exponent for the mean square separation of pairs of particles and the preferential concentration of the particles, both as a function of the response time.

preprint2010arXiv

Tactical Voting in Plurality Elections

How often will elections end in landslides? What is the probability for a head-to-head race? Analyzing ballot results from several large countries rather anomalous and yet unexplained distributions have been observed. We identify tactical voting as the driving ingredient for the anomalies and introduce a model to study its effect on plurality elections, characterized by the relative strength of the feedback from polls and the pairwise interaction between individuals in the society. With this model it becomes possible to explain the polarization of votes between two candidates, understand the small margin of victories frequently observed for different elections, and analyze the polls' impact in American, Canadian, and Brazilian ballots. Moreover, the model reproduces, quantitatively, the distribution of votes obtained in the Brazilian mayor elections with two, three, and four candidates.

preprint2010arXiv

Transport on exploding percolation clusters

We propose a simple generalization of the explosive percolation process [Achlioptas et al., Science 323, 1453 (2009)], and investigate its structural and transport properties. In this model, at each step, a set of q unoccupied bonds is randomly chosen. Each of these bonds is then associated with a weight given by the product of the cluster sizes that they would potentially connect, and only that bond among the q-set which has the smallest weight becomes occupied. Our results indicate that, at criticality, all finite-size scaling exponents for the spanning cluster, the conducting backbone, the cutting bonds, and the global conductance of the system, change continuously and significantly with q. Surprisingly, we also observe that systems with intermediate values of q display the worst conductive performance. This is explained by the strong inhibition of loops in the spanning cluster, resulting in a substantially smaller associated conducting backbone.

preprint2009arXiv

Biased Percolation on Scale-free Networks

Biased (degree-dependent) percolation was recently shown to provide new strategies for turning robust networks fragile and vice versa. Here we present more detailed results for biased edge percolation on scale-free networks. We assume a network in which the probability for an edge between nodes $i$ and $j$ to be retained is proportional to $(k_ik_j)^{-α}$ with $k_i$ and $k_j$ the degrees of the nodes. We discuss two methods of network reconstruction, sequential and simultaneous, and investigate their properties by analytical and numerical means. The system is examined away from the percolation transition, where the size of the giant cluster is obtained, and close to the transition, where nonuniversal critical exponents are extracted using the generating functions method. The theory is found to agree quite well with simulations. By introducing an extension of the Fortuin-Kasteleyn construction, we find that biased percolation is well described by the $q\to 1$ limit of the $q$-state Potts model with inhomogeneous couplings.

preprint2009arXiv

Swimming in Granular Media

We study a simple model of periodic contraction and extension of large intruders in a granular bed to understand the mechanism for swimming in an otherwise solid media. Using an event-driven simulation, we find optimal conditions that idealized swimmers must use to critically fluidize a sand bed so that it is rigid enough to support a load when needed, but fluid enough to permit motion with minimal resistance. Swimmers - or other intruders - that agitate the bed too rapidly produce large voids that prevent traction from being achieved, while swimmers that move too slowly cannot travel before the bed re-solidifies around them i.e., the swimmers locally probe the fundamental time-scale in a granular packing.

preprint2008arXiv

Avalanches in anisotropic sheared granular media

We study the influence of particle shape anisotropy on the occurrence of avalanches in sheared granular media. We use molecular dynamic simulations to calculate the relative movement of two tectonic plates. % with transform boundaries. Our model considers irregular polygonal particles constituting the material within the shear zone. We find that the magnitude of the avalanches is approximately independent on particle shape and in good agreement with the Gutenberg-Richter law, but the aftershock sequences are strongly influenced by the particle anisotropy yielding variations on the exponent characterizing the empirical Omori's law. Our findings enable one to identify the presence of anisotropic particles at the macro-mechanical level only by observing the avalanche sequences of real faults. In addition, we calculate the probability of occurrence of an avalanche for given values of stiffness or frictional strength and observe also a significant influence of the particle anisotropy.

preprint2008arXiv

The dune size distribution and scaling relations of barchan dune fields

Barchan dunes emerge as a collective phenomena involving the generation of thousands of them in so called barchan dune fields. By measuring the size and position of dunes in Moroccan barchan dune fields, we find that these dunes tend to distribute uniformly in space and follow an unique size distribution function. We introduce an analyticalmean-field approach to show that this empirical size distribution emerges from the interplay of dune collisions and sand flux balance, the two simplest mechanisms for size selection. The analytical model also predicts a scaling relation between the fundamental macroscopic properties characterizing a dune field, namely the inter-dune spacing and the first and second moments of the dune size distribution.

preprint2007arXiv

Analytical approach to directed sandpile models on the Apollonian network

We investigate a set of directed sandpile models on the Apollonian network, which are inspired on the work by Dhar and Ramaswamy (PRL \textbf{63}, 1659 (1989)) for Euclidian lattices. They are characterized by a single parameter $q$, that restricts the number of neighbors receiving grains from a toppling node. Due to the geometry of the network, two and three point correlation functions are amenable to exact treatment, leading to analytical results for the avalanche distributions in the limit of an infinite system, for $q=1,2$. The exact recurrence expressions for the correlation functions are numerically iterated to obtain results for finite size systems, when larger values of $q$ are considered. Finally, a detailed description of the local flux properties is provided by a multifractal scaling analysis.

preprint2007arXiv

Modeling slow deformation of polygonal particles using DEM

We introduce two improvements in the numerical scheme to simulate collision and slow shearing of irregular particles. First, we propose an alternative approach based on simple relations to compute the frictional contact forces. The approach improves efficiency and accuracy of the Discrete Element Method (DEM) when modeling the dynamics of the granular packing. We determine the proper upper limit for the integration step in the standard numerical scheme using a wide range of material parameters. To this end, we study the kinetic energy decay in a stress controlled test between two particles. Second, we show that the usual way of defining the contact plane between two polygonal particles is, in general, not unique which leads to discontinuities in the direction of the contact plane while particles move. To solve this drawback, we introduce an accurate definition for the contact plane based on the shape of the overlap area between touching particles, which evolves continuously in time.

preprint2007arXiv

New approaches to model and study social networks

We describe and develop three recent novelties in network research which are particularly useful for studying social systems. The first one concerns the discovery of some basic dynamical laws that enable the emergence of the fundamental features observed in social networks, namely the nontrivial clustering properties, the existence of positive degree correlations and the subdivision into communities. To reproduce all these features we describe a simple model of mobile colliding agents, whose collisions define the connections between the agents which are the nodes in the underlying network, and develop some analytical considerations. The second point addresses the particular feature of clustering and its relationship with global network measures, namely with the distribution of the size of cycles in the network. Since in social bipartite networks it is not possible to measure the clustering from standard procedures, we propose an alternative clustering coefficient that can be used to extract an improved normalized cycle distribution in any network. Finally, the third point addresses dynamical processes occurring on networks, namely when studying the propagation of information in them. In particular, we focus on the particular features of gossip propagation which impose some restrictions in the propagation rules. To this end we introduce a quantity, the spread factor, which measures the average maximal fraction of nearest neighbors which get in contact with the gossip, and find the striking result that there is an optimal non-trivial number of friends for which the spread factor is minimized, decreasing the danger of being gossiped.

preprint2007arXiv

Numerical Improvement of the Discrete Element Method applied to Shear of Granular Media

We present a detailed analysis of the bounds on the integration step in Discrete Element Method (DEM) for simulating collisions and shearing of granular assemblies. We show that, in the numerical scheme, the upper limit for the integration step, usually taken from the average time $t_c$ of one contact, is in fact not sufficiently small to guarantee numerical convergence of the system during relaxation. In particular, we study in detail how the kinetic energy decays during the relaxation stage and compute the correct upper limits for the integration step, which are significantly smaller than the ones commonly used. In addition, we introduce an alternative approach, based on simple relations to compute the frictional forces, that converges even for integration steps above the upper limit.

preprint2007arXiv

Obtaining the size distribution of fault gouges with polydisperse bearings

We generalize the recent study of random space-filling bearings to a more realistic situation, where the spacing offset varies randomly during the space-filling procedure, and show that it reproduces well the size-distributions observed in recent studies of real fault gouges. In particular, we show that the fractal dimensions of random polydisperse bearings sweep predominantly the low range of values in the spectrum of fractal dimensions observed along real faults, which strengthen the evidence that polydisperse bearings may explain the occurrence of seismic gaps in nature. In addition, the influence of different distributions for the offset is studied and we find that the uniform distribution is the best choice for reproducing the size-distribution of fault gouges.

preprint2007arXiv

Spreading gossip in social networks

We study a simple model of information propagation in social networks, where two quantities are introduced: the spread factor, which measures the average maximal fraction of neighbors of a given node that interchange information among each other, and the spreading time needed for the information to reach such fraction of nodes. When the information refers to a particular node at which both quantities are measured, the model can be taken as a model for gossip propagation. In this context, we apply the model to real empirical networks of social acquaintances and compare the underlying spreading dynamics with different types of scale-free and small-world networks. We find that the number of friendship connections strongly influences the probability of being gossiped. Finally, we discuss how the spread factor is able to be applied to other situations.

preprint2006arXiv

Computer Simulation of Particle Suspensions

Particle suspensions are ubiquitous in our daily life, but are not well understood due to their complexity. During the last twenty years, various simulation methods have been developed in order to model these systems. Due to varying properties of the solved particles and the solvents, one has to choose the simulation method properly in order to use the available compute resources most effectively with resolving the system as well as needed. Various techniques for the simulation of particle suspensions have been implemented at the Institute for Computational Physics allowing us to study the properties of clay-like systems, where Brownian motion is important, more macroscopic particles like glass spheres or fibers solved in liquids, or even the pneumatic transport of powders in pipes. In this paper we will present the various methods we applied and developed and discuss their individual advantages.

preprint2006arXiv

Networks based on collisions among mobile agents

We investigate in detail a recent model of colliding mobile agents [Phys. Rev. Lett.~96, 088702], used as an alternative approach to construct evolving networks of interactions formed by the collisions governed by suitable dynamical rules. The system of mobile agents evolves towards a quasi-stationary state which is, apart small fluctuations, well characterized by the density of the system and the residence time of the agents. The residence time defines a collision rate and by varying the collision rate, the system percolates at a critical value, with the emergence of a giant cluster whose critical exponents are the ones of two-dimensional percolation. Further, the degree and clustering coefficient distributions and the average path length show that the network associated with such a system presents non-trivial features which, depending on the collision rule, enables one not only to recover the main properties of standard networks, such as exponential, random and scale-free networks, but also to obtain other topological structures. Namely, we show a specific example where the obtained structure has topological features which characterize accurately the structure and evolution of social networks in different contexts, ranging from networks of acquaintances to networks of sexual contacts.

preprint2006arXiv

Plurality Voting: the statistical laws of democracy in Brazil

We explore the statistical laws behind the plurality voting system by investigating the election results for mayor held in Brazil in 2004. Our analysis indicate that the vote partition among mayor candidates of the same city tends to be "polarized" between two candidates, a phenomenon that can be closely described by means of a simple fragmentation model. Complex concepts like "government continuity" and "useful vote" can be identified and even statistically quantified through our approach.

preprint2005arXiv

Coherence in complex networks of oscillators

We study fully synchronized (coherent) states in complex networks of chaotic oscillators, reviewing the analytical approach of determining the stability conditions for synchronizability and comparing them with numerical criteria. As an example, we present detailed results for networks of chaotic logistic maps having three different scale-free topologies: random scale-free topology, deterministic pseudo-fractal scale-free network and Apollonian network. For random scale-free topology we find that the lower boundary of the synchronizability region scales approximately as $k^{-μ}$, where $k$ is the outgoing connectivity and $μ$ depends on the local nonlinearity. For deterministic scale-free networks coherence is observed only when the coupling is heterogeneous, namely when it is proportional to some power of the neighbor connectivity. In all cases, stability conditions are determined from the eigenvalue spectrum of the Laplacian matrix and agree well with numerical results based on histograms of coherent states in parameter space. Additionally, we show that almost everywhere in the synchronizability region the basin of attraction of coherent states fills the entire phase space, and that the transition to coherence is of first-order.

preprint2005arXiv

Cycles and clustering in bipartite networks

We investigate the clustering ability in bipartite networks where cycles of size three are absent and therefore the standard definition of clustering coefficient cannot be used. Instead, we use another coefficient given by the fraction of cycles with size four, showing that both coefficients yield the same clustering properties. The new coefficient is computed for two networks of sexual contacts, one monopartite and another bipartite. In both cases the clustering ability is similar. Furthermore, combining both clustering coefficients we deduce an expression for estimating cycles of larger size, which improves previous estimations and is suitable for either monopartite and multipartite networks.

preprint2004arXiv

Fragmentation of a Circular Disc by Impact on a Frictionless Plate

The break-up of a two-dimensional circular disc by normal and oblique impact on a hard frictionless plate is investigated by molecular dynamics simulations. The disc is composed of numerous unbreakable randomly shaped convex polygons connected together by simple elastic beams that break when bent or stretched beyond a certain limit. It is found that for both normal and oblique impacts the crack patterns are the same and depend solely on the normal component of the impact velocity. Analysing the pattern of breakage, amount of damage, fragment masses and velocities, we show the existence of a critical velocity which separates two regimes of the impact process: below the critical point only a damage cone is formed at the impact site (damage), cleaving of the particle occurs at the critical point, while above the critical velocity the disc breaks into several pieces (fragmentation). In the limit of very high impact velocities the disc suffers complete disintegration (shattering) into many small fragments. In agreement with experimental results, fragment masses are found to follow the Gates-Gaudin-Schuhmann distribution (power law) with an exponent independent of the velocity and angle of impact. The velocity distribution of fragments exhibit an interesting anomalous scaling behavior when changing the impact velocity and the size of the disc.

preprint2003arXiv

Modelling transverse dunes

Transverse dunes appear in regions of mainly unidirectional wind and high sand availability. A dune model is extended to two dimensional calculation of the shear stress. It is applied to simulate dynamics and morphology of transverse dunes which seem to reach translational invariance and do not stop growing. Hence, simulations of two dimensional dune fields have been performed. Characteristic laws were found for the time evolution of transverse dunes. Bagnold's law of the dune velocity is modified and reproduced. The interaction between transverse dunes led to interesting results which conclude that small dunes can pass through bigger ones.

preprint2002arXiv

Scaling laws of creep rupture of fiber bundles

We study the creep rupture of fiber composites in the framework of fiber bundle models. Two novel fiber bundle models are introduced based on different microscopic mechanisms responsible for the macroscopic creep behavior. Analytical and numerical calculations show that above a critical load the deformation of the creeping system monotonically increases in time resulting in global failure at a finite time $t_f$, while below the critical load the system suffers only partial failure and the deformation tends to a constant value giving rise to an infinite lifetime. It is found that approaching the critical load from below and above the creeping system is characterized by universal power laws when the fibers have long range interaction. The lifetime of the composite above the critical point has a universal dependence on the system size.

preprint2002arXiv

Size segregation and convection

The size segregation of granular materials in a vibrating container is investigated using Molecular Dynamics. We find that the rising of larger particles is accompanied by the existence of convection cells even in the case of the lowest possible frequencies. The convection can, however, also be triggered by the larger particle itself. The possibility of rising through this mechanism strongly depends on the depth of the larger particle.

preprint1998arXiv

Broad Histogram Method for Continuous Systems: the XY-Model

We propose a way of implementing the Broad Histogram Monte Carlo method to systems with continuous degrees of freedom, and we apply these ideas to investigate the three-dimensional XY-model with periodic boundary conditions. We have found an excellent agreement between our method and traditional Metropolis results for the energy, the magnetization, the specific heat and the magnetic susceptibility on a very large temperature range. For the calculation of these quantities in the temperature range 0.7<T<4.7 our method took less CPU time than the Metropolis simulations for 16 temperature points in that temperature range. Furthermore, it calculates the whole temperature range 1.2<T<4.7 using only 2.2 times more computer effort than the Histogram Monte Carlo method for the range 2.1<T<2.2. Our way of treatment is general, it can also be applied to other systems with continuous degrees of freedom.

preprint1995arXiv

A study of fragmentation processes using a discrete element method

We present a model of solids made from polygonal cells connected via beams. We calculate the macroscopic elastic moduli from the beam and cell parameters. This modellisation is particularly suited for the simulation of fragmentation processes. We study the effects of an explosion inside a circular disk and the impact of a projectile and obtain the fragment size distribution. We find that if breaking only happens under tensile forces a layer on the free wall opposed to impact is first ejected. In that case the distribution follows a power-law with an exponent that in most cases is around two.