Brownian motion in a granular gas

Size: px
Start display at page:

Download "Brownian motion in a granular gas"

Transcription

1 PHYSICAL REVIEW E VOLUME 60, NUMBER 6 DECEMBER 1999 Brownin motion in grnulr gs J. Jvier Brey, M. J. Ruiz-Montero, nd R. Grcí-Rojo Físic Teóric, Universidd de Sevill, E Sevill, Spin Jmes W. Dufty Deprtment of Physics, University of Florid, Ginesville, Florid Received 2 June 1999 The dynmics of hevy prticle in gs of much lighter prticles is studied vi the Boltzmnn-Lorentz eqution with inelstic collisions mong ll prticles. A forml expnsion in the rtio of gs to tgged prticle mss trnsforms the Boltzmnn-Lorentz eqution into Fokker Plnck eqution. The predictions of the ltter re tested here using direct Monte Crlo simultion of the Boltzmnn-Lorentz eqution. Excellent greement is obtined for the pproch to homogeneous cooling stte, the temperture of tht stte, pproch to diffusion, nd the dependence of the diffusion constnt on dissiption prmeters. Some results from moleculr-dynmics simultions re lso presented nd shown to gree with the theoreticl predictions. S X PACS numbers: Rm, Dd, I. INTRODUCTION The description of low-density, rpid grnulr flow by Boltzmnn kinetic theory hs received much ttention in recent yers 1,2. A primry objective hs been the derivtion of hydrodynmic equtions nd expressions for the trnsport coefficients ppering in them s functions of the coefficient of restitution. The most ccurte pproch uses generliztion of the Chpmn-Enskog method known from the corresponding nlysis of the Boltzmnn eqution for elstic collisions 3,4. This method ssumes the existence of norml solution whose spce nd time dependence occurs only through the hydrodynmic fields. Implicit in this ssumption is the rpid relxtion of nonhydrodynmic excittions so tht the hydrodynmic description domintes on longer time scle. Such seprtion of time scles hs been questioned in the cse of grnulr flow, prticulrly for lrge degrees of inelsticity 5. The bsis for this concern seems to be the dditionl time scle set by the cooling of the homogeneous reference stte, such tht the hydrodynmic time scles re not simply determined by the degree of sptil inhomogeneity. However, detiled study of self-diffusion 6 nd kinetic models bsed on the Boltzmnn eqution 7 suggests this concern is not justified. For these cses it is shown tht the microscopic excittions lwys decy on times short compred to ll hydrodynmic times, including tht defined by the cooling rte. To further reinforce the cse for hydrodynmic description, n exct nlysis of the Boltzmnn-Lorentz kinetic eqution for tgged prticle in freely evolving gs hs been performed in the limit of symptoticlly lrge reltive mss for the tgged prticle 7. In this limit, the Boltzmnn- Lorentz eqution reduces to Fokker-Plnck eqution. An exct nlysis of the spectrum for this eqution confirms the seprtion of time scles for ll degrees of tgged prticle inelsticity. It remins to confirm the vlidity of the forml symptotic nlysis leding from the Boltzmnn-Lorentz eqution to the Fokker-Plnck eqution. One objective here is to provide this confirmtion on the bsis of direct Monte Crlo simultion of the Boltzmnn-Lorentz eqution. In the next section the Boltzmnn-Lorentz nd its Fokker- Plnck limit re reclled. A chnge of spce nd time vribles provides n exct mp of the Fokker-Plnck eqution for inelstic collisions to tht for elstic collisions. Consequently, ll the known results from the ltter cse for velocity relxtion nd the pproch to diffusive stge trnslte exctly to the cse of inelstic collisions. Severl exmples re considered explicitly. The qulittive differences occur only through the chnges in spce nd time scles. Some of the most interesting differences re s follows. i The renormlized time scle is relted logrithmiclly to rel time. Accordingly, velocity relxtion nd pproch to hydrodynmics is lgebric rther thn exponentil. ii The long time limit of the tgged prticle distribution for the homogeneous stte is Gussin, lthough the bth prticle distribution is non-gussin. iii The time-dependent temperture of the tgged prticle Gussin differs from the temperture of the surrounding bth, lthough the cooling rtes of both become the sme t long times. iv The men-squre displcement pproches ln t for lrge times with coefficient tht determines the diffusion coefficient. v The diffusion eqution t long times hs the usul form, lthough its solution is qulittively different from tht for the elstic cse due to the time dependence of the temperture in the diffusion coefficient. In Sec. III the direct simultion Monte Crlo method to obtin numericl solutions of the Boltzmnn eqution is described briefly. Comprisons re mde between the simultions of the Boltzmnn-Lorentz eqution nd the predictions of the Fokker-Plnck eqution for tgged prticle whose mss is 100 times tht of the surrounding gs prticles. The Mxwellin distribution is confirmed, s is the pproch to different tempertures for the bth nd tgged prticle t the sme cooling rtes. The men-squre displcement is mesured nd the diffusion coefficient extrcted from its long time behvior to confirm the detiled dependence on inelsticity s predicted by the Fokker-Plnck eqution. The results provide strong support for the ccurcy of the Fokker-Plnck X/99/606/71748/$15.00 PRE The Americn Physicl Society

2 PRE 60 BROWNIAN MOTION IN A GRANULAR GAS 7175 eqution s representtion of the Boltzmnn-Lorentz eqution for lrge reltive mss. Of course, ll the bove discussion relies on the vlidity of the inelstic Boltzmnn-Lorentz eqution to describe mssive tgged prticle in freely evolving gs. While this vlidity is well estblished for elstic collisions, it is sometimes questioned for inelstic collisions. To provide support for this cse, we hve performed moleculr-dynmics simultions for this system in two dimensions. The results presented in Sec. IV confirm the nlysis bsed on the kinetic eqution, indicting clerly the vlidity of kinetic equtions to study rpid grnulr flows. This section lso includes summry nd discussion of the min results in the pper. II. FOKKER-PLANCK EQUATION AND ITS SOLUTION A tgged prticle immersed in low-density gs is considered. The gs is formed by hrd spheres (d3) or disks (d2) of mss m g nd dimeter g, nd the tgged prticle is lso hrd sphere or disk, but with mss m nd dimeter. All prticles re smooth nd collide inelsticlly. Collisions re chrcterized by velocity-independent coefficients of norml restitution. For the collisions between the gs prticles it will be denoted by g, while tht for collisions between the tgged prticle nd gs prticles will be represented by. The probbility density F(r,v,t) describing the dynmics of the tgged prticle obeys the Boltzmnn-Lorentz eqution t v FJr,v,tF, f. The collision opertor J is given by 7,8 Jr,v,tF, f 0 d1 dv 1 dˆ g ˆ g ˆ 2 Fr,v,tfr,v 1,t Fr,v,tfr,v 1,t, where f (r,v,t) is the corresponding distribution for the surrounding gs prticles, is the Heviside step function, ˆ is unit vector pointing from the center of the gs prticle 1 to the center of the tgged prticle t contct, nd 0 ( g )/2. The precollisionl or restituting velocities v nd v 1 re given by vv 1 1 g ˆ ˆ, v 1 v g ˆ ˆ, 2.3 with gvv 1 the reltive velocity nd m g /m the rtio of gs to tgged prticle mss. The surrounding gs is tken to be in its homogeneous cooling stte HCS s determined from the solution to the nonliner Boltzmnn eqution. Its distribution function hs the scling form 1 f H v,tn g v d v g t, 2.4 v g t where n g is the constnt number density of the gs, v g (t) 2k B T g (t)/m g 1/2 is the therml velocity of the gs prticles t time t with k B being the Boltzmnn constnt, nd T g (t) is the temperture of the gs which cools ccording to the eqution dt g t tt dt g t. 2.5 The cooling rte (t) depends on time only through T g (t) nd is determined from the second moment of the Boltzmnn eqution for the gs 8, t1 2 g n g d1 g v g t 2d d1/2 d3 2 dv dv 1 g 3 vv The solution to Eq. 2.5 is T g tt g t, 2.7 showing tht the temperture of the gs decreses s t 2 for lrge times if the system remins in the HCS. The explicit form of will be not given here, but it is known in the so-clled first Sonine pproximtion 1,9. In Ref. 7 it ws shown tht the Boltzmnn-Lorentz eqution reduces to Fokker-Plnck eqution for symptoticlly smll, t v Fr,v,t e t v v k BT g t m vfr,v,t, 2.8 where e (t) is the sme friction coefficient s for elstic collisions, except s function of the time-dependent temperture T g (t), e t 4d1/2 d1 0 n dd/2 g v g t, 2.9 where v g (t) is defined in terms of T g (t) following Eq. 2.4 bove. All effects of inelstic collisions mong the gs prticles in Eq. 2.8 pper through the time dependence of T g (t), while the inelsticity of collisions between the tgged prticle nd gs prticles mnifests itself only through the prmeter (1)/2. The Fokker-Plnck eqution 2.8 cn be mpped onto the corresponding eqution for elstic collisions using the dimensionless vribles v* v v 0 t, where we hve introduced r*1 et v 0 t r, t*1 0 t dte t, 2.10

