CO_2水合物导热特性的分子动力学模拟

使用分子动力学模拟方法在NVT系综下对结构完整CO_2水合物以及结构缺陷CO_2水合物进行了导热模拟计算.对于结构完整的CO_2水合物,在200-230 K温度区间内,体系导热系数由0.4684 W·m~(-1)·K~(-1)变化到0.4836 W·m~(-1)·K~(-1),温度相关性较弱;而在230-280 K温度区间内,体系导热系数由0.4836 W·m~(-1)·K~(-1)变化到0.7494 W·m~(-1)·K~(-1),温度相关性变强;另外,通过计算功率图谱发现主体分子对水合物体系的导热贡献更大.对于结构缺陷CO_2水合物,发现晶穴占有率和笼形结构缺陷对体系导热均有一定影响,空笼晶胞导热系数约为完整晶胞导热系数的86.67%,体系的导热能力主要取决于主体结构的性质.     

温度依变性对纳米颗粒悬浮液导热系数测量影响

本文报告了50 nm氧化铜(体积百分比φ≤0.6%)与去离子水的悬浮液有效导热系数测量结果,讨论了温度依变性对导热系数的影响,运用适用于低浓度悬浮液的能量和质量方程进行分析表明,颗粒空间分布和温度场间相互影响,进而影响导热系数测量结果.考虑低浓度悬浮液中颗粒布朗扩散和热泳作用,讨论了运用准稳态方法测量纳米颗粒悬浮液导热系数的有效性,即需同时保证合适的加热热流密度和液体试样两侧的温差.     

CO2水合物生长特性的分子动力学模拟

采用分子动力学对CO2水合物生长进行模拟并分析其变化规律,探究了CO2水合物晶体生长的微观过程以及不同温度、压力条件对水合物生长特性的影响机理. CO2水合物的生长是从已有晶胞附近向外扩散并逐具有序性的过程,水分子间逐渐形成四面体氢键与CO2共同形成完整的水合物笼. 另外, CO2水合物生长需要合适的过冷度,在压力为30MPa、温度范围265K至275K,水合物笼型晶胞可正常生成,并且温度越低,生成速率越快;在高温290K和300K时,体系分子运动加剧,水合物笼直接散开. 此外,在温度为270K、不同压力条件下,发现相较温度而言,由于CO2溶解度随压力变化的不明显性,导致CO2水合物增长速度伴随压力的变化相对不敏感.     

纳米变截面结构的导热性质

采用分子动力学方法模拟了固态氩的纳米变截面结构的导热性质,研究发现纳米变截面材料的热阻和热流的大小与方向都相关:当纳米结构沿热流的方向为渐缩时,纳米结构的热阻随热流的增加而增大,而当纳米结构沿热流的方向为渐扩时,纳米结构的热阻随热流的增加呈减小的趋势;当热流较大时,热流沿渐缩方向时的热阻明显大于热流沿渐扩方向时的热阻,但当热流较小时纳米变截面结构的热阻和热流方向的关系不大.最后依据热质的运动和传递理论的动能效应对该现象进行了分析解释.     

聚乙烯/银纳米颗粒复合物的分子动力学模拟研究

通过分子动力学模拟对聚乙烯/银纳米颗粒复合物的结构、极化率和红外光谱、热力学性质、力学特性进行计算, 分析其随模拟温度和银颗粒尺寸的变化规律. 模拟结果表明: 聚乙烯/银纳米颗粒复合物为各向同性的无定形结构, 温度升高可提高银纳米颗粒的分散均匀性; 银纳米颗粒表面多个原子层呈现无定形状态, 并在银颗粒和聚乙烯基体的界面形成电极化层, 界面区域随颗粒尺寸和温度的增加分别减小和增加; 与聚乙烯体系相比, 聚乙烯/银纳米颗粒复合物的极化率高很多, 且随温度的升高和银颗粒尺寸的减小而增大; 银颗粒尺寸直接影响界面电偶极矩的强度和振动频率, 红外光谱峰强度和峰位随颗粒尺寸发生变化; 聚乙烯/银纳米颗粒复合物具有比聚乙烯体系更高的等容热容和与聚乙烯体系相反的负值热压力系数, 热容随颗粒尺寸的变化较小, 但随温度的升高而明显减小, 具有显著的温度效应; 热压力系数随温度的变化较小, 但随颗粒尺寸的增加而减小, 具有明显的尺度效应, 温度稳定性更好; 聚乙烯/银纳米颗粒复合物的力学特性表现出各向同性材料的弹性常数张量, 具有比聚乙烯体系更高的杨氏模量和泊松比, 并且都随温度的升高和银颗粒尺寸的增大而减小, 加入银纳米颗粒可有效改善聚乙烯的力学性质. 关键词： 分子动力学模拟 聚合物纳米复合物 纳米颗粒     

