离心泵非结构化四面体网格生成及应用研究

发布时间：2020-10-12 16:20

　　 网格是CFD计算和有限元分析的关键之一,其质量直接影响CFD数值计算的精度、收敛性和效率。由于离心泵几何外形复杂、曲率变化剧烈,广泛应用的非结构化四面体网格生成算法仍然面临诸多挑战,很难生成高质量的四面体网格。因此,如何精确而高效地生成与计算区域相适应的高质量四面体网格是离心泵数值模拟亟待解决的难题之一,这不仅能为离心泵数值计算提供高质量网格,促进网格生成理论和方法的发展,也能为其它类型泵及几何形状复杂机械部件高质量网格的生成提供借鉴。 本文在国家自然科学基金项目“基于RM-DT算法的水泵CFD非结构化网格生成研究”(No.51179075)的资助下,针对离心泵非结构化四面体网格生成过程中存在的问题展开研究,在理论分析、算例验证、数值模拟和试验研究的基础上提出了一整套适用于离心泵几何形状特点的四面体网格生成算法。本文的主要工作和创造性成果有： 1.系统地总结了国内外非结构化四面体网格生成技术的研究现状,重点阐述了现有网格质量衡量准则、空间曲面和三维实体网格划分以及网格优化算法存在的主要问题。归纳了离心泵内部流动数值模拟干扰因素的研究成果。此外,还对离心泵内部流动PIV测量做了较为全面的概括。 2.在对比分析不同造型方法对网格划分难易、生成网格的质量以及效率影响的基础上,确定了离心泵全计算域各过流部件的造型方法,并基于VC 6.0软件开发平台以及Pro/TOOLKIT,通过建立尺寸变量驱动的三维实体模板,开发了离心泵全计算域三维参数化自动造型及装配软件P3D。同时,为了完成设计参数输入和数据传递,采用动态连接库通信方式调用MFC,设计了直观易用的可视化交互界面。实例表明,程序运行效果良好,能大幅提高工作效率。 3.为增强评估四面体网格质量的能力,建立一种基于四面体单元几何特性的网格质量衡量准则。选取9种常用衡量准则进行比较,结果表明新准则在对劣质单元、单元形状变化的评判以及判别网格质量所耗费时间方面都显示出了较好的优势。在此基础上,结合弯管数值计算和PIV试验,探讨了构建网格独立解的方法。同时,采用新准则和其它3种准则,对网格质量衡量准则与计算精度的关系进行了深入分析。研究表明,在计算误差为E,e0.15的区间内,误差的变化能够引起新衡量准则度量值的变化,故该衡量准则能够正确地评判出对计算精度影响较大的劣质单元。 4.针对离心泵空间曲面网格生成过程中映射畸变问题,首次提出了一种基于黎曼度量并结合AFT和Delaunay法的曲面网格生成方法,即采用AFT来生成参数域的背景网格,且在对参数域内背景网格进行加密时,为了抵消映射过程的畸变,采用插入三角形外接椭圆中心的方法代替传统Delaunay插点内核。该方法通过黎曼度量的引入以及合理的背景网格充分考虑了映射畸变和算法效率,不同比转数离心泵叶片实例表明,算法易于实施,稳定性好且生成的曲面网格质量较高。 5.针对三维实体网格Sliver单元消除问题,在归纳分析加密消除Sliver单元算法原理的基础上,对一种较优的消除算法进行了改进,优化了其判断Sliver单元的方法,并通过保证插入新点后不产生细小Sliver单元的方式,解决了原方法对于插入点位置要求过于苛刻的问题。算例表明,本文算法能够完全消除离心泵中的Sliver单元。 6.针对网格优化过程中边界网格质量难以控制的问题,提出了一种基于优化算法的光顺(Optimization-based smoothing, OPS)对比分析了不同目标函数和求最优解方法在网格数、迭代次数以及需求精度等因素变化下对OPS优化效果的影响,选用一种最优的目标函数,并将其变换为包含有考虑边界网格质量和内部网格质量的函数项,且为函数项添加了一个权重系数C,从而实现了对内部和边界处网格单元质量的同时控制。通过实例验证,该方法不但能够提高内部网格质量且能够提高边界处网格质量,权重系数C取0.7能够取得较好的优化效果且效率较高。 7.针对OPS算法本身固有的无法提高“全边界单元”质量的缺陷,首次提出了一种结合OPS和拓扑优化的网格质量提高策略,该方法通过判别“全边界单元”中边(或面)被相邻单元共享的形式,而采用不同的拓扑交换方式,且在拓扑交换过程中仅考虑优化后网格的拓扑相关性,并不考虑连通区域的质量,保证能够将此类单元转换为非“全边界单元”,从而在随后的OPS优化中获得较好的优化效果。 8.采用CFD计算与PIV试验方式验证四面体网格生成算法的有效性。选取16台不同比转数的离心泵,探讨了不同计算域、湍流模型、进出口边界条件、速度压力耦合算法以及收敛精度等干扰因素对数值计算的影响,以确定最佳数值计算方法。在此基础上,应用某商用网格划分软件和本文算法生成网格进行了数值计算,并与离心泵PIV试验进行了比较。结果表明网格数量相近的情况下,本文算法得到的最差单元质量以及网格平均质量都要好于商用网格划分软件。对于外特性而言,两套网格计算的扬程和效率误差均在3.5%以内,但采用本文算法生成的网格预测精度更高；对于内流场而言,两套网格都与试验值存在一些的差异,数值上的差异尤为明显,但本文算法生成的网格无论是在预测离心泵内部流动趋势还是在数值上都要好于商用网格划分软件。
【学位单位】：江苏大学
【学位级别】：博士
【学位年份】：2012
【中图分类】：TH311
【部分图文】：

