Source author record

H. E. Stanley

H. E. Stanley 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

29works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

29 published item(s)

preprint2019arXiv

Non-Poisson donation behaviors in virtual worlds

Similar to charitable giving in real world, donation behaviors play an important role in the complex interactions among individuals in virtual worlds. However, it is not clear if the donation process is random or not. We investigate this problem using detailed data from parallel virtual worlds adhered to a massively multiplayer online role-playing game. We find that the inter-donation durations follow power-law-tailed distributions distributed with an average tail exponent close to 1.91, have strong long-range correlations, and possess multifractal features. These findings indicate that the donation process is non-Poissonian, which has potential worth in modeling the complicated individuals behaviors in virtual worlds.

preprint2019arXiv

Scaling relations and finite-size scaling in gravitationally correlated lattice percolation models

In some systems, the connecting probability (and thus the percolation process) between two sites depends on the geometric distance between them. To understand such process, we propose gravitationally correlated percolation models for link-adding networks on the two-dimensional lattice $G$ with two strategies $S_{\rm max}$ and $S_{\rm min}$, to add a link $l_{i,j}$ to connect site $i$ and site $j$ with mass $m_i$ and $m_j$, respectively; $m_i$ and $m_j$ are sizes of the clusters which contain site $i$ and site $j$, respectively. The probability to add the link $l_{i,j}$ is related to the generalized gravity $g_{ij} \equiv m_i m_j/r_{ij}^d$, where $r_{ij}$ is the geometric distance between $i$ and $j$, and $d$ is an adjustable decaying exponent. In the beginning of the simulation, all sites of $G$ are occupied and there is no link. In the simulation process, two inter-cluster links $l_{i,j}$ and $l_{k,n}$ are randomly chosen and the generalized gravities $g_{ij}$ and $g_{kn}$ are computed. In the strategy $S_{\rm max}$, the link with larger generalized gravity is added. In the strategy $S_{\rm min}$, the link with smaller generalized gravity is added, which include percolation on the Erd\H os-Rényi random graph and the Achlioptas process of explosive percolation as the limiting cases, $d \to \infty$ and $d \to 0$, respectively. Adjustable strategies facilitate or inhibit the network percolation in a generic view. We calculate percolation thresholds $T_c$ and critical exponents $β$ by numerical simulations. We also obtain various finite-size scaling functions for the node fractions in percolating clusters or arrival of saturation length with different intervening strategies.

preprint2017arXiv

Characterizing departure delays of flights in passenger aviation network of United States

Flight delay happens every day in airports all over the world. However, systemic investigation in large scales remains a challenge. We collect primary data of domestic departure records from Bureau of Transportation Statistics of United States, and do empirical statistics with them in form of complementary cumulative distributions functions (CCDFs) and transmission function of the delays. Fourteen main airlines are characterized by two types of CCDFs: shifted power-law and exponentially truncated shifted power-law. By setting up two phenomenological models based on mean-field approximation in temporal regime, we convert effect from other delay factors into a propagation one. Three parameters meaningful in measuring airlines emerge as universal metrics. Moreover, method used here could become a novel approach to revealing practical meanings hidden in temporal big data in wide fields.

preprint2016arXiv

Cascading failures in interdependent networks with finite functional components

We present a cascading failure model of two interdependent networks in which functional nodes belong to components of size greater than or equal to $s$. We find theoretically and via simulation that in complex networks with random dependency links the transition is first-order for $s\geq 3$ and second-order for $s=2$. We find for two square lattices with a distance constraint $r$ in the dependency links that increasing $r$ moves the system from a regime without a phase transition to one with a second-order transition. As $r$ continues to increase the system collapses in a first-order transition. Each regime is associated with a different structure of domain formation of functional nodes.

preprint2016arXiv

Quantifying immediate price impact of trades based on the $k$-shell decomposition of stock trading networks

