Force networks and elasticity in granular silos

We have made experimental observations of the force networks within a two-dimensional granular silo similar to the classical system of Janssen. Models like that of Janssen predict that pressure within a silo saturates with depth as the result of vertical forces being redirected to the walls of the silo where they can then be carried by friction. We use photoelastic particles to obtain information not available in previous silo experiments --the internal force structure. We directly compare various predictions with the results obtained by averaging ensembles of experimentally obtained force networks. We identify several differences between the mean behavior in our system and that predicted by Janssen-like models: We find that the redirection parameter describing how the force network transfers vertical forces to the walls varies with depth. We find that changes in the preparation of the material can cause the pressure within the silo to either saturate or to continue building with depth. Most strikingly, we observe a nonlinear response to overloads applied to the top of the material in the silo. For larger overloads we observe the previously reported “giant overshoot” effect where overload pressure decays only after an initial increase (G. Ovarlez et al., Phys. Rev. E 67, 060302(R) (2003)). For smaller overloads we find that additional pressure propagates to great depth. Analysis of the differences between the inter-grain contact and force networks suggests that, for our system, when the load and the particle weight are comparable, particle elasticity acts to stabilize the force network, allowing deep propagation. For larger loads, the force network rearranges, resulting in the expected, Janssen-like behavior. Thus, a meso-scale network phenomenon results in an observable nonlinearity in the mean pressure profile.     

Violation of finite-size scaling in three dimensions

We reexamine the range of validity of finite-size scaling in the lattice model and the field theory below four dimensions. We show that general renormalization-group arguments based on the renormalizability of the theory do not rule out the possibility of a violation of finite-size scaling due to a finite lattice constant and a finite cutoff. For a confined geometry of linear size L with periodic boundary conditions we analyze the approach towards bulk critical behavior as at fixed for where is the bulk correlation length. We show that for this analysis ordinary renormalized perturbation theory is sufficient. On the basis of one-loop results and of exact results in the spherical limit we find that finite-size scaling is violated for both the lattice model and the field theory in the region . The non-scaling effects in the field theory and in the lattice model differ significantly from each other. Received 5 February 1999     

Emergent organization of oscillator clusters in coupled self-regulatory chaotic maps

Here we introduce a model of parametrically coupled chaotic maps on a one-dimensional lattice. In this model, each element has its internal self-regulatory dynamics, whereby at fixed intervals of time the nonlinearity parameter at each site is adjusted by feedback from its past evolution. Additionally, the maps are coupled sequentially and unidirectionally, to their nearest neighbor, through the difference of their parametric variations. Interestingly we find that this model asymptotically yields clusters of superstable oscillators with different periods. We observe that the sizes of these oscillator clusters have a power-law distribution. Moreover, we find that the transient dynamics gives rise to a 1/f power spectrum. All these characteristics indicate self-organization and emergent scaling behavior in this system. We also interpret the power-law characteristics of the proposed system from an ecological point of view.      

Modeling Evolution of Weighted Clique Networks

We propose a weighted clique network evolution model, which expands continuously by the addition of a new clique (maximal complete sub-graph) at each time step. And the cliques in the network overlap with each other. The structural expansion of the weighted clique network is combined with the edges' weight and vertices' strengths dynamical evolution. The model is based on a weight-driven dynamics and a weights' enhancement mechanism combining with the network growth. We study the network properties, which include the distribution of vertices' strength and the distribution of edges' weight, and find that both the distributions follow the scale-free distribution. At the same time, we also find that the relationship between strength and degree of a vertex are linear correlation during the growth of the network. On the basis of mean-field theory, we study the weighted network model and prove that both vertices' strength and edges' weight of this model follow the scale-free distribution. And we exploit an algorithm to forecast the network dynamics, which can be used to reckon the distributions and the corresponding scaling exponents. Furthermore, we observe that mean-field based theoretic results are
consistent with the statistical data of the model, which denotes the theoretical result in this paper is effective.     

Universal critical properties of the Eulerian bond-cubic model

We investigate the Eulerian bond-cubic model on the square lattice by means of Monte Carlo simulations,using an efficient cluster algorithm and a finite-size scaling analysis.The critical points and four critical exponents of the model are determined for several values of n.Two of the exponents are fractal dimensions,which are obtained numerically for the first time.Our results are consistent with the Coulomb gas predictions for the critical O(n) branch for n < 2 and the results obtained by previous transfer matrix calculations.For n=2,we find that the thermal exponent,the magnetic exponent and the fractal dimension of the largest critical Eulerian bond component are different from those of the critical O(2) loop model.These results confirm that the cubic anisotropy is marginal at n=2 but irrelevant for n < 2.     

