CN105912745A

CN105912745A - 袋式除尘器气流组织多参数优化方法

Info

Publication number: CN105912745A
Application number: CN201610170844.XA
Authority: CN
Inventors: 解彬; 李义新; 徐兵
Original assignee: ANHUI WEIDA ENVIRONMENTAL PROTECTION TECHNOLOGY Co Ltd
Current assignee: ANHUI WEIDA ENVIRONMENTAL PROTECTION TECHNOLOGY Co Ltd
Priority date: 2016-03-24
Filing date: 2016-03-24
Publication date: 2016-08-31

Abstract

本发明涉及到一种袋式除尘器气流组织多参数优化方法。其首先确定袋式除尘器气流组织主控设计参数及除尘器基本几何参数，利用响应面法进行多参数工况优化设计，建立各工况除尘器几何模型并确定计算区域，运用网格划分工具对计算区域划分拓扑网格，然后采用计算流体动力学原理进行离散化方程迭代求解计算，基于响应面法对求解计算结果进行多参数优化设计，最终得到最优解集。该方法不仅可以得到除尘器内任意断面处的速度、压力分布，而且还可以基于响应面法也就是多参数优化设计方法，对袋式除尘器气流组织进行多参数优化设计，该方法适用于各种类型袋式除尘器的气流组织优化设计，并且结果可靠准确。

Description

袋式除尘器气流组织多参数优化方法

技术领域

本发明涉及袋式除尘器，属于袋式除尘器气流组织优化设计领域，具体涉及到一种袋式除尘器气流组织多参数优化方法。

背景技术

随着计算机技术的飞速发展，计算流体动力学已被广泛应用于各行各业的除尘设备的气流组织数值模拟。袋式除尘器作为一种可靠性强的气固分离设备已被广泛应用于各行各业的粉尘治理，其内部流过的气流为稀相气固两相流（可简化为单向流进行气流组织数值模拟），并且袋式除尘器内部气流组织对其滤袋使用寿命、过滤效率和压力损失等有至关重要的影响，因此，需要对袋式除尘器内部气流组织进行优化设计。目前，针对袋式除尘器气流组织数值模拟多局限于单一工况单变量参数的数值模拟，只能保证特定工况下除尘器的性能，但不能确保多参数（除尘器本体结构参数、滤料特性参数及运行参数）共同作用下的除尘器最佳性能。

发明内容

本发明是为了解决传统袋式除尘器内气流组织设计方法的不足，公开了一种袋式除尘器气流组织多参数优化方法，该方法不仅可以得到除尘器内任意断面处的速度、压力分布，而且还可以基于响应面法对袋式除尘器气流组织进行多参数优化设计。

为了实现以上目的，本发明袋式除尘器气流组织多参数优化方法的技术方案为：

包括以下步骤：

1）根据烟气量、过滤风速及现场空地尺寸确定除尘器基本几何尺寸，确定主控设计参数及袋式除尘器基本几何参数；

2）以定量衡量除尘器气流组织分布的综合流量不均幅值作为目标函数，以步骤1）中所确定的主控设计参数作为自变量，利用响应曲面法进行多参数工况优化设计；

3）利用建模单元建立袋式除尘器几何模型并确定求解计算区域；

4）将步骤3）中几何模型文件导入计算网格划分单元进行离散化拓扑网格划分，对除尘器实际物理边界进行简化，设定计算边界条件，保存计算结果并导出网格文件；

5）将网格文件导入CFD求解单元，选取计算求解模型，对方程进行离散化处理，设定求解计算初始参数，选取数值计算迭代方法进行迭代求解计算，直至计算收敛，保存计算结果；

6）对计算结果进行相关后处理，得到任意截面处的速度、压力分布云图及各工况下的综合流量不均幅值，基于响应面法对计算结果进行多参数优化设计，得到最优解集。

优选的，步骤1）中所述的主控设计参数包括安全高度，进口风速，滤料厚度；所述的袋式除尘器基本几何参数包括灰斗、中箱体、上箱体具体尺寸。

优选的，步骤2）中所述的响应曲面法以滤料厚度、安全高度、进口速度为响应因素，响应值为数值模拟所获得袋式除尘器气流组织分布的综合流量不均幅值。

进一步优选的，所述的步骤3）采用Solidworks建模单元建立步骤2）中各工况的三维几何模型；再打开所述的步骤4）中的拓扑网格划分单元，读取所述的三维几何模型，对读取的模型进行拓扑计算，依据实际几何模型边界对各部分进行划分并命名，设定全局及各网格划分尺寸，选择八叉树算法对计算模型进行拓扑网格划分，并对网格进行光滑处理，再对除尘器实际物理边界进行简化，分别设定袋式除尘器进口为速度进口边界条件、出口为压力出口边界条件、花板及壳体壁面为固壁边界条件、滤袋为跳跃多孔介质边界条件，导出网格文件再将网格文件导入步骤5）中的CFD求解单元Fluent中，采用标准k-ε双方程模型，对方程进行离散化处理，设定求解计算初始参数，选取数值计算迭代方法进行迭代求解计算，直至计算收敛，保存计算结果。