Traders in a stock market exchange stock shares and form a stock trading network. Trades at different positions of the stock trading network may contain different information. We construct stock trading networks based on the limit order book data and classify traders into $k$ classes using the $k$-shell decomposition method. We investigate the influences of trading behaviors on the price impact by comparing a closed national market (A-shares) with an international market (B-shares), individuals and institutions, partially filled and filled trades, buyer-initiated and seller-initiated trades, and trades at different positions of a trading network. Institutional traders professionally use some trading strategies to reduce the price impact and individuals at the same positions in the trading network have a higher price impact than institutions. We also find that trades in the core have higher price impacts than those in the peripheral shell.

preprint2016arXiv

Recovery of Interdependent Networks

Recent network research has focused on the cascading failures in a system of interdependent networks and the necessary preconditions for system collapse. An important question that has not been addressed is how to repair a failing system before it suffers total breakdown. Here we introduce a recovery strategy of nodes and develop an analytic and numerical framework for studying the concurrent failure and recovery of a system of interdependent networks based on an efficient and practically reasonable strategy. Our strategy consists of repairing a fraction of failed nodes, with probability of recovery $γ$, that are neighbors of the largest connected component of each constituent network. We find that, for a given initial failure of a fraction $1-p$ of nodes, there is a critical probability of recovery above which the cascade is halted and the system fully restores to its initial state and below which the system abruptly collapses. As a consequence we find in the plane $γ-p$ of the phase diagram three distinct phases. A phase in which the system never collapses without being restored, another phase in which the recovery strategy avoids the breakdown, and a phase in which even the repairing process cannot avoid the system collapse.

preprint2015arXiv

Epidemic Model with Isolation in Multilayer Networks

The Susceptible-Infected-Recovered (SIR) model has successfully mimicked the propagation of such airborne diseases as influenza A (H1N1). Although the SIR model has recently been studied in a multilayer networks configuration, in almost all the research the isolation of infected individuals is disregarded. Hence we focus our study in an epidemic model in a two-layer network, and we use an isolation parameter to measure the effect of isolating infected individuals from both layers during an isolation period. We call this process the Susceptible-Infected-Isolated-Recovered ($SI_IR$) model. The isolation reduces the transmission of the disease because the time in which infection can spread is reduced. In this scenario we find that the epidemic threshold increases with the isolation period and the isolation parameter. When the isolation period is maximum there is a threshold for the isolation parameter above which the disease never becomes an epidemic. We also find that epidemic models, like $SIR$ overestimate the theoretical risk of infection. Finally, our model may provide a foundation for future research to study the temporal evolution of the disease calibrating our model with real data.

preprint2015arXiv

Estimating Tipping Points in Feedback-Driven Financial Networks

Much research has been conducted arguing that tipping points at which complex systems experience phase transitions are difficult to identify. To test the existence of tipping points in financial markets, based on the alternating offer strategic model we propose a network of bargaining agents who mutually either cooperate or where the feedback mechanism between trading and price dynamics is driven by an external "hidden" variable R that quantifies the degree of market overpricing. Due to the feedback mechanism, R fluctuates and oscillates over time, and thus periods when the market is underpriced and overpriced occur repeatedly. As the market becomes overpriced, bubbles are created that ultimately burst in a market crash. The probability that the index will drop in the next year exhibits a strong hysteresis behavior from which we calculate the tipping point. The probability distribution function of R has a bimodal shape characteristic of small systems near the tipping point. By examining the S&P500 index we illustrate the applicability of the model and demonstate that the financial data exhibits a hysteresis and a tipping point that agree with the model predictions. We report a cointegration between the returns of the S&P 500 index and its intrinsic value.

preprint2015arXiv

Predicting the extinction of Ebola spreading in Liberia due to mitigation strategies

