SPH方法对不可压自由表面流动的数值模拟.docVIP

SPH方法对不可压自由表面流动的数值模拟* 强洪夫，高巍然（西安高科技研究所201室，陕西西安 710025）ncompressible Flow with unsteady Free-Surface Qiang Hong Fu, Gao Wei Ran* (Xi’an Hi-Tech Institute No.201 Faculty, Shaanxi, Xi’an, PRC, 710025) Abstract：This paper describes the SPH (Smoothed Particle Hydrodynamics) method for simulation of incompressible flow with emphasis on the treatments of incompressible condition and rigid boundary condition, and the coefficients of those conditions are discussed. The 2D Dam-break problem was simulated with weak compressible condition and combined boundary condition of boundary force method and ghost particle method to illustrate the performance of these algorithms. Keyword：Smoothed Particle Hydrodynamics; Free-surface flow; Incompressible; Rigid boundary 引言 自由表面流动问题的研究在很多工业和环境工程中有着重要的意义，然而由于流体自由界面的存在给传统基于网格的N-S方程求解算法的应用带来很大障碍，特别是对流体自由表面发生复杂变形、破碎或者融合等流动问题的求解。 近年来，作为一种纯Lagrangian型的无网格算法——SPH算法[1]得到了很快发展，并越来越显示出它在解决复杂自由表面流动问题时的优势。该算法的特点是在计算中不需要生成任何辅助网格，不存在Lagrangian网格算法中网格缠绕的技术瓶颈，理论上可以处理任意的变形问题；由于算法的Lagrangian特性，SPH离散粒子自然地追踪流体物质的运动及边界的变化，不需要像Eluerian网格算法中采取复杂的界面追踪技术，不存在因数值耗散造成界面追踪精度的降低。SPH算法的特点使得它在处理自由表面流动问题有着得天独厚的优势，并且SPH算法很容易向高维问题扩展，特别适合于求解三维算例。 从目前的情况来看，SPH算法在自由表面流动问题中的应用主要遇到不可压缩条件和固壁边界条件难以处理的两个困难。本文将阐述SPH求解自由表面流动问题算法，特别是不可压缩条件和固壁边界条件的处理方法，利用作者所编写的SPH求解程序，对溃坝问题进行数值模拟，讨论所选用方法的优缺点。文章最后展望了SPH算法在这一领域的研究发展方向。 SPH基本方程组 SPH的核心思想是核函数插值对于任意连续场函数，， (1) 式(1)中为插值核函数,通常选取具有紧支性的偶函数具有，并在时趋向一个强尖峰函数。是插值核宽度的一种度量，称为光滑长度，表示不显著为零时的取值范围，通常设定时。离散的求和插值近似处场函数 (2) 式(2)中下标为核函数紧支域内所有其它粒子的编号。粒子的质量为、坐标为、密度为。如果选取的核函数是可微的，导数可以 (3) 采用核函数近似，可将连续的NS方程组离散为一组常微分方程组，即SPH基本方程组，再对这组方程采用相应的常微分方程组求解方法来推进时间进程的求解。下面是一组笛卡尔坐标系下无粘流体动力学的SPH离散方程组： (4) (5) (6) (7) 在式(4)~(7)中，上标表示矢量沿坐标轴方向上的分量，、、、和分别为粒子的密度、速度、内能和坐标，表示粒子受到的体积力，是人为粘性项，，，。 SPH算法在模拟含激波间断问题时，需要采用人为粘性项以获得稳定解，其最常用的是J. J. Monaghan提出的形式[1] ，如式(8) (8) 式(10)中，，，，为粒子的当地声速，。，，为可调参数，在含激波间断问题中选取，，。在无激波间断问题的数值模拟中，为了使计算更加稳定，也采用较小参数的人为粘性项。 不可压缩条件的处理 传统SPH方法是基于密度的算法，和其它基于密度算法一样，在处理不可压问题时，难以保持密度场绝对不变。通常SPH算法在处理不可压条件时，将不可压流体近似处理为弱可压缩流体， 如选取流体的物态方程为[3] (9