CN109214053A - 一种外热式回转窑污泥热解***的传热模型计算方法 - Google Patents

Abstract

本发明公开了一种外热式回转窑污泥热解***的传热模型计算方法，包含以下步骤：S1、计算烟气通道内的对流换热系数h₁；S2、计算烟气通道内的辐射传热系数h₂；S3、计算回转窑内壁面的污泥颗粒上表面的对流换热系数h₃和辐射传热系数h₄；S4、计算回转窑内壁与污泥颗粒之间的辐射传热系数h₅；S5、计算回转窑内壁面污泥颗粒下表面的导热系数h₆；S6、计算外热式回转窑污泥热解***的传热量ΔQ。本发明切实可行的提出了一种外热式回转窑污泥热解***的传热模型计算方法，对外热式回转窑应用于污泥热解的设计改造、燃料利用率均有重要意义。

Description

一种外热式回转窑污泥热解***的传热模型计算方法

技术领域

本发明涉及回转窑技术领域，具体涉及一种外热式回转窑污泥热解***的传热模型计算方法。

背景技术

回转窑是一种广泛应用于食品、建材、冶金、化工、环保等行业，对固体物料进行机械、物理或化学处理的回转圆筒类热工设备。根据对窑内物料加热方式的不同，回转窑可分外热式和内热式。

目前，外热式回转窑正在诸如食品干燥等加工环节、催化剂焙烧等工业生产以及固体废物和塑料的热解处理中获得越来越广泛的应用。市政污泥作为城市的主要生活垃圾之一，利用外热式回转窑进行热解的工艺也开始投入实际生产。

但是，现有技术中在利用外热式回转窑进行污泥热解的运行过程中，回转窑内部存在着热效率较低的缺陷，导致污泥热解所需燃料量较高且污泥热解不充分等问题。

发明内容

为了解决现有技术中的上述不足，有必要针对上述外热式回转窑内污泥热解进行传热量计算，提供一种外热式回转窑污泥热解***的传热模型计算方法，利用此计算方法分析污泥颗粒在外热式回转窑内进行热解的热量传递影响因素及变化规律，对外热式回转窑应用于污泥热解的设计改造、燃料利用率均有重要意义。

本发明的目的可以通过采取如下技术方案达到：

一种外热式回转窑污泥热解***的传热模型计算方法，所述的外热式回转窑污泥热解***由双通道组成，其中，外通道为环形烟气通道，作为热源；内通道为圆柱形污泥通道，作为热解对象，整个热解***的传热过程分解为：窑外换热、窑壁换热和窑内换热，所述的传热模型计算方法包括以下步骤，

S1、计算烟气通道内的对流换热系数h₁；

S2、计算烟气通道内的辐射传热系数h₂；

S3、计算回转窑内壁面的污泥颗粒上表面的对流换热系数h₃和辐射传热系数h₄；

S4、计算回转窑内壁与污泥颗粒之间的辐射传热系数h₅；

S5、计算回转窑内壁面污泥颗粒下表面的导热系数h₆；

S6、计算外热式回转窑污泥热解***的传热量ΔQ。

进一步地，所述的烟气通道内的对流换热系数h₁的计算公式如下，

h₁＝d/Nu₁λ_g (1)

式中，d为当量直径，λ_g为烟气通道内换热系数，Nu₁为烟气通道内烟气的努塞尔数，其值由下式计算，

Nu₁＝0.023Re^0.8Pr^0.3 (2)

式中，Re为雷诺数，Pr为普朗特数。

进一步地，所述的烟气通道内热辐射的传热系数h₂的计算公式如下，

式中，σ为黑体辐射常数，T_g为烟气温度，T_o为回转窑外壁温度，α_g为烟气平均吸收比，ε_g为烟气平均发射率，α_g和ε_g由下式决定，

式中，C_H2O为H₂O气体修正系数，C_CO2为CO₂气体修正系数，ε_H2O为H₂O气体发射率，ε_CO2为CO₂气体发射率，α_H2O为H₂O气体吸收比，α_CO2为CO₂气体吸收比，Δε为H₂O气体与CO₂气体共存时的发射率修正量，Δα为H₂O气体与CO₂气体共存时的修正量，其中，α_H2O和α_CO2通过下式确定，