3 7176 BREY, RUIZ-MONTERO, GARCÍA-ROJO, AND DUFTY PRE 60 v 0 t 1 1/2 2k BT g t m 1/2, t 2 e t Since both (t) nd e (t) re proportionl to T 1/2 g (t), it follows tht is time-independent quntity. The bove definitions pply only for 1, point discussed further in the finl section. Here we only note tht it is necessry condition in the derivtion of Eq In terms of the new vribles, the Fokker Plnck eqution 2.8 becomes t * v* F*r*,v*,t* r* v* v* 1 F*r*,v*,t*, v* with the scled probbility density F* given by v 2d 0 t F*r*,v*,t* 1 e t Fr,v,t, 2.13 d which is lso normlized to unity. Eqution 2.12 is the sme s the dimensionless form of the Fokker-Plnck eqution for elstic collisions ( g 1). Consequently, the physicl properties of mssive tgged prticle with inelstic collisions moving in gs in the HCS re the sme s those for n elstic prticle in n equilibrium gs; the only differences re the relevnt spce nd time scles. For exmple, t* e t for elstic collisions while for inelstic collisions the reltionship is t* 1 ln t, 2.14 nd the time scle is stretched logrithmiclly. The generl solution of Eq to the initil vlue problem for n unbounded system is well known 10: F*r*,v*,t* dr* dv* G*r*,v*,t*;r*,v*,0 with G*r*,v*,t*;r*,v*,0 et*bt* 2 F*r*,v*,0, d/2 expet*r*r* ct* v*v* 2 bt*v*e t* v* Consider first sptilly homogeneous initil stte for the tgged prticle F*(r*,v*,0)F*(v*,0). Then Eq simplifies to F*v*,t* bt* d/2 dv* expbt* v*e t* v* 2 F*v*,0. For t*1, b(t*) 1 nd, consequently, F*v*,t* F M * v* 1 * d/2 e v*2, where * is the volume of the system mesured in the reduced length scle defined in Eq Thus, for generl homogeneous initil conditions, the distribution function pproches exponentilly fst sttionry Mxwellin distribution in the reduced units. In terms of the originl vribles this result is Fv,t F M v,t 1 d/2 m e 2k B T t mv2 /2k B T (t), 2.20 with the temperture prmeter T (t) given by T tt g t Interestingly, the tgged prticle pproches homogeneous cooling stte with scling form similr to Eq. 2.4 for the gs, but with two importnt differences. First, the cooling stte for the tgged prticle is Gussin while tht for the gs is not. Second, the cooling temperture for the tgged prticle is different from tht for the gs, lthough the cooling rtes re the sme. To elborte on this ltter feture, it is convenient to define more generlly the kinetic temperture for the tgged prticle by d 2 k BTt dr dv 1 2 mv2 Fr,v,t k BT g t 1 dr* dv* v* 2 F*r*,v*,t. Evluting the integrl using Eq gives Tt T g t 1 T0 e2t* T g In the bove expressions we hve defined 1 et* 2t*2ct*, bt* 1 1e 2t*, ct* 1et* 1e t* The pproch of the tgged prticle temperture to the symptotic vlue T (t) is exponentilly fst on the reduced time scle. Note tht T (t)t g (t) even in the cse of elstic collisions for the tgged prticle (1), if the gs prticle collisions re inelstic. Conversely, if gs prticle collisions re elstic but the tgged prticle collisions re inelstic, then 0 nd the tgged prticle pproches the constnt temperture of the gs, but still differs from it by fctor of.

4 PRE 60 BROWNIAN MOTION IN A GRANULAR GAS 7177 The result in Eq pplies even for rbitrry inhomogeneous sttes of the Brownin prticle. Consider next n initil stte whose velocity distribution is given by the homogeneous cooling Mxwellin but with sptil inhomogeneity, F*r*,v*,0n*r*,0 d/2 e v*2, 2.24 where n*(r*,0) is the probbility density for finding the tgged prticle t position r* t time t0. The corresponding quntity t time t* is obtined by integrtion of F*(r*,v*,t*), s given in Eq. 2.15, with respect to v*, giving d/2 d n*r*,t* dr* where 2l* 2 t* expdr*r* 2 /2l* 2 t*n*r*,0, 2.25 l* 2 t*dt*1e t* It is esily verified tht l* 2 (t*) is the men-squre displcement of the tgged prticle in terms of the dimensionless vribles. Moreover, it follows directly from Eq tht n*(r*,t*) obeys the extended diffusion eqution t * n*r*,t*1et* D* * 2 n*r*,t*, 2.27 with D*1/2. This eqution is exct for ll times if the initil condition hs the ssumed form, nd shows most directly the pproch to hydrodynmic stge, i.e., the usul diffusion eqution pplies exponentilly fst for t*1. The vlue D*1/2 is consistent with the Einstein result D* 1 2d lim t* 1 l* 2 t* t* In terms of the originl vribles, the symptotic diffusion eqution is with t nr,tdt 2 nr,t, Dt D et Here D e (t)k B T g (t)/m e (t) is the sme s the diffusion coefficient for elstic collisions, except s function of the time-dependent gs temperture. For given temperture, the diffusion coefficient is seen to be enhnced by fctor of (1) 2, depending only on the inelsticity of gs collisions. For more generl initil conditions, the dynmics is qulittively similr with rpid velocity relxtion followed by sptil diffusion. In the next section, severl predictions from this Fokker-Plnck nlysis re tested ginst direct simultion of the corresponding property from the Boltzmnn- Lorentz eqution. III. MONTE CARLO SIMULATION OF THE BOLTZMANN-LORENTZ EQUATION To test the forml nlysis leding from the Boltzmnn- Lorentz eqution to the Fokker-Plnck, the direct simultion Monte Crlo DSMC method 11 hs been pplied to the former. Specificlly, the predictions of velocity relxtion to Mxwellin t new temperture nd the pproch to diffusion with n enhnced diffusion coefficient hve been investigted for system of hrd spheres with mss rtio 10 2, equl dimeters for ll prticles, nd coefficient of norml restitution for collisions between gs prticles g These vlues ssure the necessry condition 1 for 01, s shown in the Appendix. The DSMC for the Boltzmnn eqution with elstic collisions is well described in Ref. 11. Its dpttion to the Boltzmnn-Lorentz eqution with inelstic collisions is strightforwrd. As indicted in Eq. 2.1, the distribution function for the gs, f s, is required for input, nd is tken here to be the homogeneous cooling solution to the Boltzmnn eqution f H. It is known tht this solution is unstble to long-wvelength sptil perturbtions t sufficiently long times, so ny possible interference of relxtion processes by this instbility is not ddressed. However, for the chosen vlue of g the homogeneous cooling stte is stble on the time scles studied 12. In the simultions 10 5 trjectories of tgged prticle hve been generted. Independently of its position, collisions of the tgged prticle lwys took plce with prticles of homogeneous gs whose velocity distribution ws given by Eq As lredy mentioned, the exct form of is not known, nd the expression obtined in the first Sonine pproximtion 1,9 ws used to generte the velocity of the colliding gs prticles. Thus, the trjectories of the gs prticles re not required, which increses gretly the efficiency of the numericl simultion nd voids the introduction of specific boundry conditions the system is formlly considered s infinite. A similr method lredy hs been used to study self-difussion in low-density grnulr flow 6. In typicl run, collisions between the tgged prticle nd the fluid prticles were considered s uncoupled during time step t 0, chosen much smller thn the initil verge collision time of the tgged prticle. This mens tht the position of the tgged prticle ws chnged t constnt velocity for n intervl t 0 between every two pplictions of the collision lgorithm. For the collisions, the velocities of the prticles were generted from the HCS distribution, s mentioned bove, with gs temperture determined from the lw given by Eq. 2.5 nd considered to be constnt in the time intervl t 0. This implies tht t 0 hs to be chosen lso much smller thn the time chrcterizing the cooling of the gs, 2/(0). In the simultions we report here we used t /(0). The initil condition in ll simultions ws Gussin velocity distribution with the sme temperture s the surrounding fluid. As first test of the Fokker-Plnck limit, the predicted pproch to Gussin cooling stte is studied. The ssoci-