粮仓内颗粒压力的测量:Janssen行为及其偏差

基于侧壁分为上下两部分的装置,测量了圆柱对称颗粒系统(竖直粮仓)的压力随高度的变化,并讨论了测得的仓内压力与Janssen理论模型的偏离情况.实验结果显示无论是否对粮仓进行Vanel等建议的让"侧壁摩擦力充分动员"的底座缓慢下沉操作处理,测得的仓内压力都有不同程度的偏离Janssen现象.另外如果没有沉降处理,在填充颗粒的荷载下上侧壁会发生微弱下沉,从而显著增大压力偏离的程度,但沉降处理可将其消除.这些结果表明,粮仓的颗粒重量在仓底和沿侧壁边界上的分配情况受制备方式和装置变形的影响明显,具体分布情况复杂多样.因此颗粒物质的静应力弹性理论分析应该取应力边界条件.     

On the scaling behavior of the chiral phase transition in QCD in finite and infinite volume

We study the scaling behavior of the two-flavor chiral phase transition using an effective quark–meson model. We investigate the transition between infinite-volume and finite-volume scaling behavior when the system is placed in a finite box. We can estimate effects that the finite volume and the explicit symmetry breaking by the current quark masses have on the scaling behavior which is observed in full QCD lattice simulations. The model allows us to explore large quark masses as well as the chiral limit in a wide range of volumes, and extract information about the scaling regimes. In particular, we find large scaling deviations for physical pion masses and significant finite-volume effects for pion masses that are used in current lattice simulations.     

The phase diagram of the Flory-Huggins-de Gennes model of a binary polymer blend

We undertake a numerical study of the Flory-Huggins-de Gennes functional ind=3 dimensions describing a polymer blend. By discretising the functional on a three-dimensional lattice and employing the hybrid Monte Carlo simulation algorithm, we investigate to what extent the inclusion of the term describing fluctuations in local polymer concentration alters the phase diagram of the model. We find that, despite the relatively small weight of the fluctuation term, the coexistence curve is shifted by an appreciable amount from that predicted by naive mean-field theory, which ignores such spatial fluctuations. The direction of the shift is consistent with that already observed in experiment and in simulations of microscopic models of polymer blends. A finite-size scaling analysis indicates that the critical behavior of the model seems to belong to the 3D Ising universality class rather than being mean-field in nature.It is a pleasure to dedicate this paper to Oliver Penrose on the occasion of his 65th birthday.     

Spin twists, domain walls, and the cluster spin-glass phase of weakly doped cuprates

We examine the role of spin twists in the formation of domain walls, often called stripes, by focusing on the spin textures found in the cluster spin glass phases of and . To this end, we derive improved analytic expressions for the spin distortions produced by a frustrating bond, both near the core region of the bond and in the far field, and then derive an improved expression for interaction energies between such bonds. We critique our analytical theory by comparison to numerical solutions of this problem and find excellent agreement. By looking at collections of small numbers of such bonds localized in some region of a lattice, we demonstrate the stability of small “clusters” of spins, each cluster having its own orientation of its antiferromagnetic order parameter. Then, we display a domain wall corresponding to spin twists between clusters of locally ordered spins showing how spin twists can serve as a mechanism for stripe formation. Since the charges are localized in this model, we emphasize that these domain walls are produced in a situation for which no kinetic energy is present in the problem. Received 10 January 2000     

Magnetic ordering and phase transitions in planar antiferromagnetic systems with a Kagomé lattice

We study the process of magnetic ordering in planar antiferromagnetic systems with a Kagomé lattice. It is found that if the interaction between next-nearest-neighbor spins is taken into account, the heat capacity of such systems has a singularity at a finite temperature T. On the basis of a scaling analysis of finite-size systems we study the behavior of thermodynamic quantities in the neighborhood of a phase transition. We find that the phase transition at the critical point is due to discrete-and continuous-symmetry breaking, in which the long-range chiral order and the power-law translational spin order emerge simultaneously. Finally, we calculate the temperatures of the transition to different (with three and nine spins per unit cell) ordered states. Zh. éksp. Teor. Fiz. 113, 2209–2220 (June 1998)     

Supersolid hard-core bosons on the triangular lattice

We determine the phase diagram of hard-core bosons on a triangular lattice with nearest-neighbor repulsion, paying special attention to the stability of the supersolid phase. Similar to the same model on a square lattice we find that for densities rho<1/3 or rho>2/3 a supersolid phase is unstable and the transition between a commensurate solid and the superfluid is of first order. At intermediate fillings 1/3相似文献   

