2.1.問題描述

本次研究擬采用LJ體系模擬二維Couette flow，Couette flow（庫愛特流）指的是粘性流體在相對運動著的兩平行平板之間的層流流動。這個流動是由作用在流體上的粘性力和與平板平行的外部壓力推動的。本次研究通過固定底端，移動頂端來制造Couette flow。

2.2.模型描述

具體模型如圖2.1所示。本次模擬采用LJ約化單位，初始晶體模型為六方最密堆積結構，晶格參數為0.7，沿x（100）方向為20倍晶格長度，y方向（010）為20倍晶格長度。此次模型為2維模型，x為流動方向，因此設置為周期性邊界。y方向采用收縮邊界，以模擬平板移動。采用OVTIO進行模型可視化處理。在模擬流動前先設置流動區域和平板區域。具體方式為采用velocity和fix setforce命令固定底端和頂端1倍晶格長度的區域作為平板。對中間的流體區域采用速度標定法進行控溫。在進行流動模擬時，為頂端的固定區域設置沿x方向的初速度為5.0，其他方向速度為0。底端繼續保持固定。流動模擬一共運行100000步。

圖2.1：模型示意圖

2.3結果整理與分析

圖2.2展示了在初始1000步，50000步和100000步時流體原子沿著y方向的x方向速度（vx）的分布情況。從圖中可以看出有平板與流體之間粘性力帶起的流體運動存在著明顯的滯后現象。這樣的滯后體現在空間和時間尺度上。在空間尺度表現為從固定端到移動端存在著明顯的速度梯度。從時間尺度上表現為流體的速度隨著時間逐漸增加。同時還利用OVITO分析了沿y軸不同位置原子的移動軌跡，如圖2.3所示。這里也可以清楚的看的靠近頂端移動平板的原子在相同時間內有著更長的移動距離。這樣的可視化可為Couette flow的直觀分析帶來更大的便利。

圖2.2：在不同時刻流體原子沿y軸的vx分布情況

圖2.3：在不同時刻流體原子的軌跡線

最后, 歡迎通過公眾號"320科技工作室"與我們聯系