5 7178 BREY, RUIZ-MONTERO, GARCÍA-ROJO, AND DUFTY PRE 60 FIG. 1. Plot of the rtio of the tgged prticle temperture to the surrounding gs temperture s function of the dimensionless reduced time t* defined in Eq The rtio of gs to tgged prticle mss is 10 2 nd the coefficient of restitution for the gs collisions is g The solid lines re the predictions from the Fokker-Plnck eqution nd the symbols from the DSMC method of the Boltzmnn-Lorentz eqution. ted scling velocity is defined in terms of the temperture T(t) of the tgged prticle, chrcterizing its kinetic energy. It is expected tht this temperture pproches n symptotic vlue proportionl to the gs temperture T g (t), ccording to Eq Figure 1 shows the time evolution of the rtio T(t)/T g (t) s obtined from the numericl simultion of the Boltzmnn-Lorentz eqution for 0.5, 0.7, nd In ll cses it is seen tht the rtio evolves from the unit initil condition until reching stedy vlue. Also plotted is the theoreticl prediction. Figure 2 shows detiled comprison of the numericl symptotic vlues with Eq over this rnge of vlues of. The greement is very good considering tht corrections to the Fokker Plnck limit re of order FIG. 2. Asymptotic vlue of the rtio between the tgged prticle temperture T nd the gs temperture T g s function of the coefficient of norml restitution for collisions between the tgged prticle nd the fluid prticles. The mss rtio nd the coefficient of restitution for gs collisions re the sme s in Fig. 1. The solid line is from the Fokker-Plnck eqution nd the dots from the numericl simultion of the Boltzmnn-Lorentz eqution. FIG. 3. Time evolution of the reduced fourth velocity moment of the tgged prticle distribution. Here time is mesured in units of /v 0, where is the men free pth. The continuous line corresponds to 0.95 nd the dshed line to 0.5. The stright lines represent the fourth moment of the HCS distribution for the fluid t 0.95 continuous nd 0.5 dshed. 1/2. Interestingly, we hve found tht the greement between theory nd simultions improves significntly if the mss rtio nd g re chnged to reduce the vlue of. The Gussin chrcter of the tgged prticle distribution cn be studied vi its fourth moment, or the normlized expression 3v 4 /5v 2 2, which hs the vlue unity for Gussin. Here it is v n 1 dr dv v n Fr,v,t. 3.1 Figure 3 shows the simultion vlues for the normlized moment s function of time for 0.5 nd Also shown re the corresponding results for the gs distribution, i.e., for the velocity distribution in Eq The ltter shows significnt devitions from unity wheres the tgged prticle results confirm the Gussin even t strong dissiption. Although we hve plotted in the figure the results corresponding to two extreme vlues of, similr behvior hs been obtined lso for severl intermedite vlues. It could be climed tht the Gussin distribution observed in the simultions is in some wy influenced by the initil Gussin distribution. To clrify this point we hve lso considered initil distributions fr from the Gussin, nmely uniform distribution with zero men. It ws observed tht the distribution evolved towrds Gussin very fst, indicting tht the use of n initil Gussin does not limit t ll the results presented here. The pproch to diffusion cn be studied vi the mensqure displcement of the tgged prticle. In Fig. 4, comprison of the simultions with the universl form given by Eq for 0.5, 0.7, nd 0.99 is presented. The results confirm both the pproch to diffusion nd the predicted limiting form with D*1/2. Figure 5 shows more detiled comprison with the predicted diffusion coefficient over wide rnge of vlues. Agin, the greement is very good, since the discrepncies re smller thn 2%.

6 PRE 60 BROWNIAN MOTION IN A GRANULAR GAS 7179 FIG. 4. Time evolution of the men-squre displcement of the tgged prticle for three different vlues of the coefficient of restitution for collisions between the tgged prticle nd the gs prticles. Time nd length re mesured in the reduced dimensionless units defined in the min text. The solid line is the theoreticl prediction given by Eq IV. DISCUSSION The Fokker-Plnck eqution, obtined from the Boltzmnn-Lorentz eqution in the limit of smll rtio of gs to tgged prticle mss, llows n exct nlysis of the tgged prticle dynmics. The surrounding gs is ssumed to be in the homogeneous cooling stte. In prticulr, the study shows the cler seprtion of time scles required for hydrodynmic description. For tgged prticle dynmics the ltter refers to diffusion. The objective here hs been to confirm the vlidity of the symptotic nlysis leding to the Fokker-Plnck limit by direct Monte Crlo simultion of the Boltzmnn-Lorentz eqution. The excellent greement found for the vrious properties studied provides convincing evidence for the vlidity of the Fokker-Plnck eqution nd its exct consequences. The seprtion of time scles between microscopic kinetic nd diffusive modes is quite cler, once the complictions due to cooling hve been suppressed by the chnge of FIG. 5. Comprison between the numericl vlue of the reduced self-diffusion coefficient obtined from the Boltzmnn-Lorentz eqution dots nd the predicted vlue D*1/2 solid line. FIG. 6. Distribution function of the tgged prticle normlized by the Gussin for three different vlues of the restitution coefficient. For the three of them it is 1, nd the Fokker-Plnck eqution is not expected to hold. The velocity is mesured in the reduced units defined in Eq vribles in Eq In the new vribles, the description is independent of nd, therefore, pplies for rbitrry degree of dissiption for the tgged prticle. However, the gs is implicitly restricted to wek dissiption due to the condition 1, where is defined by Eq The explicit form for given in the Appendix shows tht this condition requires 1 2 g. The origin of this condition is requirement tht the rtio of tgged prticle temperture to gs temperture does not grow in time, since it ppers s fctor of in the symptotic nlysis 7. The relevnce of this condition for the Fokker-Plnck limit is estblished in Fig. 6 showing the distribution function F normlized by the Gussin F M for stronger gs dissiption, g 0.95, nd for severl vlues of the restitution coefficient for collisions between the tgged prticle nd the fluid, 0.99, 0.9, nd 0.8. The corresponding vlues for re 1.73, 1.81, nd 1.915, respectively. The devitions from Mxwellin re now lrge nd the Fokker-Plnck description, which predicts Mxwellin velocity distribution, is no longer vlid. Let us stress tht the restriction 1 ffects only the Fokker-Plnck limit nd not the more generl picture of pproch to diffusion, s described by the Boltzmnn-Lorentz eqution. This hs been demonstrted in Ref. 6, where the men-squre displcement hs been simulted for mechniclly identicl prticles for 0.6 g 1, which includes conditions of strong gs prticle dissiption nd consequently 1. All our previous nlysis is bsed on the ccurcy of the Boltzmnn-Lorentz eqution to describe the time evolution of tgged prticle in low-density gs, when ll the collisions re inelstic. For the prticulr cse of self-diffusion the vlidity of the kinetic eqution hs been explicitly shown in Ref. 6, by compring its predictions with the results obtined from moleculr-dynmics simultion of system of inelstic disks. Here we present some dditionl evidence when the Brownin limit of mssive prticle is considered. From computtionl point of view, the study of time evolution of mssive prticle in bth, whose stte must not be perturbed by the motion of the prticle, is very demnding simultion, requiring lot of computing time in order to get

