Transport Coefficients Of Ar—Kr Mixtures by Molecular Dynamics Computer Simulation

Abstract

Equilibrium molecular dynamics computer simulations have been used to determine the transport coefficients of model Ar—Kr mixtures, which are represented by Lennard-Jones pair potentials with Lorentz—Berthelot rules for the cross-species interactions. The component self-diffusion and mutual-diffusion coefficients are calculated from time correlation functions and mean square displacements. Time correlation functions are used to evaluate the shear and bulk viscosity, thermal conductivity and the thermal diffusion coefficient (Soret/Dufour coefficient). In the case of the thermal transport coefficients, the partial enthalpy of the two species is required at each state point to define the heat flux rigorously. We obtain this and the partial volume (and species resolved chemical potential) using particle-exchange (and particle insertion) techniques implemented in separate [NPT] simulations at the same state point.

The viscoelasticity of the fluids is characterised by the relaxation times for bulk and shear stress relaxation. The results are for dense liquids close to the triple point temperature and density. Agreement with experiment and previous simulation is particularly good for the density of the mixtures, the shear modulus, shear viscosity, shear stress relaxation time and thermal conductivity. As for the single component noble gas fluids (simulated and experiment) there is a significant qualitative difference in the temperature and, for mixtures, composition dependence of the bulk viscosity.     

Thermal diffusion in Lennard–Jones fluids in the frame of the law of the corresponding states

This work is related to the definition of a reduced thermal diffusion coefficient thanks to numerical microscale molecular dynamics simulations. This cross transport process, also called Soret effect, couples mass flux and thermal gradient and is still largely misunderstood. For this study, we have applied a boundary driven non-equilibrium molecular dynamics algorithm on Lennard–Jones spheres mixtures. Simulations have been performed at a constant reduced supercritical state, using a van der Waals’ one fluid approximation in order to fulfil the law of the corresponding states. In binary mixtures, we have studied the molecular parameters and the molar fraction influences on thermal diffusion separately and then combined. It is shown that, on pressure and on thermal conductivity, the corresponding states law is fulfilled for a wide range of molecular parameters ratios. In this frame, we have then constructed simple correlations which relate thermal diffusion factor to the mixture parameters. Combining the relations obtained, a reduced thermal diffusion factor taking into account all the various contributions has been defined. Finally, it is shown that this relation enables us to estimate thermal diffusion in various binary and ternary mixtures of Lennard–Jones spheres representing alkanes with a maximum deviation of 15%.     

On the accuracy of non-equilibrium transport coefficients calculation

The algorithms of non-equilibrium transport coefficients calculation in reacting gas mixtures are discussed. The influence of the molecular interaction potential on the transport properties is estimated in the various temperature ranges.     

Molecular dynamics simulations of a chemical reaction; conditions for local equilibrium in a temperature gradient

We have examined a simple chemical reaction in a temperature gradient; 2F <==> F2. A mechanical model was used, based on Stillinger and Weber's 2- and 3-body potentials. Equilibrium and non-equilibrium molecular dynamics simulations showed that the chemical reaction is in local thermodynamic as well as in local chemical equilibrium (delta(r)G = 0) in the supercritical fluid, for temperature gradients up to 10(12) K m(-1). The reaction is thus diffusion-controlled. The velocity distributions of both components were everywhere close to being Maxwellian. The peak distributions were shifted slightly up or down from the average velocity of all particles. The shift depended on the magnitude of the temperature gradient. The results support the assumption that the entropy production of the reacting mixture can be written as a product sum of fluxes and forces. The temperature gradient promotes interdiffusion of components in the stationary state, a small reaction rate and an accumulation of the molecule in the cold region and the atom in the hot region.     

Thermodynamic and Transport Properties of Two-temperature Oxygen Plasmas

