氩晶体薄膜法向热导率的分子动力学模拟

结合卫星“微型核”的特点，研究电介质薄膜中的导热机理以及薄膜厚度对导热系数的影响.以结构较为简单、具有可靠势能函数，实验数据较为丰富和可靠的氩的(fcc)晶体为模型,采用平衡分子动力学方法(EMD)和各向异性非平衡分子动力学方法(NEMD)计算了氩晶体及其法向薄膜的热导率,并与实验结果进行比较.模拟结果表明，氩晶体纳米薄膜的热导率显著小于对应大体积晶体的实验值，具有明显的尺寸效应.在氩薄膜厚度为2.124—5.310nm的模拟范围内,薄膜的法向热导率随着薄膜厚度的增加而呈近似线性增加. 关键词： 热导率 纳米薄膜 尺寸效应 平衡分子动力学 非平衡分子动力学     

一种模拟热导率的非平衡分子动力学方法

提出一种计算热导率的非平衡分子动力学(NEMD)方法,通过构造均匀内热源获得抛物线形温度分布,并基于Fourier导热定律计算热导率,与Müller-Plathe发展的反扰动非平衡分子动力学(RNEMD)方法相比,不仅具有能量动量守恒和收敛性好的优点,还克服了常规NEMD方法中热冷源区域存在局域热力学非平衡的问题,并有模拟系统温差影响小的特点.对液态氩的热导率进行模拟并与RNEMD方法的模拟结果进行对比.     

高剪切率的Couette流下液氩流体热导率的分子动力学模拟

基于平衡分子动力学(EMD)模拟方法,对Lees-Edwards周期性边界条件作一修正并成功地应用其形成具有高恒定剪切率的Couette流动,最后通过Green-Kubo公式计算得到了剪切流下液氩流体的热导率。计算结果,表明在剪切流下液氩流体的热导率仍然表现为各向同性,同时还发现流动方向的热导率不随流动剪切率的增加而发生变化,并且该值跟静态下液氩的热导率(λ=0.132 Wm~(-1)K~(-1))基本一致。     

水中自然超空泡机理及减阻效应的非平衡分子动力学研究

应用非平衡分子动力学方法,对水中超空泡流形成机理及减阻效应进行了模拟研究.计算得到了流体密度分布、局部空化数分布、阻力系数及含气百分比等流场细节数据,结果显示空化数判据在分子层面仍然成立,局部低空化数区域与超空泡形成区域在空间上分离;超空泡形成和稳定主要受物体运动速度影响;空化器构型对空泡内含气率有较大影响;从云雾空化状态过渡到超空泡,物体表面摩擦力可以减小50%—90%.与数值模拟结果的对比表明非平衡分子动力学模拟适用于研究微观超空泡机理,能够经济有效地探讨超高速流体运动的一些自然规律. 关键词： 非平衡分子动力学 超空泡 高速流体 摩擦系数     

LJ流体非牛顿现象的分子动力学模拟

使用非平衡分子动力学模拟方法研究了单原子LJ流体的非牛顿流变行为,并在系统中分别施加稳态Couette流场和振荡剪切流场.在Couette流场的模拟中,流体出现剪切变稀和法向应力差效应,不同剪切率下的径向分布函数反映了流体分子由于剪切所导致的微观结构变化,通过分析势能函数发现当剪切率增大时,分子间排斥作用增强,吸引作用减弱.在振荡剪切流场的模拟中,发现剪切应力和剪切率之间的相位差随频率增加而增加,随频率增加复数粘度的实部先增大再减小,虚部单调增加,导致虚部粘度相对实部粘度比例增大,弹性模量和粘性模量之比也随频率增加而增加.这三点现象表明LJ流体出现粘弹性行为,且在高频率下,弹性所占比重增大.     

纳米流体热导率和粘度的分子动力学模拟计算

本文采用分子动力学(MD)模拟来计算纳米流体比较重要的热物性:热导率和粘度,与已有实验结果比较符合 较好,为进一步研究纳米流体传热效率提供了依据。     

流场环境下复杂囊泡的动力学行为

采用非平衡态分子动力学方法研究了二维复杂囊泡在剪切流中的动力学行为. 模拟发现了复杂囊泡经典的翻滚(tumbling)、摇摆(trembling)和坦克履(tank-treading)行为, 还观察到由坦克履行为向平动行为(translating)的转变. 囊泡的平动行为与剪切率大小、复杂囊泡的形状密切相关. 当大囊泡均匀嫁接较多数目的小囊泡后, 其平动方式消失. 该研究有益于加深对囊泡在剪切流场中复杂性行为的理解.     