The Ebola virus is spreading throughout West Africa and is causing thousands of deaths. In order to quantify the effectiveness of different strategies for controlling the spread, we develop a mathematical model in which the propagation of the Ebola virus through Liberia is caused by travel between counties. For the initial months in which the Ebola virus spreads, we find that the arrival times of the disease into the counties predicted by our model are compatible with World Health Organization data, but we also find that reducing mobility is insufficient to contain the epidemic because it delays the arrival of Ebola virus in each county by only a few weeks. We study the effect of a strategy in which safe burials are increased and effective hospitalisation instituted under two scenarios: (i) one implemented in mid-July 2014 and (ii) one in mid-August---which was the actual time that strong interventions began in Liberia. We find that if scenario (i) had been pursued the lifetime of the epidemic would have been three months shorter and the total number of infected individuals 80\% less than in scenario (ii). Our projection under scenario (ii) is that the spreading will stop by mid-spring 2015.

preprint2014arXiv

Network Risk and Forecasting Power in Phase-Flipping Dynamical Networks

In order to model volatile real-world network behavior, we analyze phase-flipping dynamical scale-free network in which nodes and links fail and recover. We investigate how stochasticity in a parameter governing the recovery process affects phase-flipping dynamics, and find the probability that no more than q% of nodes and links fail. We derive higher moments of the fractions of active nodes and active links, $f_n(t)$ and $f_{\ell}(t)$, and define two estimators to quantify the level of risk in a network. We find hysteresis in the correlations of $f_n(t)$ due to failures at the node level, and derive conditional probabilities for phase-flipping in networks. We apply our model to economic and traffic networks.

preprint2014arXiv

Predicting Lifetime of Dynamical Networks Experiencing Persistent Random Attacks

Empirical estimation of critical points at which complex systems abruptly flip from one state to another is among the remaining challenges in network science. However, due to the stochastic nature of critical transitions it is widely believed that critical points are difficult to estimate, and it is even more difficult, if not impossible, to predict the time such transitions occur [1-4]. We analyze a class of decaying dynamical networks experiencing persistent attacks in which the magnitude of the attack is quantified by the probability of an internal failure, and there is some chance that an internal failure will be permanent. When the fraction of active neighbors declines to a critical threshold, cascading failures trigger a network breakdown. For this class of network we find both numerically and analytically that the time to the network breakdown, equivalent to the network lifetime, is inversely dependent upon the magnitude of the attack and logarithmically dependent on the threshold. We analyze how permanent attacks affect dynamical network robustness and use the network lifetime as a measure of dynamical network robustness offering new methodological insight into system dynamics.

preprint2013arXiv

The Boson peak in supercooled water

We perform extensive molecular dynamics simulations of the TIP4P/2005 model of water to investigate the origin of the Boson peak reported in experiments on supercooled water in nanoconfined pores, and in hydration water around proteins. We find that the onset of the Boson peak in supercooled bulk water coincides with the crossover to a predominantly low-density-like liquid below the Widom line $T_W$. The frequency and onset temperature of the Boson peak in our simulations of bulk water agree well with the results from experiments on nanoconfined water. Our results suggest that the Boson peak in water is not an exclusive effect of confinement. We further find that, similar to other glass-forming liquids, the vibrational modes corresponding to the Boson peak are spatially extended and are related to transverse phonons found in the parent crystal, here ice Ih.

preprint2013arXiv

Triple Point in Correlated Interdependent Networks

Many real-world networks depend on other networks, often in non-trivial ways, to maintain their functionality. These interdependent "networks of networks" are often extremely fragile. When a fraction $1-p$ of nodes in one network randomly fails, the damage propagates to nodes in networks that are interdependent and a dynamic failure cascade occurs that affects the entire system. We present dynamic equations for two interdependent networks that allow us to reproduce the failure cascade for an arbitrary pattern of interdependency. We study the "rich club" effect found in many real interdependent network systems in which the high-degree nodes are extremely interdependent, correlating a fraction $α$ of the higher degree nodes on each network. We find a rich phase diagram in the plane $p-α$, with a triple point reminiscent of the triple point of liquids that separates a non-functional phase from two functional phases.

preprint2012arXiv

Critical Points of Correlated Percolation in a Gravitational Link-adding Network Model

