计算机数值仿真逐渐成为解决现代工程和科学问题的一条重要途径。数值仿真能为理论提供测试和检验,有助于对复杂的物理问题加深认识,甚至还能帮助解释和发现新现象。基于网格的数值方法虽然已经有广泛的应用,但是在很多方面仍存在不足之处,比如在计算流体动力学中的大变形、运动物质交界面、自由表面等问题时,由于网格产生畸变导致计算误差过大或无法进行,从而使其在许多问题的应用上受到限制。近年来,无网格法倍受关注,这种方法在许多应用中都优于传统的基于网格的有限元法、有限差分法以及有限体积法等数值方法。本文依托山东大学虚拟工程研究中心和新加坡国立大学ACES实验室,系统地研究了新一代无网格方法—光滑粒子流体动力学方法(SPH)在应用及扩展过程中的相关关键技术,实现了SPH方法的两相耦合应用和三维应用,在此基础上对三维管道气力输送过程进行了数值仿真,扩展了SPH方法的工程应用领域,为SPH方法真正辅助试验创造了条件,具有重要的理论和应用价值。论文的主要内容如下：系统总结了基于网格的数值方法,指出其在很多方面存在不足之处,由于网格畸变而使其在许多问题的应用上受到限制,在此基础上,详细分析了新一代无网格法—光滑粒子动力学方法的基本思想和求解问题的过程,由于问题域粒子之间没有任何连接,运算中使用的粒子取决于当前局部分布的粒子,得出了SPH方法的无网格属性和自适应性。系统分析了拉格朗日型的Navier-Stokes方程,应用SPH粒子近似法推导出了Navier-Stokes的密度方程、动量方程和能量方程的SPH表达式。研究了流体动力学SPH数值计算中的关键技术,指出了为了使算法适合模拟不同的流体特性问题,必须对算法进行特别地处理。在SPH方法中有两种方式对密度进行展开：一种是密度求和法,另一种是连续密度法。虽然密度求和法所需要的计算量大,但因为密度求和法体现了SPH近似法的本质,所以使用得较为广泛。连续性密度法主要用于仿真具有强间断的问题。在SPH方法中,用核函数来近似狄拉克δ函数,分析了在实际应用中常用的两种核函数：高斯核函数和三次样条核函数。在SPH的应用中,边界条件的处理既是该方法的优点,也是目前的薄弱环节。研究了处理固定边界条件的两种类型的虚粒子,第一种类型的虚粒子设置在固定边界上,第二种类型的虚粒子分布在边界的外部,通常在边界条件不断变化的场合下使用。第二种类型的虚粒子按以下的方式构造,即给定一个实粒子i,则在边界外与实粒子对称处分布一个虚粒子,这些虚粒子具有与相对应实粒子相同的压力和密度,但速度方向相反。为了避免数值震荡,提高算法的稳定性,防止粒子间相互接近时的非物理穿透,在SPH方法的动量方程中引入了人工粘度来进行修正。本文使用跳蛙法进行时间积分,跳蛙法的优点是计算时所需要的存储量低,而且在每一次计算中只需要进行一次优化估值。给出了SPH算法的程序结构,SPH的基本方法和SPH中其它数值方面的伴随算法使得SPH程序具有一些特殊性质。这些特殊性质都包含在时间积分过程的主循环中。对腔内剪切流动问题、冲击管问题相关算例进行了研究,测试了SPH方法在不同的流体动力学问题中的实用性,对于以上问题,SPH方法均可得到令人满意的结果。在分析SPH方法单相流基本原理和离散思想的基础上,提出了SPH的两相耦合流动仿真方法。由于SPH方法不需要背景网格,是一种纯无网格方法,在对问题进行粒子离散化的过程中,布置的粒子本身具有物理属性,因此,只要能够正确处理具有相同或者不同属性的粒子之间的相互作用,通过整体粒子的运动分布就能够辨析出不同物质之间的交界面,这样也就能够描述具有不同属性的两相流体的运动情况。在SPH两相耦合流动仿真中,可以通过将固体相或液体相离散为与气体相一样的SPH粒子,但其上携带的是固体相或液体相的物理属性来实现。研究了保证SPH两相耦合流动数值仿真过程正常进行需要做的技术改进。将密度正则化引入SPH两相耦合流动数值仿真方程,解决了两相流体交界面附近由于密度差异过大导致的边界效应,提高了密度不同或不连续的两相流体交界面处的精度。研究了SPH两相耦合流动中的人工状态方程,通过在两相耦合流动中密度小的物质的状态方程中添加气体内聚力项,解决了密度较小的粒子的逃逸问题。研究了SPH两相耦合流动中的速度修正,引入XSPH中的速度修正方法,防止了粒子间的相互穿透,使得到的两相耦合流动交界面更加清晰和光滑。综合使用SPH两相耦合流动方法,实现了气泡上浮和溃坝两相流动的数值仿真,得到的结果验证了本文提出的SPH两相耦合流动方法的正确性和可行性。实现了SPH方法在三维流体流动中的技术应用。分析了全配对搜索法和链表搜索法,将Bucket搜索优化算法和树形搜索优化算法引入到SPH方法中,并给出了程序实现,有效地提高了光滑粒子法的计算效率。深入研究了粒子对的相互作用,由于粒子间的相互作用是基于点对点方式的,因此可在SPH仿真分析中应用成对相互作用法来提高计算效率和节省存储空间。成对相互作用法是通过应用最近相邻粒子搜索法来实现的,并为之后进行的SPH求和法存储必须的数据。研究了提高计算结果精度的几种方法：在对核近似连续性理论分析的基础上,提出了初始粒子的配置方法,须尽量将粒子均匀分布,并且使所有粒子质量相同或者质量呈连续变化；采用一种简单高效的基于密度变换的可变光滑长度技术；对光滑长度进行优化和松弛更新,给出了优化和松弛过程的步骤和参数,有效解决了现有SPH方法在处理三维流体流动过程中计算误差大的问题。应用SPH方法实现了三维管道气力输送过程的数值仿真,分析了影响管道气力输送过程的主要因素,推导了粒子速度和输送气体速度的关系,提出了拟合方法。研究了气力输送过程中管道边界的实现方法,给出了使用应变率及应力的迭代近似法求解粘性力的方法,并分析了管道中输送物料的分布及压力情况。提出了管道边界层流效应的仿真方法：通过设置边界虚粒子和实粒子的粘度参数,而在管道中部的粒子并不设置粘度参数,按无粘处理,能够模拟边界层流效应。输送过程中气流碰到刚性壁面后将产生壁面热量,研究了壁面热量能量转换,给出了考虑人工热量条件下的SPH能量方程。合理简化影响要素,综合运用前面的相关技术,实现了管道输送问题的SPH程序原型。研究了管道气力输送SPH数值仿真方法的VC++实现技术,由于Fortran具有接近数学公式的自然描述,并且计算精度高,在计算机里具有很高的执行效率,是目前流行较广的一种面向过程的适用于科学计算的高级语言,因此整个SPH程序实现选用Fortran来编写,但是它在实现人机对话、交互处理方面不很理想,界面也不够美观,图形处理功能也比较欠缺。VC++是目前个人计算机上深受欢迎的面向对象的程序设计环境。VC++具有强大的功能,但在科学计算方面实现起来却比较复杂。因此,用VC++和Fortran进行混合编程可以编制出具有友好界面和良好计算功能的程序。本文通过综合运用SPH基本原理,SPH两相耦合流动方法,SPH三维仿真技术,实现了气力输送过程的SPH数值仿真,在保持计算效率的基础上,通过VC++调用Fortran生成的动态链接库(.dll)文件,实现了SPH程序的交互功能。综上所述,本文提出的SPH两相耦合流动仿真方法,将SPH方法延伸到了两相耦合流动中,完善了现有的光滑粒子流体动力学应用理论。本文实现了SPH方法在三维仿真流动中的应用,将树形搜索优化算法引入到SPH方法中,有效地提高了光滑粒子法的计算效率,特别是在粒子数量较多的时候,计算效率优势更加明显。本文给出了提高计算结果精度的几种方法。首次将SPH方法应用在三维管道气力输送仿真中,扩展了SPH方法的工程应用领域,给出了考虑粘性和壁面热量转换的计算方法,分析了管道中输送物料的分布及压力情况,提出了管道边界层流效应的仿真方法。在保持计算效率的基础上,通过VC++调用Fortran生成的动态链接库文件,实现了SPH程序的交互功能,为SPH方法真正辅助试验创造了条件。这些工作具有重要的理论和工程应用价值。本课题得到了国家建设高水平大学公派研究生项目(留金出[2007]3020)的支持。