纳米颗粒聚集形态对纳米流体导热系数的影响

纳米流体中悬浮的纳米颗粒可以增强其导热性能已经得到广泛认可,然而纳米流体颗粒增强传热的机理目前尚不清楚.研究表明,纳米颗粒的聚集是纳米流体导热系数增大的重要机制,而且纳米颗粒聚集的形态对纳米流体的导热系数有重要影响,但是目前的导热系数模型大多是建立在Maxwell有效介质理论的"静态"和"均匀分散"假设基础上.本文用平衡分子动力学模拟Cu-Ar纳米流体,采用Green-Kubo公式计算导热系数,采用Schmidt-Ott关系式计算不同聚集形态下的分形维数.对比导热系数与分形维数可以发现:在相同体积分数下,较低的分形维数会有更高的导热系数,分析了分形维数与导热系数的定量关系.此外,通过径向分布函数可以看出纳米颗粒紧密聚集与松散聚集的差异,基液分子在纳米颗粒附近的纳米薄层中处于动态平衡状态.研究结果有助于理解纳米颗粒聚集形态对导热系数的影响机理.     

碳纳米管表面金纳米颗粒的形成与结构转变

利用分子动力学模拟研究了室温下金纳米颗粒在碳纳米管表面的结构和作用能.研究结果表明，金纳米颗粒随着尺寸的增大会发生不同于孤立状态下的结构转变.当原子数小于130时，颗粒属于无序结构；当原子数大于140时，呈现面心立方晶体结构.小金纳米颗粒和碳纳米管结合紧密，相互作用能正比于面对碳纳米管的颗粒表面面积.     

碳纳米管表面金纳米颗粒的形成与结构转变

利用分子动力学模拟研究了室温下金纳米颗粒在碳纳米管表面的结构和作用能.研究结果表明,金纳米颗粒随着尺寸的增大会发生不同于孤立状态下的结构转变.当原子数小于130时,颗粒属于无序结构;当原子数大于140时,呈现面心立方晶体结构.小金纳米颗粒和碳纳米管结合紧密,相互作用能正比于面对碳纳米管的颗粒表面面积. 关键词： 金纳米颗粒 碳纳米管 分子动力学模拟     

硅纳米颗粒在碳纳米管表面生长的分子动力学模拟

采用分子动力学模拟方法研究了硅纳米颗粒在碳纳米管上的生长，并分析了这种复合材料的基本结构.研究表明，由于硅原子和碳纳米管之间的相互作用以及碳纳米管的巨大的表面曲率，硅原子在碳纳米管表面不是形成覆盖碳纳米管的二维薄膜，而是生成具有三维结构的硅纳米颗粒.小纳米颗粒的结构和无基底条件下生成的颗粒结构基本一致.对于大纳米颗粒，不同于无基底条件下形成的球状纳米晶体硅结构，硅纳米颗粒沿管轴方向伸长，其结构为类似于硅晶体的无定形网络结构. 关键词： 纳米颗粒 碳纳米管 硅 分子动力学模拟     

聚酰亚胺/铜纳米颗粒复合物的分子动力学模拟研究

通过分子动力学模拟对聚酰亚胺/铜纳米颗粒复合物的形态结构、 热力学性质、力学特性进行计算, 分析其随模拟温度和纳米颗粒尺寸的变化规律. 模拟结果表明, 聚酰亚胺/铜纳米颗粒复合物为各向同性的无定形态结构, 铜纳米颗粒与聚酰亚胺基体之间通过较强的范德华作用结合在一起使结构更加稳定, 铜纳米颗粒表面多个原子层呈现无定形状态, 在铜颗粒和聚酰亚胺基体之间形成界面层, 界面区域随颗粒尺寸和温度的增加分别减小和增加. 聚酰亚胺/铜纳米颗粒复合物的等容热容随着颗粒尺寸增大而明显增高, 随温度变化比聚酰亚胺体系更为缓慢, 在较低温度下较小颗粒尺寸复合物的热容比聚酰亚胺体系更低. 聚酰亚胺/铜纳米颗粒复合物的热压力系数随颗粒尺寸增加而显著增大, 比聚酰亚胺体系的热压力系数更小, 且随温度升高而减小的程度要小得多. 聚酰亚胺/铜纳米颗粒复合物的热力学性质表现出明显的尺度效应, 温度稳定性明显高于聚酰亚胺体系. 聚酰亚胺/铜纳米颗粒复合物的力学特性表现出各向同性材料的弹性常数张量, 具有比聚酰亚胺体系更低的杨氏模量和泊松比, 随温度升高分别减小和增大, 与聚酰亚胺体系随温度的变化趋势相反, 且杨氏模量的温度稳定性显著提高, 同时泊松比随纳米颗粒尺寸增大而减小, 具有明显的尺度效应. 加入铜纳米颗粒形成复合物可获得与聚酰亚胺体系显著不同的力学新特性. 关键词： 分子动力学模拟 聚合物纳米复合物 聚酰亚胺 纳米颗粒     

