CN101822957B - 基于统计学及遍历性理论的判断多相混合均匀性的方法 - Google Patents
基于统计学及遍历性理论的判断多相混合均匀性的方法 Download PDFInfo
- Publication number
- CN101822957B CN101822957B CN201010140233.3A CN201010140233A CN101822957B CN 101822957 B CN101822957 B CN 101822957B CN 201010140233 A CN201010140233 A CN 201010140233A CN 101822957 B CN101822957 B CN 101822957B
- Authority
- CN
- China
- Prior art keywords
- time
- tomography
- mixing
- mixing uniformity
- tracer grain
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000003325 tomography Methods 0.000 claims abstract description 26
- 239000000700 radioactive tracer Substances 0.000 claims abstract description 25
- 239000000463 material Substances 0.000 claims description 17
- 238000003756 stirring Methods 0.000 claims description 10
- 238000002591 computed tomography Methods 0.000 claims description 9
- 238000003384 imaging method Methods 0.000 claims description 9
- 238000011160 research Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000010348 incorporation Methods 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000002603 single-photon emission computed tomography Methods 0.000 claims 3
- 238000013461 design Methods 0.000 abstract description 6
- 239000002245 particle Substances 0.000 abstract description 4
- 239000000126 substance Substances 0.000 abstract description 4
- 238000013019 agitation Methods 0.000 abstract 4
- 239000012530 fluid Substances 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000002146 bilateral effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 1
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 1
- 239000005864 Sulphur Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000008367 deionised water Substances 0.000 description 1
- 229910021641 deionized water Inorganic materials 0.000 description 1
- 238000006477 desulfuration reaction Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000006194 liquid suspension Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000010310 metallurgical process Methods 0.000 description 1
- 238000005272 metallurgy Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Chemical compound O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种基于统计学及遍历性理论的判断多相混合均匀性的方法。具体是：(1)在混合过程中放入搅拌槽内一个示踪颗粒作为研究对象；(2)利用电子断层等成像技术来获得任一个断层的连续的搅拌实时图样；(3)根据断层面的尺寸大小，将获得的多相混合过程中该断层的连续的搅拌实时图样进行网格划分；(4)运用统计学的方法计算示踪颗粒在一定时间内经过1～n个区域中每一个区域的次数，然后根据穿过1～n个区域的总次数计算该示踪颗粒穿过每一个区域中的概率P；(5)借助混合遍历性理论，得出多相混合在t时刻之后达到混合均匀的状态。本发明方法简单方便，对化工、冶金搅拌设备中判断混合均匀性及理论上指导搅拌反应器的设计，提供了一种可靠实用的方法。
Description
技术领域
本发明涉及属于化学工程技术领域，具体地说是一种适应于冶金、化工领域所有多相混合过程混合匀性的判断方法。
背景技术
搅拌是化工、冶金过程中最常见的操作之一。大约在60到80年代期间流体混合技术得到了飞速的发展，其研究的重点主要是针对常规搅拌桨在低粘和高粘非牛顿均相体系、固液悬浮和气液分散等非均相体系中的搅拌功耗、混合时间等宏观量进行实验研究。虽然有大量的设计经验和关联式的可用于分析和预测混合体系，但将搅拌反应器从实验室规模直接放大到大规模的工业生产中，仍是没有把握的，至今仍需要通过逐级放大来达到搅拌设备所要求的传质、传热和混合。这种方法不但耗费财力和大量的人力和物力，而且设计周期很长，据相关统计显示美国的化学工业由于搅拌反应器设计不合理所造成的损失每年约为10-100亿美元。因此采用先进的测试手段和建立合理的数学模型获取搅拌槽中的速度场、温度场和浓度场，以及采取有效的衡量流体混合效果的方法，不仅对混合设备的优化设计具有十分重要的经济意义，而且对放大和混合的基础研究具有现实的理论意义。近年来，随着科技的发展，出现了激光多普勒测速仪LDV和计算流体力学CFD模拟技术以及电子断层成像技术等的广泛应用，促进了流体混合技术飞快的发展，就目前看，衡量流体混合效果的方法众多，主要有电导率法、热电偶法、光学法、脱色法等，其中电导率法在低粘性流体的搅拌混合中应用很广，但对搅拌介质的要求高，如要求搅拌介质为去离子水；热电偶法则对流场会产生破坏；一般的光学法由于装置复杂，应用也较少；脱色法用来测高粘流体和粘弹性流体非常有效，但由于采用肉眼判断，因而带有较强的主观性。
发明内容
本发明克服了上述方法的不足，提供一种具有较高的应用价值的、简便可行的用于判断多相混合过程混合均匀性的新方法。
本发明判断多相混合均匀性的方法的技术方案是：采用统计学及遍历性理论来判断多相混合均匀性，具体步骤为：
(1)在混合过程中放入搅拌槽内一个区别于周围其它物质且不跟任何物质发生化学反应的示踪颗粒作为研究对象；
(2)利用电子断层成像技术EPT，来获得多相混合过程中任一个断层的连续的搅拌实时图样；
(3)根据断层面的尺寸大小，将获得的多相混合过程中该断层的连续的搅拌实时图样进行网格划分，一般划分为n2个区域，n为自然数，并对其编号1～n；
(4)运用统计学的方法计算示踪颗粒在一定时间内经过1～n个区域中每一个区域的次数，然后根据穿过1～n个区域的总次数计算该示踪颗粒穿过每一个区域中的概率P，记作Pi，i＝1，2，3……n；
(5)借助混合遍历性理论，若经历充分的混合时间t后，Pi的数值大小相等或者Pi的标准差充分的小即趋近于0，若画出标准差的变化曲线图，则在t时刻后，标准差就几乎为零且往后的时间几乎保持不变，则表示多相混合在t时刻之后达到混合均匀的状态。
示踪颗粒的选择必须满足异于周围的任何物体，其中密度必须不同于周围物体，形状最好是球形便于捕捉。
所述电子断层成像技术EPT，英文为(Electrical Process Tomography，包括：电阻断层成像系统ERT即Electrical Resistance Tomography)、电容断层成像系统ECT即Electrical Capacitance Tomography)和电磁断层成像系统EMT，即Electro-magnetic Tomography)三类。
混合过程中获得断连续的搅拌实时图样进行网格划分时，区域的划分原则是根据图样和示踪颗粒大小、混合均匀性判断精度要求来确定n的取值大小。
在系统科学里，直观地讲，遍历变换意味着从给定点出发，经过多次迭代之后可接近任意点，即各态历经混合变换意味着从任何集合A出发的点，经过多次迭代之后所得之点属于集合B的可能性为两个集合测度之积，且与A，B的位置无关，即混合均匀，正合变换则意味着从任何非零测度集合出发，经过多次迭代之后所得集合的点将充满整个空间，即充分扩散，正合蕴含遍历性。简而言之，这种遍历性使得空间中的定点可接近任意点，即到达任意点任意区域的概率是几乎一样的，我们称之为混合均匀。
具体方法是：(1)在混合过程中放入搅拌槽内一个区别于周围其它物质且不跟任何物质发生化学反应的示踪颗粒作为研究对象；(2)利用电子断层成像技术EPT(Electrical Process Tomography)，其中包括电阻断层成像系统ERT(Electrical Resistance Tomography)、电容断层成像系统ECT(ElectricalCapacitance Tomography)和电磁断层成像系统EMT(Electro-magneticTomography)三类，来获得多相混合过程中任一个断层的连续的搅拌实时图样；(3)根据断层面的尺寸大小，将获得的多相混合过程中该断层的连续的搅拌实时图样进行网格划分，一般划分为n2个区域(n为自然数)，并对其编号1～n；(4)运用统计学的方法计算示踪颗粒在一定时间内经过1～n个区域中每一个区域的次数，然后根据穿过1～n个区域的总次数计算该示踪颗粒穿过每一个区域中的概率P，记作Pi，i＝1，2，3……n；(5)借助混合遍历性理论，若经历充分的混合时间t后，Pi的数值大小相等或者Pi的标准差充分的小(趋近于0)，若画出标准差的变化曲线图，则在t时刻后，标准差就几乎为零且往后的时间几乎保持不变，则表示多相混合在t时刻之后达到混合均匀的状态。
本发明的有益效果是：
1、该方法简单可行，有系统科学之遍历性理论做支撑；
2、解决了电导率法、热电偶法、光学法、脱色法等的不足之处；
3、可以减少搅拌反应器因不合理的设计造成的经济损失。
附图说明
图1为本发明的电子断层成像系统搅拌设备混合图样采集装置图。
图2为本发明的断层图样网格划分图。
图3为本发明的示踪颗粒穿过区域12时刻的断层图。
图4为本发明的记录示踪颗粒穿过每一个区域的概率Pi的标准差随时间变化图。
图5为本发明的示踪颗粒穿过区域29时刻的断层图。
图6为本发明的记录示踪颗粒穿过每一个区域的概率Pi的标准差随时间变化图。
在图1中搅拌设备两侧的黑色点代表了八只传感器，用于信号的接受和发送；其中，1为浆，2为搅拌设备，3为断层面，4为数据获取设备，5为图像重构设备。
具体实施方式
实施例1：某食品加工厂，运用机械搅拌的方式进行混合物料的操作，运用如图1所示装置，在搅拌设备两侧对称安装有8只传感器，用于搅拌设备内部竖断面连续图像的采集，通过判断出不同质量物料在不同的搅拌速度下的混合均匀时间，从而针对不同质量的物料可以选出最佳的搅拌速度，即在最短的时间内让搅拌槽内达到均匀状态，有利于提高食品加工效率。
实施步骤如下：1.在搅拌设备中添加一个异于物料密度的球状材料作为示踪颗粒。2.运用图1中所示电子断层成像系统采集竖直断面的连续图样，记录下每一时刻的断层图样，根据搅拌设备竖直断面实际尺寸，我们划分为42个区域，n＝4。3.计算并统计示踪颗粒穿过1～16个网格区域的次数，在t＝5s的时候，P1＝9/71，P2＝3/71，P3＝0/71，P4＝4/71，P5＝14/71，P6＝4/71，P7＝3/71，P8＝7/71，P9＝0/71，P10＝6/71，P11＝3/71，P12＝2/71，P13＝12/71，P14＝0/71，P15＝0/71，P16＝4/71，计算标准差得0.0597，继续记录穿过每个区域的次数，在t＝25s后，P1＝90/1455，P2＝93/1455，P3＝89/1455，P4＝89/1455，P5＝94/1455，P6＝94/1455，P7＝93/1455，P8＝91/1455，P9＝88/1455，P10＝87/1455，P11＝93/1455，P12＝92/1455，P13＝90/1455，P14＝86/1455，P15＝90/1455，P16＝96/1455，计算标准差得0.0019，一直记录下去，往后的数据显示标准差已经几乎不发生变化，所以，我们分析当t＝25s以后，混合进入均匀状态。如图4所示。
实施例2：某钢铁厂，运用机械搅拌的方式进行除硫的操作，运用如图1所示装置，在搅拌设备两侧对称安装有8只传感器，用于搅拌设备内部竖断面连续图像的采集，通过判断出不同质量物料在不同的搅拌速度下的混合均匀时间，从而针对不同质量的物料可以选出最佳的搅拌速度，即在最短的时间内让搅拌槽内达到均匀状态，有利于提高脱硫效率。具体实施步骤如下：1.在搅拌设备中添加一个异于物料密度的球状材料作为示踪颗粒，该颗粒不熔于物料内部。2.运用图1中所示电子断层成像系统采集竖直断面的连续图样，记录下每一时刻的断层图样，根据搅拌设备竖直断面实际尺寸，我们划分为82个区域，n＝8如图5所示，t时刻示踪颗粒位于区域29处。3.计算并统计示踪颗粒穿过1～64个网格区域的次数，在t＝1s的时候分别计算P1～P64数值大小，计算标准差得0.089，继续记录穿过每个区域的次数，在t＝250s后，再次计算P1～P64的数值大小，计算标准差得0.019，在t＝300s后，再次计算P1～P64的数值大小，计算标准差得0.009一直记录下去，往后的数据显示标准差已经几乎不发生变化，所以，我们分析当t＝300s以后，混合进入均匀状态。如图6所示标准差的演化曲线图。
Claims (3)
1.一种判断多相混合均匀性的方法，其特征在于：该方法采用统计学及遍历性理论来判断多相混合均匀性，具体步骤为：
(1)在混合过程中放入搅拌槽内一个区别于周围其它物质且不跟任何物质发生化学反应的示踪颗粒作为研究对象；
(2)利用电子断层成像技术EPT，来获得多相混合过程中任一个断层的连续的搅拌实时图样；
(3)根据断层面的尺寸大小，将获得的多相混合过程中该断层的连续的搅拌实时图样进行网格划分，一般划分为n2个区域，n为自然数，并对其编号1-n2；
(4)运用统计学的方法计算示踪颗粒在一定时间内经过1～n个区域中每一个区域的次数，然后根据穿过1-n2个区域的总次数计算该示踪颗粒穿过每一个区域中的概率P，记作Pi，i＝1，2，3……n；
(5)借助混合遍历性理论，若经历充分的混合时间t后，Pi的数值大小相等或者Pi的标准差充分的小即趋近于0，若画出标准差的变化曲线图，则在t时刻后，标准差就几乎为零且往后的时间几乎保持不变，则表示多相混合在t时刻之后达到混合均匀的状态；
示踪颗粒的选择必须满足异于周围的任何物体，其中密度必须不同于周围物体，形状是球形便于捕捉。
2.根据权利要求1所述的判断多相混合均匀性的方法，其特征在于：所述电子断层成像技术EPT，包括：电阻断层成像系统ERT、电容断层成像系统ECT和电磁断层成像系统EMT三类。
3.根据权利要求1所述的判断多相混合均匀性的方法，其特征在于：混合过程中获得断连续的搅拌实时图样进行网格划分时，区域的划分原则是根据图样和示踪颗粒大小、混合均匀性判断精度要求来确定n的取值大小。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201010140233.3A CN101822957B (zh) | 2010-04-07 | 2010-04-07 | 基于统计学及遍历性理论的判断多相混合均匀性的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201010140233.3A CN101822957B (zh) | 2010-04-07 | 2010-04-07 | 基于统计学及遍历性理论的判断多相混合均匀性的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN101822957A CN101822957A (zh) | 2010-09-08 |
| CN101822957B true CN101822957B (zh) | 2012-12-19 |
Family
ID=42687247
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201010140233.3A Expired - Fee Related CN101822957B (zh) | 2010-04-07 | 2010-04-07 | 基于统计学及遍历性理论的判断多相混合均匀性的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN101822957B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI472371B (zh) * | 2012-07-09 | 2015-02-11 | China Steel Corp | 物料混合方法 |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102087220B (zh) * | 2010-12-24 | 2012-08-15 | 江苏大学 | 液体混合浓度场检测方法 |
| CZ306017B6 (cs) * | 2014-12-15 | 2016-06-22 | Vysoká Škola Báňská-Technická Univerzita Ostrava | Validační skladovací zařízení pro měření tokových procesů sypkých hmot metodou elektrické kapacitní tomografie |
| CN105067621A (zh) * | 2015-08-17 | 2015-11-18 | 云南财经大学 | 一种判断多相混合均匀性的方法 |
| CN105910996A (zh) * | 2016-04-07 | 2016-08-31 | 昆明理工大学 | 一种利用rgb颜色模型确定多相流混合时间的方法 |
| CN112669268B (zh) * | 2020-12-21 | 2023-06-06 | 昆明理工大学 | 一种多相混合浓度场分布均匀性评价方法、系统和终端 |
| CN112871015B (zh) * | 2021-03-26 | 2022-03-08 | 福州大学 | 立式行星搅拌机搅拌臂长度与安装角的优化方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1416781A (zh) * | 2002-12-30 | 2003-05-14 | 北京质子科技开发有限公司 | 正电子断层扫描中热源高分辨快速图像迭代重建方法 |
| CN1640361A (zh) * | 2005-01-06 | 2005-07-20 | 东南大学 | 多相水平集的正电子断层扫描重建方法 |
| CN1963403A (zh) * | 2006-11-30 | 2007-05-16 | 天津大学 | 基于截面测量的气液两相流测量方法及装置 |
| CN101650292A (zh) * | 2009-09-02 | 2010-02-17 | 昆明理工大学 | 一种衡量流体混合效果的方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH07270227A (ja) * | 1994-03-31 | 1995-10-20 | Hitachi Ltd | 流体・構造連成振動解析方法 |
| JPH0933420A (ja) * | 1995-07-19 | 1997-02-07 | Honda Motor Co Ltd | 無機粉末と有機バインダーの混練状態または顆粒状態の評価方法、および混練方法 |
-
2010
- 2010-04-07 CN CN201010140233.3A patent/CN101822957B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1416781A (zh) * | 2002-12-30 | 2003-05-14 | 北京质子科技开发有限公司 | 正电子断层扫描中热源高分辨快速图像迭代重建方法 |
| CN1640361A (zh) * | 2005-01-06 | 2005-07-20 | 东南大学 | 多相水平集的正电子断层扫描重建方法 |
| CN1963403A (zh) * | 2006-11-30 | 2007-05-16 | 天津大学 | 基于截面测量的气液两相流测量方法及装置 |
| CN101650292A (zh) * | 2009-09-02 | 2010-02-17 | 昆明理工大学 | 一种衡量流体混合效果的方法 |
Non-Patent Citations (1)
| Title |
|---|
| 冯连芳，等.《流体混合技术新进展》.《化学工程》.2002,第30卷(第2期),70-74. * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI472371B (zh) * | 2012-07-09 | 2015-02-11 | China Steel Corp | 物料混合方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN101822957A (zh) | 2010-09-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101822957B (zh) | 基于统计学及遍历性理论的判断多相混合均匀性的方法 | |
| Laurent et al. | A protocol for the use of computational fluid dynamics as a supportive tool for wastewater treatment plant modelling | |
| CN101813641B (zh) | 用于验证多相搅拌混合均匀状态及程度的方法 | |
| CN105067621A (zh) | 一种判断多相混合均匀性的方法 | |
| CN101650292B (zh) | 一种衡量流体混合效果的方法 | |
| Parvareh et al. | CFD and ERT investigations on two-phase flow regimes in vertical and horizontal tubes | |
| CN104050330B (zh) | 升流式厌氧发酵生物制氢反应器的优化设计方法及其应用 | |
| Pakzad et al. | Evaluation of the mixing of non-Newtonian biopolymer solutions in the reactors equipped with the coaxial mixers through tomography and CFD | |
| Szalai et al. | Mixing analysis in a tank stirred with Ekato Intermig® impellers | |
| Bitog et al. | Optimised hydrodynamic parameters for the design of photobioreactors using computational fluid dynamics and experimental validation | |
| CN108504548A (zh) | 一种全混式厌氧发酵反应器的设计优化方法 | |
| Yang et al. | Decrease in light/dark cycle of microalgal cells with computational fluid dynamics simulation to improve microalgal growth in a raceway pond | |
| CN108760972A (zh) | Y型交汇河道污染物输移掺混模拟实验装置及试验方法 | |
| Li et al. | A critical review of experimental and CFD techniques to characterize the mixing performance of anaerobic digesters for biogas production | |
| Sardeshpande et al. | Electrical resistance tomography for gas holdup in a gas-liquid stirred tank reactor | |
| Mohammadrezaei et al. | Optimum mixing rate in biogas reactors: Energy balance calculations and computational fluid dynamics simulation | |
| CN110728423A (zh) | 一种长江流域水系统综合模拟方法及装置 | |
| Liu et al. | A novel 2-D ring-tubes model and numerical investigation of heat and moisture transfer around helix ground heat exchanger | |
| Hao et al. | An analytical model for micro-droplet steady movement on the hydrophobic wall of a micro-channel | |
| Fan et al. | CFD simulation of flow pattern in a bubble column reactor for forming aerobic granules and its development | |
| CN103279433A (zh) | 一种表征颗粒团聚体的方法及装置 | |
| Yang et al. | Topological approach for the measurement of mixing state quality in a vertical rectangular channel | |
| Xia et al. | Computational investigation of fluid dynamics in a recently developed centrifugal impeller bioreactor | |
| Zhang et al. | Bubble size modeling approach for the simulation of bubble columns | |
| Liu et al. | Lagrangian particle tracking in mechanically agitated polydisperse suspensions: Multi-component hydrodynamics and spatial distribution |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20121219 Termination date: 20150407 |
|
| EXPY | Termination of patent right or utility model |