式中，P_H2O为H₂O气体分压力，P_CO2为CO₂气体分压力，中括号[]的下角码为确定方括号内的量时所用的参量。

进一步地，所述的回转窑内壁面污泥上表面的对流换热系数h₃的计算公式如下，

h₃＝d/Nu₂λ₁ (6)

式中，λ₁为回转窑壁面导热系数，Nu₂为窑内壁面污泥上表面热解气体的努塞尔数，其值可由下式确定，

Nu₂＝0.664Re^0.5Pr^0.33 (7)

式中，Re为雷诺数，Pr为普朗特数。

进一步地，所述的回转窑内壁面污泥颗粒上表面的辐射传热系数h₄为，

式中，T_s1为污泥颗粒上表面温度，T_i为回转窑内壁表面温度，ε_s为污泥颗粒上表面处的***发射率，其值可通过下式确定：

式中，ε_s1为污泥颗粒上表面发射率，ε_i为回转窑内壁发射率，A_s1为污泥颗粒上表面面积，A_i为窑壁表面面积。

进一步地，所述的回转窑内壁与污泥颗粒之间的辐射传热系数h₅的计算公式如下，

式中，T_s2为污泥颗粒下表面温度，ε'_s为污泥颗粒下表面处的***发射率，其值可由下式确定，

式中，ε_s2为污泥颗粒下表面发射率，ε_i为回转窑内壁发射率。

进一步地，所述的回转窑内壁面污泥颗粒下表面的导热系数h₆的计算公式如下，

式中，h₆₁为回转窑内壁与气膜层的对流换热系数，h₆₂为污泥颗粒与气膜层之间对流换热系数，其中，h₆₁和h₆₂根据下式确定，

式中，λ_g’为气膜层导热系数，r_s为颗粒半径，n为回转窑转速，d_i为回转窑内径，d_s为污泥颗粒直径，v为气膜层粘度。

进一步地，所述的外热式回转窑污泥热解***的传热量Q的计算公式如下，

Q＝Q₁+Q₂+Q₃+Q₄ (14)

式中，Q₁为烟气通道总传热量，Q₂为回转窑内壁污泥上表面处的传热量，Q₃为回转窑内壁污泥下表面处的辐射传热量，Q₄为回转窑内壁污泥下表面处的导热传热量，其中，Q₁-Q₄通过下列确定，

式中，d_o为回转窑外壁直径，l_w为回转窑窑长，l_s为回转窑内壁污泥颗粒截面弦长，l_arc为回转窑内壁与单个污泥颗粒接触的平均弧长。

本发明相对于现有技术具有如下的优点及效果：

本发明提供了一种外热式回转窑污泥热解***的传热模型计算方法，利用此传热模型分析污泥颗粒在外热式回转窑内进行热解的热量传递影响因素及变化规律，对外热式回转窑应用于污泥热解的设计改造、燃料利用率均有重要意义。

附图说明

图1是本发明中外热式回转窑污泥热解***的结构图；

图2是本发明中外热式回转窑截面图；

图3是本发明中回转窑内壁污泥上表面换热模型；

图4是本发明中回转窑内壁污泥下表面换热模型；

图5是本发明中公开的一种外热式回转窑污泥热解***的传热模型计算方法的流程图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例

本实施例公开一种外热式回转窑污泥热解***的传热模型计算方法，如图1所示为一种典型的外热式回转窑污泥热解***的结构，本发明的计算模型及计算方法适用于此类回转窑。如图2所示为外热式回转窑截面图，图中清楚阐述了烟气通道、窑壁和污泥之间的对在空间上的对应关系，本发明在求解回转窑污泥热解***传热量的过程中，所有使用的参数和建立的公式均是围绕着此对应关系展开，其具体过程包括下列步骤：

1)计算烟气通道内的对流换热系数h₁，参照下列计算公式，

h₁＝d/Nu₁λ_g (1)

式中，d为当量直径，λ_g为烟气通道内换热系数，Nu₁为烟气通道内烟气的努塞尔数，其值可由下式计算，

Nu₁＝0.023Re^0.8Pr^0.3 (2)

式中，Re为雷诺数，Pr为普朗特数。