Si表面间水平碳纳米管束的分子动力学模拟研究

本文采用分子动力学模拟方法研究了Si表面间单壁水平碳纳米管束SWCNT (10,10)的变形和摩擦特性.系统在弛豫平衡后,首先对碳纳米管束施加压力至碳纳米管或Si表面结构破坏.之后在无压力和高压力两种情况下使上表面沿水平方向做剪切运动以研究碳纳米管束的摩擦特性.结果表明,由于碳纳米管的柔韧性,碳纳米管束在加载过程中出现明显变形,但直至3.8 GPa高压下并无结构破坏.系统无压力时SWCNT (10,10)在原地轻微随机滚动,压力为3.8 GPa时,碳纳米管束出现了整体的轻微滑动,同时伴随无规律的轻微滚动, 关键词： 碳纳米管束 摩擦 分子动力学模拟     

SiO_2气凝胶导热机理的分子动力学模拟研究

本文采用分子动力学方法,选用Lennard-Jones势函数,对氩分子在SiO2气凝胶骨架结构间空隙中的运动进行了平衡分子动力学模拟(EMD)和非平衡分子动力学模拟(NEMD),分析其能量输运特性.结果表明:在骨架空隙尺寸为15nm条件下,气凝胶内部氩分子的扩散系数仅为同等条件下自由空间内氩分子扩散系数的50.35%,...     

Cu-Co双金属团簇结构演化及其性质的分子动力学模拟

采用分子动力学结合嵌入原子方法对比研究了Co分布于Cu-Co团簇不同层的结构和性质.研究表明:Co原子分层掺杂可对团簇的结构转变点和熔点进行诱导控制;分层掺杂的Cu-Co团簇第一相变是一种扩散度较小的由立方八面体转变为二十面体的相变;Co原子易于向低能态团簇的亚表层(111)面偏析,从而诱导团簇结构紊乱,造成其熔点差异.     

Molecular Dynamics Simulations for the Shear Viscosity of the One‐Component Plasma

We discuss two methods for determining the shear viscosity of a fluid of particles with Yukawa interaction potential (a one‐component plasma). Both methods are based on computing the equilibrium dynamics using large‐scale molecular dynamics (MD) simulations. Our MD results illustrate that the hydrodynamic method for computing the shear viscosity is feasible and therefore complements the more widely used method based on the Green‐Kubo relation. We expect that in the future our shear viscosity calculations will be used to assist with the interpretation and analysis of x‐ray scattering experiments, which could in principle measure this fundamental dynamical quantity (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Many body forces between long conducting molecules

Non-equilibrium molecular dynamics is used to calculate the wavevector and strain rate dependence of shear viscosity in a soft sphere fluid. The calculations are consistent with a non-analytic functional dependence of viscosity with wavevector. A consequence of such non-analyticity would be that the linear and non-linear Burnett coefficients for viscosity would not exist for fluids in three dimensions. The calculations also show the numerical consistency of three different non-equilibrium simulation methods for calculating shear viscosity.     

Nonequilibrium molecular dynamics simulation of shear viscosity by a uniform momentum source-and-sink scheme

A uniform momentum source-and-sink scheme of nonequilibrium molecular dynamics (NEMD) is developed to calculate the shear viscosity of fluids in this paper. The uniform momentum source and sink are realized by momentum exchanges of individual atoms in the left and right half systems, like the reverse nonequilibrium molecular dynamics (RNEMD) method [20] [Müller-Plathe, Phys. Rev. E, 49 (359), 1999]. This method has all features of RNEMD. In addition, the present momentum swap strategy maximizes the perturbation relaxation and eliminates the boundary jumps, which often harm other NEMD methods greatly. With periodic boundary conditions quadratic velocity profiles can be constructed and from the mean velocities of the right and left half systems the shear viscosity can be easily extracted. The scheme is tested on Lennard-Jones fluids over a wide range of state points (temperature and density), momentum exchange intervals and system sizes. It is demonstrated that the present approach can give reliable results with fast convergence by properly selecting the simulation parameters, i.e. particle number and exchange interval.     

The three body interaction energy between optically active molecules

An account is given of a steady-state molecular dynamics method for the calculation of the shear viscosity of a dense fluid. The method is applied to the case of a system described by a Lennard-Jones potential with parameters appropriate to argon. Satisfactory agreement with experimental data is found over a range of temperature and density for which the viscosity varies by a factor of approximately 8. The results are interpreted in terms of a modified hard-sphere model and it is shown that the Lennard-Jones fluid obeys a Stokes-type relation in which the effective molecular diameter is close to that deduced from equilibrium structural properties.     

Molecular origin of shear thickening in transient polymer networks: A molecular dynamics study

This paper proposes a simple model of transient networks of telechelic associating polymers for molecular simulations and reports the main results obtained by molecular dynamics on the rheological properties of the transient networks. The steady shear viscosity obtained by the non-equilibrium molecular dynamics simulation exhibits shear thickening at moderate shear rates and shear thinning at larger shear rates. The behavior is similar to that observed in experiments of telechelic associating polymers. By analyzing the distribution function of the end-to-end vector of bridge chains as a function of the shear rate, we find that shear thickening is mainly caused by the stress from the bridge chains highly stretched by shear flow. We also find that fracture of the transient network occurs in the shear-thinning regime at high shear rates.     

Dynamic Shear Viscosity in a 2D Yukawa System

We present non‐equilibrium molecular dynamics simulation studies on the dynamic (complex) shear viscosity of a 2D many‐particle system interacting via a Yukawa (Debye‐Hückel) type inter‐particle potential. Our investigations reveal the complex interplay of dissipative and elastic processes, as well as the effect of single particle resonances and enhanced collective excitations, and the influence of the external forces on the structural correlations in the system (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Shear viscosity of strongly coupled Yukawa systems on finite length scales

The Yukawa shear viscosity has been calculated using nonequilibrium molecular dynamics. Near the viscosity minimum, we find exponential decay consistent with the Navier-Stokes equation, with significant deviations on finite length scales for larger viscosity values. The viscosity is determined to be nonlocal on a scale length consistent with the correlation length, revealing the length scales necessary for obtaining transport coefficients in the hydrodynamic limit by nonequilibrium molecular dynamics methods. Our results are quasiuniversal with respect to excess entropy for excess entropies well below unity.     

Structure,thermodynamic and transport properties of imidazolium-based bis(trifluoromethylsulfonyl)imide ionic liquids from molecular dynamics simulations

Molecular dynamics (MD) simulations have been performed in order to investigate the properties of [C _n mim⁺][Tf₂N^?] (n?=?4,?8,?12) ionic liquids (ILs) in a wide temperature range (298.15?498.15?K) and at atmospheric pressure (1 bar). A previously developed methodology for the calculation of the charge distribution that incorporates ab initio quantum mechanical calculations based on density functional theory (DFT) was used to calculate the partial charges for the classical molecular simulations. The wide range of time scales that characterize the segmental dynamics of these ILs, especially at low temperatures, required very long MD simulations, on the order of several tens of nanoseconds, to calculate the thermodynamic (density, thermal expansion, isothermal compressibility), structural (radial distribution functions between the centers of mass of ions and between individual sites, radial-angular distribution functions) and dynamic (relaxation times of the reorientation of the bonds and the torsion angles, self-diffusion coefficients, shear viscosity) properties. The influence of the temperature and the cation's alkyl chain length on the above-mentioned properties was thoroughly investigated. The calculated thermodynamic (primary and derivative) and structural properties are in good agreement with the experimental data, while the extremely sluggish dynamics of the ILs under study renders the calculation of their transport properties a very complicated and challenging task, especially at low temperatures.     

Flow, ordering, and jamming of sheared granular suspensions

We study the rheological properties of a granular suspension subject to constant shear stress by constant volume molecular dynamics simulations. We derive the system "flow diagram" in the volume fraction or stress plane (phi, F): at low phi the flow is disordered, with the viscosity obeying a Bagnold-like scaling only at small F and diverging as the jamming point is approached; if the shear stress is strong enough, at higher phi an ordered flow regime is found, the order-disorder transition being marked by a sharp drop of the viscosity. A broad jamming region is also observed where, in analogy with the glassy region of thermal systems, slow dynamics followed by kinetic arrest occurs when the ordering transition is prevented.     

Transient time correlation function calculation of the viscosity of a molecular fluid at low shear rates: a comparison of stress tensors

Shear viscosity of n-decane was computed using the molecular stress transient time correlation function (TTCF) formalism for the wide range of shear rates from 1.7 × 10¹⁰ s^?1 to 2.13 × 10⁴ s^?1. In earlier work calculations were presented for the shear viscosity of n-decane using the atomic stress formalism of the TTCF method (G. Pan and C. McCabe, J. Chem. Phys. 125(19), 4527 (2006)) in which we were able to close the gap between the lowest shear rates accessible by direct nonequilibrium molecular dynamics (NEMD) simulations and the highest shear rates possible in experimental studies. Here it is shown that the application of the molecular stress approach within the TTCF formalism, as an alternative to the atomic stress method, significantly reduces the number of NEMD trajectories necessary to obtain the shear viscosity.