7 7180 BREY, RUIZ-MONTERO, GARCÍA-ROJO, AND DUFTY PRE 60 FIG. 7. Time evolution of the men-squre displcement of tgged disk in low-density gs of inelstic prticles, using moleculr dynmics. The mss rtio is nd two vlues for the coefficient of restitution for collisions of the tgged prticles hve been considered: 0.9 squres nd 0.8 circles. The coefficient of restitution for collisions mong fluid prticles is g Quntities re mesured in the dimensionless reduced units defined in the min text. The solid line is the theoreticl prediction from the Fokker-Plnck description, Eq significnt level of sttistics. We hve simulted system of hrd disks in squre domin with periodic boundry conditions. The gs ws composed of 3025 prticles nd the results we present re verges over 500 trjectories of the Brownin prticle, ech of them generted by different computer run. Agin the tgged nd gs prticles were of the sme size. The initil condition ws generted by running the gs, without tgged prticle, for period of time lrge enough to llow the system to rech the homogeneous cooling stte. Then, one of the gs prticles, rndomly selected, ws substituted by the tgged prticle, i.e., its mss ws chnged to the new vlue. The simultion technique ws bsed on the event driven lgorithm 13. Figure 7 shows the results obtined for the time evolution of the men-squre displcement of the tgged prticle in system with g 0.99 nd mss rtio The density of the gs is n , which corresponds to solid frction Two vlues of the restitution coefficient for collisions of the tgged prticle hve been considered, 0.9 nd 0.8. The solid line is the prediction given by Eq Agin good greement is obtined, providing very strong test on the vlidity of the kinetic eqution even before the system reches the hydrodynmic regime. Similr results were obtined for other vlues of the prmeters of the system, nmely those considered in the DSMC results reported in the preceding section. In the figure it is observed tht the discrepncy between theory nd simultion increses with time. This is sttistics effect which decreses s the number of trjectories increses. This is more clerly seen in Fig. 8, where we hve represented the time evolution of the temperture of the tgged prticle for the simultions corresponding to 0.9 in Fig. 7. Although the results re consistent with the predictions from the Fokker-Plnck eqution, the fluctutions re too lrge to mke detiled quntittive comprison. This is not the cse FIG. 8. Time evolution of the temperture of the Brownin prticle for the sme system s considered in Fig. 7 with 0.9. The temperture is scled with the time-dependent temperture of the gs. The solid line is the theoreticl prediction given by Eq Time is mesured in the reduce unit defined in the text. for the men-squre displcement, s seen in Fig. 7. In this sense, when compring the results obtined by using the DSMC method with those from MD, we must tke into ccount tht the number of trjectories we hve considered in the former is 200 times lrger thn in the ltter. Incresing the number of MD trjectories by such fctor would imply time too long for prcticl computer simultions. In summry, forml limit of the Boltzmnn-Lorentz eqution for tgged prticle with mss lrge compred to prticles of the surrounding gs leds to Fokker-Plnck kinetic eqution. This eqution cn be mpped onto the corresponding eqution for elstic collisions by chnge of spce nd time scles whose exct solution is known. In prticulr, this eqution demonstrtes the seprtion of time scles ssocited with the rpid trnsition from complex initil trnsients to hydrodynmics diffusion. The property occurs independent of the degree of dissiption in the collisions between the hevy prticle nd surrounding gs prticles. In the present pper, the forml nlysis leding to the Fokker- Plnck eqution hs been verified by direct simultion of the Boltzmnn-Lorentz eqution. Finlly, the vlidity of the Boltzmnn-Lorentz eqution itself hs been confirmed by selected moleculr-dynmics simultions. ACKNOWLEDGMENTS The reserch of J.J.B., M.J.R-M., nd R.G-R. ws prtilly supported by the Dirección Generl de Investigción CientíficyTécnic Spin through the grnt PB The reserch of J.W.D ws supported in prt by NSF Grnt No. PHY APPENDIX: DETAILED FORM FOR The Fokker-Plnck limit leding to Eq. 2.8 is restricted to smll mss rtio nd 1. The prmeter is defined by Eq. 2.11, with the cooling rte nd friction coefficient e given by Eqs. 2.6 nd 2.9, respectively,

8 PRE 60 BROWNIAN MOTION IN A GRANULAR GAS 7181 t 2 e t 1 g 2 d/2 16 d3 2 1 g 2 d1 g 42 0 d1 g 0 dv dv 1 g 3 vv c* g. A1 In the third line the integrl hs been performed using n pproximte solution to the Boltzmnn eqution obtined by expnding (v) in Sonine polynomils nd retining the first correction to the Gussin 1,9. This pproximtion hs been tested vi Monte Crlo simultion nd is ccurte for the relevnt velocities within few percent 12. The contribution c*( g ) bove is due to this first correction nd is given by 321 g 12 2 g c* g 924d8d41 g 30 2 g 1 g. A2 The condition 1 restricts the rnge of g to tht for wek dissiption in the gs since must be smll. The detiled vlues cn be djusted by considering prticles of very different sizes. For simplicity, the sizes hve been tken equl in the simultions reported in Sec. III. Then, for the chosen vlues in the DSMC method in Sec. III, g 0.99 nd 10 2, the bove expression for simplifies to , A3 where we hve lso prticulrized for d3. Clerly 1 for ll 01 nd the Fokker-Plnck eqution is expected to hold for rbitrry inelsticity of the collisions of the tgged prticle with the gs prticles. In the reported MD simultions, we used hrd disks (d 2), 210 2, nd g This leds to , A4 nd, gin, the condition 1 is fulfilled for rbitrry vlues of. When crrying out the simultions, we hve observed tht the greement with the theoreticl predictions inproves s the vlue of decreses, lthough the vlidity of the theory requires just 1. 1 A. Goldshtein nd M. Shpiro, J. Fluid Mech. 282, N. Sel nd I. Goldhirsch, J. Fluid Mech. 361, J.J. Brey, J.W. Dufty, C-S. Kim, nd A. Sntos, Phys. Rev. E 58, V. Grzo nd J.W. Dufty, Phys. Rev. E 59, M-L. Tn nd I. Goldhirsch, Phys. Rev. Lett. 81, ; L.P. Kdnoff, Rev. Mod. Phys. 71, J.J. Brey, M.J. Ruiz-Montero, D. Cubero, nd R. Grci-Rojo, Phys. Fluids to be published. 7 J.J. Brey, J.W. Dufty, nd A. Sntos, J. Stt. Phys. 97, J.J. Brey, J.W. Dufty, nd A. Sntos, J. Stt. Phys. 87, T.P.C. vn Noije nd M.H. Ernst, Grnulr Mtter 1, J. A. McLennn, Introduction to Nonequilibrium Sttisticl Mechnics Prentice-Hll, Englewood Cliffs, NJ, G. Bird, Moleculr Gs Dynmics nd the Direct Simultion of Gs Flows Clrendon Press, Oxford, J.J. Brey, M.J. Ruiz-Montero, nd D. Cubero, Phys. Rev. E 54, M.P. Allen nd D.J. Tildesley, Computer Simultion of Liquids Oxford University Press, Oxford, 1987.

CHM Physical Chemistry I Chapter 1 - Supplementary Material

CHM Physical Chemistry I Chapter 1 - Supplementary Material CHM 3410 - Physicl Chemistry I Chpter 1 - Supplementry Mteril For review of some bsic concepts in mth, see Atkins "Mthemticl Bckground 1 (pp 59-6), nd "Mthemticl Bckground " (pp 109-111). 1. Derivtion

More information

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0)