对于不含水分蒸发的污泥热解段：取烟气平均温度550℃，烟气平均流速0.85m/s，烟气密度0.431kg/Nm³，烟气运动黏度8.50×10^-5m²/s，烟气的普朗特数Pr为0.61，烟气的导热系数0.0699W/(m·K)，通过计算得到窑外壁面对流换热系数h₁为3.55W/(m²·K)。

对于含水分蒸发的污泥热解段：取烟气平均温度375℃，烟气平均流速0.67m/s，烟气密度0.548kg/Nm³，烟气运动黏度5.67×10^-5m²/s，烟气的普朗特数Pr为0.65，烟气的导热系数0.05W/(m·K)，通过计算得到窑外壁面对流换热系数h₁为4.07W/(m²·K)。

2)计算烟气通道内的辐射传热系数h₂，参照下列计算公式，

式中，σ为黑体辐射常数，T_g为烟气温度，T_o为回转窑外壁温度，α_g为烟气平均吸收比，ε_g为烟气平均发射率，α_g和ε_g可由下式决定，

式中，C_H2O为H₂O气体修正系数，C_CO2为CO₂气体修正系数，ε_H2O为H₂O气体发射率，ε_CO2为CO₂气体发射率，α_H2O为H₂O气体吸收比，α_CO2为CO₂气体吸收比，Δε为H₂O气体与CO₂气体共存时的发射率修正量，Δα为H₂O气体与CO₂气体共存时的修正量，其中，α_H2O和α_CO2可通过下式确定，

式中，P_H2O为H₂O气体分压力，P_CO2为CO₂气体分压力，中括号[]的下角码为确定方括号内的量时所用的参量。

烟气成分按标准烟气计算，其中二氧化碳占13％，水蒸气占11％，氮气占76％。烟气总压力为10⁵Pa。O₂、H₂和N₂等分子结构对称的双原子气体无发射与吸收辐射能的能力，为辐射透明体，不参与辐射换热。根据烟气温度T_g为843.15K，可得，C_H2O为1.11，C_CO2为0.98，ε_H2O为0.095，ε_CO2为0.098，Δε为0.005，最终计算可得ε_g为0.196。T_g和T_o分别为843.15K和763.15K，通过可得射率；α_H2O为0.102，α_CO2为0.098，Δα为0.004，最终计算可得α_g为0.206。由此可得，污泥热解段烟气侧的总传热系数h₁+h₂为20.68W/(m²·K)，水分蒸发段烟气侧的总传热系数h₁+h₂为12.30W/(m²·K)。

3)计算回转窑内壁面的污泥颗粒上表面的对流换热系数h₃，参照下列计算公式，

h₃＝d/Nu₂λ₁ (6)

式中，λ₁为回转窑壁面导热系数，Nu₂为窑内壁面污泥上表面的努塞尔数，其值可由下式确定，

Nu₂＝0.664Re^0.5Pr^0.33 (7)

对于不含水分蒸发的污泥热解段取热解气温度431.65℃，热解气体流速0.38m/s，求得换热系数h₃为1.41W/(m²·K)，含水分蒸发的污泥热解段不存在热解气体。

4)计算回转窑内壁面的污泥颗粒上表面的辐射传热系数h₄，参照下列计算公式，

式中，T_s1为污泥颗粒上表面温度，T_i为回转窑内壁表面温度，ε_s为污泥颗粒上表面处的***发射率,其值可通过下列确定：

式中，ε_s1为污泥颗粒上表面发射率，ε_i为回转窑内壁发射率，A_s1为污泥颗粒上表面面积，A_i为窑壁表面面积；

ε_s1取黑体发射率1，ε_i取钢材发射率0.52，计算可得ε_s为0.819。对于污泥热解段，取窑壁温度为366℃，污泥上表面温度260℃，计算得到h₄为34.31W/(m²·K)。对于水分蒸发段，取窑壁温度为222℃，污泥上表面温度100℃，计算得到h₄为14.10W/(m²·K)。因此，对于不含水分蒸发的污泥热解段通过污泥上表面的总换热系数h₃+h₄为35.72W/(m²·K)；对于含水分蒸发的污泥热解段通过污泥上表面的总换热系数h₃+h₄为14.10W/(m²·K)。

5)计算回转窑内壁与污泥颗粒之间的辐射传热系数h₅，参照下列计算公式，

式中，T_s2为污泥颗粒下表面温度，ε'_s为污泥颗粒下表面处的***发射率，其值可由下式确定，

