搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种新型光滑粒子动力学固壁边界施加模型

刘虎 强洪夫 陈福振 韩亚伟 范树佳

引用本文:
Citation:

一种新型光滑粒子动力学固壁边界施加模型

刘虎, 强洪夫, 陈福振, 韩亚伟, 范树佳

A new boundary treatment method in smoothed particle hydrodynamics

Liu Hu, Qiang Hong-Fu, Chen Fu-Zhen, Han Ya-Wei, Fan Shu-Jia
PDF
导出引用
  • 由于Lagrange粒子法的本质, 固壁边界条件的施加一直是光滑粒子动力学方法的难点之一. 本文从固壁边界的物理原理出发, 应用多层虚粒子表征固壁边界, 提出了一种新型固壁边界施加模型. 将虚粒子看作流体的扩展, 计算中虚粒子密度保持不变, 压力、速度等参数通过对流体粒子的插值获得, 虚粒子有条件的参与控制方程的计算, 对流体的密度/压力产生影响, 通过压力梯度隐式地表征壁面与流体之间的作用强度并对流体粒子施加沿壁面法线方向的斥力作用, 防止流体粒子对壁面的穿透. 数值算例测试结果表明, 与现有固壁边界施加方法相比, 本文方法更加符合流体与固壁边界作用的物理原理, 可以简单、有效地施加固壁边界条件, 方便地应用于具有复杂几何边界的问题, 获得稳定的流场形态、规则的粒子秩序及良好的速度、压力等参量的分布.
    As the smoothed particle hydrodynamics (SPH) is a truly Lagrangian meshfree method, the implementation of solid boundary condition has been one of the key problems that hinder SPH from applying to lots of engineering problems. In order to treat the boundary conditions efficiently, based on the boundary-fluid interaction principles, a new boundary treatment method is proposed. In this method, the solid boundary is repreflented implicitly by several layers of dummy particles along the boundary line. During the simulation, the dummy particles are treated as an extension of the fluid phase. The densities of dummy particles are kept constant, and the pressures and velocities are interpolated from the nearby fluid particles at each time step. Dummy particles can be involved in the calculation of the continuity equation conditionally and exert influences on the density/pressure field of the fluid phase. Then, for the fluid particles that approach the solid boundary, local pressure gradients are used to repreflent the dummy-fluid particle pair’s interaction strength and act as the boundary force term implicitly, which is tuned to be repulsive only and normal to the boundary. Thus, large pressure gradients mean strong boundary-fluid interaction strength, and the boundary force from the dummy particles should also be large enough to preflent the fluid particles from penetrating the solid boundary; and on the contrary, small pressure gradients mean weak boundary-fluid interaction strength and the boundary force becomes soft and little disturbs the flow field. Results of numerical tests demonstrate that, compared with the existing boundary treatment methods, the new method is in better accordance with the physical principles of the fluid-boundary interaction, and is able to treat arbitrary solid boundaries with limited modeling and computational costs. With the help of this new boundary treatment method, the stable flow field, well-ordered particle distribution, smooth velocity and pressure fields could be obtained. Theoretically, this new boundary treatment method could be directly used in three-dimensional multi-phase problems. Further tests are planned to be carried out; meanwhile, expanding the new boundary treatment method to rigid-fluid interaction problems is also a work in the future.
    • 基金项目: 国家自然科学基金(批准号: 51276192)资助的课题.
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 51276192).
    [1]

    Lucy L B 1977 Astron. J. 82 1013

    [2]

    Gingold R A, Monaghan J J 1977 Mon. Not. R. Astron. Soc. 181 375

    [3]

    Zhang A M 2008 Chin. Phys. B 22 927

    [4]

    Sun Z H, Han R J 2008 Chin. Phys. B 17 3185

    [5]

    Qiang H F, Shi C, Chen F Z, Han Y W 2013 Acta Phys. Sin. 62 214701 (in Chinese) [强洪夫, 石超, 陈福振, 韩亚伟 2013 物理学报 62 214701]

    [6]

    Monaghan J J 2005 Rep. Prog. Phys. 68 1703

    [7]

    Monaghan J J 1994 J. Comput. Phys. 110 399

    [8]

    Monaghan J J, Kajtar J B 2009 Comput. Phys. Commun. 180 1811

    [9]

    Liu M B, Shao J R, Chang J Z 2012 Sci. China Technol. Sc. 55 244

    [10]

    Han Y W, Qiang H F, Zhao J L, Gao W R 2013 Acta Phys. Sin. 62 044702 (in Chinese) [韩亚伟, 强洪夫, 赵玖玲, 高巍然 2013 物理学报 62 044702]

    [11]

    Han Y W, Qiang H F, Wang K P, Gao W R 2011 Eng. Mech. 28 245 (in Chinese) [强洪夫, 韩亚伟, 王坤鹏, 高巍然 2011 工程力学 28 245]

    [12]

    Morris J P, Fox P J, Zhu Y 1997 J. Comput. Phys. 136 214

    [13]

    Liu M B, Liu G R, Lam K Y 2002 Shock Waves 12 181

    [14]

    Colagrossi A, Landrini M 2003 J. Comput. Phys. 191 448

    [15]

    Colagrossi A, Lugni C, Brocchini M 2010 J. Hydraul. Res. 48 94

    [16]

    Marrone S, Antuono M, Colagrossi A, Colicchio G, Touzé D L, Graziani G 2011 Comput. Methods Appl. Mech. Engrg. 200 1526

    [17]

    Adami S, Hu X Y, Adams N A 2012 J. Comput. Phys. 231 7057

    [18]

    Liu M B, Chang J Z 2010 Acta Phys. Sin. 59 3654 (in Chinese) [刘谋斌, 常建忠 2010 物理学报 59 3654]

    [19]

    Liu G R, Liu M B (translate by Han X, Yang G, Qiang H F) 2005 Smoothed Particle Hydrodynamics: A Meshfree Particle Method (Changsha: Hunan University Press) pp58-67 (in Chinese) [Liu G R, Liu M B (韩旭, 杨刚, 强洪夫 译) 2005 光滑粒子流体动力学-一种无网格粒子法 (长沙:湖南大学出版社)第58-67页]

    [20]

    Bonet J, Lok T S L 1999 Comput. Method. Appl. M. 180 97

    [21]

    Marrone S, Bouscasse B, Colagrossi A, Antuono M 2012 Comput. Fluids 69 54

    [22]

    Chen J K, Beraun J E, Carney T C 1999 Int. J. Numer. Meth. Eng. 46 231

    [23]

    Schmid M, Klein F 1995 NADCA 18. International Die Casting Congress and Exposition, Indianapolis, 1995, p93

    [24]

    He Y 2012 M.S. Thesis (Guangzhou: South China University of Technology) (in Chinese) [何毅 2012 硕士学位论文 (广州: 华南理工大学)]

    [25]

    Cleary P W, Ha J 2000 Int. J. Cast Metal. Res. 12 409

    [26]

    Koshizuka S, Oka Y 1996 Nucl. Sci. Eng. 123 421

    [27]

    Lobovsky L, Botia-Vera E, Castellana F, Mas-Soler J, Souto-Iglesias A 2014 J. Fluid. Struct. 48 407

    [28]

    Monaghan J J 2012 Annu. Rev. Fluid Mech. 44 323

  • [1]

    Lucy L B 1977 Astron. J. 82 1013

    [2]

    Gingold R A, Monaghan J J 1977 Mon. Not. R. Astron. Soc. 181 375

    [3]

    Zhang A M 2008 Chin. Phys. B 22 927

    [4]

    Sun Z H, Han R J 2008 Chin. Phys. B 17 3185

    [5]

    Qiang H F, Shi C, Chen F Z, Han Y W 2013 Acta Phys. Sin. 62 214701 (in Chinese) [强洪夫, 石超, 陈福振, 韩亚伟 2013 物理学报 62 214701]

    [6]

    Monaghan J J 2005 Rep. Prog. Phys. 68 1703

    [7]

    Monaghan J J 1994 J. Comput. Phys. 110 399

    [8]

    Monaghan J J, Kajtar J B 2009 Comput. Phys. Commun. 180 1811

    [9]

    Liu M B, Shao J R, Chang J Z 2012 Sci. China Technol. Sc. 55 244

    [10]

    Han Y W, Qiang H F, Zhao J L, Gao W R 2013 Acta Phys. Sin. 62 044702 (in Chinese) [韩亚伟, 强洪夫, 赵玖玲, 高巍然 2013 物理学报 62 044702]

    [11]

    Han Y W, Qiang H F, Wang K P, Gao W R 2011 Eng. Mech. 28 245 (in Chinese) [强洪夫, 韩亚伟, 王坤鹏, 高巍然 2011 工程力学 28 245]

    [12]

    Morris J P, Fox P J, Zhu Y 1997 J. Comput. Phys. 136 214

    [13]

    Liu M B, Liu G R, Lam K Y 2002 Shock Waves 12 181

    [14]

    Colagrossi A, Landrini M 2003 J. Comput. Phys. 191 448

    [15]

    Colagrossi A, Lugni C, Brocchini M 2010 J. Hydraul. Res. 48 94

    [16]

    Marrone S, Antuono M, Colagrossi A, Colicchio G, Touzé D L, Graziani G 2011 Comput. Methods Appl. Mech. Engrg. 200 1526

    [17]

    Adami S, Hu X Y, Adams N A 2012 J. Comput. Phys. 231 7057

    [18]

    Liu M B, Chang J Z 2010 Acta Phys. Sin. 59 3654 (in Chinese) [刘谋斌, 常建忠 2010 物理学报 59 3654]

    [19]

    Liu G R, Liu M B (translate by Han X, Yang G, Qiang H F) 2005 Smoothed Particle Hydrodynamics: A Meshfree Particle Method (Changsha: Hunan University Press) pp58-67 (in Chinese) [Liu G R, Liu M B (韩旭, 杨刚, 强洪夫 译) 2005 光滑粒子流体动力学-一种无网格粒子法 (长沙:湖南大学出版社)第58-67页]

    [20]

    Bonet J, Lok T S L 1999 Comput. Method. Appl. M. 180 97

    [21]

    Marrone S, Bouscasse B, Colagrossi A, Antuono M 2012 Comput. Fluids 69 54

    [22]

    Chen J K, Beraun J E, Carney T C 1999 Int. J. Numer. Meth. Eng. 46 231

    [23]

    Schmid M, Klein F 1995 NADCA 18. International Die Casting Congress and Exposition, Indianapolis, 1995, p93

    [24]

    He Y 2012 M.S. Thesis (Guangzhou: South China University of Technology) (in Chinese) [何毅 2012 硕士学位论文 (广州: 华南理工大学)]

    [25]

    Cleary P W, Ha J 2000 Int. J. Cast Metal. Res. 12 409

    [26]

    Koshizuka S, Oka Y 1996 Nucl. Sci. Eng. 123 421

    [27]

    Lobovsky L, Botia-Vera E, Castellana F, Mas-Soler J, Souto-Iglesias A 2014 J. Fluid. Struct. 48 407

    [28]

    Monaghan J J 2012 Annu. Rev. Fluid Mech. 44 323

  • [1] 许晓阳, 周亚丽, 余鹏. eXtended Pom-Pom黏弹性流体的改进光滑粒子动力学模拟. 物理学报, 2023, 72(3): 034701. doi: 10.7498/aps.72.20221922
    [2] 邵绪强, 梅鹏, 陈文新. 基于稳定性SPH-SWE数值模型的真实感流体动画实时模拟. 物理学报, 2021, 70(23): 234701. doi: 10.7498/aps.70.20211251
    [3] 林晨森, 陈硕, 肖兰兰. 适用复杂几何壁面的耗散粒子动力学边界条件. 物理学报, 2019, 68(14): 140204. doi: 10.7498/aps.68.20190533
    [4] 蒋涛, 黄金晶, 陆林广, 任金莲. 非线性薛定谔方程的高阶分裂改进光滑粒子动力学算法. 物理学报, 2019, 68(9): 090203. doi: 10.7498/aps.68.20190169
    [5] 蒋涛, 陈振超, 任金莲, 李刚. 基于修正并行光滑粒子动力学方法三维变系数瞬态热传导问题的模拟. 物理学报, 2017, 66(13): 130201. doi: 10.7498/aps.66.130201
    [6] 马理强, 苏铁熊, 刘汉涛, 孟青. 微液滴振荡过程的光滑粒子动力学方法数值模拟. 物理学报, 2015, 64(13): 134702. doi: 10.7498/aps.64.134702
    [7] 雷娟棉, 黄灿. 一种改进的光滑粒子流体动力学前处理方法. 物理学报, 2014, 63(14): 144702. doi: 10.7498/aps.63.144702
    [8] 雷娟棉, 杨浩, 黄灿. 基于弱可压与不可压光滑粒子动力学方法的封闭方腔自然对流数值模拟及算法对比. 物理学报, 2014, 63(22): 224701. doi: 10.7498/aps.63.224701
    [9] 蒋涛, 任金莲, 徐磊, 陆林广. 非等温非牛顿黏性流体流动问题的修正光滑粒子动力学方法模拟. 物理学报, 2014, 63(21): 210203. doi: 10.7498/aps.63.210203
    [10] 邱流潮. 基于不可压缩光滑粒子动力学的黏性液滴变形过程仿真. 物理学报, 2013, 62(12): 124702. doi: 10.7498/aps.62.124702
    [11] 蒋涛, 陆林广, 陆伟刚. 等直径微液滴碰撞过程的改进光滑粒子动力学模拟. 物理学报, 2013, 62(22): 224701. doi: 10.7498/aps.62.224701
    [12] 苏铁熊, 马理强, 刘谋斌, 常建忠. 基于光滑粒子动力学方法的液滴冲击固壁面问题数值模拟. 物理学报, 2013, 62(6): 064702. doi: 10.7498/aps.62.064702
    [13] 韩亚伟, 强洪夫, 赵玖玲, 高巍然. 光滑粒子流体动力学方法固壁处理的一种新型排斥力模型. 物理学报, 2013, 62(4): 044702. doi: 10.7498/aps.62.044702
    [14] 强洪夫, 刘开, 陈福振. 液滴在气固交界面变形移动问题的光滑粒子流体动力学模拟. 物理学报, 2012, 61(20): 204701. doi: 10.7498/aps.61.204701
    [15] 杨秀峰, 刘谋斌. 光滑粒子动力学SPH方法应力不稳定性的一种改进方案. 物理学报, 2012, 61(22): 224701. doi: 10.7498/aps.61.224701
    [16] 马理强, 刘谋斌, 常建忠, 苏铁熊, 刘汉涛. 液滴冲击液膜问题的光滑粒子动力学模拟. 物理学报, 2012, 61(24): 244701. doi: 10.7498/aps.61.244701
    [17] 马理强, 常建忠, 刘汉涛, 刘谋斌. 液滴溅落问题的光滑粒子动力学模拟. 物理学报, 2012, 61(5): 054701. doi: 10.7498/aps.61.054701
    [18] 蒋涛, 欧阳洁, 栗雪娟, 张林, 任金莲. 瞬态热传导问题的一阶对称SPH方法模拟. 物理学报, 2011, 60(9): 090206. doi: 10.7498/aps.60.090206
    [19] 蒋涛, 欧阳洁, 赵晓凯, 任金莲. 黏性液滴变形过程的核梯度修正光滑粒子动力学模拟. 物理学报, 2011, 60(5): 054701. doi: 10.7498/aps.60.054701
    [20] 刘谋斌, 常建忠. 光滑粒子动力学方法中粒子分布与数值稳定性分析. 物理学报, 2010, 59(6): 3654-3662. doi: 10.7498/aps.59.3654
计量
  • 文章访问数:  5258
  • PDF下载量:  545
  • 被引次数: 0
出版历程
  • 收稿日期:  2014-10-09
  • 修回日期:  2014-11-07
  • 刊出日期:  2015-05-05

/

返回文章
返回