Anomalous Collective Diffusion in One-Dimensional Driven Granular Media
|
|
- Sharyl Holland
- 5 years ago
- Views:
Transcription
1 Typeset with jpsj2.cls <ver.1.2> Anomalous Collective Diffusion in One-Dimensional Driven Granular Media Yasuaki Kobayashi and Masaki Sano Department of Physics, University of Tokyo, Hongo, Tokyo We numerically investigate the nonequilibrium steady state of a one-dimensional rarefied granular system in which each particle is driven by Gaussian white noise and subject to inelastic mutual collisions. We find that the dynamic structure factor S(q, t) decays as exp( Dq 2 t 2 ), which is in marked contrast to the usual Fickian behavior exp( Dq 2 t). Employing newly obtained S(q, t), we demonstrate that the static spatial correlation function and the one-point time correlation function obey the same scaling behavior. KEYWORDS: granular media, clustering, dissipation, power law, dynamic structure factor, collective diffusion 1. Introduction Simple driven dissipative systems have been widely investigated to gain insights into a variety of phenomena found in systems far from equilibrium. For example, systems that are subjected to a driving field such as a shear flow or an electric field often exhibit various self-organized structures, which never occur in equilibrium states. Above all, granular media have been the subject of much interest. Among many intriguing properties of a granular gas, one of the most striking features is the clustering; Goldhirsch and Zanetti 1) showed that a gas composed of inelastically colliding particles loses its stability and forms high density clusters. Another example is dynamic properties of granular media passing in a narrow pipe. In the experiment of granular pipe flow, the power law decay in density power spectra was obtained. 2 7) Lattice Gas Automata have also confirmed such power law behavior. 8, 9) Granular media with continuous energy input evolve into a nonequilibrium steady state in which much interesting properties have been studied such as an equation of state, 10) static spatial correlations, 10 12) velocity correlations, 11, 12) and a large scale spatial structure. 13) It is worth noting that in such a nonequilibrium system clusters are always unstable and the system repeats ceaseless creation and destruction of the clusters. To fully understand the nonequilibrium properties in such a sistem, we should pay attention to the time-dependent structures of fluctuations, which are absent in most of the previous studies. In this paper, we numerically study the one-dimensional driven granular system in detail from the dynamical viewpoint. Unexpectedly, we find that the collective diffusion is anomalous address: yasuaki@daisy.phys.s.u-tokyo.ac.jp 1/11
2 where the usual Fickian diffusion mechanism breaks down. However, as we will demonstrate below, the unusual finding will be reasonably understood in terms of a relation between the static spatial correlation function and the one-point time correlation function, both of which are also calculated in the present simulations. 2. Model Oursystem iscomposed ofn identical particles of diameter d distributedon alineoflength L. Mutual passages of the particles are prohibited so that their order remains unchanged. Each particle is driven by Gaussian random forces which are uncorrelated spatially and temporally. Energy input by random forces can be balanced with the dissipation due to inelastic mutual collisions and frictional forces. The equation of motions for the particles are given by ẍ i = γẋ i + Θ ξ i (t) [ ( θ(d r ij ) κ r ij r ) ] ij r ij d ηṙ ij, (1) j=i±1 where x i is the position of i-th particle, θ(x) is the Heaviside step function, r ij = x i x j,θis the strength of the noise and γ is the drag force coefficient. The random driving force ξ i (t) is assumed to be the Gaussian white noise with zero mean and satisfy ξ i (t)ξ j (t ) = δ ij δ(t t ). It should be emphasized that ξ i (t) is not a source of fluctuation which ensures the fluctuationdissipation relation around the thermal equilibrium state but just an external force in a nonequilibrium steady state. This way of energy input was originally introduced by Williams and MacKintosh. 10) The last term in eq. (1) represents inelastic collisions between particles. 14) κ is the elastic constant and η is the friction coefficient. κ and η are determined from the coefficient of restitution r and the contact time t col during a collision by κ =(π 2 +(lnr) 2 )/2t 2 col and η = ln r/t col. A prohibition rule of mutual passages of the particles is not explicitly included in our model. By choosing proper values of κ and η, however, we can avoid the mutual passages. We chose d as the unit length and t col as the unit time. As for the other parameters, we set N = 100, L = 10000, Θ = and γ =10 4, whereby we limited our concern to rarefied gas with low drag force. We have verified that the mutual passage do not occur in the range of r used in our simulations. To integrate eq. (1) we use the explicit Euler scheme with the time step t =0.01. This value of the time step provides the good accuracy and efficiency of our calculation, and we have also confirmed that using smaller values of t do not change our results. Throughout our simulations, periodic boundary conditions were imposed. In order to decide that the system is in the steady state, we calculated the total kinetic energy E = N j=1 ẋ2 j /2 and verified that it conveged to the constant value. We also confirmed that the two-point static spatial correlation function (see below) converged at the same time. 2/11
3 All of our simulations use the data after this transient period. 3. Simulation Results Fig. 1. Time evolution of a distribution of particles in the steady state. (a)r =0.9; (b)r =0.01.The horizontal axis and the vertical one denote the space and the time respectively. The direction of thetimeaxesisfromtop(t = t 0 ) to bottom (t = t ). Figure 1 shows the typical space-time plots in the steady state for r = 0.9(a) and r =0.01(b). As is well known, and also supported by our simulations, for small r the density profile shows strong inhomogeneity and the particles form clusters due to the inelastic collisions. In addition, in Fig. 1(b) the system exhibits a much more dynamic behavior; each cluster moves around and repeats coalition and breakdown. In order to study such highly inhomogeneous spatial structures more quantitatively, we consider the correlation properties of this system. The spatial distribution can be characterized by the two-point static spatial 3/11
4 0.1 r=0.001 r=0.5 r=0.9 g(x,0) 0.01 x x Fig. 2. Plots of the two-point correlation function g(x, 0) versus distance for several values of r. The fitting lines represent g(x, 0) x α with α =0.49 for r =0.001, α =0.32 for r =0.5, and α =0.009 for r =0.9. correlation function g(x, 0), which is defined by g(x, 0) = n(x 0,t 0 )n(x 0 + x, t 0 ), (2) where n(x, t) = 1 N N j=1 δ(x x j(t)) is the number density of the particle at position x at time t, andthebrackets denote the averages with respect to the random force ξ i and the reference position x 0. Figure 2 shows that g(x, 0) follows the power law with its exponent close to 0.5 as r approaches 0. A similar power law behavior has also been obtained for an inviscid case. 10) This finding suggest that the density fluctuation becomes long-ranged and the system is strongly correlated. Nevertheless, the equal-time spatial correlation is not sufficient to describe the dynamical aspect of the system, which is our main concern. We thus calculated the one-point time correlation function g(0,t) defined by g(0,t)= n(x 0,t 0 )n(x 0,t 0 + t). (3) Figure 3 shows that g(0,t) also clearly exhibits the power law behavior for small r. Asr 0, the exponent approaches 0.5, which is the same value as obtained for g(x, 0). Letting g(x, 0) x α and g(0,t) t β, we show the dependence of the exponents on r showninfig.4. To further characterize the dynamics of the density fluctuation we calculated the dynam- 4/11
5 1 r=0.001 r=0.1 r=0.9 g(0,t) t Fig. 3. Plots of the one-point time correlation function g(0,t) for several values of r. The fitting line represents g(0,t) t β with β =0.45 for r = ical structure factor S(q,t), which determines the relaxation of density fluctuations. S(q,t) is defined as S(q,t) = n(q,t)n( q,0) 1 = N j 0 e iq(x j(t) x 0 (0)), (4) where n(q,t) denotes the spatial Fourier transform of n(x, t), and the brackets are the same as that of eq.(2). In the following we set r =0.01. Figure 5 shows the time dependence of S(q,t) for several values of q; q = 10q min, q = 15q min,an = 20q min,whereq min = 2π/L. The plotted data are well fitted by exp( D(q)t 2 )atleastuntils(q,t) diminishes by half. Therefore the main process of the relaxation of fluctuations can be described by S(q,t) exp( D(q)t 2 ). D(q) are plotted in Fig. 6 as a function of q 2. The plotted data were well fitted by the line D(q) =Dq 2.From these results we can determine the form of S(q,t) ass(q,t) =S(q,0) exp( Dq 2 t 2 ). This result is in marked contrast to the usual Fickian behavior S(q,t) exp( D c q 2 t), where D q is the collective diffusion constant. This implies anomalous collective diffusion. The characteristic relaxation time τ of a cluster with size l is defined as τ = l/ D, so the relaxation of 5/11
6 Fig. 4. Dependence of α and β on r, whereα and β are the exponent of the correlation functions g(x, 0) and g(0,t) respectively. As r 0, the two exponents approach the same value, 1/2. fluctuations is not diffusive but ballistic. 4. Discussion It is surprising that the collective motions of the particles are not diffusive but rather ballistic in the steady state. We also studied an individual motion of a tagged particle in the steady state with r =0.01. The mean square distance defined as x(t) 2 = 1 N N i=1 (x i(t + t 0 ) x(t 0 )) 2,where is the same as before, is shown in Fig. 7. This shows that particles lose their ballistic property at τ c 600, because of the prohibition of mutual passages in one dimensional system. 16, 17) But the anomaly of the collective diffusion continues beyond τ c (see Fig. 5). Therefore, the ballistic collective behavior cannot be attributed to the individual particle motions. The origin of this anomaly is still under investigation. The finding successfully explains why g(x, 0) and g(0,t) exhibits the same scaling behavior. By definition, g(x, t) = Using the form of S(q,t) obtained above, we obtain g(0,t)= 2π S(q,t) = 2π S(q,t)e iqx. (5) t 2 2π S(q,0)e Dq2 6/11
7 S(q,t) t Fig. 5. Plots of dynamical structure factor S(q, t) averaged over configurations in the steady state. Each curve corresponds to q =10q min, q =15q min,an =20q min (q min =2π/L) from top to bottom. Fitting curves have the function form exp( D(q)t 2 ). = 2π dx g(x, 0)e iqx e Dq2 t 2. (6) Substituting the asymptotic form g(x, 0) x α into the above equation, we obtain g(0,t) 2π dx x α e iqx e Dq2 t 2 2π qα 1 e Dq2 t 2 const. t α const. (7) Comparing the asymptotic form of the time correlation function, g(0,t) t β,weobtain α = β. This indicates that the spatial correlations has a close relation to the time correlations via the anomalous collective diffusion. We believe that such a consideration is quite reasonable, because nonequilibrium steady states generally involve various relaxational processes of their internal degrees of freedom, and cannot be characterized only by time-independent properties. The ballistic collective diffusion has been further confirmed by another simulation result. Figure 8 shows the spatial distribution of particles initially close-packed in a small region with their center of mass at the origin. All the parameters are the same as before, except that the system size is large enough for the particles not to reach an end of the system. 7/11
8 -6 8x 6x D(q) 4x 2x (q/qmin) 2 Fig. 6. Plots of D(q) onq. The fitting line represents D(q) =Dq 2 with D = Rescaled as x/t x and tn n, all the three plots at the different moments sit on the same curve n(x, t) = 1 2D t exp( x /D t). This indicates again that the lifetime of the cluster is proportional to its size, or the relaxation process is ballistic. In conclusion, we numerically investigated the nonequilibrium steady state of a onedimensional driven granular system and found that the dynamic structure factor S(q,t) decayed as exp( Dq 2 t 2 ), as opposed to the usual Fickian behavior exp( Dq 2 t). This causes the coincidence of the scaling exponent between the spatial density correlation and the temporal density correlation. Acknowledgement We would like to thank Yoshihiro Murayama and Hirofumi Wada for fruitful discussions. 8/11
9 Fig. 7. The mean square distance of a particle. The two straight lines represent t and t 2, respectively. 9/11
10 Fig. 8. Density distribution of particles initially close-packed in a region with their center at the origin, rescaled as x/t x and tp (x/t) P (x). +:t = 1000; :t = 2000; + :t = The fitting curve is A exp( x /b) witha = and b = /11
11 References 1) I. Goldhirsch and G. Zanetti: Phys. Rev. Lett. 70 (1993) ) G. W. Baxter and R. P. Behringer: Phys. Rev. Lett. 62 (1989) ) S. Horikawa, A. Nakahara, T. Nakayama and M. Matsushita: J. Phys. Soc. Jpn. 64 (1995) ) O. Moriyama, N. Kuroiwa, M. Matsushita and H. Hayakawa: Phys. Rev. Lett. 80 (1998) ) A. Nakahara and T. Isoda: Phys. Rev. E 55 (1997) ) S. Horikawa, T. Isoda, T. Nakayama, A. Nakahara and M. Matsushita: Physica A 233 (1996) ) O. Moriyama, N. Kuroiwa, S. Tateda, T. Arai, A. Awazu, Y. Yamazaki and M. Matsushita: Prog. Theor. Phys. 150 (2003) ) G.PengandH.J.Herrmann:Phys.Rev.E49 (1994) R ) G.PengandH.J.Herrmann:Phys.Rev.E51 (1995) ) D. R. M. Williams and F. C. MacKintosh: Phys. Rev. E. 54 (1996) R9. 11) A. Puglisi, V. Loreto, U. M. B. Marconi, A. Petri, and A. Vulpiani: Phys. Rev. Lett. 81 (1998) ) G. Peng and T. Ohta: Phys. Rev. E 58 (1998) ) T. P. C. van Noije, M. H. Ernst, E. Trizac and I. Pagonabarraga: Phys. Rev. E 59 (1999) ) Y-h. Taguchi: Phys. Rev. Lett. 69 (1992) ) Y. Murayama and M. Sano: J. Phys. Soc. Jpn. 67 (1998) ) D. G. Levitt, Phys. Rev. A 8 (1973) ) Q.-H. Wei, C. Bechinger, and P. Leiderer, Science, 287 (2000) /11
Collision Statistics of Driven Polydisperse Granular Gases
Commun. Theor. Phys. (Beijing, China) 49 (2008) pp. 1333 1338 c Chinese Physical Society Vol. 49, No. 5, May 15, 2008 Collision Statistics of Driven Polydisperse Granular Gases CHEN Zhi-Yuan, 1,2, ZHANG
More informationVIII.B Equilibrium Dynamics of a Field
VIII.B Equilibrium Dynamics of a Field The next step is to generalize the Langevin formalism to a collection of degrees of freedom, most conveniently described by a continuous field. Let us consider the
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 informationarxiv: v2 [cond-mat.soft] 8 Jun 2010
Model for Density Waves in Gravity-Driven Granular Flow in Narrow Pipes Simen Å. Ellingsen Department of Energy and Process Engineering, Norwegian University of Science and Technology, N-7491 Trondheim,
More informationVelocity 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 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 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 informationLow-Resistant Band-Passing Noise and Its Dynamical Effects
Commun. Theor. Phys. (Beijing, China) 48 (7) pp. 11 116 c International Academic Publishers Vol. 48, No. 1, July 15, 7 Low-Resistant Band-Passing Noise and Its Dynamical Effects BAI Zhan-Wu Department
More informationOn self-organised criticality in one dimension
On self-organised criticality in one dimension Kim Christensen Imperial College ondon Department of Physics Prince Consort Road SW7 2BW ondon United Kingdom Abstract In critical phenomena, many of the
More informationDiffusion in Reduced Dimensions. Clemens Bechinger 2. Physikalisches Institut, Universität Stuttgart
Diffusion in Reduced Dimensions Clemens Bechinger 2. Physikalisches Institut, Universität Stuttgart Diffusion in narrow Channels t = 0 t = t 1 3D, 2D: mixing 1D: sequence unchanged 1D diffusion entirely
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 informationAggregated spots and Waving Loops in a Reaction-Advection-Diffusion System with a Global Coupling. Abstract
Aggregated spots and Waving Loops in a Reaction-Advection-Diffusion System with a Global Coupling Yasuaki Kobayashi and Masaki Sano Department of Physics, The University of Tokyo, Hongo, Tokyo 113-0033,
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 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 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 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 informationarxiv:chao-dyn/ v1 30 Jan 1997
Universal Scaling Properties in Large Assemblies of Simple Dynamical Units Driven by Long-Wave Random Forcing Yoshiki Kuramoto and Hiroya Nakao arxiv:chao-dyn/9701027v1 30 Jan 1997 Department of Physics,
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 informationPhysics 562: Statistical Mechanics Spring 2003, James P. Sethna Homework 5, due Wednesday, April 2 Latest revision: April 4, 2003, 8:53 am
Physics 562: Statistical Mechanics Spring 2003, James P. Sethna Homework 5, due Wednesday, April 2 Latest revision: April 4, 2003, 8:53 am Reading David Chandler, Introduction to Modern Statistical Mechanics,
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 informationPhase Transitions in Nonequilibrium Steady States and Power Laws and Scaling Functions of Surface Growth Processes
Phase Transitions in Nonequilibrium Steady States and Power Laws and Scaling Functions of Surface Growth Processes Term-paper for PHY563 Xianfeng Rui, UIUC Physics Abstract: Three models of surface growth
More informationarxiv:cond-mat/ v1 17 Aug 1994
Universality in the One-Dimensional Self-Organized Critical Forest-Fire Model Barbara Drossel, Siegfried Clar, and Franz Schwabl Institut für Theoretische Physik, arxiv:cond-mat/9408046v1 17 Aug 1994 Physik-Department
More informationNon-equilibrium phase transitions
Non-equilibrium phase transitions An Introduction Lecture III Haye Hinrichsen University of Würzburg, Germany March 2006 Third Lecture: Outline 1 Directed Percolation Scaling Theory Langevin Equation 2
More information2. Molecules in Motion
2. Molecules in Motion Kinetic Theory of Gases (microscopic viewpoint) assumptions (1) particles of mass m and diameter d; ceaseless random motion (2) dilute gas: d λ, λ = mean free path = average distance
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 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 informationLangevin Methods. Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 10 D Mainz Germany
Langevin Methods Burkhard Dünweg Max Planck Institute for Polymer Research Ackermannweg 1 D 55128 Mainz Germany Motivation Original idea: Fast and slow degrees of freedom Example: Brownian motion Replace
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 informationBSc/MSci EXAMINATION. Vibrations and Waves. Date: 4 th May, Time: 14:30-17:00
BSc/MSci EXAMINATION PHY-217 Vibrations and Waves Time Allowed: 2 hours 30 minutes Date: 4 th May, 2011 Time: 14:30-17:00 Instructions: Answer ALL questions in section A. Answer ONLY TWO questions from
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 informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
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 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 information1. Introduction to Chemical Kinetics
1. Introduction to Chemical Kinetics objectives of chemical kinetics 1) Determine empirical rate laws H 2 + I 2 2HI How does the concentration of H 2, I 2, and HI change with time? 2) Determine the mechanism
More informationK. Pyragas* Semiconductor Physics Institute, LT-2600 Vilnius, Lithuania Received 19 March 1998
PHYSICAL REVIEW E VOLUME 58, NUMBER 3 SEPTEMBER 998 Synchronization of coupled time-delay systems: Analytical estimations K. Pyragas* Semiconductor Physics Institute, LT-26 Vilnius, Lithuania Received
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 informationarxiv:nlin/ v1 [nlin.ps] 4 Sep 2004
arxiv:nlin/0409008v1 [nlin.ps] 4 Sep 2004 Classification of KPZQ and BDP models by multiaffine analysis Hiroaki Katsuragi and Haruo Honjo Department of Applied Science for Electronics and Materials, Interdisciplinary
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 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 informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationDynamics of net-baryon density correlations near the QCD critical point
Dynamics of net-baryon density correlations near the QCD critical point Marcus Bluhm and Marlene Nahrgang The work of M.B. is funded by the European Union s Horizon 22 research and innovation programme
More informationKnots and Random Walks in Vibrated Granular Chains
Los Alamos National Laboratory From the SelectedWorks of Eli Ben-Naim February, 2 Knots and Random Walks in Vibrated Granular Chains E. Ben-Naim, Los Alamos National Laboratory Z.A. Daya, Los Alamos National
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 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 informationClassical Statistical Mechanics: Part 1
Classical Statistical Mechanics: Part 1 January 16, 2013 Classical Mechanics 1-Dimensional system with 1 particle of mass m Newton s equations of motion for position x(t) and momentum p(t): ẋ(t) dx p =
More informationClusters and Percolation
Chapter 6 Clusters and Percolation c 2012 by W. Klein, Harvey Gould, and Jan Tobochnik 5 November 2012 6.1 Introduction In this chapter we continue our investigation of nucleation near the spinodal. We
More informationMelting of Ultrathin Lubricant Film Due to Dissipative Heating of Friction Surfaces
ISSN 6-78, Technical Physics, 7, Vol. 5, No. 9, pp. 9. Pleiades Publishing, Ltd., 7. Original Russian Text.V. Khomenko, I.. Lyashenko, 7, published in Zhurnal Tekhnicheskoœ Fiziki, 7, Vol. 77, No. 9, pp.
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 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 informationThe effect of mutual angular misalignment in the quantized sliding of solid lubricants
Facoltà di Scienze e Tecnologie Laurea Triennale in Fisica The effect of mutual angular misalignment in the quantized sliding of solid lubricants Relatore: Prof. Nicola Manini Correlatore: Prof. Rosario
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: v1 [nlin.ps] 9 May 2015
Scaling properties of generalized two-dimensional Kuramoto-Sivashinsky equations V. Juknevičius Institute of Theoretical Physics and Astronomy, Vilnius University, A. Goštauto 2, LT-008 Vilnius, Lithuania
More informationCombined Influence of Off-diagonal System Tensors and Potential Valley Returning of Optimal Path
Commun. Theor. Phys. (Beijing, China) 54 (2010) pp. 866 870 c Chinese Physical Society and IOP Publishing Ltd Vol. 54, No. 5, November 15, 2010 Combined Influence of Off-diagonal System Tensors and Potential
More informationarxiv:cond-mat/ v2 [cond-mat.soft] 28 Mar 2007
Universal Anisotropy in Force Networks under Shear Srdjan Ostojic, 1, Thijs J. H. Vlugt, 2 and Bernard Nienhuis 1 arxiv:cond-mat/0610483v2 [cond-mat.soft] 28 Mar 2007 1 Institute for Theoretical Physics,
More informationFluid-Particles Interaction Models Asymptotics, Theory and Numerics I
Fluid-Particles Interaction Models Asymptotics, Theory and Numerics I J. A. Carrillo collaborators: T. Goudon (Lille), P. Lafitte (Lille) and F. Vecil (UAB) (CPDE 2005), (JCP, 2008), (JSC, 2008) ICREA
More informationCoarsening process in the 2d voter model
Alessandro Tartaglia (LPTHE) Coarsening in the 2d voter model May 8, 2015 1 / 34 Coarsening process in the 2d voter model Alessandro Tartaglia LPTHE, Université Pierre et Marie Curie alessandro.tartaglia91@gmail.com
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 informationActive and Driven Soft Matter Lecture 3: Self-Propelled Hard Rods
A Tutorial: From Langevin equation to Active and Driven Soft Matter Lecture 3: Self-Propelled Hard Rods M. Cristina Marchetti Syracuse University Boulder School 2009 A Tutorial: From Langevin equation
More informationcooperative motion in sheared granular matter Takahiro Hatano
cooperative motion in sheared granular matter Takahiro Hatano (Earthquake Research Institute, University of Tokyo) amorphous particulate systems: structure? 2D granular matter close to jamming spontaneous
More informationTheory of bifurcation amplifiers utilizing the nonlinear dynamical response of an optically damped mechanical oscillator
Theory of bifurcation amplifiers utilizing the nonlinear dynamical response of an optically damped mechanical oscillator Research on optomechanical systems is of relevance to gravitational wave detection
More informationLinear and Nonlinear Oscillators (Lecture 2)
Linear and Nonlinear Oscillators (Lecture 2) January 25, 2016 7/441 Lecture outline A simple model of a linear oscillator lies in the foundation of many physical phenomena in accelerator dynamics. A typical
More informationIntroduction to a few basic concepts in thermoelectricity
Introduction to a few basic concepts in thermoelectricity Giuliano Benenti Center for Nonlinear and Complex Systems Univ. Insubria, Como, Italy 1 Irreversible thermodynamic Irreversible thermodynamics
More informationPHYSICAL REVIEW LETTERS
PHYSICAL REVIEW LETTERS VOLUME 80 1 JUNE 1998 NUMBER 22 Field-Induced Stabilization of Activation Processes N. G. Stocks* and R. Mannella Dipartimento di Fisica, Università di Pisa, and Istituto Nazionale
More information73 L. L. Gonοcalves et al. ments, each one containing n sites occupied by isotopes characterized by n dierent spin relaxation times. The Hamiltonian f
Brazilian Journal of Physics, vol. 3, no. 4, December, 73 Nagel Scaling and Relaxation in the Kinetic Ising Model on an n-isotopic Chain L. L. Gonοcalves, Λ Departamento de Fisica, Universidade Federal
More informationarxiv:cond-mat/ v4 [cond-mat.stat-mech] 19 Jun 2007
arxiv:cond-mat/060065v4 [cond-mat.stat-mech] 9 Jun 007 Restoration of Isotropy in the Ising Model on the Sierpiński Gasket Naoto Yajima Graduate School of Human and Environmental Studies, Kyoto University,
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 informationNonlinear Analysis: Modelling and Control, Vilnius, IMI, 1998, No 3 KINK-EXCITATION OF N-SYSTEM UNDER SPATIO -TEMPORAL NOISE. R.
Nonlinear Analysis: Modelling and Control Vilnius IMI 1998 No 3 KINK-EXCITATION OF N-SYSTEM UNDER SPATIO -TEMPORAL NOISE R. Bakanas Semiconductor Physics Institute Go štauto 11 6 Vilnius Lithuania Vilnius
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 informationwhere r n = dn+1 x(t)
Random Variables Overview Probability Random variables Transforms of pdfs Moments and cumulants Useful distributions Random vectors Linear transformations of random vectors The multivariate normal distribution
More informationModified models of polymer phase separation
Modified models of polymer phase separation Douglas Zhou, Pingwen Zhang,, * and Weinan E,2, LMAM, Peking University, Beijing 0087, People s Republic of China and School of Mathematical Science, Peking
More informationStochastic equations for thermodynamics
J. Chem. Soc., Faraday Trans. 93 (1997) 1751-1753 [arxiv 1503.09171] Stochastic equations for thermodynamics Roumen Tsekov Department of Physical Chemistry, University of Sofia, 1164 Sofia, ulgaria The
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 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 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 informationarxiv:cond-mat/ v2 [cond-mat.stat-mech] 8 Sep 1999
BARI-TH 347/99 arxiv:cond-mat/9907149v2 [cond-mat.stat-mech] 8 Sep 1999 PHASE ORDERING IN CHAOTIC MAP LATTICES WITH CONSERVED DYNAMICS Leonardo Angelini, Mario Pellicoro, and Sebastiano Stramaglia Dipartimento
More informationModeling of Micro-Fluidics by a Dissipative Particle Dynamics Method. Justyna Czerwinska
Modeling of Micro-Fluidics by a Dissipative Particle Dynamics Method Justyna Czerwinska Scales and Physical Models years Time hours Engineering Design Limit Process Design minutes Continious Mechanics
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 informationNetwork formation in viscoelastic phase separation
INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICS: CONDENSED MATTER J. Phys.: Condens. Matter 15 (2003) S387 S393 PII: S0953-8984(03)54761-0 Network formation in viscoelastic phase separation Hajime Tanaka,
More informationStable One-Dimensional Dissipative Solitons in Complex Cubic-Quintic Ginzburg Landau Equation
Vol. 112 (2007) ACTA PHYSICA POLONICA A No. 5 Proceedings of the International School and Conference on Optics and Optical Materials, ISCOM07, Belgrade, Serbia, September 3 7, 2007 Stable One-Dimensional
More informationClusters in granular flows : elements of a non-local rheology
CEA-Saclay/SPEC; Group Instabilities and Turbulence Clusters in granular flows : elements of a non-local rheology Olivier Dauchot together with many contributors: D. Bonamy, E. Bertin, S. Deboeuf, B. Andreotti,
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 informationthe renormalization group (RG) idea
the renormalization group (RG) idea Block Spin Partition function Z =Tr s e H. block spin transformation (majority rule) T (s, if s i ; s,...,s 9 )= s i > 0; 0, otherwise. b Block Spin (block-)transformed
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 informationKibble-Zurek dynamics and off-equilibrium scaling of critical cumulants in the QCD phase diagram
Kibble-Zurek dynamics and off-equilibrium scaling of critical cumulants in the QCD phase diagram Raju Venugopalan BNL/Heidelberg arxiv:1310.1600 [cond-mat.stat-mech] Heidelberg Seminar, June 8 th, 2016
More informationDimensionality influence on energy, enstrophy and passive scalar transport.
Dimensionality influence on energy, enstrophy and passive scalar transport. M. Iovieno, L. Ducasse, S. Di Savino, L. Gallana, D. Tordella 1 The advection of a passive substance by a turbulent flow is important
More informationCover times of random searches
Cover times of random searches by Chupeau et al. NATURE PHYSICS www.nature.com/naturephysics 1 I. DEFINITIONS AND DERIVATION OF THE COVER TIME DISTRIBUTION We consider a random walker moving on a network
More informationMetropolis, 2D Ising model
Metropolis, 2D Ising model You can get visual understanding from the java applets available, like: http://physics.ucsc.edu/~peter/ising/ising.html Average value of spin is magnetization. Abs of this as
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 informationInhomogeneous elastic response of amorphous solids
Inhomogeneous elastic response of amorphous solids Jean-Louis Barrat Université de Lyon Institut Universitaire de France Acknowledgements: Anne Tanguy, Fabien Chay Goldenberg, Léonforte, Michel Tsamados
More informationChapter 2 Ensemble Theory in Statistical Physics: Free Energy Potential
Chapter Ensemble Theory in Statistical Physics: Free Energy Potential Abstract In this chapter, we discuss the basic formalism of statistical physics Also, we consider in detail the concept of the free
More informationThe Brownian Oscillator
Chapter The Brownian Oscillator Contents.One-DimensionalDiffusioninaHarmonicPotential...38 The one-dimensional Smoluchowski equation, in case that a stationary flux-free equilibrium state p o (x) exists,
More informationThe behaviour of high Reynolds flows in a driven cavity
The behaviour of high Reynolds flows in a driven cavity Charles-Henri BRUNEAU and Mazen SAAD Mathématiques Appliquées de Bordeaux, Université Bordeaux 1 CNRS UMR 5466, INRIA team MC 351 cours de la Libération,
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 informationCompacton-like solutions in some nonlocal hydrodynamic-type models
Compacton-like solutions in some nonlocal hydrodynamic-type models Vsevolod Vladimirov AGH University of Science and technology, Faculty of Applied Mathematics Protaras, October 26, 2008 WMS AGH Compactons
More informationEVOLUTION OF COMPLEX FOOD WEB STRUCTURE BASED ON MASS EXTINCTION
EVOLUTION OF COMPLEX FOOD WEB STRUCTURE BASED ON MASS EXTINCTION Kenichi Nakazato Nagoya University Graduate School of Human Informatics nakazato@create.human.nagoya-u.ac.jp Takaya Arita Nagoya University
More informationSynchronization of Limit Cycle Oscillators by Telegraph Noise. arxiv: v1 [cond-mat.stat-mech] 5 Aug 2014
Synchronization of Limit Cycle Oscillators by Telegraph Noise Denis S. Goldobin arxiv:148.135v1 [cond-mat.stat-mech] 5 Aug 214 Department of Physics, University of Potsdam, Postfach 61553, D-14415 Potsdam,
More informationSimple approach to the creation of a strange nonchaotic attractor in any chaotic system
PHYSICAL REVIEW E VOLUME 59, NUMBER 5 MAY 1999 Simple approach to the creation of a strange nonchaotic attractor in any chaotic system J. W. Shuai 1, * and K. W. Wong 2, 1 Department of Biomedical Engineering,
More informationGranular materials (Assemblies of particles with dissipation )
Granular materials (Assemblies of particles with dissipation ) Saturn ring Sand mustard seed Ginkaku-ji temple Sheared granular materials packing fraction : Φ Inhomogeneous flow Gas (Φ = 012) Homogeneous
More informationThis is a Gaussian probability centered around m = 0 (the most probable and mean position is the origin) and the mean square displacement m 2 = n,or
Physics 7b: Statistical Mechanics Brownian Motion Brownian motion is the motion of a particle due to the buffeting by the molecules in a gas or liquid. The particle must be small enough that the effects
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More information