Source author record

H. J. Herrmann

H. J. Herrmann appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

106works

36topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2020arXiv

Dynamics of extended Schelling models

We explore extensions of Schelling's model of social dynamics, in which two types of agents live on a checkerboard lattice and move in order to optimize their own satisfaction, which depends on how many agents among their neighbors are of their same type. For each number $n$ of same-type nearest neighbors we independently assign a binary satisfaction variable $s_{k}$ which is equal to one only if the agent is satisfied with that condition, and is equal to zero otherwise. This defines 32 different satisfaction rules, which we investigate in detail, focusing on pattern formation and measuring segregation with the help of an "energy" function which is related to the number of neighboring agents of different types and plays no role in the dynamics. We consider the checkerboard lattice to be fully occupied and the dynamics consists of switching the locations of randomly selected unsatisfied agents of opposite types. We show that, starting from a random distribution of agents, only a small number of rules lead to (nearly) fully segregated patterns in the long run, with many rules leading to chaotic steady-state behavior. Nevertheless, other interesting patterns may also be dynamically generated, such as "anti-segregate d" patterns as well as patterns resembling sponges.

preprint2020arXiv

Frustrated bearings

In a bearing state, touching spheres (disks in two dimensions) roll on each other without slip. Here we frustrate a system of touching spheres by imposing two different bearing states on opposite sides and search for the configurations of lowest energy dissipation. If the dissipation between contacts of spheres is viscous (with random damping constants), the angular momentum continuously changes from one bearing state to the other. For Coulomb friction (with random friction coefficients) in two dimensions, a sharp line separates the two bearing states and we show that this line corresponds to the minimum cut. Astonishingly however, in three dimensions, intermediate bearing domains, that are not synchronized with either side, are energetically more favorable than the minimum-cut surface. Instead of a sharp cut, the steady state displays a fragmented structure. This novel type of state of minimum dissipation is characterized by a spanning network of slipless contacts that reaches every sphere. Such a situation becomes possible because in three dimensions bearing states have four degrees of freedom.

preprint2020arXiv

Geometry-induced non-equilibrium phase transition in sandpiles

We study the sandpile model on three-dimensional spanning Ising clusters with the temperature $T$ treated as the control parameter. By analyzing the three dimensional avalanches and their two-dimensional projections (which show scale-invariant behavior for all temperatures), we uncover two universality classes with different exponents (an ordinary BTW class, and SOC$_{T=\infty}$), along with a tricritical point (at $T_c$, the critical temperature of the host) between them. The transition between these two criticalities is induced by the transition in the support. The SOC$_{T=\infty}$ universality class is characterized by the exponent of the avalanche size distribution $τ^{T=\infty}=1.27\pm 0.03$, consistent with the exponent of the size distribution of the Barkhausen avalanches in amorphous ferromagnets (Phys. Rev. L 84, 4705 (2000)). The tricritical point is characterized by its own critical exponents. In addition to the avalanche exponents, some other quantities like the average height, the spanning avalanche probability (SAP) and the average coordination number of the Ising clusters change significantly the behavior at this point, and also exhibit power-law behavior in terms of $ε\equiv \frac{T-T_c}{T_c}$, defining further critical exponents. Importantly the finite size analysis for the activity (number of topplings) per site shows the scaling behavior with exponents $β=0.19\pm 0.02$ and $ν=0.75\pm 0.05$. A similar behavior is also seen for the SAP and the average avalanche height. The fractal dimension of the external perimeter of the two-dimensional projections of avalanches is shown to be robust against $T$ with the numerical value $D_f=1.25\pm 0.01$.

preprint2020arXiv

The elastic backbone phase transition in the Ising model

The two-dimensional (zero magnetic field) Ising model is known to undergo a second order para-ferromagnetic phase transition, which is accompanied by a correlated percolation transition for the Fortuin-Kasteleyn (FK) clusters. In this paper we uncover that there exists also a second temperature $T_{\text{eb}}<T_c$ at which the elastic backbone of FK clusters undergoes a second order phase transition to a dense phase. The corresponding universality class, which is characterized by determining various percolation exponents, is shown to be completely different from directed percolation, proposing a new anisotropic universality class with $β=0.54\pm 0.02$, $ν_{||}=1.86\pm 0.01$, $ν_{\perp}=1.21\pm 0.04$ and $d_f=1.53\pm 0.03$. All tested hyper-scaling relations are shown to be valid.

preprint2016arXiv

Conformal Invariance of Graphene Sheets

Suspended graphene sheets exhibit correlated random deformations that can be studied under the framework of rough surfaces with a Hurst (roughness) exponent $0.72 \pm 0.01$. Here, we show that, independent of the temperature, the iso-height lines at the percolation threshold have a well-defined fractal dimension and are conformally invariant, sharing the same statistical properties as Schramm-Loewner evolution (SLE$_κ$) curves with $κ=2.24\pm0.07$. Interestingly, iso-height lines of other rough surfaces are not necessarily conformally invariant even if they have the same Hurst exponent, e.g. random Gaussian surfaces. We have found that the distribution of the modulus of the Fourier coefficients plays an important role on this property. Our results not only introduce a new universality class and place the study of suspended graphene membranes within the theory of critical phenomena, but also provide hints on the long-standing question about the origin of conformal invariance in iso-height lines of rough surfaces.

preprint2016arXiv

Connectivity disruption sparks explosive epidemic spreading

We investigate the spread of an infection or other malfunction of cascading nature when a system component can recover only if it remains reachable from a functioning central component. We consider the susceptible-infected-susceptible model, typical of mathematical epidemiology, on a network. Infection spreads from infected to healthy nodes, with the addition that infected nodes can only recover when they remain connected to a predefined central node, through a path that contains only healthy nodes. In this system, clusters of infected nodes will absorb their noninfected interior because no path exists between the central node and encapsulated nodes. This gives rise to the simultaneous infection of multiple nodes. Interestingly, the system converges to only one of two stationary states: either the whole population is healthy or it becomes completely infected. This simultaneous cluster infection can give rise to discontinuous jumps of different sizes in the number of failed nodes. Larger jumps emerge at lower infection rates. The network topology has an important effect on the nature of the transition: we observed hysteresis for networks with dominating local interactions. Our model shows how local spread can abruptly turn uncontrollable when it disrupts connectivity at a larger spatial scale.

preprint2016arXiv

Critical fragmentation properties of random drilling: How many random holes need to be drilled to collapse a wooden cube?

A solid wooden cube fragments into pieces as we sequentially drill holes through it randomly. This seemingly straightforward observation encompasses deep and nontrivial geometrical and probabilistic behavior that is discussed here. Combining numerical simulations and rigorous results, we find off-critical scale-free behavior and a continuous transition at a critical density of holes that significantly differs from classical percolation.

preprint2016arXiv

Crumpling Damaged Graphene

Through molecular mechanics we find that non-covalent interactions modify the fractality of crumpled damaged graphene. Pristine graphene membranes are damaged by adding random vacancies and carbon-hydrogen bonds. Crumpled membranes exhibit a fractal dimension of $ 2.71 \pm 0.02$ when all interactions between carbon atoms are considered, and $2.30 \pm 0.05$ when non-covalent interactions are suppressed. The transition between these two values, obtained by switching on/off the non-covalent interactions of equilibrium configurations, is shown to be reversible and independent on thermalisation. In order to explain this transition, we propose a theoretical model that is compatible with our numerical findings. Finally, we also compare damaged graphene membranes with other crumpled structures, as for instance, polymerised membranes and paper sheets, that share similar scaling properties.

preprint2016arXiv

Cutting self-similar space-filling sphere packings

Any space-filling packing of spheres can be cut by a plane to obtain a space-filling packing of disks. Here, we deal with space-filling packings generated using inversive geometry leading to exactly self-similar fractal packings. First, we prove that cutting along a random hyperplane leads in general to a packing with a fractal dimension of the one of the uncut packing minus one. Second, we find special cuts which can be constructed themselves by inversive geometry. Such special cuts have specific fractal dimensions, which we demonstrate by cutting a three- and a four-dimensional packing. The increase in the number of found special cuts with respect to a cutoff parameter suggests the existence of infinitely many topologies with distinct fractal dimensions.

preprint2016arXiv

Emergence of core-peripheries in networks

A number of important transport networks, such as the airline and trade networks of the world, exhibit a characteristic core-periphery structure, wherein a few nodes are highly interconnected and the rest of the network frays into a tree. Mechanisms underlying the emergence of core-peripheries, however, remain elusive. Here, we demonstrate that a simple pruning process based on removal of underutilized links and redistribution of loads can lead to the emergence of core-peripheries. Links are assumed beneficial if they either carry a sufficiently large load or are essential for global connectivity. This incentivized redistribution process is controlled by a single parameter which balances connectivity and profit. The obtained networks exhibit a highly resilient and connected core with a frayed periphery. The balanced network shows a higher resilience than the World Airline Network or the World Trade Network, revealing a pathway towards robust structural features through pruning.

preprint2016arXiv

Energy dissipation in flows through curved spaces

Fluid dynamics in intrinsically curved geometries is encountered in many physical systems in nature, ranging from microscopic bio-membranes all the way up to general relativity at cosmological scales. Despite the diversity of applications, all of these systems share a common feature: the free motion of particles is affected by inertial forces originating from the curvature of the embedding space. Here we reveal a fundamental process underlying fluid dynamics in curved space: the free motion of fluids, in the complete absence of solid walls or obstacles, exhibits loss of energy due exclusively to the intrinsic curvature of space. We find that local sources of curvature generate viscous stresses as a result of the inertial forces. The curvature-induced viscous forces are shown to cause hitherto unnoticed and yet appreciable energy dissipation, which might play a significant role for a variety of physical systems involving fluid dynamics in curved spaces.

preprint2016arXiv

Field Driven Charging Dynamics of a Fluidized Granular Bed

A simplified model has previously described the inductive charging of colliding identical grains in the presence of an external electric field. Here we extend that model by including heterogeneous surface charge distributions, grain rotations and electrostatic interactions between grains. We find from this more realistic model that strong heterogeneities in charging can occur in agitated granular beds, and we predict that shielding due to these heterogeneities can dramatically alter the charging rate in such beds.

preprint2016arXiv

Gender gap in the ERASMUS mobility program

Studying abroad has become very popular among students. The ERASMUS mobility program is one of the largest international student exchange programs in the world, which has supported already more than three million participants since 1987. We analyzed the mobility pattern within this program in 2011-12 and found a gender gap across countries and subject areas. Namely, for almost all participating countries, female students are over-represented in the ERASMUS program when compared to the entire population of tertiary students. The same tendency is observed across different subject areas. We also found a gender asymmetry in the geographical distribution of hosting institutions, with a bias of male students in Scandinavian countries. However, a detailed analysis reveals that this latter asymmetry is rather driven by subject and consistent with the distribution of gender ratios among subject areas.

