重构核粒子法对光滑粒子法的改进效果

计算物理 ›› 2009, Vol. 26 ›› Issue (4): 553-558.

重构核粒子法对光滑粒子法的改进效果

殷建伟¹, 马智博²

1. 中国工程物理研究院北京研究生部, 北京 100088;
2. 北京应用物理与计算数学研究所, 北京 100088

收稿日期:2008-03-25 修回日期:2008-05-23 出版日期:2009-07-25 发布日期:2009-07-25
作者简介:殷建伟(1983-),江苏泰兴,硕士生,从事计算流体力学方面的研究,北京市2101信箱研究生部100088.

Reproducing Kernel Particle Method in Smoothed Particle Hydrodynamics

YIN Jianwei¹, MA Zhibo²

1. Graduate School, China Academy of Engineering Physics, Beijing 100088, China;
2. Institute of Applied Physics and Computational Mathematics, Beijing 100088, China

Received:2008-03-25 Revised:2008-05-23 Online:2009-07-25 Published:2009-07-25

摘要/Abstract

摘要： 光滑粒子法通过核函数进行近似估计,在计算域边界附近,核估计的精度明显下降.重构核粒子法通过校正函数对核函数进行重新构造,提高核估计方法在边界点和内点上对函数的估计精度以及计算的稳定性.研究发现,虽然校正函数的构造立足于对函数的精确估计,但这个优势同样能在对导数的估计中继续保持.通过理论研究及数学、物理模型的模拟,展示重构核粒子法的改进效果,揭示其能够提高精度的原因.

关键词: 重构核粒子法, 光滑粒子法, 无网格方法, 高速碰撞

Abstract: Kernel functions approximate physical quantities in smoothed particle hydrodynamics (SPH) method.However,approximation accuracy decreases near boundaries.By reconstruction of a kernel with correcting function,reproducing kernel particle method (RKPM) improves approximation accuracy and computational stability at internal or boundary points.Though correction is proposed to approximate functions accurately,it also benefits in approximation of derivatives.

Key words: reproducing kernel particle method, smoothed particle hydrodynamics, meshless method, high-velocity impact

中图分类号:

O347
O389

殷建伟, 马智博. 重构核粒子法对光滑粒子法的改进效果[J]. 计算物理, 2009, 26(4): 553-558.

YIN Jianwei, MA Zhibo. Reproducing Kernel Particle Method in Smoothed Particle Hydrodynamics[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2009, 26(4): 553-558.

[1]	马智博. 物理质团法:一个普适的数值方法体系[J]. 计算物理, 2017, 34(3): 261-272.
[2]	郑兴, 马庆位, 段文洋. K2_SPH方法及二维破碎波的模拟[J]. 计算物理, 2012, 29(3): 317-325.
[3]	王晓东, 欧阳洁, 王玉龙, 蒋涛. 积分节点迎风偏移的节点积分无单元Galerkin方法[J]. 计算物理, 2012, 29(2): 183-190.
[4]	郑兴, 段文洋. K2_SPH方法及其对二维非线性水波的模拟[J]. 计算物理, 2011, 28(5): 659-666.
[5]	杨秀峰, 彭世镠, 刘谋斌. 物体入水的光滑粒子法模拟[J]. 计算物理, 2011, 28(4): 523-528.
[6]	孙顺凯. 二维可压缩流体力学Lagrange有限点方法[J]. 计算物理, 2011, 28(2): 159-166.
[7]	卞梁, 王肖钧, 章杰, 赵凯. 高速碰撞数值计算中的SPH分区算法[J]. 计算物理, 2011, 28(2): 207-212.
[8]	孙顺凯, 沈隆钧, 沈智军. 耦合界面计算格式的一维流体力学Lagrange有限点方法[J]. 计算物理, 2010, 27(3): 317-325.
[9]	杨秀峰, 彭世镠. 溃坝流的光滑粒子法模拟[J]. 计算物理, 2010, 27(2): 173-180.
[10]	徐金中, 汤文辉. 高速碰撞SPH方法模拟中的初始光滑长度和粒子间距[J]. 计算物理, 2009, 26(4): 548-552.
[11]	温万治, 董师舜, 杭旭登, 胡晓燕. 二维MLSPH无网格方法[J]. 计算物理, 2008, 25(2): 157-162.
[12]	李寿佛, 张瑗, 刘玉珍. 热传导方程的一类无网格方法[J]. 计算物理, 2007, 24(5): 573-580.
[13]	曾亿山, 曾清红. 一种改进的无单元伽辽金方法[J]. 计算物理, 2007, 24(1): 35-41.
[14]	曾清红, 卢德唐. 耦合径向基函数与多项式基函数的无网格方法[J]. 计算物理, 2005, 22(1): 43-50.
[15]	张刚明, 王肖钧, 王元博, 王吉, 王峰. 高速碰撞数值计算中的光滑粒子法[J]. 计算物理, 2003, 20(5): 447-454.