Motivated by the importance of geometric information in real systems, a new model for long-range correlated percolation in link-adding networks is proposed with the connecting probability decaying with a power-law of the distance on the two-dimensional(2D) plane. By overlapping it with Achlioptas process, it serves as a gravity model which can be tuned to facilitate or inhibit the network percolation in a generic view, cover a broad range of thresholds. Moreover, it yields a set of new scaling relations. In the present work, we develop an approach to determine critical points for them by simulating the temporal evolutions of type-I, type-II and type-III links(chosen from both inter-cluster links, an intra-cluster link compared with an inter-cluster one, and both intra-cluster ones, respectively) and corresponding average lengths. Numerical results have revealed objective competition between fractions, average lengths of three types of links, verified the balance happened at critical points. The variation of decay exponents $a$ or transmission radius $R$ always shifts the temporal pace of the evolution, while the steady average lengths and the fractions of links always keep unchanged just as the values in Achlioptas process. Strategy with maximum gravity can keep steady average length, while that with minimum one can surpass it. Without the confinement of transmission range, $\bar{l} \to \infty$ in thermodynamic limit, while $\bar{l}$ does not when with it. However, both mechanisms support critical points. In two-dimensional free space, the relevance of correlated percolation in link-adding process is verified by validation of new scaling relations with various exponent $a$, which violates the scaling law of Weinrib's.

preprint2012arXiv

Epidemics on Interconnected Networks

Populations are seldom completely isolated from their environment. Individuals in a particular geographic or social region may be considered a distinct network due to strong local ties, but will also interact with individuals in other networks. We study the susceptible-infected-recovered (SIR) process on interconnected network systems, and find two distinct regimes. In strongly-coupled network systems, epidemics occur simultaneously across the entire system at a critical infection strength $β_c$, below which the disease does not spread. In contrast, in weakly-coupled network systems, a mixed phase exists below $β_c$ of the coupled network system, where an epidemic occurs in one network but does not spread to the coupled network. We derive an expression for the network and disease parameters that allow this mixed phase and verify it numerically. Public health implications of communities comprising these two classes of network systems are also mentioned.

preprint2012arXiv

Fragile to strong crossover coupled to liquid-liquid transition in hydrophobic solutions

Using discrete molecular dynamics simulations we study the relation between the thermodynamic and diffusive behaviors of a primitive model of aqueous solutions of hydrophobic solutes consisting of hard spheres in the Jagla particles solvent, close to the liquid-liquid critical point of the solvent. We find that the fragile-to-strong dynamic transition in the diffusive behavior is always coupled to the low-density/high-density liquid transition. Above the liquid-liquid critical pressure, the diffusivity crossover occurs at the Widom line, the line along which the thermodynamic response functions show maxima. Below the liquid-liquid critical pressure, the diffusivity crossover occurs when the limit of mechanical stability lines are crossed, as indicated by the hysteresis observed when going from high to low temperature and vice versa. These findings show that the strong connection between dynamics and thermodynamics found in bulk water persists in hydrophobic solutions for concentrations from low to moderate, indicating that experiments measuring the relaxation time in aqueous solutions represent a viable route for solving the open questions in the field of supercooled water.

preprint2012arXiv

Identifying States of a Financial Market

The understanding of complex systems has become a central issue because complex systems exist in a wide range of scientific disciplines. Time series are typical experimental results we have about complex systems. In the analysis of such time series, stationary situations have been extensively studied and correlations have been found to be a very powerful tool. Yet most natural processes are non-stationary. In particular, in times of crisis, accident or trouble, stationarity is lost. As examples we may think of financial markets, biological systems, reactors or the weather. In non-stationary situations analysis becomes very difficult and noise is a severe problem. Following a natural urge to search for order in the system, we endeavor to define states through which systems pass and in which they remain for short times. Success in this respect would allow to get a better understanding of the system and might even lead to methods for controlling the system in more efficient ways. We here concentrate on financial markets because of the easy access we have to good data and because of the strong non-stationary effects recently seen. We analyze the S&P 500 stocks in the 19-year period 1992-2010. Here, we propose such an above mentioned definition of state for a financial market and use it to identify points of drastic change in the correlation structure. These points are mapped to occurrences of financial crises. We find that a wide variety of characteristic correlation structure patterns exist in the observation time window, and that these characteristic correlation structure patterns can be classified into several typical "market states". Using this classification we recognize transitions between different market states. A similarity measure we develop thus affords means of understanding changes in states and of recognizing developments not previously seen.