preprint2016arXiv

Hybrid Percolation Transition in Cluster Merging Processes: Continuously Varying Exponents

Consider growing a network, in which every new connection is made between two disconnected nodes. At least one node is chosen randomly from a subset consisting of $g$ fraction of the entire population in the smallest clusters. Here we show that this simple strategy for improving connection exhibits a phase transition barely studied before, namely a hybrid percolation transition exhibiting the properties of both first-order and second-order phase transitions. The cluster size distribution of finite clusters at a transition point exhibits power-law behavior with a continuously varying exponent $τ$ in the range $2 < τ(g) \le 2.5$. This pattern reveals a necessary condition for a hybrid transition in cluster aggregation processes, which is comparable to the power-law behavior of the avalanche size distribution arising in models with link-deleting processes in interdependent networks.

preprint2016arXiv

Model for the growth of the World Airline Network

We propose a probabilistic growth model for transport networks which employs a balance between popularity of nodes and the physical distance between nodes. By comparing the degree of each node in the model network and the WAN, we observe that the difference between the two is minimized for $α\approx 2$. Interestingly, this is the value obtained for the node-node correlation function in the WAN. This suggests that our model explains quite well the growth of airline networks.

preprint2016arXiv

Poiseuille flow in curved spaces

We investigate Poiseuille channel flow through intrinsically curved media, equipped with localized metric perturbations. To this end, we study the flux of a fluid driven through the curved channel in dependence of the spatial deformation, characterized by the parameters of the metric perturbations (amplitude, range and density). We find that the flux depends only on a specific combination of parameters, which we identify as the average metric perturbation, and derive a universal flux law for the Poiseuille flow. For the purpose of this study, we have improved and validated our recently developed lattice Boltzmann model in curved space by considerably reducing discrete lattice effects.

preprint2016arXiv

Prediction and control of slip-free rotation states in sphere assemblies

We study fixed assemblies of touching spheres that can individually rotate. From any initial state, sliding friction drives an assembly toward a slip-free rotation state. For bipartite assemblies, which have only even loops, this state has at least four degrees of freedom. For exactly four degrees of freedom, we analytically predict the final state, which we prove to be independent of the strength of sliding friction, from an arbitrary initial one. With a tabletop experiment, we show how to impose any slip-free rotation state by only controlling two spheres, regardless of the total number.

preprint2016arXiv

Self-similar space-filling sphere packings in three and four dimensions

Inversive geometry can be used to generate exactly self-similar space-filling sphere packings. We present a construction method in two dimensions and generalize it to search for packings in higher dimensions. We newly discover 29 three-dimensional and 13 four-dimensional topologies of which 10 and 5, respectively, are bearings. To distinguish and characterize the packing topologies, we numerically estimate their fractal dimensions and we analyze their contact networks.

preprint2016arXiv

Synaptic plasticity and neuronal refractory time cause scaling behaviour of neuronal avalanches

Neuronal avalanches measured in vitro and in vivo in different cortical networks consistently exhibit power law behaviour for the size and duration distributions with exponents typical for a mean field self-organized branching process. These exponents are also recovered in neuronal network simulations implementing various neuronal dynamics on different network topologies. They can therefore be considered a very robust feature of spontaneous neuronal activity. Interestingly, this scaling behaviour is also observed on regular lattices in finite dimensions, which raises the question about the origin of the mean field behaviour observed experimentally. In this study we provide an answer to this open question by investigating the effect of activity dependent plasticity in combination with the neuronal refractory time in a neuronal network. Results show that the refractory time hinders backward avalanches forcing a directed propagation. Hebbian plastic adaptation plays the role of sculpting these directed avalanche patterns into the topology of the network slowly changing it into a branched structure where loops are marginal.

preprint2015arXiv

Analysis of the velocity field of granular hopper flow

We report the analysis of radial characteristics of the flow of granular material through a conical hopper. The discharge is simulated for various orifice sizes and hopper opening angles. Velocity profiles are measured along two radial lines from the hopper cone vertex: along the main axis of the cone and along its wall. An approximate power law dependence on the distance from the orifice is observed for both profiles, although differences between them can be noted. In order to quantify these differences, we propose a Local Mass Flow index that is a promising tool in the direction of a more reliable classification of the flow regimes in hoppers.

preprint2015arXiv

Bursts in discontinuous Aeolian saltation

Close to the onset of Aeolian particle transport through saltation we find in wind tunnel experiments a regime of discontinuous flux characterized by bursts of activity. Scaling laws are observed in the time delay between each burst and in the measurements of the wind fluctuations at the fluid threshold Shields number $θ_c$. The time delay between each burst decreases on average with the increase of the Shields number until sand flux becomes continuous. A numerical model for saltation including the wind-entrainment from the turbulent fluctuations can reproduce these observations and gives insight about their origin. We present here also for the first time measurements showing that with feeding it becomes possible to sustain discontinuous flux even below the fluid threshold.

preprint2015arXiv

Complex networks from space-filling bearings

Two dimensional space-filling bearings are dense packings of disks that can rotate without slip. We consider the entire first family of bearings for loops of size four and propose a hierarchical construction of their contact network. We provide analytic expressions for the clustering coefficient and degree distribution, revealing bipartite scale-free behavior with tunable degree exponent depending on the bearing parameters. We also analyze their average shortest path and percolation properties.

preprint2015arXiv

Cooling Effect of the Richtmyer-Meshkov Instability

We provide numerical evidence that the Richtmyer-Meshkov (RM) instability contributes to the cooling of a relativistic fluid. Due to the presence of jet particles traveling throughout the medium, shock waves are generated in the form of Mach cones. The interaction of multiple shock waves can trigger the RM instability, and we have found that this process leads to a down-cooling of the relativistic fluid. To confirm the cooling effect of the instability, shock tube Richtmyer-Meshkov instability simulations are performed. Additionally, in order to provide an experimental observable of the RM instability resulting from the Mach cone interaction, we measure the two particle correlation function and highlight the effects of the interaction. The simulations have been performed with an improved version of the relativistic lattice Boltzmann model, including general equations of state and external forces.

preprint2015arXiv

Coupled DEM-LBM method for the free-surface simulation of heterogeneous suspensions

The complexity of the interactions between the constituent granular and liquid phases of a suspension requires an adequate treatment of the constituents themselves. A promising way for numerical simulations of such systems is given by hybrid computational frameworks. This is naturally done, when the Lagrangian description of particle dynamics of the granular phase finds a correspondence in the fluid description. In this work we employ extensions of the Lattice-Boltzmann Method for non-Newtonian rheology, free surfaces, and moving boundaries. The models allows for a full coupling of the phases, but in a simplified way. An experimental validation is given by an example of gravity driven flow of a particle suspension.

preprint2015arXiv

Critical Cooperation Range to Improve Spatial Network Robustness

A robust worldwide air-transportation network (WAN) is one that minimizes the number of stranded passengers under a sequence of airport closures. Building on top of this realistic example, here we address how spatial network robustness can profit from cooperation between local actors. We swap a series of links within a certain distance, a cooperation range, while following typical constraints of spatially embedded networks. We find that the network robustness is only improved above a critical cooperation range. Such improvement can be described in the framework of a continuum transition, where the critical exponents depend on the spatial correlation of connected nodes. For the WAN we show that, except for Australia, all continental networks fall into the same universality class. Practical implications of this result are also discussed.

preprint2015arXiv

En route to a multi-model scheme for clinker comminution with chemical grinding aids

We present a multi-model simulation approach, targeted at understanding the behavior of comminution and the effect of grinding aids (GAs) in industrial cement mills. On the atomistic scale we use Molecular Dynamics (MD) simulations with validated force field models to quantify elastic and structural properties, cleavage energies as well as the organic interactions with mineral surfaces. Simulations based on the Discrete Element Method (DEM) are used to integrate the information gained from MD simulations into the clinker particle behavior at larger scales. Computed impact energy distributions from DEM mill simulations can serve as a link between large-scale industrial and laboratory sized mills. They also provide the required input for particle impact fragmentation models. Such a multi-scale, multi-model methodology paves the way for a structured approach to the design of chemical additives aimed at improving mill performance.

preprint2015arXiv

Flow through randomly curved manifolds

We have found that the relation between the flow through campylotic (generically curved) media, consisting of randomly located curvature perturbations, and the average Ricci scalar of the system exhibits two distinct functional expressions (hysteresis), depending on whether the typical spatial extent of the curvature perturbation lies above or below the critical value maximizing the overall Ricci curvature. Furthermore, the flow through such systems as a function of the number of curvature perturbations presents a sublinear behavior for large concentrations due to the interference between curvature perturbations that, consequently, produces a less curved space. For the purpose of this study, we have developed and validated a lattice kinetic model capable of describing fluid flow in arbitrarily curved manifolds, which allows to deal with highly complex spaces in a very compact and efficient way.

preprint2015arXiv

Granular front formation in free-surface flow of concentrated suspensions

Granular fronts are a common yet unexplained phenomenon emerging during the gravity driven free-surface flow of concentrated suspensions. They are usually believed to be the result of fluid convection in combination with particle size segregation. However, suspensions composed of uniformly sized particles also develop a granular front. Within a large rotating drum, a stationary recirculating avalanche is generated. The flowing material is a mixture of a visco-plastic fluid obtained from a kaolin-water dispersion, with spherical ceramic particles denser than the fluid. The goal is to mimic the composition of many common granular-fluid materials, like fresh concrete or debris flow. In these materials, granular and fluid phases have the natural tendency to segregate due to particle settling. However, through the shearing caused by the rotation of the drum, a reorganization of the phases is induced, leading to the formation of a granular front. By tuning the material properties and the drum velocity, it is possible to control this phenomenon. The setting is reproduced in a numerical environment, where the fluid is solved by a Lattice-Boltzmann Method, and the particles are explicitly represented using the Discrete Element Method. The simulations confirm the findings of the experiments, and provide insight into the internal mechanisms. Comparing the time-scale of particle settling with the one of particle recirculation, a non-dimensional number is defined, and is found to be effective in predicting the formation of a granular front.

preprint2015arXiv

High-Order Kinetic Relaxation Schemes as High-Accuracy Poisson Solvers

We present a new approach to find accurate solutions to the Poisson equation, as obtained from the steady-state limit of a diffusion equation with strong source terms. For this purpose, we start from Boltzmann's kinetic theory and investigate the influence of higher order terms on the resulting macroscopic equations. By performing an appropriate expansion of the equilibrium distribution, we provide a method to remove the unnecessary terms up to a desired order and show that it is possible to find, with high level of accuracy, the steady-state solution of the diffusion equation for sizeable Knudsen numbers. In order to test our kinetic approach, we discretise the Boltzmann equation and solve the Poisson equation, spending up to six order of magnitude less computational time for a given precision than standard lattice Boltzmann methods.

