[論文レビュー] 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) で表される。
• 熱伝導度 $ \mathcal{L} $ を、さまざまな摩擦強度 $ L_0 $、相互作用数 $ n $、および粒子密度の下で計算するためにシミュレーションが実施された。
• 理論的予測とシミュレーションデータを比較し、特に $ 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 $ が非常に大きい領域では、相互作用数 $ n $ に依存してシミュレーションデータセット間で $ \mathcal{L} $ に約2倍の乖離が観察されたが、その明確な説明は得られなかった。
• 観察された乖離は、等温DPDの粘性度シミュレーションでも同様に見られるものと類似しており、低 $ n $ 領域における非マルコフ的効果や有限サイズ補正の影響が原因である可能性がある。