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

 光滑质点流体动力学( SP H) 方法( 综述) 张锁春 ( 中国科学院应用数学研究所 ,北京 100080) 要 综合介绍了一种新的纯 (Lagrangian 方法 ———光滑质点流体动力学 ( SPH) 方法 ,由于它计算 空间导数时不需要使用网格 ,从而避免了高维拉氏差分网格法中的网格缠结和扭曲的麻烦 ,它对非 对称的和内含真空区域的三维问题特别有效 。重点介绍了该方法的理论基础 ,流体动力学方程组的 推导 ,人为粘性和热传导 ,自引力 、汇和磁场 ,光滑核的选取 ,以及 SPH 执行过程等有关问题 。 SP H 光滑核 积分插值 O35112 引 言 1 大多数流体动力学问题 ,由于其复杂性都要求数值计算 ,为此发展了多种数值方法 。其中 有一种名叫“光滑质点流体动力学”( SP H2Smoot hed Particle Hydro dynamics) 方法 ,它对天体物 理问题尤为合适 ,是近廿多年来发展起来的一种新的纯 L agrange 方法 。描述 SP H 技巧的第 一篇论 文 是 L ucy 于 1977 年 提 出 的 , 经 过 Gingold 和 Moanaghan ( GM77 , GM82 , M82 , GM83 ,M G83 ,ML 85) 等 人 的 工 作 , 奠 定 了 扎 实 的 基 础 , 再 经 过 Benz (B T90 , HB91 , HBC92 , FBD95) 等人的改进和完善 ,现已发展成为比较成熟的计算高维 (尤其是三维) 的天体物理问题 的有效方法了 。本文的目的就是对 SP H 方法作综述性介绍 ,以促进该方法在我国国内得到应 用和发展 。 SP H 是一种质点方法 ,有点类似质点网格法 ( P IC2Particle in Cell ) (见 Harlow 1957 ,1974 , 1988) ,但根本的不同点在于 SP H 方法中计算空间导数时不需要使用任何网格 ,而被插值公式 中的解析微分式子所替代 ,从而避免了高维拉氏差分网格法中的网格缠结 ( tangling) 和扭曲 ( disto rtio n) 等最令人头痛的问题 。这种方法另一个突出的优点是表现在对缺乏对称性 (lack symmet ry) 和内含真空 (large voids) 的三维系统的计算特别有效 。因为利用网格建立的所有传 统的差分方法 ,随着网格数目的增加 ,到了三维情形变得难以招架 ,真空区域的网格剖分造成 内存贮量的大量浪费又实在令人可惜 ,而 SP H 恰能克服上述的缺点 。由于受到来自旋转和磁 场的非球对称力的影响 , 天体物理领域内要研究的大多数问题都是非对称的三维问题 , 而 SP H 方法正是适应这种需要而发展起来的 ( M85 ,M88 ,M92 ,B88 ,B90 ,B91) ,但也在其它领域 内得到广泛的应用 ,例如高速碰撞的材料动载响应 (见 Libersky 等人 (1993) ) 的数值模拟 。 2 SP H 方法的核心实为是一种插值 (interpolatio n) 。在 SP H 中任一宏观变量 ( 如密度 、压力 1995 年 7 月 6 日收到原稿 ,1996 年 8 月 9 日收到修改稿 。 温度等) A ( r ) 能方便地借助于一组无序点 ( diso rdered point s) 上的值表示成积分插值 (integral interpolant ) 计算得到 , 其形式为 〈 A ( r ) 〉= ∫A ( r ) W ( r - r, h) d r′, ( 211) D 其中 D 为整个求解区域 。W 是插值核 (interpolating ker nel) ,它具有两个特性 : (i) ∫W ( r - r, h) d r′= 1 , ( 212) D (ii) limW ( r - r′, h) = δ( r - ( 213) r′) . h →0 这里 h 是核宽度的一种度量 , 暂时假定它取为定常数 。显然 , ( i) 是一个规一化条件 , 而 ( ii) 实 际上是考虑 W 在 r= r 处 , 它是一个强尖峰函数 ( St ro ngly peaked f unctio n) . 在| r - r| h , 它零 ; 而当 h →0 , 它是一个 (Delta) 函数 , 即有 lim〈 A ( r ) 〉= A ( r ) 。通常称〈 A ( r ) 〉 A h →0 ( r) 的一个核估计 , 它亦等价于围绕场 A ( r ) 带有一个光滑器 ( Smoot hing) 或过滤器 (fiter) 函 数