preprint2015arXiv

Lattice Boltzmann Model for Electronic Structure Simulations

Recently, a new connection between density functional theory and kinetic theory has been proposed. In particular, it was shown that the Kohn-Sham (KS) equations can be reformulated as a macroscopic limit of the steady-state solution of a suitable single-particle kinetic equation. By using a discrete version of this new formalism, the exchange and correlation energies of simple atoms and the geometrical configuration of the methane molecule were calculated accurately. Here, we discuss the main ideas behind the lattice kinetic approach to electronic structure computations, offer some considerations for prospective extensions, and also show additional numerical results, namely the geometrical configuration of the water molecule.

preprint2015arXiv

Lattice Boltzmann model for resistive relativistic magnetohydrodynamics

In this paper, we develop a lattice Boltzmann model for relativistic magnetohydrodynamics (MHD). Even though the model is derived for resistive MHD, it is shown that it is numerically robust even in the high conductivity (ideal MHD) limit. In order to validate the numerical method, test simulations are carried out for both ideal and resistive limits, namely the propagation of Alfvén waves in the ideal MHD and the evolution of current sheets in the resistive regime, where very good agreement is observed comparing to the analytical results. Additionally, two-dimensional magnetic reconnection driven by Kelvin-Helmholtz instability is studied and the effects of different parameters on the reconnection rate are investigated. It is shown that the density ratio has negligible effect on the magnetic reconnection rate, while an increase in shear velocity decreases the reconnection rate. Additionally, it is found that the reconnection rate is proportional to $σ^{-\frac{1}{2}}$, $σ$ being the conductivity, which is in agreement with the scaling law of the Sweet-Parker model. Finally, the numerical model is used to study the magnetic reconnection in a stellar flare. Three-dimensional simulation suggests that the reconnection between the background and flux rope magnetic lines in a stellar flare can take place as a result of a shear velocity in the photosphere.

preprint2015arXiv

Schramm-Loewner Evolution and isoheight lines of correlated landscapes

Real landscapes are usually characterized by long-range height-height correlations, which are quantified by the Hurst exponent $H$. We analyze the statistical properties of the isoheight lines for correlated landscapes of $H\in [-1,1]$. We show numerically that, for $H\leq 0$ the statistics of these lines is compatible with $SLE$ and that established analytic results are recovered for $H=-1$ and $H=0$. This result suggests that for negative $H$, in spite of the long-range nature of correlations, the statistics of isolines is fully encoded in a Brownian motion with a single parameter in the continuum limit. By contrast, for positive $H$ we find that the one-dimensional time series encoding the isoheight lines is not Markovian and therefore not consistent with $SLE$.

preprint2015arXiv

Stochastic Loewner Evolution Relates Anomalous Diffusion and Anisotropic Percolation

We disclose the origin of anisotropic percolation perimeters in terms of the Stochastic Loewner Evolution (SLE) process. Precisely, our results from extensive numerical simulations indicate that the perimeters of multi-layered and directed percolation clusters at criticality are the scaling limits of the Loewner evolution of an anomalous Brownian motion, being subdiffusive and superdiffusive, respectively. The connection between anomalous diffusion and fractal anisotropy is further tested by using long-range power-law correlated time series (fractional Brownian motion) as driving functions in the evolution process. The fact that the resulting traces are distinctively anisotropic corroborates our hypothesis. Under the conceptual framework of SLE, our study therefore reveals new perspectives for mathematical and physical interpretations of non-Markovian processes in terms of anisotropic paths at criticality and vice-versa.

preprint2014arXiv

Dean Instability in Double-Curved Channels

We study the Dean instability in curved channels using the lattice Boltzmann model for generalized metrics. For this purpose, we first improve and validate the method by measuring the critical Dean number at the transition from laminar to vortex flow for a streamwise curved rectangular channel, obtaining very good agreement with the literature values. Taking advantage of the easy implementation of arbitrary metrics within our model, we study the fluid flow through a double-curved channel, using ellipsoidal coordinates, and study the transition to vortex flow in dependence of the two perpendicular curvature radii of the channel. We observe not only transitions to 2-cell vortex flow, but also to 4-cell and even 6-cell vortex flow, and we find that the critical Dean number at the transition to 2-cell vortex flow exhibits a minimum when the two curvature radii are approximately equal.

preprint2014arXiv

Flood avalanches in a semiarid basin with a dense reservoir network

This study investigates flood avalanches in a dense reservoir network in the semiarid north-eastern Brazil. The population living in this area strongly depends on the availability of the water from this network. Water is stored during intense wet-season rainfall events and evaporates from the reservoir surface during the dry season. These seasonal changes are the driving forces behind the water dynamics in the network. The reservoir network and its connectivity properties during flood avalanches are investigated with a model called ResNetM, which simulates each reservoir explicitly. It runs on the basis of daily calculated water balances for each reservoir. A spilling reservoir contributes with water to the reservoir downstream, which can trigger avalanches affecting, in some cases, large fractions of the network. The main focus is on the study of the relation between the total amount of water stored and the largest observable cluster of connected reservoirs that overspill in the same day. It is shown that the thousands of small and middle-sized reservoirs are eminent for the retention of water upstream the large ones. Therefore, they prevent large clusters at a low level of water. Concerning connectivity measures, the actual reservoir network, which evolved without an integrated plan, performed better (i.e., generated smaller avalanches for similar amount of stored water) than numerous stochastically generated artificial reservoir networks on the same river network.

preprint2014arXiv

Mandala Networks: ultra-robust, ultra-small-world and highly sparse graphs

The increasing demands in security and reliability of infrastructures call for the optimal design of their embedded complex networks topologies. The following question then arises: what is the optimal layout to fulfill best all the demands? Here we present a general solution for this problem with scale-free networks, like the Internet and airline networks. Precisely, we disclose a way to systematically construct networks which are 100$\%$ robust against random failures as well as to malicious attacks. Furthermore, as the sizes of these networks increase, their shortest paths become asymptotically invariant and densities of links go to zero, making them ultra-small worlds and highly sparse, respectively. The first property is ideal for communication and navigation purposes, while the second is interesting economically.

preprint2014arXiv

Modelling the influence of photospheric turbulence on solar flare statistics

Solar flares stem from the reconnection of twisted magnetic field lines in the solar photosphere. The energy and waiting time distributions of these events follow complex patterns that have been carefully considered in the past and that bear some resemblance with earthquakes and stockmarkets. Here we explore in detail the tangling motion of interacting flux tubes anchored in the plasma and the energy ejections resulting when they recombine. The mechanism for energy accumulation and release in the flow is reminiscent of self-organized criticality. From this model we suggest the origin for two important and widely studied properties of solar flare statistics, including the time-energy correlations. We first propose that the scale-free energy distribution of solar flares is largely due to the twist exerted by the vorticity of the turbulent photosphere. Second, the long-range temporal and time-energy correlations appear to arise from the tube-tube interactions. The agreement with satellite measurements is encouraging.

preprint2014arXiv

Retention capacity of correlated surfaces

We extend the water retention model [C. L. Knecht et al., Phys. Rev. Lett. 108, 045703 (2012)] to correlated random surfaces. We find that the retention capacity of discrete random landscapes is strongly affected by spatial correlations among the heights. This phenomenon is related to the emergence of power-law scaling in the lake volume distribution. We also solve the uncorrelated case exactly for a small lattice and present bounds on the retention of uncorrelated landscapes.

preprint2014arXiv

Shortest path and Schramm-Loewner Evolution

We numerically show that the statistical properties of the shortest path on critical percolation clusters are consistent with the ones predicted for Schramm-Loewner evolution (SLE) curves for $κ=1.04\pm0.02$. The shortest path results from a global optimization process. To identify it, one needs to explore an entire area. Establishing a relation with SLE permits to generate curves statistically equivalent to the shortest path from a Brownian motion. We numerically analyze the winding angle, the left passage probability, and the driving function of the shortest path and compare them to the distributions predicted for SLE curves with the same fractal dimension. The consistency with SLE opens the possibility of using a solid theoretical framework to describe the shortest path and it raises relevant questions regarding conformal invariance and domain Markov properties, which we also discuss.

preprint2014arXiv

Usage leading to an abrupt collapse of connectivity

Network infrastructures are essential for the distribution of resources such as electricity and water. Typical strategies to assess their resilience focus on the impact of a sequence of random or targeted failures of network nodes or links. Here we consider a more realistic scenario, where elements fail based on their usage. We propose a dynamic model of transport based on the Bak-Tang-Wiesenfeld sandpile model where links fail after they have transported more than an amount $μ$ (threshold) of the resource and we investigate it on the square lattice. As we deal with a new model, we provide insight on its fundamental behavior and dependence on parameters. We observe that for low values of the threshold due to a positive feedback of link failure, an avalanche develops that leads to an abrupt collapse of the lattice. By contrast, for high thresholds the lattice breaks down in an uncorrelated fashion. We determine the critical threshold $μ^*$ separating these two regimes and show how it depends on the toppling threshold of the nodes and the mass increment added stepwise to the system. We find that the time of major disconnection is well described with a linear dependence on $μ$. Furthermore, we propose a lower bound for $μ^*$ by measuring the strength of the dynamics leading to abrupt collapses.

preprint2014arXiv

Watersheds in disordered media

What is the best way to divide a rugged landscape? Since ancient times, watersheds separating adjacent water systems that flow, for example, toward different seas, have been used to delimit boundaries. Interestingly, serious and even tense border disputes between countries have relied on the subtle geometrical properties of these tortuous lines. For instance, slight and even anthropogenic modifications of landscapes can produce large changes in a watershed, and the effects can be highly nonlocal. Although the watershed concept arises naturally in geomorphology, where it plays a fundamental role in water management, landslide, and flood prevention, it also has important applications in seemingly unrelated fields such as image processing and medicine. Despite the far-reaching consequences of the scaling properties on watershed-related hydrological and political issues, it was only recently that a more profound and revealing connection has been disclosed between the concept of watershed and statistical physics of disordered systems. This review initially surveys the origin and definition of a watershed line in a geomorphological framework to subsequently introduce its basic geometrical and physical properties. Results on statistical properties of watersheds obtained from artificial model landscapes generated with long-range correlations are presented and shown to be in good qualitative and quantitative agreement with real landscapes.

preprint2013arXiv

Avoiding a Spanning Cluster in Percolation Models

