arxiv:cond-mat/ v3 [cond-mat.stat-mech] 23 Oct 1997

Transcription

1 Dynamcs of nelastcally collng rough spheres: Relaxaton of translatonal an rotatonal energy arxv:con-mat/ v3 [con-mat.stat-mech] 23 Oct 1997 Martn Huthmann an Annette Zppelus Insttut für Theoretsche Physk, Georg August Unverstät, Göttngen, Germany (January 31, 2018) We stuy the exchange of knetc energy between translatonal an rotatonal egrees of freeom for nelastc collsons of rough spheres. Even f equpartton hols n the ntal state t s mmeately estroye by collsons. The smplest generalsaton of the homogeneous coolng state allows for two temperatures, characterzng translatonal an rotatonal egrees of freeom separately. For tmes larger than a crossover frequency, whch s etermne by the Enskog frequency an the ntal temperature, both energes ecay algebracally lke t 2 wth a fxe rato of ampltues, fferent from one. Knetc theory of nelastcally collng partcles has become a subject of growng research actvty, motvate partly by renewe nterest n granular materals. A Boltzmann equaton an the Enskog varant of t have been formulate for nelastcally collng partcles wth normal an tangental resttuton [1 6]. Whereas the ervaton of the knetc equaton s relatvely straghtforwar, the methos of soluton whch were eveloppe for elastc collsons cannot be taken over to the nelastc case because there s no smple statonary or local equlbrum strbuton, aroun whch one coul expan the nonlnear knetc equaton. One smple strbuton s the homogeneous coolng state (HCS) [7,8], whch epens on tme only mplctly va the average knetc energy T(t). The latter s precte to ecay lke t 2 for large tmes. The homogenety assumpton s certanly volate when clusterng occurs an, n partcular, f nelastc collapse happens. Nevertheless we have foun recently [9] n a moel of nelastcally collng ros that the knetc energy follows on average a t 2 behavour, even f clusterng occurs. Ths suggests that the above scalng law may be a useful approxmate escrpton, even when the assumptons of HCS break own. In ths note we show that rotatonal an translatonal energy relax wth fferent rates n general. Once frcton s nclue, HCS wth a sngle tme epenent temperature s no longer consstent wth the tme evoluton of translatonal an rotatonal energy separately. Instea one has to ntrouce two temperatures whch characterze translatonal an rotatonal egress of freeom separately. Both are foun to fall off lke t 2 wth the rato approachng a constant value, whch s etermne by the coeffcents of normal an tangental resttuton. We brefly recall the collson ynamcs of har spheres wth normal an tangental resttuton. These results can for example be foun n Cergnan [6]. We conser two spheres of equal ameter a, mass M an moment of nerta I. The unt-vector from the center of the frst sphere to the center of the secon s enote by ˆn an veloctes an angular veloctes before collson by v 1, v 2, ω 1 an ω 2. The relatve velocty of the contactpont before collson s gven by V = v 2 + a 2 ˆn ω 2 v 1 + a 2 ˆn ω 1. Normal an tangental resttuton etermne the relatve velocty after collson accorng to ˆn V = ǫ(ˆn V) wth ǫ [0,1], (1) ˆn V = β(ˆn V ) wth β [ 1,1]. (2) Usng the property of conserve lnear an angular momenta one obtans for the veloctes after collson v 1 = v 1 η t v 12 (η n η t )(ˆn v 12 )ˆn+η t a 2ˆn (ω 1 +ω 2 ), v 2 = v 2 +η t v 12 +(η n η t )(ˆn v 12 )ˆn η t a 2ˆn (ω 1 +ω 2 ), ω 1 = ω 1 2 ak η tˆn v 12 + η t k ˆn (ˆn (ω 1 +ω 2 )), ω 2 = ω 2 2 ak η tˆn v 12 + η t k ˆn (ˆn (ω 1 +ω 2 )). (3) wth v 12 = v 1 v 2 an parameters k := 4I Ma 2 (k=0.4 for homogeneous spheres), η n := 1+ǫ 2, an η t := k(1+β) 2(1+k). We conser a system of N classcal partcles, confne to a 3-mensonal volume V an nteractng va a har core potental. Each partcle s characterze by ts poston r (t), ts lnear momentum p (t) = Mv (t) an angular velocty ω (t). The tme evelopment of a ynamcal varable A = A({r (t),p (t),ω (t)}) s etermne by a Pseuo Louvlle operator L + A({r,p,ω },t) = exp(l + t)a({r,p,ω },0) for t > 0. (4) 1