Molecular dynamics simulation of decomposition and thermal conductivity of methane hydrate in porous media

The hydrate has characteristics of low thermal conductivity and temperature sensitivity. To further analysis the mechanism of thermal conductivity and provide method for the exploitation, transportation and utilization of hydrate, the effect of decomposition and thermal conductivity of methane hydrate in porous media has been studied by using the molecular dynamics simulation. In this study, the simulation is carried out under the condition of temperature 253.15 K-273.15 K and pressure 1 MPa. The results show that the thermal conductivity of methane hydrate increases with the increase of temperature and has a faster growth near freezing. With the addition of porous media, the thermal conductivity of the methane hydrate improves significantly. The methane hydrate-porous media system also has the characteristics of vitreous body.With the decrease of the pore size of the porous media, thermal conductivity of the system increases gradually at the same temperature. It can be ascertained that the porous media of different pore sizes have strengthened the role of the thermal conductivity of hydrates.     

四丁基溴化铵水合物浆体导热系数测量研究

四丁基溴化铵(TBAB)水合物浆体在常压下的相变温度介于0-12℃之间,作为蓄冷材料使用时由于相变过程的存在使得其蓄冷能力较高,而且在管道中具有良好的流动特性,因而是一种理想的蓄冷和冷量输送材料。对比分析了传统导热系数计算公式和基于一维非稳态导热模型导出的导热系数计算公式的区别。利用热线装置分别测量了TBAB溶液和水合物浆体的导热系数。实验得出5-30 wt%TBAB溶液的导热系数在0.4-0.6 W.m-1.K-1之间,并随浓度的增加而减小;10-40vol%的水合物浆体的导热系数在0.5-0.6 W.m-1.K-1之间,并随体积分数的增加而增大;相同体积浓度时A型水合物浆体的导热系数大于B型水合物浆体的导热系数。     

Effect of triangular vacancy defect on thermal conductivity and thermal rectification in graphene nanoribbons

We investigate the thermal transport properties of armchair graphene nanoribbons (AGNRs) possessing various sizes of triangular vacancy defect within a temperature range of 200–600 K by using classical molecular dynamics simulation. The results show that the thermal conductivities of the graphene nanoribbons decrease with increasing sizes of triangular vacancy defects in both directions across the whole temperature range tested, and the presence of the defect can decrease the thermal conductivity by more than 40% as the number of removed cluster atoms is increased to 25 (1.56% for vacancy concentration) owing to the effect of phonon–defect scattering. In the meantime, we find the thermal conductivity of defective graphene nanoribbons is insensitive to the temperature change at higher vacancy concentrations. Furthermore, the dependence of temperatures and various sizes of triangular vacancy defect for the thermal rectification ration are also detected. This work implies a possible route to achieve thermal rectifier for 2D materials by defect engineering.     

Investigation on thermal conductivity of graphene/Si heterostructure with different defect ratios and sizes

The thermal conductivity (TC) of graphene/Si heterostructures with different defect ratios and sizes was investigated using the molecular dynamics method. As the defect ratio of heterostructure increased, the TC decreased first sharply and then slowly under a high temperature stage. The TC of heterostructure also showed a significant size effect. This phenomenon was explained by phonon dispersion and flip competition. The phonon density of states for the graphene heterostructure with different defect ratios and sizes was obtained to understand the thermal transport mechanism. Analysis showed that with the increase in the defect ratio and when the flexural modes of the heterostructure became weak, the longitudinal and transverse modes gradually dominated the phonon transport. This phenomenon can be explained that the Si atom vibration was harder in the vertical plane than that of graphene. The vibration mode hindered the heat carrier of graphene and affected heat transport to the heterostructure.     

纳米悬浮液的有效导热系数

