Dense fluid transport for inelastic hard spheres
|
|
- Iris Dixon
- 5 years ago
- Views:
Transcription
1 PHYSICAL REVIEW E VOLUME 59, NUMBER 5 MAY 999 Dense fluid transport for inelastic hard spheres V. Garzó* and J. W. Dufty Department of Physics, University of Florida, Gainesville, Florida 36 Received December 998 The revised Enskog theory for inelastic hard spheres is considered as a model for rapid flow granular media at finite densities. A normal solution is obtained via the Chapman-Enskog method for states near the local homogeneous cooling state. The analysis is performed to first order in the spatial gradients, allowing identification of the Navier-Stokes order transport coefficients associated with the heat and momentum fluxes. In addition, the cooling rate is calculated to first order in the gradients and expressed in terms of the transport coefficients. The transport coefficients are determined from linear integral equations analogous to those for elastic collisions. The solubility conditions for these equations are confirmed and the transport coefficients are calculated as explicit functions of the density and restitution coefficient using a Sonine polynomial expansion. The results are not limited to small dissipation. Finally, the analysis is repeated using a simpler kinetic model. Excellent agreement is obtained with the results from the revised Enskog equation. S063-65X PACS numbers: 8.05.Rm, 05.0.Dd, 5.0.y, 47.0.k I. INTRODUCTION The analysis of granular media under rapid flow frequently exploits a hydrodynamic description. The essential difference from that for normal fluids is the absence of energy conservation, leading to both obvious and subtle modifications of the usual Navier-Stokes equations for states with small spatial gradients of the hydrodynamic fields. The justification for a hydrodynamic description and a detailed derivation of the form of the transport coefficients remains a topic of interest and controversy. For a low density system these problems have been addressed using the Boltzmann kinetic equation modified to account for inelastic binary collisions. Recently, the Chapman-Enskog solution to the Boltzmann equation, familiar from transport theory for normal gases, has been extended to the inelastic case to obtain the Navier-Stokes equations and detailed expressions for the transport coefficients as functions of the restitution coefficient. The case considered was an idealized gas of smooth, spherical hard spheres with inelastic binary collisions. Preliminary comparisons with direct Monte Carlo simulation of the Boltzmann equation indicate the results are quite accurate, even for strong dissipation. The objective here is to extend this analysis to the revised Enskog kinetic theory RET for a description of hydrodynamics and transport at higher densities. The RET for elastic collisions 3 is known to be an accurate kinetic theory over the entire fluid domain. In fact it describes the crystal phase as well, which may suggest its relevance for granular media undergoing cluster formation. Its generalization to inelastic collisions is straightforward see, for example, Ref. 4 and the Chapman-Enskog method can be applied to obtain the Navier-Stokes hydrodynamic equations and the associated transport coefficients. Similar attempts to apply the Boltzmann and Enskog kinetic equations to derive transport coefficients began fifteen years *Permanent address: Departamento de Física, Universidad de Extremadura, 0607 Badajoz, Spain. ago 5 8, but the technical difficulties of the analysis entailed approximations that limited their accuracy. The analysis here provides formally exact results for the distribution function and transport coefficients in terms of the solutions to linear integral equations. In this respect it parallels closely the analysis for elastic collisions. As in the latter case, a first Sonine polynomial approximation is used to obtain practical results from this formulation. The approximation technique is known to be accurate for elastic collisions and is expected to have comparable accuracy here as well. The results apply to all degrees of dissipation and over the entire fluid density range. Consequently, it provides the basis for an unambiguous basis for the application of hydrodynamics to granular flow under realistic conditions. There are several motivations for this study. First, a quantitative test of the hydrodynamic description within the context of kinetic theory is possible by comparing predicted transport properties with those obtained directly from the RET via Monte Carlo simulation. Such simulation methods have been developed and applied in recent years 9 and the results obtained here provide the basis for a test of the standard method for derivation of hydrodynamics. Second, by extending the Boltzmann analysis to high densities comparisons with molecular dynamics simulations become practical. For example, a comparison of the calculated viscosity dependence on both the density and restitution coefficient with that from molecular dynamics would determine the validity and any limitations of the kinetic theory and the hydrodynamic descriptions for granular flow. Such a test is essential to address the frequently made speculation that a hydrodynamic and/or kinetic theory description of granular flow is limited to weak dissipation. None of the analysis here supports this latter view, beyond complications of the well-known hydrodynamic instability of the reference cooling state. Finally, accurate predictions from the RET allow tests of simpler kinetic models that have been proposed 0, for application to more complex boundary driven states. In summary, the results presented here provide the basis for both practical quantitative applications and fundamental benchmark tests. The RET and associated macroscopic conservation laws X/99/595/58957/$5.00 PRE The American Physical Society
2 5896 V. GARZÓ AND J.W. DUFTY PRE 59 are reviewed in the next section. An overview of the Chapman-Enskog method for solving this equation is given in Sec. III. The results for the distribution function, heat flux, and momentum flux to first order in the spatial gradients are provided, with the details of the calculation appearing in a series of Appendixes. As in the case of elastic collisions these properties are expressed in terms of solutions to a set of linear inhomogeneous integral equations. It is shown explicitly that the necessary solubility conditions are satisfied. In Sec. IV these equations are solved using a lowest order expansion in Sonine polynomials. Explicit expressions for the distribution function and transport coefficients are given and compared to the elastic limit as a function of density and restitution coefficient. In Sec. V the corresponding results from a simple kinetic model are obtained and compared with those from the RET. Excellent agreement is found for a wide range of densities and restitution coefficients. This analysis is summarized and discussed in the last section. II. RET AND CONSERVATION LAWS In this section a brief review of the RET kinetic equation and associated conservation laws are presented. The oneparticle distribution function f (r,v,t) is determined from 4 t v fr,v,t,j E r,v ft, where J E is the Enskog collision operator J E r,v f t dv dˆ ˆ gˆ g f () r,r,v,v,t f () r,r,v,v,t, f () r,r,v,v,t r,r nt fr,v,t,fr,v,t. 3 Here, is the hard sphere diameter, ˆ, ˆ being a unit vector, is the Heaviside step function, and gv v. The primes on the velocities denote the initial values v,v that lead to v,v following a binary collision, v v ( )(ˆ g)ˆ, v v ( )(ˆ g)ˆ. The parameter is the coefficient of restitution with 0. Finally, r,rn(t) is the equilibrium pair correlation function at contact as a functional of the nonequilibrium density field n(r,t) defined by n(r,t)dvf (r,v,t). The macroscopic balance equations for density n, momentum mu, and energy 3 nk B T follow directly from Eq. by multiplying with, mv, and mv / and integrating over v, D t nn U0, D t U i mn j P ij 0, D t T 3nk B P ij j U i qt, where D t t U is the material derivative, m is the mass of a particle, k B is the Boltzmann s constant, T(r,t) is the temperature, and U(r,t) is the flow velocity. The cooling rate is proportional to ( ) and is due to dissipative collisions. The pressure tensor P(r,t) and the heat flux q(r,t) have both kinetic and collisional transfer contributions, i.e., PP k P c and qq k q c. The kinetic contributions are given by P k dv mvvf r,v,t, q k dv mv Vfr,v,t, and the collisional transfer contributions are 4 P c m 3 dv 4 dv dˆ ˆ g 7 ˆ g ˆ ˆ 0 d f () r,r,v,v,t, q c m 3 dv 4 dv dˆ ˆ gˆ g G ˆ ˆ 0 df () r,r,v,v,t. Here, VvU(r,t) is the peculiar velocity and G (V V ) is the velocity of center of mass. Finally, the cooling rate in Eq. 6 is given by m n dv dv dˆ ˆ g ˆ g 3 f () r,r,v,v,t, with /(k B T). The macroscopic balance equations 4 6 are exact consequences of the RET and provide the basis for developing a hydrodynamic description. The latter follows from a representation of P(r,t), q(r,t), and (r,t) as explicit functionals of the hydrodynamic fields and their gradients. The resulting closed set of equations constitutes the hydrodynamics for the RET. Since P(r,t), q(r,t), and (r,t) are given above as explicit functionals of f (r,v,t) such a representation can be obtained if a corresponding solution to the RET can be obtained as a function of the fields and their gradients. This is done in the next section using the Chapman-Enskog method. III. CHAPMAN-ENSKOG EXPANSION The Chapman-Enskog method assumes the existence of a normal solution such that all space and time dependence of the distribution function occurs through the hydrodynamic fields f r,v,t f vnt,tt,ut. The notation on the right side indicates a functional dependence on the density, temperature, and flow velocity. When this form is substituted into the RET all space and time de-
3 PRE 59 DENSE FLUID TRANSPORT FOR INELASTIC HARD SPHERES 5897 rivatives occur through the corresponding derivatives of the hydrodynamic fields. The time derivatives can be expressed in terms of space derivatives using the macroscopic balance equations 4 6. In this way the RET becomes an equation to determine the functional dependence of f vn(t),t(t),u(t) on the fields which are then determined self-consistently from the macroscopic balance equations. This procedure becomes more practical for states with small spatial gradients. Then the functional dependence can be made local by an expansion of the fields in terms of their gradients, with the distribution function determined perturbatively. To implement this expansion it is convenient to introduce a formal expansion parameter referred to as the uniformity parameter. It measures the dependence on spatial gradients e.g., a term of order is of first order in a hydrodynamic gradient, is either a product of two first order hydrodynamic gradients or one second degree hydrodynamic gradient. Accordingly the distribution function, collision operator, and time derivative are given by the representations f f (0) f (), t t (0) t (). J E J (0) J (), The coefficients in the time derivative expansion are identified by a representation of the fluxes and the cooling rate in the macroscopic balance equations as a similar series through their definitions as functionals of f. To zeroth order in the RET becomes t (0) f (0) J (0) f (0),f (0), J (0) f (0),f (0) dv dˆ ˆ gˆ g 3 f (0) v f (0) v f (0) v f (0) v. 4 Here r,rn(t) nn(r,t) is the pair functional evaluated with all density fields at the local point r. The collision operator 4 can be recognized as the Boltzmann operator for inelastic collisions multiplied by this factor. The macroscopic balance equations to this order read t (0) n0, t (0) U0, T t (0) T (0), 5 where the cooling rate (0) is determined by Eq. 0 to zeroth order (0) m 4n dv dv f (0) v f (0) v v v 3. 6 Since f (0) has the form the time derivative can be evaluated as f (0) nv 0 3 V/v 0, v 0 k B Tr,t/m. 8 The dependence on the magnitude of V follows from the requirement that to zeroth order in the gradients the distribution function must be isotropic with respect to the peculiar velocity. The RET to this order is therefore the solution to (0) v Vf (0) J (0) f (0),f (0). 9 To date an exact solution has not been found. However, a very good approximation can be obtained from an expansion in Sonine polynomials 3. In all of the following it is presumed that f (0) is known. Since the distribution function is isotropic the zeroth order pressure tensor and heat flux are found from Eqs. 7 0 to be where p is the hydrostatic pressure P (0) ij p ij, q (0) 0, 0 pnk B T m3 dv dv f (0) V f (0) V dˆ ˆ gˆ g nk B T n*. 3 The dimensionless density n*n 3 has been introduced. The analysis is similar to first order in, although technically more complex. Only the results are presented in this section with the details given in the appendixes. The RET equation to first order is found to be t (0) Lf () (,) f () Q v Vf (0) AV ln TBV ln nc ij V j U i i U j 3 ij UDV U, where L is the linearized collision operator LXQJ (0) f (0),XJ (0) X,f (0). 3 Here QP, where P is the projection operator onto the space spanned by, V, and V see Eq. A3 of Appendix A. The velocity dependence on the right side of Eq. is AVn,U,T 5 f (0) V vf (0)V t (0) f (0) (0) T T f (0) (0) v Vf (0), 7 p p k Vf(0) m v f (0) where v /v. The second equality follows from dimensional analysis which requires that the temperature dependence of f (0) must occur in the form K v Vf (0), 4
4 5898 V. GARZÓ AND J.W. DUFTY PRE 59 BVn,U,T p p kn n ln p p k Vf (0) m v f (0) n n ln K f (0), 5 P A B C ij D0. 36 C ij Vn,U,TV i vj f (0) K i vj f (0), 6 DV 3 K i vi f (0) ppk 3Q 3p k v Vf (0). 7 Here, vi /v i, p k nk B T and the linear operator KX is defined by KXQ dv dˆ ˆ g ˆ g f (0) v Xv f (0) v Xv. 8 The functional (,) ( f () ) is related to () by see, Eq. A3 () ppk p k U (,) f (), 9 (,) f () m n dv dv f (0) V f () V g 3. The solution to Eq. is of the form f () A ln TB ln n 30 C ij j U i i U j 3 ij UD U. 3 Substitution of this into Eq. gives the integral equations to determine A(V), B(V), C ij (V), and D(V) (0) T T L (0) AA, 3 (0) T T LBBn n ln (0) A, (0) T T LC ij C ij, (0) T T LD (,) DQ v Vf (0) D. 35 Use has been made of the fact that (,) ( f () ) is a scalar and therefore couples only to D. The lowest order distribution f (0) is a functional of the exact nonequilibrium fields. Consequently, its moments with respect to, V, and V are the same as those for the exact distribution, i.e., Pf (0) Pf. This implies that f () lies entirely in the orthogonal subspace so acceptable solutions to Eqs must obey There are additional necessary conditions for the solution to these equations to exist, the Fredholm alternative 4. These solubility conditions are proved in Appendix B. With the distribution function f () determined, the heat and momentum fluxes can be calculated to first order in the spatial gradients from Eqs The details of the calculations are given in Appendixes C and D with the results P () ij j U i i U j 3 ij U ij U, q () T n, 37 where is the shear viscosity, is the bulk viscosity, is the thermal conductivity, and is an additional transport coefficient not present for elastic collisions k n* 5 3, 38 5 mn*nv 0 I 8, 39 k n* 5 mn*nv 0 3 I 48T, 40 k n* 5. 4 Here the superscript k denotes the contributions from the kinetic parts of the fluxes. The kinetic parts of the transport coefficients in these results are expressed in terms of the solution to the integral equations as k 0 dv D ij VC ij V, k 0, 4 where D ij (V) and S(V) are D ij VmV i V j 3 V ij, k 3T dvsv AV, k 3n dvsv BV, 43 SVV mv 5 k BT. 44 The dimensionless integrals I () and I () in Eqs. 39 and 40 are
5 PRE 59 DENSE FLUID TRANSPORT FOR INELASTIC HARD SPHERES 5899 I n v 0 dv dv f (0) V f (0) V g, I n v 0 3 dv dv f (0) V f (0) V g g G gg 3 gg G 4 g In addition to the fluxes, the cooling rate can be calculated to first order in the gradients from Eqs. 9 and 30 () ppk p k f (0) V DV g 3 U. m 4n dv dv 47 This completes the results of the Chapman-Enskog solution to first order in the spatial gradients. The distribution function f f (0) f () with f (0) given by the solution to Eq. 9 and f () given by Eq. 3 is exact to this order in the gradients, as are the fluxes 37. These results require solution to the linear equations An approximate solution, known to be accurate for elastic collisions, is given in the next section. For that purpose it is convenient to use an equivalent exact expression for the kinetic parts of the transport coefficients see Appendix C k n (0) 5 3n*, 48 k 5nk B m (0) dvd ij VLC ij V, dvd ij VC ij V dvsv LAV dvsv LBV,. dvsv AV dvsv BV 5 Further, 0 5k B 0 /4m and 0 5m / /6 () / are the low density values of the thermal conductivity and the shear viscosity in the elastic limit, respectively. The constant c*() is related to the fourth moment of f (0) by c* 5 8 m dvf n (0) VV It is seen that c*() vanishes if f (0) is replaced by the Maxwellian distribution, and hence it is a measure of the deviation of the reference state from that for a gas with elastic collisions. IV. SONINE POLYNOMIAL APPROXIMATION More explicit results require determination of f (0) and the solutions A, B, and C ij to the linear integral equations It is useful to represent these quantities as an expansion in a complete set of polynomials and generate approximations by truncating the expansion. In practice the leading term in these expansions provides a very accurate description over the full range of dissipation and density. The determination of f (0) to leading order in the Sonine polynomial expansion has been discussed elsewhere 3 so only the result is quoted here f (0) Vf M V 4 c*ev, 53 p p kc* 0 n* 3 c* 5, k 5 m 3 (0) n n ln (0) k / 0 p 3p kn n ln p p k c* 5 n* n n ln c*, with the definitions EV mv 5 mv The constant c* is the same as that of Eq. 5. It is determined by substituting Eq. 53 into the lowest order RET equation 9 and retaining linear terms in the first Sonine polynomial to get c* Comparison of this prediction with direct Monte Carlo simulation of the kinetic equation shows a very good agreement even for relatively large dissipation 5. The tensorial properties of A, B, C ij, and D are obvious and traditionally the scalar magnitudes have been expanded in terms of Sonine polynomials. For elastic collisions truncation of the series at the leading first Sonine polynomial approximation, taking
6 5900 V. GARZÓ AND J.W. DUFTY PRE 59 into account the requirement 36, gives very good results accurate to within a few percent. In the present case this approximation is AV BV c B C ij c C D ij DV fmvcasv c D EV, f M Vn v 3/ 0 e (V/v 0 ). 56 The factor f M (V) occurs since these polynomials are defined relative to a Gaussian scalar product. The coefficients are the projections of A, B, C ij, and D along S(V), D ij (V), and E(V) c A c B m m3 5n dv AV SV k 5nk B BV SV c C 0n dvc ij VD ij V n k, c D 5n dvdvev, k, mn3 5n where use has been made of the definitions 4 and 43. The distribution function f () in this approximation is obtained from Eq. 3 f () f M m3 5n k T k n SV n k i U j j U i 3 ij UD ij V c D EV U. 60 To evaluate the kinetic parts of the transport coefficients the forms can be used, with the frequencies (), (), and () determined from Eq. 5 using the approximation 56 dvd ij Lf M D ij dvs Lf M S,. dvf M D ij D ij dv f M S S 6 The calculation of all properties now has been reduced to quadratures. An outline of the integrations involved are given in Appendix E while the contributions to the cooling rate are obtained in Appendix F. The results are summarized in Table I. The transport coefficients have been reduced by 0 and 0. The cooling rate and the frequencies and TABLE I. Summary of the main results of the paper. ff M 4c*EV m3 5n k T k n SV n k i U j j U i 3 ij UD ij Vc D EV U, * k * 5n*() 5*, 3 k *(* (0) *) 5()(3)n*, * 45n* 3 ()( 3c*), * k * 5n*() 5n* 64 ()( 3c*), 7 k * 3(* (0) *) (p*)c* 0n*() 5 () 3()c*, * k * 5n*(), k * *3 (0) * n n ln (0) * k * p* 3 n n lnp*c* c D 0 5n* n n ln c*, (0) * * 5c* c* 45n*p* 3c*, p* n*, ( 3 n*)( 6n*) 3, * (0) * () *, (0) * 5 ( )( 3c*), 3 () * 0 p* c*c D U, c*3()( )8730 (), c* , * 4 () 64c*, * * c*, * c* have been reduced by the characteristic collision frequency 0 p k / 0. The resulting dimensionless variables are */ 0, */ 0, */ 0, *n/t 0, 6 p*p/p k, */ 0, * / 0, * * / Also shown in this table is the Carnahan-Starling approximation for as a function of n*. It is easy to check that all results presented in Table I have the correct elastic limit for.
7 PRE 59 DENSE FLUID TRANSPORT FOR INELASTIC HARD SPHERES 590 FIG.. Plot of the low density transport coefficients as a function of the restitution coefficient : shear viscosity *(,0) solid line, thermal conductivity *(,0) dashed line, and *(,0) dotted line. FIG. 3. Plot of the reduced thermal conductivity *(,n*)/*(,0) as a function of the reduced density n* as obtained from the Sonine expansion of the RET solid line and from the kinetic model equation discussed in Sec. V dashed line. From the top to the bottom,, 0.8, and 0.6. The transport coefficients for granular flow of dense, weakly inelastic hard spheres were studied by Lun et al. 5 by using an approximate moment method to solve the Enskog kinetic equation. However, the method does not include all contributions to leading order in (), and neglects completely those from *() and c*(). It also gives ()0 in contrast to the finite result found here. In general, Lun s results show important quantitative discrepancies with those derived here. For example, at 0.6 and n* 0., the discrepancies are % for the reduced shear viscosity * and 73% for the reduced thermal conductivity *. Figure shows *(,n*0), *(,n*0), and *(,n*0) as functions of in the low density limit. These are the same as those given in Ref.. Recent Monte Carlo simulations of the dilute shear viscosity confirm the accuracy of this calculation. Figures 4 show *(,n*)/*(,0), *(,n*)/*(,0), and *(,n*)/*(,0) as functions of the density for, 0.8, and 0.6. Also shown in Figs. 4 are the results from the kinetic model discussed in Sec. V 6. In general, all transport coefficients are increasing functions of the density at sufficiently large, whose slopes decrease with decreasing i.e., greater dissipation. V. KINETIC MODEL The Chapman-Enskog solution to the RET described here for states with small spatial gradients is technically difficult but straightforward. For more complex states driven far from equilibrium the RET becomes intractable. In these cases it is useful to have kinetic models with the same qualitative features as the RET but with a mathematically simpler structure. The prototype example for this approach is the Bhatnagar- Gross-Krook BGK kinetic model for the Boltzmann equation 7. Generalizations to the RET for both elastic 8 and inelastic 4,0, cases have been discussed recently and applied with excellent success to shear flow far from equilibrium. The objective is to test the most sophisticated of these models for the prediction of the transport proper- FIG.. Plot of the reduced shear viscosity *(,n*)/*(,0) as a function of the reduced density n* as obtained from the Sonine expansion of the RET. From the top to the bottom, we have considered the values of the restitution coefficient, 0.8, and 0.6. FIG. 4. Plot of the reduced transport coefficient *(,n*)/*(,0) as a function of the reduced density n* as obtained from the Sonine expansion of the RET solid line and from the kinetic model equation dashed line. From the top to the bottom, 0.8 and 0.6.
8 590 V. GARZÓ AND J.W. DUFTY PRE 59 ties calculated from the RET in the previous sections. The detailed analysis leading to the model kinetic equation can be found in Ref. and will not be repeated here. The model kinetic equation is t v fpj E fqff (0) X S Y ij D ij f (0) v Vf. 64 Equation 64 has the same form as Eq. except that the RET collision operator has been replaced by the terms on the right side of Eq. 64. The first term is the exact projection of the collision operator into the subspace spanned by, V, and V see Eq. A of Appendix A. The second term represent an approximation to the collision operator in the orthogonal subspace, QP. It differs from the corresponding kinetic model for elastic collisions 8 only by the last term of the brackets which accounts for the cooling. The functions S(V) and D ij (V) are given by Eq. 44 while X and Y ij are defined by X dvs J E f (0) dvd ij J E f (0), dv f (0) S S Y ij. 65 dv f (0) D kl D kl The Chapman-Enskog solution to the kinetic model equation proceeds in the same way as described in Sec. III. To lowest order in the gradients f (0) is again the solution, and f () is determined from Eq. with only the replacement Lf () f () (0) v Vf (). 66 All results of Table I are the same except now and ** (0) * * ** 3 (0) *, with */ 0. Since the kinetic model has only one free parameter * only one of these frequencies can be reproduced exactly by the model. For example, the choice * equal to * (0) * will yield the exact shear and bulk viscosities, but the thermal conductivity and will be approximate the Prandtl number at low density is instead of /3). This limitation of the one parameter kinetic model is well known and occurs for the case of elastic collisions as well. It is remarkable that the kinetic model reproduces exactly all other features of transport, both the dependence on dissipation and on the density. This dependence is shown in Figs. 4 for different values of the restitution coefficient. Excellent agreement with the Enskog results is found over a wide range of values for n* and. VI. DISCUSSION The RET kinetic theory for hard spheres with elastic collisions provides the basis for an accurate description of transport over a wide range of densities, as confirmed by both Monte Carlo and molecular dynamics simulations. The objective of the present work has been to extend the application of this equation, with appropriate modifications, to the case of inelastic collisions. The transport processes considered are those for a fluid with small spatial gradients of the hydrodynamic fields. The macroscopic balance equations for mass, momentum, and energy or, equivalently, for number density, flow velocity, and temperature, as given in Eqs. 4 6 follow directly from moments of the RET with respect to, v, and v. These equations include the unknown pressure tensor, heat flux, and cooling rate which are given as functionals of the distribution function f. The closed set of hydrodynamic equations are obtained when a solution to the RET is known as a functional of the hydrodynamic fields. Here, this solution is obtained perturbatively for small spatial gradients using the Chapman-Enskog method well known from the case of elastic collisions. The distribution function has been calculated through first order in the gradients. Use of this in the functionals for the pressure tensor, heat flux, and cooling rate provides a representation of these as linear combinations of the gradients. The coefficients in these expressions are the transport coefficients which are functions of the density and restitution coefficient. Use of these first order results in the macroscopic balance equations provides the hydrodynamic equations which are the basis for a fluid dynamics description of granular flow. This work extends a recent corresponding analysis for a low density granular gas. The following comments are offered for context and clarification. The distribution function f is given in terms of a reference function f (0) which is the solution to Eq. 9. While the exact solution to this equation has not been obtained to date, it is easy to verify that it differs from the Maxwellian for elastic collisions. This reference function depends on space and time only through the hydrodynamic fields. In particular, for both elastic and inelastic cases, the reference state is time dependent. However, for the inelastic case the time scale for the temperature is set by the cooling rate (0) rather than the spatial gradients. This means that within the hydrodynamic description there can be two well-separated time scales, for weak gradients and strong dissipation. This should not be confused with the separation of both hydrodynamic time scales from the shorter kinetic excitations in the solution to RET. Investigations to date suggest that the kinetic excitations remain isolated from the hydrodynamic excitations so that the later dominate on sufficiently large time scales, just as in the elastic case. The deviations of f from f (0) are proportional to spatial gradients of the temperature, density, and flow velocity. The coefficients of each gradient are functions of the velocity which are determined from linear inhomogeneous integral equations. In the elastic case these coefficients are nonzero only for the temperature gradient and the traceless part of the velocity gradient. For inelastic collisions there are two additional terms, one proportional to the density gradient and one proportional to the divergence of the velocity field. The
9 PRE 59 DENSE FLUID TRANSPORT FOR INELASTIC HARD SPHERES 5903 former gives rise to an additional term in the heat flux, while the latter provides an additional term in the cooling rate. The latter provides a correction to the results of Ref. 7. The additional term of the heat flux due to the density gradient has been recognized in earlier work. The form of f (0) and the solutions to the linear integral equations to determine f and the transport coefficients were determined approximately using a first Sonine polynomial approximation. This is known to be accurate in the elastic limit and to give an excellent approximation for f (0). Preliminary comparisons of the resulting viscosity with Monte Carlo simulations at low density also confirm the accuracy of this method. The method does not make any a priori limitation on the degree of dissipation so it is expected that the dependence of the shear viscosity, bulk viscosity, thermal conductivity,, and cooling rate on both n* and is well described. A primary outcome of these calculations is a form for the hydrodynamic equations with all parameters given explicitly. The heat and momentum fluxes are determined to first order in the gradients. Since these occur as divergences in the hydrodynamic equations their contributions are of second order in the gradients. In contrast, the cooling rate is determined only to first order, showing a contribution proportional to U and given by () * in Table I. The second term in the brackets of that expression does not appear to have predicted before. In principle, the Navier-Stokes order hydrodynamics for inelastic collisions requires going one order further in the Chapman-Enskog expansion Burnett order to obtain the cooling rate to second order. This has been done in the low density limit where it is found that the corrections are very small. Consequently, it is likely that the hydrodynamics determined here provides a good basis for applications of fluid dynamics to granular flows. The Chapman-Enskog method does not make explicit use of the form of the collision operator until a late stage in the analysis. Therefore, the results can be extended in several directions without new conceptual difficulties. One generalization of importance for practical applications is to include hard objects of different shapes, a degree of roughness, and a possible velocity dependence of. In the opposite direction, it is of interest to consider simpler collision operators that allow access to solutions for more complex nonequilibrium states far from equilibrium. One such example has been explored in Sec. V. It is shown there that a simple kinetic model has the capacity to reproduce most of the n* and dependence of the transport coefficients, giving credibility to its application to more complex states. One example of the latter is a recent study of the rheology of a granular medium under shear flow 9. Good agreement with Monte Carlo simulations of the RET are obtained for a wide range of n*,, and the shear rate. As noted in the Introduction there are concerns that the qualitative changes introduced by inelastic collisions may invalidate the necessary conditions for a hydrodynamic or kinetic theory description. If the RET is a good description, the validity of hydrodynamics can be tested by Monte Carlo simulation of the evolution of initial states with small spatial gradients to check the n* and dependence of the transport properties predicted here. Further, a test of the validity of the RET is possible using molecular dynamics simulations. The calculations here, appropriate at high densities, make comparisons of transport properties in this manner as well. It is hoped that such simulations will be performed in the near future. ACKNOWLEDGMENTS We are grateful to Dr. Andrés Santos for a critical reading of the manuscript. This research was supported by National Science Foundation under Grant No. PHY 9733, and by the Dirección General de Investigación Científica y Técnica Spain through Grant Nos. PR and PB97-50 V.G.. APPENDIX A: CHAPMAN-ENSKOG EXPANSION To initiate the Chapman-Enskog expansion it is useful to make a formal transformation of the RET to expose the effects of cooling. First, define the operator J E f by J E f JE f v Vf, and decompose it according to J E f PJE fqje f, A A where P is the projection operator onto an orthonormal set constructed from,v,v: Pgvn vf (0) v dv vgv, A3 m V,c/ V 3,m V /. A4 The normalization coefficient cn dvf v m (0) V c*, A5 8 with c* given by Eq. 5. Also Q(P) is the orthogonal projection. The contributions from PJ E f can be calculated directly in terms of the collisional parts of the fluxes and the RET equation becomes t v f v Vf QJ Ef n f(0)v :Pc c 3 m 3 V q c P c : U. To zeroth order in the gradients this equation becomes Q J (0) f (0),f (0) (0) v Vf (0) 0, which is the same as Eq. 9 of Sec. III. A6 A7
10 5904 V. GARZÓ AND J.W. DUFTY PRE 59 To first order in the gradients A6 becomes f (0) r,v,t t (0) f () P (0) v Vf () QJ () f f (0) r,v,t ln n vi ln f (0) U i Q () v Vf (0) t () v f (0) 3V vln f (0) ln T, A6 () v Vf (0) n f(0)v p c 3pc c m 3 V U. A8 The first term on the right side vanishes, i.e., P v Vf () 0, A9 which follows from the condition Pf () 0. The fourth term is calculated using the form 8 t () v f (0) () v Vf (0) f (0) t () v ln n vi f (0) t () v U i where is obtained from the functional (r,r,n) by evaluating all density fields at n(r,t). The collision operator to first order then becomes QJ E f QJ (0) f (0),f (0) Lf () n n ln K f (0) ln n K j vi f (0) j U i K v Vf (0) ln T, J (0) f (0),f (0) dv dˆ ˆ gˆ g A7 f (0) v f (0) v f (0) v f (0) v, A8 v Vf (0) t () v ln T () v Vf (0). LXQJ (0) f (0),XJ (0) X,f (0), A9 A0 The macroscopic balance equations to first order in the gradients are t () U nn U0, t () U U i mn i p0, t () U T p 3nk B UT (). Use of these in Eq. A0 gives finally t () v f (0) () v Vf (0) f (0) V ln nmn v f (0) p A A A3 v Vf (0) V ln T vi f (0) V U i f (0) p 3p k v Vf (0) U. A4 Finally, the collision operator must be expanded to first order in the gradients. To do so the following results are needed: r,r n n n ln ln n, A5 KXQ dv dˆ ˆ g ˆ gˆ f (0) v Xv f (0) v Xv. A0 Substitution of Eqs. A4 and A7 into Eq. A8 gives t (0) Lf () Q () v Vf (0) p p kn n ln p p k f (0) Vm v f (0) n n ln K f (0) ln n f (0) Vm v f (0) p p k f(0) V v Vf (0) V K v Vf (0) ln T V j vi f (0) 3 ijv v f (0) K j vi f (0) ij pp k 3p k Q v Vf (0) j U i. In the last line use has been made of the identity A
11 PRE 59 DENSE FLUID TRANSPORT FOR INELASTIC HARD SPHERES 5905 P v Vf (0) f (0) 45c* m 3 V. A To expose the functional dependence of () on f () Eq. 0 is evaluated directly to first order in the gradients m n dv dv f (0) V f (0) V dˆ ˆ gˆ g 3 m dv n dv f (0) V f (0) V dˆ ˆ gˆ g 3 vi ln f (0) U i m n dv dv f (0) V f () V f () V f (0) V dˆ ˆ gˆ g 3 m 4n dv dv f (0) V f (0) V g 3 m dv dv f (0) V f (0) V g U 8n m n dv dv f (0) V f () V g 3 (0) 3 n* U (,) f (), A3 (0) m 4n dv dv f (0) V f (0) V g 3, A4 (,) f () m n dv dv f (0) V f () V g 3. A5 With these results A becomes t (0) Lf () (,) f () Q v Vf (0) p p kn n ln p p k f (0) Vm v f (0) n n ln K f (0) ln n f (0) Vm v f (0) p p k f(0) V v Vf (0) V K v Vf (0) ln T V j vi f (0) K j vi f (0) j U i i U j 3 ij U 3 K i vi f (0) ppk 3p k 3Q v Vf (0) U. A6 This is the result used in the text.
12 5906 V. GARZÓ AND J.W. DUFTY PRE 59 APPENDIX B: SOLUBILITY CONDITIONS The solubility conditions Fredholm alternative 4 for Eqs can be identified by using the property 36 and the definitions 3 35 to write these equations in the equivalent form Q (0) T T L (0) AA, B Q (0) T T LBBn n ln (0) A, B k dvd 0 / (0) ij VC ij V. C The integral on the right side can be performed using the definition 6 dvd 0 ij VC ij V Q (0) T T LC ij C ij, B3 0 dvm V i V j 3 V ijv i vj f (0) Q (0) T T LD (,) DQ v Vf (0) D. B4 Use has been made of the property 0 dvd ij VK i vj f (0) p k 0 dvd ij VK i vj f (0), C3 T D A B QT C ij v VA DQ B B B5 C ij C ij DTTA as follows from Eq. A9. Since the left side of these equations lies in the orthogonal subspace it is necessary that the right sides must as well, or equivalently P A B B6 C ij D0. These are the solubility conditions. It is straightforward to verify that they are satisfied by direct integration using the definitions 4 7 and Eq. A3. APPENDIX C: KINETIC CONTRIBUTIONS Equations for the kinetic contributions to the transport coefficients follow directly from the integral equations. Consider the viscosity. Multiply Eq. 34 by D ij and integrate over the velocity to get (0) T T k 0 dvd ij VC ij, C where is defined by Eq. 5. From dimensional analysis k T / so and Eq. C becomes k (0) dvd 0p k ij VK i vj f (0). p k C4 The remaining integral is performed using the definition of K in Eq. 8 dvd ij VK i vj f (0) dv D ij V QdV dˆ ˆ g ˆ g i f (0) V vj f(0) V f (0) V vj f (0) V. C5 The projection operator Q can be replaced by the identity operator because D ij (V ) is traceless. The scattering law for V and V is given in Sec. II. However, a simpler form is obtained by changing variables to integrate over V and V instead of V and V in the first term of Eq. C5. The Jacobian of the transformation is and ˆ gˆ g. Also, V (V,V )V V ()ˆ(ˆ g). The integral then becomes
13 PRE 59 DENSE FLUID TRANSPORT FOR INELASTIC HARD SPHERES 5907 dvd ij K i vj f (0) dv dv f (0) V vj f (0) V dˆ ˆ gˆ g i D ij V D ij V 6 m3 dv dv f (0) V vj f (0) V dˆ ˆ gˆ g ˆ jˆ V ˆ g3v j 6 m3 dv dv f (0) V f (0) V dˆˆ g4ˆ V 3 ˆ g 3 m3 dv dv f (0) V f (0) V 8 3 V gg 3 pk n*3. C6 Use of this in Eq. C4 gives the final result ˆ g i Q*S i V S i V. C9 k (0) 5 n*3. C7 p k The kinetic part of the thermal conductivity is obtained in a similar way. Multiplication of Eq. 3 by S(V) and integration over the velocity leads to k 3T (0) dvsv AV 3T (0) dvsv 5 f (0) V vf (0)V p p k Vf(0) m v f (0) 6T (0) dvsv K v Vf (0) 5k Bp k m (0) p p kc* Here Q* is the adjoint projection operator with the effect Q*S i V S i V 3m V i dvf (0) VV SV V i mv 5 c*. With this result C9 becomes dvsv K v Vf (0) C0 m3 dv dv f (0) V f (0) V V v d ˆˆ gˆ g ˆ V m 5p k dvsv K v Vf (0). The last term on the right side is, more explicitly, dvsv K v Vf (0) dv dv f (0) V C8 3 4 ˆ gˆ V 8 ˆ g V m 5 c* 3 m pk n* 5 3c*. C v V f (0) V dˆ ˆ g Use of this in Eq. C8 gives the desired result
14 5908 V. GARZÓ AND J.W. DUFTY PRE 59 k 5k Bp k m (0) p p kc* 0 n* 5 3c*. C To evaluate k, multiply Eq. 33 by S(V) and integrate over the velocity to get k 3 (0) 3n dvsv n n ln 3n dvs B 3 (0) n n ln (0) T n k 5 4 p p k m n n ln p p k c* n n ln 3n dvsv K f (0). C3 (0) AB The last integral is 3 (0) n n ln (0) T n k dvsv K f (0) dv dv f (0) V f (0) V dˆ ˆ gˆ g i Q*S i V S i V Substitution of Eq. C4 into Eq. C3 gives m3 dv dv f (0) V f (0) V dˆ ˆ gˆ g ˆ V 3 4 ˆ gˆ V 8 ˆ g V m 5 c* 3 m pk n* 4 4 3c*. k T n 3 (0) n n ln (0) k 5 4 k Bp k m n* c*. pk B m n n ln p p k c* n n ln C4 C5 APPENDIX D: COLLISIONAL TRANSFER CONTRIBUTIONS The collisional transfer contributions to the pressure tensor and heat flux are determined from Eqs. 8 and 9. Consider first the pressure tensor which becomes to first order in the gradients P c ij 4 m3 dv dv dˆ ˆ gˆ g ˆiˆ jf (0) V f () V f () V f (0) V f(0) V ˆ f (0) V f(0) V ˆ f (0) V m3 dv dv dˆ ˆ gˆ g ˆ iˆ jf () V f (0) V f(0) V ˆ f (0) V 5 m3 dv dv f () V f (0) V g i g j g ij k U l m 4 dv 4 dv f (0) V vl f (0) V dˆ g ˆ g ˆ ˆ iˆ jˆ k
15 PRE 59 DENSE FLUID TRANSPORT FOR INELASTIC HARD SPHERES n*p ij k k U l 48 m4 dv dv f (0) V f (0) V vl g g i g j g k g g j ik g i jk g k ij 5 mn*p ij k k U l 48 m4 dv dv f (0) V f (0) V g 3 g l g i g j g k g g j g k il g i g k jl g j g i kl g l g j ik g l g i jk g l g k ij g jl ik il jk kl ij 5 n*p ij k mn* dv 8 n dv f (0) V f (0) V g 3 5 j U i i U j 3 ij U ij U. D From this equation one can easily identify the collisional transfer contributions to the shear viscosity and the bulk viscosity given by Eqs. 38, 39, and 45. The collisional transfer contribution to the heat flux to first order in the gradients can be obtained in a similar way. The result is q i c m3 dv dv dˆ ˆ gˆ g ˆ Gˆ i f () V f (0) V f(0) V ˆ f (0) V 5 n*q i k i T mn* 48 nt dv dv f (0) V f (0) V g g G gg 3 gg G 4 g3. From Eq. D one may identify the collisional contributions to the thermal conductivity and the coefficient given in Eqs. 40, 4, and 46. D APPENDIX E: FIRST SONINE APPROXIMATION In this appendix an outline of the integrations appearing in the evaluation of the transport coefficients is given using the first Sonine approximations for the zeroth-order distribution 53 and for the first-order distribution 56. The frequencies * and * * are given by Eq. 6 in the first approximation. It is straightforward to show that the projection operator Q in the definition of L can be replaced by the identity in both integrations, and these integrals become the same as those appearing in the Boltzmann limit except for a factor ). The details will not be repeated here and only the results are quoted in Table I. The collisional contributions are given in terms of the dimensionless integrals I and I defined by Eqs. 45 and 46, respectively. In both cases, the integrations over the relative and center of mass variables g and G are performed and nonlinear terms in c* are neglected. The integral I is given by I 3 v 0 7 dg dge g /v 0 e G /v 0 g c* 4 5 5v 0 g G v 4 0 g G v 0 4 G g 4 4 / 3 c*. E In a similar way, the integral I is given by I 3 v 0 9 dg dge g /v 0 e G /v 0 c* 4 5 5v 0 g G v 4 0 g G v 0 4 G g 4 g3 g g G gg 3 gg G 4 8 / 7 3 c*. APPENDIX F: CONTRIBUTIONS TO THE COOLING RATE Up to the first order in gradients, the cooling rate is (0) 3 n* U (,), E F where (0) and (,) are defined by Eqs. 6 and 30, respectively. The contribution (0) can be evaluated easily by using the first Sonine approximation to f (0). Its expression appears in Table I. The two first contributions to the cooling rate were calculated by Goldshtein and Shapiro 7,
16 590 V. GARZÓ AND J.W. DUFTY PRE 59 although an error was made in the expression for c*. The corrected value for the latter was obtained in Ref. 3 and is given by Eq. 55. With this change the results of Ref. 7 agree with the first two terms on the right side of Eq. F. In the first Sonine approximation, (,) is given by (,) D m c n D dv dv f (0) V f M V EV g c* / m n c D. F In order to determine c D, multiply the integral equation 35 by E(V) and integrate over V to get (0) T T c D 5n dvevdv 5n (,) dvevq v Vf (0), with the definition F3 dvevl f M VEV. F4 dvf M VEVEV Dimensional analysis requires that (0) and are proportional to T / while the integrals on the right-hand side are independent of T. Therefore c D T / and Eq. F3 gives To convert Eq. F5 in an explicit expression for c D the two integrals appearing in the right-hand side of Eq. F5 and the quantity must be evaluated. The calculation is straightforward, leading to the results dvevq v (Vf (0) ) 5 nc*, F6 dvevk i vi f (0) c* n*nn*n, * and consequently c D 0 F7 c* , F8 /45n*pp k /3p k 3c* (0) **5c*/ c* F9 c D 5n dvevdv (0) (,) dvevq v Vf (0). F5 Finally, (,) is obtained from Eqs. F and F9. Weare not aware of any previous calculation of (,) for comparison. It is noted that (,) vanishes in the limits of elastic spheres (, arbitrary density and of dilute inelastic spheres (n*0, arbitrary inelasticity. J.J. Brey, J.W. Dufty, C-S Kim, and A. Santos, Phys. Rev. E 58, J.J. Brey, M.J. Ruiz-Montero, and D. Cubero unpublished. 3 H. van Beijeren and M.H. Ernst, Physica A 68, ; 70, 5 973; J. Stat. Phys., These references provide the statistical mechanical basis for the RET. It corrects the earlier phenomenological kinetic theory of Enskog see Ref. below. While the corrections are important for mixtures and states far from equilibrium, they are not important for the calculations done here. 4 J.J. Brey, J.W. Dufty, and A. Santos, J. Stat. Phys. 87, C.K.W. Lun, S.B. Savage, D.J. Jeffrey, and N. Chepurniy, J. Fluid Mech. 40, 3 984; C.K.K. Lun, ibid. 33, J.T. Jenkins and M.W. Richman, Arch. Ration. Mech. Anal. 87, ; Phys. Fluids 8, A. Goldshtein and M. Shapiro, J. Fluid Mech. 8, N. Sela and I. Goldhirsch, J. Fluid Mech. 36, J.M. Montanero and A. Santos, Phys. Rev. E 54, ; Phys. Fluids 9, J.W. Dufty, J.J. Brey, and A. Santos, Physica A 40, 997. J.J. Brey, J.W. Dufty, and A. Santos, J. Stat. Phys. to be published. J. Ferziger and H. Kaper, Mathematical Theory of Transport
17 PRE 59 DENSE FLUID TRANSPORT FOR INELASTIC HARD SPHERES 59 Processes in Gases North-Holland, Amsterdam, T.P.C. van Noije and M.H. Ernst, Granular Matter, N. Dunford and J. Schwartz, Linear Operators Interscience, New York, J.J. Brey, M.J. Ruiz-Montero, and D. Cubero, Phys. Rev. E 54, A comparison of the low density results in Fig. with the kinetic model was presented in Ref., showing good agreement. The application of the kinetic model in Ref. has an additional approximation, not made here, of neglecting c*. However, the results are quantitatively almost the same. 7 C. Cercignani, Theory and Applications of the Boltzmann Equation Elsevier, New York, A. Santos, J.M. Montanero, J.W. Dufty, and J.J. Brey, Phys. Rev. E 57, J.M. Montanero, V. Garzó, A. Santos, and J. J. Brey, J. Fluid Mech. to be published.
Hydrodynamics for granular flow at low density
PHYSICAL REVIEW E VOLUME 58, NUMBER 4 OCTOBER 998 Hydrodynamics for granular flow at low density J. Javier Brey Física Teórica, Universidad de Sevilla, E-4080 Sevilla, Spain James W. Dufty Department of
More informationBOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES
BOLTZMANN KINETIC THEORY FOR INELASTIC MAXWELL MIXTURES Vicente Garzó Departamento de Física, Universidad de Extremadura Badajoz, SPAIN Collaborations Antonio Astillero, Universidad de Extremadura José
More informationKinetic Models for Granular Flow
Journal of Statistical Physics, Vol. 97, Nos. 12, 1999 Kinetic Models for Granular Flow J. Javier Brey, 1 James W. Dufty, 2 and Andre s Santos 3 Received September 28, 1998; final March 29, 1999 The generalization
More informationKinetic model for the hard-sphere fluid and solid
PHYSICAL REVIEW E VOLUME 57, NUMBER FEBRUARY 1998 Kinetic model for the hard-sphere fluid and solid Andrés Santos Departamento de Física, Universidad de Extremadura, E-671 Badajoz, Spain José M. Montanero
More informationComparison between the Boltzmann and BGK equations for uniform shear flow. Vicente Garz6, Andr6s Santos
ELSEVIER Physica A 213 (1995) 426-434 PHYSICA Comparison between the Boltzmann and BGK equations for uniform shear flow Vicente Garz6, Andr6s Santos Departamento de Ffsica, Universidad de Extremadura,
More informationarxiv:cond-mat/ v2 [cond-mat.stat-mech] 21 Nov 2003
Diffusion of impurities in a granular gas Vicente Garzó Departamento de Física, Universidad de Extremadura, E-0607 Badajoz, Spain arxiv:cond-mat/0309639v [cond-mat.stat-mech] Nov 003 José María Montanero
More informationAndrés Santos* Badajoz (Spain)
Granular Poiseuille flow Andrés Santos* University of Extremadura Badajoz (Spain) *In collaboration with Mohamed Tij, Université Moulay Ismaïl, Meknès (Morocco) Departament de Física, Universitat Autònoma
More informationarxiv: v3 [cond-mat.stat-mech] 2 Dec 2015
Inelastic Maxwell models for monodisperse gas-solid flows Aleksander Kubicki Departamento de Física, Universidad de Extremadura, Universidad de Extremadura, E-0607 Badajoz, Spain arxiv:4.345v3 [cond-mat.stat-mech]
More informationKinetic Theory for Binary Granular Mixtures at Low Density
10 Kinetic Theory for Binary Granular Mixtures at Low Density V. Garzó Departamento de Física, Universidad de Extremadura, E-06071 Badajoz, Spain vicenteg@unex.es Many features of granular media can be
More informationAsymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks
Physica A 275 (2000) 483 504 www.elsevier.com/locate/physa Asymptotic solution of the Boltzmann equation for the shear ow of smooth inelastic disks V. Kumaran Department of Chemical Engineering, Indian
More informationAndrés Santos Universidad de Extremadura, Badajoz (Spain)
Andrés Santos Universidad de Extremadura, Badajoz (Spain) Outline Moment equations molecules for Maxwell Some solvable states: Planar Fourier flow Planar Fourier flow with gravity Planar Couette flow Force-driven
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
More informationEPJ E. Soft Matter and Biological Physics. EPJ.org. Segregation by thermal diffusion in moderately dense granular mixtures
EPJ E Soft Matter and Biological Physics EPJ.org your physics journal Eur. Phys. J. E 9, 61 74 009) DOI: 10.1140/epje/i009-10488-4 Segregation by thermal diffusion in moderately dense granular mixtures
More informationKINETIC THEORY OF GRANULAR SOME APPLICATIONS
KINETIC THEORY OF GRANULAR BINARY MIXTURES: SOME APPLICATIONS Vicente Garzó Departamento t de Física, Universidad id dde Extremadura, Badajoz, Spain OUTLINE 1. Introduction 2. Granular binary mixtures.
More informationRheology of two- and three-dimensional granular mixtures under uniform shear flow: Enskog kinetic theory versus molecular dynamics simulations
Rheology of two- and three-dimensional granular mixtures under uniform shear flow: Enskog kinetic theory versus molecular dynamics simulations José María Montanero Departamento de Electrónica e Ingeniería
More informationCHAPTER 9. Microscopic Approach: from Boltzmann to Navier-Stokes. In the previous chapter we derived the closed Boltzmann equation:
CHAPTER 9 Microscopic Approach: from Boltzmann to Navier-Stokes In the previous chapter we derived the closed Boltzmann equation: df dt = f +{f,h} = I[f] where I[f] is the collision integral, and we have
More informationHydrodynamic Limits for the Boltzmann Equation
Hydrodynamic Limits for the Boltzmann Equation François Golse Université Paris 7 & Laboratoire J.-L. Lions golse@math.jussieu.fr Academia Sinica, Taipei, December 2004 LECTURE 2 FORMAL INCOMPRESSIBLE HYDRODYNAMIC
More informationREGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS
1 REGULARIZATION AND BOUNDARY CONDITIONS FOR THE 13 MOMENT EQUATIONS HENNING STRUCHTRUP ETH Zürich, Department of Materials, Polymer Physics, CH-8093 Zürich, Switzerland (on leave from University of Victoria,
More informationRheology of two- and three-dimensional granular mixtures under uniform shear flow: Enskog kinetic theory versus molecular dynamics simulations
Granular Matter (2006) 8: 103 115 DOI 10.1007/s10035-006-0001-7 ORIGINAL PAPER José María Montanero Vicente Garzó MeheboobAlam Stefan Luding Rheology of two- and three-dimensional granular mixtures under
More informationarxiv:comp-gas/ v1 28 Apr 1993
Lattice Boltzmann Thermohydrodynamics arxiv:comp-gas/9304006v1 28 Apr 1993 F. J. Alexander, S. Chen and J. D. Sterling Center for Nonlinear Studies and Theoretical Division Los Alamos National Laboratory
More informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More informationLecture 5: Kinetic theory of fluids
Lecture 5: Kinetic theory of fluids September 21, 2015 1 Goal 2 From atoms to probabilities Fluid dynamics descrines fluids as continnum media (fields); however under conditions of strong inhomogeneities
More informationComparison of the DSMC Method with an Exact Solution of the Boltzmann Equation
Comparison of the DSMC Method with an Exact Solution of the Boltzmann Equation J M Montanero and A Santos Departamento de Fisica, Universidad de Extremadura, Spain Abstract Results obtained from the Direct
More informationComputational test of kinetic theory of granular media
Physica A 274 (1999) 158 170 www.elsevier.com/locate/physa Computational test of kinetic theory of granular media M.D. Shattuck a, C. Bizon b, J.B. Swift a, Harry L. Swinney a; a Center for Nonlinear Dynamics
More informationHydrodynamics of Inelastic Maxwell Models
Math. Model. Nat. Phenom. Vol. 6, No. 4, 2011, pp. 37-76 DOI: 10.1051/mmnp/20116403 Hydrodynamics of Inelastic Maxwell Models V. Garzó and A. Santos Departamento de Física, Universidad de Extremadura,
More informationModelling and numerical methods for the diffusion of impurities in a gas
INERNAIONAL JOURNAL FOR NUMERICAL MEHODS IN FLUIDS Int. J. Numer. Meth. Fluids 6; : 6 [Version: /9/8 v.] Modelling and numerical methods for the diffusion of impurities in a gas E. Ferrari, L. Pareschi
More informationChapman-Enskog expansion about nonequilibrium states with application to the sheared granular fluid
PHYSICAL REVIEW E 73 030 006 Chapman-Enskog expansion about nonequilibrium states with application to the sheared granular fluid James F. Lutsko Center for Nonlinear Phenomena and Complex Systems Université
More informationDerivation of BGK models.
Derivation of BGK models. Stéphane BRULL, Vincent PAVAN, Jacques SCHNEIDER. University of Bordeaux, University of Marseille, University of Toulon September 27 th 2014 Stéphane BRULL (Bordeaux) Derivation
More informationEnergy production rates in fluid mixtures of inelastic rough hard spheres Andrés Santos, 1 Gilberto M. Kremer, 2 and Vicente Garzó 1 1 Universidad de Extremadura, Badajoz (Spain) 2 Universidade id d Federal
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationFree cooling of a binary granular fluid of inelastic rough hard spheres Andrés Santos Universidad de Extremadura, Badajoz (Spain) In collaboration with Gilberto M. Kremer and Vicente Garzó Minimal model
More informationEl fluido de esferas penetrables. Teoría y simulación Andrés Santos Universidad de Extremadura, Badajoz, España Colaboradores: Luis Acedo (Universidad
Influence of particle roughness on some properties of granular gases Andrés Santos Universidad de Extremadura, Badajoz (Spain) In collaboration with Gilberto M. Kremer and Vicente Garzó Instituto de Ciencias
More informationarxiv: v1 [nucl-th] 9 Jun 2008
Dissipative effects from transport and viscous hydrodynamics arxiv:0806.1367v1 [nucl-th] 9 Jun 2008 1. Introduction Denes Molnar 1,2 and Pasi Huovinen 1 1 Purdue University, Physics Department, 525 Northwestern
More informationBulk equations and Knudsen layers for the regularized 13 moment equations
Continuum Mech. Thermon. 7 9: 77 89 DOI.7/s6-7-5- ORIGINAL ARTICLE Henning Struchtrup Toby Thatcher Bulk equations and Knudsen layers for the regularized 3 moment equations Received: 9 December 6 / Accepted:
More informationVelocity Distributions in a Cooling Granular Gas
School of Physical Sciences, Jawaharlal Nehru University, New Delhi 1167, India E-mail: puri@mail.jnu.ac.in Syed Rashid Ahmad School of Physical Sciences, Jawaharlal Nehru University, New Delhi 1167, India
More informationThe existence of Burnett coefficients in the periodic Lorentz gas
The existence of Burnett coefficients in the periodic Lorentz gas N. I. Chernov and C. P. Dettmann September 14, 2006 Abstract The linear super-burnett coefficient gives corrections to the diffusion equation
More informationFUNDAMENTALS OF CHEMISTRY Vol. II - Irreversible Processes: Phenomenological and Statistical Approach - Carlo Cercignani
IRREVERSIBLE PROCESSES: PHENOMENOLOGICAL AND STATISTICAL APPROACH Carlo Dipartimento di Matematica, Politecnico di Milano, Milano, Italy Keywords: Kinetic theory, thermodynamics, Boltzmann equation, Macroscopic
More informationRotational-translational relaxation effects in diatomic-gas flows
Rotational-translational relaxation effects in diatomic-gas flows V.V. Riabov Department of Computer Science, Rivier College, Nashua, New Hampshire 03060 USA 1 Introduction The problem of deriving the
More informationThe Methods of Chapman-Enskog and Grad and Applications
The Methods of Chapman-Enskog and Grad and Applications Gilberto Medeiros Kremer Departamento de Física Universidade Federal do Paraná 81531-980 Curitiba, BRAZIL Contents 1 Introduction 3 Boltzmann and
More informationKinetic relaxation models for reacting gas mixtures
Kinetic relaxation models for reacting gas mixtures M. Groppi Department of Mathematics and Computer Science University of Parma - ITALY Main collaborators: Giampiero Spiga, Giuseppe Stracquadanio, Univ.
More informationLattice Boltzmann Method
3 Lattice Boltzmann Method 3.1 Introduction The lattice Boltzmann method is a discrete computational method based upon the lattice gas automata - a simplified, fictitious molecular model. It consists of
More informationGas flow around a longitudinally oscillating plate at arbitrary ratio of collision frequency to oscillation frequency
Gas flow around a longitudinally oscillating plate at arbitrary ratio of collision frequency to oscillation frequency Feli Sharipov and Denize Kalempa Departamento de Física, Universidade Federal do Paraná
More informationPlasmas as fluids. S.M.Lea. January 2007
Plasmas as fluids S.M.Lea January 2007 So far we have considered a plasma as a set of non intereacting particles, each following its own path in the electric and magnetic fields. Now we want to consider
More informationMechanisms of Cluster Formation in Force-Free Granular Gases
Math. Model. Nat. Phenom. Vol. 6, No. 4, 011, pp. 175-190 DOI: 1051/mmnp/017108 Mechanisms of Cluster Formation in Force-Free Granular Gases C. Salueña 1, L. Almazán 1, and N. V. Brilliantov 3 1 Department
More informationTransport Properties of a Kinetic Model for Chemical Reactions without Barriers
Transport Properties of a Kinetic Model for Chemical Reactions without arriers Giselle M. Alves, Gilberto M. Kremer and Ana Jacinta Soares Escola Técnica, Universidade Federal do Paraná, Curitiba, razil
More informationarxiv:cond-mat/ v3 [cond-mat.stat-mech] 27 Mar 2001
Velocity correlations and the structure of nonequilibrium hard core fluids arxiv:cond-mat/8359v3 [cond-mat.stat-mech] 27 Mar 21 James F. Lutsko Physics Division, Starlab Rue Engelandstraat 555 B-118 Brussels,
More informationHIGH FRICTION LIMIT OF THE KRAMERS EQUATION : THE MULTIPLE TIME SCALE APPROACH. Lydéric Bocquet
HIGH FRICTION LIMIT OF THE KRAMERS EQUATION : THE MULTIPLE TIME SCALE APPROACH Lydéric Bocquet arxiv:cond-mat/9605186v1 30 May 1996 Laboratoire de Physique, Ecole Normale Supérieure de Lyon (URA CNRS 1325),
More informationarxiv: v1 [cond-mat.stat-mech] 26 Jan 2008
Long Time Tail of the Velocity Autocorrelation Function in a Two-Dimensional Moderately Dense Hard Disk Fluid Masaharu Isobe Graduate School of Engineering, Nagoya Institute of Technology, Nagoya 466-8555,
More information12. MHD Approximation.
Phys780: Plasma Physics Lecture 12. MHD approximation. 1 12. MHD Approximation. ([3], p. 169-183) The kinetic equation for the distribution function f( v, r, t) provides the most complete and universal
More informationFluid equations, magnetohydrodynamics
Fluid equations, magnetohydrodynamics Multi-fluid theory Equation of state Single-fluid theory Generalised Ohm s law Magnetic tension and plasma beta Stationarity and equilibria Validity of magnetohydrodynamics
More informationTransport coefficients for dense hard-disk systems. Abstract
Transport coefficients for dense hard-disk systems Ramón García-Rojo Institute for Computational Physics, University of Stuttgart, Pfaffenwaldring 27, 70569 Stuttgart, Germany Stefan Luding Technische
More informationSimple shear flow of collisional granular-fluid mixtures
Manuscript Click here to download Manuscript: Manuscript_r1.docx 1 Simple shear flow of collisional granular-fluid mixtures 2 3 4 5 D. Berzi 1 1 Department of Environmental, Hydraulic, Infrastructure,
More informationA Unified Gas-kinetic Scheme for Continuum and Rarefied Flows
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows K. Xu and J.C. Huang Mathematics Department, Hong Kong University of Science and Technology, Hong Kong Department of Merchant Marine, National
More informationThe Source of 2nd Moment in Dilute Granular Flows of Highly Inelastic Spheres
Worcester Polytechnic Institute Digital WPI Mechanical Engineering Faculty Publications Department of Mechanical Engineering 11-1989 The Source of 2nd Moment in Dilute Granular Flows of Highly Inelastic
More informationRevisit to Grad s Closure and Development of Physically Motivated Closure for Phenomenological High-Order Moment Model
Revisit to Grad s Closure and Development of Physically Motivated Closure for Phenomenological High-Order Moment Model R. S. Myong a and S. P. Nagdewe a a Dept. of Mechanical and Aerospace Engineering
More informationKinetic theory of the ideal gas
Appendix H Kinetic theory of the ideal gas This Appendix contains sketchy notes, summarizing the main results of elementary kinetic theory. The students who are not familiar with these topics should refer
More informationDry and wet granular flows. Diego Berzi
Dry and wet granular flows Diego Berzi Outline 2 What? Why? How? When? Who? Where? Then? What? Granular flows many solid moving particles 3 particle mass is large (at least 10 20 molecular masses). Hence,
More informationPhysical Modeling of Multiphase flow. Boltzmann method
with lattice Boltzmann method Exa Corp., Burlington, MA, USA Feburary, 2011 Scope Scope Re-examine the non-ideal gas model in [Shan & Chen, Phys. Rev. E, (1993)] from the perspective of kinetic theory
More informationFrom Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray
From Boltzmann Equations to Gas Dynamics: From DiPerna-Lions to Leray C. David Levermore Department of Mathematics and Institute for Physical Science and Technology University of Maryland, College Park
More informationLattice Boltzmann Methods for Fluid Dynamics
Lattice Boltzmann Methods for Fluid Dynamics Steven Orszag Department of Mathematics Yale University In collaboration with Hudong Chen, Isaac Goldhirsch, and Rick Shock Transient flow around a car LBGK
More informationMulticomponent diffusion in gases and plasma mixtures
High Temperatures ^ High Pressures, 2002, volume 34, pages 109 ^ 116 15 ECTP Proceedings pages 1337 ^ 1344 DOI:10.1068/htwu73 Multicomponent diffusion in gases and plasma mixtures Irina A Sokolova Institute
More informationThe dynamics of small particles whose size is roughly 1 µmt or. smaller, in a fluid at room temperature, is extremely erratic, and is
1 I. BROWNIAN MOTION The dynamics of small particles whose size is roughly 1 µmt or smaller, in a fluid at room temperature, is extremely erratic, and is called Brownian motion. The velocity of such particles
More informationEstimations of Rotational Relaxation Parameters in Diatomic Gases
Estimations of Rotational Relaxation Parameters in Diatomic Gases Vladimir V. Riabov Department of Mathematics and Computer Science, Rivier College, 420 S. Main St., Nashua, NH 03060, USA Abstract: The
More informationLong-Time-Tail Effects on Lyapunov Exponents of a Random, Two-Dimensional Field-Driven Lorentz Gas
Journal of Statistical Physics, Vol. 00, Nos., 000 Long-Time-Tail Effects on Lyapunov Exponents of a Random, Two-Dimensional Field-Driven Lorentz Gas D. Panja, J. R. Dorfman, and Henk van Beijeren Received
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationDSMC Simulation of Binary Rarefied Gas Flows between Parallel Plates and Comparison to Other Methods
Simulation of Binary Rarefied Gas Flows between Parallel Plates and Comparison to Other Methods L. Szalmas Department of Mechanical Engineering, University of Thessaly, Pedion Areos, Volos 38334, Greece
More informationDivergence of the gradient expansion and the applicability of fluid dynamics Gabriel S. Denicol (IF-UFF)
Divergence of the gradient expansion and the applicability of fluid dynamics Gabriel S. Denicol (IF-UFF) arxiv:1608.07869, arxiv:1711.01657, arxiv:1709.06644 Frankfurt University 1.February.2018 Preview
More informationRate of Heating and Cooling
Rate of Heating and Cooling 35 T [ o C] Example: Heating and cooling of Water E 30 Cooling S 25 Heating exponential decay 20 0 100 200 300 400 t [sec] Newton s Law of Cooling T S > T E : System S cools
More informationPage 52. Lecture 3: Inner Product Spaces Dual Spaces, Dirac Notation, and Adjoints Date Revised: 2008/10/03 Date Given: 2008/10/03
Page 5 Lecture : Inner Product Spaces Dual Spaces, Dirac Notation, and Adjoints Date Revised: 008/10/0 Date Given: 008/10/0 Inner Product Spaces: Definitions Section. Mathematical Preliminaries: Inner
More informationRegularization of the Chapman-Enskog Expansion and Its Description of Shock Structure
NASA/CR-2001-211268 ICASE Report No. 2001-39 Regularization of the Chapman-Enskog Expansion and Its Description of Shock Structure Kun Xu Hong Kong University, Kowloon, Hong Kong ICASE NASA Langley Research
More informationGlobal Weak Solution to the Boltzmann-Enskog equation
Global Weak Solution to the Boltzmann-Enskog equation Seung-Yeal Ha 1 and Se Eun Noh 2 1) Department of Mathematical Science, Seoul National University, Seoul 151-742, KOREA 2) Department of Mathematical
More informationExtended kinetic theory applied to dense, granular, simple shear flows
Noname manuscript No. (will be inserted by the editor) Diego Berzi Extended kinetic theory applied to dense, granular, simple shear flows Received: date / Accepted: date Abstract We apply the Extended
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationTheorem on the Distribution of Short Time Single Particle Displacements
Theorem on the Distribution of Short Time Single Particle Displacements R. van Zon and E. G. D. Cohen The Rockefeller University, 1230 York Avenue, New York, NY 10021, USA September 30, 2005 Abstract The
More informationPhysics 221A Fall 2018 Notes 22 Bound-State Perturbation Theory
Copyright c 2018 by Robert G. Littlejohn Physics 221A Fall 2018 Notes 22 Bound-State Perturbation Theory 1. Introduction Bound state perturbation theory applies to the bound states of perturbed systems,
More informationLow Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model
Low Variance Particle Simulations of the Boltzmann Transport Equation for the Variable Hard Sphere Collision Model G. A. Radtke, N. G. Hadjiconstantinou and W. Wagner Massachusetts Institute of Technology,
More informationLattice Bhatnagar Gross Krook model for the Lorenz attractor
Physica D 154 (2001) 43 50 Lattice Bhatnagar Gross Krook model for the Lorenz attractor Guangwu Yan a,b,,liyuan a a LSEC, Institute of Computational Mathematics, Academy of Mathematics and System Sciences,
More informationForcing independent velocity distributions in an experimental granular fluid
Forcing independent velocity distributions in an experimental granular fluid P. M. Reis,* R. A. Ingale, and M. D. Shattuck Benjamin Levich Institute, The City College of the City University of New York,
More informationAnalysis of the Reaction Rate of Binary Gaseous Mixtures near from the Chemical Equilibrium
Analysis of the Reaction Rate of Binary Gaseous Mixtures near from the Chemical Equilibrium Adriano W. da Silva May 30, 2013 Abstract A binary gaseous mixture with reversible reaction of type A + A = B
More informationNon-equilibrium mixtures of gases: modeling and computation
Non-equilibrium mixtures of gases: modeling and computation Lecture 1: and kinetic theory of gases Srboljub Simić University of Novi Sad, Serbia Aim and outline of the course Aim of the course To present
More informationCitation PHYSICS OF FLUIDS (2003), 15(3): 68.
Vapor flows condensing at incidence Titlephase in the presence of a nonconde condensation Author(s) Taguchi, S; Aoki, K; Takata, S Citation PHYSICS OF FLUIDS (2003), 15(3): 68 Issue Date 2003-03 URL http://hdl.handle.net/2433/39782
More informationMonte Carlo methods for kinetic equations
Monte Carlo methods for kinetic equations Lecture 4: Hybrid methods and variance reduction Lorenzo Pareschi Department of Mathematics & CMCS University of Ferrara Italy http://utenti.unife.it/lorenzo.pareschi/
More informationSteady Flow and its Instability of Gravitational Granular Flow
Steady Flow and its Instability of Gravitational Granular Flow Namiko Mitarai Department of Chemistry and Physics of Condensed Matter, Graduate School of Science, Kyushu University, Japan. A thesis submitted
More informationGeneralized Gas Dynamic Equations
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 2009 Orlando Florida AIAA 2009-672 Generalized Gas Dynamic Equations Kun Xu Mathematics Department
More informationPhysical models for plasmas II
Physical models for plasmas II Dr. L. Conde Dr. José M. Donoso Departamento de Física Aplicada. E.T.S. Ingenieros Aeronáuticos Universidad Politécnica de Madrid Physical models,... Plasma Kinetic Theory
More informationLattice Boltzmann Method
Chapter 5 Lattice Boltzmann Method 5.1 Pseudo-kinetic computational methods Understanding the flow properties of complex fluids in complex topologies is of importance to technological applications and
More informationDynamics of sheared gases
Computer Physics Communications 121 122 (1999) 225 230 www.elsevier.nl/locate/cpc Dynamics of sheared gases Patricio Cordero a,1, Dino Risso b a Departamento de Física, Facultad de Ciencias Físicas y Matemáticas,
More informationDiscontinuous shear thickening on dense non-brownian suspensions via lattice Boltzmann method
Discontinuous shear thickening on dense non-brownian suspensions via lattice Boltzmann method Pradipto and Hisao Hayakawa Yukawa Insitute for Theoretical Physics Kyoto University Rheology of disordered
More informationPerfect entropy functions of the Lattice Boltzmann method
EUROPHYSICS LETTERS 15 July 1999 Europhys. Lett., 47 (2), pp. 182-188 (1999) Perfect entropy functions of the Lattice Boltzmann method I. V. Karlin( ), A. Ferrante and H. C. Öttinger ETH Zürich, Department
More information1.3 Molecular Level Presentation
1.3.1 Introduction A molecule is the smallest chemical unit of a substance that is capable of stable, independent existence. Not all substances are composed of molecules. Some substances are composed of
More informationSolution of time-dependent Boltzmann equation for electrons in non-thermal plasma
Solution of time-dependent Boltzmann equation for electrons in non-thermal plasma Z. Bonaventura, D. Trunec Department of Physical Electronics Faculty of Science Masaryk University Kotlářská 2, Brno, CZ-61137,
More informationDynamics and Patterns in Sheared Granular Fluid : Order Parameter Description and Bifurcation Scenario
Dynamics and Patterns in Sheared Granular Fluid : Order Parameter Description and Bifurcation Scenario NDAMS Workshop @ YITP 1 st November 2011 Meheboob Alam and Priyanka Shukla Engineering Mechanics Unit
More informationEquation of state of additive hard-disk fluid mixtures: A critical analysis of two recent proposals
PHYSICAL REVIEW E 66, 0310 00 Equation of state of additive hard-disk fluid mixtures: A critical analysis of two recent proposals M. López de Haro* Centro de Investigación en Energía, UNAM, Temixco, Morelos
More informationKinetic Models and Gas-Kinetic Schemes with Rotational Degrees of Freedom for Hybrid Continuum/Kinetic Boltzmann Methods
Kinetic Models and Gas-Kinetic Schemes with Rotational Degrees of Freedom for Hybrid Continuum/Kinetic Boltzmann Methods Simone Colonia, René Steijl and George N. Barakos CFD Laboratory - School of Engineering
More informationOn Entropy Flux Heat Flux Relation in Thermodynamics with Lagrange Multipliers
Continuum Mech. Thermodyn. (1996) 8: 247 256 On Entropy Flux Heat Flux Relation in Thermodynamics with Lagrange Multipliers I-Shih Liu Instituto de Matemática Universidade do rio de Janeiro, Caixa Postal
More informationKinetic Theory of a Confined Quasi-Two-Dimensional Gas of Hard Spheres. J. Javier Brey *, Vicente Buzón, Maria Isabel García de Soria and Pablo Maynar
entropy Article Kinetic Theory of a Confined Quasi-Two-Dimensional Gas of Hard Spheres J. Javier Brey *, Vicente Buzón, Maria Isabel García de Soria and Pablo Maynar Física Teórica, Universidad de Sevilla,
More informationInvestigation of Kinetic Theory for Rapid Flows of Polydisperse and/or Continuous Particle Size Distributions
University of Colorado, Boulder CU Scholar Chemical & Biological Engineering Graduate Theses & Dissertations Chemical & Biological Engineering Spring 1-1-2012 Investigation of Kinetic Theory for Rapid
More informationStudies on transport phenomena of cohesive granular particles. Satoshi Takada
Studies on transport phenomena of cohesive granular particles Satoshi Takada January 7, 216 c Copyright by Satoshi Takada 216 Abstract In this thesis, we numerically and theoretically study the transport
More informationA Solution of the Two-dimensional Boltzmann Transport Equation
Applied Mathematical Sciences, Vol. 9, 2015, no. 147, 7347-7355 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2015.510665 A Solution of the Two-dimensional Boltzmann Transport Equation F.
More information