Effect of the ratios of diameter of silo to bead on the pressure screening in granular columns

We present the apparent mass measurements at the bottom of granular packings for different bead and silo sizes. The redirection parameter K in Janssen theory is found to increase with the ratios of the diameter of the silo to the bead. We attribute this feature to the friction between the beads and the confining wall of silo; it is the role of friction that leads to variations in the shielding of vertical stresses as well as pressure screening.     

Damage-spreading dynamic scaling for the ising model on the sierpinski gasket fractal

We study the relaxation towards equilibrium of the ferromagnetic Ising model on the Sierpinski gasket, which is a fractal lattice. We do this by performing Monte Carlo simulations, based on the heat-bath dynamics, and investigating the time evolution of the Hamming distance between two different configurations of the model. Starting with an initial damage created in all lattice sites, we calculate the average values of two quantities that characterize the relaxation process: the nonlinear damage relaxation time (tau), and the time for all sites to be undamaged at least once (tau(c)). We find that tau diverges, at low temperatures, with a dynamical exponent z which depends linearly on the inverse of temperature, as predicted by a generalized scaling theory developed by Henley. There is a complete breakdown of scaling for tau(c).     

二维颗粒堆积中压力问题的格点系统模型

为了便于从理论上探究粮仓效应产生的机理,处理筒仓中颗粒介质的应力分布等问题,将二维颗粒堆积简化为格点系统,并且以随机堆积为物理背景提出了一个由吸收系数p及侧向传递系数q决定的力传递模型.给出了矩阵形式的力传递方程,提出基于二阶差分方程的方法同时求解传递系数矩阵的特征值和特征向量,从理论上导出了一种典型情况下容器底部压力分布与顶部压力分布的关系式.对有效质量随总质量变化关系的理论分析表明,该模型可以给出与Janssen模型类似的结果.对无负载情况下的底部应力分布进行了理论计算,结果表明容器底部中央应力最大,离中央越远应力越小.运用数值计算讨论了p与q对容器底部压力随堆积高度变化曲线的影响.     

Two-threshold model for scaling laws of noninteracting snow avalanches

The sizes of snow slab failure that trigger snow avalanches are power-law distributed. Such a power-law probability distribution function has also been proposed to characterize different landslide types. In order to understand this scaling for gravity-driven systems, we introduce a two-threshold 2D cellular automaton, in which failure occurs irreversibly. Taking snow slab avalanches as a model system, we find that the sizes of the largest avalanches just preceding the lattice system breakdown are power-law distributed. By tuning the maximum value of the ratio of the two failure thresholds our model reproduces the range of power-law exponents observed for land, rock, or snow avalanches. We suggest this control parameter represents the material cohesion anisotropy.     

The Length of an SLE—Monte Carlo Studies

The scaling limits of a variety of critical two-dimensional lattice models are equal to the Schramm–Loewner evolution (SLE) for a suitable value of the parameter κ. These lattice models have a natural parametrization of their random curves given by the length of the curve. This parametrization (with suitable scaling) should provide a natural parametrization for the curves in the scaling limit. We conjecture that this parametrization is also given by a type of fractal variation along the curve, and present Monte Carlo simulations to support this conjecture. Then we show by simulations that if this fractal variation is used to parametrize the SLE, then the parametrized curves have the same distribution as the curves in the scaling limit of the lattice models with their natural parametrization.     

Minimal path on the hierarchical diamond lattice

We consider the minimal paths on a hierarchical diamond lattice, where bonds are assigned a random weight. Depending on the initial distribution of weights, we find all possible asymptotic scaling properties. The different cases found are the small-disorder case, the analog of Lévy's distributions with a power-law decay at-, and finally a limit of large disorder which can be identified as a percolation problem. The asymptotic shape of the stable distributions of weights of the minimal path are obtained, as well as their scaling properties. As a side result, we obtain the asymptotic form of the distribution of effective percolation thresholds for finite-size hierarchical lattices.     

颗粒固体应力转向比的光弹法探测

利用压敏双折射光学特性材料,实验测量了在自然堆积和密堆积两种制样方式下颗粒仓轴向荷载在仓壁上随深度的分布.发现在填充颗粒总质量相同、容器不变的情况下,颗粒仓轴向荷载在仓壁上的分布不随颗粒深度单调变化,而是随深度呈单峰结构,数值模拟与实验观察定性符合,并且发现峰值依赖于荷载大小和样品的制备方式.另外,我们也测量了在不同填充高度下的颗粒底部平均应力随轴向荷载的变化,将其与边壁应力对比,得到颗粒体系不同深度处的转向比.