Heat conduction and thermomolecular orientation in diatomic fluids: a non-equilibrium molecular dynamics study

2012 ◽  
Vol 38 (14-15) ◽  
pp. 1198-1208 ◽  
Author(s):  
Frank Römer ◽  
Fernando Bresme
Author(s):  
Quan-Wen Hou ◽  
Bing-Yang Cao ◽  
Zeng-Yuan Guo

The phonon relaxation and heat conduction of the Femi-Pasta-Ulam β lattice are studied via molecular dynamics simulations. The phonon relaxation rate is calculated from the energy autocorrelation function for different modes at various temperatures through equilibrium molecular dynamics simulations. The relaxation rate as a function of wave vector k is estimated to be proportional to k1.688, which leads to a N0.41 divergence of the thermal conductivity in the framework of Green-Kubo relation. This result is in agreement with that obtained by non-equilibrium molecular dynamics simulations which estimate the length dependence exponent of thermal conductivity as 0.415. Our results confirm the N2/5 divergence in one-dimensional FPU β lattice. The effect of the heat flux on the thermal conductivity is also studied by imposing large temperature differences on the two ends of the lattice in non-equilibrium molecular dynamics simulations. The results indicate that the thermal conductivity is insensitive to the heat flux under our simulation conditions, and the linear response theory is widely applicable.


2020 ◽  
Vol 45 (4) ◽  
pp. 319-332
Author(s):  
Xiaoyu Chen ◽  
Ruquan Liang ◽  
Yong Wang ◽  
Ziqi Xia ◽  
Lichun Wu ◽  
...  

AbstractThe effect of the temperature gradient on the Soret coefficient in n-pentane/n-decane (n-C5/n-C10) mixtures was investigated using non-equilibrium molecular dynamics (NEMD) with the heat exchange (eHEX) algorithm. n-Pentane/n-decane mixtures with three different compositions (0.25, 0.5, and 0.75 mole fractions, respectively) and the TraPPE-UA force field were used in computing the Soret coefficient ({S_{T}}) at 300 K and 1 atm. Added/removed heat quantities (ΔQ) of 0.002, 0.004, 0.006, 0.008, and 0.01 kcal/mol were employed in eHEX processes in order to study the effect of different thermal gradients on the Soret coefficient. Moreover, a phenomenological description was applied to discuss the mechanism of this effect. Present results show that the Soret coefficient values firstly fluctuate violently and then become increasingly stable with increasing ΔQ (especially in the mixture with a mole fraction of 0.75), which means that ΔQ has a smaller effect on the Soret coefficient when the temperature gradient is higher than a certain thermal gradient. Thus, a high temperature gradient is recommended for calculating the Soret coefficient under the conditions that a linear response and constant phase are ensured in the system. In addition, the simulated Soret coefficient obtained at the highest ΔQ within three different compositions is in great agreement with experimental data.


2015 ◽  
Vol 32 (3) ◽  
pp. 683-698 ◽  
Author(s):  
F. A. Furtado ◽  
◽  
A. J. Silveira ◽  
C. R. A. Abreu ◽  
F. W. Tavares ◽  
...  

Sign in / Sign up

Export Citation Format

Share Document