本发明的技术效果在于：本发明基于CFD多参数共同作用，该方法不仅可以得到除尘器内任意断面处的速度、压力分布，而且还可以基于响应面法也就是多参数优化设计方法，对袋式除尘器气流组织进行多参数优化设计，该方法适用于各种类型袋式除尘器的气流组织优化设计，并且结果可靠准确。其首先确定袋式除尘器气流组织主控设计参数及除尘器基本几何参数，利用响应面法进行多参数工况优化设计，建立各工况除尘器几何模型并确定计算区域，运用网格划分工具对计算区域划分拓扑网格，然后采用计算流体动力学（CFD）原理进行离散化方程迭代求解计算，基于响应面法对求解计算结果进行多参数优化设计，最终得到最优解集（气流组织最佳时，各参数实际对应取值）。

附图说明

图1为本发明实施例的流程图。

具体实施方式

为了实现在多参数共同作用下袋式除尘器气流组织优化设计，该方法具体实施方案如下：

（1）根据烟气量、过滤风速及现场空地尺寸确定除尘器基本几何尺寸，确定主控设计参数及袋式除尘器基本几何参数；

（2）以定量衡量除尘器气流组织分布的综合流量不均幅值作为目标函数，以上述所确定的主控设计参数作为自变量，利用响应面法进行多参数工况优化设计；

响应曲面法(Response surface methodology, RSM)是一种数学和统计学相结合的数据处理方法,用来对目标函数响应值受多个变量影响的问题进行建模和分析,并对响应结果进行优化。选取二次响应曲面方程，考虑所有的一次项，二次项和两两交叉项，响应曲面方程可表示为：

(1)

其中，Y为目标函数或称响应；X_i为自变量；β_i，β_ii，β_ij代表一次、二次、交互作用项的回归系数；k为影响因素的数量；e为误差。上式中的回归系数可由最小二乘法拟合得到。为方便起见，将所有变量按下式规格化：

(2)

其中X_iH和X_iL分别为变量的最大值和最小值，为变量的平均值。

（3）利用建模单元（如： Solidworks、ProE）建立袋式除尘器几何模型并确定求解计算区域；

（4）将上述几何模型文件导入计算网格划分单元（如：ICEM-cfd、Gambit）进行离散化拓扑网格划分，对除尘器实际物理边界进行合理简化，设定计算边界条件，保存计算结果并导出网格文件；

（5）将网格文件导入CFD求解单元（如：Fluent、OpenFoam），选取计算求解模型，对方程进行离散化处理，设定求解计算初始参数，选取数值计算迭代方法进行迭代求解计算，直至计算收敛，保存计算结果；

其中计算求解所涉及的三大守恒方程分别为：

a）连续性方程

由质量守恒定律可知，单位时间内流进、流出控制体的流体质量差等于控制体内流体面密度发生变化所引起的质量增量。由此可推导出连续性方程，其表达式为：

(3)

其中，ρ为流体的密度，kg/m³； u为流体速度，m/s。

b）动量守恒方程

根据动量定理，单位时间内，控制体的动量的变化率等于作用在该体积上的质量力和表面力之和，由此得出动量守恒方程如下：

(4)

其中P, τ分别为单位体积上流体所受压力及粘性应力，ρf表示单位体积上的质量力。

c）能量守恒方程

由能量守恒定律可知，单位时间内控制体内流体的动能和内能的变化率等于作用于控制体上的质量力和面力所作的功加上单位时间内给予该体积的热流量，其表达式为：

其中e，u²/2分别是单位质量流体的内能和动能；k为热传导系数，W/m²·K； T是温度；q为由于辐射或其它原因在单位时间内传入单位质量的热量分布函数。

（6）对计算结果进行相关后处理，得到任意截面处的速度、压力分布云图及各工况下的综合流量不均幅值，基于响应面法对计算结果进行多参数优化设计，得到最优解集。

为了更加形象至关地阐述本发明的方案及优点，下述实施例将结合附图对本发明的实施步骤进行详细深入的说明，具体实施步骤流程如图1所示（图示中响应面法即所述的响应曲面法）。

实施例：

以一种气箱脉冲喷吹袋式除尘器LFGM-32-2为例，基于响应面法对其气流分布进行多参数优化设计，详细实施步骤如下所示：

（1）根据烟气量、过滤风速及现场空地尺寸确定除尘器基本几何尺寸，确定主控设计参数（安全高度、进口风速、滤料厚度）及袋式除尘器基本几何参数（灰斗、中箱体、上箱体等具体尺寸）；

