具体实施方式

为进一步说明各实施例，本发明提供有附图，这些附图为本发明揭露内容的一部分，其主要用以说明实施例，并可配合说明书的相关描述来解释实施例的运作原理，配合参考这些内容，本领域普通技术人员应能理解其他可能的实施方式以及本发明的优点，图中的组件并未按比例绘制，而类似的组件符号通常用来表示类似的组件。

根据本发明的实施例，提供了一种适用于流程阀门的整流效果判定方法。

实施例一；

如图1-7所示，根据本发明实施例的适用于流程阀门的整流效果判定方法，包括以下步骤：

步骤S101、建立流道模型；对实际工况下的阀门和整流器利用Solidworks软件进行三维建模，将整流器模型装配在阀门出口处，对模型流体域进行抽取，获得流道模型，流道模型特征在于：所述流道模型阀门上游为五倍管径，整流器下游为十二倍管径；

步骤S103、划分网格；进行网格无关性验证，选取符合要求的网格数，对流道模型进行网格划分；

步骤S105、进行数值计算；将流道网格导入ANASYS FLUENT软件，设置相关参数，并进行数值计算；

步骤S107、计算下游截面的速度偏心率η；提取整流器下游管道横截面上的最大速度，根据公式，计算下游每个横截面的速度偏心率η，计算公式如下：

其中，Vmax表示下游管道横截面上的最大速度，单位为m/s；V0表示管道内的平均速度，单位为m/s。

步骤S109，判断整流效果；根据计算出来的速度偏心率，利用Origin软件绘制出速度偏心率随流动发展的曲线图，对比分析不同结构整流器的整流效果。

所述划分网格具体实施包括以下步骤：

步骤S1031、在Solidworks软件里面，将流道模型的每一个部分分别导出；

步骤S1032、利用ANSYS ICEM CFD软件，对流道模型的每一个部分分别划分结构性网格，并对阀芯和整流器的部分进行网格加密，最后利用ANSYS ICEM CFD软件的Merge功能，将每个部分的网格合并成一个整体。

所述网格无关性验证具体实施包括以下步骤：

步骤S1033、针对阀门全开、安装某一种整流器时的流道模型，进行网格划分，并划分成多套具有不同网格数量的流道网格模型；

步骤S1034、比较不同网格数下的流量系数和阻力系数，当流量系数和阻力系数随网格数量的增加而基本保持不变时，选择流量系数和阻力系数基本保持不变的最少网格数。

所述步骤S105中设置相关参数具体实施包括以下步骤：

步骤S1051、在General模块中定义求解器，求解器类型采用基于压力的求解器，时间类型设置为稳态计算方式；

步骤S1052、在Models模块中选择计算模型，选择Standard k-ε模型进行计算；

步骤S1053、在Materials模块中设置流体介质，根据实际工况，选择流体介质为水或者气；

步骤S1054、在Boundary Conditions模块中设置边界条件，入口选择velocity-inlet，出口选择pressure-outlet，壁面选择wall，交界面选择interface；当实际工作条件明确时，根据实际工作条件设置边界条件；当实际工作条件不明确时，自行设定合理的边界条件；

步骤S1055、在Mesh Interfaces模块中定义interface；

步骤S1056、在Residual模块中设置收敛残差；

步骤S1057、在Initialization模块中点击Hybrid Initialization进行初始化；

步骤S1058、在Run Calculation模块中设置迭代步数后点击Calculation模块开始进行数值计算。

所述的步骤S107中提取整流器下游管道横截面上的最大速度具体实施包括以下步骤：

步骤S1071、在Surface模块中选择Iso-Surface，在整流器下游每间隔1D创建一个横截面，共创建12个横截面；

步骤S1072、在Reports模块中的Surface Integrals模块，Report Type选择FacetMaximum，Field Variable选择Velocity Magnitude，Surfaces选择S1071步骤中在整流器下游创建的横截面，即可得到横截面上的最大速度。

所述的判断整流效果的具体方法是：横截面的速度偏心率越小，说明截面上的速度分布越均匀，流场越稳定，整流器的整流效果越好。根据Origin软件绘制出速度偏心率随流动发展的曲线图可以直观地看到不同结构整流器下游截面速度偏心率的大小，以此来判定不同结构整流器整流效果的好坏。

