CN107644132B - 一种用于模拟带电除尘器除尘效率的计算方法 - Google Patents
一种用于模拟带电除尘器除尘效率的计算方法 Download PDFInfo
- Publication number
- CN107644132B CN107644132B CN201710839402.4A CN201710839402A CN107644132B CN 107644132 B CN107644132 B CN 107644132B CN 201710839402 A CN201710839402 A CN 201710839402A CN 107644132 B CN107644132 B CN 107644132B
- Authority
- CN
- China
- Prior art keywords
- dust
- electric field
- distribution
- removal efficiency
- simulating
- 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
Links
- 239000000428 dust Substances 0.000 title claims abstract description 160
- 238000004364 calculation method Methods 0.000 title claims abstract description 20
- 239000002245 particle Substances 0.000 claims abstract description 88
- 230000005684 electric field Effects 0.000 claims abstract description 51
- 230000005484 gravity Effects 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 8
- 230000009471 action Effects 0.000 claims abstract description 7
- 239000012530 fluid Substances 0.000 claims abstract description 6
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims description 11
- 239000003546 flue gas Substances 0.000 claims description 11
- 230000001133 acceleration Effects 0.000 claims description 4
- 239000000779 smoke Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 14
- 239000012071 phase Substances 0.000 description 5
- 239000012717 electrostatic precipitator Substances 0.000 description 4
- 230000005611 electricity Effects 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000007790 solid phase Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010410 dusting Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
Images
Landscapes
- Electrostatic Separation (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种用于模拟带电除尘器除尘效率的计算方法,包括以下步骤:建立静电除尘器的全尺寸几何模型,模拟静电力作用的电压分布;将步骤1中所述的电压分布的二维/三维的结果以文件的形式导出,且输出为数据流;将所述数据流导入到CFD仿真软件中,得到整个除尘通道的电压值分布;通过编写程序将步骤3得到的电压值分布转换为电场力分布;计算全尺寸几何模型中粉尘粒子的受力情况,包括重力,曳力和电场力;根据计算流体力学方法,根据粉尘粒子的数量和运动轨迹终结位置来计算带电除尘器的除尘效率。本发明将电场力分布的模拟引入到电场除尘效率模拟中,从而能够准确模拟带电除尘器除尘效率的计算方法。
Description
技术领域
本发明涉及一种用于模拟带电除尘器除尘效率的计算方法。
背景技术
目前对电除尘器的仿真中,一般都是根据计算流体力学原理,自行编程或者利用成熟的商业软件仅对流场进行模拟,通过对电除尘器内部空气动力场的模拟,得出气流分布方案。这种方法应用较多、较成熟。另外更深入仿真的还有利用多相流/离散项模型,在电除尘器内部空气动力场中,对气固两相进行仿真,同时考虑重力的影响,通过对不同粒径的固相粉尘颗粒在烟气中的运动轨迹的计算,得出粉尘颗粒的去除效率。而带电除尘器中除尘原理为:(1)施加高电压产生强场使气体电离,即产生电晕放电;(2)悬浮粒子荷电;(3)荷电尘粒在电场力作用下向电极运行;(4)荷电尘粒在电场中被捕集;(5)振打清灰。由此可见,粉尘被去除的关键是电场力的作用,因此,带电除尘器除尘效率的模拟中,如果不能对电场力进行模拟,这样的仿真结果是不准确的,不能为工程应用提供准确依据。
发明内容
为解决上述问题,本发明的目的在于提供一种用于模拟带电除尘器除尘效率的计算方法,为带电除尘器提供一种能够准确模拟粉尘粒子在电场力、重力、惯性力、曳力共同作用下的运动过程,得出除尘效率的仿真方法。
本发明的目的是通过以下技术方案实现的:
一种用于模拟带电除尘器除尘效率的计算方法,包括以下步骤:
步骤1,建立静电除尘器的全尺寸几何模型,模拟静电力作用的电压分布;
步骤2,将步骤1中所述的电压分布的二维/三维的结果以文件的形式导出,且输出为数据流;
步骤3,将所述数据流导入到CFD仿真软件中,得到整个除尘通道的电压值分布;
步骤4,通过编写程序将步骤3得到的电压值分布转换为电场力分布;
步骤5,计算全尺寸几何模型中粉尘粒子的受力情况,包括重力,曳力和电场力;
步骤6,根据计算流体力学方法,根据粉尘粒子的数量和运动轨迹终结位置来计算带电除尘器的除尘效率。
进一步的,步骤1中,所述全尺寸几何模型为二维/三维模型,且所述全尺寸几何模型包含阳极线和阴极;将所述全尺寸几何模型划分网格,将网格文件导入电场模拟程序中,根据电晕放电公式计算出二维/三维空间中的电压分布。
进一步的,所述电晕放电公式考虑阴极线半径、阴极线的粗糙系数、静电除尘器结构和空气相对密度。
进一步的,步骤2中,将静电除尘器单个除尘通道内电压分布以固定的格式输出为一个数据流,输出格式为:网格X坐标,网格Y坐标,网格Z坐标,网格中心点电压值;通过所述数据流表征整个除尘通道内的电压分布。
进一步的,步骤3中,在CFD仿真软件中,划分与所述全尺寸几何模型相同的网格,导入所述数据流,在CFD仿真软件中将每个网格中心点的电压值表征出来,得到整个除尘通道的电压值分布。
进一步的,步骤4包括:
步骤4.1,将步骤3中的电压值分布作为CFD仿真软件中的标量值,通过求得所述标量值的梯度得到静电除尘器内电场强度E的分布;
步骤4.2,由电场强度E的分布转换为电场力的分布,其中电场力F的计算公式为:
F=q·E
式中,E为电场强度,q为粉尘粒子的荷电量;
粉尘粒子的荷电量q按下式计算:
式中,ε0为真空介电常数,dp为粉尘粒子粒径,E为电场强度,εp为粉尘粒子的相对介电常数。
进一步的,步骤5中重力G、曳力FD和电场力F的计算方法分别如下:
式中,gx为重力加速度,ρp为粉尘粒子密度,ρ为烟气密度;
式中,μ为烟气粘度,CD为曳力系数,ρp为粉尘粒子密度,dp为粉尘粒子粒径,Re为相对雷诺数;
up为粉尘粒子速度,u为烟气速度,ρ为烟气密度,dp为粉尘粒子粒径;
F=q·E
式中,E为电场强度,q为粉尘粒子的荷电量。
进一步的,步骤6中:根据多相流/离散相模型,对粉尘粒子的边界条件进行设定,根据粉尘粒子的粒度分布对粉尘粒子直径进行输入;假设粒子运动到极板上就被捕集,运动轨迹结束;根据初始粉尘粒子总数和出口粉尘粒子总数,计算除尘效率,同时计算出不同粒径粉尘粒子的分级除尘效率。
本发明的有益效果为:
本发明将电场力分布的模拟引入到电场除尘效率模拟中,从而能够准确模拟带电除尘器除尘效率的计算方法。具体的,通过电场仿真软件对阴极线放电后的二维/三维空间中的电压分布进行模拟,将电压分布的二维/三维的结果以数据流的形式导出,然后将其导入到CFD仿真软件中,通过编写程序将电压分布转换为电场力分布,计算模型中粉尘粒子的受力情况,包括重力、曳力、电场力;根据计算流体力学方法,根据粉尘粒子的数量和运动轨迹终结位置来计算带电除尘器的除尘效率。
本发明不仅仅适合于静电除尘器的除尘效率模拟,也应用于电-袋复合式除尘器、湿式电除尘器除尘效率的模拟,以及所有利用静电去除粉尘的设备的除尘效率的模拟。
附图说明
图1为本发明所述计算方法的流程示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
一种用于模拟带电除尘器除尘效率的计算方法,包括以下步骤:
步骤1,建立静电除尘器的全尺寸几何模型,模拟静电力作用的电压分布;
因为每个除尘通道具有相似性,因此可以对静电除尘器一个通道内的电场进行模拟。建立二维/三维全尺寸几何模型,包含阳极线和阴极。将所述全尺寸几何模型划分网格,将网格文件导入已编号的电场模拟程序或者商业软件中,根据电晕放电公式计算出二维/三维空间中的电压分布。其中,所述电晕放电公式考虑阴极线半径、阴极线的粗糙系数、静电除尘器结构和空气相对密度。
步骤2,将步骤1中所述的电压分布的二维/三维的结果以文件的形式导出,且输出为数据流;
将静电除尘器单个除尘通道内电压分布以固定的格式输出为一个数据流,输出格式为:网格X坐标,网格Y坐标,网格Z坐标,网格中心点电压值;通过所述数据流表征整个除尘通道内的电压分布。
步骤3,将所述数据流导入到CFD仿真软件中,得到整个除尘通道的电压值分布;
在CFD仿真软件中,划分与所述全尺寸几何模型相同的网格,导入所述数据流,在CFD仿真软件中将每个网格中心点的电压值表征出来,得到整个除尘通道的电压值分布。
步骤4,通过编写程序将步骤3得到的电压值分布转换为电场力分布,包括:
步骤4.1,将步骤3中的电压值分布作为CFD仿真软件中的标量值,通过求得所述标量值的梯度得到静电除尘器内电场强度E的分布;
步骤4.2,由电场强度E的分布转换为电场力的分布,其中电场力F的计算公式为:
F=q·E
式中,E为电场强度,单位V/m;q为粉尘粒子的荷电量;
其中饱和状态下粉尘粒子的荷电量q按下式计算:
式中,ε0为真空介电常数,ε0=8.85×10-12C/N·m2;dp为粉尘粒子粒径,单位m;E为电场强度,单位V/m;εp为粉尘粒子的相对介电常数。
步骤5,计算全尺寸几何模型中粉尘粒子的受力情况,包括重力,曳力和电场力;
重力G:因为重力加速度而产生的力,受力方向为竖直向下。
式中,gx为重力加速度,ρp为粉尘粒子密度,ρ为烟气密度;
曳力FD:因为粒子与烟气的相对速度而产生的力,受力方向与相对速度方向相同。
式中,μ为烟气粘度,CD为曳力系数,ρp为粉尘粒子密度,dp为粉尘粒子粒径,Re为相对雷诺数;
up为粉尘粒子速度,u为烟气速度,ρ为烟气密度,dp为粉尘粒子粒径。
电场力F:因为电场强度而产生的力,受力方向与电场强度方向相同。
F=q·E
式中,E为电场强度,q为粉尘粒子的荷电量。
步骤6,根据计算流体力学方法,根据粉尘粒子的数量和运动轨迹终结位置来计算带电除尘器的除尘效率。
根据多相流/离散相模型,对粉尘粒子的边界条件进行设定,根据粉尘粒子的粒度分布对粉尘粒子直径进行输入,粉尘粒子粒径很重要,影响粉尘粒子所受的电场力,因此粒径是计算除尘效率的关键参数;假设粒子运动到极板上就被捕集,运动轨迹结束;根据初始粉尘粒子总数和出口粉尘粒子总数,计算除尘效率,同时计算出不同粒径粉尘粒子的分级除尘效率。
本发明不仅仅适合于静电除尘器的除尘效率模拟,也应用于电-袋复合式除尘器、湿式电除尘器除尘效率的模拟,以及所有利用静电去除粉尘的设备的除尘效率的模拟。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种用于模拟带电除尘器除尘效率的计算方法,其特征在于,包括以下步骤:
步骤1,建立带电除尘器的全尺寸几何模型,模拟带电力作用的电压分布;
步骤2,将步骤1中所述的电压分布的二维/三维的结果以文件的形式导出,且输出为数据流;
步骤3,将所述数据流导入到CFD仿真软件中,得到整个除尘通道的电压值分布;
步骤4,通过编写程序将步骤3得到的电压值分布转换为电场力分布;
步骤5,计算全尺寸几何模型中粉尘粒子的受力情况,包括重力,曳力和电场力;重力G、曳力FD和电场力F的计算方法分别如下:
式中,gx为重力加速度,ρp为粉尘粒子密度,ρ为烟气密度;
式中,μ为烟气粘度,CD为曳力系数,ρp为粉尘粒子密度,dp为粉尘粒子粒径,Re为相对雷诺数;
up为粉尘粒子速度,u为烟气速度,ρ为烟气密度,dp为粉尘粒子粒径;
F=q·E
式中,E为电场强度,q为粉尘粒子的荷电量;
步骤6,根据计算流体力学方法,根据粉尘粒子的数量和运动轨迹终结位置来计算带电除尘器的除尘效率。
2.根据权利要求1所述的用于模拟带电除尘器除尘效率的计算方法,其特征在于,步骤1中,所述全尺寸几何模型为二维/三维模型,且所述全尺寸几何模型包含阳极线和阴极;将所述全尺寸几何模型划分网格,将网格文件导入电场模拟程序中,根据电晕放电公式计算出二维/三维空间中的电压分布。
3.根据权利要求2所述的用于模拟带电除尘器除尘效率的计算方法,其特征在于,所述电晕放电公式考虑阴极线半径、阴极线的粗糙系数、带电除尘器结构和空气相对密度。
4.根据权利要求2所述的用于模拟带电除尘器除尘效率的计算方法,其特征在于,步骤2中,将带电除尘器单个除尘通道内电压分布以固定的格式输出为一个数据流,输出格式为:网格X坐标,网格Y坐标,网格Z坐标,网格中心点电压值;通过所述数据流表征整个除尘通道内的电压分布。
5.根据权利要求4所述的用于模拟带电除尘器除尘效率的计算方法,其特征在于,步骤3中,在CFD仿真软件中,划分与所述全尺寸几何模型相同的网格,导入所述数据流,在CFD仿真软件中将每个网格中心点的电压值表征出来,得到整个除尘通道的电压值分布。
7.根据权利要求1所述的用于模拟带电除尘器除尘效率的计算方法,其特征在于,步骤6中:根据多相流/离散相模型,对粉尘粒子的边界条件进行设定,根据粉尘粒子的粒度分布对粉尘粒子直径进行输入;假设粒子运动到极板上就被捕集,运动轨迹结束;根据初始粉尘粒子总数和出口粉尘粒子总数,计算除尘效率,同时计算出不同粒径粉尘粒子的分级除尘效率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710839402.4A CN107644132B (zh) | 2017-09-18 | 2017-09-18 | 一种用于模拟带电除尘器除尘效率的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710839402.4A CN107644132B (zh) | 2017-09-18 | 2017-09-18 | 一种用于模拟带电除尘器除尘效率的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107644132A CN107644132A (zh) | 2018-01-30 |
CN107644132B true CN107644132B (zh) | 2021-06-11 |
Family
ID=61111904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710839402.4A Active CN107644132B (zh) | 2017-09-18 | 2017-09-18 | 一种用于模拟带电除尘器除尘效率的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107644132B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109759235B (zh) * | 2019-01-25 | 2020-05-12 | 国网江苏省电力有限公司检修分公司 | 一种变压器防尘装置及其集尘极电压设计方法 |
CN109883772B (zh) * | 2019-01-31 | 2021-02-26 | 中煤科工集团重庆研究院有限公司 | 一种呼吸性粉尘连续分离装置分级效率优化方法 |
CN112760125B (zh) * | 2019-11-06 | 2022-07-26 | 中国石油化工股份有限公司 | 存储器、电脱盐脱水效率的计算方法、系统和装置 |
CN112016808B (zh) * | 2020-07-31 | 2022-05-03 | 山西晋浙环保科技有限公司 | 高效静电除尘器 |
CN113033117B (zh) * | 2021-03-09 | 2024-03-19 | 江苏大学 | 一种运动荷电液滴诱导电场强度及电场力计算方法和系统 |
CN115294839B (zh) * | 2021-12-22 | 2024-06-11 | 浙江安防职业技术学院 | 模拟异极性荷电粒子在磁电凝并器中的运动情况的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102284360A (zh) * | 2011-08-18 | 2011-12-21 | 上海电力学院 | 静电除尘器可吸入颗粒物pm10捕集效率的确定方法 |
CN104785374A (zh) * | 2015-03-24 | 2015-07-22 | 福建龙净环保股份有限公司 | 一种电除尘器性能检测方法及装置 |
CN106951603A (zh) * | 2017-02-28 | 2017-07-14 | 上海电力学院 | 基于power函数的pm2.5捕集效率贡献率曲线拟合方法 |
CN106980710A (zh) * | 2017-02-28 | 2017-07-25 | 上海电力学院 | Esp工作电压贡献率与磁感应强度曲线拟合方法 |
CN106980712A (zh) * | 2017-03-06 | 2017-07-25 | 三峡大学 | 一种山火条件下输电线路间隙空间合成电场计算方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102631990B (zh) * | 2012-04-17 | 2015-09-23 | 江苏中科睿赛污染控制工程有限公司 | 一种纳微颗粒高效捕集方法及装置 |
KR101373720B1 (ko) * | 2013-05-13 | 2014-03-13 | 한국기계연구원 | 활성탄 섬유필터를 이용하는 전기집진기 |
CN107138016A (zh) * | 2017-05-08 | 2017-09-08 | 浙江大学 | 基于多场促进的颗粒物凝并团聚强化装置及方法 |
-
2017
- 2017-09-18 CN CN201710839402.4A patent/CN107644132B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102284360A (zh) * | 2011-08-18 | 2011-12-21 | 上海电力学院 | 静电除尘器可吸入颗粒物pm10捕集效率的确定方法 |
CN104785374A (zh) * | 2015-03-24 | 2015-07-22 | 福建龙净环保股份有限公司 | 一种电除尘器性能检测方法及装置 |
CN106951603A (zh) * | 2017-02-28 | 2017-07-14 | 上海电力学院 | 基于power函数的pm2.5捕集效率贡献率曲线拟合方法 |
CN106980710A (zh) * | 2017-02-28 | 2017-07-25 | 上海电力学院 | Esp工作电压贡献率与磁感应强度曲线拟合方法 |
CN106980712A (zh) * | 2017-03-06 | 2017-07-25 | 三峡大学 | 一种山火条件下输电线路间隙空间合成电场计算方法 |
Non-Patent Citations (3)
Title |
---|
太阳能静电除尘器设计与性能研究;张华;《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》;20170215(第2期);B027-1846 * |
脉冲放电电凝并结合碱液吸收烟气多种污染物协同脱除研究;徐飞;《中国博士学位论文全文数据库 工程科技Ⅰ辑》;20120115(第1期);B027-102 * |
范德华力和静电力下的细颗粒离散动力学研究;柳冠青;《中国博士学位论文全文数据库 基础科学辑》;20121115(第11期);A004-3 * |
Also Published As
Publication number | Publication date |
---|---|
CN107644132A (zh) | 2018-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107644132B (zh) | 一种用于模拟带电除尘器除尘效率的计算方法 | |
Farnoosh et al. | 3-D numerical analysis of EHD turbulent flow and mono-disperse charged particle transport and collection in a wire-plate ESP | |
Long et al. | Evaluation of various particle charging models for simulating particle dynamics in electrostatic precipitators | |
Nikas et al. | Numerical simulation of the flow and the collection mechanisms inside a laboratory scale electrostatic precipitator | |
Long et al. | Numerical simulation of the flow and the collection mechanism inside a scale hybrid particulate collector | |
Varonos et al. | Prediction of the cleaning efficiency of an electrostatic precipitator | |
CN106980710A (zh) | Esp工作电压贡献率与磁感应强度曲线拟合方法 | |
CN104785374B (zh) | 一种电除尘器性能检测方法及装置 | |
CN106951603A (zh) | 基于power函数的pm2.5捕集效率贡献率曲线拟合方法 | |
Liu et al. | Numerical analysis of charged particle collection in wire-plate ESP | |
JPWO2016021063A1 (ja) | 粒子荷電装置 | |
Zhang et al. | Theoretical and experimental study of trapping PM2. 5 particles via magnetic confinement effect in a multi-electric field ESP | |
Ning et al. | Electrode configurations inside an electrostatic precipitator and their impact on collection efficiency and flow pattern | |
Zhang et al. | A numerical simulation of diffusion charging effect on collection efficiency in wire-plate electrostatic precipitators | |
Zhang et al. | Modeling and simulation of PM 2.5 collection efficiency in a wire‐plate ESP subjected to magnetic field and diffusion charging | |
Bacchiega et al. | Back-corona model for prediction of ESP efficiency and voltage-current characteristics | |
Feng et al. | A critical review of models used in numerical simulation of electrostatic precipitators | |
Haque et al. | A numerical model of an electrostatic precipitator | |
Kantouna et al. | Analysis of a cylinder-wire-cylinder electrode configuration during corona discharge | |
González et al. | Study and Numerical Simulation of Particle Collection and Trajectory of Particle Motion in a Duct-type Electrostatic Precipitator with Seven Circular Corona Wires | |
Shen et al. | Electrical Characteristics and Electrohydrodynamic Flows in Electrostatic Precipitator of Six Shaped Discharge Electrodes | |
Sayed et al. | A new approach for accurate prediction of optimal ion mobility in dust loaded electrostatic precipitator using FDM-FMG | |
Khare et al. | Computer-aided simulation of efficiency of an electrostatic precipitator | |
Vereshchagin et al. | A new approach to calculation of flue gas scrubbing efficiency in electrostatic precipitators | |
CN110909469B (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 |