SPH后处理研究

计算物理 ›› 2011, Vol. 28 ›› Issue (2): 213-218.

SPH后处理研究

郑俊, 于开平, 魏英杰, 张嘉钟

哈尔滨工业大学航天学院, 黑龙江哈尔滨 150001

收稿日期:2009-12-08 修回日期:2010-05-25 出版日期:2011-03-25 发布日期:2011-03-25
作者简介:郑俊(1983-),男,浙江,博士生,主要从事力学、计算力学及计算流体力学的研究.
基金资助:
国家自然科学重点基金(10832007)资助项目

A General Study on Post-processing of Smoothed Particle Hydrodynamics

ZHENG Jun, YU Kaiping, WEI Yingjie, ZHANG Jiazhong

School of Astronautics, Harbin Institute of Technology, Harbin 150001, China

Received:2009-12-08 Revised:2010-05-25 Online:2011-03-25 Published:2011-03-25

摘要/Abstract

摘要： 提出一套相对完整的光滑粒子流体动力学(SPH)后处理方法.凸区域和严重变形的非凸区域上的SPH计算结果都可用该方法处理.首先,对SPH粒子集进行Delaunay三角化,可得到由粒子作为节点的三角单元集;根据每个单元中三个节点是否彼此为粒子作用对,决定是否将该单元从单元集中删除.将保留下的单元作为有限单元并利用它们节点上的函数值,根据有限元插值方法即可得单元内部任何一点的函数值.根据该方法,可提取介质的自由表面.数值算例表明方法可行.对含固体壁面、严重粒子飞溅、多重介质互相作用的情况,提出相应对策;为SPH以至其它无网格方法形成较通用的后处理软件提供可行的途径.

关键词: SPH, 非凸区域, Delaunay三角化, 后处理

Abstract: We propose a method for data-visualization of smoothed particle hydrodynamics(SPH).Data calculated with SPH on convex or non convex continuum domain are processed.Particle set is discretized into a set of triangular elements by Delaunay triangulation,and only those elements on which nodes interact with each other are kept.With kept elements and data on their nodes a finite element interpolation is employed to derive values on interior points of the elements.A free surface is easily extracted with the method.Feasibility of the method is tested.Methods dealing with walls,severe particle splashing or multi-materials are provided,respectively.They supply feasible ways for realizing a general post-processing software package for SPH and even for gridless methods.

Key words: SPH, non convex domain, Delaunay triangulation, post processing

中图分类号:

O302

郑俊, 于开平, 魏英杰, 张嘉钟. SPH后处理研究[J]. 计算物理, 2011, 28(2): 213-218.

ZHENG Jun, YU Kaiping, WEI Yingjie, ZHANG Jiazhong. A General Study on Post-processing of Smoothed Particle Hydrodynamics[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2011, 28(2): 213-218.

[1]	徐丞君, 徐胜利. SPH二阶粒子近似光滑函数的充分条件及数值验证[J]. 计算物理, 2019, 36(1): 25-38.
[2]	上官子柠, 周秀丽, 宋鑫, 陈谦. 典型自由面流动问题的SPH-ALE数值模拟[J]. 计算物理, 2017, 34(6): 641-650.
[3]	周杰, 徐胜利. SPH方法气-液界面边界条件研究[J]. 计算物理, 2017, 34(4): 409-416.
[4]	赵亚洲, 马智博. SPH核函数光滑长度最优选取和自适应准则研究[J]. 计算物理, 2017, 34(1): 29-38.
[5]	周杰, 徐胜利. SPH方法进出口数值边界条件研究[J]. 计算物理, 2016, 33(5): 516-522.
[6]	李付鹏, 汪继文, 欧莽. 一种光滑粒子流体动力学数据的后处理方法[J]. 计算物理, 2015, 32(1): 33-39.
[7]	郑俊, 于开平, 王军锋, 李昌烽. SPH在水下高速物体空泡发展模拟中的应用[J]. 计算物理, 2014, 31(1): 27-32.
[8]	王庄, 李遇春, 王立时. 不同截面水槽流体参数晃动的SPH模拟[J]. 计算物理, 2013, 30(5): 642-648.
[9]	郑兴, 马庆位, 段文洋. K2_SPH方法及二维破碎波的模拟[J]. 计算物理, 2012, 29(3): 317-325.
[10]	张志春, 强洪夫, 傅学金, 夏薇. 基于SPH-FEM转换算法的7.62mm步枪弹冲击30CrMnSiA钢板的数值计算[J]. 计算物理, 2012, 29(1): 73-81.
[11]	郑兴, 段文洋. K2_SPH方法及其对二维非线性水波的模拟[J]. 计算物理, 2011, 28(5): 659-666.
[12]	任金莲, 欧阳洁, 蒋涛. 二维管道充填过程的修正SPH模拟[J]. 计算物理, 2011, 28(4): 515-522.
[13]	肖毅华, 胡德安, 韩旭. 一种有限元-光滑粒子流体动力学耦合算法[J]. 计算物理, 2011, 28(2): 219-224.
[14]	姚熊亮, 杨文山, 初文华, 明付仁. 舰船水下接触爆炸数值研究[J]. 计算物理, 2010, 27(6): 840-846.
[15]	陈炎, 曹树良, 梁开洪, 祝宝山. 结合前沿推进的Delaunay三角化网格生成及应用[J]. 计算物理, 2009, 26(4): 527-533.