Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Segregation of an intruder in a heated granular dense gas

Vicente Garzó and Francisco Vega Reyes
Phys. Rev. E 85, 021308 – Published 22 February 2012

Abstract

A recent segregation criterion [Phys. Rev. E 78, 020301(R) (2008)] based on the thermal diffusion factor Λ of an intruder in a heated granular gas described by the inelastic Enskog equation is revisited. The sign of Λ provides a criterion for the transition between the Brazil-nut effect (BNE) and the reverse Brazil-nut effect (RBNE). The present theory incorporates two extra ingredients not accounted for by the previous theoretical attempt. First, the theory is based upon the second Sonine approximation to the transport coefficients of the mass flux of the intruder. Second, the dependence of the temperature ratio (intruder temperature over that of the host granular gas) on the solid volume fraction is taken into account in the first and second Sonine approximations. In order to check the accuracy of the Sonine approximation considered, the Enskog equation is also numerically solved by means of the direct simulation Monte Carlo method to get the kinetic diffusion coefficient D0. The comparison between theory and simulation shows that the second Sonine approximation to D0 yields an improvement over the first Sonine approximation when the intruder is lighter than the gas particles in the range of large inelasticity. With respect to the form of the phase diagrams for the BNE-RBNE transition, the kinetic theory results for the factor Λ indicate that while the form of these diagrams depends sensitively on the order of the Sonine approximation considered when gravity is absent, no significant differences between both Sonine solutions appear in the opposite limit (gravity dominates the thermal gradient). In the former case (no gravity), the first Sonine approximation overestimates both the RBNE region and the influence of dissipation on thermal diffusion segregation.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
2 More
  • Received 23 March 2011

DOI:https://doi.org/10.1103/PhysRevE.85.021308

©2012 American Physical Society

Authors & Affiliations

Vicente Garzó* and Francisco Vega Reyes

  • Departamento de Física, Universidad de Extremadura, E-06071 Badajoz, Spain

  • *vicenteg@unex.es; http://www.unex.es/eweb/fisteor/vicente/
  • fvega@unex.es

Article Text (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 85, Iss. 2 — February 2012

Reuse & Permissions
Access Options
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×

Images

  • Figure 1
    Figure 1
    Plot of Δφγφ versus the solid volume fraction φ for a hard-sphere gas in the case m0/m=8 and σ0/σ=2. Three different cases are considered: (a) α=α0=0.9, (b) α=α0=0.8, and (c) α=α0=0.5.Reuse & Permissions
  • Figure 2
    Figure 2
    Temperature ratio γT0/T versus the (common) coefficient of restitution α0=α for a hard-sphere system (d=3) in the case σ0/σ=2, φ=0.2 and two values of the mass ratio m0/m: m0/m=2 (solid line and circles) and m0/m=18 (dashed line and squares). The lines are the theoretical results, and the symbols refer to the numerical results obtained from the DSMC method.Reuse & Permissions
  • Figure 3
    Figure 3
    Reduced kinetic diffusion coefficient D0(α)/D0(1) as a function of the (common) coefficient of restitution α=α0 for a system of hard spheres with ω=2 and φ=0.2. The left panel is for M=1/8 while the right panel is for M=2. The solid lines correspond to the second Sonine approximation and the dashed lines refer to the first Sonine approximation. The symbols are the results obtained from Monte Carlo simulations. Here, D0(1) is the elastic value of the kinetic diffusion coefficient consistently obtained in each approximation.Reuse & Permissions
  • Figure 4
    Figure 4
    Reduced kinetic diffusion coefficient D0(α)/D0(1) as a function of the (common) coefficient of restitution α=α0 for a system of hard spheres with ω=2 and φ=0.2. The left panel is for M=1/5 while the right panel is for M=1/2. The solid lines correspond to the second Sonine approximation, and the dashed lines refer to the first Sonine approximation. The symbols are the results obtained from Monte Carlo simulations. Here, D0(1) is the elastic value of the kinetic diffusion coefficient consistently obtained in each approximation.Reuse & Permissions
  • Figure 5
    Figure 5
    Plot of the thermal diffusion factor Λ obtained from the second Sonine approximation as a function of the diameter ratio σ0/σ for a dilute (φ=0) hard-sphere gas in the absence of gravity (g*=0) when the intruder and gas particles have the same mass density [m0/m=(σ0/σ)3]. Two different cases are considered: α=α0=0.7 (solid line) and α=0.5, α0=0.9 (dashed line).Reuse & Permissions
  • Figure 6
    Figure 6
    Phase diagram for BNE-RBNE at φ=0.1 in the absence of gravity with α=0.8. Points above the curve correspond to Λ>0 (BNE), while points below the curve correspond to Λ<0 (RBNE). The dotted line is the result obtained in Ref. [5], the dashed line is the result obtained here from the first Sonine approximation, while the solid line is the result derived from the second Sonine approximation.Reuse & Permissions
  • Figure 7
    Figure 7
    The same as in Fig. 6 but at φ=0.5.Reuse & Permissions
  • Figure 8
    Figure 8
    Phase diagram for BNE-RBNE for φ=0.2 in the absence of a thermal gradient (|g*|) with α=0.9. The dotted and dashed lines refer to the results obtained in Ref. [5] and here, respectively, from the first Sonine approximation while the solid line is the result derived from the second Sonine approximation.Reuse & Permissions
  • Figure 9
    Figure 9
    Phase diagram for BNE-RBNE at φ=0.1 with g*=0 and two different values of α. The solid lines are the results from the second Sonine approximation, while the dashed line is the result derived from the latter approach for α=0.5 but assuming T0=T.Reuse & Permissions
×

Sign up to receive regular email alerts from Physical Review E

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×