（2）以定量衡量除尘器气流组织分布的综合流量不均幅值作为目标函数，以上述所确定的主控设计参数作为自变量，利用响应面法进行多参数工况优化设计。以滤料厚度t(X1)、安全高度h (X2)、进口速度v (X3)为响应因素，响应值为数值模拟所获得袋式除尘器气流组织分布的综合流量不均幅值ΔK (Y)。三种响应因素的水平值如表1所示。采用三因素三水平Box-Behnken法进行设计，本文所需数值模拟次数为15次，其中析因部分实验次数12次，保证均一精密性的中心点重复次数为3。具体设计工况见表2。

表1 RSM模型中各因素水平值设置

袋式除尘器气流分布综合流量不均幅值ΔK为：

(6)

其中k_qi为气流通过单个滤袋出口的流量分配系数；n为滤袋总数量。

表2 响应面试验设计

工况	X1 (t, mm)	X2 (h, mm)	X3 (v, m/s)
				1	2.5	400	10.19
2	2.5	300	9.34
				3	2.0	400	9.34
4	2.0	300	8.49
				5	2.0	400	9.34
6	1.5	400	8.49
				7	2.0	300	10.19
8	2.0	500	10.19
				9	2.5	500	9.34
10	1.5	400	10.19
				11	2.0	400	9.34
12	2.0	500	8.49
				13	1.5	300	9.34
14	2.5	400	8.49
				15	1.5	500	9.34

（3）基于气箱脉冲喷吹袋式除尘器LFGM-32-2的基本参数，确定求解计算区域，采用Solidworks建立（2）中各工况的三维几何模型，并将所建立的三维模型保存。

（4）打开拓扑网格划分单元ICEM-cfd，读取（3）中所建立模型的三维模型，对读取的模型进行拓扑计算，依据实际几何模型边界对各部分进行part划分并命名。设定全局及各part网格划分尺寸，选择八叉树算法对计算模型进行拓扑网格划分，并对网格进行光滑等处理。然后，对除尘器实际物理边界进行合理简化，分别设定袋式除尘器进口为速度进口边界条件、出口为压力出口边界条件、花板及壳体壁面为固壁边界条件、滤袋为跳跃多孔介质边界条件。最终，将划分好的网格文件保存；

多孔介质模型方程为：

(7)

其中，Si为动量耗散量；μ为气体粘度；ρ为气体密度；α为透气系数；C为内部阻力系数；v_i为垂直于介质表面的速度。

（5）将网格文件导入求解单元Fluent中，采用标准k-ε双方程模型，对方程进行离散化处理，设定求解计算初始参数，选取数值计算迭代方法进行迭代求解计算，直至计算收敛，保存计算结果。为了确保计算的准确性，采用二阶迎风格式Simple算法对离散化动量方程进行压力速度耦合求解。并将连续性方程及动量方程在直角坐标系x，y及z方向上的收敛残差设定在10-6以内。

（6）对计算结果进行相关后处理，得到任意截面处的速度、压力分布云图及各工况下的综合流量不均幅值。

（7）采用响应面法优化算法对（6）中计算结果进行多参数优化设计，最终得到最优解集。

Claims

1.一种袋式除尘器气流组织多参数优化方法，其特征在于：包括以下步骤：

2.根据权利要求1所述的袋式除尘器气流组织多参数优化方法，其特征在于：步骤1）中所述的主控设计参数包括安全高度，进口风速，滤料厚度；所述的袋式除尘器基本几何参数包括灰斗、中箱体、上箱体具体尺寸。

3.根据权利要求2所述的袋式除尘器气流组织多参数优化方法，其特征在于：步骤2）中所述的响应曲面法以滤料厚度、安全高度、进口速度为响应因素，响应值为数值模拟所获得袋式除尘器气流组织分布的综合流量不均幅值。

4.根据权利要求3所述的袋式除尘器气流组织多参数优化方法，其特征在于：所述的步骤3）采用Solidworks建模单元建立步骤2）中各工况的三维几何模型；再打开所述的步骤4）中的网格划分单元，所述的网格划分单元采用ICEM-cfd拓扑网格划分，读取步骤3）中所述的三维几何模型，对读取的模型进行拓扑计算，依据实际几何模型边界对各部分进行part划分并命名，设定全局及各part网格划分尺寸，选择八叉树算法对计算模型进行拓扑网格划分，并对网格进行光滑处理，再对除尘器实际物理边界进行简化，分别设定袋式除尘器进口为速度进口边界条件、出口为压力出口边界条件、花板及壳体壁面为固壁边界条件、滤袋为跳跃多孔介质边界条件，导出网格文件再将网格文件导入步骤5）中的CFD求解单元Fluent中，采用标准k-ε双方程模型，对方程进行离散化处理，设定求解计算初始参数，选取数值计算迭代方法进行迭代求解计算，直至计算收敛，保存计算结果。