Collision Statistics of Driven Polydisperse Granular Gases
|
|
- Teresa McLaughlin
- 6 years ago
- Views:
Transcription
1 Commun. Theor. Phys. (Beijing, China) 49 (2008) pp c Chinese Physical Society Vol. 49, No. 5, May 15, 2008 Collision Statistics of Driven Polydisperse Granular Gases CHEN Zhi-Yuan, 1,2, ZHANG Duan-Ming, 1 LI Zhong-Ming, 2 YANG Feng-Xia, 1 and GUO Xin-Ping 1 1 Department of Physics, Huazhong University of Science and Technology, Wuhan , China 2 Department of Physics, Xianning College, Xianning , China (Received June 4, 2007) Abstract We present a dynamical model of two-dimensional polydisperse granular gases with fractal size distribution, in which the disks are subject to inelastic mutual collisions and driven by standard white noise. The inhomogeneity of the disk size distribution can be measured by a fractal dimension d f. By Monte Carlo simulations, we have mainly investigated the effect of the inhomogeneity on the statistical properties of the system in the same inelasticity case. Some novel results are found that the average energy of the system decays exponentially with a tendency to achieve a stable asymptotic value, and the system finally reaches a nonequilibrium steady state after a long evolution time. Furthermore, the inhomogeneity has great influence on the steady-state statistical properties. With the increase of the fractal dimension d f, the distributions of path lengths and free times between collisions deviate more obviously from expected theoretical forms for elastic spheres and have an overpopulation of short distances and time bins. The collision rate increases with d f, but it is independent of time. Meanwhile, the velocity distribution deviates more strongly from the Gaussian one, but does not demonstrate any apparent universal behavior. PACS numbers: Bn, Hv Key words: fractal dimension d f of the disk size distribution, restitution coefficient e, standard white noise, nonequilibrium steady state, statistical properties 1 Introduction Granular materials represent a type of matter not well defined by conventional means. Although each granular particle is obviously solid, an assemblage of these particles shows distinctly nonsolid behavior when subjected to external forces. [1] Their collective behavior is peculiar and also different from other forms of matter, such as liquids, or gases. Although the ordinary statistical mechanical approach successfully deals with large assemblies of microscopic particles, it is not adequate for them. Generally speaking, granular materials cannot be described as equilibrium systems either from the configurational point of view or from the dynamical point of view. It is known that in fact these systems remain easily trapped in some metastable configurations which can last for long time intervals unless they are shaken or perturbed. [1] The understanding of the granular state still represents an open challenge. But, a granular system consisted of fluidizing grains behaves as a nonideal gas, which is relatively easier to be studied. The granular system can be regarded as granular gases. One crucial difference between ordinary gases and granular gases is represented by the intrinsic inelasticity of the interactions among the grains, which makes any theory based on energy conservation, e.g., for ideal gases, not suitable. However, several theoretical methods have been employed to deal with granular gases ranging from hydrodynamic equations, kinetic theories, to molecular dynamics. Usually researchers investigate the statistical properties of granular gases by experiments and computer simulations. They all show that the dynamics of granular gases stems from the inelastic nature of their collisions which leads to a variety of very peculiar phenomena, such as clustering, [2 5] non-maxwellian velocity distributions, [3,6 12] breakdown of energy equipartition in mixtures, [13] and so on, all of which have become one of the most active research topics in nonequilibrium statistical mechanics and fluid dynamics. It is well known that most models describing the statistical properties of granular gases are developed for monodisperse or bidisperse granular materials. In practice, however, actual granular materials are almost always polydisperse and thus the effect of it on the statistical properties should be of interest. Because the experimental study of the statistical properties in a real granular system presents a great challenge, the increasing power of computers awakens an interest in actual granular gases simulations. Recently, in the homogenous cooling state, Lambiotte and Brenig [13] have studied the nonequipartition of energy in granular fluids composed of an arbitrarily large number of components, and have discussed the asymptotic temperature ratio varying with the inelasticity of each species. Zhang et al. [14 17] presented a fractal model of polydisprse granular system. They studied that the fractal dimension d f of the particle size distribution influences the effective thermal conductivity due to the stochastic movement of the particles and the steady-state dynamic behaviors of the driven granular system, such as the granular temperature, kinetic pressure, velocity distribution and distribution of distances between nearest neighbors. While, in order to solve problems in an easier way, they only studied a one-dimensional case. Similar studies of polydisperse granular mixtures in higher dimensional systems are very scarce. The project supported by National Natural Science Foundation of China under Grant No and the Natural Science Foundation of Education Department of Hubei Province of China under Grant Nos. D and kz0627 Corresponding author, czy2004hust@163.com
2 1334 CHEN Zhi-Yuan, ZHANG Duan-Ming, LI Zhong-Ming, YANG Feng-Xia, and GUO Xin-Ping Vol. 49 In this paper, we present a dynamical model of twodimensional polydisperse inelastic granular gases with fractal size distribution. As granular gases are dissipative, it is necessary to feed energy into the system to keep the disks agitated. To thermalize the system, we choose a random external force, which gives frequent kicks to each disk between collisions. Such a driving mechanism has been studied by many authors. [2,3,6,7,12,18] The advantage of this thermalization mechanism, besides its potential physical realizations, lies in the fact that the nonequilibrium stationary state is linearly stable against spatial inhomogeneities. [18] By Monte Carlo simulations, we mainly study the evolution of average energy in transient state and the effect of the inhomogeneity of the disk size distribution on the steady-state statistical properties, focusing on the distributions of path lengths and free times between collisions, collision rate and velocity distribution. 2 Model Let us consider a polydisperse granular mixture of N inelastic hard-disks with radius r i and mass m i (where 1 i N), in which the size distribution of the disks has the fractal characteristic. Fortunately, experimental results show that the particle size distributions of many granular materials exhibit the fractal characteristic, i.e. the size of particles following the fractal distribution. [19 21] Based on the fractal theory, [22] there is a fractal dimension d f as a measurement of the inhomogeneity in the particle size distribution, [20,21] the higher d f implies greater dispersion of the finer particles, which leads to more inhomogeneity in the particle size distribution, and the fractal dimension d f can be obtained, 1 < d f < 2, in two-dimensional case. So, it is reasonable that we consider the polydisperse granular mixture with the fractal characteristic of the disk size distribution, and use the fractal dimension d f to describe the inhomogeneity of the disk size distribution. Therefore, the disk size distribution in the granular system will be more and more inhomogeneous as the fractal dimension d f increases, and a more inhomogeneous disk size distribution indicates a more asymmetric mass distribution of disks. When n 0 /N 1, where n 0 is the number of disks with the maximum size r max, the disk size distribution of the granular system satisfies the size-frequency character by fractal theory [14] ( r ) Y nr (r) = 1 N 1 df n 0, (1) r max where Y nr is the ratio of n r to N, n r is the number of disks whose size is smaller than r, N is the total number of disks, and d f is the fractal dimension of disk size distribution, here 1 < d f < 2. For simplicity, we assume that the surface of disks is smooth and the mass density of disks is identical. Then, the mass of a disk with radius r is m = πr 2 ρ m, (2) where ρ m is the mass density of the disks. From Eqs. (1) and (2), the mass and the radius of any disk in the mixture can be respectively expressed as m = m max [ N n 0 (1 Y nr )] 2/df, (3) r = r max [ N n 0 (1 Y nr )] 1/df. (4) If the values of N, m max, r max, n 0 and d f are given, according to Eqs. (3) and (4) we can randomly evaluate the mass and the radius of every disk in the polydisperse granular mixture. In this work, an inelastic hard disk model is implemented. We simulate an ensemble of N polydisperse inelastic hard-disks in a two-dimensional horizontal box of side L with periodic boundary conditions. The area packing fraction φ of the system may be defined as φ = N / i=1 πr2 i L 2. The disk collisions are inelastic: the total momentum is conserved, while the component of the relative velocity parallel to the direction joining the center of the disks is reduced to a fraction e (0 < e < 1) of its initial value, lowering in this way the kinetic energy of the pair: m j v [ i = v i (1 + e) (vi v j ) ˆn ]ˆn, (5) m i + m j m i v [ j = v j + (1 + e) (vi v j ) ˆn ]ˆn, (6) m i + m j where v i and v j are the velocities of the colliding disks after the collision, ˆn = (x i (t) x j (t))/(r i + r j ) is the unit vector along the line of centers x i (t) and x j (t) of the colliding disks at contact. In the interval between two subsequent collisions, the motion of each disk i is governed by the following Langevin equation: [3,12] m i t v i (t) = m i ηv i (t) + 2m i ηt F ξ i (t), (7) t x i (t) = v i (t), (8) where the function ξ i (t) is a stochastic process with average ξ i (t) = 0 and correlation ξi α(t 1)ξ β j (t 2) = δ ij δ(t 1 t 2 )δ αβ (α and β being component indices), i.e., a standard white noise. This means that each disk feels a hot fluid with a temperature T F and a viscosity coefficient η. Therefore, the disks obtain the kinetic energy from the white noise forcing, execute Brownian motion between collisions, and dissipate the kinetic energy through inelastic collisions and dump. From Eqs. (5) and (6), the energy dissipation due to the inelastic collisions during one collision is E = (1 e2 )[(v i v j ) ˆn] 2 m i m j 2 m i + m j = (1 e2 )[(v i v j ) ˆn] 2 m i (m i + m ij ), (9) 2 2m i + m ij where m ij = m j m i. Obviously, both the mass difference m ij between the two colliding disks and the restitution coefficient e can influence energy dissipation. We know that, a larger value of the fractal dimension d f represents more inhomogeneity of size distribution, which makes the mass difference between any two colliding disks greater. Therefore, when the restitution coefficient e is
3 No. 5 Collision Statistics of Driven Polydisperse Granular Gases 1335 fixed (0 < e < 1), the greater the fractal dimension d f is, the more the energy dissipation E is. Likewise, if the fractal dimension d f remains the same, the smaller the restitution coefficient e is, the more the energy dissipation E is. 3 Simulations and Results The system considered here is a dilute inelastic gas with an area packing fraction φ < 0.4. By Monte Carlo method, we perform the simulations of the disk colliding and Brownian motion between collisions. In the actual simulations, we reduce any variable to dimensionless one. For definiteness of our discussion, we use m max = 20, r max = 1, n 0 = 1, T F = 1, and η = 0.1 in all simulations performed. The simulation scheme consists of a discrete time integration of the motion of the disks. At each time step of length t (where t τ c, and τ c is average collision time), the following operations are performed: [12] (i) Free streaming: equations (7) and (8) are integrated for a time step t disregarding possible interactions among the disks; (ii) Collisions: Every disk has a probability p = t/τ c of undergoing a collision. Its collision mate is chosen among the disks in a circle of fixed radius r B with a probability proportional to its relative velocity. Initially, the disks are randomly placed on the two-dimensional horizontal square box which is divided into L 2 identical square lattices. They are then randomly given an initial velocity which is consistent with their position by a Gaussian random generator. The disks subsequently move, if necessary, to eliminate disks overlap. The simulations proceed from collision to collision using an event-driven algorithm, [23] and the disk motion between two consecutive collisions is governed by Eqs. (7) and (8). After a long evolution time, the system reaches a nonequilibrium steady state which is a result of the balance between the dissipation of the kinetic energy due to inelastic collisions and dump, and the energy injection due to the driving mechanism. 3.1 Approach to Steady State Starting from any initial configuration, the system with dissipation (the restitution coefficient 0 < e < 1) reaches a steady state after a large enough number of collisions. To ensure that the system has reached a steady state, we relax the system until the average energy does not drift. Figure 1 shows how the average granular energy per disk relaxes with time toward a steady state for a system with L = 10, N = , e = 0.6, and d f = 1.5. It is found that the relaxation is exponential as E(t) = E( ) + [E(0) E( )] e t/τ B. (10) E(0) is the initial value of the average energy per disk, E( ) is the asymptotic value of the energy curve when the time approaches infinity, where E(0) = , E( ) = We find that the relaxation time τ B (= 50.1) is independent of the initial state, and is controlled by the model parameters (i.e., the restitution coefficient e, the fractal dimension d f, the total number of disks N etc.). Fig. 1 Average energy relaxation starting from a random initial configuration for a system with L = 10, N = , e = 0.6, and d f = 1.5. The vertical axis is the average kinetic energy per particle (1/2) N m i=1 iυi 2 /N. The best fit using Eq. (11) is also shown with E( ) = , A = , and B = In Fig. 1, the dotted line represents the first-order exponential fit. Therefore, the time-dependent relaxation is also given by E(t) = E( ) + A e Bt, (11) where A is a constant depending on the initial condition, B = 1/τ B = Equation (11) further indicates that the average energy of the system evolves from the initial state into a stable state with an exponentially decaying form. When the system reaches a steady state, we start the simulations with the length of collisions per disk. Thereafter, we can obtain the data with small statistical error. 3.2 Distributions of Path Lengths and Free Times Between Collisions By basic kinetic theory arguments, [24] the distribution of path lengths between collision events for an elastic hardsphere gas (and by a similar treatment the distribution of free times between collisions) is given by P(l) = (2 2φ) e 2 2φl. (12) Apparently, the distribution P(l) should follow a simple exponential form depending only on the packing fraction φ of the assembly. As in Ref. [23], the distribution of path lengths from the geometric distance between collision events is calculated and shown in Fig. 2(a) for all systems with different values of the fractal dimension d f at the same packing fraction φ. For comparison, we plot in Fig. 2(a) the corresponding distribution of an elastic gas which has the same fractal dimension d f and area packing fraction φ of our inelastic mixture. However, it is clear from the dashed line in Fig. 2(a) that the simple form given by Eq. (12) does not describe the behavior over all l in the inelastic case. Therefore, the presence of the inelasticity of collisions modifies
4 1336 CHEN Zhi-Yuan, ZHANG Duan-Ming, LI Zhong-Ming, YANG Feng-Xia, and GUO Xin-Ping Vol. 49 such a simple exponential law in the way from the symbol shown in Fig. 2(a). It is seen from Fig. 2(a) that, in the inelastic case, the distribution of path lengths deviates from the theoretical prediction for elastic hard spheres. In particular, P(l) shows a peak in the small l region, not present in elastic systems, demonstrating the inherent clustering in this system. From Fig. 2(a) it can also be found that, as the fractal dimension d f increases with respect to the same restitution coefficient e (0 < e < 1), the probability of the short-path increases. This indicates that when the value of the fractal dimension d f augments in the same inelastic case, the system becomes more and more clusterized. This phenomenon is consistent with our recent simulations in one-dimensional case. [17] the fractal dimension d f increases at the same restitution coefficient e (0 < e < 1), the probability of the short-time increases. Likewise, for comparison we plot in Fig. 2(a) also the corresponding distribution of an elastic gas. The effect of the shorter average collision time is also clearly visible. Such a behavior is similar to Refs. [8] and [9]. As we know, the higher d f implies greater dispersion of the finer particles, which leads to the particle number density increasing at the same packing fraction φ. Thus, more energy dissipation is caused by the inhomogeneity of the disk size distribution and the higher collision rates due to increased density in the same inelasticity. The effect of energy dissipation leads to strong spatial correlations of density and velocities, which reduce disk relative motion. This smaller value of disk relative velocity implies that the disks tend to cluster and stream when driven. Therefore, the physical reason for the distributions of paths and free times between collisions lies in the existence of strong spatial correlations of density and velocities, which lead to the presence of clusters where the path lengths and the free times are shorter than they are in a uniform system. 3.3 Collision Rate Fig. 2 (a) The probability distributions of path lengths P(l) vs. l, on a log-linear scale. The dashed line shows the theoretical form given by Eq. (12) derived for elastic particles, while the symbols represent the simulations. (b) The probability distributions of free times P(τ) vs. τ, on a log-linear scale. The symbols represent the simulations in the inelastic case, and the dashed line is the corresponding distributions for an elastic system. In both figures φ = 0.2, L = 10, η = 0.1, T F = 1, and e = 0.6. The distribution of free times between collisions P(τ) is also computed and shown in Fig. 2(b). It displays similar behavior to the path length distributions, namely an overpopulation of the short-time bins. One clearly sees from Fig. 2(b) that the probability density for which a disk suffers a collision in a short interval is enhanced with respect to the elastic case. It is also found that, when Fig. 3 (a) Real time vs. total number of collisions for a system with different values of the fractal dimension d f. (b) Collision rate ζ vs. d f. In both figures N = , L = 10, η = 0.1, T F = 1, and e = 0.6. In our simulations, we count the total number of collisions starting from a configuration in the steady state. Figure 3(a) shows that the total collision number (C) is linear with time for all systems with different values of the fractal dimension d f in the same inelasticity case
5 No. 5 Collision Statistics of Driven Polydisperse Granular Gases 1337 (0 < e < 1). The reciprocal slopes of the linear-fit lines are the collision rates for different fractal dimension d f. In the steady state the collision rate is constant, independent of time. This is in contrast to the time-dependent collision rate in the freely evolving monodisperse granular medium where the total number of collisions increases as ln(1 + t/t e ) for homogeneous cooling. [25,26] One would expect that the collision rate should be an increasing function with decreasing d f, as the velocity is larger and the collision is more likely to happen in the smaller d f system. However, we find from Fig. 3(a) that the collision rate increases as the energy dissipation increases. Figure 3(b) shows the dependence of the collision rate ζ as a function of the fractal dimension d f. The increase of collision rate with increasing d f may be linked with the clustering mechanism. [4,5] If the velocity distribution is a Gaussian then γ x = 3, and if the distribution is non-gaussian then γ x exceeds 3. In Fig. 4(b), γ x has been plotted as a function of d f for e = 1.0 and e = 0.6, respectively. In the elastic case, γ x would be 3 over all d f, this demonstrates the velocity distribution is a Gaussian. However, in the inelastic case, the increase in γ x shown in Fig. 4(b) implies that velocity distribution function deviates more strongly from a Gaussian as d f is increased. We find that the calculated values of the kurtosis changing with the fractal dimension d f have the same evolution as the above observation in Fig. 4(a). Thus, we can conclude that the dependence of γ x on d f may be due to the dissipation in the inelasticity. This analysis is similar to MD simulations of Brey and Ruiz Montero [11] and experimental investigations of D.L. Blair and A. Kudrolli. [9,10] 3.4 Velocity Distributions The distribution of velocities is normally displayed in a log-linear fashion to accentuate the tails of the velocity distribution function; however, this suppresses the deviations at low velocities. For comparison, all the disk velocities are rescaled by their mean square values. The distribution of the x-component of the rescaled velocities is plotted in Fig. 4(a) when the system with different values of the fractal dimension d f reaches a stationary state. In the elastic case, the velocity distribution is very well fitted by Gaussian. However, in the inelastic case, we observe that the velocity distribution displays strong deviations from Gaussian one both at low and high velocities. In fact, the kurtosis of the rescaled velocity distribution function rises higher distinctly with the fractal dimension d f increasing at the same restitution coefficient e (0 < e < 1), whereas, the fractal dimension d f is larger, the velocity distribution function has more extended tails. This phenomenon is consistent with our previous simulation results in one-dimensional case. [16,17] From our velocity distribution functions, whose corresponding inelasticity of collisions remains the same (0 < e < 1) and where the fractal dimension d f varies, we cannot find universality of the velocity distribution functions that describes the overall form. Such a lack of universality was also noticed in Refs. [3] and [8]: which is shown that the tails of the velocity distribution function become broader when the dissipation rate is increased. Here the mechanism is similar: the fractal dimension d f is larger, the dissipation rate is higher. The lack of universality was also noted experimentally by Blair and Kudrolli. [9] It is well known that the degree of the deviation from a Gaussian can be measured by the kurtosis of the velocity distribution function. To further show the deviation of the velocity distribution functions, we quantitatively calculate the kurtosis of the velocity distribution functions. The kurtosis is obtained by the following equation γ x = υ4 x υ 2 x 2. (13) Fig. 4 (a) The rescaled velocity distribution P(υ x /σ x ) vs. υ x /σ x (σ x = υ 2 x 1/2 ) in two different regimes (an elastic case with e = 1.0, and an inelastic case with e = 0.6) on a log-linear scale for a system with N = , L = 10, d f = 1.2, 1.5, and 1.8 respectively. In both cases, the symbols represent the simulations. The dotted line represents the Gaussian fit. (b) The kurtosis γ x calculated from P(υ x /σ x ) as a function of d f for e = 1.0 and 0.6, respectively. The dashed line is the values given by a Gaussian in the elastic case, and the solid line represents the values of non-gaussian in the inelastic case. The velocity distribution for the y-component P(υ y /σ y ) versus υ y /σ y for each d f is also calculated. We have found that the distribution of the y-component of the disk velocities resembles the velocity distribution function for the x-component. Therefore, our analysis shows that the present consensus emerging from various studies tends
6 1338 CHEN Zhi-Yuan, ZHANG Duan-Ming, LI Zhong-Ming, YANG Feng-Xia, and GUO Xin-Ping Vol. 49 to reach the conclusion of the absence of universality in the velocity distribution: [27] various experimental conditions and various energy injection modes lead to different distributions. 4 Summary and Conclusion In this paper, we present a dynamical model of polydisperse inelastic hard disks with the fractal size distribution, in which the disks are driven by standard white noise and constrained to move in a two-dimensional horizontal square box with periodic boundary conditions. The inhomogeneity of the disk size distribution is measured by a fractal dimension d f, and the higher d f implies greater dispersion of the finer disks, which leads to more inhomogeneity in the disk size distribution. By Monte Carlo simulations, we perform long-time simulations of the disk colliding and Brownian motion between collisions in the same inelasticity case (0 < e < 1). The simulation results indicate that the average energy of the system decays exponentially, with a tendency to achieve a stable asymptotic value, and the system finally reaches a nonequilibrium steady state after a long evolution time. Moreover, the inhomogeneity of the disk size distribution has great effect on the steady-state statistical properties of the system, focusing on the distributions of path lengths and free times between collisions, collision rate and velocity distribution. With the increase of the fractal dimension d f, the distributions of path lengths and free times between collisions are shown to deviate more obviously from expected theoretical forms for elastic spheres and have an overpopulation of short distances and time bins, and the probability of the short-path and short time bins increases. The collision rate increases with d f, but it is independent of time. As d f is increased, the velocity distribution deviates more strongly from the Gaussian one, but does not demonstrate any apparent universal behavior. Of course, these conclusions are based on the numerical studies of simplified models, so the universality is still to be considered as an open question from a point of experiment. Acknowledgments The authors would like to thank Prof. YU Bo-Ming for helpful discussions. References [1] H.M. Jaeger, S.R. Nagel, and R.P. Behringer, Rev. Mod. Phys. 68 (1996) [2] D.R.M. Williams and F.C. MacKintosh, Phys. Rev. E 54 (1996) R9. [3] A. Puglisi, V. Loreto, U.M.B. Marconi, A. Petri, and A. Vulpiani, Phys. Rev. Lett. 81 (1998) 3848; Phys. Rev. E 59 (1999) [4] M.A. Hopkins and M.Y. Louge, Phys. Fluids A 3 (1991) 47. [5] I. Goldhirsch and G. Zanetti, Phys. Rev. Lett. 70 (1993) [6] T.P.C. Van Noije and M.H. Ernst, Granular Matter 1 (1998) 57. [7] S.J. Moon, M.D. Shattuck, and J.B. Swift, Phys. Rev. E 64 (2001) [8] D. Paolotti, C. Cattuto, U.M.B. Marconi, and A. Puglisi, Granular Matter 5 (2003) 75. [9] D.L. Blair and A. Kudrolli, Phys. Rev. E 67 (2003) [10] D.L. Blair and A. Kudrolli, Phys. Rev. E 64 (2001) [11] J.J. Brey and M.J. Ruiz-Montero, Phys. Rev. E 67 (2003) [12] A. Puglisi, A. Baldassarri, and V. Loreto, Phys. Rev. E 66 (2002) [13] R. Lambiotte and L. Brenig, Phys. Rev. E 72 (2005) [14] D.M. Zhang, Z. Zhang, and B.M. Yu, Commun. Theor. Phys. (Beijing, China) 31 (1999) 373; D.M. Zhang, Y.J. Lei, B.M. Yu, et al., Commun. Theor. Phys. (Beijing, China) 37 (2002) 231; D.M. Zhang, Y.J. Lei, B.M. Yu, et al., Commun. Theor. Phys. (Beijing, China) 40 (2003) 491. [15] D.M. Zhang, Y.J. Lei, G.J. Pan, and B.M. Yu, Chin. Phys. Lett. 20 (2003) 2221; D.M. Zhang, X.Y. Su, B.M. Yu, et al., Commun. Theor. Phys. (Beijing, China) 44 (2005) 551. [16] D.M. Zhang, R. Li, X.Y. Su, G.J. Pan, and B.M. Yu, J. Phys. A: Math. Gen. 38 (2005) [17] D.M. Zhang, Z.Y. Chen, Y.P. Yin, et al., Physica A 374 (2007) 187; Z.Y. Chen, D.M. Zhang, Z.C. Zhong, et al., Commun. Theor. Phys. (Beijing, China) 47 (2007) [18] I. Pagonabarraga, E. Trizac, T.P.C. Van Noije, and M.H. Ernst, Phys. Rev. E 65 (2001) [19] H. Xie, R. Bhasker, and J. Li, Minerals and Metallurgical Processing 2 (1993) 36; M.D. Normand and M. Peleg, Powder Technol. 45 (1986) 271. [20] L.M. Zhang, Coll. Surf. A: Physicochem. Eng. Aspects 202 (2002) 1. [21] J.P. Hyslip and L.E. Vallejo, Eng. Geol. 48 (1997) 231. [22] B.B. Mandelbrot, The Fractal Geometry of Nature, W.H. Freeman, New York (1982). [23] M.P. Allen and D.J. Tildesley, Computer Simulation of Liquids, Oxford University Press, New York (1989). [24] C.E. Hecht, Statistical Thermodynamics and Kinetic Theory, Dover Publications, New York (1998). [25] P. Deltour and J.L. Barrat, J. Phys. I 7 (1997) 137. [26] J.A.G. Orza, R. Brito, T.P.C. Van Noije, and M.H. Ernst, Int. J. Mod. Phys. C 8 (1997) 953. [27] A. Barrat, E. Trizac, and M.H. Ernst, J. Phys.: Condens. Matter 17 (2005) S2429.
Velocity Distributions in a Cooling Granular Gas
School of Physical Sciences, Jawaharlal Nehru University, New Delhi 1167, India E-mail: puri@mail.jnu.ac.in Syed Rashid Ahmad School of Physical Sciences, Jawaharlal Nehru University, New Delhi 1167, India
More informationAnomalous Collective Diffusion in One-Dimensional Driven Granular Media
Typeset with jpsj2.cls Anomalous Collective Diffusion in One-Dimensional Driven Granular Media Yasuaki Kobayashi and Masaki Sano Department of Physics, University of Tokyo, Hongo, Tokyo 113-0033
More informationOn high energy tails in inelastic gases arxiv:cond-mat/ v1 [cond-mat.stat-mech] 5 Oct 2005
On high energy tails in inelastic gases arxiv:cond-mat/0510108v1 [cond-mat.stat-mech] 5 Oct 2005 R. Lambiotte a,b, L. Brenig a J.M. Salazar c a Physique Statistique, Plasmas et Optique Non-linéaire, Université
More informationBOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES
BOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES Vicente Garzó Departamento de Física, Universidad de Extremadura Badajoz, SPAIN Collaborations Antonio Astillero, Universidad de Extremadura José
More informationKinetic Theory for Binary Granular Mixtures at Low Density
10 Kinetic Theory for Binary Granular Mixtures at Low Density V. Garzó Departamento de Física, Universidad de Extremadura, E-06071 Badajoz, Spain vicenteg@unex.es Many features of granular media can be
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationOn the velocity distributions of the one-dimensional inelastic gas
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICSA: MATHEMATICAL AND GENERAL J. Phys. A: Math. Gen. 35 (2002) 463 480 PII: S0305-4470(02)30071-4 On the velocity distributions of the one-dimensional inelastic
More informationSpatial and Temporal Behaviors in a Modified Evolution Model Based on Small World Network
Commun. Theor. Phys. (Beijing, China) 42 (2004) pp. 242 246 c International Academic Publishers Vol. 42, No. 2, August 15, 2004 Spatial and Temporal Behaviors in a Modified Evolution Model Based on Small
More informationCHAPTER V. Brownian motion. V.1 Langevin dynamics
CHAPTER V Brownian motion In this chapter, we study the very general paradigm provided by Brownian motion. Originally, this motion is that a heavy particle, called Brownian particle, immersed in a fluid
More informationKinetic Theory of a Confined Quasi-Two-Dimensional Gas of Hard Spheres. J. Javier Brey *, Vicente Buzón, Maria Isabel García de Soria and Pablo Maynar
entropy Article Kinetic Theory of a Confined Quasi-Two-Dimensional Gas of Hard Spheres J. Javier Brey *, Vicente Buzón, Maria Isabel García de Soria and Pablo Maynar Física Teórica, Universidad de Sevilla,
More informationarxiv: v1 [cond-mat.stat-mech] 26 Jan 2008
Long Time Tail of the Velocity Autocorrelation Function in a Two-Dimensional Moderately Dense Hard Disk Fluid Masaharu Isobe Graduate School of Engineering, Nagoya Institute of Technology, Nagoya 466-8555,
More informationThe Truth about diffusion (in liquids)
The Truth about diffusion (in liquids) Aleksandar Donev Courant Institute, New York University & Eric Vanden-Eijnden, Courant In honor of Berni Julian Alder LLNL, August 20th 2015 A. Donev (CIMS) Diffusion
More informationVelocity statistics in excited granular media
CHAOS VOLUME 9, NUMBER 3 SEPTEMBER 1999 Velocity statistics in excited granular media W. Losert, D. G. W. Cooper, and J. Delour Department of Physics, Haverford College, Haverford, Pennsylvania 19041 A.
More information1714 A. Barrat and E. Trizac free pressure tensor (r P¼0) is obeyed, taking into account anisotropies, boundaries and corrections to the ideal gas equ
MOLECULAR PHYSICS, 10JUNE 2003, VOL. 101, NO. 11, 1713 1719 A molecular dynamics Maxwell Demon experiment for granular mixtures ALAIN BARRAT and EMMANUEL TRIZAC* Laboratoire de Physique The orique (UMR
More informationVelocity Distribution in a Viscous Granular Gas
Velocity Distribution in a Viscous Granular Gas Alexandre Rosas,, Daniel ben-avraham 3, and Katja Lindenberg,4 Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA
More informationBoundary Dissipation in a Driven Hard Disk System
Boundary Dissipation in a Driven Hard Disk System P.L. Garrido () and G. Gallavotti (2) () Institute Carlos I for Theoretical and Computational Physics, and Departamento de lectromagnetismo y Física de
More informationStatistical Properties of a Ring Laser with Injected Signal and Backscattering
Commun. Theor. Phys. (Beijing, China) 35 (2001) pp. 87 92 c International Academic Publishers Vol. 35, No. 1, January 15, 2001 Statistical Properties of a Ring Laser with Injected Signal and Backscattering
More informationAsymptotic solutions of the nonlinear Boltzmann equation for dissipative systems
Contents Part I Kinetic Theory Asymptotic solutions of the nonlinear Boltzmann equation for dissipative systems Matthieu Ernst, Ricardo Brito..................................... 3 The Homogeneous Cooling
More informationForcing independent velocity distributions in an experimental granular fluid
Forcing independent velocity distributions in an experimental granular fluid P. M. Reis,* R. A. Ingale, and M. D. Shattuck Benjamin Levich Institute, The City College of the City University of New York,
More informationRheology of two- and three-dimensional granular mixtures under uniform shear flow: Enskog kinetic theory versus molecular dynamics simulations
Rheology of two- and three-dimensional granular mixtures under uniform shear flow: Enskog kinetic theory versus molecular dynamics simulations José María Montanero Departamento de Electrónica e Ingeniería
More informationMonte Carlo Study of Planar Rotator Model with Weak Dzyaloshinsky Moriya Interaction
Commun. Theor. Phys. (Beijing, China) 46 (2006) pp. 663 667 c International Academic Publishers Vol. 46, No. 4, October 15, 2006 Monte Carlo Study of Planar Rotator Model with Weak Dzyaloshinsky Moriya
More informationEPJ E. Soft Matter and Biological Physics. EPJ.org. Segregation by thermal diffusion in moderately dense granular mixtures
EPJ E Soft Matter and Biological Physics EPJ.org your physics journal Eur. Phys. J. E 9, 61 74 009) DOI: 10.1140/epje/i009-10488-4 Segregation by thermal diffusion in moderately dense granular mixtures
More informationAndrés Santos* Badajoz (Spain)
Granular Poiseuille flow Andrés Santos* University of Extremadura Badajoz (Spain) *In collaboration with Mohamed Tij, Université Moulay Ismaïl, Meknès (Morocco) Departament de Física, Universitat Autònoma
More informationarxiv:cond-mat/ v2 12 Jul 1999
Simple Model with Facilitated Dynamics for Granular Compaction J. Javier Brey, A. Prados Física Teórica, Facultad de Física, Universidad de Sevilla, Apdo. de Correos 1065, E-41080 Sevilla, Spain arxiv:cond-mat/9907144v2
More informationFree cooling of a binary granular fluid of inelastic rough hard spheres Andrés Santos Universidad de Extremadura, Badajoz (Spain) In collaboration with Gilberto M. Kremer and Vicente Garzó Minimal model
More informationEffects of Atomic Coherence and Injected Classical Field on Chaotic Dynamics of Non-degenerate Cascade Two-Photon Lasers
Commun. Theor. Phys. Beijing China) 48 2007) pp. 288 294 c International Academic Publishers Vol. 48 No. 2 August 15 2007 Effects of Atomic Coherence and Injected Classical Field on Chaotic Dynamics of
More informationarxiv: v1 [cond-mat.stat-mech] 6 Mar 2008
CD2dBS-v2 Convergence dynamics of 2-dimensional isotropic and anisotropic Bak-Sneppen models Burhan Bakar and Ugur Tirnakli Department of Physics, Faculty of Science, Ege University, 35100 Izmir, Turkey
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 6 Oct 2002
Competitive Clustering in a Bi-disperse Granular Gas René Mikkelsen, Devaraj van der Meer, Ko van der Weele, and Detlef Lohse Department of Applied Physics and J.M. Burgers Centre for Fluid Dynamics, University
More informationarxiv: v3 [cond-mat.stat-mech] 2 Dec 2015
Inelastic Maxwell models for monodisperse gas-solid flows Aleksander Kubicki Departamento de Física, Universidad de Extremadura, Universidad de Extremadura, E-0607 Badajoz, Spain arxiv:4.345v3 [cond-mat.stat-mech]
More informationarxiv: v1 [cond-mat.dis-nn] 25 Oct 2017
Slow and Long-ranged Dynamical Heterogeneities in Dissipative Fluids arxiv:70.09488v [cond-mat.dis-nn] 25 Oct 207 Introduction Supercooled liquids, colloidal suspensions, and granular systems show evidence
More informationAdiabatic piston - Draft for MEMPHYS
diabatic piston - Draft for MEMPHYS 4th ovember 2006 Christian GRUBER (a) and nnick LESE (b) (a) Institut de héorie des Phénomènes Physiques École Polytechnique Fédérale de Lausanne, CH-05 Lausanne, Switzerland
More informationContinuum Model of Avalanches in Granular Media
Continuum Model of Avalanches in Granular Media David Chen May 13, 2010 Abstract A continuum description of avalanches in granular systems is presented. The model is based on hydrodynamic equations coupled
More informationComputational test of kinetic theory of granular media
Physica A 274 (1999) 158 170 www.elsevier.com/locate/physa Computational test of kinetic theory of granular media M.D. Shattuck a, C. Bizon b, J.B. Swift a, Harry L. Swinney a; a Center for Nonlinear Dynamics
More informationarxiv:cond-mat/ v1 [cond-mat.soft] 13 Nov 2003
Europhysics Letters PREPRINT arxiv:cond-mat/31133v1 cond-mat.soft] 13 Nov 23 Oscillatory instability in a driven granular gas Evgeniy Khain and Baruch Meerson Racah Institute of Physics, Hebrew University
More informationThe First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
Commun. Theor. Phys. (Beijing, China 35 (2 pp. 42 46 c International Academic Publishers Vol. 35, No. 4, April 5, 2 The First Principle Calculation of Green Kubo Formula with the Two-Time Ensemble Technique
More informationSliding Friction in the Frenkel-Kontorova Model
arxiv:cond-mat/9510058v1 12 Oct 1995 to appear in "The Physics of Sliding Friction", Ed. B.N.J. Persson (Kluwer Academic Publishers), 1995 Sliding Friction in the Frenkel-Kontorova Model E. Granato I.N.P.E.,
More informationAn alternative approach to dissipative particle dynamics
EUROPHYSICS LETTERS 15 July 1999 Europhys. Lett., 47 (2), pp. 145-151 (1999) An alternative approach to dissipative particle dynamics C. P. Lowe Computational Physics, Delft University of Technology Lorentzweg
More informationEquivalence of the freely cooling granular gas to the sticky gas
Equivalence of the freely cooling granular gas to the sticky gas Mahendra Shinde, 1, * Dibyendu Das, 1, and R. Rajesh 2, 1 Department of Physics, Indian Institute of Technology, Bombay, Powai, Mumbai-400
More informationComputer Modeling of Binary Dipolar Monolayers
Proceedings of the 8 th International Conference on Applied Informatics Eger, Hungary, January 27 30, 2010. Vol. 1. pp. 329 336. Computer Modeling of Binary Dipolar Monolayers Imre Varga a, Ferenc Kun
More informationComputer generation of dense polydisperse sphere packings
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Computer generation of dense polydisperse sphere packings Anuraag R. Kansal Department of Chemical Engineering, Princeton University, Princeton,
More informationA Brownian ratchet driven by Coulomb friction
XCIX Congresso Nazionale SIF Trieste, 23 September 2013 A Brownian ratchet driven by Coulomb friction Alberto Petri CNR - ISC Roma Tor Vergata and Roma Sapienza Andrea Gnoli Fergal Dalton Giorgio Pontuale
More informationRandom Averaging. Eli Ben-Naim Los Alamos National Laboratory. Paul Krapivsky (Boston University) John Machta (University of Massachusetts)
Random Averaging Eli Ben-Naim Los Alamos National Laboratory Paul Krapivsky (Boston University) John Machta (University of Massachusetts) Talk, papers available from: http://cnls.lanl.gov/~ebn Plan I.
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationLattice Boltzmann Simulation of One Particle Migrating in a Pulsating Flow in Microvessel
Commun. Theor. Phys. 56 (2011) 756 760 Vol. 56, No. 4, October 15, 2011 Lattice Boltzmann Simulation of One Particle Migrating in a Pulsating Flow in Microvessel QIU Bing ( ), 1, TAN Hui-Li ( Û), 2 and
More informationProjective synchronization of a complex network with different fractional order chaos nodes
Projective synchronization of a complex network with different fractional order chaos nodes Wang Ming-Jun( ) a)b), Wang Xing-Yuan( ) a), and Niu Yu-Jun( ) a) a) School of Electronic and Information Engineering,
More informationStationary states and energy cascades in inelastic gases
Los Alamos National Laboratory From the SelectedWorks of Eli Ben-Naim May, 005 Stationary states and energy cascades in inelastic gases E. Ben-Naim, Los Alamos National Laboratory J. Machta, Uniersity
More informationEnergy production rates in fluid mixtures of inelastic rough hard spheres Andrés Santos, 1 Gilberto M. Kremer, 2 and Vicente Garzó 1 1 Universidad de Extremadura, Badajoz (Spain) 2 Universidade id d Federal
More informationOpen boundary conditions in stochastic transport processes with pair-factorized steady states
Open boundary conditions in stochastic transport processes with pair-factorized steady states Hannes Nagel a, Darka Labavić b, Hildegard Meyer-Ortmanns b, Wolfhard Janke a a Institut für Theoretische Physik,
More informationHow Different is the Dynamics of a Lennard-Jones Binary Fluid from One-Component Lennard-Jones Fluid? 1
How Different is the Dynamics of a Lennard-Jones Binary Fluid from One-Component Lennard-Jones Fluid? Takayuki NARUMI and Michio TOKUYAMA Summary We investigate the dynamics of liquids and supercooled
More informationDynamics that trigger/inhibit cluster formation in a one-dimensional granular gas
Physica A 342 (24) 62 68 www.elsevier.com/locate/physa Dynamics that trigger/inhibit cluster formation in a one-dimensional granular gas Jose Miguel Pasini a;1, Patricio Cordero b; ;2 a Department of Theoretical
More informationarxiv: v2 [cond-mat.stat-mech] 20 Sep 2016
Clustering and velocity distributions in granular gases cooling by solid friction arxiv:1609.01114v2 [cond-mat.stat-mech] 20 Sep 2016 by Prasenjit Das, 1 Sanjay Puri, 1 and Moshe Schwartz 2,3 1 School
More informationResearch of Micro-Rectangular-Channel Flow Based on Lattice Boltzmann Method
Research Journal of Applied Sciences, Engineering and Technology 6(14): 50-55, 013 ISSN: 040-7459; e-issn: 040-7467 Maxwell Scientific Organization, 013 Submitted: November 08, 01 Accepted: December 8,
More informationarxiv: v2 [cond-mat.stat-mech] 1 Oct 2009
arxiv:0909.5563v2 [cond-mat.stat-mech] 1 Oct 2009 Granular gas of viscoelastic particles in a homogeneous cooling state. A. S. Bodrova,a, N. V. Brilliantov b,a a Faculty of Physics, M.V.Lomonosov Moscow
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationThe existence of Burnett coefficients in the periodic Lorentz gas
The existence of Burnett coefficients in the periodic Lorentz gas N. I. Chernov and C. P. Dettmann September 14, 2006 Abstract The linear super-burnett coefficient gives corrections to the diffusion equation
More informationKinetic Models for Granular Flow
Journal of Statistical Physics, Vol. 97, Nos. 12, 1999 Kinetic Models for Granular Flow J. Javier Brey, 1 James W. Dufty, 2 and Andre s Santos 3 Received September 28, 1998; final March 29, 1999 The generalization
More informationPhase Diagram of One-Dimensional Bosons in an Array of Local Nonlinear Potentials at Zero Temperature
Commun. Theor. Phys. (Beijing, China) 36 (001) pp. 375 380 c International Academic Publishers Vol. 36, No. 3, September 15, 001 Phase Diagram of One-Dimensional Bosons in an Array of Local Nonlinear Potentials
More informationAvailable online at ScienceDirect. Physics Procedia 57 (2014 ) 77 81
Available online at www.sciencedirect.com ScienceDirect Physics Procedia 57 (204 ) 77 8 27th Annual CSP Workshops on Recent Developments in Computer Simulation Studies in Condensed Matter Physics, CSP
More informationarxiv: v1 [nlin.cd] 28 Nov 2008
Extremum complexity in the monodimensional ideal gas: the piecewise uniform density distribution approximation Xavier Calbet BIFI, Universidad de Zaragoza, E-50009 Zaragoza, Spain. arxiv:0811.4749v1 [nlin.cd]
More informationSynchronization and Bifurcation Analysis in Coupled Networks of Discrete-Time Systems
Commun. Theor. Phys. (Beijing, China) 48 (2007) pp. 871 876 c International Academic Publishers Vol. 48, No. 5, November 15, 2007 Synchronization and Bifurcation Analysis in Coupled Networks of Discrete-Time
More informationDiffusive Transport Enhanced by Thermal Velocity Fluctuations
Diffusive Transport Enhanced by Thermal Velocity Fluctuations Aleksandar Donev 1 Courant Institute, New York University & Alejandro L. Garcia, San Jose State University John B. Bell, Lawrence Berkeley
More informationA MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE
A MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE Shigeo Maruyama and Tatsuto Kimura Department of Mechanical Engineering The University of Tokyo 7-- Hongo, Bunkyo-ku, Tokyo -866,
More informationStatistical Mechanics and Thermodynamics of Small Systems
Statistical Mechanics and Thermodynamics of Small Systems Luca Cerino Advisors: A. Puglisi and A. Vulpiani Final Seminar of PhD course in Physics Cycle XXIX Rome, October, 26 2016 Outline of the talk 1.
More informationNon-equilibrium phenomena and fluctuation relations
Non-equilibrium phenomena and fluctuation relations Lamberto Rondoni Politecnico di Torino Beijing 16 March 2012 http://www.rarenoise.lnl.infn.it/ Outline 1 Background: Local Thermodyamic Equilibrium 2
More informationDissipation of a two-mode squeezed vacuum state in the single-mode amplitude damping channel
Dissipation of a two-mode squeezed vacuum state in the single-mode amplitude damping channel Zhou Nan-Run( ) a), Hu Li-Yun( ) b), and Fan Hong-Yi( ) c) a) Department of Electronic Information Engineering,
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 6 Feb 2004
A statistical model with a standard Gamma distribution Marco Patriarca and Kimmo Kaski Laboratory of Computational Engineering, Helsinki University of Technology, P.O.Box 923, 25 HUT, Finland arxiv:cond-mat/422v
More informationOn the Asymptotic Convergence. of the Transient and Steady State Fluctuation Theorems. Gary Ayton and Denis J. Evans. Research School Of Chemistry
1 On the Asymptotic Convergence of the Transient and Steady State Fluctuation Theorems. Gary Ayton and Denis J. Evans Research School Of Chemistry Australian National University Canberra, ACT 0200 Australia
More informationApplication of the Dissipative Particle Dynamics Method to Ferromagnetic Colloidal Dispersions
Application of the Dissipative Particle Dynamics Method to Ferromagnetic Colloidal Dispersions AKIRA SATOH Professor, Akita Prefectural University, Japan Roy W. Chantrell Professor, York University, U.
More informationTransport properties through double-magnetic-barrier structures in graphene
Chin. Phys. B Vol. 20, No. 7 (20) 077305 Transport properties through double-magnetic-barrier structures in graphene Wang Su-Xin( ) a)b), Li Zhi-Wen( ) a)b), Liu Jian-Jun( ) c), and Li Yu-Xian( ) c) a)
More information19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, 2-7 SEPTEMBER 2007
19 th INTERNATIONAL CONGRESS ON ACOUSTICS MADRID, -7 SEPTEMBER 007 INVESTIGATION OF AMPLITUDE DEPENDENCE ON NONLINEAR ACOUSTICS USING THE DIRECT SIMULATION MONTE CARLO METHOD PACS: 43.5.Ed Hanford, Amanda
More informationMomentum Distribution of a Fragment and Nucleon Removal Cross Section in the Reaction of Halo Nuclei
Commun. Theor. Phys. Beijing, China) 40 2003) pp. 693 698 c International Academic Publishers Vol. 40, No. 6, December 5, 2003 Momentum Distribution of a ragment and Nucleon Removal Cross Section in the
More informationEffects of Interactive Function Forms and Refractoryperiod in a Self-Organized Critical Model Based on Neural Networks
Commun. Theor. Phys. (Beijing, China) 42 (2004) pp. 121 125 c International Academic Publishers Vol. 42, No. 1, July 15, 2004 Effects of Interactive Function Forms and Refractoryperiod in a Self-Organized
More informationDilatancy Transition in a Granular Model. David Aristoff and Charles Radin * Mathematics Department, University of Texas, Austin, TX 78712
Dilatancy Transition in a Granular Model by David Aristoff and Charles Radin * Mathematics Department, University of Texas, Austin, TX 78712 Abstract We introduce a model of granular matter and use a stress
More informationColloidal Suspension Rheology Chapter 1 Study Questions
Colloidal Suspension Rheology Chapter 1 Study Questions 1. What forces act on a single colloidal particle suspended in a flowing fluid? Discuss the dependence of these forces on particle radius. 2. What
More informationAvalanches, transport, and local equilibrium in self-organized criticality
PHYSICAL REVIEW E VOLUME 58, NUMBER 5 NOVEMBER 998 Avalanches, transport, and local equilibrium in self-organized criticality Afshin Montakhab and J. M. Carlson Department of Physics, University of California,
More informationSelf-organized Criticality and Synchronization in a Pulse-coupled Integrate-and-Fire Neuron Model Based on Small World Networks
Commun. Theor. Phys. (Beijing, China) 43 (2005) pp. 466 470 c International Academic Publishers Vol. 43, No. 3, March 15, 2005 Self-organized Criticality and Synchronization in a Pulse-coupled Integrate-and-Fire
More informationSegregation in a noninteracting binary mixture
Segregation in a noninteracting binary mixture Filip Krzyżewski and Magdalena A. Załuska-Kotur,2 Institute of Physics, Polish Academy of Sciences, Al. Lotników 32/46, 02-668 Warsaw, Poland 2 Faculty of
More information2 Abstract In this thesis, the numerical study on the dynamics of granular fluidization in a simple model is presented. The model consists of two-dime
Phase Changes in an Inelastic Hard Disk System with a Heat Bath under Weak Gravity for Granular Fluidization A DISSERTATION SUBMITTED TO KYUSHU UNIVERSITY FOR THE DEGREE OF DOCTOR OF SCIENCE BY MASAHARU
More informationRheology of two- and three-dimensional granular mixtures under uniform shear flow: Enskog kinetic theory versus molecular dynamics simulations
Granular Matter (2006) 8: 103 115 DOI 10.1007/s10035-006-0001-7 ORIGINAL PAPER José María Montanero Vicente Garzó MeheboobAlam Stefan Luding Rheology of two- and three-dimensional granular mixtures under
More informationControl and synchronization of Julia sets of the complex dissipative standard system
Nonlinear Analysis: Modelling and Control, Vol. 21, No. 4, 465 476 ISSN 1392-5113 http://dx.doi.org/10.15388/na.2016.4.3 Control and synchronization of Julia sets of the complex dissipative standard system
More informationA kinetic study of the ordering process in ternary III V semiconductor alloys
J. Phys.: Condens. Matter 9 (1997) 5737 5749. Printed in the UK PII: S0953-8984(97)78738-1 A kinetic study of the ordering process in ternary III V semiconductor alloys Zhi-Feng Huang and Bing-Lin Gu CCAST
More informationAsymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks
Physica A 275 (2000) 483 504 www.elsevier.com/locate/physa Asymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks V. Kumaran Department of Chemical Engineering, Indian
More informationInvaded cluster dynamics for frustrated models
PHYSICAL REVIEW E VOLUME 57, NUMBER 1 JANUARY 1998 Invaded cluster dynamics for frustrated models Giancarlo Franzese, 1, * Vittorio Cataudella, 1, * and Antonio Coniglio 1,2, * 1 INFM, Unità di Napoli,
More informationA comparison of different methods in the study of dynamical fluctuations in high energy e + e collisions
A comparison of different methods in the study of dynamical fluctuations in high energy e + e collisions Chen Gang, and Liu Lianshou Institute of Particle Physics, Huazhong Normal University, Wuhan 79
More informationStatistical Mechanics of Active Matter
Statistical Mechanics of Active Matter Umberto Marini Bettolo Marconi University of Camerino, Italy Naples, 24 May,2017 Umberto Marini Bettolo Marconi (2017) Statistical Mechanics of Active Matter 2017
More information16. Working with the Langevin and Fokker-Planck equations
16. Working with the Langevin and Fokker-Planck equations In the preceding Lecture, we have shown that given a Langevin equation (LE), it is possible to write down an equivalent Fokker-Planck equation
More informationVici Progress Report, March 2013.
Vici Progress Report, March 2013. Nicolás Rivas Project description Two general approaches are used to model granular materials: microscopic, where information is taken from individual particles, and macroscopic,
More informationarxiv:cond-mat/ v2 [cond-mat.stat-mech] 21 Nov 2003
Diffusion of impurities in a granular gas Vicente Garzó Departamento de Física, Universidad de Extremadura, E-0607 Badajoz, Spain arxiv:cond-mat/0309639v [cond-mat.stat-mech] Nov 003 José María Montanero
More informationKINETIC THEORY OF GRANULAR SOME APPLICATIONS
KINETIC THEORY OF GRANULAR BINARY MIXTURES: SOME APPLICATIONS Vicente Garzó Departamento t de Física, Universidad id dde Extremadura, Badajoz, Spain OUTLINE 1. Introduction 2. Granular binary mixtures.
More informationCellular Automaton Simulation of Evacuation Process in Story
Commun. Theor. Phys. (Beijing, China) 49 (2008) pp. 166 170 c Chinese Physical Society Vol. 49, No. 1, January 15, 2008 Cellular Automaton Simulation of Evacuation Process in Story ZHENG Rong-Sen, 1 QIU
More informationMOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE
MOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE Tatsuto Kimura* and Shigeo Maruyama** *Department of Mechanical Engineering, The University of Tokyo 7-- Hongo,
More informationQuantum measurement theory and micro-macro consistency in nonequilibrium statistical mechanics
Nagoya Winter Workshop on Quantum Information, Measurement, and Quantum Foundations (Nagoya, February 18-23, 2010) Quantum measurement theory and micro-macro consistency in nonequilibrium statistical mechanics
More informationLattice Bhatnagar Gross Krook model for the Lorenz attractor
Physica D 154 (2001) 43 50 Lattice Bhatnagar Gross Krook model for the Lorenz attractor Guangwu Yan a,b,,liyuan a a LSEC, Institute of Computational Mathematics, Academy of Mathematics and System Sciences,
More informationDSMC Simulation of Binary Rarefied Gas Flows between Parallel Plates and Comparison to Other Methods
Simulation of Binary Rarefied Gas Flows between Parallel Plates and Comparison to Other Methods L. Szalmas Department of Mechanical Engineering, University of Thessaly, Pedion Areos, Volos 38334, Greece
More informationEffects of Interactive Function Forms in a Self-Organized Critical Model Based on Neural Networks
Commun. Theor. Phys. (Beijing, China) 40 (2003) pp. 607 613 c International Academic Publishers Vol. 40, No. 5, November 15, 2003 Effects of Interactive Function Forms in a Self-Organized Critical Model
More informationFrom time series to superstatistics
From time series to superstatistics Christian Beck School of Mathematical Sciences, Queen Mary, University of London, Mile End Road, London E 4NS, United Kingdom Ezechiel G. D. Cohen The Rockefeller University,
More informationQuasi-Stationary Simulation: the Subcritical Contact Process
Brazilian Journal of Physics, vol. 36, no. 3A, September, 6 685 Quasi-Stationary Simulation: the Subcritical Contact Process Marcelo Martins de Oliveira and Ronald Dickman Departamento de Física, ICEx,
More informationarxiv:cond-mat/ v1 22 Sep 1998
Scaling properties of the cluster distribution of a critical nonequilibrium model Marta Chaves and Maria Augusta Santos Departamento de Física and Centro de Física do Porto, Faculdade de Ciências, Universidade
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
More informationGeneralized continuum theory for ferroelectric thin films
Generalized continuum theory for ferroelectric thin films Tianquan Lü* and Wenwu Cao Department of Physics and Materials Science, City University of Hong Kong, 83 Tat Chee Avenue, Kowloon, Hong Kong, China
More information