preprint2011arXiv

Robustness of a Tree-like Network of Interdependent Networks

In reality, many real-world networks interact with and depend on other networks. We develop an analytical framework for studying interacting networks and present an exact percolation law for a network of $n$ interdependent networks (NON). We present a general framework to study the dynamics of the cascading failures process at each step caused by an initial failure occurring in the NON system. We study and compare both $n$ coupled Erdős-Rényi (ER) graphs and $n$ coupled random regular (RR) graphs. We found recently [Gao et. al. arXive:1010.5829] that for an NON composed of $n$ ER networks each of average degree $k$, the giant component, $P_{\infty}$, is given by $P_{\infty}=p[1-\exp(-kP_{\infty})]^n$ where $1-p$ is the initial fraction of removed nodes. Our general result coincides for $n=1$ with the known Erdős-Rényi second-order phase transition at a threshold, $p=p_c$, for a single network. For $n=2$ the general result for $P_{\infty}$ corresponds to the $n=2$ result [Buldyrev et. al., Nature, 464, (2010)]. Similar to the ER NON, for $n=1$ the percolation transition at $p_c$, is of second order while for any $n>1$ it is of first order. The first order percolation transition in both ER and RR (for $n>1$) is accompanied by cascading failures between the networks due to their interdependencies. However, we find that the robustness of $n$ coupled RR networks of degree $k$ is dramatically higher compared to the $n$ coupled ER networks of average degree $k$. While for ER NON there exists a critical minimum average degree $k=k_{\min}$, that increases with $n$, below which the system collapses, there is no such analogous $k_{\min}$ for RR NON system.

preprint2010arXiv

Catastrophic Cascade of Failures in Interdependent Networks

Modern network-like systems are usually coupled in such a way that failures in one network can affect the entire system. In infrastructures, biology, sociology, and economy, systems are interconnected and events taking place in one system can propagate to any other coupled system. Recent studies on such coupled systems show that the coupling increases their vulnerability to random failure. Properties for interdependent networks differ significantly from those of single-network systems. In this article, these results are reviewed and the main properties discussed.

preprint2010arXiv

Quarantine generated phase transition in epidemic spreading

We study the critical effect of quarantine on the propagation of epidemics on an adaptive network of social contacts. For this purpose, we analyze the susceptible-infected-recovered (SIR) model in the presence of quarantine, where susceptible individuals protect themselves by disconnecting their links to infected neighbors with probability w, and reconnecting them to other susceptible individuals chosen at random. Starting from a single infected individual, we show by an analytical approach and simulations that there is a phase transition at a critical rewiring (quarantine) threshold w_c separating a phase (w<w_c) where the disease reaches a large fraction of the population, from a phase (w >= w_c) where the disease does not spread out. We find that in our model the topology of the network strongly affects the size of the propagation, and that w_c increases with the mean degree and heterogeneity of the network. We also find that w_c is reduced if we perform a preferential rewiring, in which the rewiring probability is proportional to the degree of infected nodes.

preprint2009arXiv

Designing optimal transport networks

We investigate the optimal design of networks for a general transport system. Our network is built from a regular two-dimensional ($d=2$) square lattice to be improved by adding long-range connections (shortcuts) with probability $P_{ij} \sim r_{ij}^{-α}$, where $r_{ij}$ is the Euclidean distance between sites $i$ and $j$, and $α$ is a variable exponent. We introduce a cost constraint on the total length of the additional links and find optimal transport in the system for $α=d+1$. Remarkably, this condition remains optimal, regardless of the strategy used for navigation, being based on local or global knowledge of the network structure, in sharp contrast with the results obtained for unconstrained navigation using global or local information, where the optimal conditions are $α=0$ and $α=d$, respectively. The validity of our theoretical results is supported by data on the US airport network, for which $α\approx 3.0$ was recently found [Bianconi {\it et al.}, arXiv:0810.4412 (2008)].

