[论文解读] Dissipative Particle Dynamics with Energy Conservation: Dynamic and Transport Properties

Josep Bonet Àvalos, Allan D. Mackie|arXiv (Cornell University)|Feb 16, 2000

Material Dynamics and Properties被引用 25

一句话总结

本文提出了一种能量守恒的耗散粒子动力学模型（DPDE），可实现对动量与热量输运的精确模拟。通过结合局域平衡近似与直接模拟，推导出输运系数，并发现当摩擦系数较小时，热导率由耗散主导转变为扩散主导；模拟结果验证了理论标度形式，尤其在粒子相互作用数较大时表现显著。

ABSTRACT

Simulation results of the thermal conductivity ${\cal L}$ of Dissipative Particle Dynamics model with Energy Conservation (DPDE) are reported. We also present an analysis of the transport equations and the transport coefficients for DPDE based on a local equilibrium approximation. This approach is valid when the particle-particle thermal conductivity $λ$ and the friction coefficient $ζ$ are large. A qualitative derivation of the scaling form of the kinetic contribution of the transport of energy is derived, yielding two different forms for the kinetic contribution to the heat transport, depending on the value of $λ$. We find agreement between the theoretically predicted value for ${\cal L}$ and the simulation results, for large $λ$ and many particles interacting at one time. Significant differences are found for small number of interacting particles, even with large $λ$. For smaller values of $λ$, the obtained macroscopic thermal conductivity is dominated by diffusive transport, in agreement with the proposed scaling form.

研究动机与目标

• 开发一种在每个时间步均保持能量守恒的 DPDE 模型，以实现对复杂流体中热量输运的真实模拟。
• 利用局域平衡近似推导 DPDE 的输运方程与系数，尤其关注热导率。
• 研究在不同相互作用条件下，动能与耗散项对能量输运的标度行为。
• 将理论预测与直接模拟结果进行比较，特别是在局域平衡假设失效的区域。
• 识别理论与模拟之间差异的根源，尤其在相互作用数较少的区域。

提出的方法

• DPDE 模型通过引入粒子内能，并在粒子相互作用过程中通过耗散力与随机力中的功交换实现能量守恒，从而扩展了标准 DPD 模型。
• 应用局域平衡近似，假设动量与能量分布能快速弛豫至具有特征时间 $ t_p $ 和 $ t_q $ 的局域麦克斯韦-玻尔兹曼分布。
• 输运系数以耗散系数 $ L_0 $ 和 $ ζ_0 $ 的负幂次形式展开，其中动能贡献由公式 (3.36) 和 (3.39) 表达。
• 通过模拟计算不同摩擦强度 $ L_0 $、相互作用数 $ n $ 与粒子密度下的热导率 $ \mathcal{L} $。
• 将理论预测与模拟数据进行比较，尤其关注大 $ L_0 $ 与小 $ L_0 $ 区域，以评估局域平衡假设的有效性。
• 分析从耗散输运到扩散输运的交叉行为，预测在小 $ L_0 $ 时扩散主导，且该结果与 $ L_0 $ 无关。

实验结果

研究问题

• RQ1 在局域平衡假设下，DPDE 模型的热导率 $ \mathcal{L} $ 如何随耗散系数 $ L_0 $ 变化？
• RQ2 DPDE 中能量输运的动能贡献具有何种性质，其如何依赖于 $ \lambda $ 与 $ \zeta $ ？
• RQ3 在低 $ L_0 $ 与小相互作用数条件下，$ \mathcal{L} $ 的模拟结果与理论预测相比如何？
• RQ4 为何即使在 $ \lambda $ 较大时，低相互作用数区域仍会出现理论与模拟之间的差异？
• RQ5 当 $ L_0 $ 较小时，哪种物理机制主导热量输运，其对 $ \mathcal{L} $ 的标度行为有何影响？

主要发现

• 当相互作用粒子数较大且 $ \lambda $ 较大时，基于局域平衡的理论预测与模拟结果高度吻合。
• 当相互作用粒子数较少时，即使 $ \lambda $ 较大，理论与模拟之间仍出现显著偏差，表明局域平衡假设失效。
• 在低 $ L_0 $ 条件下，热导率 $ \mathcal{L} $ 与 $ L_0 $ 无关，表明发生向扩散输运主导的转变，与标度公式 (3.39) 的预测一致。
• 模拟结果证实 $ \mathcal{L} $ 随 $ L_0 $ 减小而增加，与低摩擦区域扩散输运占主导的预测一致。
• 在大 $ L_0 $ 条件下，不同模拟数据集之间 $ \mathcal{L} $ 存在约两倍的差异，且该差异与相互作用数 $ n $ 有关，但尚无明确解释。
• 该观测到的偏差可能源于与等温 DPD 粘度模拟中观察到的相同机制，可能与低 $ n $ 区域的非马尔可夫效应或有限尺寸修正有关。