发明名称 一种验证离心泵网格质量与计算精度关系的方法
摘要 本发明公开了一种确定网格质量与计算精度之间关系的方法,设计离心泵及计算机辅助设计制造领域,本发明可以准确地验证所提出网格质量衡量准则的有效性,即网格质量越高,其计算精度越高。该方法通过多套不完善网格CFD计算外特性与试验值对比,并比较相邻两套网格所有单元的速度误差值,从而确定问题的网格无关解,并在此基础上,确定其他多套网格与网格无关解之间的各个单元的质量及误差,进而准确地获得网格质量与计算精度之间的关系。
申请公布号 CN102750414A 申请公布日期 2012.10.24
申请号 CN201210208537.8 申请日期 2012.06.25
申请人 江苏大学 发明人 刘厚林;董亮;代;王勇;王凯;谈明高
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 南京知识律师事务所 32207 代理人 汪旭东
主权项 1.一种验证离心泵网格质量与计算精度关系的方法,其特征在于,具体步骤为:(A)通过对离心泵进行外特性试验,获得离心泵的扬程、效率数据;(B)生成试验泵模型的四面体网格,定义为第<i>k</i>套网格,将该套网格进行CFD数值计算并记录网格中各节点的速度值,定义为第<i>k</i>套网格数值解;(C)按照1.2比值增加网格单元数,定义为第<i>k</i>+1套网格,将该套网格进行CFD数值计算并记录网格中各节点的速度值,定义为第<i>k</i>+1套网格数值解;(D)采用三维线性插值算法将<i>k</i>+1套网格数值解映射到<i>k</i>套网格数值解中,获得第<i>k</i>+1套网格在第<i>k</i>套网格中相应的数值解;(E)通过下面公式计算第<i>k</i>+1套网格数值解与第<i>k</i>套网格数值解的速度相对误差<img file="2012102085378100001DEST_PATH_IMAGE002.GIF" wi="26" he="26" />;<img file="2012102085378100001DEST_PATH_IMAGE004.GIF" wi="196" he="81" />,其中<i>V</i><sub><i>k</i>,<i>i</i></sub>为第<i>k</i>套网格节点<i>i</i>的速度,<i>V</i><sub><i>k</i>+1,<i>i</i></sub>为第<i>k+</i>1套网格映射后在节点<i>i</i>处的速度,<i>N</i>为<i>k</i>套网格节点总数;(F)如果<img file="703335DEST_PATH_IMAGE002.GIF" wi="26" he="26" />小于给定的误差值<i>ε=2%</i>,那么就认为第<i>k</i>套网格为“候选”网格,否则重复前面e步骤;(G)将<i>k</i>套网格计算得到的扬程和效率结果与试验的扬程和效率进行对比,其中杨程采用相对误差而效率采用绝对误差,扬程相对误差:<img file="2012102085378100001DEST_PATH_IMAGE006.GIF" wi="182" he="25" />, 式中:<i>H</i>为计算扬程<i>H</i><sub><i>e</i></sub>为试验扬程,效率绝对误差:<img file="2012102085378100001DEST_PATH_IMAGE008.GIF" wi="78" he="25" />,式中:<i>η</i>为计算效率<i>η</i><sub><i>e</i></sub>为试验效率当扬程相对误差和效率绝对误差均小于3%时,则认为该套网格的计算结果为该问题的网格独立解;(H)通过公式<img file="2012102085378100001DEST_PATH_IMAGE010.GIF" wi="469" he="80" />,获得不完善网格中各单元质量, 式中:<i>V</i>表示四面体的体积;<i>L</i><sub><i>ij</i></sub>表示顶点标号为<i>i</i>、<i>j</i>的四面体的边长;<i>S</i><sub><i>i</i></sub>表示三角形的面积;<img file="2012102085378100001DEST_PATH_IMAGE012.GIF" wi="116" he="57" />;(I)采用三维线性插值算法将不完善网格数值解映射到独立网格数值解中,获得网格独立解网格在该套套网格中相应的数值解,获得不完善网格数值解中各单元相对于独立解的速度相对误差值<img file="765617DEST_PATH_IMAGE002.GIF" wi="26" he="26" />;(J)将不完善网格中所有单元质量以及相对于独立解的速度相对误差值<img file="545354DEST_PATH_IMAGE002.GIF" wi="26" he="26" />建立一个数据库,且单元质量与速度相对误差值<img file="2880DEST_PATH_IMAGE002.GIF" wi="26" he="26" />一一对应;(K)按照单元的速度相对误差由小到大顺序方式对数据库中的数据重新排列,且单元质量与其相对应的速度相对误差值<img file="817253DEST_PATH_IMAGE002.GIF" wi="26" he="26" />重新排列;(L)将该数据库分为15-30个区间,计算各个区间的单元质量以及速度相对误差值<img file="639715DEST_PATH_IMAGE002.GIF" wi="26" he="26" />的平均值,即采用平均值形式统计区间的数据,速度相对误差<img file="969065DEST_PATH_IMAGE002.GIF" wi="26" he="26" />与单元质量都取区间的平均值,接着以各个区间内速度相对误差值<img file="167965DEST_PATH_IMAGE002.GIF" wi="26" he="26" />的平均值为横坐标,各个区间内单元质量的平均值为纵坐标,获得网格质量与计算精度之间的关系。
地址 212013 江苏省镇江市京口区学府路301号