Thermodynamic and transport properties of two-temperature oxygen plasmas are presented. Variation of species densities, mass densities, specific heat, enthalpy, viscosity, thermal conductivity, collision frequency and electrical conductivity as a function of temperature, pressure and different degree of temperature non-equilibrium are computed. Reactional, electronic and heavy particle components of the total thermal conductivity are discussed. To meet practical needs of fluid-dynamic simulations, temperatures included in the computation range from 300 K to 45,000 K, the ratio of electron temperature (T _e) to the heavy particle temperature (T _h) ranges from 1 to 30 and the pressure ranges from 0.1 to 7 atmospheres. Results obtained for thermodynamic equilibrium (T _e = T _h) under atmospheric pressure are compared with published results obtained for similar conditions. Observed overall agreement is reasonable. Slight deviations in some properties may be attributed to the values used for collision integral data and for the two temperature formulations used. An approach for computing properties under chemical non-equilibrium and associated deviations from two-temperature results under similar conditions are discussed.     

氮化硼纳米管热输运性能的分子动力学模拟

采用基于声子散射理论的Boltzmann-Peierls声子传输方程(BTE)和非平衡态分子动力学模拟(NEMD)方法研究了氮化硼纳米管(BNNT)的热输运性能.分析了BNNT的热力耦合效应,通过BTE与NEMD两种方法相结合,分析了温度和长度对BNNT热输运性能的影响,并应用量子修正扩大了NEMD的研究范围.结果表明:随着拉伸或压缩应变的增加,BNNT热输运性能均呈降低的趋势.通过计算声子态密度(PDOS)在理论上分析了以上结果,发现在拉伸状态下,声子模式的变化是决定BNNT热输运性能变化的主要因素;在压缩状态下,热导率变化是由于模型发生明显的屈曲变形引起的.在低温段,BNNT的热输运性能受量子效应影响最初有一个线性增加的过程,当温度超过一定值时,其开始显著地降低;当BNNT长度小于120nm时,随着长度的增加,其弹道性能逐渐减弱,但仍主要体现为弹道特征,其热导率(κ)与长度(L)基本满足κ∝Lα这一关系.     

Dependence of thermal diffusion factor of binary mixtures to the thermodynamic state by NEMD simulation

In this paper, we investigate the dependence of thermal diffusion factor and thermal conductivity to the temperature, density and mole fraction in Lennard–Jones binary mixtures of isotopes, noble gases and SF₆–noble gases by non-equilibrium molecular dynamics simulations.The results for the isotopic mixtures indicated that the density has a crucial effect on the dependence of thermal diffusion factor on the temperature. For isotope system at low density, thermal diffusion factor increased with temperature then remains constant at higher temperatures and the slope of thermal diffusion factor vs. temperature is positive while at higher density, thermal diffusion factor decreased with temperature and then fluctuate. For noble gas mixtures, thermal diffusion factor reduces with increasing of temperature and remain constant at high temperatures. For SF₆–Ar system, thermal diffusion factor has a negative slope and reduced with increasing of temperature, but remain nearly constant at high temperatures. For Xe–SF₆ thermal diffusion factor changed sign and the slope of thermal diffusion factor vs. temperature was negative. The results also show that thermal conductivity increases with temperature for all systems.The dependence of thermal diffusion factor to mole fraction of heavier component also investigated. The inverse of thermal diffusion factor versus mole fraction of heavier component is linear for isotope mixtures at thermodynamic conditions: (a) Low temperature, large mass ratio and all densities. (b) High temperature, large mass ratio and low densities. For Ne–Kr mixture, the inverse of thermal diffusion factor shows a linear dependence to the mole fraction of heavier component in moderate temperatures and all densities. For SF₆–Ar and Xe–SF₆ mixtures, the inverse of thermal diffusion factor has linear behaviour at moderate temperatures and low density and high temperature and low density, respectively.     

High thermal conductivity of graphene and structure defects:Prospects for thermal applications in graphene sheets