When dynamics in a system proceeds under suppressive external bias, the system can undergo an abrupt phase transition, as it occurs for example in the epidemic spreading. Recently, an explosive percolation (EP) model was introduced in line with such phenomena. The order of the EP transition has not been clarified in a unified framework covering low dimensional systems and the mean-field limit. We introduce a stochastic model, in which a rule for dynamics is designed to avoid the formation of a spanning cluster through competitive selection in Euclidean space. We show by heuristic arguments that, in the thermodynamic limit and depending on a control parameter, the EP transition can be either continuous or discontinuous if $d < d_c$ and is always continuous if $d \geq d_c$, where $d$ is the spatial dimension and $d_c$ the upper critical dimension.

preprint2013arXiv

Breathing synchronization in interconnected networks

Global synchronization in a complex network of oscillators emerges from the interplay between its topology and the dynamics of the pairwise interactions among its numerous components. When oscillators are spatially separated, however, a time delay appears in the interaction which might obstruct synchronization. Here we study the synchronization properties of interconnected networks of oscillators with a time delay between networks and analyze the dynamics as a function of the couplings and communication lag. We discover a new breathing synchronization regime, where two groups appear in each network synchronized at different frequencies. Each group has a counterpart in the opposite network, one group is in phase and the other in anti-phase with their counterpart. For strong couplings, instead, networks are internally synchronized but a phase shift between them might occur. The implications of our findings on several socio-technical and biological systems are discussed.

preprint2013arXiv

Exact evaluation of the cutting path length in a percolation model on a hierarchical network

This work presents an approach to evaluate the exact value of the fractal dimension of the cutting path $d^{CP}_f$ on hierarchical structures with finite order of ramification. This represents the first renormalization group treatment of the universality class of watersheds. By making use of the self-similar property, we show that $d^{CP}_f$ depends only on the average cutting path (CP) of the first generation of the structure. For the simplest Wheastone hierarchical lattice (WHL), we present a mathematical proof. For a larger WHL structure, the exact value of $d^{CP}_f$ is derived based on an computer algorithm that identifies the length of all possible CP's of the first generation.

preprint2013arXiv

Formal analogy between the Dirac equation in its Majorana form and the discrete-velocity version of the Boltzmann kinetic equation

We point out a formal analogy between the Dirac equation in Majorana form and the discrete-velocity version of the Boltzmann kinetic equation. By a systematic analysis based on the theory of operator splitting, this analogy is shown to turn into a concrete and efficient computational method, providing a unified treatment of relativistic and non-relativistic quantum mechanics. This might have potentially far-reaching implications for both classical and quantum computing, because it shows that, by splitting time along the three spatial directions, quantum information (Dirac-Majorana wavefunction) propagates in space-time as a classical statistical process (Boltzmann distribution).

preprint2013arXiv

Gaussian Quadrature and Lattice Discretization of the Fermi-Dirac Distribution for Graphene

We construct a lattice kinetic scheme to study electronic flow in graphene. For this purpose, we first derive a basis of orthogonal polynomials, using as weight function the ultrarelativistic Fermi-Dirac distribution at rest. Later, we use these polynomials to expand the respective distribution in a moving frame, for both cases, undoped and doped graphene. In order to discretize the Boltzmann equation and make feasible the numerical implementation, we reduce the number of discrete points in momentum space to 18 by applying a Gaussian quadrature, finding that the family of representative wave (2+1)-vectors, that satisfies the quadrature, reconstructs a honeycomb lattice. The procedure and discrete model are validated by solving the Riemann problem, finding excellent agreement with other numerical models. In addition, we have extended the Riemann problem to the case of different dopings, finding that by increasing the chemical potential, the electronic fluid behaves as if it increases its effective viscosity.

preprint2013arXiv

Geometrical and transport properties of sequential adsorption clusters

We investigate transport properties of percolating clusters generated by irreversible cooperative sequential adsorption (CSA) on square lattices with Arrhenius rates given by ki= q^(ni), where ni is the number of occupied neighbors of the site i, and q a controlling parameter. Our results show a dependence of the prefactors on q and a strong finite size effect for small values of this parameter, both impacting the size of the backbone and the global conductance of the system. These results might be pertinent to practical applications in processes involving adsorption of particles.

preprint2013arXiv

Hydrodynamic Model for Conductivity in Graphene

Based on the recently developed picture of an electronic ideal relativistic fluid at the Dirac point, we present an analytical model for the conductivity in graphene that is able to describe the linear dependence on the carrier density and the existence of a minimum conductivity. The model treats impurities as submerged rigid obstacles, forming a disordered medium through which graphene electrons flow, in close analogy with classical fluid dynamics. To describe the minimum conductivity, we take into account the additional carrier density induced by the impurities in the sample. The model, which predicts the conductivity as a function of the impurity fraction of the sample, is supported by extensive simulations for different values of ${\cal E}$, the dimensionless strength of the electric field, and provides excellent agreement with experimental data.

preprint2013arXiv

Kinetic Formulation of the Kohn-Sham Equations for ab initio Electronic Structure Calculations

We introduce a new approach to density functional theory based on kinetic theory, showing that the Kohn-Sham equations can be derived as a macroscopic limit of a suitable Boltzmann kinetic equation in the limit of small mean free path versus the typical scale of density gradients (Chapman-Enskog expansion). To derive the approach, we first write the Schrödinger equation as a special case of a Boltzmann equation for a gas of quasi-particles, with the potential playing the role of an external source that generates and destroys particles, so as to drive the system towards the ground state. The ions are treated as classical particles, using the Born-Oppenheimer dynamics, or by imposing concurrent evolution with the electronic orbitals. In order to provide quantitative support to our approach, we implement a discrete (lattice) model and compute, the exchange and correlation energies of simple atoms, and the geometrical configuration of the methane molecule. Excellent agreement with values in the literature is found.

preprint2013arXiv

Lattice Boltzmann model for ultra-relativistic flows

We develop a relativistic lattice Boltzmann model capable of describing relativistic fluid dynamics at ultra-high velocities, with Lorentz factors up to $γ\sim 10$. To this purpose, we first build a new lattice kinetic scheme by expanding the Maxwell-Jüttner distribution function in an orthogonal basis of polynomials and applying an appropriate quadrature, providing the discrete versions of the relativistic Boltzmann equation and the equilibrium distribution. To achieve ultra-high velocities, we include a flux limiter scheme, and introduce the bulk viscosity by a suitable extension of the discrete relativistic Boltzmann equation. The model is validated by performing simulations of shock waves in viscous quark-gluon plasmas and comparing with existing models, finding very good agreement. To the best of our knowledge, we for the first time successfully simulate viscous shock waves in the highly relativistic regime. Moreover, we show that our model can also be used for near-inviscid flows even at very high velocities. Finally, as an astrophysical application, we simulate a relativistic shock wave, generated by, say, a supernova explosion, colliding with a massive interstellar cloud, e.g. molecular gas.

preprint2013arXiv

Optimal synchronizability of bearings

Bearings are mechanical dissipative systems that, when perturbed, relax toward a synchronized (bearing) state. Here we find that bearings can be perceived as physical realizations of complex networks of oscillators with asymmetrically weighted couplings. Accordingly, these networks can exhibit optimal synchronization properties through fine tuning of the local interaction strength as a function of node degree [Motter, Zhou, and Kurths, Phys. Rev. E 71, 016116 (2005)]. We show that, in analogy, the synchronizability of bearings can be maximized by counterbalancing the number of contacts and the inertia of their constituting rotor disks through the mass-radius relation, $m\sim r^α$, with an optimal exponent $α=α_{\times}$ which converges to unity for a large number of rotors. Under this condition, and regardless of the presence of a long-tailed distribution of disk radii composing the mechanical system, the average participation per disk is maximized and the energy dissipation rate is homogeneously distributed among elementary rotors.

preprint2013arXiv

Percolation with long-range correlated disorder

Long-range power-law correlated percolation is investigated using Monte Carlo simulations. We obtain several static and dynamic critical exponents as function of the Hurst exponent $H$ which characterizes the degree of spatial correlation among the occupation of sites. In particular, we study the fractal dimension of the largest cluster and the scaling behavior of the second moment of the cluster size distribution, as well as the complete and accessible perimeters of the largest cluster. Concerning the inner structure and transport properties of the largest cluster, we analyze its shortest path, backbone, red sites, and conductivity. Finally, bridge site growth is also considered. We propose expressions for the functional dependence of the critical exponents on $H$.

preprint2013arXiv

Relativistic effects on the Richtmyer-Meshkov instability

Theoretical and numerical analysis of the relativistic effects on the Richtmyer-Meshkov (RM) instability reveals new and potentially very useful effects. We find that, in contrast with the non- relativistic case, the growth rate of the RM instability depends strongly on the equation of state of the fluid, opening up the possibility to infer equations of state from experimental observations of the RM instability. As opposed to the non-relativistic case, we also discover that, above a critical value of the fluid velocity, the growth rate of the instability counter-intuitively decreases due to the Lorentz's factor, and vanishes in the ultrarelativistic limit, as the speed of the particles approaches the speed of light. Both effects might prove very useful for leading-edge applications, such as the study of the equation of state of quark-gluon matter, and the design of fast ignition inertial confinement fusion (ICF) schemes. We perform a linear stability analysis to characterize the instability, for an arbitrary equation of state, and implement numerical simulations to study the instability in the non-linear regime, using the equation of state of an ideal gas. Furthermore, based on the numerical results, we propose a general expression that characterizes the long term evolution of the instability.

preprint2013arXiv

Relativistic Lattice Boltzmann Model with Improved Dissipation

We develop a relativistic lattice Boltzmann (LB) model, providing a more accurate description of dissipative phenomena in relativistic hydrodynamics than previously available with existing LB schemes. The procedure applies to the ultra-relativistic regime, in which the kinetic energy (temperature) far exceeds the rest mass energy, although the extension to massive particles and/or low temperatures is conceptually straightforward. In order to improve the description of dissipative effects, the Maxwell-Juettner distribution is expanded in a basis of orthonormal polynomials, so as to correctly recover the third order moment of the distribution function. In addition, a time dilatation is also applied, in order to preserve the compatibility of the scheme with a cartesian cubic lattice. To the purpose of comparing the present LB model with previous ones, the time transformation is also applied to a lattice model which recovers terms up to second order, namely up to energy-momentum tensor. The approach is validated through quantitative comparison between the second and third order schemes with BAMPS (the solution of the full relativistic Boltzmann equation), for moderately high viscosity and velocities, and also with previous LB models in the literature. Excellent agreement with BAMPS and more accurate results than previous relativistic lattice Boltzmann models are reported.

preprint2013arXiv

Smart Rewiring for Network Robustness

