一种基于cfd仿真和网格自适应的阀门流量系数计算方法
【专利说明】一种基于CFD仿真和网格自适应的阀门流量系数计算方法 所属技术领域
[0001] 本发明涉及到计算阀门流通能力的领域，具体表现为应用CFD与网格自适应相结 合的方法计算阀门流量系数。
【背景技术】
[0002] 阀门的流量系数是指单位时间内，在保持恒定的压力条件下管道内介质流经阀门 的总流量。流量系数反映了阀门的流通能力和节能环保性能，是阀门重要的工艺参数和技 术指标。目前国内的很多阀门生产厂家以实验方法测得的流量系数作为标准衡量阀门的流 通能力，但是实验成本高，周期长，耗费大量的人力物力。
[0003] 随着计算机技术和计算流体动力学（ComputationalFluidDynamics，CFD)的发 展，基于CH)的数值仿真可以准确、快速而且可视化地反映复杂流场的细节特征。运用CH)方 法进行阀门流通能力的分析将是现有实验方法的有效补充，同时也有助于设计人员了解阀 门结构对内部流场的影响规律。
[0004] 计算规模大、精度要求高是阀门CFD仿真的两大主要特点。众所周知，网格是影响 CFD仿真精度和计算效率的重要因素。随着网格数量的增加，结果精度一般也会提高，但是 计算时间也会随之增加，所以在划分网格时需要综合考虑精度和效率两个方面。在阀门内 流场的CFD仿真中，通常需要对模型的不同局部设置不同的网格密度，并进行大量手工加密 处理。这些经验性的设置和处理不仅工作量较大，而且会对计算结果产生较大的影响，导致 计算精度存在一定的不确定性，限制了CFD方法在阀门制造行业的推广应用。
[0005] 本发明旨在利用网格自适应技术对阀门流场中敏感的区域进行有针对性的自动 加密，从而提高仿真精度。该方法相比人工调整流场网格更加简单易行，而且对操作者的专 业程度依赖性不高，有利于在工业界推广应用。
[0006] 经检索，目前有关应用CFD仿真方法衡量阀门流通能力的申报专利有 CN103729505A，这篇专利提出了基于CFD仿真的阀门当量长度计算方法。但是这篇专利只是 利用通用的方法求解阀门两端的压差，进而计算得到阀门当量长度，并没有考虑到网格对 仿真精度的影响，更未提到应用网格自适应的方法降低人工优化网格的难度和工作量。
【发明内容】
[0007] 本发明的目的在于克服运用实验方法计算阀门流量系数成本高、周期长、对阀门 尺寸规格有限制以及现有CH)仿真方法人工加密网格、计算精度受限于操作人员专业程度 等问题，提出了一种基于CH)仿真和网格自适应的阀门流量系数计算方法。此方法可以在阀 门设计阶段起到指导作用，提高效率，节约资源，大大的提高仿真精度，降低了人工加密网 格的工作量和难度。
[0008] 实现本发明目的的技术方案为： 利用三维软件建立阀门的装配体模型，并且抽取"阀门-管道"流场模型和"直管道"流 场模型，对两个模型分别进行混合网格和结构化网格划分，然后利用商用CFD软件FLUENT在 多种不同边界条件下结合网格自适应技术进行仿真计算，得到模型进出口的压差值，再通 过公式计算得到阀门的流量系数。
[0009] 为提高计算精度，降低了人工加密网格的工作量和难度，本发明在仿真过程中采 用以下自适应操作： (1)先利用FLUENT软件对初始网格进行试运算； (2 )进行第一次网格自适应，包括y+自适应和速度梯度自适应，其中在y+自适应中根据 y +取值范围设置粗化阀值和加密阀值，在速度梯度自适应中选取的自适应方法为 Gradient，标准化方式为Standard，根据自适应函数等值线云图确定速度梯度变化较大区 域，得到加密阀值； (3)若前一次计算结果的y+值不满足30〈y+〈300范围，则继续采用y+自适应对近壁面网 格进行不同次数的优化。
[0010] 由于自适应后网格会大量增加，故选取y+的取值范围要根据流场模型的大小和计 算机性能综合考虑，一般在工程应用中y+值的取值范围不应超出30〈 7+〈300，若模型比较 小，计算机性能又足够高，使11.5〈y+〈60，则可以获得更高的仿真精度。
[0011] 本发明具体实现步骤如下： 步骤（1):建立流场模型和网格划分 建立阀门装配体三维模型，并且分别在阀门进、出口连接5倍公称直径和10倍公称直径 长度的上、下游管道，抽取"阀门-管道"流场模型，并建立15倍公称直径长度的"直管道"流 场模型。
[0012] 对"阀门-管道"流场模型进行多种不同规模的混合网格划分，即两边管道区域划 分为结构网格，中间阀门区域划分为非结构网格;将"直管道"流场模型全部划分为结构网 格。这种划分网格的方式可以获得质量更高的初始网格。
[0013] 步骤(2):判断最优初始网格 针对"阀门-管道"流场模型，分别应用不同规模的初始网格进行仿真试运算，并分析所 得的"阀门-管道"模型进出口压差随网格规模的变化趋势，当在某初始网格基础上网格数 量增加10%而计算所得的压差变化量不超过1.5%时，将该初始网格选取为最优初始网格。
[0014] 步骤(3):不同边界条件下"阀门-管道"压差 利用步骤(2)中得到的最优初始网格，考虑不同的流动边界条件，结合多次网格自适应 操作，仿真计算各边界条件下的"阀门-管道"模型的进出口压差。
[0015] 步骤(4):不同边界条件下"直管道"压差 针对"直管道"流场模型，设置与步骤(3)中相对应的边界条件，结合网格自适应操作， 仿真计算各边界条件下的"直管道"模型的进出口压差。
[0016] 步骤(5):计算阀门流量系数 将"阀门-管道"压差与"直管道"压差的差值代入阀门流量系数计算公式心 中进行计算，得到相应的阀门流量系数，式中为流量（单位：），ΔΡν为阀门的净压差（单 位:），为密度(单位:），常温下。
[0017] 步骤(6):计算流量系数平均值 将各组边界条件下最后一次自适应计算得到的流量系数取算术平均值，并将此均值作 为阀门流量系数的预测值。
[0018] 本发明的有益效果是： (1)相比于实验方法测量阀门流量系数，本发明所述方法能够在设计初期准确预测阀 门的流通能力，而且能够提供远比实验更加丰富的流场细节信息，向设计人向直观地展示 阀门结构对流场的影响，因此可以在阀门的设计阶段起到很好的指导作用，缩短研发周期， 节约成本，避免人力物力的浪费。
[0019] (2)相比于现有应用CH)仿真对阀门流通能力进行检测的方法，本发明通过应用网 格自适应技术大大提高了仿真精度。利用y+自适应技术可以有效的改善边界层网格，利用 速度梯度自适应技术则可对流场中急变流区域网格进行自动加密，避免了人工加密网格的 盲目性，同时也降低了工作量和难度。该方法对操作者优化网格的专业程度依赖性不高，有 利于在工业界普遍推广。
[0020] (3)本发明在阀门进出口分别延长了5倍和10倍公称直径长度的管道区域，这种做 法可以让流体介质充分发展，且可以有效避免回流现象的发生。在仿真过程中用"阀门-管 道"压差Δ Pi和"直管道"压差Δ P2的差值Δ Pv对阀门流量系数进行计算，这种方法有效地避 免了长直管道产生的压差对计算结果的影响。
【附图说明】
[0021 ]图1是本发明所述流量系数计算方法流程图； 图2是同一初始网格不同边界条件自适应前后相对误差变化折线图。
【具体实施方式】 [0022] 图1是本发明所述流量系数计算方法流程图，下面提供本发明所述的一种基于CH)仿真 和网格自适应的阀门流量系数计算方法的【具体实施方式】。
[0023] 实施例1: 本实施例的操作步骤如下： (1)建立流场模型和网格划分 以DN500偏心蝶阀为例，利用Solidworks建立装配体模型，对此模型中的倒圆、倒角和 螺纹孔等细小结构进行简化后在进出口分别连接5倍公称直径和10倍公称直径长度的管 道;然后抽取"阀门-管道"流场模型，总长度为7857毫米，并且建立15倍公称直径长度的"直 管道"流场模型，总长度为7500毫米。
[0024]利用ANSYSICEMCFD前处理软件将"阀门-管道"流场模型划分为混合网格，将"直管 道"流场模型全部划分为结构网格。
[0025] (2)判断最优初始网格 针对"阀门-管道"流场模型设置2.3m/s的进口流速，在45万、140万、311万、383万、427 万五种不同规模的网格下进行仿真运算。仿真得到的压差如表1所示。
[0026]表1不同网格规模下仿真结果
每一次增加网格数量仿真得到压差之间的相对误差分别为2.464%、2.110%、1.629%和 0.358%，利用427万网格计算得到的压差相对误差低于1.5%，故选取最优网格数量为427万。 [0027] (3)不同边界条件下"阀门-管道"压差 针对"阀门-管道"流场模型，在427万网格数量下分别设置入口速度为2.3m/s、2.88m/ s、3.17m/s、3.33m/s，并且对每一组实验进行2次网格自适应，共进行12次仿真计算。不同边 界条件下仿真得到的"阀门-管道"压差如表2所示，单位为Pa。
[0028] 表2不同进口速度下"阀门-管道"压差仿真结果_
(4)不同边界条件下"直管道"压差 针对"直管道"流场模型，在373万网格数量下设置与步