Heat transfer simulation using energy conservative dissipative particle dynamics

Dissipative particle dynamics with energy conservation (eDPD) was used to investigate conduction heat transfer in two dimensions under steady-state condition. Various types of boundary condition were implemented to the conduction domain. Besides, 2D conduction with internal heat generation was studied and the heat generation term was used to measure the thermal conductivity and diffusivity of the eDPD system. The boundary conditions used include both the Neumann and Dirichlet boundary conditions. The Neumann boundary condition was applied via adiabatic surfaces and surfaces exposed to convection heat transfer. The DPD simulations were compared to analytical solutions and finite-difference techniques. It was found that DPD appropriately predicts the temperature distribution in the conduction regime. Details of boundary condition implementation and thermal diffusivity measurement are also described in this paper.