Sampling properties of directed networks

For many real-world networks only a small "sampled" version of the original network may be investigated; those results are then used to draw conclusions about the actual system. Variants of breadth-first search (BFS) sampling, which are based on epidemic processes, are widely used. Although it is well established that BFS sampling fails, in most cases, to capture the IN-component(s) of directed networks, a description of the effects of BFS sampling on other topological properties are all but absent from the literature. To systematically study the effects of sampling biases on directed networks, we compare BFS sampling to random sampling on complete large-scale directed networks. We present new results and a thorough analysis of the topological properties of seven different complete directed networks (prior to sampling), including three versions of Wikipedia, three different sources of sampled World Wide Web data, and an Internet-based social network. We detail the differences that sampling method and coverage can make to the structural properties of sampled versions of these seven networks. Most notably, we find that sampling method and coverage affect both the bow-tie structure, as well as the number and structure of strongly connected components in sampled networks. In addition, at low sampling coverage (i.e. less than 40%), the values of average degree, variance of out-degree, degree auto-correlation, and link reciprocity are overestimated by 30% or more in BFS-sampled networks, and only attain values within 10% of the corresponding values in the complete networks when sampling coverage is in excess of 65%. These results may cause us to rethink what we know about the structure, function, and evolution of real-world directed networks.Comment: 21 pages, 11 figure

A Novel Approach to Discontinuous Bond Percolation Transition

We introduce a bond percolation procedure on a $D$-dimensional lattice where two neighbouring sites are connected by $N$ channels, each operated by valves at both ends. Out of a total of $N$, randomly chosen $n$ valves are open at every site. A bond is said to connect two sites if there is at least one channel between them, which has open valves at both ends. We show analytically that in all spatial dimensions, this system undergoes a discontinuous percolation transition in the $N\to \infty$ limit when $\gamma =\frac{\ln n}{\ln N}$ crosses a threshold. It must be emphasized that, in contrast to the ordinary percolation models, here the transition occurs even in one dimensional systems, albeit discontinuously. We also show that a special kind of discontinuous percolation occurs only in one dimension when $N$ depends on the system size.Comment: 6 pages, 6 eps figure

Random elastic networks : strong disorder renormalization approach

For arbitrary networks of random masses connected by random springs, we define a general strong disorder real-space renormalization (RG) approach that generalizes the procedures introduced previously by Hastings [Phys. Rev. Lett. 90, 148702 (2003)] and by Amir, Oreg and Imry [Phys. Rev. Lett. 105, 070601 (2010)] respectively. The principle is to eliminate iteratively the elementary oscillating mode of highest frequency associated with either a mass or a spring constant. To explain the accuracy of the strong disorder RG rules, we compare with the Aoki RG rules that are exact at fixed frequency.Comment: 8 pages, v2=final versio

Two-dimensional SIR epidemics with long range infection

We extend a recent study of susceptible-infected-removed epidemic processes with long range infection (referred to as I in the following) from 1-dimensional lattices to lattices in two dimensions. As in I we use hashing to simulate very large lattices for which finite size effects can be neglected, in spite of the assumed power law $p({\bf x})\sim |{\bf x}|^{-\sigma-2}$ for the probability that a site can infect another site a distance vector ${\bf x}$ apart. As in I we present detailed results for the critical case, for the supercritical case with $\sigma = 2$, and for the supercritical case with $0< \sigma < 2$. For the latter we verify the stretched exponential growth of the infected cluster with time predicted by M. Biskup. For $\sigma=2$ we find generic power laws with $\sigma-$dependent exponents in the supercritical phase, but no Kosterlitz-Thouless (KT) like critical point as in 1-d. Instead of diverging exponentially with the distance from the critical point, the correlation length increases with an inverse power, as in an ordinary critical point. Finally we study the dependence of the critical exponents on $\sigma$ in the regime $0<\sigma <2$, and compare with field theoretic predictions. In particular we discuss in detail whether the critical behavior for $\sigma$ slightly less than 2 is in the short range universality class, as conjectured recently by F. Linder {\it et al.}. As in I we also consider a modified version of the model where only some of the contacts are long range, the others being between nearest neighbors. If the number of the latter reaches the percolation threshold, the critical behavior is changed but the supercritical behavior stays qualitatively the same.Comment: 14 pages, including 29 figure

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

Exact solutions for mass-dependent irreversible aggregations

We consider the mass-dependent aggregation process (k + 1) X -> X, given a fixed number of unit mass particles in the initial state. One cluster is chosen proportional to its mass and is merged into one, either with k neighbors in one dimension, or-in the well-mixed case-with k other clusters picked randomly. We find the same combinatorial exact solutions for the probability to find any given configuration of particles on a ring or line, and in the well-mixed case. The mass distribution of a single cluster exhibits scaling laws and the finite-size scaling form is given. The relation to the classical sum kernel of irreversible aggregation is discussed

Hole Cleaning Performance of Water vs. Polymer-Based Fluids Under Turbulent Flow Conditions

Abstract The problem of solid cleanout in horizontal wellbores is studied experimentally. The special case of drilling fluid circulation with no inner pipe rotation is considered. This case is similar to Coiled tubing drilling where often hole cleanout must be performed. Sand sized cuttings (ranging from 260 micron to 1240 micron) are used. Critical velocity and wall shear stress required for initiating bed erosion are measured. Water and viscous polymer base fluids with 3 different polymer concentrations are used. Results have shown that water always initiate cuttings movement at lower flow rates than polymer solutions. Fluids with higher polymer concentration and hence with higher viscosity required higher flow rates to start eroding the bed. Critical wall shear stress is also determined from pressure loss measurements. Analyzing the data revealed that water initiate cuttings removal at lower pressure loss than more viscous fluids. Higher viscosity fluids always showed higher pressure loss at the initiation of bed erosion. For the range of cuttings size studied, results show that intermediate cuttings size is easier to remove. Smallest and largest cuttings are more difficult to move requiring higher flow rates and higher pressure losses. Dimensionless analysis of relevant parameters to the process of cuttings movement is performed. It is shown that dimensionless wall shear stress (in the forms of Shield's stress and also ratio of shear velocity to settling velocity) correlates well with particles Reynolds number. Based on this finding two correlations are developed to predict wall shear stress required for initiating cuttings movement under different conditions.</jats:p