While new forms of attacks are developed every day to compromise essential infrastructures, service providers are also expected to develop strategies to mitigate the risk of extreme failures. In this context, tools of Network Science have been used to evaluate network robustness and propose resilient topologies against attacks. We present here a new rewiring method to modify the network topology improving its robustness, based on the evolution of the network largest component during a sequence of targeted attacks. In comparison to previous strategies, our method lowers by several orders of magnitude the computational effort necessary to improve robustness. Our rewiring also drives the formation of layers of nodes with similar degree while keeping a highly modular structure. This "modular onion-like structure" is a particular class of the onion-like structure previously described in the literature. We apply our rewiring strategy to an unweighted representation of the World Air Transportation network and show that an improvement of 30% in its overall robustness can be achieved through smart swaps of around 9% of its links.

preprint2013arXiv

Stacked triangular lattice: Percolation properties

The stacked triangular lattice has the shape of a triangular prism. In spite of being considered frequently in solid state physics and materials science, its percolation properties have received few attention. We investigate several non-universal percolation properties on this lattice using Monte Carlo simulation. We show that the percolation threshold is $p_c^\text{bond}=0.186\;02\pm0.000\;02$ for bonds and $p_c^\text{site}=0.262\;40\pm0.000\;05$ for sites. The number of clusters at the threshold per site is $n_c^\text{bond}=0.284\;58\pm0.000\;05$ and $n_c^\text{site}=0.039\;98\pm0.000\;05$. The stacked triangular lattice is a convenient choice to study the RGB model [Sci. Rep. {\bf 2}, 751 (2012)]. We present results on this model and its scaling behavior at the percolation threshold.

preprint2013arXiv

Ultrarelativistic Transport Coefficients in Two Dimensions

We compute the shear and bulk viscosities, as well as the thermal conductivity of an ultrarelativistic fluid obeying the relativistic Boltzmann equation in 2+1 space-time dimensions. The relativistic Boltzmann equation is taken in the single relaxation time approximation, based on two approaches, the first, due to Marle and using the Eckart decomposition, and the second, proposed by Anderson and Witting and using the Landau-Lifshitz decomposition. In both cases, the local equilibrium is given by a Maxwell-Juettner distribution. It is shown that, apart from slightly different numerical prefactors, the two models lead to a different dependence of the transport coefficients on the fluid temperature, quadratic and linear, for the case of Marle and Anderson-Witting, respectively. However, by modifying the Marle model according to the prescriptions given in Ref.[1], it is found that the temperature dependence becomes the same as for the Anderson-Witting model.

preprint2012arXiv

Bohman-Frieze-Wormald model on the lattice, yielding a discontinuous percolation transition

The BFW model introduced by Bohman, Frieze, and Wormald [Random Struct. Algorithms, 25, 432 (2004)] and recently investigated in the framework of discontinuous percolation by Chen and D'Souza [Phys. Rev. Lett., 106, 115701 (2011)], is studied on the square and simple-cubic lattices. In two and three dimensions, we find numerical evidence for a strongly discontinuous transition. In two dimensions, the clusters at the threshold are compact with a fractal surface of fractal dimension $d_f=1.49\pm0.02$. On the simple-cubic lattice, distinct jumps in the size of the largest cluster are observed. We proceed to analyze the tree-like version of the model, where only merging bonds are sampled, for dimension two to seven. The transition is again discontinuous in any considered dimension. Finally, the dependence of the cluster-size distribution at the threshold on the spatial dimension is also investigated.

preprint2012arXiv

Corrections to Scaling for Watersheds, Optimal Path Cracks, and Bridge Lines

We study the corrections to scaling for the mass of the watershed, the bridge line, and the optimal path crack in two and three dimensions. We disclose that these models have numerically equivalent fractal dimensions and leading correction-to-scaling exponents. We conjecture all three models to possess the same fractal dimension, namely, $d_f=1.2168\pm0.0005$ in 2D and $d_f=2.487\pm0.003$ in 3D, and the same exponent of the leading correction, $Ω=0.9\pm0.1$ and $Ω=1.0\pm0.1$, respectively. The close relations between watersheds, optimal path cracks in the strong disorder limit, and bridge lines are further supported by either heuristic or exact arguments.

preprint2012arXiv

Fracturing highly disordered materials

We investigate the role of disorder on the fracturing process of heterogeneous materials by means of a two-dimensional fuse network model. Our results in the extreme disorder limit reveal that the backbone of the fracture at collapse, namely the subset of the largest fracture that effectively halts the global current, has a fractal dimension of $1.22 \pm 0.01$. This exponent value is compatible with the universality class of several other physical models, including optimal paths under strong disorder, disordered polymers, watersheds and optimal path cracks on uncorrelated substrates, hulls of explosive percolation clusters, and strands of invasion percolation fronts. Moreover, we find that the fractal dimension of the largest fracture under extreme disorder, $d_f=1.86 \pm 0.01$, is outside the statistical error bar of standard percolation. This discrepancy is due to the appearance of trapped regions or cavities of all sizes that remain intact till the entire collapse of the fuse network, but are always accessible in the case of standard percolation. Finally, we quantify the role of disorder on the structure of the largest cluster, as well as on the backbone of the fracture, in terms of a distinctive transition from weak to strong disorder characterized by a new crossover exponent.

preprint2012arXiv

Fracturing ranked surfaces

Discretized landscapes can be mapped onto ranked surfaces, where every element (site or bond) has a unique rank associated with its corresponding relative height. By sequentially allocating these elements according to their ranks and systematically preventing the occupation of bridges, namely elements that, if occupied, would provide global connectivity, we disclose that bridges hide a new tricritical point at an occupation fraction $p=p_{c}$, where $p_{c}$ is the percolation threshold of random percolation. For any value of $p$ in the interval $p_{c}< p \leq 1$, our results show that the set of bridges has a fractal dimension $d_{BB} \approx 1.22$ in two dimensions. In the limit $p \rightarrow 1$, a self-similar fracture is revealed as a singly connected line that divides the system in two domains. We then unveil how several seemingly unrelated physical models tumble into the same universality class and also present results for higher dimensions.

preprint2012arXiv

How to share underground reservoirs

Many resources, such as oil, gas, or water, are extracted from porous soils and their exploration is often shared among different companies or nations. We show that the effective shares can be obtained by invading the porous medium simultaneously with various fluids. Partitioning a volume in two parts requires one division surface while the simultaneous boundary between three parts consists of lines. We identify and characterize these lines, showing that they form a fractal set consisting of a single thread spanning the medium and a surrounding cloud of loops. While the spanning thread has fractal dimension ${1.55\pm0.03}$, the set of all lines has dimension ${1.69\pm0.02}$. The size distribution of the loops follows a power law and the evolution of the set of lines exhibits a tricritical point described by a crossover with a negative dimension at criticality.

preprint2012arXiv

How to suppress undesired synchronization

It is delightful to observe the emergence of synchronization in the blinking of fireflies to attract partners and preys. Other charming examples of synchronization can also be found in a wide range of phenomena such as, e.g., neurons firing, lasers cascades, chemical reactions, and opinion formation. However, in many situations the formation of a coherent state is not pleasant and should be mitigated. For example, the onset of synchronization can be the root of epileptic seizures, traffic congestion in communication networks, and the collapse of constructions. Here we propose the use of contrarians to suppress undesired synchronization. We perform a comparative study of different strategies, either requiring local or total knowledge of the system, and show that the most efficient one solely requires local information. Our results also reveal that, even when the distribution of neighboring interactions is narrow, significant improvement in mitigation is observed when contrarians sit at the highly connected elements. The same qualitative results are obtained for artificially generated networks as well as two real ones, namely, the Routers of the Internet and a neuronal network.

preprint2012arXiv

Klein Tunneling in the presence of random impurities

In this paper, we study Klein tunneling in random media. To this purpose, we simulate the propagation of a relativistic Gaussian wavepacket through a graphene sample with randomly distributed potential barriers (impurities). The simulations, based on a relativistic quantum lattice Boltzmann method, permit to compute the transmission coefficient across the sample, thereby providing an estimate for the conductivity as a function of impurity concentration and strength of the potentials. It is found that the conductivity loss due to impurities is significantly higher for wave-packets of massive particles, as compared to massless ones. A general expression for the loss of conductivity as a function of the impurity percentage is presented and successfully compared with the Kozeny-Carman law for disordered media in classical fluid dynamics.

preprint2012arXiv

Overspill avalanching in a dense reservoir network

Sustainability of communities, agriculture, and industry is strongly dependent on an effective storage and supply of water resources. In some regions the economic growth has led to a level of water demand which can only be accomplished through efficient reservoir networks. Such infrastructures are not always planned at larger scale but rather made by farmers according to their local needs of irrigation during droughts. Based on extensive data from the upper Jaguaribe basin, one of the world's largest system of reservoirs, located in the Brazilian semiarid northeast, we reveal that surprisingly it self-organizes into a scale-free network exhibiting also a power-law in the distribution of the lakes and avalanches of discharges. With a new self-organized-criticality-type model we manage to explain the novel critical exponents. Implementing a flow model we are able to reproduce the measured overspill evolution providing a tool for catastrophe mitigation and future planning.

preprint2012arXiv

Pre-Turbulent Regimes in Graphene Flows

We provide numerical evidence that electronic pre-turbulent phenomena in graphene could be observed, under current experimental conditions, through detectable current fluctuations, echoing the detachment of vortices past localized micron-sized impurities. Vortex generation, due to micron sized constriction, is also explored with special focus on the effects of relativistic corrections to the normal Navier-Stokes equations. These corrections are found to cause a delay in the stability breakout of the fluid as well as a small shift in the vortex shedding frequency. Finally, a relation between the Strouhal number, a dimensionless measure of the vortex shedding frequency, and the Reynolds number is provided under conditions of interest for future experiments.

preprint2012arXiv

Single to Double Hump Transition in the Equilibrium Distribution Function of Relativistic Particles

We unveil a transition from single peaked to bimodal velocity distribution in a relativistic fluid under increasing temperature, in contrast with a non-relativistic gas, where only a monotonic broadening of the bell-shaped distribution is observed. Such transition results from the interplay between the raise in thermal energy and the constraint of maximum velocity imposed by the speed of light. We study the Bose-Einstein, the Fermi-Dirac, and the Maxwell-Jüttner distributions, all exhibiting the same qualitative behavior. We characterize the nature of the transition in the framework of critical phenomena and show that it is either continuous or discontinuous, depending on the group velocity. We analyze the transition in one, two, and three dimensions, with special emphasis on two-dimensions, for which a possible experiment in graphene, based on the measurement of the Johnson-Nyquist noise, is proposed.

preprint2012arXiv

Taylor-Couette Instability in General Manifolds: A Lattice Kinetic Approach