The utilization of thermal energy from different sources is an important development direction for conserving energy.With the development of technology,refined and rapid utilization of thermal energy is required.Traditional thermal conductive materials cannot meet the growing needs of human beings.Therefore,people pay attention to two-dimensional graphene film materials for their thermal conductivity.This review collects current modeling group of thermal transport on graphene,including non-equilibrium Green function(NEGF) theory,molecular dynamics(MD) simulations modeling and Boltzmann transport equation method.These models can well explain several phenomena of phonon transport in graphene.Furthe r,structural de fects were discussed and expounded the effect for graphene thermal conductivity,including doping,grain bounda ry and defects.Deeply understanding of defects on graphene,we can better grasp the thermal conductivity of graphene from the microscopic point of view.     

Study of the thermal diffusion behavior of alkane/benzene mixtures by thermal diffusion forced rayleigh scattering experiments and lattice model calculations

In this work the thermal diffusion behavior of binary mixtures of linear alkanes (heptane, nonane, undecane, tridecane, pentadecane, heptadecane) in benzene has been investigated by thermal diffusion forced Rayleigh scattering (TDFRS) for a range of concentrations and temperatures. The Soret coefficient ST of the alkane was found to be negative for these n-alkane/benzene mixtures indicating that the alkanes are enriched in the warmer regions of the liquid mixtures. For the compositions investigated in this work, the magnitude of the Soret coefficient decreases with increasing chain length and increasing alkane content of the mixtures. The temperature dependence of the Soret coefficient depends on mixture composition and alkane chain length; the slope of ST versus temperature changes from positive to negative with increasing chain length at intermediate compositions. To study the influence of molecular architecture on the Soret effect, mixtures of branched alkanes (2-methylhexane, 3-methylhexane, 2,3-dimethylpentane, 2,4-dimethylpentane, 2,2,3-trimethylbutane, and 2,2,4-trimethylpentane) in benzene were also investigated. Our results for the Soret coefficients show that the tendency for the alkanes to move to the warmer regions of the fluid decreases with increasing degree of branching. The branching effect is so strong that for 2,2,4-trimethylpentane/benzene mixtures the Soret coefficient changes sign at high alkane content and that equimolar 2,2,3-trimethylbutane/benzene mixtures have positive Soret coefficients in the investigated temperature range. In order to investigate the effect of molecular interactions on thermal diffusion, we adapted a recently developed two-chamber lattice model to n-alkane/benzene mixtures. The model includes the effects of chain-length, compressibility, and orientation dependence of benzene-benzene interactions and yields good qualitative predictions for the Soret effect in n-alkane/benzene mixtures. For the branched isomers, we find some correlations between the moments of inertia of the molecules and the Soret coefficients. PACS numbers: 66.10.Cb, 61.25.Hq.     

Thermophysical Properties of High Temperature Reacting Mixtures of Carbon and Water in the Range 400–30,000?K and 0.1–10?atm. Part 2: Transport Coefficients

The present contribution is continuation of Part 1: Equilibrium composition and thermodynamic properties. This paper is devoted to the calculation of transport properties of mixtures of water and carbon at high temperature. The transport properties, including electron diffusion coefficient, viscosity, thermal conductivity, and electrical conductivity are obtained by using the Chapman?CEnskog method expanded to the third-order approximation (second-order for viscosity), taking only elastic processes into account. The calculations, which assume local thermodynamic equilibrium, are performed for atmospheric pressure plasmas in the temperature range from 400 to 30,000?K for pressures of 0. 10, 1.0, 3.0, 5.0 and 10.0?atm. with the results obtained are compared to those of previously published studies, and the reasons for discrepancies are analyzed. The results provide reliable reference data for simulation of plasmas in mixtures of carbon and water.     

Thermodynamic and transport properties of argon/carbon and helium/carbon mixtures in fullerene synthesis

