CN113553737B - 一种基于阀门压差的阀门流量预测方法 - Google Patents

一种基于阀门压差的阀门流量预测方法 Download PDF

Info

Publication number
CN113553737B
CN113553737B CN202110723676.3A CN202110723676A CN113553737B CN 113553737 B CN113553737 B CN 113553737B CN 202110723676 A CN202110723676 A CN 202110723676A CN 113553737 B CN113553737 B CN 113553737B
Authority
CN
China
Prior art keywords
valve
flow
pressure
model
grid
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110723676.3A
Other languages
English (en)
Other versions
CN113553737A (zh
Inventor
林哲
杨林
陶俊宇
田帅
朱祖超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang Sci Tech University ZSTU
Original Assignee
Zhejiang Sci Tech University ZSTU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang Sci Tech University ZSTU filed Critical Zhejiang Sci Tech University ZSTU
Priority to CN202110723676.3A priority Critical patent/CN113553737B/zh
Publication of CN113553737A publication Critical patent/CN113553737A/zh
Application granted granted Critical
Publication of CN113553737B publication Critical patent/CN113553737B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/14—Pipes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Measuring Volume Flow (AREA)

Abstract

本发明公开了一种基于阀门压差的阀门流量预测方法，包括以下步骤：确定流体物性；根据所选流经阀门的流体，确定流体的物性参数，如密度(Kg/m3)、动力粘性(m2/s)，确定阀门流量系数曲线方程；阀门流量系数曲线方程一般通过两种方法获得，一种是查阀门手册，可以快速读取阀门流量系数曲线方程；另外一种可以通过数值模拟的方法得到阀门流量系数曲线方程。有益效果：转换为流经阀门流量的变化，与传统的流量计测流量相比，本方法不用在流场内部放置检测装置，对流体的扰动很小，不会改变检测区域流体的流量，且压力传感器对压力的变化比流量计对流量的变化敏感。因此，该方法对流量的预测具有干扰小、响应快的特点。

Description

一种基于阀门压差的阀门流量预测方法
技术领域
本发明涉及流体流量预测技术领域,具体来说，涉及一种基于阀门压差的阀门流量预测方法。
背景技术
流量计是计量科学技术的组成部分之一，它与国民经济、国防建设、科学研究有密切的关系。流量测量工作，对保证产品质量、提高生产效率、促进科学技术的发展都具有重要作用。传统的流量测量，如差压式流量计、转子流量计、节流式流量计等，都需要将装置伸入流场内部，对于一些流体流动剧烈复杂的流场，这些测量装置会对流场造成剧烈影响，破坏流场内部的流动特性，不仅得不到准确的流量信息，还改变了流场的速度分布。
特别地，在一些狭窄微小阀门管道空间，不足以放置传统的流量计，或者放置的流量计强烈干扰测量区域的流场流量，导致测量的流量与真实值的偏差很大。因此运用不伸入流场、较小扰动的基于阀门压差的阀门流量预测方法，预测阀门流道流体流量显得较为合适可靠.
针对相关技术中的问题，目前尚未提出有效的解决方案。
发明内容
针对相关技术中的问题，本发明提出一种基于阀门压差的阀门流量预测方法，以克服现有相关技术所存在的上述技术问题。
为此，本发明采用的具体技术方案如下：
一种基于阀门压差的阀门流量预测方法，包括以下步骤：
S101、确定流体物性；根据所选流经阀门的流体，确定流体的物性参数，如密度(Kg/m3)、动力粘性(m2/s)；
S103、确定阀门流量系数曲线方程；阀门流量系数曲线方程一般通过两种方法获得，一种是查阀门手册，可以快速读取阀门流量系数曲线方程；另外一种通过数值模拟的方法得到阀门流量系数曲线方程；
S1031、建立流道模型；建立所要模拟阀门模型，设置阀门相对开度M为5％，对模型流体区域进行抽取，获得流道模型，流道模型特征在于：所述流道模型为进口管长为五倍管径、出口管长为十倍管径的蝶阀；
S1032、网格无关性验证；对流道模型进行网格划分，比较不同网格模型下流道中的流量系数，选取符合要求的网格模型；
S105、测量阀前压力P1 *和阀后压力后P2 *，从而计算阀门前后的压力差，计算公式如下：
其中ΔP表示流体流经阀门的压力差，单位为百帕斯卡；P1 *表示测得的流体流经阀前五倍管径处的压力，单位为帕斯卡；P2 *表示测得的流体流经阀后十倍管径处的压力，单位为帕斯卡；
S107、预测阀门流量；确定阀门相对开度M，根据阀门流量系数曲线方程确定流量系数Kv *，从而计算得到阀门的流量Q*，计算公式如下：
其中Q*表示预测的阀门流量，单位为立方米每小时；Kv *表示根据所得阀门流量系数曲线方程得到阀门流量系数；ρ*表示流体密度，单位为千克每立方米，ΔP表示流体流经阀前五倍管径和阀后十倍管径处的压差，单位为百帕斯卡。
作为优选的，所述步骤S103中选取符合要求的网格模型具体实施步骤如下：
S1033、根据流道模型的局部特征，绘制不同数量网格数的流道网格，确保网格质量大于0.3，将所有流道网格按网格数按升序进行排序，比较相邻网格数下的流量系数，相邻网格数下流量系数的偏差在以内的网格即为符合要求的网格模型；
S1034、数值模拟计算，将符合要求的网格模型导入ANSYS-FLUENT软件进行数值模拟计算；
S1035、设置边界条件；当阀门实际工作条件明确时，根据实际工作条件设置边界条件；当阀门实际工作条件不明确时，设置进口边界条件为压力进口，压力值为0.1MPa，出口边界条件为压力出口，压力值为0MPa，其余边界设置为无滑移的壁面边界条件；
S1036、求解设置；求解方法选用SIMPLE算法，其中SIMPLE算法中的梯度离散方式采用基于单元体的最小二乘法，SIMPLE算法中的压力、动量、能量离散方式均采用二阶迎风格式，初始化设置为全局初始化，设置求解残差和迭代步数，确保求解收敛有效；
S1037、输出计算结果；在ANSYS-FLUENT软件后处理中，导出阀门流量系数和进出口压力。
作为优选的，所述步骤S1034具体实施步骤如下:
S10341、定义求解器；求解器类型采用基于压力的求解器，时间类型设置为稳态计算方式，求解器的速度方程采用相对速度，操作压力设置为当地大气压；
S10342、选择计算模型；根据流体的流动状态，选择不同的计算模型，其步骤如下：
S103421、当雷诺数低于2300，选择层流模型，使用Laminar层流模型；
S103422、当雷诺数大于2300，选择Standard湍流模型。
作为优选的，所述步骤S1037中输出计算结果方法步骤如下：
S10371、选择Results模块中的Reports模块，其中ReportType选择VolumeFlowRate，Surfaces选择出口面，即可显示其体积流量Q；
S10372、选择Results模块中的Reports模块，其中Surface Integrals选择Surface Integrals，在ReportType模块选择FacetAverage，在Surface Integrals模块右侧的FieldVariable模块选择StaticPressure,在Surfaces模块分别选择出口面和进口面，单击Compute在ANSYS-FLUENT的Console窗口即可显示其进口面压力P1和出口面压力P2；
S1037、计算阀门流量系数；通过ANSYS-FLUENT软件导出的体积流量、进出口压力和流体密度计算得到该阀门开度下的流量系数，计算公式如下：
Kv表示计算所得流量系数；Q表示ANSYS-FLUENT软件导出的阀门流量，单位为立方米每秒；ρ表示流体密度，单位为千克每立方米；P1表示测得的流体流经阀前五倍管径处的压力，单位为帕斯卡；P2表示测得的流体流经阀后十倍管径处的压力，单位为帕斯卡；
S1038、得到不同开度下阀门流量系数曲线方程；在已有相对开度的基础上，第一次将阀门相对开度M增加5％，之后在前者相对开度M的基础上增加10％，不断重复所述S103中S1031～S1034的步骤直到阀门相对开度M为100％，获取5％、10％、20％、30％、40％、50％、60％、70％、80％、90％和100％相对开度M下的阀门流量系数，将所得相对开度下的流量系数进行五阶函数拟合，从而得到不同开度下阀门流量系数曲线方程。
作为优选的，通过对管道以及管道上的阀门进行三维建模，其特征在于：阀前选取长度五倍于管径的管道，阀后十倍于管径的管道，根据所建模型，使用ANSYS Workbench中的DesigenModeler抽取阀门流道模型。
作为优选的，其压差为阀门前五倍管径处与阀门后十倍管径处的压力差值，即上述模型进口和出口的压力差值。
作为优选的，对流道进行网格划分，其特征在于：所述S1032的网格划分，使用ANSYS Workbench的Meshing对流道进行网格划分，网格类型采用非结构网格，且保证最差网格的网格质量大于0.3。
本发明的有益效果为：通过监测阀门进出口的压差变化，转换为流经阀门流量的变化，与传统的流量计测流量相比，本方法不用在流场内部放置检测装置，对流体的扰动很小，不会改变检测区域流体的流量，且压力传感器对压力的变化比流量计对流量的变化敏感。因此，该方法对流量的预测具有干扰小、响应快的特点。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是根据本发明实施例的一种基于阀门压差的阀门流量预测方法的流程图；
图2是根据本发明实施例的一种基于阀门压差的阀门流量预测方法中蝶阀流道模型图；
图3是根据本发明实施例的一种基于阀门压差的阀门流量预测方法中蝶阀流道模型网格图；
图4是根据本发明实施例的一种基于阀门压差的阀门流量预测方法中蝶阀流量系数网格无关性折线图；
图5是根据本发明实施例的一种基于阀门压差的阀门流量预测方法中蝶阀流量系数曲线方程图。
具体实施方式
为进一步说明各实施例，本发明提供有附图，这些附图为本发明揭露内容的一部分，其主要用以说明实施例，并可配合说明书的相关描述来解释实施例的运作原理，配合参考这些内容，本领域普通技术人员应能理解其他可能的实施方式以及本发明的优点，图中的组件并未按比例绘制，而类似的组件符号通常用来表示类似的组件。
根据本发明的实施例，提供了一种基于阀门压差的阀门流量预测方法。
实施例一；
如图1-5所示，根据本发明实施例的基于阀门压差的阀门流量预测方法，包括以下步骤：
步骤一：确定流体物性；水被广泛使用在阀门管道中，是较为常见的流体。因此，选择水作为流过阀门的流量预测流体，则流体的密度ρ＝998Kg/m3、动力粘度为v＝1.01×10- 3Pa·s；
步骤二：确定阀门流量系数；阀门流量系数一般运用两种方法，一种是查阀门手册，可以快速读取阀门流量系数；另外一种可以通过数值模拟的方法得到流量系数，其特征在于：方法步骤如下:
①建立流道模型；建立所要模拟阀门模型，本例为蝶阀(不仅限于蝶阀)，其阀门公称直径为1600mm进口管长8000mm内径1600mm，出口管长16000mm内径1600mm，设置阀门相对开度为5％，对模型流体区域进行抽取，获得流道模型，流道模型如图2所示；
②网格无关性验证；对流道模型进行网格划分，比较不同网格模型下流道中的流量系数，选取符合要求的网格模型，其网格无关性验证结果如图4所示；
所述步骤②中选取符合要求的网格模型具体实施步骤如下：
S1:根据流道模型的局部特征，绘制不同数量网格数的流道网格，确保网格质量大于0.3。将所有流道网格按网格数按升序进行排序，比较相邻网格数下的流量系数，相邻网格数下流量系数的偏差在以内的网格即为符合要求的网格模型；
③数值模拟计算
将符合要求的网格模型如图3所示，导入ANSYS-FLUENT软件进行数值模拟计算；
所述步骤③具体实施步骤如下:
S1:定义求解器；
求解器类型采用基于压力的求解器，时间类型设置为稳态计算方式，求解器的速度方程采用相对速度，操作压力设置为当地大气压。；
S2:选择计算模型；
根据流体的流动状态，选择不同的计算模型，其步骤如下：
①当雷诺数低于2300，选择层流模型，使用Laminar层流模型；
②当雷诺数大于2300，选择湍流模型，由于Standard湍流模型的稳定性，因此将Standard湍流模型作为计算模型；
S3:设置边界条件；当阀门实际工作条件明确时，根据实际工作条件设置边界条件；当阀门实际工作条件不明确时，设置进口边界条件为压力进口(压力值为0.1MPa)，出口边界条件为压力出口(压力值为0MPa)，其余边界设置为无滑移的壁面边界条件；
S4:求解设置；求解方法旋转SIMPLE算法，其中SIMPLE算法中的梯度离散方式采用基于单元体的最小二乘法，SIMPLE算法中的压力、动量、能量离散方式均采用二阶迎风格式，初始化设置为全局初始化，设置求解残差为，迭代步数为10000，确保求解收敛有效；
S5:输出计算结果；在ANSYS-FLUENT软件后处理中，导出阀门流量系数和进出口压力，其特征在于：方法步骤如下：
I.选择Results模块中的Reports模块，其中Report Type选择Volume Flow Rate，Surfaces选择出口面，即可显示其体积流量Q＝0.8689m3/s；
II.选择Results模块中的Reports模块，其中Surface Integrals选择SurfaceIntegrals，在Report Type模块选择Facet Average，在Surface Integrals模块右侧的Field Variable模块选择Static Pressure,在Surfaces模块分别选择出口面和进口面，单击Compute在ANSYS-FLUENT的Console窗口即可显示其进口面压力P1＝100000Pa和出口面压力P2＝0Pa；
④计算阀门流量系数；通过ANSYS-FLUENT软件导出的体积流量、进出口压力和流体密度计算得到该阀门开度下的流量系数，计算公式如下：
Kv表示计算所得流量系数；Q表示ANSYS-FLUENT软件导出的阀门流量，单位为立方米每秒；ρ表示流体密度，单位为千克每立方米；P2表示测得的流体流经阀前五倍管径处的压力，单位为帕斯卡；表示测得的流体流经阀后十倍管径处的压力，单位为帕斯卡；
⑤统计不同开度下的阀门流量系数
在已有相对开度M的基础上，第一次将阀门相对开度M增加5％，之后在前者相对开度M的基础上增加10％，不断重复所述步骤二中①～④的步骤直到阀门相对开度M为100％，获取M为5％、10％、20％、30％、40％、50％、60％、70％、80％、90％和100％的阀门流量系数，将所得相对开度下的流量系数进行五阶函数拟合，从而得到不同开度下阀门流量系数曲线方程，如图5所示；2.测量阀门前后压差ΔP；当阀门相对开度为100％时，测得阀前五倍管径处的压力P1 *为0.09934MPa、阀后十倍管径处的压力P2 *为0.0007017MPa，则流体在阀门处的压降可按照如下公式计算，公式如下：
其中ΔP表示流体流经阀门的压力差，单位为百帕斯卡；P1 *表示测得的流体流经阀前五倍管径处的压力，单位为帕斯卡；P2 *表示测得的流体流经阀后十倍管径处的压力，单位为帕斯卡；
3.预测阀门流量；在实际M＝100％工况下，通过在所得阀门流量系数曲线方程得到阀门流量系数Kv *＝179711.44，通过实际工作流体密度以及测得的流经阀门前后流体的压差，计算的得到阀门的流量Q*，计算公式如下：
其中Q*表示预测的阀门流量，单位为立方米每小时；Kv *表示根据所得流量特性曲线或流量特性曲线的五阶拟合方程得到阀门流量系数；ρ*表示流体密度，单位为千克每立方米，ΔP表示流体流经阀前五倍管径和阀后十倍管径处的压差，单位为百帕斯卡。
4.当改变相对阀门开度，可以读取不同阀门相对开度下的流量系数。
按照相同计算方式，当相对开度为20％、40％、60％、80％时，在所得阀门流量系数曲线上读取所述蝶阀流量系数为10000、24000、49800、101000，测量阀的前后压差，通过计算可以对该阀门的流量进行预测，各相对开度M下阀门的预测流量如下表1所示。
表1不同相对开度下流量预测表；
相对开度M 20％ 40％ 60％ 80％
流量系数Kv * 10000 24000 49800 101000
阀门前后压差(Pa×102) 999.9 999.62 998.73 993.84
预测流量(m3/h) 10009.51 24019.47 49818.2 100789.3
综上所述，借助于本发明的上述技术方案，通过监测阀门进出口的压差变化，转换为流经阀门流量的变化，与传统的流量计测流量相比，本方法不用在流场内部放置检测装置，对流体的扰动很小，不会改变检测区域流体的流量，且压力传感器对压力的变化比流量计对流量的变化敏感。因此，该方法对流量的预测具有干扰小、响应快的特点。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (7)

1.一种基于阀门压差的阀门流量预测方法，其特征在于，包括以下步骤：
S101、确定流体物性；根据所选流经阀门的流体，确定流体的物性参数，如密度(Kg/m3)、动力粘性(m2/s)；
S103、确定阀门流量系数曲线方程；阀门流量系数曲线方程一般通过两种方法获得，一种是查阀门手册，可以快速读取阀门流量系数曲线方程；另外一种通过数值模拟的方法得到阀门流量系数曲线方程；
S1031、建立流道模型；建立所要模拟阀门模型，设置阀门相对开度M为5％，对模型流体区域进行抽取，获得流道模型，流道模型特征在于：所述流道模型为进口管长为五倍管径、出口管长为十倍管径的蝶阀；
S1032、网格无关性验证；对流道模型进行网格划分，比较不同网格模型下流道中的流量系数，选取符合要求的网格模型；
S105、测量阀前压力P1 *和阀后压力后P2 *，从而计算阀门前后的压力差，计算公式如下：
其中ΔP表示流体流经阀门的压力差，单位为百帕斯卡；P1 *表示测得的流体流经阀前五倍管径处的压力，单位为帕斯卡；P2 *表示测得的流体流经阀后十倍管径处的压力，单位为帕斯卡；
S107、预测阀门流量；确定阀门相对开度M，根据阀门流量系数曲线方程确定流量系数Kv *，从而计算得到阀门的流量Q*，计算公式如下：
其中Q*表示预测的阀门流量，单位为立方米每小时；Kv *表示根据所得阀门流量系数曲线方程得到阀门流量系数；ρ*表示流体密度，单位为千克每立方米，ΔP表示流体流经阀前五倍管径和阀后十倍管径处的压差，单位为百帕斯卡。
2.根据权利要求1所述的一种基于阀门压差的阀门流量预测方法，其特征在于，选取符合要求的网格模型具体实施包括以下步骤；
S1033、根据流道模型的局部特征，绘制不同数量网格数的流道网格，确保网格质量大于0.3，将所有流道网格按网格数按升序进行排序，比较相邻网格数下的流量系数，相邻网格数下流量系数的偏差在以内的网格即为符合要求的网格模型；
S1034、数值模拟计算，将符合要求的网格模型导入ANSYS-FLUENT软件进行数值模拟计算；
S1032、网格无关性验证；对流道模型进行网格划分，比较不同网格模型下流道中的流量系数，选取符合要求的网格模型；
S1035、设置边界条件；当阀门实际工作条件明确时，根据实际工作条件设置边界条件；当阀门实际工作条件不明确时，设置进口边界条件为压力进口，压力值为0.1MPa，出口边界条件为压力出口，压力值为0MPa，其余边界设置为无滑移的壁面边界条件；
S1036、求解设置；求解方法选用SIMPLE算法，其中SIMPLE算法中的梯度离散方式采用基于单元体的最小二乘法，SIMPLE算法中的压力、动量、能量离散方式均采用二阶迎风格式，初始化设置为全局初始化，设置求解残差和迭代步数，确保求解收敛有效；
S1037、输出计算结果；在ANSYS-FLUENT软件后处理中，导出阀门流量系数和进出口压力。
3.根据权利要求2所述的一种基于阀门压差的阀门流量预测方法，其特征在于，所述步骤S1034具体实施步骤如下:
S10341、定义求解器；求解器类型采用基于压力的求解器，时间类型设置为稳态计算方式，求解器的速度方程采用相对速度，操作压力设置为当地大气压；
S10342、选择计算模型；根据流体的流动状态，选择不同的计算模型，其步骤如下：
S103421、当雷诺数低于2300，选择层流模型，使用Laminar层流模型；
S103422、当雷诺数大于2300，选择Standard湍流模型。
4.根据权利要求2所述的一种基于阀门压差的阀门流量预测方法，其特征在于，所述步骤S1037中输出计算结果方法步骤如下：
S10371、选择Results模块中的Reports模块，其中Report Type选择Volume FlowRate，Surfaces选择出口面，即可显示其体积流量Q；
S10372、选择Results模块中的Reports模块，其中Surface Integrals选择SurfaceIntegrals，在Report Type模块选择Facet Average，在Surface Integrals模块右侧的Field Variable模块选择Static Pressure,在Surfaces模块分别选择出口面和进口面，单击Compute在ANSYS-FLUENT的Console窗口即可显示其进口面压力P1和出口面压力P2；
S1037、计算阀门流量系数；通过ANSYS-FLUENT软件导出的体积流量、进出口压力和流体密度计算得到该阀门开度下的流量系数，计算公式如下：
Kv表示计算所得流量系数；Q表示ANSYS-FLUENT软件导出的阀门流量，单位为立方米每秒；ρ表示流体密度，单位为千克每立方米；P1表示测得的流体流经阀前五倍管径处的压力，单位为帕斯卡；P2表示测得的流体流经阀后十倍管径处的压力，单位为帕斯卡；
S1038、得到不同开度下阀门流量系数曲线方程；在已有相对开度的基础上，第一次将阀门相对开度M增加5％，之后在前者相对开度M的基础上增加10％，不断重复所述步骤S103中S1031～S1034的步骤直到阀门相对开度M为100％，获取5％、10％、20％、30％、40％、50％、60％、70％、80％、90％和100％相对开度M下的阀门流量系数，将所得相对开度下的流量系数进行五阶函数拟合，从而得到不同开度下阀门流量系数曲线方程。
5.根据权利要求4所述的一种基于阀门压差的阀门流量预测方法，其特征在于，通过对管道以及管道上的阀门进行三维建模，其特征在于：阀前选取长度五倍于管径的管道，阀后十倍于管径的管道，根据所建模型，使用ANSYS Workbench中的DesigenModeler抽取阀门流道模型。
6.根据权利要求5所述的一种基于阀门压差的阀门流量预测方法，其特征在于，其压差为阀门前五倍管径处与阀门后十倍管径处的压力差值，即上述模型进口和出口的压力差值。
7.根据权利要求6所述的一种基于阀门压差的阀门流量预测方法，其特征在于：所述S1032的网格划分，使用ANSYS Workbench的Meshing对流道进行网格划分，网格类型采用非结构网格，且保证最差网格的网格质量大于0.3。
CN202110723676.3A 2021-06-28 2021-06-28 一种基于阀门压差的阀门流量预测方法 Active CN113553737B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110723676.3A CN113553737B (zh) 2021-06-28 2021-06-28 一种基于阀门压差的阀门流量预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110723676.3A CN113553737B (zh) 2021-06-28 2021-06-28 一种基于阀门压差的阀门流量预测方法

Publications (2)

Publication Number Publication Date
CN113553737A CN113553737A (zh) 2021-10-26
CN113553737B true CN113553737B (zh) 2023-10-13

Family

ID=78102458

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110723676.3A Active CN113553737B (zh) 2021-06-28 2021-06-28 一种基于阀门压差的阀门流量预测方法

Country Status (1)

Country Link
CN (1) CN113553737B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113970365A (zh) * 2021-10-29 2022-01-25 华能重庆两江燃机发电有限责任公司 一种大管径大流量测量***及方法
CN114323435B (zh) * 2021-12-13 2023-10-20 潍柴动力股份有限公司 压差传感器可信性检测方法、装置、电子设备及存储介质
CN115193342B (zh) * 2022-07-04 2024-04-05 中国科学院过程工程研究所 一种稳定浆料态物料流量的方法及实现其的装置与应用

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105677964A (zh) * 2016-01-07 2016-06-15 江苏神通阀门股份有限公司 一种基于cfd仿真和网格自适应的阀门流量系数计算方法
CN111578371A (zh) * 2020-05-22 2020-08-25 浙江大学 一种数据驱动的城市集中供热***精准调控方法
CN111767663A (zh) * 2020-05-29 2020-10-13 江苏神通阀门股份有限公司 一种基于cfd仿真的阀门流量系数便捷计算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9715237B2 (en) * 2013-08-26 2017-07-25 Spraying Systems Co. Flow-control valve system and method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105677964A (zh) * 2016-01-07 2016-06-15 江苏神通阀门股份有限公司 一种基于cfd仿真和网格自适应的阀门流量系数计算方法
CN111578371A (zh) * 2020-05-22 2020-08-25 浙江大学 一种数据驱动的城市集中供热***精准调控方法
CN111767663A (zh) * 2020-05-29 2020-10-13 江苏神通阀门股份有限公司 一种基于cfd仿真的阀门流量系数便捷计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
崔宝玲 ; 马光飞 ; 王慧杰 ; 林哲 ; 尚照辉.阀芯结构对节流截止阀流阻特性和内部流动特性的影响.机械工程学报.2015,(第012期),178-184. *
许洪斌 ; 刁富强 ; 杨长辉 ; 秦飞龙 ; .滑板式调节阀流场的数值模拟及性能预测.机床与液压.2015,(第19期),187+197. *

Also Published As

Publication number Publication date
CN113553737A (zh) 2021-10-26

Similar Documents

Publication Publication Date Title
CN113553737B (zh) 一种基于阀门压差的阀门流量预测方法
CN105042339B (zh) 一种基于无量纲的成品油管道泄漏量估计***及方法
CN105181040B (zh) 一种差压式流量计的数字化标定及优化方法
CN103729505B (zh) 一种基于cfd的阀门当量长度计算方法
CN110197049B (zh) 一种基于瞬变反问题的非金属管道泄漏定位方法
Aly et al. Experimental study of the pressure drop after fractal-shaped orifices in turbulent pipe flows
CN103074873B (zh) 一种非理想流场下的多声道超声流量计声道布置方法
CN111767663A (zh) 一种基于cfd仿真的阀门流量系数便捷计算方法
CN108019622A (zh) 一种基于压力差的管道泄露定位的计算方法
CN106503396A (zh) 基于有限差分法与有限体积法耦合的多维水力***瞬变模拟方法
El Drainy et al. CFD analysis of incompressible turbulent swirling flow through zanker plate
Sreekala et al. Study of flow performance of a globe valve and design optimisation
Meng et al. A flow stability evaluation method based on flow-pressure correlation
CN108280300B (zh) 基于计算流体力学的大数据流量计开发方法
CN116698141A (zh) 一种不同工况下超声波流量计的测速误差修正方法及***
CN114811448B (zh) 一种流动条件下管道泄漏检测、泄露流速估计与泄露定位的方法
CN108506622B (zh) 一种基于弧形导流片的低阻力三通构件
Hao et al. Influence of the hole chamfer on the characteristics of a multi-hole orifice flowmeter
CN113536711B (zh) 一种基于田口法的橄榄形流量计结构优化设计方法
Hariri Asli et al. Isolated pressure zones based on GIS as a solution for water network problems
CN103323066A (zh) 一种低含液率气液两相流测量方法及测量***
CN112836447A (zh) 一体化交互式管路瞬变流动演示***及方法
Menge Analysis of turbulent flow in a pipe at constant reynolds number using computational fluid dynamics
Sravani et al. Transient behavior analysis for orifice flow meter performance for variation in sensor dimension
CN106202743A (zh) 一种复杂汽、水管道压力取样位置的确定方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant