论文部分内容阅读
大规模科学与工程计算已经成为科学研究中与试验和理论研究方法同等重要的第三种方法。单核处理器受硬件发展的的限制,计算能力很难再有大的提高,其现有计算能力远远不能满足大规模科学与工程计算的要求,高性能计算成为满足这种要求的唯一有效途径。高性能计算除了需要有良好的并行计算机系统结构做支撑,更需要有适合并行计算的算法才能很好解决大规模科学与工程问题的计算。复杂流体运动的数值模拟一直是大规模科学与工程计算最重要且具有挑战性的研究领域之一。传统求解Navier-Stokes方程的方法需要求解大规模的代数方程组,Navier-Stokes方程显式解法并行性较好,但是数值稳定性差,收敛速度慢;隐式解法数值稳定性好,收敛速度快,但是其并行可扩展性较差,由于这些因素,极大限制了该方法求解大规模问题的能力。建立在Boltzmann方程介观模型的Lattice Boltzmann Methods(LBM)是近年来发展起来的一种模拟流体流动新的计算方法。LBM方法具有计算简单,能够处理复杂边界问题,天然并行等很多优点,但是LBM方法本身的局限使得求解高速可压缩问题非常困难,大大的限制了该方法的应用范围,本文探讨新的LBM方法理论,开展对可压缩高速流动研究,并应用新方法进行了大量数值实验,验证了可压缩LBM方法的正确性。实现了可压缩高速LBM方法并行算法,并对其并行性能进行了测试分析。除了LBM方法,Gas-Kinetic方法也是建立在Boltzmann方程模型上模拟流体流动新的数值方法,该方法借鉴有限体积方法的思想,通过建立控制体积剖分计算区域,应用Boltzmann-BGK方程模型求解控制体积界面上宏观物理量通量,该方法可以求解高速可压缩流体问题。本文提出了将非结构网格应用到Gas-Kinetic方法求解复杂几何形状、复杂边界条件的可压缩高速流体问题,数值实验验证了该方法的正确性,并对该方法并行计算的可扩展性进行了分析。最后本文对这两种基于Boltzmann-BGK模型的方法进行了讨论,对比了两种方法的优缺点,并提出融合两种方法的Lattice Boltzmann Gas-Kinetic Methods(LBGKM)方法的设想。本文的创新性工作体现在以下几个方面:1.探讨新LBM理论用于可压缩高速流动问题研究,提出了适用于LBM方法的TVD(Total Variation Diminishing)格式模拟激波现象,实现该算法并得到比较好的计算结果;2.进行了可压缩LBM方法算法并行化研究,对并行算法的性能影响因素进行了深入分析,提出了改进并行性能的方法;3.通过大量的数值实验对本文的新LBM理论进行了验证,并实现了三维可压缩流动问题的并行计算,得到了三维流动的结果;4.提出对Gas-Kinetic方法采用非结构网格处理复杂边界问题,对结构和非结构网格的Gas-Kinetic方法进行了深入研究,提出了基于Gas-Kinetic方法的并行算法,并对算法进行了程序实现,初步得到了一些三维高速流动的计算结果;5.针对这两种方法编制了LBM3D和GAS3D两个规模较大的三维高速流体并行计算程序。