We present a new lattice kinetic method to simulate fluid dynamics in curvilinear geometries. A suitable discrete Boltzmann equation is solved in contravariant coordinates, and the equilibrium distribution function is obtained by a Hermite polynomials expansion of the Maxwell-Boltzmann distribution, expressed in terms of the contravariant coordinates and the metric tensor. To validate the model, we calculate the critical Reynolds number for the onset of the Taylor-Couette instability between two concentric cylinders, obtaining excellent agreement with the theory. In order to extend this study to more general geometries, we also calculate the critical Reynolds number for the case of two concentric spheres, finding good agreement with experimental data. In the case of two concentric tori, we have found that the critical Reynolds is about 10% larger than the respective value for the two concentric cylinders.

preprint2012arXiv

The balance between excitation and inhibition controls the temporal organization of neuronal avalanches

Neuronal avalanches, measured in vitro and in vivo, exhibit a robust critical behaviour. Their temporal organization hides the presence of correlations. Here we present experimental measurements of the waiting time distribution between successive avalanches in the rat cortex in vitro. This exhibits a non-monotonic behaviour, not usually found in other natural processes. Numerical simulations provide evidence that this behaviour is a consequence of the alternation between states of high and low activity, named up and down states, leading to a balance between excitation and inhibition controlled by a single parameter. During these periods both the single neuron state and the network excitability level, keeping memory of past activity, are tuned by homeostatic mechanisms.

preprint2012arXiv

Watersheds are Schramm-Loewner Evolution curves

We show that in the continuum limit watersheds dividing drainage basins are Schramm-Loewner Evolution (SLE) curves, being described by one single parameter $κ$. Several numerical evaluations are applied to ascertain this. All calculations are consistent with SLE$_κ$, with $κ=1.734\pm0.005$, being the only known physical example of an SLE with $κ<2$. This lies outside the well-known duality conjecture, bringing up new questions regarding the existence and reversibility of dual models. Furthermore it constitutes a strong indication for conformal invariance in random landscapes and suggests that watersheds likely correspond to a logarithmic Conformal Field Theory (CFT) with central charge $c\approx-7/2$.

preprint2011arXiv

Fractality of eroded coastlines of correlated landscapes

Using numerical simulations of a simple sea-coast mechanical erosion model, we investigate the effect of spatial long-range correlations in the lithology of coastal landscapes on the fractal behavior of the corresponding coastlines. In the model, the resistance of a coast section to erosion depends on the local lithology configuration as well as on the number of neighboring sea sides. For weak sea forces, the sea is trapped by the coastline and the eroding process stops after some time. For strong sea forces erosion is perpetual. The transition between these two regimes takes place at a critical sea force, characterized by a fractal coastline front. For uncorrelated landscapes, we obtain, at the critical value, a fractal dimension D=1.33, which is consistent with the dimension of the accessible external perimeter of the spanning cluster in two-dimensional percolation. For sea forces above the critical value, our results indicate that the coastline is self-affine and belongs to the Kardar-Parisi-Zhang universality class. In the case of landscapes generated with power-law spatial long-range correlations, the coastline fractal dimension changes continuously with the Hurst exponent H, decreasing from D=1.34 to 1.04, for H=0 and 1, respectively. This nonuniversal behavior is compatible with the multitude of fractal dimensions found for real coastlines.

preprint2011arXiv

Gaussian model of explosive percolation in three and higher dimensions

The Gaussian model of discontinuous percolation, recently introduced by Araújo and Herrmann [Phys. Rev. Lett., 105, 035701 (2010)], is numerically investigated in three dimensions, disclosing a discontinuous transition. For the simple-cubic lattice, in the thermodynamic limit, we report a finite jump of the order parameter, $J=0.415 \pm 0.005$. The largest cluster at the threshold is compact, but its external perimeter is fractal with fractal dimension $d_A = 2.5 \pm 0.2$. The study is extended to hypercubic lattices up to six dimensions and to the mean-field limit (infinite dimension). We find that, in all considered dimensions, the percolation transition is discontinuous. The value of the jump in the order parameter, the maximum of the second moment, and the percolation threshold are analyzed, revealing interesting features of the transition and corroborating its discontinuous nature in all considered dimensions. We also show that the fractal dimension of the external perimeter, for any dimension, is consistent with the one from bridge percolation and establish a lower bound for the percolation threshold of discontinuous models with finite number of clusters at the threshold.

preprint2011arXiv

Granular packings of cohesive elongated particles

We report numerical results of effective attractive forces on the packing properties of two-dimensional elongated grains. In deposits of non-cohesive rods in 2D, the topology of the packing is mainly dominated by the formation of ordered structures of aligned rods. Elongated particles tend to align horizontally and the stress is mainly transmitted from top to bottom, revealing an asymmetric distribution of local stress. However, for deposits of cohesive particles, the preferred horizontal orientation disappears. Very elongated particles with strong attractive forces form extremely loose structures, characterized by an orientation distribution, which tends to a uniform behavior when increasing the Bond number. As a result of these changes, the pressure distribution in the deposits changes qualitatively. The isotropic part of the local stress is notably enhanced with respect to the deviatoric part, which is related to the gravity direction. Consequently, the lateral stress transmission is dominated by the enhanced disorder and leads to a faster pressure saturation with depth.

preprint2011arXiv

Impact of Perturbations on Watersheds

We find that watersheds in real and artificial landscapes can be strongly affected by small, local perturbations like landslides or tectonic motions. We observe power-law scaling behavior for both the distribution of areas enclosed by the original and the displaced watershed as well as the probability density to induce, after perturbation, a change at a given distance. Scaling exponents for real and artificial landscapes are determined, where in the latter case the exponents depend linearly on the Hurst exponent of the applied fractional Brownian noise. The obtained power-laws are shown to be independent on the strength of perturbation. Theoretical arguments relate our scaling laws for uncorrelated landscapes to properties of invasion percolation.

preprint2011arXiv

Jump at the onset of saltation

We reveal a discontinuous transition in the saturated flux for aeolian saltation by simulating explicitly particle motion in turbulent flow. The discontinuity is followed by a coexistence interval with two metastable solutions. The modification of the wind profile due to momentum exchange exhibits a second maximum at high shear strength. The saturated flux depends on the strength of the wind as $q_s=q_0+A(u_*-u_t)(u_*^2+u_t^2)$.

preprint2011arXiv

Modelling the Hyphal Growth of the Wood-decay Fungus Physisporinus vitreus

The white-rot fungus, Physisporinus vitreus, degrades the membranes of bordered pits in tracheids and consequently increases the permeability of wood, which is a process that can be used by the wood industry to improve the uptake of wood preservatives and environmentally benign wood modification substances to enhance the use and sustainability of native conifer wood species. To understand and apply this process requires an understanding of how a complex system (fungus-wood) interacts under defined conditions. We present a three-dimensional fungal growth model (FGM) of the hyphal growth of P. vitreus in the heartwood of Norway spruce. The model considers hyphae and nutrients as discrete structures and links the microscopic interactions between fungus and wood (e.g. degradation rate and degree of opening of pits) with macroscopic system properties, such penetration depth of the fungus, biomass and distribution of destroyed pits in early- and latewood. Simulations were compared with experimental data. The growth of P. vitreus is characterized by a stepwise capture of the substrate and the effect of this on wood according to different model parameters is discussed.

preprint2011arXiv

Oil displacement through a porous medium with a temperature gradient

We investigate the effect of a temperature gradient on oil recovery in a two-dimensional pore-network model. The oil viscosity depends on temperature as, $μ_o=exp(B/T)$, where $B$ is a physico-chemical parameter depending on the type of oil, and $T$ is the temperature. A temperature gradient is applied across the medium in the flow direction. Initially, the porous medium is saturated with oil and, then, another fluid is injected. We have considered two cases representing different injection strategies. In the first case, the invading fluid viscosity is constant (finite viscosity ratio) while in the second one, the invading fluid is inviscid (infinite viscosity ratio). Our results show that, for the case of finite viscosity ratio, recovery increases with $ΔT$ independently on strength or sign of the gradient. For an infinite viscosity ratio, a positive temperature gradient is necessary to enhance recovery. Moreover, we show that, for $ΔT>0$, the percentage of oil recovery generally decreases (increases) with $B$ for a finite (infinite) viscosity ratio. Finally, we also extend our results for infinite viscosity ratio to a three-dimensional porous media geometry.

preprint2011arXiv

On the relativistic Lattice Boltzmann method for quark-gluon plasma simulations

In this paper, we investigate the recently developed lattice Boltzmann model for relativistic hydrodynamics. To this purpose, we perform simulations of shock waves in quark-gluon plasma in the low and high viscosities regime, using three different computational models, the relativistic lattice Boltzmann (RLB), the Boltzmann Approach Multi-Parton Scattering (BAMPS), and the viscous sharp and smooth transport algorithm (vSHASTA). From the results, we conclude that the RLB model departs from BAMPS in the case of high speeds and high temperature(viscosities), the departure being due to the fact that the RLB is based on a quadratic approximation of the Maxwell-Jüttner distribution, which is only valid for sufficiently low temperature and velocity. Furthermore, we have investigated the influence of the lattice speed on the results, and shown that inclusion of quadratic terms in the equilibrium distribution improves the stability of the method within its domain of applicability. Finally, we assess the viability of the RLB model in the various parameter regimes relevant to ultra-relativistic fluid dynamics.

preprint2011arXiv

Optimal path cracks in correlated and uncorrelated lattices

The optimal path crack model on uncorrelated surfaces, recently introduced by Andrade et al. (Phys. Rev. Lett. 103, 225503, 2009), is studied in detail and its main percolation exponents computed. In addition to beta/nu = 0.46 \pm 0.03 we report, for the first time, gamma/nu = 1.3 \pm 0.2 and tau = 2.3 \pm 0.2. The analysis is extended to surfaces with spatial long-range power-law correlations, where non-universal fractal dimensions are obtained when the degree of correlation is varied. The model is also considered on a three-dimensional lattice, where the main crack is found to be a surface with a fractal dimension of 2.46 \pm 0.05.

preprint2011arXiv

Packing of elastic wires in spherical cavities

We investigate the morphologies and maximum packing density of thin wires packed into spherical cavities. Using simulations and experiments, we find that ordered as well as disordered structures emerge, depending on the amount of internal torsion. We find that the highest packing densities are achieved in low torsion packings for large systems, but in high torsion packings for small systems. An analysis of both situations is given in terms of energetics and comparison is made to analytical models of DNA packing in viral capsids.

preprint2011arXiv

Scaling Relations for Watersheds