preprint2006arXiv

Optimal Path and Minimal Spanning Trees in Random Weighted Networks

We review results on the scaling of the optimal path length in random networks with weighted links or nodes. In strong disorder we find that the length of the optimal path increases dramatically compared to the known small world result for the minimum distance. For Erdős-Rényi (ER) and scale free networks (SF), with parameter $λ$ ($λ>3$), we find that the small-world nature is destroyed. We also find numerically that for weak disorder the length of the optimal path scales logaritmically with the size of the networks studied. We also review the transition between the strong and weak disorder regimes in the scaling properties of the length of the optimal path for ER and SF networks and for a general distribution of weights, and suggest that for any distribution of weigths, the distribution of optimal path lengths has a universal form which is controlled by the scaling parameter $Z=\ell_{\infty}/A$ where $A$ plays the role of the disorder strength, and $\ell_{\infty}$ is the length of the optimal path in strong disorder. The relation for $A$ is derived analytically and supported by numerical simulations. We then study the minimum spanning tree (MST) and show that it is composed of percolation clusters, which we regard as "super-nodes", connected by a scale-free tree. We furthermore show that the MST can be partitioned into two distinct components. One component the {\it superhighways}, for which the nodes with high centrality dominate, corresponds to the largest cluster at the percolation threshold which is a subset of the MST. In the other component, {\it roads}, low centrality nodes dominate. We demonstrate the significance identifying the superhighways by showing that one can improve significantly the global transport by improving a very small fraction of the network.

preprint2005arXiv

Relation Between the Widom line and the Strong-Fragile Dynamic Crossover in Systems with a Liquid-Liquid Phase Transition

We investigate, for two water models displaying a liquid-liquid critical point, the relation between changes in dynamic and thermodynamic anomalies arising from the presence of the liquid-liquid critical point. We find a correlation between the dynamic fragility transition and the locus of specific heat maxima $C_P^{\rm max}$ (``Widom line'') emanating from the critical point. Our findings are consistent with a possible relation between the previously hypothesized liquid-liquid phase transition and the transition in the dynamics recently observed in neutron scattering experiments on confined water. More generally, we argue that this connection between $C_P^{\rm max}$ and dynamic crossover is not limited to the case of water, a hydrogen bond network forming liquid, but is a more general feature of crossing the Widom line. Specifically, we also study the Jagla potential, a spherically-symmetric two-scale potential known to possess a liquid-liquid critical point, in which the competition between two liquid structures is generated by repulsive and attractive ramp interactions.

preprint2004arXiv

Optimization of Network Robustness to Waves of Targeted and Random Attack

We study the robustness of complex networks to multiple waves of simultaneous (i) targeted attacks in which the highest degree nodes are removed and (ii) random attacks (or failures) in which fractions $p_t$ and $p_r$ respectively of the nodes are removed until the network collapses. We find that the network design which optimizes network robustness has a bimodal degree distribution, with a fraction $r$ of the nodes having degree $k_2= (\kav - 1 +r)/r$ and the remainder of the nodes having degree $k_1=1$, where $\kav$ is the average degree of all the nodes. We find that the optimal value of $r$ is of the order of $p_t/p_r$ for $p_t/p_r\ll 1$.

preprint2003arXiv

Time correlations and 1/f behavior in backscattering radar reflectivity measurements from cirrus cloud ice fluctuations