图1一1间接法曲面网格生成RM算法在消除映射畸变方面具有很好的效果，能够准确地度量映射畸变，并且在合理定义黎曼度量后可以很好地控制待划分域网格的生成必一叨，如图1一1所示为文献口7]中间接法曲面网格生成过程。RM算法通过曲面本身信息和网格划分技术对生成的网格节点进行控制，其关键是如何获得当前节点的控制度量「28]。曲面上每一点的控制度量都是不一样的，网格划分过程中时刻需要得到曲面上某一点的控制度量。若每一次求取点的控制度量时都调用相关接口函数，那么算法的运行效率将会变得很低。一个有效的解决方法是采用背景网格法，即首先给定背景网格，接着将曲面上若干点的控制度量存储到背景网格中，而对于其它点的控制度量，可根据背景网格中已有点的控制度量插值得到卿]。对于复杂几何模型，这样的工作通常困难且十分费时，对用户的经验和专业素养要求也很高。同时，在离心泵曲面网格划分过程中

网的同时保证边界边在最终网格中存在，而更通用的做法是在Deaunay三角化后执行边界恢复算法来保证边界一致。相比二维边界恢复问题，三维边界恢复则存在理论上的困难，Sch6nhardt多面体(如图1一2所示)的存在从理论上否定了在不加点情形下Delaunay三角化任意多面体的可能脚]。Ruppe[40]则证明在不加点的情况下判断任意多面体能否三角化是一个Np(Non一 deterministiepolynomial)问题。因此，要实现边界恢复必定要考虑如何添加辅助点，添加辅助点不仅破坏了原有单元尺寸场，还可能在其附近形成质量较差，甚至产生退化的单元。因此，对添加辅助点的要求是[4l]:一要尽量减少辅助点的数量;二要优化辅助点附近的单元质量。c玲一子一厂处图1一 2Sch6nhardt多面体