实施例二；

如图1-7所示，根据本发明实施例的适用于流程阀门的整流效果判定方法，包括以下步骤：

步骤一：建立流道模型；选择球阀为扰流元件，对实际工况下的阀门和整流器利用Solidworks软件进行三维建模，将整流器I和整流器II模型装配在球阀出口处，整流器I和整流器II模型如图2和图3所示，设定球阀相对开度为30％，对模型流体域进行抽取，获得流道模型，流道模型特征在于：所述流道模型阀门上游为五倍管径，整流器下游为十二倍管径，流道模型如图4所示；

步骤二：划分网格；进行网格无关性验证，选取符合要求的网格数，对流道模型进行网格划分；

所述的划分网格具体实施包括以下步骤：

在Solidworks软件里面，将流道模型的每一个部分分别导出；

利用ANSYS ICEM CFD软件，对流道模型的每一个部分分别划分结构性网格，并对阀芯和整流器的部分进行网格加密，最后利用ANSYS ICEM CFD软件的Merge功能，将每个部分的网格合并成一个整体，流道网格模型如图5所示。

所述的网格无关性验证具体实施包括以下步骤：

针对阀门全开、安装整流器I时的流道模型，进行网格划分，并划分成多套具有不同网格数量的流道网格模型；

比较不同网格数下的流量系数和阻力系数，当流量系数和阻力系数随网格数量的增加而基本保持不变时，选择流量系数和阻力系数基本保持不变的最少网格数，如图6所示，选择350万网格。

步骤三：进行数值计算；将流道网格导入ANASYS FLUENT软件，设置相关参数，并进行数值计算；

所述的设置相关参数具体实施包括以下步骤：

①在General模块中定义求解器，求解器类型采用基于压力的求解器，时间类型设置为稳态计算方式；

②在Models模块中选择计算模型，选择Standard k-ε模型进行计算；

③在Materials模块中设置流体介质，根据实际工况，选择流体介质为水；

④在Boundary Conditions模块中设置边界条件，入口选择velocity-inlet，出口选择pressure-outlet，壁面选择wall，交界面选择interface；根据实际工作条件设置边界条件，入口条件设定为1m/s，安装整流器I和整流器II的出口条件分别设定为17876.8Pa和17880.9Pa；

⑤在Mesh Interfaces模块中定义interface；

⑥在Residual模块中设置收敛残差；

⑦在Solution Initialization模块卡中点击Hybrid Initialization进行初始化；

⑧在Run Calculation模块中设置迭代步数后点击Calculation模块开始进行数值计算。

步骤四：计算下游截面的速度偏心率η；提取整流器下游管道横截面上的最大速度，根据公式，计算下游每个横截面的速度偏心率η，计算公式如下：

其中，Vmax表示下游管道横截面上的最大速度，单位为m/s；V0表示管道内的平均速度，单位为m/s。

提取整流器下游管道横截面上的最大速度具体实施包括以下步骤：

①在Surface模块中选择Iso-Surface，在整流器下游每间隔1D创建一个横截面，共创建12个横截面；

②在Reports模块中的Surface Integrals模块，Report Type选择FacetMaximum，Field Variable选择Velocity Magnitude，Surfaces选择在整流器下游创建的横截面，即可得到横截面上的最大速度，整流器I和整流器II下游横截面上的最大速度如下表：

速度偏心率η计算结果如下：

步骤五：判断整流效果；根据计算出来的速度偏心率，利用Origin软件绘制出速度偏心率随流动发展的曲线图，对比分析不同结构整流器的整流效果。

所述的判断整流效果的具体方法是：横截面的速度偏心率越小，说明截面上的速度分布越均匀，流场越稳定，整流器的整流效果越好。，如图7所示的速度偏心率随流动发展的曲线图可以直观地看到整流器I下游截面速度偏心率均大于整流器II下游截面速度偏心率，整流器II下游的流场更加均匀，因此可以判定，整流器II的整流效果要比整流器I好。

综上所述，借助于本发明的上述技术方案，通过数值模拟的方法，对阀后的整流情况进行了反映，利用整流器下游横截面的速度偏心率判定不同结构整流器整流效果的好坏，不需要进行复杂的实验来判断，可以大大减少人力物力，在工程上具有实际意义。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。