P 3 (x) = f(0) + f (0)x + f (0) 2. x 2 + f (0) . In the problem set, you are asked to show, in general, the n th order term is a n = f (n) (0) 1 Tylor polynomils In Section 3.5, we discussed how to pproximte function f(x) round point in terms of its first derivtive f (x) evluted t, tht is using the liner pproximtion f() + f ()(x ). We clled this

More information

13: Diffusion in 2 Energy Groups

13: Diffusion in 2 Energy Groups 3: Diffusion in Energy Groups B. Rouben McMster University Course EP 4D3/6D3 Nucler Rector Anlysis (Rector Physics) 5 Sept.-Dec. 5 September Contents We study the diffusion eqution in two energy groups

More information

Monte Carlo method in solving numerical integration and differential equation

Monte Carlo method in solving numerical integration and differential equation Monte Crlo method in solving numericl integrtion nd differentil eqution Ye Jin Chemistry Deprtment Duke University yj66@duke.edu Abstrct: Monte Crlo method is commonly used in rel physics problem. The

More information

Lecture 14: Quadrature

Lecture 14: Quadrature Lecture 14: Qudrture This lecture is concerned with the evlution of integrls fx)dx 1) over finite intervl [, b] The integrnd fx) is ssumed to be rel-vlues nd smooth The pproximtion of n integrl by numericl

More information

Math 8 Winter 2015 Applications of Integration

Math 8 Winter 2015 Applications of Integration Mth 8 Winter 205 Applictions of Integrtion Here re few importnt pplictions of integrtion. The pplictions you my see on n exm in this course include only the Net Chnge Theorem (which is relly just the Fundmentl

More information

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM ROMAI J., v.9, no.2(2013), 173 179 THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM Alicj Piseck-Belkhyt, Ann Korczk Institute of Computtionl Mechnics nd Engineering,

More information

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by. NUMERICAL INTEGRATION 1 Introduction The inverse process to differentition in clculus is integrtion. Mthemticlly, integrtion is represented by f(x) dx which stnds for the integrl of the function f(x) with

More information

Review of Calculus, cont d

Review of Calculus, cont d Jim Lmbers MAT 460 Fll Semester 2009-10 Lecture 3 Notes These notes correspond to Section 1.1 in the text. Review of Clculus, cont d Riemnn Sums nd the Definite Integrl There re mny cses in which some

More information

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations ME 3600 Control ystems Chrcteristics of Open-Loop nd Closed-Loop ystems Importnt Control ystem Chrcteristics o ensitivity of system response to prmetric vritions cn be reduced o rnsient nd stedy-stte responses

More information

New Expansion and Infinite Series

New Expansion and Infinite Series Interntionl Mthemticl Forum, Vol. 9, 204, no. 22, 06-073 HIKARI Ltd, www.m-hikri.com http://dx.doi.org/0.2988/imf.204.4502 New Expnsion nd Infinite Series Diyun Zhng College of Computer Nnjing University

More information

Travelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing

Travelling Profile Solutions For Nonlinear Degenerate Parabolic Equation And Contour Enhancement In Image Processing Applied Mthemtics E-Notes 8(8) - c IN 67-5 Avilble free t mirror sites of http://www.mth.nthu.edu.tw/ men/ Trvelling Profile olutions For Nonliner Degenerte Prbolic Eqution And Contour Enhncement In Imge

More information

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

INTRODUCTION. The three general approaches to the solution of kinetics problems are: INTRODUCTION According to Newton s lw, prticle will ccelerte when it is subjected to unblnced forces. Kinetics is the study of the reltions between unblnced forces nd the resulting chnges in motion. The

More information

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions Physics 6C Solution of inhomogeneous ordinry differentil equtions using Green s functions Peter Young November 5, 29 Homogeneous Equtions We hve studied, especilly in long HW problem, second order liner

More information

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1 3 9. SEQUENCES AND SERIES 63. Representtion of functions s power series Consider power series x 2 + x 4 x 6 + x 8 + = ( ) n x 2n It is geometric series with q = x 2 nd therefore it converges for ll q =

More information

p-adic Egyptian Fractions

p-adic Egyptian Fractions p-adic Egyptin Frctions Contents 1 Introduction 1 2 Trditionl Egyptin Frctions nd Greedy Algorithm 2 3 Set-up 3 4 p-greedy Algorithm 5 5 p-egyptin Trditionl 10 6 Conclusion 1 Introduction An Egyptin frction

More information

221B Lecture Notes WKB Method

221B Lecture Notes WKB Method Clssicl Limit B Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using

More information

Numerical integration

Numerical integration 2 Numericl integrtion This is pge i Printer: Opque this 2. Introduction Numericl integrtion is problem tht is prt of mny problems in the economics nd econometrics literture. The orgniztion of this chpter

More information

Lecture 19: Continuous Least Squares Approximation

Lecture 19: Continuous Least Squares Approximation Lecture 19: Continuous Lest Squres Approximtion 33 Continuous lest squres pproximtion We begn 31 with the problem of pproximting some f C[, b] with polynomil p P n t the discrete points x, x 1,, x m for

More information

Part I: Basic Concepts of Thermodynamics

Part I: Basic Concepts of Thermodynamics Prt I: Bsic Concepts o Thermodynmics Lecture 4: Kinetic Theory o Gses Kinetic Theory or rel gses 4-1 Kinetic Theory or rel gses Recll tht or rel gses: (i The volume occupied by the molecules under ordinry

More information

5.7 Improper Integrals

5.7 Improper Integrals 458 pplictions of definite integrls 5.7 Improper Integrls In Section 5.4, we computed the work required to lift pylod of mss m from the surfce of moon of mss nd rdius R to height H bove the surfce of the

More information

4 The dynamical FRW universe

4 The dynamical FRW universe 4 The dynmicl FRW universe 4.1 The Einstein equtions Einstein s equtions G µν = T µν (7) relte the expnsion rte (t) to energy distribution in the universe. On the left hnd side is the Einstein tensor which

More information

Math 1B, lecture 4: Error bounds for numerical methods

Math 1B, lecture 4: Error bounds for numerical methods Mth B, lecture 4: Error bounds for numericl methods Nthn Pflueger 4 September 0 Introduction The five numericl methods descried in the previous lecture ll operte by the sme principle: they pproximte the

More information

Math& 152 Section Integration by Parts

Math& 152 Section Integration by Parts Mth& 5 Section 7. - Integrtion by Prts Integrtion by prts is rule tht trnsforms the integrl of the product of two functions into other (idelly simpler) integrls. Recll from Clculus I tht given two differentible

More information

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018

Physics 201 Lab 3: Measurement of Earth s local gravitational field I Data Acquisition and Preliminary Analysis Dr. Timothy C. Black Summer I, 2018 Physics 201 Lb 3: Mesurement of Erth s locl grvittionl field I Dt Acquisition nd Preliminry Anlysis Dr. Timothy C. Blck Summer I, 2018 Theoreticl Discussion Grvity is one of the four known fundmentl forces.

More information

1 1D heat and wave equations on a finite interval

1 1D heat and wave equations on a finite interval 1 1D het nd wve equtions on finite intervl In this section we consider generl method of seprtion of vribles nd its pplictions to solving het eqution nd wve eqution on finite intervl ( 1, 2. Since by trnsltion

More information

Lecture 20: Numerical Integration III

Lecture 20: Numerical Integration III cs4: introduction to numericl nlysis /8/0 Lecture 0: Numericl Integrtion III Instructor: Professor Amos Ron Scribes: Mrk Cowlishw, Yunpeng Li, Nthnel Fillmore For the lst few lectures we hve discussed

More information

8 Laplace s Method and Local Limit Theorems

8 Laplace s Method and Local Limit Theorems 8 Lplce s Method nd Locl Limit Theorems 8. Fourier Anlysis in Higher DImensions Most of the theorems of Fourier nlysis tht we hve proved hve nturl generliztions to higher dimensions, nd these cn be proved

More information

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

1.2. Linear Variable Coefficient Equations. y + b ! = a y + b  Remark: The case b = 0 and a non-constant can be solved with the same idea as above. 1 12 Liner Vrible Coefficient Equtions Section Objective(s): Review: Constnt Coefficient Equtions Solving Vrible Coefficient Equtions The Integrting Fctor Method The Bernoulli Eqution 121 Review: Constnt

More information

MAC-solutions of the nonexistent solutions of mathematical physics

MAC-solutions of the nonexistent solutions of mathematical physics Proceedings of the 4th WSEAS Interntionl Conference on Finite Differences - Finite Elements - Finite Volumes - Boundry Elements MAC-solutions of the nonexistent solutions of mthemticl physics IGO NEYGEBAUE

More information

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation

1.1. Linear Constant Coefficient Equations. Remark: A differential equation is an equation 1 1.1. Liner Constnt Coefficient Equtions Section Objective(s): Overview of Differentil Equtions. Liner Differentil Equtions. Solving Liner Differentil Equtions. The Initil Vlue Problem. 1.1.1. Overview

More information

Conservation Law. Chapter Goal. 5.2 Theory

Conservation Law. Chapter Goal. 5.2 Theory Chpter 5 Conservtion Lw 5.1 Gol Our long term gol is to understnd how mny mthemticl models re derived. We study how certin quntity chnges with time in given region (sptil domin). We first derive the very

More information

221A Lecture Notes WKB Method

221A Lecture Notes WKB Method A Lecture Notes WKB Method Hmilton Jcobi Eqution We strt from the Schrödinger eqution for single prticle in potentil i h t ψ x, t = [ ] h m + V x ψ x, t. We cn rewrite this eqution by using ψ x, t = e

More information

A SHORT NOTE ON THE MONOTONICITY OF THE ERLANG C FORMULA IN THE HALFIN-WHITT REGIME. Bernardo D Auria 1

A SHORT NOTE ON THE MONOTONICITY OF THE ERLANG C FORMULA IN THE HALFIN-WHITT REGIME. Bernardo D Auria 1 Working Pper 11-42 (31) Sttistics nd Econometrics Series December, 2011 Deprtmento de Estdístic Universidd Crlos III de Mdrid Clle Mdrid, 126 28903 Getfe (Spin) Fx (34) 91 624-98-49 A SHORT NOTE ON THE

More information

Review of basic calculus

Review of basic calculus Review of bsic clculus This brief review reclls some of the most importnt concepts, definitions, nd theorems from bsic clculus. It is not intended to tech bsic clculus from scrtch. If ny of the items below

More information

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a). The Fundmentl Theorems of Clculus Mth 4, Section 0, Spring 009 We now know enough bout definite integrls to give precise formultions of the Fundmentl Theorems of Clculus. We will lso look t some bsic emples

More information

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams Chpter 4 Contrvrince, Covrince, nd Spcetime Digrms 4. The Components of Vector in Skewed Coordintes We hve seen in Chpter 3; figure 3.9, tht in order to show inertil motion tht is consistent with the Lorentz

More information

The Active Universe. 1 Active Motion

The Active Universe. 1 Active Motion The Active Universe Alexnder Glück, Helmuth Hüffel, Sš Ilijić, Gerld Kelnhofer Fculty of Physics, University of Vienn helmuth.hueffel@univie.c.t Deprtment of Physics, FER, University of Zgreb ss.ilijic@fer.hr

More information

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı Mthemticl nd Computtionl Applictions, Vol., No., pp. 37-5, 5 http://dx.doi.org/.99/mc-5- SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD Mehmet Pkdemirli nd Gözde

More information

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as

Improper Integrals. Type I Improper Integrals How do we evaluate an integral such as Improper Integrls Two different types of integrls cn qulify s improper. The first type of improper integrl (which we will refer to s Type I) involves evluting n integrl over n infinite region. In the grph

More information

Sufficient condition on noise correlations for scalable quantum computing

Sufficient condition on noise correlations for scalable quantum computing Sufficient condition on noise correltions for sclble quntum computing John Presill, 2 Februry 202 Is quntum computing sclble? The ccurcy threshold theorem for quntum computtion estblishes tht sclbility

More information

Numerical Analysis: Trapezoidal and Simpson s Rule

Numerical Analysis: Trapezoidal and Simpson s Rule nd Simpson s Mthemticl question we re interested in numericlly nswering How to we evlute I = f (x) dx? Clculus tells us tht if F(x) is the ntiderivtive of function f (x) on the intervl [, b], then I =

More information

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS.

THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS. THE EXISTENCE-UNIQUENESS THEOREM FOR FIRST-ORDER DIFFERENTIAL EQUATIONS RADON ROSBOROUGH https://intuitiveexplntionscom/picrd-lindelof-theorem/ This document is proof of the existence-uniqueness theorem

More information

Numerical Integration

Numerical Integration Chpter 5 Numericl Integrtion Numericl integrtion is the study of how the numericl vlue of n integrl cn be found. Methods of function pproximtion discussed in Chpter??, i.e., function pproximtion vi the

More information

Session 13

Session 13 780.20 Session 3 (lst revised: Februry 25, 202) 3 3. 780.20 Session 3. Follow-ups to Session 2 Histogrms of Uniform Rndom Number Distributions. Here is typicl figure you might get when histogrmming uniform

More information

Theoretical foundations of Gaussian quadrature

Theoretical foundations of Gaussian quadrature Theoreticl foundtions of Gussin qudrture 1 Inner product vector spce Definition 1. A vector spce (or liner spce) is set V = {u, v, w,...} in which the following two opertions re defined: (A) Addition of

More information

Ordinary differential equations

Ordinary differential equations Ordinry differentil equtions Introduction to Synthetic Biology E Nvrro A Montgud P Fernndez de Cordob JF Urchueguí Overview Introduction-Modelling Bsic concepts to understnd n ODE. Description nd properties

More information

Name Solutions to Test 3 November 8, 2017

Name Solutions to Test 3 November 8, 2017 Nme Solutions to Test 3 November 8, 07 This test consists of three prts. Plese note tht in prts II nd III, you cn skip one question of those offered. Some possibly useful formuls cn be found below. Brrier

More information

Week 10: Line Integrals

Week 10: Line Integrals Week 10: Line Integrls Introduction In this finl week we return to prmetrised curves nd consider integrtion long such curves. We lredy sw this in Week 2 when we integrted long curve to find its length.

More information

+ x 2 dω 2 = c 2 dt 2 +a(t) [ 2 dr 2 + S 1 κx 2 /R0

+ x 2 dω 2 = c 2 dt 2 +a(t) [ 2 dr 2 + S 1 κx 2 /R0 Notes for Cosmology course, fll 2005 Cosmic Dynmics Prelude [ ds 2 = c 2 dt 2 +(t) 2 dx 2 ] + x 2 dω 2 = c 2 dt 2 +(t) [ 2 dr 2 + S 1 κx 2 /R0 2 κ (r) 2 dω 2] nd x = S κ (r) = r, R 0 sin(r/r 0 ), R 0 sinh(r/r

More information

Rel Gses 1. Gses (N, CO ) which don t obey gs lws or gs eqution P=RT t ll pressure nd tempertures re clled rel gses.. Rel gses obey gs lws t extremely low pressure nd high temperture. Rel gses devited

More information

( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that

( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that Arc Length of Curves in Three Dimensionl Spce If the vector function r(t) f(t) i + g(t) j + h(t) k trces out the curve C s t vries, we cn mesure distnces long C using formul nerly identicl to one tht we

More information

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral

f(x) dx, If one of these two conditions is not met, we call the integral improper. Our usual definition for the value for the definite integral Improper Integrls Every time tht we hve evluted definite integrl such s f(x) dx, we hve mde two implicit ssumptions bout the integrl:. The intervl [, b] is finite, nd. f(x) is continuous on [, b]. If one

More information

Review of Gaussian Quadrature method

Review of Gaussian Quadrature method Review of Gussin Qudrture method Nsser M. Asi Spring 006 compiled on Sundy Decemer 1, 017 t 09:1 PM 1 The prolem To find numericl vlue for the integrl of rel vlued function of rel vrile over specific rnge

More information

SUMMER KNOWHOW STUDY AND LEARNING CENTRE

SUMMER KNOWHOW STUDY AND LEARNING CENTRE SUMMER KNOWHOW STUDY AND LEARNING CENTRE Indices & Logrithms 2 Contents Indices.2 Frctionl Indices.4 Logrithms 6 Exponentil equtions. Simplifying Surds 13 Opertions on Surds..16 Scientific Nottion..18

More information

Math 124A October 04, 2011

Math 124A October 04, 2011 Mth 4A October 04, 0 Viktor Grigoryn 4 Vibrtions nd het flow In this lecture we will derive the wve nd het equtions from physicl principles. These re second order constnt coefficient liner PEs, which model

More information

Forces from Strings Under Tension A string under tension medites force: the mgnitude of the force from section of string is the tension T nd the direc

Forces from Strings Under Tension A string under tension medites force: the mgnitude of the force from section of string is the tension T nd the direc Physics 170 Summry of Results from Lecture Kinemticl Vribles The position vector ~r(t) cn be resolved into its Crtesin components: ~r(t) =x(t)^i + y(t)^j + z(t)^k. Rtes of Chnge Velocity ~v(t) = d~r(t)=

More information

Applicable Analysis and Discrete Mathematics available online at

Applicable Analysis and Discrete Mathematics available online at Applicble Anlysis nd Discrete Mthemtics vilble online t http://pefmth.etf.rs Appl. Anl. Discrete Mth. 4 (2010), 23 31. doi:10.2298/aadm100201012k NUMERICAL ANALYSIS MEETS NUMBER THEORY: USING ROOTFINDING

More information

10 Vector Integral Calculus

10 Vector Integral Calculus Vector Integrl lculus Vector integrl clculus extends integrls s known from clculus to integrls over curves ("line integrls"), surfces ("surfce integrls") nd solids ("volume integrls"). These integrls hve

More information

Lecture 3 Gaussian Probability Distribution

Lecture 3 Gaussian Probability Distribution Introduction Lecture 3 Gussin Probbility Distribution Gussin probbility distribution is perhps the most used distribution in ll of science. lso clled bell shped curve or norml distribution Unlike the binomil

More information

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007

A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H. Thomas Shores Department of Mathematics University of Nebraska Spring 2007 A REVIEW OF CALCULUS CONCEPTS FOR JDEP 384H Thoms Shores Deprtment of Mthemtics University of Nebrsk Spring 2007 Contents Rtes of Chnge nd Derivtives 1 Dierentils 4 Are nd Integrls 5 Multivrite Clculus

More information

CBE 291b - Computation And Optimization For Engineers

CBE 291b - Computation And Optimization For Engineers The University of Western Ontrio Fculty of Engineering Science Deprtment of Chemicl nd Biochemicl Engineering CBE 9b - Computtion And Optimiztion For Engineers Mtlb Project Introduction Prof. A. Jutn Jn

More information

and that at t = 0 the object is at position 5. Find the position of the object at t = 2.

and that at t = 0 the object is at position 5. Find the position of the object at t = 2. 7.2 The Fundmentl Theorem of Clculus 49 re mny, mny problems tht pper much different on the surfce but tht turn out to be the sme s these problems, in the sense tht when we try to pproimte solutions we

More information

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature CMDA 4604: Intermedite Topics in Mthemticl Modeling Lecture 19: Interpoltion nd Qudrture In this lecture we mke brief diversion into the res of interpoltion nd qudrture. Given function f C[, b], we sy

More information

Entropy ISSN

Entropy ISSN Entropy 006, 8[], 50-6 50 Entropy ISSN 099-4300 www.mdpi.org/entropy/ ENTROPY GENERATION IN PRESSURE GRADIENT ASSISTED COUETTE FLOW WITH DIFFERENT THERMAL BOUNDARY CONDITIONS Abdul Aziz Deprtment of Mechnicl

More information

Continuous Random Variables

Continuous Random Variables STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 217 Néhémy Lim Continuous Rndom Vribles Nottion. The indictor function of set S is rel-vlued function defined by : { 1 if x S 1 S (x) if x S Suppose tht

More information

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), )

Euler, Ioachimescu and the trapezium rule. G.J.O. Jameson (Math. Gazette 96 (2012), ) Euler, Iochimescu nd the trpezium rule G.J.O. Jmeson (Mth. Gzette 96 (0), 36 4) The following results were estblished in recent Gzette rticle [, Theorems, 3, 4]. Given > 0 nd 0 < s

More information

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION XX IMEKO World Congress Metrology for Green Growth September 9,, Busn, Republic of Kore THERMAL EXPANSION COEFFICIENT OF WATER FOR OLUMETRIC CALIBRATION Nieves Medin Hed of Mss Division, CEM, Spin, mnmedin@mityc.es

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 2013 Outline 1 Riemnn Sums 2 Riemnn Integrls 3 Properties

More information

Chapter 5 : Continuous Random Variables

Chapter 5 : Continuous Random Variables STAT/MATH 395 A - PROBABILITY II UW Winter Qurter 216 Néhémy Lim Chpter 5 : Continuous Rndom Vribles Nottions. N {, 1, 2,...}, set of nturl numbers (i.e. ll nonnegtive integers); N {1, 2,...}, set of ll

More information

Physics 202H - Introductory Quantum Physics I Homework #08 - Solutions Fall 2004 Due 5:01 PM, Monday 2004/11/15

Physics 202H - Introductory Quantum Physics I Homework #08 - Solutions Fall 2004 Due 5:01 PM, Monday 2004/11/15 Physics H - Introductory Quntum Physics I Homework #8 - Solutions Fll 4 Due 5:1 PM, Mondy 4/11/15 [55 points totl] Journl questions. Briefly shre your thoughts on the following questions: Of the mteril

More information

Reversals of Signal-Posterior Monotonicity for Any Bounded Prior

Reversals of Signal-Posterior Monotonicity for Any Bounded Prior Reversls of Signl-Posterior Monotonicity for Any Bounded Prior Christopher P. Chmbers Pul J. Hely Abstrct Pul Milgrom (The Bell Journl of Economics, 12(2): 380 391) showed tht if the strict monotone likelihood

More information

MATH 144: Business Calculus Final Review

MATH 144: Business Calculus Final Review MATH 144: Business Clculus Finl Review 1 Skills 1. Clculte severl limits. 2. Find verticl nd horizontl symptotes for given rtionl function. 3. Clculte derivtive by definition. 4. Clculte severl derivtives

More information

On the Uncertainty of Sensors Based on Magnetic Effects. E. Hristoforou, E. Kayafas, A. Ktena, DM Kepaptsoglou

On the Uncertainty of Sensors Based on Magnetic Effects. E. Hristoforou, E. Kayafas, A. Ktena, DM Kepaptsoglou On the Uncertinty of Sensors Bsed on Mgnetic Effects E. ristoforou, E. Kyfs, A. Kten, DM Kepptsoglou Ntionl Technicl University of Athens, Zogrfou Cmpus, Athens 1578, Greece Tel: +3177178, Fx: +3177119,

More information

5.3 Nonlinear stability of Rayleigh-Bénard convection

5.3 Nonlinear stability of Rayleigh-Bénard convection 118 5.3 Nonliner stbility of Ryleigh-Bénrd convection In Chpter 1, we sw tht liner stbility only tells us whether system is stble or unstble to infinitesimlly-smll perturbtions, nd tht there re cses in

More information

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac REVIEW OF ALGEBRA Here we review the bsic rules nd procedures of lgebr tht you need to know in order to be successful in clculus. ARITHMETIC OPERATIONS The rel numbers hve the following properties: b b

More information

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes Jim Lmbers MAT 169 Fll Semester 2009-10 Lecture 4 Notes These notes correspond to Section 8.2 in the text. Series Wht is Series? An infinte series, usully referred to simply s series, is n sum of ll of

More information

The Wave Equation I. MA 436 Kurt Bryan

The Wave Equation I. MA 436 Kurt Bryan 1 Introduction The Wve Eqution I MA 436 Kurt Bryn Consider string stretching long the x xis, of indeterminte (or even infinite!) length. We wnt to derive n eqution which models the motion of the string

More information

Riemann Sums and Riemann Integrals

Riemann Sums and Riemann Integrals Riemnn Sums nd Riemnn Integrls Jmes K. Peterson Deprtment of Biologicl Sciences nd Deprtment of Mthemticl Sciences Clemson University August 26, 203 Outline Riemnn Sums Riemnn Integrls Properties Abstrct

More information

Chapter 0. What is the Lebesgue integral about?

Chapter 0. What is the Lebesgue integral about? Chpter 0. Wht is the Lebesgue integrl bout? The pln is to hve tutoril sheet ech week, most often on Fridy, (to be done during the clss) where you will try to get used to the ides introduced in the previous

More information

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1 MATH34032: Green s Functions, Integrl Equtions nd the Clculus of Vritions 1 Section 1 Function spces nd opertors Here we gives some brief detils nd definitions, prticulrly relting to opertors. For further

More information

20 MATHEMATICS POLYNOMIALS

20 MATHEMATICS POLYNOMIALS 0 MATHEMATICS POLYNOMIALS.1 Introduction In Clss IX, you hve studied polynomils in one vrible nd their degrees. Recll tht if p(x) is polynomil in x, the highest power of x in p(x) is clled the degree of

More information

The steps of the hypothesis test

The steps of the hypothesis test ttisticl Methods I (EXT 7005) Pge 78 Mosquito species Time of dy A B C Mid morning 0.0088 5.4900 5.5000 Mid Afternoon.3400 0.0300 0.8700 Dusk 0.600 5.400 3.000 The Chi squre test sttistic is the sum of

More information

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0. STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA STEPHEN SCHECTER. The unit step function nd piecewise continuous functions The Heviside unit step function u(t) is given by if t

More information

Intro to Nuclear and Particle Physics (5110)

Intro to Nuclear and Particle Physics (5110) Intro to Nucler nd Prticle Physics (5110) Feb, 009 The Nucler Mss Spectrum The Liquid Drop Model //009 1 E(MeV) n n(n-1)/ E/[ n(n-1)/] (MeV/pir) 1 C 16 O 0 Ne 4 Mg 7.7 14.44 19.17 8.48 4 5 6 6 10 15.4.41

More information

Massachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 6

Massachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 6 Msschusetts Institute of Technology Quntum Mechnics I (8.) Spring 5 Solutions to Problem Set 6 By Kit Mtn. Prctice with delt functions ( points) The Dirc delt function my be defined s such tht () (b) 3

More information

ftp.fe?a:fmmmhm Quickly get policy ) long generally equilibrium independent steady P # steady E amp= : Dog steady systems Every equilibrium by density

ftp.fe?a:fmmmhm Quickly get policy ) long generally equilibrium independent steady P # steady E amp= : Dog steady systems Every equilibrium by density Introduction fle SI 3 61 Non equilibrium Sttisticl Mechnics 2520 1 Techer : Tens It Brdrsonbrdrson@kthse Office 174 : 1049 ( open door policy Lecture 1 & The second lw nd origin of irreversibility Wht

More information

Scientific notation is a way of expressing really big numbers or really small numbers.

Scientific notation is a way of expressing really big numbers or really small numbers. Scientific Nottion (Stndrd form) Scientific nottion is wy of expressing relly big numbers or relly smll numbers. It is most often used in scientific clcultions where the nlysis must be very precise. Scientific

More information

Exam 1 Solutions (1) C, D, A, B (2) C, A, D, B (3) C, B, D, A (4) A, C, D, B (5) D, C, A, B

Exam 1 Solutions (1) C, D, A, B (2) C, A, D, B (3) C, B, D, A (4) A, C, D, B (5) D, C, A, B PHY 249, Fll 216 Exm 1 Solutions nswer 1 is correct for ll problems. 1. Two uniformly chrged spheres, nd B, re plced t lrge distnce from ech other, with their centers on the x xis. The chrge on sphere

More information

Section 6.1 INTRO to LAPLACE TRANSFORMS

Section 6.1 INTRO to LAPLACE TRANSFORMS Section 6. INTRO to LAPLACE TRANSFORMS Key terms: Improper Integrl; diverge, converge A A f(t)dt lim f(t)dt Piecewise Continuous Function; jump discontinuity Function of Exponentil Order Lplce Trnsform

More information

Measuring Electron Work Function in Metal

Measuring Electron Work Function in Metal n experiment of the Electron topic Mesuring Electron Work Function in Metl Instructor: 梁生 Office: 7-318 Emil: shling@bjtu.edu.cn Purposes 1. To understnd the concept of electron work function in metl nd

More information

Discrete Mathematics and Probability Theory Summer 2014 James Cook Note 17

Discrete Mathematics and Probability Theory Summer 2014 James Cook Note 17 CS 70 Discrete Mthemtics nd Proility Theory Summer 2014 Jmes Cook Note 17 I.I.D. Rndom Vriles Estimting the is of coin Question: We wnt to estimte the proportion p of Democrts in the US popultion, y tking

More information

Math 426: Probability Final Exam Practice

Math 426: Probability Final Exam Practice Mth 46: Probbility Finl Exm Prctice. Computtionl problems 4. Let T k (n) denote the number of prtitions of the set {,..., n} into k nonempty subsets, where k n. Argue tht T k (n) kt k (n ) + T k (n ) by

More information

Section 14.3 Arc Length and Curvature

Section 14.3 Arc Length and Curvature Section 4.3 Arc Length nd Curvture Clculus on Curves in Spce In this section, we ly the foundtions for describing the movement of n object in spce.. Vector Function Bsics In Clc, formul for rc length in

More information

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives Block #6: Properties of Integrls, Indefinite Integrls Gols: Definition of the Definite Integrl Integrl Clcultions using Antiderivtives Properties of Integrls The Indefinite Integrl 1 Riemnn Sums - 1 Riemnn

More information

Series: Elementary, then Taylor & MacLaurin

Series: Elementary, then Taylor & MacLaurin Lecture 3 Series: Elementry, then Tylor & McLurin This lecture on series strts in strnge plce, by revising rithmetic nd geometric series, nd then proceeding to consider other elementry series. The relevnce

More information

F k B T = 1 2 K. θ q e iq x, (2) q 2 θq. d 2 1

F k B T = 1 2 K. θ q e iq x, (2) q 2 θq. d 2 1 Physics 7c: Sttisticl Mechnics Kosterlitz-Thouless Trnsition The trnsition to superfluidity in thin films is n exmple of the Kosterlitz-Thouless trnsition, n exotic new type of phse trnsition driven by

More information

Stuff You Need to Know From Calculus

Stuff You Need to Know From Calculus Stuff You Need to Know From Clculus For the first time in the semester, the stuff we re doing is finlly going to look like clculus (with vector slnt, of course). This mens tht in order to succeed, you

More information

First midterm topics Second midterm topics End of quarter topics. Math 3B Review. Steve. 18 March 2009

First midterm topics Second midterm topics End of quarter topics. Math 3B Review. Steve. 18 March 2009 Mth 3B Review Steve 18 Mrch 2009 About the finl Fridy Mrch 20, 3pm-6pm, Lkretz 110 No notes, no book, no clcultor Ten questions Five review questions (Chpters 6,7,8) Five new questions (Chpters 9,10) No

More information

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1

Exam 2, Mathematics 4701, Section ETY6 6:05 pm 7:40 pm, March 31, 2016, IH-1105 Instructor: Attila Máté 1 Exm, Mthemtics 471, Section ETY6 6:5 pm 7:4 pm, Mrch 1, 16, IH-115 Instructor: Attil Máté 1 17 copies 1. ) Stte the usul sufficient condition for the fixed-point itertion to converge when solving the eqution

More information