Abstract

The homogeneous non-equilibrium molecular dynamics (NEMD) heat flow algorithm due to Evans is known to be unstable for large systems (N ⩾ 896) at sufficiently large values of the applied heat field. Instead of heat being conducted homogeneously through the system, the system transfers heat by generating a solitary supersonic shock wave. As a result the effective thermal conductivity increases dramatically. This instability reduces the usefulness of the heat flow NEMD algorithm by increasing the difficulty of making the zero field extrapolation. This extrapolation is necessary in order to calculate the Navier-Stokes thermal conductivity. We present an algorithm which strongly inhibits the generation of the shock wave and results in a conductivity which is largely independent of system size allowing a more efficient extrapolation to zero field.

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call

Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.