Heat induced instability in a model liquid
- 1 September 1989
- journal article
- research article
- Published by Taylor & Francis in Molecular Physics
- Vol. 68 (1) , 97-105
- https://doi.org/10.1080/00268978900101991
Abstract
The stability of the homogeneous nonequilibrium molecular dynamics (NEMD) heat flow algorithm due to Evans is discussed with respect to a two dimensional soft disc inverse-twelve N-particle liquid. The thermal conductivity coefficient is evaluated as a function of an induced heat field, related to ∇ ln T, for N = 56, 224, 896 and 3584. We show that homogeneous heat conduction becomes unstable for large systems (N ⩾ 896) at sufficiently large values of the field; the conductivity increases and the system conducts heat by generating a solitary shock wave. The results are strongly dependent on N and differ in this respect from parallel NEMD studies on a liquid subjected to a shear.This publication has 7 references indexed in Scilit:
- Local order in a dense liquidPhysical Review A, 1988
- Shear induced anisotropy in two-dimensional liquidsPhysica A: Statistical Mechanics and its Applications, 1988
- Shear-induced angular dependence of the liquid pair correlation functionPhysical Review A, 1987
- Thermal conductivity of the Lennard-Jones fluidPhysical Review A, 1986
- Shear-induced melting and reentrant positional ordering in a system of spherical particlesInternational Journal of Thermophysics, 1985
- The calculation of thermal conductivities by perturbed molecular dynamics simulationJournal of Physics C: Solid State Physics, 1983
- Stationary nonequilibrium states by molecular dynamics. Fourier's lawPhysical Review A, 1982