式中，ε_s2为污泥颗粒下表面发射率，ε_i为回转窑内壁发射率。

ε_s2取黑体发射率1，ε_i取钢材发射率0.52，计算可得ε_s’为0.52。对于不含水分蒸发的污泥热解段，窑内壁面与外壁面温度相等(366℃)，颗粒表面平均温度取260℃，求得回转窑内壁与污泥颗粒之间的辐射传热系数h₅为23.94W/(m²·K)；对于含水分蒸发的污泥热解段，窑内壁面与外壁面温度相等(222℃)，颗粒表面平均温度取338.50℃，求得回转窑内壁与污泥颗粒之间的辐射传热系数h₅为9.7W/(m²·K)。

6)计算回转窑内壁面污泥颗粒下表面的导热系数h₆，参照下列计算公式，

式中，h₆₁为回转窑内壁与气膜层的对流换热系数；h₆₂为污泥颗粒与气膜层之间对流换热系数。其中，h₆₁和h₆₂可根据下式确定，

式中，λ_g’为气膜层导热系数；r_s为颗粒半径；n为回转窑转速；d_i为回转窑内径；d_s为污泥颗粒直径；v为气膜层粘度。

d_i为0.5m，颗粒直径为0.006m，回转窑转速0.67r/min。对于不含水分蒸发的污泥热解段，窑内气膜温度与内壁一致(366℃)，计算可得回转窑内壁与气膜层的对流换热系数h₆₁为30.6W/(m²·K)，污泥颗粒与气膜层之间对流换热系数h₆₂为6.96W/(m²·K)，回转窑内壁面污泥颗粒下表面的导热系数h₆为5.67W/(m²·K)。对于含水分蒸发的污泥热解段，窑内气膜温度与内壁一致(222℃)，计算可得回转窑内壁与气膜层的对流换热系数h₆₁为20.88W/(m²·K)，污泥颗粒与气膜层之间对流换热系数h₆₂为4.75W/(m²·K)，回转窑内壁面污泥颗粒下表面的导热系数h₆为3.87W/(m²·K)。

7)计算外热式回转窑污泥热解***的传热量ΔQ，参照下列计算公式，

ΔQ＝-Q₁+Q₂+Q₃+Q₄ (14)

式中，Q₁为烟气通道总传热量，Q₂为回转窑内壁污泥上表面处的传热量，Q₃为回转窑内壁污泥下表面处的辐射传热量，Q₄为回转窑内壁污泥下表面处的导热传热量。其中，Q₁-Q₄可通过下列确定，

式中，d_o为回转窑外壁直径，l_w为回转窑窑长，l_s为回转窑内壁污泥颗粒截面弦长，l_arc为回转窑内壁与单个污泥颗粒接触的平均弧长。

回转窑外壁直径为0.506m，回转窑窑长为3.45m，回转窑内壁污泥颗粒截面弦长为0.812m，所有颗粒弧长l_arc总和为0.9476m。计算可得，烟气通道总传热量Q₁为2.09×10⁴J，回转窑内壁污泥上表面处的传热量Q₂为1.06×10⁴J，回转窑内壁污泥下表面处的辐射传热量Q₃为8297.72J，回转窑内壁污泥下表面处的导热传热量Q₄为1965.39J。由此可得，外热式回转窑污泥热解***的传热量ΔQ为-36.89J，结果符合能量守恒定律，因此传热量Q₁-Q₄计算结果可靠。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (8)

1.一种外热式回转窑污泥热解***的传热模型计算方法，所述的外热式回转窑污泥热解***由双通道组成，其中，外通道为环形烟气通道，作为热源；内通道为圆柱形污泥通道，作为热解对象，整个热解***的传热过程分解为：窑外换热、窑壁换热和窑内换热，其特征在于，所述的传热模型计算方法包括以下步骤，

S1、计算烟气通道内的对流换热系数h₁；

S2、计算烟气通道内的辐射传热系数h₂；

S3、计算回转窑内壁面的污泥颗粒上表面的对流换热系数h₃和辐射传热系数h₄；

S4、计算回转窑内壁与污泥颗粒之间的辐射传热系数h₅；

S5、计算回转窑内壁面污泥颗粒下表面的导热系数h₆；

S6、计算外热式回转窑污泥热解***的传热量ΔQ。

