[논문 리뷰] Dissipative Particle Dynamics with Energy Conservation: Dynamic and Transport Properties

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

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 $ 및 $ \zeta_0 $ 의 역수 거듭제곱의 급수 전개로 유도되며, 운동적 기여는 식 (3.36) 및 (3.39)로 표현된다.
• 열전도도 $ \mathcal{L} $ 를 다양한 마찰 강도 $ L_0 $, 상호작용 수 $ n $ 및 입자 밀도에서 계산하기 위해 시뮬레이션을 수행한다.
• 이론적 예측을 대조하여, 특히 큰 $ L_0 $ 및 작은 $ L_0 $ 영역에서 국부 평형 가정의 타당성을 평가한다.
• 소산성에서 확산성 운반으로의 전이 행동을 분석하며, 작은 $ L_0 $ 에서는 확산성 지배가 예측되며, 이는 $ L_0 $ 와 무관하다.

실험 결과

연구 질문

• RQ1 국부 평형 가정 하에서, DPDE 모델의 열전도도 $ \mathcal{L} $ 는 소산 계수 $ L_0 $ 와 어떻게 스케일링되는가?
• RQ2 DPDE에서 에너지 운반에 대한 운동적 기여의 성격은 무엇이며, $ \lambda $ 와 $ \zeta $ 에 어떻게 의존하는가?
• RQ3 소수의 상호작용 수가 있는 영역에서, $ \mathcal{L} $ 의 시뮬레이션 결과는 이론적 예측과 어떻게 비교되는가?
• RQ4 왜 고 $ \lambda $ 에서조차도 상호작용 수가 적은 영역에서 이론과 시뮬레이션 간의 격리가 발생하는가?
• RQ5 소마찰 조건에서 열 운반을 지배하는 물리적 메커니즘은 무엇이며, 이는 $ \mathcal{L} $ 의 스케일링에 어떻게 영향을 미치는가?

주요 결과

• 대규모 상호작용 수와 큰 $ \lambda $ 조건에서, 국부 평형 근사를 기반으로 한 열전도도 $ \mathcal{L} $ 의 이론적 예측은 시뮬레이션 결과와 잘 일치한다.
• 작은 상호작용 수 조건에서는 큰 $ \lambda $ 를 가진다 해도 이론과 시뮬레이션 간에 뚜렷한 격리가 발생하며, 이는 국부 평형 가정의 붕괴를 시사한다.
• 저마찰 조건($ L_0 $ 가 작은 경우), 열전도도 $ \mathcal{L} $ 는 $ L_0 $ 에 의존하지 않게 되며, 이는 확산성 지배로의 전이가 일어남을 나타내며, 식 (3.39)의 스케일링 형태에 의해 예측된 그대로이다.
• 시뮬레이션 결과는 $ L_0 $ 가 감소할수록 $ \mathcal{L} $ 가 증가하는 경향을 보이며, 이는 저마찰 영역에서 확산성 운반의 지배가 이루어짐을 일관되게 뒷받침한다.
• 큰 $ L_0 $ 에서는 상호작용 수 $ n $ 에 따라 시뮬레이션 데이터 세트 간에 약 두 배의 격리가 관측되며, 이에 대한 명확한 설명은 아직 찾지 못했다.
• 관측된 격리는 등온 DPD 점성도 시뮬레이션에서 관찰된 것과 같은 기원에서 비롯되었을 가능성이 있으며, 저 $ n $ 영역에서의 비마르코프 효과 또는 유한 체적 보정 요인 때문일 수 있다.