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

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

单晶硅薄膜面向热导率分子动力学研究

采用非平衡分子动力学方法(NEMD)研究了室温(300 K)下厚度为2～32 nm的单晶硅薄膜的沿膜平面方向的热导率,并使用Debye-Einstein模型对模拟温度进行了量子修正。模拟表明薄膜面向热导率小于相应的大体积值,并随膜厚度减小而减小,具有显著的尺寸效应。在模拟范围内膜面向热导率略大于其法向热导率;与声子气动力论的定性结果一致。晶体的表面弛豫和表面重构现象导致了MD模拟中体系总内能的升高。     

单晶硅薄膜法向热导率分子动力学研究

采用非平衡分子动力学方法（NEMD）研究了平均温度为 500K、厚度为 2～32nm的单晶硅薄膜的法向热导率。模拟结果表明,薄膜热导率显著低于对应温度下的体硅单晶的实验值,并随膜厚度减小以接近线性的规律减小。用声子气动力论模型的分析结果与NEMD模拟相一致,表明纳米单晶硅薄膜中声子平均自由程显著减小。     

微尺度下硅锗合金热导率的分子动力学模拟

采用非平衡分子动力学(NEMD)方法研究了微尺度下硅锗合金的热导率变化情况.结果表明,硅锗合金热导率具有明显的尺寸效应,显著小于大体积晶体的实验值;并且受边界散射效应的影响,热导率值随硅锗原子百分比的不同发生变化;同时热导率随温度的升高而增大,与实验值比较大致吻合.     

不同周期结构硅锗超晶格导热性能研究

构造了均匀、梯度、随机3种不同周期分布的硅/锗(Si/Ge)超晶格结构.采用非平衡分子动力学(NEMD)方法模拟了硅/锗超晶格在3种不同周期分布下的热导率,并研究了样本总长度和温度对热导率的影响.模拟结果表明:梯度和随机周期Si/Ge超晶格的热导率明显低于均匀周期结构超晶格;在不同的周期结构下,声子分别以波动和粒子性质输运为主;均匀周期超晶格热导率具有显著的尺寸效应和温度效应,而梯度、随机周期Si/Ge超晶格的热导率对样本总长度和温度的依赖性较小.     

扭转对碳纳米管束热物性的影响

采用非平衡分子动力学方法(NEMD)研究了扭转对碳纳米管(CNT)管束热导率的影响。CNT管束由7根(10,10) CNT组成。通过对CNT管束施加不同程度扭转,结果表明CNT管束热导率随着扭转程度增大迅速下降。采用T形法测量了扭转前后CNT管束热导率随温度的变化,实验结果与模拟预测趋势吻合。声子谱和拉曼光谱的结果表明扭转对原结构的破坏,结构有序性降低造成声子软化,最终导致了热导率的下降。     

微尺度下硅锗合金热导率的分子动力学模拟

采用非平衡分子动力学（NEMD）方法研究了微尺度下硅锗合金的热导率变化情况。结果表明,硅锗合金热导率具有明显的尺寸效应,显著小于大体积晶体的实验值;并且受边界散射效应的影响,热导率值随硅锗原子百分比的不同发生变化;同时热导率随温度的升高而增大,与实验值比较大致吻合。     

液态Co物性的平衡及非平衡分子动力学模拟

本文采用EAM作用势,通过平衡分子动力学(EMD)模拟的方法计算了Co熔体的自扩散系数、剪切粘度等物理性质.同时采用非平衡分子动力学(NEMD)方法计算了Co的剪切粘度.研究表明有关传输性质的计算是可与实验比较的,能够反映出液态Co典型的动力学特性.     

TiO2/ZnO纳米薄膜界面热导率的分子动力学模拟

采用平衡分子动力学方法及Buckingham势研究了金红石型TiO₂薄膜与闪锌矿型ZnO薄膜构筑的纳米薄膜界面沿晶面[0001](z轴方向)的热导率.通过优化分子模拟初始条件中的截断半径r_c和时间步后,计算并分析了平衡温度、薄膜厚度、薄膜截面大小对热导率的影响.研究表明,薄膜热导率受薄膜温度和厚度的影响很大,当温度由300 K升高600 K时,薄膜的热导率逐渐减小;当薄膜厚度由1.8 nm增大到5 nm时,热导率会逐渐增大;并在此基础 关键词： 热导率 分子动力学 2/ZnO纳米薄膜界面')" href="#">TiO₂/ZnO纳米薄膜界面 数值模拟     

固态金属中声子热传递的分子动力学模拟研究

固态金属中的热传递是声子和自由电子共同作用的结果。自由电子引起的热导率可以通过电导率,利用Wiedemann-Franz定律得到,声子引起的热导率目前仍然不能进行实验测量,只能借助其他方法来研究。本文采用非平衡分子动力学(NEMD)方法,用镶嵌原子方法(EAM)势能模型,模拟计算了不同厚度(1.760-10.56nm)金属镍薄膜中由于声子-声子作用引起的热导率。然后根据纳米厚度金属薄膜的热导率借助关联式推到宏观尺度下由于声子-声子作用引起的热导率。结果表明,对于纳米厚度金属薄膜,由于声子-声子作用引起的热导率比块体金属镍的热导率小一个数量级;薄膜厚度越小,声子-声子作用引起的热导率越小;对于块体金属镍,由于声子-声子作用引起的热导率约占其总热导率的33.0%左右。     

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.     