文中以有效介质近似理论为基础,考虑了纳米颗粒在基液中强烈的B rown ian运动对强化传热的作用和纳米颗粒的表面吸附液体层、纳米颗粒的粒径和体积分数对纳米悬浮液有效导热系数的影响,建立了预测纳米悬浮液有效导热系数的模型,通过对纳米CuO-去离子水溶液的验证,发现该模型比几种经典模型具有更高的精度,因此具有一定的参考价值。     

Effect of electrostatics techniques on the estimation of thermal conductivity via equilibrium molecular dynamics simulation: application to methane hydrate

Equilibrium molecular dynamics (MD) simulations for three system sizes of fully occupied methane hydrate have been performed at around 265 K to estimate the thermal conductivity using the Ewald, Lekner, reaction field, shifted-force and undamped Fennell–Gezelter methods. The TIP4P water model was used in conjunction with a fully atomistic methane potential with which it had been parameterized from quantum simulation. The thermal conductivity was evaluated by integration of the heat flux autocorrelation function (ACF) derived from the Green–Kubo formalism; this approach vas validated by estimation of the average phonon mean free path. The thermal conductivities predicted by non-periodic techniques were in reasonable agreement with the experimental results of 0.62 and 0.68 W/m K, although it was found that the estimates by the non-periodic techniques were up to 25% larger than those of Lekner and Ewald estimates, particularly for larger systems. The results for the Lekner method exhibited the least variation with respect to system size. A decomposition of the heat flux vector into its respective contributions revealed the importance of electrostatic interactions, and how different electrostatic treatments affect the contribution to the thermal conductivity.     

Improvement in thermal conductivity of paraffin by adding high aspect-ratio carbon-based nano-fillers

Molecular dynamics simulations have been utilized to study thermal conductivity of liquid and solid mixtures of paraffin and carbon-based high aspect-ratio nano-additives, i.e. carbon nanotubes and graphene. In agreement with existing experimental data, we observe high enhancement in thermal conductivity through adding these graphitic nano-additives into paraffin, particularly in the solid phase. We demonstrate that this significant improvement is mainly achieved by the enhancement in thermal conductivity of the matrix itself. This is caused by carbon nanotubes and graphene promoting ordering of the matrix molecules which consequently leads to improvement of its thermal conductivity.     

Molecular dynamics study for the thermal conductivity of diatomic liquid

The thermal conductivity of diatomic liquids was analyzed using a nonequilibrium molecular dynamics (NEMD) method. Five liquids, namely, O₂, CO, CS₂, Cl₂ and Br₂, were assumed. The two-center Lennard-Jones (2CLJ) model was used to express the intermolecular potential acting on liquid molecules. First, the equation of state of each liquid was obtained using MD simulation, and the critical temperature, density and pressure of each liquid were determined. Heat conduction of each liquid at various liquid states [metastable (ρ=1.9ρ_cr), saturated (ρ=2.1ρ_cr), and stable (ρ=2.3ρ_cr)] at T=0.7T_cr was simulated and the thermal conductivity was estimated. These values were compared with experimental results and it was confirmed that the simulated results were consistent with the experimental data within 10%. Obtained thermal conductivities at saturated state were reduced by the critical temperature, density and mass of molecules and these values were compared with each other. It was found that the reduced thermal conductivity increased with the increase in the molecular elongation. Detailed analysis of the molecular contribution to the thermal conductivity revealed that the contribution of the heat flux caused by energy transport and by translational energy transfer to the thermal conductivity is independent of the molecular elongation while the contribution of the heat flux caused by rotational energy transfer to the thermal conductivity increases with the increase in the molecular elongation. Moreover, by comparing the reduced thermal conductivity at various states, it was found that the increase of thermal conductivity with the increase in the density, or pressure, was caused by the increase of the contribution of energy transfer due to molecular interaction.     

Lattice thermal conductivity of δ-graphyne — A molecular dynamics study

Thermal conductivity of δ-graphyne was investigated using reverse non-equilibrium molecular dynamics simulations. The dependence of the thermal conductivities with the temperature, acetylenic linkages, and external strain were explained by the phonon density of states. Our simulations revealed that as the temperature increased, the thermal conductivity of graphene first increased and then decreased, whereas that of δ-graphyne monotonically decreased. Owing to the presence of the acetylenic linkages, a significant reduction was found in the thermal conductivity of δ-graphyne, which resulted in a phonon vibration mismatch or weakened coupling. Moreover, the temperature profile changed from mono linear to the ladder the number of acetylenic linkages increased. These results play a guidance role in the design and application of thermoelectrics devices using 2D carbon materials.