The equilibrium composition and thermodynamic and transport properties of argon; carbon and helium/carbon mixtures are calculated in the temperature range 300–20,000 K. The curves for the composition of mixtures of 50%, carbon in argon or helium are shown fir a pressure of 1.33 × 10⁴ Pa. The calculations for the heat capacity at constant pressure (C_p) and transport coefficients are validated with other studies, for the cases or pure argon and pure helium at a pressure of 10⁵ Pa. The properties of mixtures with various proportions of carbon in argon and helium are calculated. Results are presented at pressures of 105 and 1.33 × 10⁴ Pa, typical of reactors for the synthesis of fullerenes and nanotubes. It is observed that the properties of carbon and mixtures of carbon with a buffer gas (argon or helium) are very different from those of the buffer gas, thus the need to consider this effect in simulations. In general, the mixtures follow trends intermediate to those of the pure gases from which they are composed except for the thermal conductivity which shows a deviation from this tendency in the region between 11,500 and 19,000 K for argon/carbon mixtures and between 8,000 and 12,000 K for helium/carbon mixtures. Also, the electrical conductivity of mixtures of low carbon concentration is very close to that ofpure carbon. A datafile containing the transport properties of mixtures for pressures between 10⁴ and 10⁵ Pa is available free of charge from the authors.     

Cause and effect reversed in non-equilibrium molecular dynamics: an easy route to transport coefficients

A novel non-equilibrium method for calculating transport coefficients is presented. It reverses the experimental cause-and-effect picture, e.g. for the calculation of viscosities: the effect, the momentum flux or stress, is imposed, whereas the cause, the velocity gradient or shear rates, is obtained from the simulation. It differs from other Norton-ensemble methods by the way, in which the steady-state fluxes are maintained. This method involves a simple exchange of particle momenta, which is easy to implement and to analyse. Moreover, it can be made to conserve the total energy as well as the total linear momentum, so no thermostatting is needed. The resulting raw data are robust and rapidly converging. The method is tested on the calculation of the shear viscosity, the thermal conductivity and the Soret coefficient (thermal diffusion) for the Lennard–Jones (LJ) fluid near its triple point. Possible applications to other transport coefficients and more complicated systems are discussed.     

在确定的热力学状态下测量气体导热系数与黏度

气体的导热系数和黏度是重要的热物性参数,其数值大小取决于所处的热力学状态。在目前的导热系数和黏度主要测量方法中,待测工质在测量时需经历非定常的过程或处于具有物性梯度的非平衡态之下,使得待测工质的物性在时间或者空间上不处于一个确定的热力学状态。本文利用圆柱定程干涉法,通过分析气体导热系数和黏度导致的声波能量耗散,结合气体输运理论中对稀疏气体的描述,探索了在确定的热力学状态下同时测量气体导热系数和黏度的方法,并以氩(Ar)为例进行了实验验证。测量结果与已有文献一致性较好,初步证实了方法的可行性。     

Thermophysical Properties of High-Temperature Reacting Mixtures of Carbon and Water in the Range 400–30,000 K and 0.1–10 atm. Part 1: Equilibrium Composition and Thermodynamic Properties

This paper is devoted to the calculation of the chemical equilibrium composition and thermodynamic properties of reacting mixtures of carbon and water at high temperature. Equilibrium particle concentrations and thermodynamic properties including mass density, molar weight, entropy, enthalpy and specific heat at constant pressure, sonic velocity, and heat capacity ratio are determined by the method of Gibbs free energy minimization, using species data from standard thermodynamic tables. The calculations, which assume local thermodynamic equilibrium, are performed in the temperature range from 400 to 30,000 K for pressures of 0.10, 1.0, 3.0, 5.0 and 10.0 atm. The properties of the reacting mixture are affected by the possible occurrence of solid carbon formation at low temperature, and therefore attention is paid to the influence of the carbon phase transition by comparing the results obtained with and without considering solid carbon formation. The results presented here clarify some basic chemical process and are reliable reference data for use in the simulation of plasmas in reacting carbon and water mixtures together with the need of transport coefficients computation.     

Transport in molten LiF-NaF-ZrF4 mixtures: A combined computational and experimental approach

