光滑质点流体动力学_SPH_方法_综述_.pdfVIP

第 13 卷第4 期 计 算 物 理 Vol. 13,No. 4 1996 年 12 月 CHINESE JOURNAL OF COM PUTATIONAL PHYSICS Dec. , 1996 ( SPH) ( ) 张锁春 ( 中国科学院应用数学研究所, 北京 100080) 综合介绍了一种新 纯( Lagrangian 方法) ) ) 光滑质点流体动力学( SPH) 方法, 由于它计算 空间导数时不需要使用网格, 从而避免了高维拉氏差分网格法中 网格缠结和扭曲 麻烦, 它对非 对称 和内含真空区域 三维问题特别有效。重点介绍了该方法 理论基础, 流体动力学方程组 推导, 人为粘性和热传导, 自引力、汇和磁场, 光滑核 选取, 以及SPH 执行过程等有关问题。 SPH 光滑核 积分插值 O35 11 2 1 大多数流体动力学问题, 由于其复杂性都要求数值计算, 为此发展了多种数值方法。其中 有一种名叫/ 光滑质点流体动力学0 ( SPH-Smoothed Particle H drod namics) 方法, 它对天体物 理问题尤为合适, 是近廿多年来发展起来 一种新 纯Lagrange 方法。描述 SPH 技巧 第 一篇论文是 Luc 于 1977 年提出 , 经过 Gingold 和 Moanaghan ( GM77, GM 82, M82, GM83,M G83, ML85) 等人 工作, 奠定了扎实 基础, 再经过 Benz( BT90, HB9 1, HBC92, FBD95) 等人 改进和完善, 现已发展成为比较成熟 计算高维( 尤其是三维) 天体物理问题 有效方法了。本文 目 就是对SPH 方法作综述性介绍, 以促进该方法在我国国内得到应 用和发展。 SPH 是一种质点方法, 有点类似质点网格法( PIC-Particle in Cell ) ( 见Harlow 1957, 1974, 1988) , 但根本 不同点在于SPH 方法中计算空间导数时不需要使用任何网格, 而被插值公式 中 解析微分式子所替代, 从而避免了高维拉氏差分网格法中 网格缠结( tangling) 和扭曲 ( distortion) 等最令人头痛 问题。这种方法另一个突出 优点是表现在对缺乏对称性( lack s mmetr ) 和内含真空( large voids) 三维系统 计算特别有效。因为利用网格建立 所有传 统 差分方法, 随着网格数目 增加, 到了三维情形变得难以招架, 真空区域 网格剖分造成 内存贮量 大量浪费又实在令人可惜, 而 SPH 恰能克服上述 缺点。由于受到来自旋转和磁 场 非球对称力 影响, 天体物理领域内要研究 大多数问题都是非对称 三维问题, 而 SPH 方法正是适应这种需要而发展起来 ( M 85, M88, M92, B88, B90, B91) , 但也在其它领域 内得到广泛 应用, 例如高速碰撞 材料动载响应( 见Libersk 等人( 1993) ) 数值模拟。 2 SPH 方法 核心实为是一种插值( interpolation) 。在SPH 中任一宏观变量( 如密度、压力 1995 年7 月6 日收到原稿, 1996 年8 月9 日收到修改稿。 386 计 算 物 理 第 13 卷 温度等) A( r ) 能方便地借助于一组无序点(disordered points) 上 值表示成积分插值( integral interpolant) 计算得到, 其形式为 3A (r ) 4= A (r )W(r - rc, h)drc, (211) Q