Modelling transient heat conduction in solids at multiple length and time scales: A coupled non-equilibrium molecular dynamics/continuum approach

详细信息    查看全文

• 作者：Kenny Jolley ; Simon P.A. Gill
• 关键词：Molecular dynamics ; Boundary conditions ; Heat transfer ; Multiscale modelling ; Atomistic/continuum coupling
• 刊名：Journal of Computational Physics
• 出版年：2009
• 出版时间：20 October 2009
• 年：2009
• 卷：228
• 期：19
• 页码：7412-7425
• 全文大小：749 K

文摘

A method for controlling the thermal boundary conditions of non-equilibrium molecular dynamics simulations is presented. The method is simple to implement into a conventional molecular dynamics code and independent of the atomistic model employed. It works by regulating the temperature in a thermostatted boundary region by feedback control to achieve the desired temperature at the edge of an inner region where the true atomistic dynamics are retained. This is necessary to avoid intrinsic boundary effects in non-equilibrium molecular dynamics simulations. Three thermostats are investigated: the global deterministic Nosé–Hoover thermostat and two local stochastic thermostats, Langevin and stadium damping. The latter thermostat is introduced to avoid the adverse reflection of phonons that occurs at an abrupt interface. The method is then extended to allow atomistic/continuum models to be thermally coupled concurrently for the analysis of large steady state and transient heat conduction problems. The effectiveness of the algorithm is demonstrated for the example of heat flow down a three-dimensional atomistic rod of uniform cross-section subjected to a variety of boundary conditions.