工业炉炉壁面散热模型及热环境分析.docxVIP

工业炉炉壁面散热模型及热环境分析 行业渣的外壁通常占炉渣总热损失的40%左右。其次，它在炉渣热损失项目中排名第二，占相当大的数量。据不完全统计,2005年,全国县以上企业工业炉窑有10多万台,其能耗为全国总能耗的1/4,约占工业总能耗的60%,以我国2005年能源消费总量22.2 亿t标准煤、散热损失以20%计,则我国全年通过炉壁散热向大气环境排放废热约1.11 亿t标准煤,数量巨大,是造成城市热岛效应(urban heat island)的主要元凶之一,造成我国城市热岛效应非常突出:北京1961～2000年市郊日均气温差为3.3 ℃, 上海热岛强度最高达7.4 ℃,并且伴随着工业化进程的继续,城市热岛效应强度及规模还在日益加剧。本文通过对工业炉壁面的散热建立数学模型,探索可行的壁面热量回收方法,以此减少工业炉炉壁对大气环境的热污染,实现节能减排。 1 工业设备壁热环境分布规律 1.1 炉壁热风蒸焦 工业炉设备运行时燃料燃烧产生的热量一部分加热物件,称为有效热;其余部分则为无效热(或称废热),以各种形态、不同方式散发到周围的环境中去。炉壁散热将炉体周围的冷空气加热成热空气而形成烟囱效应,最后排放到大气中去,炉壁温度愈高,则散发到空气中的热量愈多,工业炉的热利用率愈低,对环境造成的热污染愈严重。 炉体壁面散热特点:外表散热面积大、热量散发分散、散发的热量品位低,如受到外界干扰,则其流场会发生很大的影响。如图1所示为工业炉炉壁散热对周围环境造成热污染的示意图。 1.2 自然对流控制方程 以工业炉壁面为研究对象,图2为工业炉物理模型及坐标系,由于工业炉属大型加热设备,描述其壁面外空气的运动,可以不考虑三维效应,而只采用二维稳定流动控制方程。设壁面温度tw,空气环境温度t∞,近壁面处流体受热而沿壁面作上升运动。采用Boussinesq假定的工业炉壁面二维稳定流动的控制方程组如下: 连续性方程: ?u?x+?v?y=0(1) 动量方程: x方向:ρu?u?x+ρv?u?y=-?p?x+μ(?2u?x2+?2u?y2)-ρg y方向:ρu?v?x+ρv?v?y=-?p?y+μ(?2v?x2+?2v?y2)(2) 能量方程: ρucp?t?x+ρvcp?t?y-u?p?x-v?p?y=??x(λ?t?x)+??y(λ?t?y)(3) 根据流体在炉壁表面流动的物理特点,采用尺度分析法判断方程中各项的数量级别,忽略小数量级项,简化上述方程组。如图2所示,炉体高度为L,边界层厚度为δ,对于变量x,y和u规定以下的尺度:x类比于L,y类比于δ和u类比于u∞,并且δL。 由尺度分析法判断方程中各项的数量级别可知:动量方程中y方向的各项中均相应地减少了δ/L数量级,因此y方向动量方程可不预考虑。同时由于自然对流边界层外的流体是静止的,压力只是x方向的函数,则:?p/?y=0;?p/?x=dp/dx=dp∞/dx=-ρ∞g。 由流体性质知:动量扩散率(运动粘度)为υ=μ/ρ;热扩散率为a=λ/ρcp,空气体膨胀系数β=1/T∞。则二维控制方程组简化为: ?u?x+?v?y=0(4) u?u?x+v?u?y=υ?2u?y2-(t∞-t)βg(5) u?t?x+v?t?y=a?2t?y2(6) 定壁温方程边界条件: {u=v=0t=tw(当y=0时)u→0t=t∞(当y=∞时)u=0t=t∞(当x=0时) 将上述方程组无因次化。定义无因次参数为: X=x/L U=u/(gβΔtL)1/2 Y=y/L V=v/(gβΔtL)1/2 T=(t-t∞)/(tw-t∞) 则无因次控制方程组为: ?U?X+?V?Y=0(7) U?U?X+V?U?Y=1Gr1/2t?2U?Y2+Τ(8) U?Τ?X+V?Τ?Y=1Gr1/2tΡr?2Τ?Y2(9) 式中:Grt为格拉肖夫(Grashof)数,Grt=gβΔtL3υ3,Δt=tw-t∞。 无因次边界条件为: Y=0,0≤X≤1:T=1,U=V=0 Y=0.4? :Τin=0,V=0,?Τ?X|out=0,?U?X|out=0 图3为数值边界条件示意图。根据经验,自然对流主要是在炉壁面形成,它相对于炉壁高度来说比较小。所以,当Y大于0.4时,认为速度、温度值已不再发生变化,故Y方向的计算宽度取为0.4。 1.3 压力场的离散化 SIMPIE方法是求解压力关联的半隐式方法(semi-implicit method for pressure-linked equation)的简称。算法的关键是利用连续性方程和动量方程来构成一个近似的压力校正方程来求解速度场。 采用控制容积法对方程进行离散化,导入交错网格对计算区域进行划分,并可在不同的网格点上求解不同的变量,如图4所示,E、W、S、N为网格东西南北各方向。压力P在控制体(i,