2 Such a Pseuo Louvlle operator was frst ntrouce for har (perfectly smooth) spheres by Ernst et al. [10] an subsequently apple to the calculaton of transport coeffcents of a har sphere flu (see e.g. [11] an references theren). Noje an Ernst [12] have generalse the formalsm to nelastc collsons wth normal resttuton. Here we exten these results to rough spheres wth normal an tangental resttuton. The Pseuo Louvlle operators L + conssts of two parts L + = L 0 +L +. The frst one, L 0, escrbes free streamng of partcles an the secon one, L + = 1 2 L 0 = p M n rn (5) n n m T +(nm), escrbes har-core collsons of two partcles by T + (nm) = M (p nm ˆr nm )Θ( p nm ˆr nm )δ( r nm a)(b + nm 1). (6) The operatorb + nm replacesthe lnearan angularmomentaoftwo partclesnan m before collsonby the corresponng ones after collson. Θ(x) s the Heavse step functon an we have ntrouce the notaton p nm = p n p m, r nm = r n r m an ˆr = r/ r. Eq. (6) has a smple nterpretaton. The factor (p nm ˆr nm )/M gves the flux of ncomng partcles. The Θ- an δ-functons specfy the contons for a collson to take place. A collson between partcles n an m happens only f the two partcles are approachng each other, whch s taken nto account by Θ( p nm ˆr nm ). At the nstant of a collson the stance between the two partcles has to vansh, as expresse by δ( r nm a). Fnally (b nm 1) generates the change of lnear an angular momenta [13]. The ensemble average of a ynamcal varable s efne by < A > t = Γρ(0)A(t) = Γρ(t)A(0) = (r p ω ) <jθ( r j a)ρ(t)a(0). (7) Here ρ(t) = exp( L + t)ρ(0) s the N-partcle strbuton functon, whose tme evoluton s governe by the ajont L + of the tme evoluton operator L +. The quanttes of nterest are the translatonal an rotatonal energes per partcle as well as the total knetc energy E = E tr +E rot. As a frst step we compute the ntal ecay rates E tr = 1 N E rot = 1 N M 2 v2 (8) I 2 ω2 (9) t < E tr > t=0 = < L + E tr > t=0 = ν tr (10) t < E rot > t=0 = < L + E rot > t=0 = ν rot (11) assumng that the system has been prepare n a thermal equlbrum state ρ(0) exp( E/T) whch s characterze by the temperature 3T/2 =< E > t=0. The change of translatonal an rotatonal energy are gven by ν tr = 2πa2 T 3/2 ( n 0 g(a) (1 ǫ 2 )+ k ) (πm) 1/2 1+k (1 β2 ) (12) ν rot = 2πa2 T 3/2 ( ) n 0 g(a) 1 (πm) 1/2 1+k (1 β2 ). (13) Here g(a) enotes the par correlaton at contact an n 0 = N/V. Rotatonal energy s conserve n two cases, for ether perfectly smooth spheres (β = 1) or perfectly rough spheres (β = +1). Translatonal energy s only conserve f n aton ǫ = 1. For all other values of the parameters ǫ an β the translatonal an rotatonal energy ecrease lnearly wth tme but wth fferent rates. Ths mples that after a small tme nterval t equpartton among rotatonal an translatonal egrees of freeom no longer hols 2

3 < E tr > t = < E tr > 0 ν tr t < E rot > t = < E rot > 0 ν rot t. (14) Hence for collsons wth frcton (β 1) the homogeneous coolng state s not consstent wth the tme evoluton of translatonal an rotatonal energy separately. To generalse the concept of a homogeneous coolng state to collsons wth frcton we ntrouce two temperatures T tr (t) = 2/3 < E tr > t an T rot (t) = 2/3 < E rot > t. We keep the assumpton of spatal homogenety an assume that both lnear an angular momenta are normally strbute wth n general fferent tme epenent wths or temperatures ρ(t) exp 1 2 N (Mv 2 (t)+iω 2 /T rot(t)). (15) =1 The above strbuton allows for a calculaton of the acay rates of T tr (t) an T rot (t) for arbtrary tmes, resultng n two couple fferental equatons: t T tr = (1 ǫ 2 ) γ 4 T3/2 tr γη t (1 η t )T 3/2 tr + γ k η2 tt 1/2 tr T rot (16) t T rot = γ k η2 tt 3/2 tr γ η ( t 1 η ) t T rot T 1/2 tr (17) k k wth γ = (16/3)a 2 n 0 g(a) π/m. As we have seen above the energes ecrease lnearly for short tmes. For large tmes both, translatonal an rotatonal energy fall off algebracally lke t 2 wth however fferent ampltues. The ampltues can be etermne analytcally, e.g. by solvng the fferental equaton for T tr T rot = T tr ( (1 ǫ 2 )/4+η t (1 η t ) ) +T rot η 2 t/k T tr η 2 t /k T rot(1 η t /k)η t /k (18) Ths equaton s solve by a constant rato T tr /T rot = c 1 /c (c 1 /c 2 ) 2 wth c 1,c 2 gven n terms of ǫ,β an k by c 1 = (1 ǫ 2 )/4+(1 β 2 )(k 1)/(4k+4) an c 2 = 2k(1+β) 2 /(2+2k) 2 n agreement wth ref. [5]. The crossover tme t 0 between the lnear regme an the algebrac ecay s etermne by γ an T tr (0) accorng to t 1 0 γ T tr (0). The full soluton [15] has been obtane by numercal ntegraton of eqs.(16,17) an s shown n Fg. (1) for k = 0.4, γ = 1, ǫ = 0.4, β = 0.6 an T tr (0) = T rot (0) = 20. The asymptotcs s ncate by a ashe otte lne an the nset s a blow up of the short tme behavour T tr (0) T rot (0) Asymptotc FIG. 1. Tme ecay of the translatonal an rotatonal energy. Parameters were chosen accorng to: k = 0.4, ǫ = 0.4, β = 0.6 an T tr(0) = T rot(0) = 20. We have ntrouce a mensonless tme argument τ = tγ T tr(0). The ashe otte lne ncates the asymptotcs τ 3

