论文部分内容阅读
在超级润滑性、聚合电解质刷、密闭流体中的相变等实验中,通过振荡流对流体粘弹性的研究受到了人们极大的关注。但由于不同研究组的实验条件不同,得出的结论有很大差异。因此,人们开始尝试通过数值模拟的方法在计算机中进行分子动力学模拟,以避免实验中环境所带来的不利影响。然而,在计算机模拟过程中,如果想要得到与真实情况非常近似的模拟结果,就需要对数量很大的分子进行数值计算,其计算过程的时间成本很高,往往会持续几天、几周、甚至几个月,很多问题由于耗时过长而不得不对问题进行某种程度的简化,而这些简化是以牺牲结果准确性作为代价的。本文以近年来高速发展的GPU计算方法为基础,基于CUDA和MPI的混合编程方法,提出了以GPU集群为主要计算环境,进行分子动力学模拟的实验方法,相对于CPU获得了5倍以上的加速比,这表明该方法在大规模分子动力学计算机模拟实验中有一定的应用价值。在此基础上,本文研究了基于GPU的纳米级流体振荡流的数值模拟计算。在实际模拟实验的过程中,本文采用基于域分解的GPU计算方法,以及基于GPU的快速傅立叶变换计算,有效地减少了整个模拟过程所需的时间。本文采用杂化连续微粒的限制性动力学方法来引起振荡流动,首次研究了当ω=1、0.5、0.1三种情况时纳米流体的各种性质。结果表明当振荡振幅恒定不变时,速度分布对频率的响应随着振荡频率迅速减弱。尤其是ω=1时,在接近壁面区域内的流体没有明显的振荡,然而当ω=0.1时,几乎在整个区域内都出现振荡。主要原因是远离C→P区域的微粒不能快速地对高振荡频率做出响应。同时,实验结果表明振荡流的温度分布是不对称的。另外,剪切应力的快速傅立叶变换表明在给定的条件下流体存在线性粘弹性质。这些结果在纳米流体的研究中具有一定的意义。