计算共存相之间的界面张力

会计学1计算共存相之间的界面张力计算共存相之间的界面张力1设置非立方的二元混合体系 在本教程中，我们将会研究一个非常简单的油水体系。水分子将会使用一个单独的珠子来表述，油分子则表述为较短的珠子链。油珠子和水珠子的斥力将会被设置为较大值，以保证能够快速的发生相分离。从模块工具栏上选择DPD 工具 ，然后选择Calculation。选择Species 标签。对于这一计算，你需要两种不同的珠子类型。你可以将其标记为O 和W，用于反映油和水。在Bead typesBead types 部分，将A A 改为O O，将B B 改为W W。你可以自己定义油和水珠子，但是现在必须定义用什么样的珠子构造你的分子。你将会使用单个W W 珠子来定义水分子，并使用八个O O 珠子来定义油分子。在Mesoscale molecule typeMesoscale molecule type 部分将DiblockDiblock 的名称改为OilOil。在TopologyTopology 区域输入O 8O 8。在下一行，输入WaterWater 作为水分子名称，并在TopologyTopology 输入W 1W 1 作为水分子。这样定义了要计算的物种。第1页/共9页 现在你需要指定这些珠子之间是如何作用的。选择Interactions Interactions 标签。在本例中相互作用的默认值就比较合适，因此在这里不需要作任何改变。选择SystemSystem 标签。在当前体系中，并没有包括水。你必须改变两种组分的相对数量，以使油和水的数量相等。将Relative amounts Relative amounts fields 的两个数值都改为0.50.5。如上面所说，需要一个非立方的模拟单元来确定界面发生在X 方向。因此，需要增加X方向的长度，并且缩小Y、Z 方向的长度。在System extentSystem extent 部分的第一项上单击，将数值从10.0 改为15。将其他两部分改为6.0。第2页/共9页 3000 步的计算对于产生所需的数据已经是足够多的了。现在我们需要确认每过250 步就将体系的密度值和介观分子子集就被保存一次以便于后期分析。选择Setup 标签。将Number of Steps 设置为3000，并将Frame every 设置为250。确定Coordinates every 设置为1，并且Restart file every 设置为3000。3000/250=12 是指经过12步，水油分离。第3页/共9页2.执行DPD 模拟 现在需要指定输出输出文件的名称。选择Job Control 标签。去掉Automatic 前面的对勾，并在Job Description 里输入Oil_Water。计算所有产生的文档名称都将会包括oil_water 的前缀。按下Run 按钮。等待作业结束，并进入到下一部份。第4页/共9页3分析运算结果，并获得界面张力 如下图，计算结果为oil_water.mtd 文档。按下Animation 工具栏上的按钮。从下拉菜单中选Options将playback speed 设为Slow将playback style设为Loop(提示：要显示Animation 工具栏，请首先从View 菜单中选择Toolbars，并选中Animation。)第5页/共9页现在可以通过动画播放计算的结果。*确定oil-water.mtd被激活 *单击Animation 工具栏的 *在动画演示的过程中，油珠子和水珠子将发生分离。从Animation 中的 Options上可以看到每一步时水、油的混合情况。)第6页/共9页 双击oil_water.Dpd_out 文档，将屏幕滚动到文档的最后部分。Oil-water.Dpd_out 文件的最后部分将会报道体系的界面张力。第7页/共9页 表面张力的单位为质量.时间-2。该数据以DPD 单位给出。质量的单位就是珠子的质量，而时间单位则以珠子扩散一个珠子直径的时间为准。DPD 参数可以通过以下转换公式来转换为真实参数（带有上划线的符号使用DPD 单位）：珠子的质量必须从手册（参数化工程）中来获得。第8页/共9页