1、計算共存相之間的界面張力背景 界面張力是評估兩種流體能否互相混合的重要物理學量。DPD 模塊提供了計算兩種流體間界面張力的可能性。DPD 所使用的算法會計算體系單元內沿著X 軸方向的界面張力值。為了能夠達到這一目的，我們需要在X 軸方向延長體系的長度。通過這種改變，會使得X 界面會有比較小的表面積，從而導致較低的自由能。該教程包括： 設置非立方的二元混合體系 執行DPD 分析計算結果并獲得界面張力11設置非立方的二元混合體系 在本教程中，我們將會研究一個非常簡單的油水體系。水分子將會使用一個單獨的珠子來表述，油分子則表述為較短的珠子鏈。油珠子和水珠子的斥力將會被設置為較大值，以保證能夠快速的發

2、生相分離。 從模塊工具欄上選擇DPD 工具 ，然后選擇Calculation。選擇Species 標簽。 對于這一計算，你需要兩種不同的珠子類型。你可以將其標記為O 和W，用于反映油和水。在Bead types 部分，將A 改為O，將B 改為W。 你可以自己定義油和水珠子，但是現在必須定義用什么樣的珠子構造你的分子。你將會使用單個W 珠子來定義水分子，并使用八個O 珠子來定義油分子。 在Mesoscale molecule type 部分將Diblock 的名稱改為Oil。在Topology 區域輸入O 8。在下一行，輸入Water 作為水分子名稱，并在Topology 輸入W 1 作為水分子

3、。 這樣定義了要計算的物種。2 現在你需要指定這些珠子之間是如何作用的。選擇Interactions 標簽。在本例中相互作用的默認值就比較合適，因此在這里不需要作任何改變。 選擇System 標簽。在當前體系中，并沒有包括水。你必須改變兩種組分的相對數量，以使油和水的數量相等。將Relative amounts fields 的兩個數值都改為0.5。 如上面所說，需要一個非立方的模擬單元來確定界面發生在X 方向。因此，需要增加X方向的長度，并且縮小Y、Z 方向的長度。在System extent 部分的第一項上單擊，將數值從10.0 改為15。將其他兩部分改為6.0。3 3000 步的計算對于

4、產生所需的數據已經是足夠多的了。現在我們需要確認每過250 步就將體系的密度值和介觀分子子集就被保存一次以便于后期分析。 選擇Setup 標簽。將Number of Steps 設置為3000，并將Frame every 設置為250。確定Coordinates every 設置為1，并且Restart 設置為3000。3000/250=12是指經過12步，水油分離。42. 執行DPD 模擬 現在需要指定輸出輸出文件的名稱。 選擇Job Control 標簽。去掉Automatic 前面的對勾，并在Job Description 里輸入Oil_Water。計算所有產生的文檔名稱都將會包括oil

5、_water 的前綴。 按下Run 按鈕。等待作業結束，并進入到下一部份。53分析運算結果，并獲得界面張力 如下圖，計算結果為oil_water.mtd 文檔。 按下Animation 工具欄上的按鈕。從下拉菜單中選Options將playback speed 設為Slow將playback style設為Loop(提示：要顯示Animation 工具欄，請首先從View 菜單中選擇Toolbars，并選中Animation。)6現在可以通過動畫播放計算的結果。 * 確定oil-water.mtd被激活 * 單擊Animation 工具欄的 * 在動畫演示的過程中，油珠子和水珠子將發生分離。從Animation 中的 Options上可以看到每一步時水、油的混合情況。)7 表面張力的單位為質量.時間-2。該數據以DPD 單位給出。質量的單位就是珠子的質量，而時間單位則以