2.根据权利要求1所述的一种外热式回转窑污泥热解***的传热模型计算方法，其特征在于，所述的烟气通道内的对流换热系数h₁的计算公式如下，

h₁＝d/Nu₁λ_g (1)

式中，d为当量直径，λ_g为烟气通道内换热系数，Nu₁为烟气通道内烟气的努塞尔数，其值由下式计算，

Nu₁＝0.023Re^0.8Pr^0.3 (2)

式中，Re为雷诺数，Pr为普朗特数。

3.根据权利要求1所述的一种外热式回转窑污泥热解***的传热模型计算方法，其特征在于，所述的烟气通道内热辐射的传热系数h₂的计算公式如下，

式中，σ为黑体辐射常数，T_g为烟气温度，T_o为回转窑外壁温度，α_g为烟气平均吸收比，ε_g为烟气平均发射率，α_g和ε_g由下式决定，

式中，为H₂O气体修正系数，为CO₂气体修正系数，为H₂O气体发射率，为CO₂气体发射率，为H₂O气体吸收比，为CO₂气体吸收比，Δε为H₂O气体与CO₂气体共存时的发射率修正量，Δα为H₂O气体与CO₂气体共存时的修正量，其中，和通过下式确定，

式中，为H₂O气体分压力，为CO₂气体分压力，中括号[]的下角码为确定方括号内的量时所用的参量。

4.根据权利要求1所述的一种外热式回转窑污泥热解***的传热模型计算方法，其特征在于，所述的回转窑内壁面污泥上表面的对流换热系数h₃的计算公式如下，

h₃＝d/Nu₂λ₁ (6)

式中，λ₁为回转窑壁面导热系数，Nu₂为窑内壁面污泥上表面热解气体的努塞尔数，其值可由下式确定，

Nu₂＝0.664Re^0.5Pr^0.33 (7)

式中，Re为雷诺数，Pr为普朗特数。

5.根据权利要求1所述的一种外热式回转窑污泥热解***的传热模型计算方法，其特征在于，所述的回转窑内壁面污泥颗粒上表面的辐射传热系数h₄的计算公式如下，

式中，T_s1为污泥颗粒上表面温度，T_i为回转窑内壁表面温度，ε_s为污泥颗粒上表面处的***发射率，其值通过下式确定：

式中，ε_s1为污泥颗粒上表面发射率，ε_i为回转窑内壁发射率，A_s1为污泥颗粒上表面面积，A_i为窑壁表面面积。

6.根据权利要求1所述的一种外热式回转窑污泥热解***的传热模型计算方法，其特征在于，所述的回转窑内壁与污泥颗粒之间的辐射传热系数h₅的计算公式如下，

式中，T_s2为污泥颗粒下表面温度，ε'_s为污泥颗粒下表面处的***发射率，其值由下式确定，

式中，ε_s2为污泥颗粒下表面发射率，ε_i为回转窑内壁发射率。

7.根据权利要求1所述的一种外热式回转窑污泥热解***的传热模型计算方法，其特征在于，所述的回转窑内壁面污泥颗粒下表面的导热系数h₆的计算公式如下，

式中，h₆₁为回转窑内壁与气膜层的对流换热系数，h₆₂为污泥颗粒与气膜层之间对流换热系数，其中，h₆₁和h₆₂根据下式确定，

式中，λ_g’为气膜层导热系数，r_s为颗粒半径，n为回转窑转速，d_i为回转窑内径，d_s为污泥颗粒直径，v为气膜层粘度。

8.根据权利要求1所述的一种外热式回转窑污泥热解***的传热模型计算方法，其特征在于，所述的外热式回转窑污泥热解***的传热量Q的计算公式如下，

Q＝Q₁+Q₂+Q₃+Q₄ (14)

式中，Q₁为烟气通道总传热量，Q₂为回转窑内壁污泥上表面处的传热量，Q₃为回转窑内壁污泥下表面处的辐射传热量，Q₄为回转窑内壁污泥下表面处的导热传热量，其中，Q₁、Q₂、Q₃、Q₄通过下列确定，

式中，d_o为回转窑外壁直径，l_w为回转窑窑长，l_s为回转窑内壁污泥颗粒截面弦长，l_arc为回转窑内壁与单个污泥颗粒接触的平均弧长。