[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN113591414B - Rectifying effect judging method suitable for flow valve - Google Patents

Rectifying effect judging method suitable for flow valve Download PDF

Info

Publication number
CN113591414B
CN113591414B CN202110912459.9A CN202110912459A CN113591414B CN 113591414 B CN113591414 B CN 113591414B CN 202110912459 A CN202110912459 A CN 202110912459A CN 113591414 B CN113591414 B CN 113591414B
Authority
CN
China
Prior art keywords
rectifier
flow
module
downstream
section
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
CN202110912459.9A
Other languages
Chinese (zh)
Other versions
CN113591414A (en
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 CN202110912459.9A priority Critical patent/CN113591414B/en
Publication of CN113591414A publication Critical patent/CN113591414A/en
Application granted granted Critical
Publication of CN113591414B publication Critical patent/CN113591414B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/14Pipes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force 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)
  • Control Of Positive-Displacement Air Blowers (AREA)

Abstract

The invention discloses a rectifying effect judging method suitable for a flow valve, which comprises the following steps: establishing a flow channel model; three-dimensional modeling is carried out on a valve and a rectifier under actual working conditions by utilizing Solidworks software, a rectifier model is assembled at an outlet of the valve, a model fluid domain is extracted, a flow channel model is obtained, and the flow channel model is characterized in that: the upstream of the flow channel model valve is five-pipe diameter, and the downstream of the rectifier is twelve-pipe diameter. The beneficial effects are that: the rectifying condition after the valve is reflected by a numerical simulation method, the speed eccentricity of the downstream cross section of the rectifier is utilized to judge whether the rectifying effect of the rectifier with different structures is good or not, and complex experiments are not needed to judge, so that manpower and material resources can be greatly reduced, and the method has practical significance in engineering.

Description

一种适用于流程阀门的整流效果判定方法A method for determining the rectification effect suitable for process valves

技术领域Technical field

本发明涉及整流效果判定技术领域,具体来说,涉及一种适用于流程阀门的整流效果判定方法。The present invention relates to the technical field of rectification effect determination, and specifically, to a rectification effect determination method suitable for process valves.

背景技术Background technique

流程阀门在流体传输过程中起着越来越重要的作用,在石油化工、水利水电等行业得到了广泛的应用。在不同的工况下具备调节、分流、启闭等功能。然而,流程阀门内部的急变流与高剪切流引起的漩涡等复杂不稳定流动,会沿下游管道发展,影响流程阀门运行可靠性及管道流体准确计量。Process valves play an increasingly important role in the fluid transmission process and are widely used in petrochemical, water conservancy and hydropower industries. It has functions such as adjustment, diversion, opening and closing under different working conditions. However, complex and unstable flows such as rapid flow changes and vortices caused by high shear flow inside the process valve will develop along the downstream pipeline, affecting the operational reliability of the process valve and the accurate measurement of pipeline fluids.

工程中往往会在管道流体计量设备前安装整流器来获得较为稳定的计量信号。整流器是加速不规则流体发展、缓减流场畸变对流量计影响的设备,可在较小的空间内采集稳定的流量计量信号,大大提高流量计量结果的准确性。安装整流器是当前流场优化和流量计性能提升的重要途经。In engineering projects, rectifiers are often installed in front of pipeline fluid metering equipment to obtain a more stable metering signal. The rectifier is a device that accelerates the development of irregular fluids and mitigates the impact of flow field distortion on the flow meter. It can collect stable flow measurement signals in a small space and greatly improve the accuracy of flow measurement results. Installing a rectifier is an important way to optimize the current flow field and improve the performance of the flow meter.

但是,在实际工程中,整流器整流效果的好坏通常根据实验测量得到的结果来判定,这花费了大量的人力和物力,所得到的实验结果还存在误差,不一定准确。However, in actual engineering, the rectification effect of the rectifier is usually judged based on the results obtained from experimental measurements, which consumes a lot of manpower and material resources, and the experimental results obtained still have errors and are not necessarily accurate.

针对相关技术中的问题,目前尚未提出有效的解决方案。No effective solutions have yet been proposed for the problems in related technologies.

发明内容Contents of the invention

针对相关技术中的问题,本发明提出一种适用于流程阀门的整流效果判定方法,以克服现有相关技术所存在的上述技术问题。In view of the problems in the related technology, the present invention proposes a rectification effect determination method suitable for process valves to overcome the above technical problems existing in the existing related technology.

为此,本发明采用的具体技术方案如下:To this end, the specific technical solutions adopted by the present invention are as follows:

一种适用于流程阀门的整流效果判定方法,包括以下步骤:A method for determining the rectification effect suitable for process valves, including the following steps:

S101、建立流道模型;对实际工况下的阀门和整流器利用Solidworks软件进行三维建模,将整流器模型装配在阀门出口处,对模型流体域进行抽取,获得流道模型,流道模型特征在于:所述流道模型阀门上游为五倍管径,整流器下游为十二倍管径;S101. Establish a flow channel model; use Solidworks software to conduct three-dimensional modeling of the valve and rectifier under actual working conditions, assemble the rectifier model at the valve outlet, extract the model fluid domain, and obtain the flow channel model. The characteristics of the flow channel model are: : The upstream of the flow channel model valve is five times the pipe diameter, and the downstream of the rectifier is twelve times the pipe diameter;

S103、划分网格;进行网格无关性验证,选取符合要求的网格数,对流道模型进行网格划分;S103. Mesh division; perform grid independence verification, select the number of meshes that meet the requirements, and mesh the flow channel model;

S105、进行数值计算;将流道网格导入ANASYS FLUENT软件,设置相关参数,并进行数值计算;S105. Perform numerical calculations; import the flow channel grid into the ANASYS FLUENT software, set relevant parameters, and perform numerical calculations;

S107、计算下游截面的速度偏心率η;提取整流器下游管道横截面上的最大速度,根据公式,计算下游每个横截面的速度偏心率η,计算公式如下: S107. Calculate the velocity eccentricity η of the downstream section; extract the maximum velocity on the cross-section of the pipeline downstream of the rectifier, and calculate the velocity eccentricity eta of each downstream cross-section according to the formula. The calculation formula is as follows:

其中,Vmax表示下游管道横截面上的最大速度,单位为m/s;V0表示管道内的平均速度,单位为m/s。Among them, Vmax represents the maximum velocity on the cross section of the downstream pipeline, in m/s; V0 represents the average velocity in the pipeline, in m/s.

S109,判断整流效果;根据计算出来的速度偏心率,利用Origin软件绘制出速度偏心率随流动发展的曲线图,对比分析不同结构整流器的整流效果。S109, determine the rectification effect; based on the calculated speed eccentricity, use the Origin software to draw a curve chart of the speed eccentricity developing with the flow, and compare and analyze the rectification effects of rectifiers with different structures.

作为优选的,所述划分网格具体实施包括以下步骤:Preferably, the specific implementation of meshing includes the following steps:

S1031、在Solidworks软件里面,将流道模型的每一个部分分别导出;S1031. In the Solidworks software, export each part of the flow channel model separately;

S1032、利用ANSYS ICEM CFD软件,对流道模型的每一个部分分别划分结构性网格,并对阀芯和整流器的部分进行网格加密,最后利用ANSYS ICEM CFD软件的Merge功能,将每个部分的网格合并成一个整体。S1032. Use ANSYS ICEM CFD software to divide each part of the flow channel model into a structural mesh, and refine the mesh for the valve core and rectifier parts. Finally, use the Merge function of ANSYS ICEM CFD software to merge each part. The grids merge into a whole.

作为优选的,所述网格无关性验证具体实施包括以下步骤:Preferably, the specific implementation of the grid independence verification includes the following steps:

S1033、针对阀门全开、安装某一种整流器时的流道模型,进行网格划分,并划分成多套具有不同网格数量的流道网格模型;S1033. For the flow channel model when the valve is fully open and a certain rectifier is installed, mesh it and divide it into multiple sets of flow channel mesh models with different mesh numbers;

S1034、比较不同网格数下的流量系数和阻力系数,当流量系数和阻力系数随网格数量的增加而基本保持不变时,选择流量系数和阻力系数基本保持不变的最少网格数。S1034. Compare the flow coefficient and resistance coefficient under different grid numbers. When the flow coefficient and resistance coefficient remain basically unchanged as the number of grids increases, select the minimum number of grids where the flow coefficient and resistance coefficient remain basically unchanged.

作为优选的,所述步骤S105中设置相关参数具体实施包括以下步骤:Preferably, the specific implementation of setting relevant parameters in step S105 includes the following steps:

S1051、在General模块中定义求解器,求解器类型采用基于压力的求解器,时间类型设置为稳态计算方式;S1051. Define a solver in the General module. The solver type is a pressure-based solver, and the time type is set to the steady-state calculation method;

S1052、在Models模块中选择计算模型,选择Standard k-ε模型进行计算;S1052. Select the calculation model in the Models module and select the Standard k-ε model for calculation;

S1053、在Materials模块中设置流体介质,根据实际工况,选择流体介质为水或者气;S1053. Set the fluid medium in the Materials module, and select the fluid medium as water or gas according to the actual working conditions;

S1054、在Boundary Conditions模块中设置边界条件,入口选择velocity-inlet,出口选择pressure-outlet,壁面选择wall,交界面选择interface;当实际工作条件明确时,根据实际工作条件设置边界条件;当实际工作条件不明确时,自行设定合理的边界条件;S1054. Set boundary conditions in the Boundary Conditions module. Select velocity-inlet for the inlet, pressure-outlet for the outlet, wall for the wall, and interface for the interface. When the actual working conditions are clear, set the boundary conditions according to the actual working conditions; when the actual working conditions When the conditions are unclear, set reasonable boundary conditions by yourself;

S1055、在Mesh Interfaces模块中定义interface;S1055. Define interface in the Mesh Interfaces module;

S1056、在Residual模块中设置收敛残差;S1056. Set the convergence residual in the Residual module;

S1057、在Initialization模块中点击Hybrid Initialization进行初始化;S1057. Click Hybrid Initialization in the Initialization module to initialize;

S1058、在Run Calculation模块中设置迭代步数后点击Calculation模块开始进行数值计算。S1058. After setting the number of iteration steps in the Run Calculation module, click the Calculation module to start numerical calculations.

作为优选的,所述的步骤S107中提取整流器下游管道横截面上的最大速度具体实施包括以下步骤:Preferably, the specific implementation of extracting the maximum speed on the cross-section of the pipeline downstream of the rectifier in step S107 includes the following steps:

S1071、在Surface模块中选择Iso-Surface,在整流器下游每间隔1D创建一个横截面,共创建12个横截面;S1071. Select Iso-Surface in the Surface module and create a cross section every 1D downstream of the rectifier, creating a total of 12 cross sections;

S1072、在Reports模块中的Surface Integrals模块,Report Type选择FacetMaximum,Field Variable选择Velocity Magnitude,Surfaces选择S1071步骤中在整流器下游创建的横截面,即可得到横截面上的最大速度。S1072. In the Surface Integrals module in the Reports module, select FacetMaximum for Report Type, select Velocity Magnitude for Field Variable, and select the cross section created downstream of the rectifier in step S1071 for Surfaces to get the maximum velocity on the cross section.

作为优选的,所述的判断整流效果的具体方法是:横截面的速度偏心率越小,说明截面上的速度分布越均匀,流场越稳定,整流器的整流效果越好。根据Origin软件绘制出速度偏心率随流动发展的曲线图可以直观地看到不同结构整流器下游截面速度偏心率的大小,以此来判定不同结构整流器整流效果的好坏。Preferably, the specific method for judging the rectification effect is: the smaller the velocity eccentricity of the cross-section, the more uniform the velocity distribution on the cross-section, the more stable the flow field, and the better the rectification effect of the rectifier. Drawing the curve of velocity eccentricity with flow development based on the Origin software can visually see the velocity eccentricity of the downstream cross-section of rectifiers with different structures, so as to determine the rectification effect of rectifiers with different structures.

本发明的有益效果为:通过数值模拟的方法,对阀后的整流情况进行了反映,利用整流器下游横截面的速度偏心率判定不同结构整流器整流效果的好坏,不需要进行复杂的实验来判断,可以大大减少人力物力,在工程上具有实际意义。The beneficial effects of the present invention are: through the method of numerical simulation, the rectification situation behind the valve is reflected, and the speed eccentricity of the downstream cross-section of the rectifier is used to determine the rectification effect of the rectifiers with different structures, without the need to conduct complex experiments to judge. , can greatly reduce manpower and material resources, and has practical significance in engineering.

附图说明Description of the drawings

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to explain the embodiments of the present invention or the technical solutions in the prior art more clearly, the drawings needed to be used in the embodiments will be briefly introduced below. Obviously, the drawings in the following description are only some of the drawings of the present invention. Embodiments, for those of ordinary skill in the art, other drawings can also be obtained based on these drawings without exerting creative efforts.

图1是根据本发明实施例的一种适用于流程阀门的整流效果判定方法的流程图;Figure 1 is a flow chart of a method for determining the rectification effect of a process valve according to an embodiment of the present invention;

图2是根据本发明实施例的一种适用于流程阀门的整流效果判定方法中整流器I的二维结构图;Figure 2 is a two-dimensional structural diagram of the rectifier I in a method for determining the rectification effect of a process valve according to an embodiment of the present invention;

图3是根据本发明实施例的一种适用于流程阀门的整流效果判定方法中整流器II的二维结构图;Figure 3 is a two-dimensional structural diagram of the rectifier II in a method for determining the rectification effect of a process valve according to an embodiment of the present invention;

图4是根据本发明实施例的一种适用于流程阀门的整流效果判定方法中球阀及整流器流道模型图;Figure 4 is a model diagram of the ball valve and rectifier flow channel in a rectification effect determination method suitable for process valves according to an embodiment of the present invention;

图5是根据本发明实施例的一种适用于流程阀门的整流效果判定方法中球阀及整流器流道模型网格图;Figure 5 is a grid diagram of the ball valve and rectifier flow channel model in a rectification effect determination method suitable for process valves according to an embodiment of the present invention;

图6是根据本发明实施例的一种适用于流程阀门的整流效果判定方法中流量系数和阻力系数的网格无关性柱形图;Figure 6 is a grid-independent column chart of the flow coefficient and the resistance coefficient in a method for determining the rectification effect of a process valve according to an embodiment of the present invention;

图7是根据本发明实施例的一种适用于流程阀门的整流效果判定方法中两种整流器下游截面的速度偏心率随流动发展的曲线图。Figure 7 is a graph showing the velocity eccentricity of the downstream sections of two rectifiers as the flow develops in a rectification effect determination method suitable for process valves according to an embodiment of the present invention.

具体实施方式Detailed ways

为进一步说明各实施例,本发明提供有附图,这些附图为本发明揭露内容的一部分,其主要用以说明实施例,并可配合说明书的相关描述来解释实施例的运作原理,配合参考这些内容,本领域普通技术人员应能理解其他可能的实施方式以及本发明的优点,图中的组件并未按比例绘制,而类似的组件符号通常用来表示类似的组件。In order to further explain each embodiment, the present invention provides drawings. These drawings are part of the disclosure of the present invention. They are mainly used to illustrate the embodiments and can be used with the relevant descriptions in the specification to explain the operating principles of the embodiments. For reference From this, those of ordinary skill in the art will be able to understand other possible implementations and advantages of the present invention. The components in the figures are not drawn to scale, and similar component symbols are generally used to represent similar components.

根据本发明的实施例,提供了一种适用于流程阀门的整流效果判定方法。According to an embodiment of the present invention, a method for determining the rectification effect suitable for process valves is provided.

实施例一;Embodiment 1;

如图1-7所示,根据本发明实施例的适用于流程阀门的整流效果判定方法,包括以下步骤:As shown in Figures 1-7, the rectification effect determination method suitable for process valves according to the embodiment of the present invention includes the following steps:

步骤S101、建立流道模型;对实际工况下的阀门和整流器利用Solidworks软件进行三维建模,将整流器模型装配在阀门出口处,对模型流体域进行抽取,获得流道模型,流道模型特征在于:所述流道模型阀门上游为五倍管径,整流器下游为十二倍管径;Step S101. Establish a flow channel model; use Solidworks software to conduct three-dimensional modeling of the valve and rectifier under actual working conditions, assemble the rectifier model at the valve outlet, extract the model fluid domain, and obtain the flow channel model and flow channel model characteristics. The method is: the upstream of the flow channel model valve is five times the pipe diameter, and the downstream of the rectifier is twelve times the pipe diameter;

步骤S103、划分网格;进行网格无关性验证,选取符合要求的网格数,对流道模型进行网格划分;Step S103: Mesh division; perform mesh independence verification, select the number of meshes that meet the requirements, and mesh the flow channel model;

步骤S105、进行数值计算;将流道网格导入ANASYS FLUENT软件,设置相关参数,并进行数值计算;Step S105: Perform numerical calculations; import the flow channel grid into the ANASYS FLUENT software, set relevant parameters, and perform numerical calculations;

步骤S107、计算下游截面的速度偏心率η;提取整流器下游管道横截面上的最大速度,根据公式,计算下游每个横截面的速度偏心率η,计算公式如下: Step S107: Calculate the velocity eccentricity η of the downstream section; extract the maximum velocity on the cross section of the pipeline downstream of the rectifier, and calculate the velocity eccentricity eta of each downstream cross section according to the formula. The calculation formula is as follows:

其中,Vmax表示下游管道横截面上的最大速度,单位为m/s;V0表示管道内的平均速度,单位为m/s。Among them, Vmax represents the maximum velocity on the cross section of the downstream pipeline, in m/s; V0 represents the average velocity in the pipeline, in m/s.

步骤S109,判断整流效果;根据计算出来的速度偏心率,利用Origin软件绘制出速度偏心率随流动发展的曲线图,对比分析不同结构整流器的整流效果。Step S109, determine the rectification effect; based on the calculated speed eccentricity, use the Origin software to draw a curve chart of the speed eccentricity developing with the flow, and compare and analyze the rectification effects of rectifiers with different structures.

所述划分网格具体实施包括以下步骤:The specific implementation of the meshing includes the following steps:

步骤S1031、在Solidworks软件里面,将流道模型的每一个部分分别导出;Step S1031. In the Solidworks software, export each part of the flow channel model separately;

步骤S1032、利用ANSYS ICEM CFD软件,对流道模型的每一个部分分别划分结构性网格,并对阀芯和整流器的部分进行网格加密,最后利用ANSYS ICEM CFD软件的Merge功能,将每个部分的网格合并成一个整体。Step S1032: Use ANSYS ICEM CFD software to divide each part of the flow channel model into a structural mesh, and perform mesh densification on the valve core and rectifier parts. Finally, use the Merge function of ANSYS ICEM CFD software to merge each part. The grids are merged into a whole.

所述网格无关性验证具体实施包括以下步骤:The specific implementation of the grid independence verification includes the following steps:

步骤S1033、针对阀门全开、安装某一种整流器时的流道模型,进行网格划分,并划分成多套具有不同网格数量的流道网格模型;Step S1033. For the flow channel model when the valve is fully open and a certain type of rectifier is installed, meshing is performed and divided into multiple sets of flow channel mesh models with different mesh numbers;

步骤S1034、比较不同网格数下的流量系数和阻力系数,当流量系数和阻力系数随网格数量的增加而基本保持不变时,选择流量系数和阻力系数基本保持不变的最少网格数。Step S1034: Compare the flow coefficient and resistance coefficient under different grid numbers. When the flow coefficient and resistance coefficient remain basically unchanged as the number of grids increases, select the minimum number of grids where the flow coefficient and resistance coefficient remain basically unchanged. .

所述步骤S105中设置相关参数具体实施包括以下步骤:The specific implementation of setting relevant parameters in step S105 includes the following steps:

步骤S1051、在General模块中定义求解器,求解器类型采用基于压力的求解器,时间类型设置为稳态计算方式;Step S1051: Define a solver in the General module. The solver type is a pressure-based solver, and the time type is set to the steady-state calculation method;

步骤S1052、在Models模块中选择计算模型,选择Standard k-ε模型进行计算;Step S1052: Select the calculation model in the Models module and select the Standard k-ε model for calculation;

步骤S1053、在Materials模块中设置流体介质,根据实际工况,选择流体介质为水或者气;Step S1053: Set the fluid medium in the Materials module, and select the fluid medium as water or gas according to the actual working conditions;

步骤S1054、在Boundary Conditions模块中设置边界条件,入口选择velocity-inlet,出口选择pressure-outlet,壁面选择wall,交界面选择interface;当实际工作条件明确时,根据实际工作条件设置边界条件;当实际工作条件不明确时,自行设定合理的边界条件;Step S1054. Set boundary conditions in the Boundary Conditions module. Select velocity-inlet for the inlet, pressure-outlet for the outlet, wall for the wall, and interface for the interface. When the actual working conditions are clear, set the boundary conditions according to the actual working conditions. When working conditions are unclear, set reasonable boundary conditions by yourself;

步骤S1055、在Mesh Interfaces模块中定义interface;Step S1055: Define interface in the Mesh Interfaces module;

步骤S1056、在Residual模块中设置收敛残差;Step S1056: Set the convergence residual in the Residual module;

步骤S1057、在Initialization模块中点击Hybrid Initialization进行初始化;Step S1057: Click Hybrid Initialization in the Initialization module to initialize;

步骤S1058、在Run Calculation模块中设置迭代步数后点击Calculation模块开始进行数值计算。Step S1058: Set the number of iteration steps in the Run Calculation module and click the Calculation module to start numerical calculation.

所述的步骤S107中提取整流器下游管道横截面上的最大速度具体实施包括以下步骤:The specific implementation of extracting the maximum speed on the cross-section of the downstream pipe of the rectifier in step S107 includes the following steps:

步骤S1071、在Surface模块中选择Iso-Surface,在整流器下游每间隔1D创建一个横截面,共创建12个横截面;Step S1071. Select Iso-Surface in the Surface module and create a cross section every 1D downstream of the rectifier, creating a total of 12 cross sections;

步骤S1072、在Reports模块中的Surface Integrals模块,Report Type选择FacetMaximum,Field Variable选择Velocity Magnitude,Surfaces选择S1071步骤中在整流器下游创建的横截面,即可得到横截面上的最大速度。Step S1072. In the Surface Integrals module in the Reports module, select FacetMaximum for Report Type, select Velocity Magnitude for Field Variable, and select the cross section created downstream of the rectifier in Step S1071 for Surfaces to obtain the maximum velocity on the cross section.

所述的判断整流效果的具体方法是:横截面的速度偏心率越小,说明截面上的速度分布越均匀,流场越稳定,整流器的整流效果越好。根据Origin软件绘制出速度偏心率随流动发展的曲线图可以直观地看到不同结构整流器下游截面速度偏心率的大小,以此来判定不同结构整流器整流效果的好坏。The specific method for judging the rectification effect is: the smaller the velocity eccentricity of the cross section, the more uniform the velocity distribution on the cross section, the more stable the flow field, and the better the rectification effect of the rectifier. Drawing the curve of velocity eccentricity with flow development based on the Origin software can visually see the velocity eccentricity of the downstream cross-section of rectifiers with different structures, so as to determine the rectification effect of rectifiers with different structures.

实施例二;Embodiment 2;

如图1-7所示,根据本发明实施例的适用于流程阀门的整流效果判定方法,包括以下步骤:As shown in Figures 1-7, the rectification effect determination method suitable for process valves according to the embodiment of the present invention includes the following steps:

步骤一:建立流道模型;选择球阀为扰流元件,对实际工况下的阀门和整流器利用Solidworks软件进行三维建模,将整流器I和整流器II模型装配在球阀出口处,整流器I和整流器II模型如图2和图3所示,设定球阀相对开度为30%,对模型流体域进行抽取,获得流道模型,流道模型特征在于:所述流道模型阀门上游为五倍管径,整流器下游为十二倍管径,流道模型如图4所示;Step 1: Establish the flow channel model; select the ball valve as the spoiler element, use Solidworks software to conduct three-dimensional modeling of the valve and rectifier under actual working conditions, and assemble the rectifier I and rectifier II models at the outlet of the ball valve. Rectifier I and rectifier II The model is shown in Figures 2 and 3. The relative opening of the ball valve is set to 30%. The model fluid domain is extracted to obtain a flow channel model. The characteristics of the flow channel model are: the upstream of the valve in the flow channel model is five times the pipe diameter. , the downstream of the rectifier is twelve times the pipe diameter, and the flow channel model is shown in Figure 4;

步骤二:划分网格;进行网格无关性验证,选取符合要求的网格数,对流道模型进行网格划分;Step 2: Divide the mesh; verify the grid independence, select the number of meshes that meet the requirements, and mesh the flow channel model;

所述的划分网格具体实施包括以下步骤:The specific implementation of meshing includes the following steps:

在Solidworks软件里面,将流道模型的每一个部分分别导出;In the Solidworks software, export each part of the flow channel model separately;

利用ANSYS ICEM CFD软件,对流道模型的每一个部分分别划分结构性网格,并对阀芯和整流器的部分进行网格加密,最后利用ANSYS ICEM CFD软件的Merge功能,将每个部分的网格合并成一个整体,流道网格模型如图5所示。Use ANSYS ICEM CFD software to divide each part of the flow channel model into structural meshes, and refine the mesh for the valve core and rectifier parts. Finally, use the Merge function of ANSYS ICEM CFD software to merge the meshes of each part. Merged into a whole, the flow channel grid model is shown in Figure 5.

所述的网格无关性验证具体实施包括以下步骤:The specific implementation of the grid independence verification includes the following steps:

针对阀门全开、安装整流器I时的流道模型,进行网格划分,并划分成多套具有不同网格数量的流道网格模型;For the flow channel model when the valve is fully open and the rectifier I is installed, the mesh is divided into multiple sets of flow channel mesh models with different mesh numbers;

比较不同网格数下的流量系数和阻力系数,当流量系数和阻力系数随网格数量的增加而基本保持不变时,选择流量系数和阻力系数基本保持不变的最少网格数,如图6所示,选择350万网格。Compare the flow coefficient and resistance coefficient under different grid numbers. When the flow coefficient and resistance coefficient remain basically unchanged as the number of grids increases, select the minimum number of grids where the flow coefficient and resistance coefficient remain basically unchanged, as shown in the figure As shown in 6, select 3.5 million grid.

步骤三:进行数值计算;将流道网格导入ANASYS FLUENT软件,设置相关参数,并进行数值计算;Step 3: Perform numerical calculations; import the flow channel grid into ANASYS FLUENT software, set relevant parameters, and perform numerical calculations;

所述的设置相关参数具体实施包括以下步骤:The specific implementation of setting relevant parameters includes the following steps:

①在General模块中定义求解器,求解器类型采用基于压力的求解器,时间类型设置为稳态计算方式;① Define the solver in the General module. The solver type is a pressure-based solver, and the time type is set to the steady-state calculation method;

②在Models模块中选择计算模型,选择Standard k-ε模型进行计算;②Select the calculation model in the Models module and select the Standard k-ε model for calculation;

③在Materials模块中设置流体介质,根据实际工况,选择流体介质为水;③Set the fluid medium in the Materials module, and select water as the fluid medium according to the actual working conditions;

④在Boundary Conditions模块中设置边界条件,入口选择velocity-inlet,出口选择pressure-outlet,壁面选择wall,交界面选择interface;根据实际工作条件设置边界条件,入口条件设定为1m/s,安装整流器I和整流器II的出口条件分别设定为17876.8Pa和17880.9Pa;④Set boundary conditions in the Boundary Conditions module, select velocity-inlet for the inlet, pressure-outlet for the outlet, wall for the wall, and interface for the interface; set the boundary conditions according to the actual working conditions, set the inlet condition to 1m/s, and install the rectifier The outlet conditions of I and rectifier II are set to 17876.8Pa and 17880.9Pa respectively;

⑤在Mesh Interfaces模块中定义interface;⑤Define interface in the Mesh Interfaces module;

⑥在Residual模块中设置收敛残差;⑥Set the convergence residual in the Residual module;

⑦在Solution Initialization模块卡中点击Hybrid Initialization进行初始化;⑦Click Hybrid Initialization in the Solution Initialization module card to initialize;

⑧在Run Calculation模块中设置迭代步数后点击Calculation模块开始进行数值计算。⑧After setting the number of iteration steps in the Run Calculation module, click the Calculation module to start numerical calculations.

步骤四:计算下游截面的速度偏心率η;提取整流器下游管道横截面上的最大速度,根据公式,计算下游每个横截面的速度偏心率η,计算公式如下: Step 4: Calculate the velocity eccentricity η of the downstream section; extract the maximum velocity on the cross-section of the pipeline downstream of the rectifier, and calculate the velocity eccentricity eta of each downstream cross-section according to the formula. The calculation formula is as follows:

其中,Vmax表示下游管道横截面上的最大速度,单位为m/s;V0表示管道内的平均速度,单位为m/s。Among them, Vmax represents the maximum velocity on the cross section of the downstream pipeline, in m/s; V0 represents the average velocity in the pipeline, in m/s.

提取整流器下游管道横截面上的最大速度具体实施包括以下步骤:The specific implementation of extracting the maximum velocity on the cross-section of the pipeline downstream of the rectifier includes the following steps:

①在Surface模块中选择Iso-Surface,在整流器下游每间隔1D创建一个横截面,共创建12个横截面;① Select Iso-Surface in the Surface module and create a cross section every 1D downstream of the rectifier, creating a total of 12 cross sections;

②在Reports模块中的Surface Integrals模块,Report Type选择FacetMaximum,Field Variable选择Velocity Magnitude,Surfaces选择在整流器下游创建的横截面,即可得到横截面上的最大速度,整流器I和整流器II下游横截面上的最大速度如下表:② In the Surface Integrals module in the Reports module, select FacetMaximum for Report Type, select Velocity Magnitude for Field Variable, and select the cross section created downstream of the rectifier for Surfaces to get the maximum speed on the cross section, on the cross sections downstream of rectifier I and rectifier II. The maximum speed is as follows:

速度偏心率η计算结果如下:The calculation results of velocity eccentricity eta are as follows:

步骤五:判断整流效果;根据计算出来的速度偏心率,利用Origin软件绘制出速度偏心率随流动发展的曲线图,对比分析不同结构整流器的整流效果。Step 5: Determine the rectification effect; based on the calculated speed eccentricity, use Origin software to draw a curve chart of the speed eccentricity developing with the flow, and compare and analyze the rectification effects of rectifiers with different structures.

所述的判断整流效果的具体方法是:横截面的速度偏心率越小,说明截面上的速度分布越均匀,流场越稳定,整流器的整流效果越好。,如图7所示的速度偏心率随流动发展的曲线图可以直观地看到整流器I下游截面速度偏心率均大于整流器II下游截面速度偏心率,整流器II下游的流场更加均匀,因此可以判定,整流器II的整流效果要比整流器I好。The specific method for judging the rectification effect is: the smaller the velocity eccentricity of the cross section, the more uniform the velocity distribution on the cross section, the more stable the flow field, and the better the rectification effect of the rectifier. , as shown in Figure 7, the velocity eccentricity curve develops with the flow. It can be intuitively seen that the velocity eccentricity of the downstream section of rectifier I is greater than the velocity eccentricity of the downstream section of rectifier II. The flow field downstream of rectifier II is more uniform, so it can be determined , the rectification effect of rectifier II is better than that of rectifier I.

综上所述,借助于本发明的上述技术方案,通过数值模拟的方法,对阀后的整流情况进行了反映,利用整流器下游横截面的速度偏心率判定不同结构整流器整流效果的好坏,不需要进行复杂的实验来判断,可以大大减少人力物力,在工程上具有实际意义。To sum up, with the help of the above technical solution of the present invention, the rectification situation behind the valve is reflected through the method of numerical simulation, and the speed eccentricity of the downstream cross section of the rectifier is used to determine the rectification effect of the rectifiers with different structures. Complex experiments are required to judge, which can greatly reduce manpower and material resources, and has practical significance in engineering.

以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The above descriptions are only preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention shall be included in the present invention. within the scope of protection.

Claims (4)

1. The rectifying effect judging method suitable for the flow valve is characterized by comprising the following steps of:
s101, establishing a flow channel model; three-dimensional modeling is carried out on a valve and a rectifier under actual working conditions by utilizing Solidworks software, a rectifier model is assembled at an outlet of the valve, a model fluid domain is extracted, a flow channel model is obtained, and the flow channel model is characterized in that: the upstream of the flow channel model valve is five-pipe diameter, and the downstream of the rectifier is twelve-pipe diameter;
s103, dividing grids; performing grid independence verification, selecting grid numbers meeting requirements, and performing grid division on the flow model;
s105, performing numerical calculation; importing the runner grid into ANASYS FLUENT software, setting related parameters, and calculating a digital value;
s107, calculating the velocity eccentricity eta of the downstream section; the maximum speed on the cross section of the pipeline downstream of the rectifier is extracted, and the speed eccentricity eta of each cross section downstream is calculated according to the formula, wherein the calculation formula is as follows:the method comprises the steps of carrying out a first treatment on the surface of the Wherein Vmax represents the maximum velocity in m/s over the downstream pipe cross section; v0 represents the average velocity in the pipeline in m/s;
s109, judging the rectifying effect; drawing a graph of the speed eccentricity along with the flow development by using Origin software according to the calculated speed eccentricity, and comparing and analyzing the rectifying effects of rectifiers with different structures;
the implementation of extracting the maximum speed of the cross section of the pipeline at the downstream of the rectifier in the step S107 comprises the following steps:
s1071, selecting an Iso-Surface in a Surface module, and creating a cross section at intervals of 1D at the downstream of the rectifier to create 12 cross sections in total;
s1072, a Surface Integrals module in the Reports module, wherein the Report Type selects Facet Maximum, the Field Variable selects Velocity Magnitude, and the Surfaces select the cross section created at the downstream of the rectifier in the step S1071, so that the Maximum speed on the cross section can be obtained;
the specific method for judging the rectifying effect is as follows: the smaller the speed eccentricity of the cross section is, the more uniform the speed distribution on the cross section is, the more stable the flow field is, and the better the rectifying effect of the rectifier is; and drawing a graph of the speed eccentricity along with the flow development according to Origin software, and intuitively seeing the speed eccentricity of the downstream section of the rectifiers with different structures, so as to judge the rectifying effect of the rectifiers with different structures.
2. The rectifying effect determining method for a flow valve according to claim 1, wherein the meshing implementation comprises the steps of:
s1031, respectively exporting each part of the flow channel model in Solidworks software;
s1032, dividing each part of the flow model into structural grids by using ANSYS ICEM CFD software, encrypting the grids of the valve core and the rectifier part, and finally merging the grids of each part into a whole by using the Merge function of ANSYS ICEM CFD software.
3. The rectifying effect determining method for a flow valve according to claim 1, wherein the grid independence verification implementation comprises the following steps:
s1033, dividing grids aiming at the runner model when the valve is fully opened and a certain rectifier is installed, and dividing the runner model into a plurality of runner grid models with different grid numbers;
s1034, comparing the flow coefficient and the resistance coefficient under different grid numbers, and selecting the minimum grid number with the flow coefficient and the resistance coefficient unchanged when the flow coefficient and the resistance coefficient keep unchanged along with the increase of the grid number.
4. The rectifying effect determining method for a flow valve according to claim 1, wherein the setting of the relevant parameters in the step S105 comprises the following steps:
s1051, defining a solver in a General module, wherein the solver adopts a solver based on pressure, and the time type is set as a steady-state calculation mode;
s1052, selecting a calculation model in a Models module, and selecting a Standard k-epsilon model for calculation;
s1053, setting a fluid medium in the Materials module, and selecting the fluid medium as water or gas according to actual working conditions;
s1054, setting boundary conditions in a Boundary Conditions module, selecting a property-inlet as an inlet, selecting a pressure-outlet as an outlet, selecting wall surfaces, and selecting interfaces as interfaces; when the actual working conditions are clear, setting boundary conditions according to the actual working conditions; when the actual working condition is ambiguous, setting reasonable boundary conditions by oneself;
s1055, defining Interfaces in the Mesh Interfaces module;
s1056, setting convergence Residual in a Residual module;
s1057, clicking Hybrid Initialization in the Initialization module for Initialization;
s1058, after the number of iteration steps is set in the Run Calculation module, clicking the Calculation module to start numerical Calculation.
CN202110912459.9A 2021-08-10 2021-08-10 Rectifying effect judging method suitable for flow valve Active CN113591414B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110912459.9A CN113591414B (en) 2021-08-10 2021-08-10 Rectifying effect judging method suitable for flow valve

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110912459.9A CN113591414B (en) 2021-08-10 2021-08-10 Rectifying effect judging method suitable for flow valve

Publications (2)

Publication Number Publication Date
CN113591414A CN113591414A (en) 2021-11-02
CN113591414B true CN113591414B (en) 2023-10-13

Family

ID=78256624

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110912459.9A Active CN113591414B (en) 2021-08-10 2021-08-10 Rectifying effect judging method suitable for flow valve

Country Status (1)

Country Link
CN (1) CN113591414B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101097849A (en) * 2004-07-09 2008-01-02 积水化学工业株式会社 Method for processing outer periphery of substrate and apparatus thereof
JP2012202249A (en) * 2011-03-24 2012-10-22 Toshiba Corp Steam valve
WO2013148021A1 (en) * 2012-03-28 2013-10-03 Exxonmobil Upstream Research Company Method for mutiphase flow upscaling
CN109325265A (en) * 2018-08-29 2019-02-12 浙江理工大学 A design and construction method of pipeline valve noise reduction model
CN111241757A (en) * 2020-01-10 2020-06-05 中核第四研究设计工程有限公司 Dam break three-dimensional numerical simulation method for uranium tailing pond based on computational fluid mechanics

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101097849A (en) * 2004-07-09 2008-01-02 积水化学工业株式会社 Method for processing outer periphery of substrate and apparatus thereof
JP2012202249A (en) * 2011-03-24 2012-10-22 Toshiba Corp Steam valve
WO2013148021A1 (en) * 2012-03-28 2013-10-03 Exxonmobil Upstream Research Company Method for mutiphase flow upscaling
CN109325265A (en) * 2018-08-29 2019-02-12 浙江理工大学 A design and construction method of pipeline valve noise reduction model
CN111241757A (en) * 2020-01-10 2020-06-05 中核第四研究设计工程有限公司 Dam break three-dimensional numerical simulation method for uranium tailing pond based on computational fluid mechanics

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
屈天鹏 ; 郭瑞琪 ; 郝月莹 ; 张文帝 ; 袁子凯 ; 王德永 ; .连铸钢包浇注过程中旋涡的形成机理.钢铁研究学报.2016,(第08期),全文. *
崔宝玲 ; 尚照辉 ; 石柯 ; 陈德胜 ; 林哲 ; .基于CFD的蝶板结构改进设计及数值分析.排灌机械工程学报.2013,(第06期),全文. *

Also Published As

Publication number Publication date
CN113591414A (en) 2021-11-02

Similar Documents

Publication Publication Date Title
CN105677964B (en) A kind of valve flow coefficient computational methods based on CFD emulation and mesh adaption
CN103729505B (en) CFD (computational fluid dynamics) based method for computing equivalent length of valve
CN111767663A (en) A convenient calculation method of valve flow coefficient based on CFD simulation
Nikpour et al. Experimental and numerical simulation of water hammer
CN111400941A (en) A Numerical Prediction Method of Internal Backflow and Backflow Vortex Cavitation in Vane Pumps
CN111460699A (en) Design method of microtexture with drag reduction function on flat wall surface
CN112597709B (en) A CFD Calculation Method of Turbine Cascade Flow Channel Considering Transition Disturbance Factors
CN113553737B (en) Valve flow prediction method based on valve pressure difference
Cao et al. 3D CFD simulation and analysis of transient flow in a water pipeline
CN113591414B (en) Rectifying effect judging method suitable for flow valve
Wang et al. Inter-stage performance and energy characteristics analysis of electric submersible pump based on entropy production theory
Jiang et al. Large eddy simulation on the vortex evolution in a squirrel-cage fan based on a slice computational model
CN102750414B (en) A kind of method verifying centrifugal pump mesh quality and computational accuracy relation
CN113536640B (en) Optimization design method for internal flow channel structure of air distributor based on orthogonal test
CN117807893A (en) Multi-objective optimization design method for impeller of high-speed centrifugal pump
Adamkowski et al. The comparative analysis of the current-meter method and the pressure-time method used for discharge measurements in the Kaplan turbine penstocks
CN110390118A (en) A Partially Cut Pipeline Network Hydraulic Model Construction Method
CN114036873A (en) Calculation method of super-hydrophobic surface flat laminar boundary layer with leading edge area correction considered
Yang Flow patterns causing saddle instability in the performance curve of a centrifugal pump with vaned diffuser
CN115659689B (en) Pump device inflow and performance numerical value iterative calculation method reflecting actual boundary conditions
CN113673133B (en) An Analytical Method for Energy Conversion of Centrifugal Pumps Based on Energy Transport Theory
CN113536711A (en) Olive-shaped flow meter structure optimization design method based on Taguchi method
Sun et al. Numerical simulation of turbine flowmeter's three-dimensional flow fields
Zheng et al. Study on internal flow field simulation accuracy of centrifugal impellers based on different meshing types
CN113536705A (en) Evaluation method of flow field distribution based on computational fluid dynamics

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