Layer and size dependence of thermal conductivity in multilayer graphene nanoribbons

Using nonequilibrium molecular dynamics method (NEMD), we have found that the thermal conductivity of multilayer graphene nanoribbons monotonously decreases with the increase of the number of layers which can be attributed to the phonon resonance effect of out-of-plane phonon modes. The reduction of thermal conductivity is proportional to the layer size, which is caused by the increase of phonon resonance. The results clearly show the dimensional evolution of thermal conductivity from quasi-one dimension to higher dimensions in graphene nanoribbons.     

Unexpected large thermal rectification in asymmetric grain boundary of graphene

We have investigated the lattice thermal transport across the asymmetric tilt grain boundary between armchair and zigzag graphene by nonequilibrium molecular dynamics (NEMD). We have observed significant temperature drop and ultra-low temperature-dependent thermal boundary resistance. More importantly, we find an unexpected thermal rectification phenomenon. The thermal conductivity and Kapitza conductance is direction-dependent. The effect of thermal rectification could be amplified by increasing the difference of temperature imposed on two sides. Our results propose a promising kind of thermal rectifier and phonon diodes based on polycrystalline graphene without delicate manipulation of the atomic structure.     

Thermal conductivity of periodic array of intermolecular junctions of silicon nanowires

Recent studies on intramolecular junctions of silicon nanowires and carbon nanotubes have revealed a wealth of intriguing phenomena. However, the thermal properties of the intramolecular junctions of silicon nanowires (SiNWs) are not yet well understood. In this study periodic arrays of intramolecular junctions with different lattice orientations are investigated, and their thermal conductivities are calculated using nonequilibrium molecular dynamics (NEMD) simulations. Different from the X-shaped and Y-shaped junctions of carbon nanotubes, no distinct jump is found in the temperature profile at the junctions. Compared with straight pristine SiNWs of the same length, the thermal conductivity of the periodic array of intramolecular junctions is reduced. The underlying mechanism of the observed behavior is analyzed by the phonon spectral density of the atomic velocities. The dependence of temperature on the thermal conductivity of this junction array structure is discussed.     

Atomistic hybrid DSMC/NEMD method for nonequilibrium multiscale simulations

A multiscale hybrid method for coupling the direct simulation Monte Carlo (DSMC) method to the nonequilibrium molecular dynamics (NEMD) method is introduced. The method addresses Knudsen layer type gas flows within a few mean free paths of an interface or about an object with dimensions of the order of a few mean free paths. It employs the NEMD method to resolve nanoscale phenomena closest to the interface along with coupled DSMC simulation of the remainder of the Knudsen layer. The hybrid DSMC/NEMD method is a particle based algorithm without a buffer zone. It incorporates a new, modified generalized soft sphere (MGSS) molecular collision model to improve the poor computational efficiency of the traditional generalized soft sphere GSS model and to achieve DSMC compatibility with Lennard-Jones NEMD molecular interactions. An equilibrium gas, a Fourier thermal flow, and an oscillatory Couette flow, are simulated to validate the method. The method shows good agreement with Maxwell–Boltzmann theory for the equilibrium system, Chapman–Enskog theory for Fourier flow, and pure DSMC simulations for oscillatory Couette flow. Speedup in CPU time of the hybrid solver is benchmarked against a pure NEMD solver baseline for different system sizes and solver domain partitions. Finally, the hybrid method is applied to investigate interaction of argon gas with solid surface molecules in a parametric study of the influence of wetting effects and solid molecular mass on energy transfer and thermal accommodation coefficients. It is determined that wetting effect strength and solid molecular mass have a significant impact on the energy transfer between gas and solid phases and thermal accommodation coefficient.     

Molecular dynamics studies on the thermal conductivity of single-walled carbon nanotubes

We have studied the thermal conductivity of single-walled carbon nanotubes (SWCNTs) using the NEMD method. The results indicate that the thermal conductivity values are not profoundly influenced by the specific simulation-technique used in the MD simulations. Some possible reasons, which could be responsible for the discrepancy on thermal conductivity values of SWCNTs in the literatures, are discussed.      

Influence of chirality on the thermal conductivity of single-walled carbon nanotubes

The influence of chirality on the thermal conductivity of single-walled carbon nanotubes(SWNTs) is discussed in this paper, using a non-equilibrium molecular dynamics(NEMD) method. The tube lengths of the SWNTs studied here are 20, 50, and 100 nm, respectively, and at each length the relationship between chiral angle and thermal conductivity of a SWNT is revealed. We find that if the tube length is relatively short, the influence of chirality on the thermal conductivity of a SWNT is more obvious and that a SWNT with a larger chiral angle has a greater thermal conductivity. Moreover, the thermal conductivity of a zigzag SWNT is smaller than that of an armchair one. As the tube length becomes longer, the thermal conductivity increases while the influence of chirality on the thermal conductivity decreases.