出现将导致在CFD计算时出现刚度矩阵病态，严重影响计算结果的精度，甚至导致无法收敛。Slive:单元也称薄元，是指体积趋近于零的四面体单元，其4个顶点几乎共面，如图1一3所示。当A趋近于A‘时就会形成Sliver单元。它是唯一一种同时具有四面体外接球半径与最短边的比值(Redius一 edge)147]较小但四面体外接球半径与内接球半径的比值 (aspeetratio)较大的畸形单元。通常情况下，Radius一ed罗比值小于2且最小二面角小于5“的四面体，就认为是Sliver单元。图1一 3Slive:单元定义Slive;单元的消除是非结构化网格生成中一个非常重要的问题，国内外学者提出了多种消除方法，主要分为通过优化消除Slive;单元和通过加密消除sliver单元两种方法。第一种方法是对己生成的网格进行优化处理，从而消除网格中存在的Slive;单元。oiaehin、smith、ohadyani等￡48一50]通过对s一iver单元的相邻单元进行局部变换来消除网格中的Sliver单元。Escobars、vartziotis等〔5’一52]提出采用助plaeian光顺和广义薄元分解相结合的优化算法来消
【相似文献】

相关期刊论文 前10条

1 李世堃;;离心泵蜗壳蜗线绘形方法的研究[J];机械;1987年03期

2 李海生,杨钦,陈其明;根据STL文件生成三维四面体网格[J];中国机械工程;2003年11期

3 赵建军,王启付;边界一致的Delaunay四面体网格稳定生成算法[J];机械工程学报;2004年06期

4 胡于进,赵虎跃,赵建军;基于Delaunay准则的三维网格自动插点算法[J];华中理工大学学报;2000年05期

5 刘建国,习华勇,王玉昆;基于MDT的离心泵蜗壳三维造型研究与实现[J];机械设计与制造;2005年05期

6 李世煌;许建中;;离心泵蜗壳内流动的高速摄影测试[J];北京农业工程大学学报;1990年01期

7 李锡夔，方玉玲;三维域四面体网格生成的推进网阵方法和网格光顺[J];计算力学学报;1994年01期

8 王德生,万水;与波前法相结合的Delaunay四面体网格生成方法(英文)[J];Journal of Southeast University;2002年02期

9 刘明晖;;四面体网格修匀算法研究[J];建筑技术开发;2009年10期

10 李安虎 ,李红 ,刘炜巍;离心泵蜗壳内流场的数值模拟[J];通用机械;2004年11期

相关博士学位论文 前5条

1 董亮;离心泵非结构化四面体网格生成及应用研究[D];江苏大学;2012年

2 杨孙圣;离心泵作透平的理论分析数值计算与实验研究[D];江苏大学;2012年

3 王春林;旋流自吸泵设计及内部流动研究[D];江苏大学;2010年

4 陈斌;低比转速无过载排污泵优化设计及实验研究[D];江苏大学;2012年

5 李伟;斜流泵启动过程瞬态非定常内流特性及实验研究[D];江苏大学;2012年

相关硕士学位论文 前10条

1 薛冰军;离心泵内部流场的数值模拟与实验分析[D];河北工程大学;2012年

2 孟瑞锋;基于UG的离心泵蜗壳参数化三维造型软件二次开发的研究[D];兰州理工大学;2012年

3 田素垒;从表面模型生成四面体网格的算法研究[D];西北农林科技大学;2012年

4 周剑松;基于CFD分析的CQB65-50-160F离心泵的性能预测及优化[D];浙江工业大学;2012年

5 许文静;低比转速复合叶轮离心泵非定常流场的数值模拟[D];浙江理工大学;2012年

6 陈延华;面向有限元的四面体网格生成算法研究[D];济南大学;2011年

7 马鹏飞;转速变化对离心泵性能影响的数值模拟研究及调速调节节能分析[D];兰州理工大学;2011年

8 狄驰;KMAS系统中三维四面体网格生成器及相关理论研究[D];吉林大学;2006年

9 万忠;低比转速复合叶轮离心泵的优化设计和试验研究[D];浙江理工大学;2011年

10 崔凌国;约束Delaunay四面体剖分及其相关算法的研究[D];西北工业大学;2006年

本文编号：2838005