具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述。
图1为逆流式自然通风湿式冷却塔的结构示意图。如图1所示,自然塔淋水面积是指冷却塔内“填料区顶部”的断面面积,按淋水面积A冷却塔可初步划分为:
小型塔A<4000m2
中型塔4000m2≤A<7000m2
大型塔7000m2≤A<12000m2
超大型塔A≥12000m2(塔底部直径D>110m、进风口高度h>11m)
自然通风逆流式冷却塔的散热包含了淋水区(图中配水区与填料之间的区域)、填料区和雨区(图中填料区以下、塔底水池水面以上的部分称为雨区)三个部分。
目前,冷却塔热力设计计算方法有一维仿真方法、二维仿真方法、三维仿真方法。一维仿真方法计算时将淋水区、填料区、雨区等三个区作为一个区处理,即将室内淋水填料实测的热力特性作为冷却塔的热力特性进行冷却塔的设计;二维仿真方法可以计算喷淋及填料区、雨区的每一点的热交换,可以对塔内的非均匀配水配风进行优化计算,但是二维仿真方法比较复杂,一般适用于冷却塔的局部优化;三维仿真方法尚处在初级阶段,就是对冷却塔进行三维数值建模,生成冷却塔整体计算域的三维计算网格,一般适用于冷却塔特定工况的影响分析。
随着冷却塔的淋水面积、塔高、塔底部直径增大,塔的进风口高度也响应增高,雨区的高度也随之增高,特别是超大型冷却塔底部直径一般大于140m、进风口高度大于12.5m时,若仍沿用现有的修正一维热力计算(拟二维计算)公式进行冷却塔的设计计算,将会使冷却塔的设计工况与实际运行工况偏离增大,严重的会造成设计不合理甚至。
图2为本发明超大型湿式冷却塔雨区热力特性的仿真计算方法的雨区模拟模型示意图。由于超大型冷却塔的体积庞大,内部结构复杂,本发明的超大型湿式冷却塔雨区热力特性的仿真计算方法,建立超大型冷却塔内部雨区的物理模型进行仿真研究。
图3为本发明超大型湿式冷却塔雨区热力特性的仿真计算方法的流程图。下面结合图2、图3对本发明模拟仿真方法做进一步说明。
S101:建立冷却塔的雨区模拟模型,其中,设置该模型的底部直径大于110米、进风口高度大于11米;
S102:在预设不同工况的条件下,计算所述雨区模拟模型的雨滴当量直径;
S103:根据所述雨滴当量直径进行运动仿真,获取所述雨区模拟模型的热力特性数学模型。
对于“底部直径一般大于110m、进风口高度大于11m”的超大型逆流式自然通风湿式冷却塔,要获得冷却塔的雨区热力特性,通过原型观测的方法显然是不现实,须通过一套模拟仿真的方法来取得冷却塔的雨区热力特性。
雨区热力特性模拟模型如图2所示。该模型主要由水循环和空气循环两个模拟系统组成,各模拟系统分别介绍如下:
水循环模拟系统:水从集水池经水泵提升后进入电加热器进行加热,然后送至配水池内,配水池底部设有孔,孔下面有一层淋水填料将水散开形成淋水区,填料内不能通风没有热交换,淋水区的水滴回落在装置的下方集水池中循环利用。
空气循环模拟系统:由装置尾部的抽风机从进风口吸入外界空气,进入塔的测试段与水进行气、热交换后,再从风机出口排出。进塔风量通过变频器改变风机的电机的转速来控制。
在装置上布置进塔空气的干、湿球温度测点、配水池上布置进塔水温测点、集水池中布置出塔水温的测点、在出口处布置出塔气温测点。温度采用多点温度测量系统进行测量,温度测量精度为0.2℃;循环水流量在进塔竖管上布置电磁流量计,进塔风量在进风口采用风速表进行测量。
建立冷却塔的雨区模拟模型的步骤,还包括:
设置所述雨区模拟模型的平面尺寸为1m×1m、高度15m、淋水高度为5m~15m、淋水密度在7-13m3/m2·h、风速为0.6m/s~2.8m/s。
测试时不安装填料只淋水,在模拟装置上测出:循环水量、进水温度、出水温度、进气量、进气干球温度、湿球温度、大气压力,根据这些测试数据,可以算出冷却数及容积散质系数。
(1)冷却数N
式中:i″—与水温相应的饱和空气比焓(kJ/kg);
i—空气比焓(kJ/kg);
Q—循环水量,Kg/s;
Cw—水的比热(kJ/(kg·℃));
t1、t2—进、出塔水温(℃);
h—淋水高度(即模拟进风口高度)(m);
A—淋水面积(m2)。
(2)容积散质系数Ka(kg/(m3·h))
式中:q—淋水密度(kg/(m2·h));
h-淋水高度(m)。
计算所述雨区模拟模型的雨滴当量直径的步骤之前,包括:
预设数个淋水高度,在各个淋水高度下进行测试;
采集循环水量、进水温度、出水温度、进气量、进气干球温度、湿球温度、大气压力进行雨滴当量直径的计算;
将计算的结果存储为雨滴当量直径表,该表包含了淋水高度、冷却数、气水比、风速、淋水密度以及对应的雨滴当量直径。
优选地,分别对淋水高度5.0m、10.0m和14.5m进行测试,按式(2-1)和式(2-2)可整理出不同淋水高度的雨区的热力特性关系为:
淋水高度为5.0m时:
N=0.246λ0.227 (2-3)
淋水高度为10.0m时:
N=0.422λ0.214 (2-4)
淋水高度为14.5m时:
N=0.560λ0.311 (2-5)
假定淋水的雨滴为不变型的刚性球,刚性球在空气中运动的表面面积蒸发系数为:
水滴容积散质系数为:
式中:v—气流流过水滴的速度,m/s;
d—水滴直径,m;
ρw—水的密度,kg/m3;
vf-水滴下降速度,m/s;
T-空气的绝对温度,K。
雨滴在空气中的运动受到两个力的作用,自身重力和空气的阻力,空气的阻力为:
其中:
式中:μl--空气的动力粘性系数(kg/(m2.s));
ρw—空气的密度,kg/m3;
u—空气的速度,m/s;
uw—雨滴的的速度,m/s。
雨滴的运动方程为:
其余符号同前;
淋水区中的空气流速由风机的风量决定,为已知量,雨滴的流速可由式(2.2-5)求出,雨区的各点的容积散质系数可由式(2.2-1)至式(2.2-3)求得。
有了雨区各点的水滴容积散质系数,可通过下式计算雨区的平均容积散质系数。
式中:ΔV--小的单元体体积,m3;
V--淋水区的体积,m3。
最后雨区的平均容积散质系数,可按式(2-2)求出淋水区的冷却数。
给定不同的雨滴当量直径可以计算出不同条件和淋水高度(模拟进风口高度)的雨区冷却数,根据测试数据,利用计算机的数据拟合程序,分析不同淋水高度雨滴当量直径与冷却数的关系,归纳为下式:
淋水高度为5m时:
淋水高度为10m时:
淋水高度为14.5m时:
对比雨区的测试结果与不同当量直径的热力特性可以找出不同淋水高度、不同运行工况条件下的雨滴当量直径,在计算机内存储为雨滴当量直径表2.2-1。
表2.2-1不同运行工况条件下的雨滴当量直径表
根据所述雨滴当量直径进行运动仿真的步骤,包括:
构建所述雨区模拟模型的塔内气流运动数学模型;
将雨滴的属性设置为相同直径的刚性球,在所述塔内气流运动数学模型中计算所述雨滴的运动方程。
自然通风冷却塔在无自然风时,塔内空气流场为轴对称的二维流动。冷却塔在稳定运行状态下,塔内空气流动为定常流,空气为不可压缩。塔内填料顶面气流的雷诺数可达107,为湍流状态。流动符合轴对称定常不可压二维雷诺时均方程,采用k-ε双方程湍流模型进行方程封闭。
连续方程:
动量方程:
κ方程:
ε方程
式中:
u x方向的气流速度(m/s);
v r方向的气流速度(m/s);
ρ--气流密度(kg/m3);
p--压力(Pa);
g--重力加速度(m/s2);
Fx x方向的填料及水滴阻力(N/m3);
Fr r方向的填料及水滴阻力(N/m3);
μe--有效粘性系数(N/(s·m2));
Cμ、C1、C2、σk、σε--湍流模型常数,其值如下表:
表2.3-1湍流模型常数表
雨滴运动的数学模型
雨滴在冷却塔中的运动是一个十分复杂的过程,计算时将雨滴简化为相同直径的刚性球,在此基础上按牛顿第二定律即可列出雨滴的运动方程。
式中:m--雨滴的质量(kg);
uw,ur--雨滴x、r方向速度(m/s);
t--时间(s);
fx--空气在x方向对雨滴的作用力(N);
fr--空气在r方向对雨滴的作用力(N);
ρ——空气密度(kg/m3);
其余符号同前。
其中:
式中:μl--空气的动力粘性系数(kg/(m2.s));
其余符号同前
雨区的热力特性
在自然通风冷却塔的淋雨区由于风速分布是二维的,沿高度与深度变化,淋水水滴与空气之间的相对速度也是变化的,可以通过式(2.2-1)~式(2.2-3)计算特定雨滴当量直径时雨区各空间点的水滴容积散质系数,采用式(2.1-6)和式(2-2)可求得雨区的冷却数。
在边界条件给定后,利用式(2.3-1)~式(2.3-11)即可进行冷却塔的塔内流场与雨滴流速场的计算。
上式(2.3-2)至式(2.3-5)可以归纳成一个统一的方程形式,设求解参数为φ:
式(2.3-12)中的各参数意义见下表2.3-2。
将式(2.3-12)在图4-3所示的交错网格内积分可以得到式(2.3-13)的离散化方程:
aPφP=aEφE+aWφW+aSφS+aNφN+b (2.3-13)
采用交错网格时将标量放在网格节点上,而将速度分量分别放在w、s面上,如图2.3-3所示。并将方程的源项S线性化为:
S=SC+SPφP (2.3-14)
那么:b=SCr·Δr·Δx
aP=aE+aW+aS+aN-SPr·Δr·Δx
aE=DeA(|Pe|)+[|-Fe,0|]
aW=DwA(|Pw|)+[|Fw,0|]
aS=DsA(|Ps|)+[|-Fs,0|]
aN=DnA(|Pn|)+[|Fn,0|]
Fe=(ρv)ereΔx,Fw=(ρv)wrwΔx,Fs=(ρu)srsΔr,Fn=(ρu)nrnΔr
Pe=Fe/De,Pw=Fw/Dw,Ps=Fs/Ds,Pn=Fn/Dn
A(|P|)=[|0,1-0.5|P||],符号[|A,B|]表示取A、B中较大者。
对于式(2.3-1)至式(2.3-5),按上式(2.3-12)离散方程采用SIMPLE方法迭代求解,迭代采用欠松驰,通过调试松驰因子使迭代收敛,收敛的控制条件是以所有计算点所有方程的余额的和与塔的流入量比小于0.05%。而雨滴运动方程式(2.3-8)、式(2.3-9)系常微分方程,采用四阶龙格库塔法求解。具体求解程序如下:
1)、给定进塔风量;
2)、用SIMPLE法求解空气流场;
3)、用龙格库塔法求解雨滴速度;
4)、判别流场是否收敛,若不收敛返回第2)步;
5)、计算确定雨区平均容积散质系数;
6)、计算雨区冷却数;
7)、输出计算结果。
采用上文的方法和表2.2-1的雨区雨滴当量直径可求出不同淋水面积不同雨区高度的雨区冷却数。取三种淋水面积的冷却塔,对其各种工况进行组合如表2.3-3示。
表2.3-3雨区冷却数计算工况表
将各计算工况进行分析整理,获取所述雨区模拟模型的热力特性数学模型,结果如下:
塔底直径为130m,进风口面积/淋水面积为0.30时,淋水高度为10.25m时:
N=0.329λ0.201 (2.3-15)
塔底直径为130m,进风口面积/淋水面积为0.35时,淋水高度为11.88m时:N=0.369λ0.223 (2.3-16)
塔底直径为130m,进风口面积/淋水面积为0.40时,淋水高度为13.50m时:N=0.409λ0.256 (2.3-17)
塔底直径为150m,进风口面积/淋水面积为0.30时,淋水高度为11.75m时:
N=0.368λ0.226 (2.3-18)
塔底直径为150m,进风口面积/淋水面积为0.35时,淋水高度为13.63m时:
N=0.414λ0.263 (2.3-19)
塔底直径为150m,进风口面积/淋水面积为0.40时,淋水高度为15.50m时:
N=0.456λ0.291 (2.3-20)
塔底直径为170m,进风口面积/淋水面积为0.30时,淋水高度为13.25m时:
N=0.405λ0.260 (2.3-21)
塔底直径为170m,进风口面积/淋水面积为0.35时,淋水高度为15.38m时:
N=0.456λ0.295 (2.3-22)
塔底直径为170m,进风口面积/淋水面积为0.40时,淋水高度为17.50m时:
N=0.474λ0.285 (2.3-23)
根据式(2.3-15)~式(2.3-23),按计算机的数据拟合程序,分析出雨区的热力特性,即雨区冷却数为(底部直径大于110m、进风口高度适合于11~16m):
N=(-0.0014h2+0.061h-0.1496)λ0.0182h+0.0124(2.3-24)
上式即为适用于“底部直径一般大于110m、进风口高度大于11m、淋水面积大于12000m2”的超大型冷却塔雨区热力特性的计算公式。
引用上式雨区测试结果,则式(1-3)变化为:
上式即为“底部直径大于110m、进风口高度11~16m”的超大型冷却塔的拟二维热力计算修正公式。
一种超大型湿式冷却塔雨区热力特性的仿真计算应用方法,包括:
应用如前所述的热力特性数学模型进行冷却塔的选型。
其中,所述热力特性数学模型为:
N=(-0.0014h2+0.061h-0.1496)λ0.0182h+0.0124。
本发明的仿真应用方法,应用所述热力数值计算公式进行冷却塔的建造设计,选型优化,能够节省工程成本,使预期设计工况与实际运行工况相符,误差较小。
图4为本发明超大型湿式冷却塔雨区热力特性的仿真计算模拟装置的示意图,包括:
模型建立模块,用于建立冷却塔的雨区模拟模型;
参数采集模块,用于在预设不同工况的条件下,采集所述雨区模拟模型的雨滴当量直径;
数学模型转换模块,用于根据所述雨滴当量直径进行模拟计算,获取所述雨区模拟模型的热力特性数学模型。
在其中一个实施例当中,所述模型建立模块包括:
模型设置单元,用于设置该模型的底部直径大于110米、进风口高度大于11米,还用于设置所述雨区模拟模型的平面尺寸为1m×1m、高度15m、淋水高度为5m~15m、淋水密度在7-13m3/m2h、风速为0.6m/s~2.8m/s。
在其中一个实施例当中,参数采集模块包括:
测试单元,用于在预设的数个淋水高度下进行测试,采集循环水量、进水温度、出水温度、进气量、进气干球温度、湿球温度、大气压力;
计算单元,用于进行雨滴当量直径的计算;
存储单元,用于将计算的结果存储为雨滴当量直径表,该表包含了淋水高度、冷却数、气水比、风速、淋水密度以及对应的雨滴当量直径。
在其中一个实施例当中,数学模型转换模块包括:
气流仿真单元,用于构建所述雨区模拟模型的塔内气流运动数学模型;
运动仿真单元,用于将雨滴的属性设置为相同直径的刚性球,在所述塔内气流运动数学模型中计算所述雨滴的运动方程。
本模拟仿真装置还包括显示模块,用于显示仿真测试的结果。
图4与图3相对应,图中各个单元的运行方式与方法中的相同。
以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。