The state of the atmosphere is governed by the classical laws of fluid motion and exhibits correlations in various spatial and temporal scales. These correlations are crucial to understand the short and long term trends in climate. Cirrus clouds are important ingredients of the atmospheric boundary layer. To improve future parameterization of cirrus clouds in climate models, it is important to understand the cloud properties and how they change within the cloud. We study correlations in the fluctuations of radar signals obtained at isodepths of winter and fall cirrus clouds. In particular we focus on three quantities: (i) the backscattering cross-section, (ii) the Doppler velocity and (iii) the Doppler spectral width. They correspond to the physical coefficients used in Navier Stokes equations to describe flows, i.e. bulk modulus, viscosity, and thermal conductivity. In all cases we find that power-law time correlations exist with a crossover between regimes at about 3 to 5 min. We also find that different type of correlations, including 1/f behavior, characterize the top and the bottom layers and the bulk of the clouds. The underlying mechanisms for such correlations are suggested to originate in ice nucleation and crystal growth processes.

preprint1999arXiv

Applications of the Stell-Hemmer Potential to Understanding Second Critical Points in Real Systems

We consider the novel properties of the Stell-Hemmer core-softened potentials. First we explore how the theoretically predicted second critical point for these potentials is related to the occurrence of the experimentally observed solid-solid isostructural critical point. We then discuss how this class of potentials can generate anomalies analogous to those found experimentally in liquid water.

preprint1998arXiv

Inertial Effects on Fluid Flow through Disordered Porous Media

We study the fluid flow through disordered porous media by numerically solving the complete set of the Navier-Stokes equations in a two dimensional lattice with a spatially random distribution of solid obstacles (plaquettes). We simulate viscous and non-viscous flow through these idealized pore spaces to determine the origin of the deviations from the classical Darcy's law behavior. Due to the non-linear contribution of inertia to the transport of momentum at the pore scale, we observe a typical departure from Darcy's law at sufficiently high Reynolds numbers. Moreover, we show that the classical Forchheimer equation provides a valid phenomenological model to correlate the variations of the friction factor of the porous media over a wide range of Reynolds conditions.

preprint1997arXiv

Scaling behavior in economics: I. Empirical results for company growth

We address the question of the growth of firm size. To this end, we analyze the Compustat data base comprising all publicly-traded United States manufacturing firms within the years 1974-1993. We find that the distribution of firm sizes remains stable for the 20 years we study, i.e., the mean value and standard deviation remain approximately constant. We study the distribution of sizes of the ``new'' companies in each year and find it to be well approximated by a log-normal. We find (i) the distribution of the logarithm of the growth rates, for a fixed growth period of one year, and for companies with approximately the same size $S$ displays an exponential form, and (ii) the fluctuations in the growth rates -- measured by the width of this distribution $σ_1$ -- scale as a power law with $S$, $σ_1\sim S^{-β}$. We find that the exponent $β$ takes the same value, within the error bars, for several measures of the size of a company. In particular, we obtain: $β=0.20\pm0.03$ for sales, $β=0.18\pm0.03$ for number of employees, $β=0.18\pm0.03$ for assets, $β=0.18\pm0.03$ for cost of goods sold, and $β=0.20\pm0.03$ for property, plant, & equipment.

preprint1995arXiv

Scaling in Late Stage Spinodal Decomposition with Quenched Disorder

We study the late stages of spinodal decomposition in a Ginzburg-Landau mean field model with quenched disorder. Random spatial dependence in the coupling constants is introduced to model the quenched disorder. The effect of the disorder on the scaling of the structure factor and on the domain growth is investigated in both the zero temperature limit and at finite temperature. In particular, we find that at zero temperature the domain size, $R(t)$, scales with the amplitude, $A$, of the quenched disorder as $R(t) = A^{-β} f(t/A^{-γ})$ with $β\simeq 1.0$ and $γ\simeq 3.0$ in two dimensions. We show that $β/γ= α$, where $α$ is the Lifshitz-Slyosov exponent. At finite temperature, this simple scaling is not observed and we suggest that the scaling also depends on temperature and $A$. We discuss these results in the context of Monte Carlo and cell dynamical models for phase separation in systems with quenched disorder, and propose that in a Monte Carlo simulation the concentration of impurities, $c$, is related to $A$ by $A \sim c^{1/d}$.