从lammps自带例子中 找到viscosity文件夹 包含有五种不同的计算粘度的方法 我使用过的是G-K方法 现将例子中的in文件进行解析: 例子中采用的是二维的,#表示注释 # settings #设置 variable x equal 20 # 二维空间的xy表示的值 variable y equal 20 # variable rho equal 0.6 #晶格长度 variable t equal 1.0 ...
Green-Kubo方法的计算热导率的原理是利用系统热流的时间自相关函数进行计算,其公式为 由于体系是各向同性的所以xyz三个方向上的热流自相关函数计算出来的热导率应该是相同的,所以可以计算出来三个方向上的热导率系数然后取平均。在LAMMPS中利用compute heat/flux来计算系统的热流。这个compute会输出6个全局矢量。该矢量...
从lammps自带例子中 找到viscosity文件夹 包含有五种不同的计算粘度的方法 我使用过的是G-K方法 现将例子中的in文件进行解析: 例子中采用的是二维的,#表示注释 # settings #设置 variable x equal 20 # 二维空间的xy表示的值 variable y equal 20 # variable rho equal 0.6 #晶格长度 variable t equal 1.0 ...
平衡分子动力学(EMD)Green-Kubo方法,适用于均质体系及原子电荷较小的体系,适当大小的初始模型便可以得到准确的结果。 10.2 MedeA® Viscosity 预测纯液体以及液体混合物在不同温度和压力下的粘度。 平衡分子动力学(EMD)采用Green-Kubo方法,适用于粘度中等的体系。 反向非平衡分子动力学方法(RNEMD),适用于粘度较大的...
在运行模拟之后,使用Green-Kubo方法计算了液体的平均平方位移(msd),并通过一些计算得到了液体的粘度。 在lammps的输入文件准备好后,我们可以通过运行lammps软件来进行粘度计算。命令行中输入以下命令: ``` lammps -in viscosity.in ``` 其中,viscosity.in是前面准备好的lammps输入文件。执行上述命令后,lammps将进行...
MedeA基本建模工具中包含了纳米团簇,纳米管等模型的创建功能,通过形状、键长、键角、悬挂键、对称性等建模参数的设置,可以轻松创建各种不同元素的团簇,包裹,管状,线状的周期/非周期结构。 Nanoparticles能够创建团簇及基础壳层结构。 Nanotube能够创建出各种纳米管结构,比如碳纳米管、硼氮纳米管。
第四种方法是使用基于 Green-Kubo(GK)公式,其将应力张量的自相关的 整体平均值与 eta 相关联。该方法在完全平衡的模拟中完成,这与第一第二个方 法完全不同,其中动量连续地流过模拟盒子。 以下是通过 GK 公式计算液体 Ar 的粘度的示例输入脚本: # Sample LAMMPS input script for viscosity of liquid Ar ...
#采用Green-Kubo方法计算沥青分子粘度 #采用real单位系统（记住：长度：A；时间：fs）#定义几个参数 var...
#采用Green-Kubo方法计算沥青分子粘度 #采用real单位系统（记住：长度：A；时间：fs）#定义几个参数 var...
4. Calculate viscosity: Once the system reaches steady state, calculate the viscosity using the Green-Kubo relation, which relates the viscosity to the autocorrelation function of the stress tensor. Advantages of Reverse Dynamics Method: The reverse dynamics method offers several advantages for viscosity...