The transport properties of several LiF-NaF-ZrF₄ mixtures have been determined. Our work primarily consisted in the determination of the electrical conductivity from experimental measurements and from computer simulations. A good agreement was observed between both approaches. The simulations are based on the molecular dynamics technique and they employ a polarizable interaction potential, which was parameterized from first-principles calculations only. The diffusion coefficients were also determined from the simulations, which allowed us to understand the mechanisms responsible for the variations of electrical conductivity with temperature and composition of the melt.     

双笼氟化富勒烯分子C_(20)F_(18)(CO)_2C_(20)F_(18)电子输运性能的第一性原理研究

以双笼氟化富勒烯C_(20)F_(18)(CO)_2C_(20)F_(18)为中心分子,与Ag(100)纳米线电极连接构筑分子电子器件,通过第一性原理和非平衡格林函数相结合的方法,对器件的电子输运特性进行了研究.结果显示,在外加偏压的作用下,中心分子的前线轨道逐渐定域在分子的左侧,电子透射通道被阻塞,所对应的共振隧穿峰被压制,器件的电流-电压特性曲线在0.3~0.8V区间内表现出明显的负微分电阻(NDR)现象.     

Prediction of the transport properties of a polyatomic gas

An ab initio molecular potential model is employed in this paper to show its excellent predictability for the transport properties of a polyatomic gas from molecular dynamics simulations. A quantum mechanical treatment of molecular vibrational energies is included in the Green and Kubo integral formulas for the calculation of the thermal conductivity by the Metropolis Monte Carlo method. Using CO₂ gas as an example, the fluid transport properties in the temperature range of 300–1000 K are calculated without using any experimental data. The accuracy of the calculated transport properties is significantly improved by the present model, especially for the thermal conductivity. The average deviations of the calculated results from the experimental data for self-diffusion coefficient, shear viscosity, thermal conductivity are, respectively, 2.32%, 0.71% and 2.30%.     

Transport Coefficients of Two-temperature Lithium Plasma for Space Propulsion Applications

Lithium has been proposed as an attractive metal propellant for advanced electric propulsion. In our current work, transport coefficients including the viscosity, thermal conductivity, and electrical conductivity of lithium plasma under both the equilibrium and non-equilibrium conditions are calculated based on a two-temperature model. The collision integrals used in calculating the transport coefficients are significantly more accurate than values used in previous theoretical studies, resulting in more reliable values of the transport coefficients. Results are computed for different degrees of thermal non-equilibrium, i.e. the ratio of electron to heavy particle temperatures, from 1 to 15, with the electron temperature ranging from 300 to 60,000 K in a wide pressure range from 0.0001 to 100 atm. We compare our calculated results with existing published results and discrepancies are found and explained.     

Heats of transfer in the diffusion layer before the surface and the surface temperature for a catalytic hydrogen oxidation (H2 + (1/2)O2 --> H2O) reaction

The surface temperature and surface mole fractions are calculated for a catalytic hydrogen oxidation reaction over a Pt/Al2O3 catalyst pellet. The thermodynamics of irreversible processes was used in order to ensure the correct introduction of coupled heat and mass transfer. Two pathways, one using the 4 x 4 resistivity matrix and the other using a simplified effective conductivity matrix, were proven to yield equivalent results. By using expressions for the thermal diffusion coefficients, heats of transfer, and the Maxwell-Stefan diffusion coefficients given in the literature, available experimental data could be reproduced. The Dufour effect was found to be negligible for the prediction of the surface temperature. Neglecting the Soret effect would increase the predicted value of the surface temperature significantly-more than 30 K out of an average of about 400 K. It is found that the reaction rate can be used to predict the surface temperature.     

F+CH2CO的反应机理和动力学研究

用G3(MP2)方法对F与CH₂CO的反应进行研究,揭示了该反应的加成-消除机理.F原子首先与CH₂CO作用形成富能的中间体CH₂FCO^*,此加成反应为无势垒过程.富能的CH₂FCO^*可进一步发生解离或异构化反应生成各种可能的产物.其中CO和CH₂F可能为反应的主要产物.根据从头算的结果,用RRKM-TST理论计算该反应的速率常数.总包反应速率常数与温度存在弱的依赖关系,与总压力无关.