Dynamical approach to weakly dissipative granular collisions


 Belinda Reeves
 1 years ago
 Views:
Transcription
1 Dyamical approach to weakly dissipative graular collisios Italo Ivo Lima Dias Pito ad Alexadre Rosas Departameto de Física, CCEN, Uiversidade Federal da Paraíba, Caixa Postal 58, , João Pessoa, Brazil Katja Lideberg Departmet of Chemistry ad Biochemistry ad BioBircuits Istitute, Uiversity of Califoria Sa Diego, La Jolla, Califoria , USA (Dated: Jue 3, 25) Graular systems preset surprisigly complicated dyamics. I particular, oliear iteractios ad eergy dissipatio play importat roles i these dyamics. Usually (but admittedly ot always), costat coefficiets of restitutio are itroduced pheomeologically to accout for eergy dissipatio whe grais collide. The collisios are assumed to be istataeous ad to coserve mometum. Here, we itroduce the dissipatio through a viscous (velocity depedet) term i the equatios of motio for two collidig grais. Usig a first order approximatio, we solve the equatios of motio i the low viscosity regime. This approach allows us to calculate the collisio time, the fial velocity of each grai, ad a coefficiet of restitutio that depeds o the relative velocity of the grais. We compare our aalytic results with those obtaied by umerical itegratio of the equatios of motio ad with exact oes obtaied by other methods for some geometries. I. INTRODUCTION The characterizatio of graular matter is extremely broad, ad icludes essetially ay coglomeratio of discrete macroscopic particles. These ca be as small as grais of sad ad as large as asteroids, they ca be i a codesed or gaslike phase. The codesed phases may exhibit characteristics of solids or fluids or gases or various combiatios thereof. Graular matter is importat i more idustrial applicatios tha ca be listed here, ad exhibits a huge variety of iterestig behaviors that have provided food for thought over ceturies of time. Behaviors such as the socalled jammig trasitio ad the formatio of patters are frequet subjects of curret research, as is the propagatio of eergy i graular materials. A feature commo to graular matter is the fact that eergy is lost every time grais collide. Ideed, the grais have to be very hard ad difficult to compress for a collisio ot to lead to a loss of eergy. Yet it is usually the case that mometum is coserved i these ielastic collisios. The coservative limit, where oly elastic collisios are ivolved, is famously illustrated by Newto s cradle, cosistig of a row of very hard balls that just touch, each hagig o a strig of the same legth attached to a commo support. Whe the ball at oe ed is picked up ad released so that it collides with the ext ball, the eergy passes dow the row util the last ball flies up to the same height as the first ball before it was released (the other balls remai at rest). The last ball the flies back, the eergy is trasferred across the row agai, ad the first ball flies up to the same height [ 3]. This cotiues, although ot forever because of course some small amout of eergy must be lost at each collisio evet. The prototypical pheomeological descriptio of eergy loss ivolves the coefficiet of restitutio ε i the equatio that describes a collisio betwee two grais, v f2 v f = ε (v i v i2 ). () Here the v s represet the velocities, the subscripts i ad f stad for iitial (before collisio) ad fial (after collisio) ad the umbers label particles ad 2. This descriptio leads to a eergy loss at each collisio of ε 2 of the kietic eergy of the ceter of mass before the collisio. For a successfully built Newto s cradle, ε is exceedigly small. The coefficiet of restitutio is most ofte treated as a parameter idepedet of the velocities. Ad yet, it is broadly recogized that this ca ot be totally correct because it leads to problematic features i the asymptotic behavior such as the socalled ielastic collapse i a graular gas because there may be a ifiite umber of collisios i a fiite time [4, 5]. Ideed, whe oe cosiders realistic iteractio models, it is i fact uiversally the case that iteractios of ay two compressible grais are oliear. For istace, the iteractios betwee two spherical objects obey Hertz s law, where the repulsio is proportioal to the compressio to the power 3/2 rather tha the more familiar Hook s law where the repulsio is simply proportioal to the compressio. The cosequece of this oliearity is that the duratio of a collisio depeds o the iitial velocities of the particles before the collisio. Therefore, whe a dissipative collisio occurs, the mechaism resposible for the dissipatio of eergy acts for differet legths of time depedig o the iitial velocities, leadig to distict eergy losses, ad cosequetly to a velocitydepedet coefficiet of restitutio. It is iterestig to ote that Hertz s law is frequetly used together with a velocityidepedet coefficiet of restitutio. The history of the aalysis of the effect of velocitydepedet frictioal forces o the coefficiet of restitutio is log ad varied. Here we oly summarize some of
2 2 its saliet poits. Hertz s law assumes that there are o attractive surface forces i lightly loaded spherical graules. Perhaps the earliest work to recogize that such forces lead to a fiite cotact area betwee surfaces uder zero load (adhesio), ad that this i tur modifies the exteral force required to separate two bodies of give surface eergy ad geometry, is that of Johso et al. (commoly kow as JKR theory) [6]. This i tur modifies the velocity depedece of the coefficiet of restitutio, as aalyzed i detail by Brilliatov et al. [7]. While this correctio may be small, eve egligible, betwee graules of high elastic modulus such as metals or glass, it is cosiderable i soft graules such as rubber ad agrees well with experimets i these cases. Much of the literature o this topic is based o spherical graules, but other shapes are also discussed i this cotext. For istace, the work of Walto ad Brau o frictioal disks is oteworthy [8], as is the work of Herbold ad Nestereko [9]. I [] the coefficiet of restitutio for spherical graules obeyig Hertz s law for the elastic portio of the iteractios plus a velocitydepedet frictioal force was obtaied exactly as a ifiite series. The series ca ot be summed aalytically ad coverges very slowly, ad trucatio of the series leads to uphysical divergeces. A compact Padé approximatio to the series makes it possible to perform much more efficiet evetdrive molecular dyamics simulatios as well as direct Mote Carlo methods tha was possible with earlier methods, ad provides excellet results whe compared to those obtaied from umerical itegratios of the equatios of motio []. Determiatio of the velocity depedece of the restitutio coefficiet has recetly become icreasigly detailed ad addresses a more varied rage of particle sizes ad compositios. For istace, experimetal determiatio [2] ad molecular dyamics calculatios [3] of the velocity depedece of the coefficiet of restitutio of argo aoparticles have show that these aoparticles are hard ad highly elastic at collisio velocities smaller tha the sizedepedet yield velocity, while they progressively softe as the collisio velocity icreases beyod the yield velocity. I this cotributio, our goal is ot to arrive at a more efficiet approach for umerical simulatios, but rather to preset a perturbative but very broadly applicable methodology to aalytically aalyze the cosequeces of a geeral viscoelastic model of dissipatio o the outcome of a collisio i the low dissipatio regime. Our study provides a alterative approach to the problem. We arrive at approximate but accurate aalytic expressios for the iitial velocity depedece ot oly for the coefficiet of restitutio but also for the duratio of a collisio ad for the velocities of the graules at the ed of the collisio for large rages of parameter values. Our results have iterestig implicatios for the uderstadig of a full rage of realistic collisios, which should be useful for the study of graular gases ad for pulse propagatio i graular chais. Examples of graular systems with low itergraular frictio iclude sad, glass beads, ad steel beads. We accomplish three objectives, all i the limit of low but ozero dissipative forces. The first, preseted i Sec. II, is to discuss the equatios of motio that describe the collisio of two ot ecessarily spherical grais. The sceario is this: the two grais are iitially just touchig. Grai has velocity v ad grai 2 has velocity v 2. How these grais came to have these velocities is ot importat. If our two grais are elemets i a graular gas, they may, for example, have arrived at this poit due to previous collisios with other grais. Or, oe might be preparig a experimet with just two grais, where oe grai is give a kick of some kid at time t = to cause it to start movig with velocity v while the other is iitially at rest (v 2 = ). The equatios of motio the determie the further evolutio of the two graules. I particular, we are able to calculate the relative velocities of the graules at the ed of a collisio as a fuctio of the iitial relative velocities. Our secod objective, detailed i Sec. III, is to use this result to calculate the duratio of a collisio as a fuctio of the iitial velocities. This makes explicit our assertio that collisios are i geeral ot istataeous ad that their duratio i fact depeds o the iitial velocities of the collidig graules. As poited out above, this i tur leads to a velocitydepedet coefficiet of restitutio. Our third objective is to calculate the coefficiet of restitutio ad the velocities of the collidig graules at the ed of the collisio. This is preseted i Sec. IV. I the absece of dissipative forces the coefficiet of restitutio is equal to uity. We are able to explicitly calculate the lowest order correctios to this ad thus to obtai a explicit form for the depedece of the coefficiet of restitutio o the iitial velocity differece of the graules. Fially, Sec. V cotais a summary of the paper ad commets o the possible geeralizatios of this study. II. THE MODEL Viscoelastic forces for the collisio of two spheres iclude two terms. The first is due to the elastic repulsio betwee the two particles ad has its origis i Hertz s theory [4, 5]. The secod term stads for the viscous dissipatio via a dashpot [6, 7]. Hece, the cotact force ca be writte as F = r (x x 2 ) γ (x x 2 ) α (ẋ ẋ 2 ). (2) Here x is the displacemet of a particle from its iitial positio at the begiig of a collisio. A dot deotes a derivative with respect to time, ad the subscripts o x label the two particles. The coefficiet r is a costat depedet o Youg s modulus ad Poisso s ratio. For istace, for collidig spheres of radii R ad R 2 this costat is give by r = 2 5D 2 2R R 2 R + R 2, (3)
3 3 where D 2 = 3 4 ( σ 2 + ) σ2 2, (4) E E 2 σ i is Poisso s ratio ad E i the Youg s modulus of sphere i. Returig to Eq. (2), γ is the coefficiet of viscosity, α is a costat that defies the specific viscoelastic model, ad depeds o the topology of the cotact betwee the particles. For spheres is equal to 5/2 [4, 5], ad α to /2 [8, 9], but we leave them as ad α for the sake of geerality. We say that particle is to the left of particle 2, ad that our system of coordiates icreases from left to right. Whe we work with two equal spheres so that R = R 2, σ = σ 2, ad E = E 2, we drop the subscript o D, that is, we set D D. The equatios of motio for two particles of mass m ad m 2 durig a collisio are m ẍ = r (x x 2 ) γ (x x 2 ) α (ẋ ẋ 2 ), m 2 ẍ 2 = r (x x 2 ) + γ (x x 2 ) α (ẋ ẋ 2 ). (5) From Eq. (5), coservatio of mometum immediately follows, m ẍ + m 2 ẍ 2 =, (6) so that m ẋ + m 2 ẋ 2 = cost. Equatio (5) also leads to a ucoupled equatio for the differece variable z = x x 2, z = r µ z γ µ zα ż, (7) where µ is the reduced mass µ = m + m 2. For the latter equatio, the iitial coditios are z() = because we deal with cofiguratios where the graules are iitially just touchig each other, ad ż() = v i v i2. Here v i ad v i2 are the iitial velocities of the two collidig graules. A aalytic solutio z(t) of Eq. (7) for arbitrary ad α seems ot to be available. However, we have bee able to obtai the velocities at the ed of the collisio as a fuctio of the iitial velocities i the low viscosity limit. We rewrite Eq. (7) as a first order differetial equatio of the velocity as a fuctio of the positio. Defiig v = ż, ad otig that we rewrite Eq. (7) as z = dv dt = dv dz dz dt = v dv dz, (8) v dv dz = r µ z γ µ zα v. (9) I the absece of dissipatio (γ = ), this equatio admits two solutios for v(z), v ± (z) = ± 2 2rz µ. () Obviously, the positive sig should be cosidered durig compressio, ad the egative oe durig decompressio. Our approximatio i the low dissipatio regime starts by writig the velocity as a perturbatio o the odissipative solutio. Cosequetly, durig compressio we have v(z) = v + (z) + γv comp (z), () where v comp (z) is a fuctio to be determied. Substitutig the trial solutio Eq. () i Eq. (9), ad collectig the terms of order γ, we have r µ z v comp + ( ) ( ) 2rz z µ 2 α µ + v comp =, (2) where a prime deotes a derivative with respect to z. Furthermore, the coditio v comp () = is ecessary to satisfy the iitial coditios. The solutio of Eq. (2) with the iitial coditios v comp () =, whe added to v + (z), the gives us the compressio velocity to first order i the dissipatio, v(z) = 2 2rz µ γ z+α [ 2( + α) + 2 F (, 2 + +α ; + +α ( + α)( α)µ ] ; 2rz µv 2 ), (3) where 2 F is a hypergeometric fuctio [2]. Next we move o to the decompressio, which starts with v = ad z = z max. To calculate the maximum compressio z max, we agai look for a first order correctio i γ, ( ) 2 / µv z max = ( γz c ), (4) 2r where z c is a costat to be determied ad ( µv 2 /2r ) / is the maximum compressio i the absece of dissipatio, obtaied by settig v ± (z) = i Eq. (). The relative velocity must vaish for the maximum compressio. This is the poit at which the collidig grais stop movig ad begi to separate. Therefore, we evaluate Eq. (3) at z = z max as give i Eq. (4),
4 4 expad to first order i γ, ad set the left had side equal to zero. Solvig the resultig equatio for z c ad isertig the result i Eq. (4), we fid for the maximum compressio ( 2 µv z max = 2r ) γ ++α 2 ( π µ ) +α r 2 +2α 2 µγ ( α Γ ( ) +α ). (5) While it might be temptig to assume that the collisio eds whe z =, that is, whe the particles lose cotact, this would eglect the fact that the particles have ot recovered their origial shape at this poit []. The effect of eglectig this cotributio leads to serious spurious attractive forces. The collisio actually eds whe the force is zero, z =, at which poit the iitial shape is recovered. Settig z =, z = z f, ad ż = v f i Eq. (7), where the subscript f stads for fial, we have ( γvf z f = r ) α. (6) I what follows we restrict to values > α +, so that z f is small for low viscosity. We write the perturbative solutio for the decompressio velocity i the low viscosity case as v(z) = v (z) + γv decomp (z). (7) Agai substitutig ito Eq. (9) ad expadig the latter to first order i γ, we fid v decomp (z) = C µv2 2rz 4rz ++α z+α ( α) µ (2rz µ2 ) [2 + 2α + 2 F (, 2 + +α ; + +α ( + α) ( α) µ )] ; 2rz µv 2. (8) Here C is a costat to be determied by the cotiuity of the solutios (8) ad (3) at z max. Rememberig that v(z max ) =, where z max is give by Eq. (5), ad expadig up to first order i γ, C is foud to be C = 2 π(µ) ( 2 + ) r v +2α+2 Γ ( ) ). Γ ( (9) To use this result to calculate the leadig cotributio to the fial relative velocity we write v f = v(z f ) = v() + v ()z f, (2) sice z f is small. A prime deotes a derivative with respect to the argumet z. Sice the force is zero at the ed of the collisio, we have z = = dż/dt = dż/dz dz/dt = v (z f )v(z f ) = [v() + v ()z f ][v () + v ()z f ] = v()v () + [v()v () + (v ()) 2 ] + O(zf 2 ). (2) This equatio has two solutios for v (): either v = v()/z f or v =. The first leads to a fial relative velocity that vaishes, which is ot physical i our perturbative approach where we expect the fial relative velocity to go to whe γ. We must thus choose v () = (or, more rigorously, v () = O(zf 2 )). Therefore, v f = v() + O(z 3 f ) = v() + O(γ 3/( α) ). (22) We choose values of ad α such that < α < 3. For spheres with = 5/2 ad α = /2 this coditio is clearly satisfied. We later also cosider = 5/2, α = (used by others), for which the coditio is also satisfied. The correctio i the fial velocity ca the be eglected i our aalysis to first order i γ. The first term, v(), is obtaied by substitutig Eq. (9) ito Eq. (8), usig this result together with Eq. (7), ad takig the limit z. We the obtai v f = + γ III. ( µv 2 2r ) +α πγ( + +α ) µγ( α ). (23) COLLISION TIME Our ext objective is to calculate the collisio time, which we do i two parts. First we calculate the com
5 5 pressio time of the collisio as the two grais compress oe aother ad the, as the grais move apart, the attedat decompressio time. The collisio time is the the sum of the two. First we calculate the compressio time as T compressio = zmax dz (24) v(z) where, to order γ, v(z) is obtaied from Eq. (3) as v(z) = 2 2rz µ z [2( +α + α) + 2 F (, 2 + +α + γ ] ; 2rz µv 2 ) ; + +α ( + α) ( α) µ(v 2 2rz µ ). (25) Substitutig this ito Eq. (24) ad itegratig leads to a cotributio of order γ, ad a cacellatio of two terms of order γ /2. Cosequetly, the compressio time does ot show ay γ depedece up to first order: T compressio = 2 + πr µ + v 2 Γ( ) Γ( 2 + ) +O(γ 3/2 ). (26) Next we move o to decompressio which starts with v = ad z = z max (maximum compressio) ad eds whe z = (force is zero), at which poit z = z f as give i Eq. (6). Hece, the decompressio time is T decompressio = = zf z max z max dz v(z) dz v(z) + zf dz v(z). (27) We first deal with the first itegral. Sice the zeroth order term of the velocity durig decompressio, Eq. (7), is the egative of the compressio velocity i Eq. (3), ad the limits of itegratio of the compressio ad the first term i the decompressio times are switched, the latter time is the same as the compressio time up to terms of order γ 3/2. This still leaves the secod term i the decompressio time give i Eq. (27). Sice z max is small, zf dz v(z) z f v(z f ) z f = ( γv ) α. (28) r For α >, this correctio is ideed small ad our perturbative approach is valid. I particular, for spherical grais = 5/2, ad the particular values of α = /2 (commoly used i the literature) ad α = (used by others, see below), the value of the expoet of γ is 2/3 ad respectively. This i tur implies that this term is a small perturbatio. We have thus established that for the parameters used herei, the decompressio time is equal to the compressio time up to order γ 2/3 or order γ depedig o the value of α. Hece the total collisio time to this order is twice the compressio time plus a additioal cotributio to the decompressio time, T = 2 + πr µ + v 2 Γ( ) Γ( 2 + ) ( γv ) α, (29) r ad we have arrived at our secod objective. We have ot foud ay other calculatio of the velocitydepedece of the collisio time i the literature. I a recet paper, the merits ad problems of differet choices of the parameter α, ad eve a geeralizatio of the above model, were discussed [7]. Here, for = 5/2 we cosider two choices of this parameter that have bee commoly used i the literature to test Eq. (2). The simplest case, α =, was proposed i [2] ad further developed i [22], but i our subsequet results below we see that this choice leads to problematic outcomes. The case α = /2, proposed idepedetly i [8, 23], is more widespread i the graular gas commuity. Yet aother combiatio of expoets was foud experimetally for chais of origs i[9] ad aalyzed theoretically i[24]. I this model dissipatio was ot icluded, but the elastic force is a double powerlaw rather tha a sigle oe, with expoets 5/2 ad 7. I the latter work it was show that depedig o the characteristics of the origs ad the experimetal setup, either the oe or the other cotributio ca be domiat. I Fig. we show the duratio of the collisio for two equal spherical grais ( = 5/2) as a fuctio of the relative velocity at the begiig of the collisio for both values of α. As ca be see i the figure, for small γ the data is very well predicted by our approximatio, idepedetly of the expoet α (the approximatio oly starts to deviate from the theoretical predictio for γ =., represeted by the plus sigs). This αidepedece of the duratio of the collisio is oe of the strikig predictios of our theory. Aother iterestig characteristic of our solutio is the powerlaw depedece of T o the iitial relative velocity, as evideced i the iset of Fig.. I Fig. 2 we agai show T as a fuctio of, but this time for spheres of differet sizes. Fixig the radius of graule ad varyig the radius of graule 2 (assumig
6 6 T log T log IV. FINAL VELOCITIES AND COEFFICIENT OF RESTITUTION We ext tur our attetio to the fial velocities. From the coservatio of mometum, we kow that the total mometum at the begiig ad ed of the collisio must be the same, m v i + m 2 v i2 = m v f + m 2 v f2, (3) where the subscripts i ad f oce agai label the iitial ad fial velocities of the grais. O the other had, we also kow from the defiitio of z(t) that ż(t ) = v f = v f v f2. (3) FIG.. Collisio duratio for several values of γ (from top to bottom,.,.,. ad.) ad two differet values of α ( for lower group ad /2 for higher group), obtaied via umerical itegratio of the equatios of motio (symbols) ad from our theory (lie). The lie represet the theoretical predictio of Eq. (29). The iset shows the same data o a loglog scale. I this figure, the parameters are as follows: m = m 2 =, R = R 2 =, ad 2/(5D) =. As predicted, for these values of γ ad α the collisio duratios are essetially idepedet of these parameters. T FIG. 2. Duratio of collisios for several values of R 2 (from bottom to top, R 2 varies from. to.9 i steps of.) obtaied via umerical itegratio of the equatios of motio. The lies represet the theoretical predictio Eq. (29). I this figure the parameters are as follows: m =, R =, α = /2, γ =., ad 2/(5D 2) =. The desities of the two collidig spheres are equal. that they have the same desity), we ca see that the collisio takes loger for larger values of R 2. Solvig the set of the two equatios give above for v f ad v f2, ad usig Eq. (7), we fid v f = m m 2 m + m 2 v + 2m 2 m + m 2 v m 2 m + m 2 ( πγ ) ( r µ rγ ( v f2 = m m 2 m + m 2 v 2 + 2m m + m 2 v 2 m m + m 2 ( πγ ) ( r µ rγ ( ) ) (v i v i2 ) 2() γ, ) (32) ) (v i v i2 ) 2() γ. As expected, the result for a elastic collisio is recovered whe γ =. Further, the ifluece of the dissipatio is greater o the lighter particle, ad the ifluece of the iitial coditio o the chage i the fial velocities due to dissipatio depeds oly o the relative velocity. We coclude this sectio by usig the above results i Eq. () to calculate the coefficiet of restitutio ad thus completig our third ad pricipal objective: ε = γ 2 ( πγ ) ( r µ ) ) (v i v i2 ) +2α+2. rγ ( (33) I Fig. 3 we show the coefficiet of restitutio for several values of α. The agreemet is equally good for all of them. A importat characteristic of ε is that its qualitative depedece o the iitial relative velocity is drastically differet for α larger tha or smaller tha ( 2)/2 (i the case of spheres, this value is /4): α smaller tha this value leads to the uphysical situatio of egative coefficiets of restitutio for very small relative velocities. For larger α, the collisio approaches the elastic case for small relative velocities. Perhaps most importatly, the power ( + 2α + 2)/ of the iitial relative velocities of the two graules is equal to /5 with the physically justified values = 5/2 ad α = /2 for spheres, see e.g. Eq. (2) i [25]. I additio to recoverig this expoet, as i Eq. (4) of that referece, we also obtai
7 ε ε FIG. 3. Coefficiet of restitutio as a fuctio of the iitial relative velocity. Each curve correspods to a differet value of α (from bottom to top o the left side, α varies from to.9 i steps of.). The other parameters are: γ =., m = m 2 =, R = R 2 =, ad 2/(5D) =. The lies correspod to the theoretical predictio of Eq. (33) (a)α = the liear depedece of the correctio of the coefficiet of restitutio o the coefficiet of viscosity γ. I Fig. 4, we show the coefficiet of restitutio as a fuctio of for particles of differet sizes (same desity) for α = ad α = /2. It is evidet from the figure that ε icreases with the radius i both cases. However, the qualitative behavior is idepedet of the sizes of the grais. V. CONCLUSIONS We have succeeded i calculatig quatities that characterize the collisio of two graules that lead to the loss of eergy (but ot mometum) to the eviromet via viscous dissipatio. We started with a equatio of motio (Newto s Law) cotaiig a kietic eergy cotributio, a force due to the elastic repulsio betwee the two graules, ad a dashpot viscous dissipatio term. I additio to parameters related to the shape ad size of the graules, the model cotais two importat parameters: a coefficiet of viscosity γ, ad a costat α that defies the specific viscoelastic model, cf. Eq (2). Our calculatios are perturbative i the coefficiet of viscosity, that is, we preset lowest order correctios to elastic (eergycoservig) collisios. A collisio begis with the two graules just touchig heado toward each other with a relative velocity. This velocity ad cofiguratio defie the collisio stregth. The collisio begis at this iitial momet with compressio of the graules util their relative velocity is zero (at which poit the compressio is a maximum). Decompressio the follows, util the force betwee the graules vaishes, at which poit the collisio eds. ε (b)α = /2 FIG. 4. Coefficiet of restitutio as a fuctio of the iitial relative velocity. Each curve correspods to a differet value of the radius of graule 2 (from bottom to top, R 2 varies from.5 to.4 i steps of.). The other parameters are: γ =., m =, R =, ad 2/(5D 2) =. The desities of the two collidig spheres are equal. The lies correspod to the theoretical predictio of Eq. (33). Itegratio of the equatios of motio leads to aalytic results for several importat quatities usually specified simply as pheomeological parameters. The first is the relative velocity of the graules durig compressio ad durig decompressio. We calculate the fial relative velocity as a fuctio of the separatio of the ceters of the graules ad fid the depedece o iitial relative velocity ad o the parameters γ ad α, cf. Eq. (23); if the collisio were elastic, the fial ad iitial relative velocities would of course just be the egatives of oe aother. A secod set of useful results are the fial velocities of each grai, for which we obtai explicit expressios as a fuctio of the parameters ad of the iitial velocities of each grai, cf. Eq. (32). These are importat for simu
8 8 latios of graular gases. The third quatity we calculate is the duratio of the collisio, cf. Eq. (29). I most pheomeologies, collisios are assumed to be istataeous. Collisios are of course ot istataeous. Fially, we fid a aalytic expressio for the coefficiet of restitutio ε defied i Eq. (). This coefficiet, usually chose pheomeologically, recogizes the ielasticity of graular collisios. We have foud the depedece of this coefficiet o particle shape (via the expoet i the force that determies the topology of the cotact betwee the graules), the coefficiet depedet o Youg s modulus ad Poisso s ratio, ad, most importatly, o iitial relative velocity ad o the parameters γ (to lowest order) ad α that defie the viscoelastic model, cf. Eq. (33). If γ = the coefficiet of restitutio is uity, that is, there is o eergy loss i the collisio. Similarly, if the iitial velocities of the two graules are equal, the coefficiet is trivially uity agai. The depedeces o these quatities are otrivial ad, we submit, essetially impossible to arrive at pheomeologically. This the provides a physical basis for the usual pheomeological choice ε <. While our results are perturbative i γ ad thus ot as geeral for elastic spheres ( = 5/2, α = /2) as is the ifiite series developed i [], our model allows differet values of ad α ad yields relatively simple explicit results for the iitial velocity depedece of the fial velocities of the collidig graules, of the duratio of a collisio, ad of the restitutio coefficiet. I this paper we have oly dealt with two collidig graules, takig ito accout the eergy loss due to a explicit viscoelastic force i the equatios of motio. This reders our results immediately applicable to graular gases where at low desities biary collisios are the most commo iteractios. The geeralizatio to a graular chai or to eve higher dimesioal graular arrays is ot trivial, cf. Ref. [25], but is ow made cosiderably easier by the fact that we have explicitly foud the pricipal igrediets of the problem. There is evertheless a great deal of work to be doe, especially toward higher dimesioal geeralizatios. Ackowledgmets I. L. D. Pito ad A. Rosas akowledge the CNPq ad BioaotecCAPES for fiacial support. [] S. Hutzler, G. Delaey, D. Weaire ad F. Macleod, Am. J. Phys. 72, 58 (24). [2] P. Glediig, Phys. Rev. E 84, 672 (2). [3] K. Sekimoto, Phys. Rev. Lett. 4, 2432 (2). [4] Sea McNamara, W. R. Youg. Ielastic collapse ad clumpig i a oedimesioal graular medium. Phys. Fluids A 4 (3), 496 (992). [5] B. Beru ad R. Mazighi J. Phys. A:Math. Ge. 23, 5745 (99). [6] K. L. Johso, K. Kedall ad A. D. Roberts, Proc. R. Soc. Lodo A 324, 3 (97). [7] N. V. Brilliatov, N. Albers, F. Spah ad T. Pöschel, Phys. Rev. E 76, 532 (27). [8] O. R. Walto ad R. L. Brau, J. Rheol. 3, 949 (986). [9] E. B. Herbold ad V. F. Nestereko, Appl. Phys. Lett. 9, 2692 (27). [] T. Schwager ad T. Pöschel, Phys. Rev. E 78, 534 (28). [] P. Mueller ad T. Pöschel, Phys. Rev. E 84, 232 (2). [2] A. I. Ayesh, S. A. Brow, A. Awasthi, S. C. Hedy, P. Y. Covers, ad K. Nichol, Phys. Rev. B 8, (2). [3] Y. Takato, S. Se, ad J. B. Lechma, Phys. Rev. E 89, 3338 (24). [4] H. Hertz ad J. Reie Agew. Math. 92, 56 (88). [5] L. D. Ladau ad E. M. Lifshitz, Theory of Elasticity, AddisoWesley, MA, (959). [6] S. Ji ad H. H. She, J. Eg. Mech. 32, 252 (26). [7] E. Alizadeh, F. Bertrad ad J. Chaouki, Powder Techology 237, 22 (23). [8] N. Brilliatov, F. Spah, J.M. Hertsch ad T. Pöschel, Phys. Rev. E 53,5382 (996). [9] W. A. M. Morgado ad I. Oppeheim, Phys. Rev. E 55, 94 (997). [2] George B. Arfke, Has J. Weber, Frak Harris, Mathematical Methods for Physicists (Academic Press, 5th editio, 2). [2] J. Lee ad H. J. Herrma, J. Phys. A: Math. Ge. 26, 373 (993). [22] A. Rosas, A. H. Romero, V. F. Nestereko, ad K. Lideberg, Phys. Rev. Lett. 98, 643 (27). [23] G. Kuwabara ad K. Koo, Japaese J. Appl. Phys. 26, 23 (987). [24] I. L. Dias Pito, A. Rosas, A. H. Romero, ad K. Lideberg, Phys. Rev. E 82, 338 (2). [25] T. Pöschel ad N. V. Brilliatov, Phys. Rev. E 63, 255 (2).
Kinetics of Complex Reactions
Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet
More informationx a x a Lecture 2 Series (See Chapter 1 in Boas)
Lecture Series (See Chapter i Boas) A basic ad very powerful (if pedestria, recall we are lazy AD smart) way to solve ay differetial (or itegral) equatio is via a series expasio of the correspodig solutio
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rateofchage quatities, such as velocity ad acceleratio). Numerical
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationSECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES
SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES Read Sectio 1.5 (pages 5 9) Overview I Sectio 1.5 we lear to work with summatio otatio ad formulas. We will also itroduce a brief overview of sequeces,
More informationProbability, Expectation Value and Uncertainty
Chapter 1 Probability, Expectatio Value ad Ucertaity We have see that the physically observable properties of a quatum system are represeted by Hermitea operators (also referred to as observables ) such
More information1 6 = 1 6 = + Factorials and Euler s Gamma function
Royal Holloway Uiversity of Lodo Departmet of Physics Factorials ad Euler s Gamma fuctio Itroductio The is a selfcotaied part of the course dealig, essetially, with the factorial fuctio ad its geeralizatio
More informationRotationally invariant integrals of arbitrary dimensions
September 1, 14 Rotatioally ivariat itegrals of arbitrary dimesios James D. Wells Physics Departmet, Uiversity of Michiga, A Arbor Abstract: I this ote itegrals over spherical volumes with rotatioally
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3337 HIKARI Ltd, www.mhikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationCEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering
CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio
More informationSECTION 2 Electrostatics
SECTION Electrostatics This sectio, based o Chapter of Griffiths, covers effects of electric fields ad forces i static (timeidepedet) situatios. The topics are: Electric field Gauss s Law Electric potetial
More informationDiscreteTime Systems, LTI Systems, and DiscreteTime Convolution
EEL5: DiscreteTime Sigals ad Systems. Itroductio I this set of otes, we begi our mathematical treatmet of discretetime s. As show i Figure, a discretetime operates or trasforms some iput sequece x [
More information4.3 Growth Rates of Solutions to Recurrences
4.3. GROWTH RATES OF SOLUTIONS TO RECURRENCES 81 4.3 Growth Rates of Solutios to Recurreces 4.3.1 Divide ad Coquer Algorithms Oe of the most basic ad powerful algorithmic techiques is divide ad coquer.
More informationAppendix: The Laplace Transform
Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio
More informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465475 HIKARI Ltd, www.mhikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More informationStochastic Matrices in a Finite Field
Stochastic Matrices i a Fiite Field Abstract: I this project we will explore the properties of stochastic matrices i both the real ad the fiite fields. We first explore what properties 2 2 stochastic matrices
More informationLinear Regression Demystified
Liear Regressio Demystified Liear regressio is a importat subject i statistics. I elemetary statistics courses, formulae related to liear regressio are ofte stated without derivatio. This ote iteds to
More informationThe Random Walk For Dummies
The Radom Walk For Dummies Richard A Mote Abstract We look at the priciples goverig the oedimesioal discrete radom walk First we review five basic cocepts of probability theory The we cosider the Beroulli
More informationChapter 9  CD companion 1. A Generic Implementation; The CommonMerge Amplifier. 1 τ is. ω ch. τ io
Chapter 9  CD compaio CHAPTER NINE CD9.2 CD9.2. Stages With Voltage ad Curret Gai A Geeric Implemetatio; The CommoMerge Amplifier The advaced method preseted i the text for approximatig cutoff frequecies
More informationTimeDomain Representations of LTI Systems
2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costatcoefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. Statevariable
More informationSNAP Centre Workshop. Basic Algebraic Manipulation
SNAP Cetre Workshop Basic Algebraic Maipulatio 8 Simplifyig Algebraic Expressios Whe a expressio is writte i the most compact maer possible, it is cosidered to be simplified. Not Simplified: x(x + 4x)
More informationMath 312 Lecture Notes One Dimensional Maps
Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 2123 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,
More informationPhysics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.
Physics 324, Fall 2002 Dirac Notatio These otes were produced by David Kapla for Phys. 324 i Autum 2001. 1 Vectors 1.1 Ier product Recall from liear algebra: we ca represet a vector V as a colum vector;
More informationStreamfunctionVorticity Formulation
StreamfuctioVorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctiovorticity formulatio was amog the first
More informationNICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =
AN INTRODUCTION TO SCHRÖDER AND UNKNOWN NUMBERS NICK DUFRESNE Abstract. I this article we will itroduce two types of lattice paths, Schröder paths ad Ukow paths. We will examie differet properties of each,
More information1 Adiabatic and diabatic representations
1 Adiabatic ad diabatic represetatios 1.1 BorOppeheimer approximatio The timeidepedet Schrödiger equatio for both electroic ad uclear degrees of freedom is Ĥ Ψ(r, R) = E Ψ(r, R), (1) where the full molecular
More informationShannon s noiseless coding theorem
18.310 lecture otes May 4, 2015 Shao s oiseless codig theorem Lecturer: Michel Goemas I these otes we discuss Shao s oiseless codig theorem, which is oe of the foudig results of the field of iformatio
More informationThe time evolution of the state of a quantum system is described by the timedependent Schrödinger equation (TDSE): ( ) ( ) 2m "2 + V ( r,t) (1.
Adrei Tokmakoff, MIT Departmet of Chemistry, 2/13/2007 11 574 TIMEDEPENDENT QUANTUM MECHANICS 1 INTRODUCTION 11 Timeevolutio for timeidepedet Hamiltoias The time evolutio of the state of a quatum system
More information18.01 Calculus Jason Starr Fall 2005
Lecture 18. October 5, 005 Homework. Problem Set 5 Part I: (c). Practice Problems. Course Reader: 3G 1, 3G, 3G 4, 3G 5. 1. Approximatig Riema itegrals. Ofte, there is o simpler expressio for the atiderivative
More informationChapter 6 Infinite Series
Chapter 6 Ifiite Series I the previous chapter we cosidered itegrals which were improper i the sese that the iterval of itegratio was ubouded. I this chapter we are goig to discuss a topic which is somewhat
More informationEndofYear Contest. ERHS Math Club. May 5, 2009
EdofYear Cotest ERHS Math Club May 5, 009 Problem 1: There are 9 cois. Oe is fake ad weighs a little less tha the others. Fid the fake coi by weighigs. Solutio: Separate the 9 cois ito 3 groups (A, B,
More informationQuantum Annealing for Heisenberg Spin Chains
LAUR #  Quatum Aealig for Heiseberg Spi Chais G.P. Berma, V.N. Gorshkov,, ad V.I.Tsifriovich Theoretical Divisio, Los Alamos Natioal Laboratory, Los Alamos, NM Istitute of Physics, Natioal Academy of
More informationCALCULATION OF FIBONACCI VECTORS
CALCULATION OF FIBONACCI VECTORS Stuart D. Aderso Departmet of Physics, Ithaca College 953 Daby Road, Ithaca NY 14850, USA email: saderso@ithaca.edu ad Dai Novak Departmet of Mathematics, Ithaca College
More informationPHY4905: NearlyFree Electron Model (NFE)
PHY4905: NearlyFree Electro Model (NFE) D. L. Maslov Departmet of Physics, Uiversity of Florida (Dated: Jauary 12, 2011) 1 I. REMINDER: QUANTUM MECHANICAL PERTURBATION THEORY A. Nodegeerate eigestates
More informationA statistical method to determine sample size to estimate characteristic value of soil parameters
A statistical method to determie sample size to estimate characteristic value of soil parameters Y. Hojo, B. Setiawa 2 ad M. Suzuki 3 Abstract Sample size is a importat factor to be cosidered i determiig
More informationMATHEMATICAL MODELLING OF ARCH FORMATION IN GRANULAR MATERIALS
6 th INTERNATIONAL MULTIDISCIPLINARY CONFERENCE MATHEMATICAL MODELLING OF ARCH FORMATION IN GRANULAR MATERIALS Istva eppler SZIE Faculty of Mechaics, H2103 Gödöllő Páter. 1., Hugary Abstract: The mathematical
More informationFinally, we show how to determine the moments of an impulse response based on the example of the dispersion model.
5.3 Determiatio of Momets Fially, we show how to determie the momets of a impulse respose based o the example of the dispersio model. For the dispersio model we have that E θ (θ ) curve is give by eq (4).
More informationMicroscopic Theory of Transport (Fall 2003) Lecture 6 (9/19/03) Static and Short Time Properties of Time Correlation Functions
.03 Microscopic Theory of Trasport (Fall 003) Lecture 6 (9/9/03) Static ad Short Time Properties of Time Correlatio Fuctios Refereces  Boo ad Yip, Chap There are a umber of properties of time correlatio
More informationPrinciple Of Superposition
ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give
More informationSeunghee Ye Ma 8: Week 5 Oct 28
Week 5 Summary I Sectio, we go over the Mea Value Theorem ad its applicatios. I Sectio 2, we will recap what we have covered so far this term. Topics Page Mea Value Theorem. Applicatios of the Mea Value
More informationThe axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.
5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =
More informationThe ztransform. 7.1 Introduction. 7.2 The ztransform Derivation of the ztransform: x[n] = z n LTI system, h[n] z = re j
The Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discretetime LTI systems. 7.
More informationCHAPTER 8 SYSTEMS OF PARTICLES
CHAPTER 8 SYSTES OF PARTICLES CHAPTER 8 COLLISIONS 45 8. CENTER OF ASS The ceter of mass of a system of particles or a rigid body is the poit at which all of the mass are cosidered to be cocetrated there
More informationR is a scalar defined as follows:
Math 8. Notes o Dot Product, Cross Product, Plaes, Area, ad Volumes This lecture focuses primarily o the dot product ad its may applicatios, especially i the measuremet of agles ad scalar projectio ad
More informationSequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018
CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisherprovided material Overview May mathematical
More information62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +
62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of
More informationRandom Models. Tusheng Zhang. February 14, 2013
Radom Models Tusheg Zhag February 14, 013 1 Radom Walks Let me describe the model. Radom walks are used to describe the motio of a movig particle (object). Suppose that a particle (object) moves alog the
More information1 1 2 = show that: over variables x and y. [2 marks] Write down necessary conditions involving first and secondorder partial derivatives for ( x0, y
Questio (a) A square matrix A= A is called positive defiite if the quadratic form waw > 0 for every ozero vector w [Note: Here (.) deotes the traspose of a matrix or a vector]. Let 0 A = 0 = show that:
More informationThe Pendulum. Purpose
The Pedulum Purpose To carry out a example illustratig how physics approaches ad solves problems. The example used here is to explore the differet factors that determie the period of motio of a pedulum.
More informationTHE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE ABSTRACT
Europea Joural of Egieerig ad Techology Vol. 3 No., 5 ISSN 56586 THE NUMERICAL SOLUTION OF THE NEWTONIAN FLUIDS FLOW DUE TO A STRETCHING CYLINDER BY SOR ITERATIVE PROCEDURE Atif Nazir, Tahir Mahmood ad
More informationUNIVERSITY OF CALIFORNIA  SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 116C. Problem Set 4. Benjamin Stahl. November 6, 2014
UNIVERSITY OF CALIFORNIA  SANTA CRUZ DEPARTMENT OF PHYSICS PHYS 6C Problem Set 4 Bejami Stahl November 6, 4 BOAS, P. 63, PROBLEM.5 The Laguerre differetial equatio, x y + ( xy + py =, will be solved
More informationMicroscopic traffic flow modeling
Chapter 34 Microscopic traffic flow modelig 34.1 Overview Macroscopic modelig looks at traffic flow from a global perspective, whereas microscopic modelig, as the term suggests, gives attetio to the details
More informationTEACHER CERTIFICATION STUDY GUIDE
COMPETENCY 1. ALGEBRA SKILL 1.1 1.1a. ALGEBRAIC STRUCTURES Kow why the real ad complex umbers are each a field, ad that particular rigs are ot fields (e.g., itegers, polyomial rigs, matrix rigs) Algebra
More informationRecursive Algorithms. Recurrences. Recursive Algorithms Analysis
Recursive Algorithms Recurreces Computer Sciece & Egieerig 35: Discrete Mathematics Christopher M Bourke cbourke@cseuledu A recursive algorithm is oe i which objects are defied i terms of other objects
More informationRecurrence Relations
Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial()); } Let t be the umber of multiplicatios eeded to calculate factorial(). The
More informationCS284A: Representations and Algorithms in Molecular Biology
CS284A: Represetatios ad Algorithms i Molecular Biology Scribe Notes o Lectures 3 & 4: Motif Discovery via Eumeratio & Motif Represetatio Usig Positio Weight Matrix Joshua Gervi Based o presetatios by
More informationMath 2784 (or 2794W) University of Connecticut
ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really
More informationFIR Filter Design: Part II
EEL335: DiscreteTime Sigals ad Systems. Itroductio I this set of otes, we cosider how we might go about desigig FIR filters with arbitrary frequecy resposes, through compositio of multiple siglepeak
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
More informationWhat is Physical Chemistry. Physical Chemistry for Chemical Engineers CHEM251. Basic Characteristics of a Gas
7/6/0 hysical Chemistry for Chemical Egieers CHEM5 What is hysical Chemistry hysical Chemistry is the study of the uderlyig physical priciples that gover the properties ad behaviour of chemical systems
More informationBertrand s Postulate
Bertrad s Postulate Lola Thompso Ross Program July 3, 2009 Lola Thompso (Ross Program Bertrad s Postulate July 3, 2009 1 / 33 Bertrad s Postulate I ve said it oce ad I ll say it agai: There s always a
More informationmx bx kx F t. dt IR I LI V t, Q LQ RQ V t,
Lecture 5 omplex Variables II (Applicatios i Physics) (See hapter i Boas) To see why complex variables are so useful cosider first the (liear) mechaics of a sigle particle described by Newto s equatio
More informationPolynomials with Rational Roots that Differ by a Nonzero Constant. Generalities
Polyomials with Ratioal Roots that Differ by a Nozero Costat Philip Gibbs The problem of fidig two polyomials P(x) ad Q(x) of a give degree i a sigle variable x that have all ratioal roots ad differ by
More informationThe standard deviation of the mean
Physics 6C Fall 20 The stadard deviatio of the mea These otes provide some clarificatio o the distictio betwee the stadard deviatio ad the stadard deviatio of the mea.. The sample mea ad variace Cosider
More informationDETERMINATION OF MECHANICAL PROPERTIES OF A NON UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.
ICSV4 Cairs Australia 9 July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics
More informationDefinition 4.2. (a) A sequence {x n } in a Banach space X is a basis for X if. unique scalars a n (x) such that x = n. a n (x) x n. (4.
4. BASES I BAACH SPACES 39 4. BASES I BAACH SPACES Sice a Baach space X is a vector space, it must possess a Hamel, or vector space, basis, i.e., a subset {x γ } γ Γ whose fiite liear spa is all of X ad
More informationA sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam will cover..9. This sheet has three sectios. The first sectio will remid you about techiques ad formulas that you should kow. The secod gives a umber of practice questios for you
More informationDamped Vibration of a Nonprismatic Beam with a Rotational Spring
Vibratios i Physical Systems Vol.6 (0) Damped Vibratio of a Noprismatic Beam with a Rotatioal Sprig Wojciech SOCHACK stitute of Mechaics ad Fudametals of Machiery Desig Uiversity of Techology, Czestochowa,
More informationECES352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECES352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More informationFinal Review for MATH 3510
Fial Review for MATH 50 Calculatio 5 Give a fairly simple probability mass fuctio or probability desity fuctio of a radom variable, you should be able to compute the expected value ad variace of the variable
More informationOptimally Sparse SVMs
A. Proof of Lemma 3. We here prove a lower boud o the umber of support vectors to achieve geeralizatio bouds of the form which we cosider. Importatly, this result holds ot oly for liear classifiers, but
More informationENGI Series Page 601
ENGI 3425 6 Series Page 601 6. Series Cotets: 6.01 Sequeces; geeral term, limits, covergece 6.02 Series; summatio otatio, covergece, divergece test 6.03 Stadard Series; telescopig series, geometric series,
More information2C09 Design for seismic and climate changes
2C09 Desig for seismic ad climate chages Lecture 02: Dyamic respose of sigledegreeoffreedom systems I Daiel Grecea, Politehica Uiversity of Timisoara 10/03/2014 Europea Erasmus Mudus Master Course Sustaiable
More informationFluid Physics 8.292J/12.330J % (1)
Fluid Physics 89J/133J Problem Set 5 Solutios 1 Cosider the flow of a Euler fluid i the x directio give by for y > d U = U y 1 d for y d U + y 1 d for y < This flow does ot vary i x or i z Determie the
More information6.867 Machine learning, lecture 7 (Jaakkola) 1
6.867 Machie learig, lecture 7 (Jaakkola) 1 Lecture topics: Kerel form of liear regressio Kerels, examples, costructio, properties Liear regressio ad kerels Cosider a slightly simpler model where we omit
More informationA Lattice Green Function Introduction. Abstract
August 5, 25 A Lattice Gree Fuctio Itroductio Stefa Hollos Exstrom Laboratories LLC, 662 Nelso Park Dr, Logmot, Colorado 853, USA Abstract We preset a itroductio to lattice Gree fuctios. Electroic address:
More informationSummary: CORRELATION & LINEAR REGRESSION. GC. Students are advised to refer to lecture notes for the GC operations to obtain scatter diagram.
Key Cocepts: 1) Sketchig of scatter diagram The scatter diagram of bivariate (i.e. cotaiig two variables) data ca be easily obtaied usig GC. Studets are advised to refer to lecture otes for the GC operatios
More informationAssignment 2 Solutions SOLUTION. ϕ 1 Â = 3 ϕ 1 4i ϕ 2. The other case can be dealt with in a similar way. { ϕ 2 Â} χ = { 4i ϕ 1 3 ϕ 2 } χ.
PHYSICS 34 QUANTUM PHYSICS II (25) Assigmet 2 Solutios 1. With respect to a pair of orthoormal vectors ϕ 1 ad ϕ 2 that spa the Hilbert space H of a certai system, the operator Â is defied by its actio
More informationThe Riemann Zeta Function
Physics 6A Witer 6 The Riema Zeta Fuctio I this ote, I will sketch some of the mai properties of the Riema zeta fuctio, ζ(x). For x >, we defie ζ(x) =, x >. () x = For x, this sum diverges. However, we
More informationSRC Technical Note June 17, Tight Thresholds for The Pure Literal Rule. Michael Mitzenmacher. d i g i t a l
SRC Techical Note 1997011 Jue 17, 1997 Tight Thresholds for The Pure Literal Rule Michael Mitzemacher d i g i t a l Systems Research Ceter 130 Lytto Aveue Palo Alto, Califoria 94301 http://www.research.digital.com/src/
More informationABOUT CHAOS AND SENSITIVITY IN TOPOLOGICAL DYNAMICS
ABOUT CHAOS AND SENSITIVITY IN TOPOLOGICAL DYNAMICS EDUARD KONTOROVICH Abstract. I this work we uify ad geeralize some results about chaos ad sesitivity. Date: March 1, 005. 1 1. Symbolic Dyamics Defiitio
More informationINTEGRATION BY PARTS (TABLE METHOD)
INTEGRATION BY PARTS (TABLE METHOD) Suppose you wat to evaluate cos d usig itegratio by parts. Usig the u dv otatio, we get So, u dv d cos du d v si cos d si si d or si si d We see that it is ecessary
More informationBasics of Probability Theory (for Theory of Computation courses)
Basics of Probability Theory (for Theory of Computatio courses) Oded Goldreich Departmet of Computer Sciece Weizma Istitute of Sciece Rehovot, Israel. oded.goldreich@weizma.ac.il November 24, 2008 Preface.
More informationOlli Simula T / Chapter 1 3. Olli Simula T / Chapter 1 5
Sigals ad Systems Sigals ad Systems Sigals are variables that carry iformatio Systemstake sigals as iputs ad produce sigals as outputs The course deals with the passage of sigals through systems T6.4
More informationThe Gamma function Michael Taylor. Abstract. This material is excerpted from 18 and Appendix J of [T].
The Gamma fuctio Michael Taylor Abstract. This material is excerpted from 8 ad Appedix J of [T]. The Gamma fuctio has bee previewed i 5.7 5.8, arisig i the computatio of a atural Laplace trasform: 8. ft
More informationCosmology Tutorial 8
Cosmology 018  Tutorial 8 Jue 5, 018 I this ote I have largely followed Ladau & Lifshitz Fluid Mechaics (1987, Weiberg s Graviatios ad Cosmology (197 ad otes by Yacie AliHaïmoud. We have mostly covered
More informationThe BornOppenheimer approximation
The BorOppeheimer approximatio 1 Rewritig the Schrödiger equatio We will begi from the full timeidepedet Schrödiger equatio for the eigestates of a molecular system: [ P 2 + ( Pm 2 + e2 1 1 2m 2m m
More informationPhysics 232 Gauge invariance of the magnetic susceptibilty
Physics 232 Gauge ivariace of the magetic susceptibilty Peter Youg (Dated: Jauary 16, 2006) I. INTRODUCTION We have see i class that the followig additioal terms appear i the Hamiltoia o addig a magetic
More informationU8L1: Sec Equations of Lines in R 2
MCVU U8L: Sec. 8.9. Equatios of Lies i R Review of Equatios of a Straight Lie (D) Cosider the lie passig through A (,) with slope, as show i the diagram below. I poit slope form, the equatio of the lie
More informationSEQUENCES AND SERIES
9 SEQUENCES AND SERIES INTRODUCTION Sequeces have may importat applicatios i several spheres of huma activities Whe a collectio of objects is arraged i a defiite order such that it has a idetified first
More informationFourier Series and the Wave Equation
Fourier Series ad the Wave Equatio We start with the oedimesioal wave equatio u u =, x u(, t) = u(, t) =, ux (,) = f( x), u ( x,) = This represets a vibratig strig, where u is the displacemet of the strig
More informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
More informationExact Solutions for a Class of Nonlinear Singular TwoPoint Boundary Value Problems: The Decomposition Method
Exact Solutios for a Class of Noliear Sigular TwoPoit Boudary Value Problems: The Decompositio Method Abd Elhalim Ebaid Departmet of Mathematics, Faculty of Sciece, Tabuk Uiversity, P O Box 741, Tabuki
More informationLecture 25 (Dec. 6, 2017)
Lecture 5 8.31 Quatum Theory I, Fall 017 106 Lecture 5 (Dec. 6, 017) 5.1 Degeerate Perturbatio Theory Previously, whe discussig perturbatio theory, we restricted ourselves to the case where the uperturbed
More informationChapter 14: Chemical Equilibrium
hapter 14: hemical Equilibrium 46 hapter 14: hemical Equilibrium Sectio 14.1: Itroductio to hemical Equilibrium hemical equilibrium is the state where the cocetratios of all reactats ad products remai
More informationHE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT. Examples:
5.6 4 Lecture #34 page HE ATOM & APPROXIMATION METHODS MORE GENERAL VARIATIONAL TREATMENT Do t restrict the wavefuctio to a sigle term! Could be a liear combiatio of several wavefuctios e.g. two terms:
More informationLecture 4 Conformal Mapping and Green s Theorem. 1. Let s try to solve the following problem by separation of variables
Lecture 4 Coformal Mappig ad Gree s Theorem Today s topics. Solvig electrostatic problems cotiued. Why separatio of variables does t always work 3. Coformal mappig 4. Gree s theorem The failure of separatio
More information