Megías Alberto, Santos Andrés
Departamento de Física, Universidad de Extremadura, E-06006 Badajoz, Spain.
Instituto de Computación Científica Avanzada (ICCAEx), Universidad de Extremadura, E-06006 Badajoz, Spain.
Phys Rev E. 2021 Sep;104(3-1):034902. doi: 10.1103/PhysRevE.104.034902.
Conditions for the stability under linear perturbations around the homogeneous cooling state are studied for dilute granular gases of inelastic and rough hard disks or spheres with constant coefficients of normal (α) and tangential (β) restitution. After a formally exact linear stability analysis of the Navier-Stokes-Fourier hydrodynamic equations in terms of the translational (d_{t}) and rotational (d_{r}) degrees of freedom, the transport coefficients derived in the companion paper [A. Megías and A. Santos, "Hydrodynamics of granular gases of inelastic and rough hard disks or spheres. I. Transport coefficients" Phys. Rev. E 104, 034901 (2021)10.1103/PhysRevE.104.034901] are employed. Known results for hard spheres [Garzó, Santos, and Kremer, Phys. Rev. E 97, 052901 (2018)10.1103/PhysRevE.97.052901] are recovered by setting d_{t}=d_{r}=3, while novel results for hard disks (d_{t}=2, d_{r}=1) are obtained. In the latter case, a high-inelasticity peculiar region in the (α,β) parameter space is found, inside which the critical wave number associated with the longitudinal modes diverges. Comparison with event-driven molecular dynamics simulations for dilute systems of hard disks at α=0.2 shows that this theoretical region of absolute instability may be an artifact of the extrapolation to high inelasticity of the approximations made in the derivation of the transport coefficients, although it signals a shrinking of the conditions for stability. In the case of moderate inelasticity (α=0.7), however, a good agreement between the theoretical predictions and the simulation results is found.
对于具有恒定法向(α)和切向(β)恢复系数的非弹性粗糙硬磁盘或球体的稀颗粒气体,研究了围绕均匀冷却状态的线性扰动下的稳定性条件。在对Navier-Stokes-Fourier流体动力学方程进行了关于平移(dₜ)和旋转(dᵣ)自由度的形式上精确的线性稳定性分析之后,采用了在配套论文[A. Megías和A. Santos,“非弹性粗糙硬磁盘或球体的颗粒气体的流体动力学。I. 输运系数”,《物理评论E》104,034901(2021)10.1103/PhysRevE.104.034901]中导出的输运系数。通过设置dₜ = dᵣ = 3,可以得到硬球的已知结果[Garzó,Santos和Kremer,《物理评论E》97,052901(2018)10.1103/PhysRevE.97.052901],同时得到硬磁盘(dₜ = 2,dᵣ = 1)的新结果。在后一种情况下,在(α,β)参数空间中发现了一个高非弹性特殊区域,在该区域内与纵向模式相关的临界波数发散。与α = 0.2时硬磁盘稀系统的事件驱动分子动力学模拟的比较表明,这个绝对不稳定的理论区域可能是在输运系数推导中所做近似向高非弹性外推的产物,尽管它表明稳定性条件在缩小。然而,在中等非弹性(α = 0.7)的情况下,理论预测与模拟结果之间发现了很好的一致性。