4 The evaton of the rato of rotatonal an translatonal energy from 1,.e. 1 T rot, s shown n Fg. (2) for two fferent values of β an ǫ ε = 0.6; β = 0.4 ε = 0.4 β = T rot FIG. 2. Rato 1 T rot as a functon of tme wth k = 0.4, γ = 1, T tr(0) = T rot(0) = 100, ǫ = 0.6, β = 0.4 an ǫ = 0.4, β = 0.6. Unt of tme was chosen as n Fg. (1). Coolng of a granular gas has been nvestgate by several groups[16,17] moellng collsons wth a normal coeffcent of resttuton. The t 2 behavour of the total energy has been confrme for a we range of parameters such that the system s stable to shear fluctuatons [17] an nelastc collapse. Lung [18] has consere a more etale moel of nelastc collsons nclung tangental resttuton an Coulomb s law of frcton. It s straghtforwar to generalse the moel to nclue Coulomb frcton, one just has to mofy the upatng rules (eq. (3)) accorngly. Drvng the system wth a vbratng wall s an open problem so far. One may also generalse the approach to nonsphercal objects, e.g. har ros or neeles [19], for whch elastc collsons have been analyse by Frenkel an Magure [20]. We wsh to thank Tmo Aspelmeer an Kurt Broerx for useful scussons an one of us (M.H.) the Lan Neersachsen for grantng a scholarshp. Note ae n proof: After ths paper was complete we learnt about relate work by S. McNamara an S. Lung who also stuy the rato of translatonal an rotatonal energy n the HCS. τ [1] J. T. Jenkns an S. B. Savage, J. Flu Mech. 130, 187 (1983). [2] C. K. K. Lun, S. B. Savage, D. J. Jeffrey an N. Chepurny, J. Flu Mech. 140, 223 (1984). [3] J. T. Jenkns an M. W. Rchman, Phys. Flus 28, 3485 (1985). [4] C. K. K. Lun, J. Flu Mech.233, 539 (1991). [5] A. Golshten an M. Shapro, J. Flu Mech. 282, 75 (1995). [6] C. Cercgnan n Fsca matematca et ngegnera elle strutture: rapport e compatblta, e. G. Ferrarese, Ptagora Etrce, Bologna 1995, p.119. [7] S. McNamara an W. R. Young, Phys. Flus A 5, 34 (1993). [8] S. E. Espov an T. Poeschel, J. Stat. Phys. 86, 1385 (1997). [9] T. Aspelmeer, G. Gese an A. Zppelus, con-mat/ [10] M. H. Ernst, J. R. Dorfmann, W. R. Hoegy an J. M. J. van Leeuwen, Physca 45, 127 (1969). [11] P. M. Furtao, G. F. Mazenko an S.Yp, Phys. Rev. A 12, 1653 (1975). [12] T. P. C. van Noje an M. H. Ernst, con-mat/ (1997). [13] In some calculatons t s helpful to know the evoluton backwars n tme. In that case the momenta after collson are gven an one wants to know the momenta before collson. These are obtane by replacng the parameters n eq. (3) accorng to β 1 an ǫ 1. Wth ths replacement the momenta before collson are gven by the rght han se of β ǫ eq.(3) n terms of the momenta after collson. These collson rules can then be use to construct an operator L, whch escrbes the evoluton backwars n tme, but n general L s not the ajont of L + (see also [12]). [14] T. P. C. van Noje, M. H. Ernst, R. Brto an J. A. G. Orza, Phys. Rev. Lett. 79, 411 (1997). [15] Wth the substtuton T tr = z 2 the two couple fferental equatons can be reuce to a sngle secon orer equaton z + azz +bz 3 = 0, whch s scusse n E. Kamke Dfferentalglechungen, Akaemsche Verlagsgesellschaft, Lepzg 1944, p

5 [16] S. McNamara an W. R. Young, Phys. Rev. E 53, 5089 (1996). [17] I. Golhrsch an G. Zanett, Phys. Rev. Lett. 70, 1619 (1993). [18] S. Lung, Phys. Rev. E 52, 4442 (1995). [19] M. Huthmann an A. Zppelus, n preparaton. [20] D. Frenkel an J. F. Magure, Mol. Phys. 49, 503 (1983). 5