一种光滑粒子流体动力学数据的后处理方法.pdfVIP

第 32卷 第 1期 计 算 物 理 V0I_32．NO．1 2015年 1月 CHINESE JOURNALOFCOMPUTAT10NALPHYSICS Jan．，2015 文章编号 ：1001—246X(2015)01-0033—07 一 种光滑粒子流体动力学数据的后处理方法 李付鹏 ， 汪继文 ， 欧 莽 (1．安徽大学计算智能与信号处理教育部重点实验室，安徽 合肥 230039 2．安徽大学计算机科学与技术学院，安徽 合肥 230039) 摘 要：提出一种光滑流体动力学数据后处理的三角网格化方法．首先确定流场中每个粒子作用域 内的配对粒子 集；再以当前粒子的配对粒子集为限定条件进行Delaunay三角剖分 ，可得到由粒子作为节点的三角形单元集；最后 对 “虚假”单元做过滤处理 ．算例表明复杂边界条件和严重变形的流体 区域上 的SPH计算结果都可用该方法处理 ． 关键词 ：光滑粒子方法 ；非 凸区域 ；Delaunay三角化 ；后处理 中图分类号 ：0302 文献标志码 ：A 0 引言 近年来 ，光滑粒子流体动力学(smoothedparticlehydrodynamics，简称为 SPH)方法…逐渐成为解决物理 问题的一个有效工具 ，与传统的基于网格的数值方法 ，如有限差分 、有 限元、有限体积方法不 同，这是一种典 型的无 网格拉格 朗 日数值方法，它利用核函数对物理问题进行近似处理 ，用离散的粒子来描述宏观连续分布 微观仍为粒子的流体，而每个粒子则携带了其所在位置的流体的各种性质 ，如质量 、密度 、速度 、能量等 ，算法 节点可以任意分布，并非有某种网格联系在一起，这种 自适应的特点使得该方法在模拟流体剧烈变形或间断 问题时具有显著优势，目前方法的应用 已扩展到气体动力学、不可压缩 、爆炸、固体力学和弹性体等多个领 域．然而 ，与该算法应用快速发展所不相适应 的是 ，方法的数据后处理技术一直没有太大的发展，甚至可以说 是一大难题 ．如果将计算后的结果直接图形化，得到 的是大量无组织的离散粒子 ，无法表现 出流体 的性 态，更无法对流体的等值线 (面)、流线图等分析，现有的比较成熟的基于网格的可视化技术也无法直接利 用．目前，就我们查阅到的文献而言，有基于 Splash、Paraview等软件开发的SPH后处理程序 ，但这类软件 较为专业 ，通用性不好 ，同时在流体力学方面的处理能力有限． 文献E4—5]提出了一种处理 SPH数据后处理的方法 ，该方法先利用 Delaunay三角剖分技术对离散数据 三角化 ，然后根据每个单元 中三个节点是否彼此为粒子作用对 ，决定是否将该单元从单元集中删除，从算例 来看，这种方法是可行的．本文根据 SPH算法的特点，对上述方法进行了改进 ，提出了一种更为简便的SPH 数据后处理方法 ，以粒子支持域 (影响域)为限定条件进行限定三角剖分 (constrainedDelaunaytriangulation， 简称为 CDT或 conformingDelaunaytriangulation，简称为RDT)，与文献 [4—5]相 比，方法减少了所生成的多 余的三角形网格的数量，删除的空白的三角形 网格也随之减少 ，提高了算法运算 的效率．这样做另一个优点 在于三角网格化输出的图形能够清晰地体现各粒子之间的影响关系 ，凡与某个节点(粒子)相关联的其它节 点(粒子)均为数据输 出时刻与该粒子有相互影响的粒子，而且方法更加简便易行，一般情况下不需要过滤 “虚假”单元 ，简化了算法流程 ，增强了算法 的适应性 ，验证结果表明方法稳定性好 ，凸区域和严重变形的非 凸区域上的SPH计算结果都可用该方法处理． 首先介绍 SPH方法中与数据三角化相关的几个概念 以及数据后处理的特点 ；接下来阐述 SPH数据三角 化的算法思想和步骤 ；最后通过算例验证 ，得出结论． 收稿El期 ：2014—0l一24；修回 日期 ：2014—06—18 基金项 目：安徽高校省级 自然科学研究 (KJ2013A009)资助项 目 作者简介 ：李付鹏 (1974一)，男，安徽阜阳，博士生，从事计算流体力学、计算机流体动画模拟研究，E-mail：lifupengan