We study the morphology of watersheds in two and three dimensional systems subjected to different degrees of spatial correlations. The response of these objects to small, local perturbations is also investigated with extensive numerical simulations. We find the fractal dimension of the watersheds to generally decrease with the Hurst exponent, which quantifies the degree of spatial correlations. Moreover, in two dimensions, our results match the range of fractal dimensions $1.10 \leq d_f \leq 1.15$ observed for natural landscapes. We report that the watershed is strongly affected by local perturbations. For perturbed two and three dimensional systems, we observe a power-law scaling behavior for the distribution of areas (volumes) enclosed by the original and the displaced watershed, and for the distribution of distances between outlets. Finite-size effects are analyzed and the resulting scaling exponents are shown to depend significantly on the Hurst exponent. The intrinsic relation between watershed and invasion percolation, as well as relations between exponents conjectured in previous studies with two dimensional systems, are now confirmed by our results in three dimensions.

preprint2011arXiv

Study of the combined effect of temperature, pH and water activity on the radial growth rate of the white-rot basidiomycete Physisporinus vitreus by using a hyphal growth model

The present work investigates environmental effects on the growth of fungal colonies of P. vitreus by using a lattice-free discrete modelling approach called FGM (Fuhr et al. (2010), arXiv:1101.1747), in which hyphae and nutrients are considered as discrete structures. A discrete modelling approach allows studying the underlying mechanistic rule concerning the basic architecture and dynamic of fungal networks on the scale of a single colony. By comparing simulations of the FGM with laboratory experiments of growing fungal colonies on malt extract agar we show that combined effect of temperature, pH and water activity on the radial growth rate of a fungal colony on a macroscopic scale may be explained by a power law for the growth costs of hyphal expansion on a microscopic scale. The information about the response of the fungal mycelium on a microscopic scale to environmental conditions is essential to simulate its behavior in complex structure substrates such as wood, where the impact of the fungus to the wood (i.e. the degradation of bordered pits or the creation of bore holes and cavities) changes the local environmental condition, e.g. the permeability of the substrate and therefore the water activity level in a colonized wood cell lumen. A combination of diffusion and moisture processes into wood with the FGM may brighten the knowledge about the colonization strategy of P. vitreus and helps to optimize its growth behavior for biotechnological application such as bioincising.

preprint2010arXiv

Adhesive Penetration in Beech Wood Part I: Experiments

A study with synchrotron radiation X-ray tomographic microscopy (SRXTM) of PUR, PVAC, and UF adhesive bond lines in beech wood, bonded under various growth ring angles is presented. After determining the hardening characteristics of the adhesives, we evaluate the bond line morphologies, and the adhesive penetration into the wood structure. We find distinct bond line imperfections for the different adhesive systems. To describe the adhesive distribution inside the bond line we propose the saturation of the pore space instead of the commonly used maximum penetration depth. The results are the basis for a penetration model of hardening fluids into hardwood (part II).

preprint2010arXiv

Adhesive penetration in Beech wood Part II: Penetration Model

We propose an analytical model to predict the adhesives penetration into hard wood. Penetration of hard wood is dominated by the vessel network which prohibits porous medium approximations. Our model considers two scales: a one dimensional capillary fluid transport of a hardening adhesive through a single, straight vessel with diffusion of solvent through its walls and a mesoscopic scale based on topological characteristics of the vessel network, where results from the single vessel scale are mapped onto a periodic network. Given an initial amount of adhesive and applied bonding pressure, we calculate the portion of the filled structure. The model is applied to beech wood samples joined with three different types of adhesive (PUR, UF, PVAC) under various growth ring angles. We evaluate adhesive properties and bond line morphologies described in part I of this work. The model contains one free parameter that can be adjusted in order to fit the experimental data.

preprint2010arXiv

Derivation of the Lattice Boltzmann Model for Relativistic Hydrodynamics

A detailed derivation of the Lattice Boltzmann (LB) scheme for relativistic fluids recently proposed in Ref. [1], is presented. The method is numerically validated and applied to the case of two quite different relativistic fluid dynamic problems, namely shock-wave propagation in quark-gluon plasmas and the impact of a supernova blast-wave on massive interstellar clouds. Close to second order convergence with the grid resolution, as well as linear dependence of computational time on the number of grid points and time-steps, are reported.

preprint2010arXiv

Explosive percolation via control of the largest cluster

We show that only considering the largest cluster suffices to obtain a first-order percolation transition. As opposed to previous realizations of explosive percolation our models obtain Gaussian cluster distributions and compact clusters as one would expect at first-order transitions. We also discover that the cluster perimeters are fractal at the transition point, yielding a fractal dimension of $1.23\pm0.03$, close to that of watersheds.

preprint2010arXiv

New universality class for the fragmentation of plastic materials

We present an experimental and theoretical study of the fragmentation of polymeric materials by impacting polypropylene particles of spherical shape against a hard wall. Experiments reveal a power law mass distribution of fragments with an exponent close to 1.2, which is significantly different from the known exponents of three-dimensional bulk materials. A 3D discrete element model is introduced which reproduces both the large permanent deformation of the polymer during impact, and the novel value of the mass distribution exponent. We demonstrate that the dominance of shear in the crack formation and the plastic response of the material are the key features which give rise to the emergence of the novel universality class of fragmentation phenomena.

preprint2009arXiv

A Hamiltonian approach for explosive percolation

We introduce a cluster growth process that provides a clear connection between equilibrium statistical mechanics and an explosive percolation model similar to the one recently proposed by Achlioptas et al. [Science 323, 1453 (2009)]. We show that the following two ingredients are essential for obtaining an abrupt (first-order) transition in the fraction of the system occupied by the largest cluster: (i) the size of all growing clusters should be kept approximately the same, and (ii) the inclusion of merging bonds (i.e., bonds connecting vertices in different clusters) should dominate with respect to the redundant bonds (i.e., bonds connecting vertices in the same cluster). Moreover, in the extreme limit where only merging bonds are present, a complete enumeration scheme based on tree-like graphs can be used to obtain an exact solution of our model that displays a first-order transition. Finally, the proposed mechanism can be viewed as a generalization of standard percolation that discloses an entirely new family of models with potential application in growth and fragmentation processes of real network systems.

preprint2009arXiv

Fracturing the optimal paths

Optimal paths play a fundamental role in numerous physical applications ranging from random polymers to brittle fracture, from the flow through porous media to information propagation. Here for the first time we explore the path that is activated once this optimal path fails and what happens when this new path also fails and so on, until the system is completely disconnected. In fact numerous applications can be found for this novel fracture problem. In the limit of strong disorder, our results show that all the cracks are located on a single self-similar connected line of fractal dimension $D_{b} \approx 1.22$. For weak disorder, the number of cracks spreads all over the entire network before global connectivity is lost. Strikingly, the disconnecting path (backbone) is, however, completely independent on the disorder.

preprint2009arXiv

Lattice Boltzmann scheme for relativistic fluids

A Lattice Boltzmann formulation for relativistic fluids is presented and numerically verified through quantitative comparison with recent hydrodynamic simulations of relativistic shock-wave propagation in viscous quark-gluon plasmas. This formulation opens up the possibility of exporting the main advantages of Lattice Boltzmann methods to the relativistic context, which seems particularly useful for the simulation of relativistic fluids in complicated geometries.

preprint2009arXiv

Simulation of Flow of Mixtures Through Anisotropic Porous Media using a Lattice Boltzmann Model

We propose a description for transient penetration simulations of miscible and immiscible fluid mixtures into anisotropic porous media, using the lattice Boltzmann (LB) method. Our model incorporates hydrodynamic flow, diffusion, surface tension, and the possibility for global and local viscosity variations to consider various types of hardening fluids. The miscible mixture consists of two fluids, one governed by the hydrodynamic equations and one by diffusion equations. We validate our model on standard problems like Poiseuille flow, the collision of a drop with an impermeable, hydrophobic interface and the deformation of the fluid due to surface tension forces. To demonstrate the applicability to complex geometries, we simulate the invasion process of mixtures into wood spruce samples.

preprint2008arXiv

Universality behind Basquin's law of fatigue

One of the most important scaling laws of time dependent fracture is Basquin's law of fatigue, namely, that the lifetime of the system increases as a power law with decreasing external load amplitude, $t_f\sim σ_0^{-α}$, where the exponent $α$ has a strong material dependence. We show that in spite of the broad scatter of the Basquin exponent $α$, the fatigue fracture of heterogeneous materials exhibits intriguing universal features. Based on stochastic fracture models we propose a generic scaling form for the macroscopic deformation and show that at the fatigue limit the system undergoes a continuous phase transition when changing the external load. On the microlevel, the fatigue fracture proceeds in bursts characterized by universal power law distributions. We demonstrate that in a range of systems, including deformation of asphalt, a realistic model of deformation, and a fiber bundle model, the system dependent details are contained in Basquin's exponent for time to failure, and once this is taken into account, remaining features of failure are universal.

preprint2007arXiv

Parabolic dunes in north-eastern Brazil

In this work we present measurements of vegetation cover over parabolic dunes with different degree of activation along the north-eastern Brazilian coast. We are able to extend the local values of the vegetation cover density to the whole dune by correlating measurements with the gray-scale levels of a high resolution satellite image of the dune field. The empirical vegetation distribution is finally used to validate the results of a recent continuous model of dune motion coupling sand erosion and vegetation growth.

preprint2006arXiv

Community structure and ethnic preferences in school friendship networks

Recently developed concepts and techniques of analyzing complex systems provide new insight into the structure of social networks. Uncovering recurrent preferences and organizational principles in such networks is a key issue to characterize them. We investigate school friendship networks from the Add Health database. Applying threshold analysis, we find that the friendship networks do not form a single connected component through mutual strong nominations within a school, while under weaker conditions such interconnectedness is present. We extract the networks of overlapping communities at the schools (c-networks) and find that they are scale free and disassortative in contrast to the direct friendship networks, which have an exponential degree distribution and are assortative. Based on the network analysis we study the ethnic preferences in friendship selection. The clique percolation method we use reveals that when in minority, the students tend to build more densely interconnected groups of friends. We also find an asymmetry in the behavior of black minorities in a white majority as compared to that of white minorities in a black majority.

preprint2006arXiv

Influence of particle shape on sheared dense granular media

We study by means of molecular dynamics simulations of periodic shear cells, the influence of particle shape on the global mechanical behavior of dense granular media. Results at macro-mechanical level show that for large shear deformation samples with elongated particles, independent of their initial orientation, reach the same stationary value for both shear force and void ratio. At the micro-mechanical level the stress, the fabric and the inertia tensors of the particles are used to study the evolution of the media. In the case of isotropic particles the direction of the principal axis of the fabric tensor is aligned with the one of the principal stress, while for elongated particles the fabric orientation is strongly dependent on the orientation of the particles. The shear band width is shown to depend on the particle shape due to the tendency of elongated particles to preferential orientations and less rotation.

preprint2005arXiv

A simple beam model for the shear failure of interfaces

We propose a novel model for the shear failure of a glued interface between two solid blocks. We model the interface as an array of elastic beams which experience stretching and bending under shear load and break if the two deformation modes exceed randomly distributed breaking thresholds. The two breaking modes can be independent or combined in the form of a von Mises type breaking criterion. Assuming global load sharing following the beam breaking, we obtain analytically the macroscopic constitutive behavior of the system and describe the microscopic process of the progressive failure of the interface. We work out an efficient simulation technique which allows for the study of large systems. The limiting case of very localized interaction of surface elements is explored by computer simulations.

preprint2005arXiv

Characterization of the material response in the granular ratcheting

The existence of a very special ratcheting regime has recently been reported in a granular packing subjected to cyclic loading \cite{alonso04}. In this state, the system accumulates a small permanent deformation after each cycle. After a short transient regime, the value of this permanent strain accumulation becomes independent on the number of cycles. We show that a characterization of the material response in this peculiar state is possible in terms of three simple macroscopic variables. They are defined that, they can be easily measured both in the experiments and in the simulations. We have carried out a thorough investigation of the micro- and macro-mechanical factors affecting these variables, by means of Molecular Dynamics simulations of a polydisperse disk packing, as a simple model system for granular material. Biaxial test boundary conditions with a periodically cycling load were implemented. The effect on the plastic response of the confining pressure, the deviatoric stress and the number of cycles has been investigated. The stiffness of the contacts and friction has been shown to play an important role in the overall response of the system. Specially elucidating is the influence of the particular hysteretical behavior in the stress-strain space on the accumulation of permanent strain and the energy dissipation.

preprint2005arXiv

Renormalizing Sznajd model on complex networks taking into account the effects of growth mechanisms

We present a renormalization approach to solve the Sznajd opinion formation model on complex networks. For the case of two opinions, we present an expression of the probability of reaching consensus for a given opinion as a function of the initial fraction of agents with that opinion. The calculations reproduce the sharp transition of the model on a fixed network, as well as the recently observed smooth function for the model when simulated on a growing complex networks.

preprint2005arXiv

Scaling behavior of fragment shapes

We present an experimental and theoretical study of the shape of fragments generated by explosive and impact loading of closed shells. Based on high speed imaging, we have determined the fragmentation mechanism of shells. Experiments have shown that the fragments vary from completely isotropic to highly anisotropic elongated shapes, depending on the microscopic cracking mechanism of the shell-material. Anisotropic fragments proved to have self-affine character described by a scaling exponent. The distribution of fragment shapes exhibits a power law decay. The robustness of the scaling laws is illustrated by a stochastic hierarchical model of fragmentation. Our results provide a possible improvement of the representation of fragment shapes in models of space debris.

preprint2004arXiv

Break-up of shells under explosion and impact

A theoretical and experimental study of the fragmentation of closed thin shells made of a disordered brittle material is presented. Experiments were performed on brown and white hen egg-shells under two different loading conditions: fragmentation due to an impact with a hard wall and explosion by a combustion mixture giving rise to power law fragment size distributions. For the theoretical investigations a three-dimensional discrete element model of shells is constructed. Molecular dynamics simulations of the two loading cases resulted in power law fragment mass distributions in satisfactory agreement with experiments. Based on large scale simulations we give evidence that power law distributions arise due to an underlying phase transition which proved to be abrupt and continuous for explosion and impact, respectively. Our results demonstrate that the fragmentation of closed shells defines a universality class different from that of two- and three-dimensional bulk systems.

preprint2004arXiv

Investigation of the incremental response of soils using a discrete element model

The incremental stress-strain relation of dense packings of polygons is investigated here by using molecular dynamics simulations. The comparison of the simulation results to the continuous theories is performed using explicit expressions for the averaged stress and strain over a representative volume element. The discussion of the incremental response raises two important questions of soil deformation: Is the incrementally non-linear theory appropriate to describe the soil mechanical response? Does a purely elastic regime exists in the deformation of granular materials?. In both cases our answer will be "no". The question of stability is also discussed in terms of the Hill condition of stability for non-associated materials.

preprint2004arXiv

Shakedown of unbound granular material

Compacted unbound granular materials are extensively used as sub-layer in pavement design. Most pavement design guides assume that they are responsible for the degradation and deformation of the roads and railways that they support. Biaxial tests are usually employed to investigate the elasto-plastic response of these materials to cyclic loading. A particularly interesting question is whether a limit load exists, below which the excitations shake down, in the sense that the material does not accumulate further deformations. We have carried out a detailed study of the elasto-plastic behavior of a simple model of unbound granular matter submitted to cyclic loading. The dissipated energy through out the simulation has been used for the characterization of the different regimes of responses.

preprint1998arXiv

Contact forces in regular 3D granular pile

We present exact results for the contact forces in a three dimensional static piling of identical, stiff and frictionless spheres. The pile studied is a pyramid of equilateral triangular base (``stack of cannonballs'') with a FCC (face centered cubic) structure. We show in particular that, as for the two dimensional case, the pressure on the base of such a pile is uniform.

preprint1998arXiv

Shock propagation in a granular chain

We numerically solve the propagation of a shock wave in a chain of elastic beads with no restoring forces under traction (no-tension elasticity). We find a sequence of peaks of decreasing amplitude and velocity. Analyzing the main peak at different times we confirm a recently proposed scaling law for its decay.

preprint1996arXiv

Fragmentation of Colliding Discs

We study the phenomena associated with the low-velocity impact of two solid discs of equal size using a cell model of brittle solids. The fragment ejection exhibits a jet-like structure the direction of which depends on the impact parameter. We obtain the velocity and the mass distribution of the debris. Varying the radius and the initial velocity of the colliding particles, the velocity components of the fragments show anomalous scaling. The mass distribution follows a power law in the region of intermediate masses.

preprint1995arXiv

On Hydrodynamic Diffusion and Velocity Fluctuations in Two--Dimensional Simulations of Sedimentation

We present a numerical method to deal efficiently with large numbers of particles in incompressible fluids. The interactions between particles and fluid are taken into account by a physically motivated ansatz based on locally defined drag forces. We demonstrate the validity of our approach by performing numerical simulations of sedimenting non-Brownian spheres in two spatial dimensions and compare our results with experiments. Our method reproduces essential aspects of the experimental findings, in particular the strong anisotropy of the hydrodynamic bulk self-diffusivities.

H. J. Herrmann

What is connected

Connect this record

See the researcher in context

Building this map preview

106 published item(s)

Dynamics of extended Schelling models

Frustrated bearings

Geometry-induced non-equilibrium phase transition in sandpiles

The elastic backbone phase transition in the Ising model

Conformal Invariance of Graphene Sheets

Connectivity disruption sparks explosive epidemic spreading

Critical fragmentation properties of random drilling: How many random holes need to be drilled to collapse a wooden cube?

Crumpling Damaged Graphene

Cutting self-similar space-filling sphere packings

Emergence of core-peripheries in networks

Energy dissipation in flows through curved spaces

Field Driven Charging Dynamics of a Fluidized Granular Bed

Gender gap in the ERASMUS mobility program

Hybrid Percolation Transition in Cluster Merging Processes: Continuously Varying Exponents

Model for the growth of the World Airline Network

Poiseuille flow in curved spaces

Prediction and control of slip-free rotation states in sphere assemblies

Self-similar space-filling sphere packings in three and four dimensions

Synaptic plasticity and neuronal refractory time cause scaling behaviour of neuronal avalanches

Analysis of the velocity field of granular hopper flow

Bursts in discontinuous Aeolian saltation

Complex networks from space-filling bearings

Cooling Effect of the Richtmyer-Meshkov Instability

Coupled DEM-LBM method for the free-surface simulation of heterogeneous suspensions

Critical Cooperation Range to Improve Spatial Network Robustness

En route to a multi-model scheme for clinker comminution with chemical grinding aids

Flow through randomly curved manifolds

Granular front formation in free-surface flow of concentrated suspensions

High-Order Kinetic Relaxation Schemes as High-Accuracy Poisson Solvers

Lattice Boltzmann Model for Electronic Structure Simulations

Lattice Boltzmann model for resistive relativistic magnetohydrodynamics

Schramm-Loewner Evolution and isoheight lines of correlated landscapes

Stochastic Loewner Evolution Relates Anomalous Diffusion and Anisotropic Percolation

Dean Instability in Double-Curved Channels

Flood avalanches in a semiarid basin with a dense reservoir network

Mandala Networks: ultra-robust, ultra-small-world and highly sparse graphs

Modelling the influence of photospheric turbulence on solar flare statistics

Retention capacity of correlated surfaces

Shortest path and Schramm-Loewner Evolution

Usage leading to an abrupt collapse of connectivity

Watersheds in disordered media

Avoiding a Spanning Cluster in Percolation Models

Breathing synchronization in interconnected networks

Exact evaluation of the cutting path length in a percolation model on a hierarchical network

Formal analogy between the Dirac equation in its Majorana form and the discrete-velocity version of the Boltzmann kinetic equation

Gaussian Quadrature and Lattice Discretization of the Fermi-Dirac Distribution for Graphene

Geometrical and transport properties of sequential adsorption clusters

Hydrodynamic Model for Conductivity in Graphene

Kinetic Formulation of the Kohn-Sham Equations for ab initio Electronic Structure Calculations

Lattice Boltzmann model for ultra-relativistic flows

Optimal synchronizability of bearings

Percolation with long-range correlated disorder

Relativistic effects on the Richtmyer-Meshkov instability

Relativistic Lattice Boltzmann Model with Improved Dissipation

Smart Rewiring for Network Robustness

Stacked triangular lattice: Percolation properties

Ultrarelativistic Transport Coefficients in Two Dimensions

Bohman-Frieze-Wormald model on the lattice, yielding a discontinuous percolation transition

Corrections to Scaling for Watersheds, Optimal Path Cracks, and Bridge Lines

Fracturing highly disordered materials

Fracturing ranked surfaces

How to share underground reservoirs

How to suppress undesired synchronization

Klein Tunneling in the presence of random impurities

Overspill avalanching in a dense reservoir network

Pre-Turbulent Regimes in Graphene Flows

Single to Double Hump Transition in the Equilibrium Distribution Function of Relativistic Particles

Taylor-Couette Instability in General Manifolds: A Lattice Kinetic Approach

The balance between excitation and inhibition controls the temporal organization of neuronal avalanches

Watersheds are Schramm-Loewner Evolution curves

Fractality of eroded coastlines of correlated landscapes

Gaussian model of explosive percolation in three and higher dimensions

Granular packings of cohesive elongated particles

Impact of Perturbations on Watersheds