CN114970255A - 一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法 - Google Patents

一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法 Download PDF

Info

Publication number
CN114970255A
CN114970255A CN202210538749.6A CN202210538749A CN114970255A CN 114970255 A CN114970255 A CN 114970255A CN 202210538749 A CN202210538749 A CN 202210538749A CN 114970255 A CN114970255 A CN 114970255A
Authority
CN
China
Prior art keywords
fire source
nozzle
height
vertical shaft
shaft
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.)
Pending
Application number
CN202210538749.6A
Other languages
English (en)
Inventor
梅丹
王晨霞
朱禹政
陈旺生
幸福堂
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University of Science and Engineering WUSE
Original Assignee
Wuhan University of Science and Engineering WUSE
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University of Science and Engineering WUSE filed Critical Wuhan University of Science and Engineering WUSE
Priority to CN202210538749.6A priority Critical patent/CN114970255A/zh
Publication of CN114970255A publication Critical patent/CN114970255A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • A—HUMAN NECESSITIES
    • A62—LIFE-SAVING; FIRE-FIGHTING
    • A62C—FIRE-FIGHTING
    • A62C31/00—Delivery of fire-extinguishing material
    • A62C31/02—Nozzles specially adapted for fire-extinguishing
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24—HEATING; RANGES; VENTILATING
    • F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F8/00—Treatment, e.g. purification, of air supplied to human living or working spaces otherwise than by heating, cooling, humidifying or drying
    • F24F8/10—Treatment, e.g. purification, of air supplied to human living or working spaces otherwise than by heating, cooling, humidifying or drying by separation, e.g. by filtering
    • F24F8/117—Treatment, e.g. purification, of air supplied to human living or working spaces otherwise than by heating, cooling, humidifying or drying by separation, e.g. by filtering using wet filtering
    • F24F8/133—Treatment, e.g. purification, of air supplied to human living or working spaces otherwise than by heating, cooling, humidifying or drying by separation, e.g. by filtering using wet filtering by direct contact with liquid, e.g. with sprayed liquid
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Emergency Management (AREA)
  • Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Business, Economics & Management (AREA)
  • Chemical & Material Sciences (AREA)
  • Mechanical Engineering (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Combustion & Propulsion (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法，旨在利用竖井内一侧墙面每层楼交界面高度处安装的喷嘴喷射的最小强制水雾来抑制热气流以及烟气的上升和扩散，阻止由于烟囱效应以及自然风的影响而导致的高层建筑火灾的加剧，同时用水雾裹挟烟颗粒，阻止其对人体的进一步伤害，起到“洗烟”作用；于此同时水雾下降运动，可降低着火部位温度，实现抑烟要求的同时又起到降温效果。

Description

一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法
技术领域
本发明涉及高层建筑竖井内水雾喷射抑烟领域，尤其涉及的是一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法。
背景技术
近年来，随着我国经济的高速发展，城市土地资源变得稀缺，高层建筑的兴建速度也飞速增长，然而相应的高层建筑物的火灾事故发生率也持续增加。高层建筑物中的火势的蔓延的重要途径之一为烟气对流现象，发挥主要作用的一般为建筑物的烟囱效应以及自然界的风力效果，因此有效控制火灾中的烟气传播是阻止火势蔓延办法。
烟囱效应的发生是造成火焰可以在较短的时间内迅速蔓延的重要原因。一旦高层建筑物的低楼层出现火灾险情时，火焰燃烧发出的热空气将会经过楼梯间、电缆管道井、电梯设备井等竖直形式的管道，快速向建筑物的上方蔓延，从而在高层建筑的上部也形成新的起火点，这样不仅会增加消防员抢险救灾工作的难度，还会对建筑物内部被困人员的生命安全带来非常严重的威胁。
因此为了抑制热气流以及烟气的上升和扩散，阻止由于烟囱效应以及自然风的影响而导致的高层建筑火灾的加剧，有必要对高层建筑竖井内的气流流动与传热的机理进行研究，构建一种水雾喷射抑烟模型，实现抑烟要求的同时又可起到降温的效果。
发明内容
本发明的目的在于，提供一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法，通过获取最小喷射速度，实现抑烟要求的同时又起到降温的效果。
本发明解决上述技术问题的技术方案如下：
步骤1：使用Ansys软件设置多种竖井火灾模型工况，所述模型是按高层建筑内竖直形式的管道几何参数和层高1：1等比例设置的y方向高为H，x方向宽为W，z方向长为L 的矩形竖井通道，矩形竖井通道内各楼层高度处设置可喷射水雾的单排喷嘴，喷射方向为Z方向，喷嘴间距以及个数根据竖井尺寸设置，火源顶部距矩形竖井通道底部的距离为C，距火源顶部上方最近的第一排喷嘴至火源顶部的距离为A，喷嘴y方向的尺寸高度为D，火源底部距矩形竖井通道底部的距离为J，所述的模型工况中还包含火源功率Q以及固定的楼层高度、传热系数、导热系数；
竖井通道内紊流流动和传热满足以下方程：
质量守恒方程：
Figure BDA0003647532940000021
式中，u、v和w是速度矢量在三维坐标下x、y和z方向的分量；
动量守恒方程：
x方向，
Figure BDA0003647532940000022
y方向，
Figure BDA0003647532940000023
z方向，
Figure BDA0003647532940000024
式中，μ为粘性系数，μt为湍动能粘性系数，
Figure BDA0003647532940000025
Cμ取值0.0845，k为湍动能，ε为耗散率，ρ为气体密度，g为重力加速度；
能量守恒方程：
Figure BDA0003647532940000026
式中α*＝kf(ρCP)f,kf是烟气的传热系数,(ρCP)f是烟气的密度及热容量；CP：J/(Kg·K)，σT为常数0.9-1.0，Pr是Prandtl数常数0.85；
纲湍动能k方程:
Figure BDA0003647532940000027
式中，αk为湍动能k对应的Prandtl数常数1.39，Gk为平均速度梯度引起的湍动能k的产生项，k是湍动能；
耗散率ε方程
Figure BDA0003647532940000028
Figure BDA0003647532940000029
Figure BDA00036475329400000210
Figure BDA00036475329400000211
Figure BDA00036475329400000212
式中，αε为耗散率ε对应的Prandtl数常数1.39；模型常数C1ε与C2ε分别取值1.42和1.68；η0为 4.377；Eij为时均应变率，i、j为非同方向的x或y或z方向的张量指标符号，ui、uj为非同一坐标方向的时均速度；
Figure BDA00036475329400000213
为
Figure BDA00036475329400000214
中的一个张量，
Figure BDA00036475329400000215
为
Figure BDA00036475329400000216
中的一个张量；β是热膨胀系数，取值为
Figure BDA00036475329400000217
T为烟气平均温度；
按以下边界条件在Ansys软件中设置边界位置处的流动参数和传热参数：
Figure BDA0003647532940000031
qw为热流密度，λ为导热系数，vs为烟气速度，TC为竖井内最低温度，T1为外界空气温度；
步骤2：基于竖井通道高宽比Ar，改变喷嘴喷出的射流的速度v1，使雷诺数Re表征的单排喷射射流产生的惯性力改变，通过Ansys模拟软件进行数据模拟计算，获得单排喷嘴喷射射流时的速度分布、射流方向和竖直方向的流函数图及喷嘴局部速度矢量图，当喷嘴喷出的射流水雾喷到矩形竖井通道内壁上的位置不高于喷嘴的高度时，喷嘴水平界面处热浮力与射流水雾的惯性力处于平衡状态，提取各工况中喷嘴水平界面处的平均速度，得到所述平衡状态时Re的临界值和对应的格拉晓夫数Gr(表征火灾热源产生的热浮力对烟气的作用)；
Ar＝d1/W，
Figure BDA0003647532940000032
ρ为气体密度，d1为喷嘴高度，即最接近火源高度顶层的喷嘴上部至底部面的距离；μ为粘性系数；β为热膨胀系数；g为重力加速度；Δt＝TH-TC：TH和TC分别为竖井内最高温度(火源温度)和竖井内最低温度；Hi为喷嘴高度，喷嘴上部至底部面的距离(Hi＝d1)；
步骤3：以宽高比Ar与所述平衡状态时对应的浮力参数Gr/Re的关系进行拟合，得到水雾喷射抑烟降温模型中Ar与Gr/Re关系：Gr/Re＝A1+B1Ar：A1与B1为拟合常数；
步骤4：改变火源位置，根据Ar与Gr/Re的关系得到最接近火源高度顶层的喷嘴最小的喷射速度。
上述的一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法，所述步骤4之后还包括，改变火源功率，按步骤1～4获得火源在竖井内任一火源位置时，最接近火源高度顶层的喷嘴最小的喷射速度。
一种高层建筑竖井内抑烟降温的水雾喷射装置，其特征在于，所述高层建筑竖井中每层楼安装有所述单排水雾喷射喷嘴，喷嘴通过管道、泵与水源连通，相邻楼层的单排水雾喷射喷嘴之间的竖井壁面设置有连通竖井内外的狭缝，所述竖井各楼层之间还安装有火灾传感器，用于检测火源顶部和底部的位置以及火源功率，并传输给高层建筑内的工控机，工控机内存储有按所述方法获得的不同火源功率和火源位置下的Ar与Gr/Re的关系，工控机根据当前火源功率和火源位置和对应的Ar与Gr/Re的关系查询到最接近火源高度顶层的喷嘴最小的喷射速度，然后控制泵将水以不小于最小喷射速度的速度，从最接近火源高度顶层的喷嘴喷出。
本发明提出的水雾喷射抑烟降温方法，旨在利用最接近建筑竖井内火源楼层上层的第一组喷嘴产生的强制喷射水雾抑制热气流与烟气的上升，阻止由于烟囱效应以及自然风的影响而导致的高层建筑火灾的加剧，同时用水雾裹挟烟颗粒，阻止其对人体的进一步伤害，起到“洗烟”作用；于此同时水雾下降运动，可降低着火部位温度，实现抑烟要求的同时又起到降温效果。
附图说明
图1为竖井喷嘴特征参数模型图；
图2为火源上方距离火源最近楼层中喷嘴喷射射流模型示意图；
图3(a)～3(e)为单排喷射冷却不同Re时速度流函数分布图；
图4(a)～4(e)为单排喷射冷却不同Re时温度分布图；
图5(a)～5(b)为Ar与浮力参数(a)Gr/Re(b)Gr/Re2的拟合线图。
具体实施方式
以下将结合附图，对本发明的具体实施方式和实施例加以详细说明，所描述的具体实施例仅用于解释本发明，并非用于限定本发明的具体实施方式。
如图2所示，图2为水雾喷射抑烟模型图，具体实施步骤如下：
实施例1
现场测量得到竖井以及多排喷嘴结构的相关几何参数，利用Ansys软件建立竖井通道内包含热源、和多排喷嘴的三维几何模型，模型是按高层建筑内竖直形式的管道几何参数和层高1：1等比例设置的，模型是高(y方向)为H，宽(x方向)为W，长(z方向)为L的矩形竖井通道，矩形竖井通道内设置于各楼层的可喷射水雾的一排喷嘴，火源顶部距矩形竖井通道底部的距离为C，距火源顶部上方最近的第一排喷嘴的距离为A，喷嘴尺寸高度为D，火源底部距矩形竖井通道底部的距离为J，v1为距火源顶部上方最近的第一排喷嘴的喷射速度，通过改变喷射水雾速度研究水雾对矩形竖井通道内竖向气流的流动和传热特性的影响，确定临界状态，喷嘴喷射水雾的方向垂直竖直方向，喷嘴底部以下为抑烟降温段。喷嘴入口为速度入口，位于矩形竖井通道底部面上方d1处，矩形竖井通道顶部壁面为压力出口，矩形竖井通道底部为压力进口，热源面为定热流密度加热壁面。模型工况中还包含火源功率Q以及固定的楼层高度、传热系数、导热系数。
步骤1.使用有限体积法建立描述竖井通道内紊流流动和传热的数学物理模型。
1)流动为定常湍流流动，流动的物理量不随时间发生变化；
2)水雾为不可压缩流体，密度满足Boussinesq近似：密度的变化不会使流体的性质发生明显的改变。即可认为，流体的流动过程中粘性耗散忽略不计，仅考虑对质量力项的影响，可理解为在浮力项中密度受温度的影响，除此之外的其他项密度均为常数；
3)喷射入口速度、温度等参数分布均匀。
矩形竖井通道内的流动所需的方程包括基本的控制方程和附加的湍流输运方程。本发明选用“矩形微通道流动换热特性的数值分析”中的湍流模型，湍流模型的具体形式如下：
(1)质量守恒方程：
Figure BDA0003647532940000051
式中，u、v和w是速度矢量在三维坐标下x、y和z方向的分量，单位m/s。
(2)动量守恒方程：
x方向，
Figure BDA0003647532940000052
y方向，
Figure BDA0003647532940000053
z方向，
Figure BDA0003647532940000054
式中，μ为粘性系数，μt为湍动能粘性系数，
Figure BDA0003647532940000055
Cμ取值0.0845，k为湍动能，ε为耗散率，ρ为气体密度，g为重力加速度。
(3)能量守恒方程：
Figure BDA0003647532940000056
式中α*＝kf(ρCP)f,kf是烟气的传热系数,(ρCP)f是烟气的密度及热容量；CP： J/(Kg·K)，σT为常数0.9-1.0，Pr是Prandtl数(普朗特数)常数0.85。
(4)纲湍动能k方程:
Figure BDA0003647532940000057
式中，αk为湍动能k对应的Prandtl数常数1.39，Gk为平均速度梯度引起的湍动能k的产生项，k是湍动能。
(5)耗散率ε方程
Figure BDA0003647532940000061
Figure BDA0003647532940000062
Figure BDA0003647532940000063
Figure BDA0003647532940000064
式中，αε为耗散率ε对应的Prandtl数常数1.39；模型常数C1ε与C2ε分别取值1.42和1.68；η0为4.377；Eij为时均应变率，i、j为非同方向的x或y或z方向的张量指标符号，ui、uj为非同一坐标方向的时均速度；
Figure BDA0003647532940000065
为
Figure BDA0003647532940000066
中的一个张量，
Figure BDA0003647532940000067
为
Figure BDA0003647532940000068
中的一个张量；β是热膨胀系数，取值为
Figure BDA0003647532940000069
T为烟气平均温度；
因此式(1-8)可展开为：
Figure BDA00036475329400000610
Figure BDA00036475329400000611
式(1-7)中的Eij·Eij表示：Exy·Exz或Exy·Eyz或Exz·Eyz；
步骤2.确定计算模型边界条件，喷嘴射流隔热的边界条件设置见表1。
表1计算边界条件
Figure BDA00036475329400000612
qw为热流密度，W·m-2；λ为导热系数，w/(m·K)。
图2为火源上方距离火源最近的楼层中喷嘴喷射射流模型示意图。
步骤3.确定无量纲方程
以雷诺数Re来确定火源上方距离火源最近的楼层中单排喷射射流中的流速，以格拉晓夫数Gr表征火灾热源产生的热浮力对烟气的作用。
Figure BDA0003647532940000071
Figure BDA0003647532940000072
其中：v1为喷嘴喷射射流的速度；d1为喷嘴高度，火源上方距离火源最近的楼层中喷嘴上部至底部面的距离；β为热膨胀系数；Δt＝TH-TC：TH和TC分别为竖井内最高温度(火源最高温度)和竖井内最低温度；H1＝d1。
选取火源上方距离火源最近的楼层中喷嘴第一排喷嘴的高度d1为特征长度，合成速度标尺为特征速度：
公式3Δv＝v1其中：合成温度标尺为特征温度Δt＝TH-TC。使用特征参数进行无量纲化，得公式4：
(X,Y,Z)＝(x,y,z)/d1
(U,V,W)＝(u,v,w)/Δv
P＝p/(ρΔv2)
θ＝(T-TC)/Δt
K＝k/Δv2
E＝ε/(Δv3/d1)
vtn＝vt/(d1Δv)
Figure BDA0003647532940000074
其中：(X,Y,Z)表示笛卡尔坐标，(U,V,W)是X、Y和Z方向上的速度分量，P，θ， K，E，vtn表示无量纲量。P为压强；K是无量纲湍动能，E是无量纲耗散率，vtn是无量纲湍流粘度(vt单位：N·s/m2)，vtn＝(Re CμK2)/E，Cμ为模型常数；Gk为平均速度梯度引起的湍动能k的产生项。
将公式4带入流动模型中的质量、能量、动量方程即可得到无量纲方程：
无量纲质量守恒方程，公式5：
Figure BDA0003647532940000073
无量纲动量守恒方程：
X方向，公式6：
Figure BDA0003647532940000081
Y方向(垂直方向)，公式7：
Figure BDA0003647532940000082
Z方向，公式8：
Figure BDA0003647532940000083
无量纲能量守恒方程，公式9：
Figure BDA0003647532940000084
无量纲湍动能k方程，公式10：
Figure BDA0003647532940000085
无量湍流耗散率ε方程，公式11：
Figure BDA0003647532940000086
其中：
Pr：Prandtl数(普朗特数)常数0.85
σT：常数0.9-1.0
αk：湍动能k对应的Prandtl数常数1.39
αε：耗散率ε对应的Prandtl数常数1.39
C1ε：常数1.42
C2ε：常数1.68
步骤4.确定无量纲边界条件，无量纲化边界条件具体见表2。
表2无量纲化边界条件
Figure BDA0003647532940000091
步骤5：根据步骤3的无量纲化方程和步骤4确定的无量纲化边界条件可以看出，影响计算结果的为烟气物理属性、温度、喷射速度以及竖井尺寸(可由竖井尺寸确定喷嘴位置楼层)(可变物理量)，则该问题中控制参数为Pr(与烟气瞬时物理属性有关)、Gr、Re和 Ar(通道高宽比Ar＝d1/W)，由于烟气的Pr除临界点附近外，几乎与温度及压力无关，可将其看为常数，则在研究中仅考虑Gr和Re对流动与传热的影响，可通过改变喷射速度控制 Re，改变喷嘴高度、温度控制Gr，改变喷嘴高度控制竖井通道高宽比Ar＝d1/W。
步骤6：以步骤5选取的竖井通道高宽比Ar作为变量，改变喷嘴喷出的介质的速度，使用Ansys软件设置多种竖井火灾模型工况。模拟矩形竖井通道内单排喷射射流时的速度分布，获得射流方向和竖直方向的流函数图及喷嘴局部速度矢量图，所述的模型工况中还包含火源位置(即A)、固定火源功率Q以及固定的矩形竖井通道的尺寸。
计算例：矩形竖井通道高20m，宽3.1m，长6.6m(约六层楼高度)，为了让流体充分发展，将火源面设置在楼层高度A＝5m处，火源功率Q设置为15KW，设置火源为稳态火，在模拟过程中放热速率恒定。喷嘴设置在最接近火源面楼层的竖井二楼高度处(约6.6m)， 10个一排均布在竖井宽度方向，喷嘴尺寸高度(y方向)为0.01m。通道入口与出口温度为外界环境温度设置为25℃，通道入口气流流速2m/s，通道出口为压力出口-100pa，其他不加热壁面为绝热壁面。
根据图3(a)～3(e)的流函数图发现，将喷嘴喷出的介质的速度换算为Re后，当喷嘴喷射水雾时会在喷嘴上方形成一个大的涡旋，喷射水雾下方的流场在Re不同时产生差异。Re较小时，浮升力引起的上升气流作用大于喷射水雾产生的惯性力，喷射水雾喷到矩形竖井通道内壁上的位置会高于喷嘴的高度，增加Re会使喷射水雾喷到竖井壁面上的位置降低；当Re 增加到临界值时，在壁面处对应喷嘴高度附近，喷射水雾产生的惯性力和上升热气流达到平衡，喷射水雾产生的惯性力会抑制烟气上升气流的上升并在此处形成一个小的涡旋；当Re 大于临界值时，喷射产生的惯性力作用大于上升气流，喷射惯性力出现向下的逆向流动，同理根据图4(a)～4(e)可见，温度分布存在相同规律。水雾气流由喷嘴高度下降，抑制烟气上升的同时，并到达火源位置进行降温活动。
在Ansys软件设置多种竖井火灾模型工况，基于Ar的改变，继而改变喷嘴喷射水雾射流的速度和雷诺数Re，通过Ansys模拟软件进行数据模拟计算获得单排喷嘴喷射射流时的速度分布、射流方向和竖直方向的流函数图及喷嘴局部速度矢量图，当喷嘴喷出的射流水雾喷到矩形竖井通道内壁上的位置不高于喷嘴的高度时，喷嘴水平界面处热浮力与射流水雾的惯性力处于平衡状态，提取各工况中喷嘴水平界面处的平均速度，得到所述平衡状态时Re 的临界值，和对应的Gr。
步骤:7：在前人如“Mei D,Xing F,Wen M,et al.Numerical simulation ofmixed convection heat transfer of galvanized steel sheets in the verticalalloying furnace[J].Applied Thermal Engineering,2016,93:500-508”和“相梦如.超临界压力CO2对流传热数值研究[D].北京:中国科学院大学,2018.”的研究中浮力参数的形式主要为Gr/Re和Gr/Re2，故以宽高比Ar与临界值对应的浮力参数Gr/Re或Gr/Re2的关系进行拟合。
使用拟合软件以Ar为x轴分别对Gr/Re和Gr/Re2进行拟合，根据其数据点分布样式选择合适的模型进行拟合。通过图5(a)～5(b)可发现Ar与Gr/Re的拟合效果较好，则得到拟合函数：Gr/Re＝A1+B1Ar：A1与B1为拟合常数。
改变火源位置，按步骤1～步骤7重新获得Ar与Gr/Re的关系，根据Ar与Gr/Re的关系得到最接近火源高度顶层的喷嘴最小的喷射速度。
实施例2
一种高层建筑竖井内抑烟降温的水雾喷射装置，高层建筑竖井中每层楼安装有实施例1 中的单排水雾喷射喷嘴，喷嘴通过管道、泵与水源连通，相邻楼层的单排水雾喷射喷嘴之间的竖井壁面设置有连通竖井内外的狭缝，高层建筑竖井各楼层之间还安装有火灾传感器，用于检测火源顶部和底部的位置以及火源功率，并传输给高层建筑内的工控机，工控机内存储有按实施例1所述的方法获得的不同火源功率和火源位置下的Ar与Gr/Re的关系，工控机根据当前火源功率和火源位置和对应的Ar与Gr/Re的关系查询到最接近火源高度顶层的喷嘴最小的喷射速度，然后控制泵将水以不小于最小喷射速度的速度，从最接近火源高度顶层的喷嘴喷出。
虽然结合附图描述了本发明的实施方式，但是本领域技术人员可以在不脱离本发明的精神和范围的情况下做出各种修改和变型，这样的修改和变型均落入由所附权利要求所限定的范围之内。

Claims (3)

1.一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法，其特征在于，包括以下步骤：
步骤1：使用Ansys软件设置多种竖井火灾模型工况，所述模型是按高层建筑内竖直形式的管道几何参数和层高1：1等比例设置的y方向高为H，x方向宽为W，z方向长为L的矩形竖井通道，矩形竖井通道内各楼层高度处设置可喷射水雾的单排喷嘴，喷射方向为Z方向，喷嘴间距以及个数根据竖井尺寸设置，火源顶部距矩形竖井通道底部的距离为C，距火源顶部上方最近的第一排喷嘴至火源顶部的距离为A，喷嘴y方向的尺寸高度为D，火源底部距矩形竖井通道底部的距离为J，所述的模型工况中还包含火源功率Q以及固定的楼层高度、传热系数、导热系数；
竖井通道内紊流流动和传热满足以下方程：
质量守恒方程：
Figure FDA0003647532930000011
式中，u、v和w是速度矢量在三维坐标下x、y和z方向的分量；
动量守恒方程：
x方向，
Figure FDA0003647532930000012
y方向，
Figure FDA0003647532930000013
z方向，
Figure FDA0003647532930000014
式中，μ为粘性系数，μt为湍动能粘性系数，
Figure FDA0003647532930000015
Cμ取值0.0845，k为湍动能，ε为耗散率，ρ为气体密度，g为重力加速度；
能量守恒方程：
Figure FDA0003647532930000016
式中α*＝kf(ρCP)f,kf是烟气的传热系数,(ρCP)f是烟气的密度及热容量；CP：J/(Kg·K)，σT为常数0.9-1.0，Pr是Prandtl数常数0.85；
纲湍动能k方程:
Figure FDA0003647532930000017
式中，αk为湍动能k对应的Prandtl数常数1.39，Gk为平均速度梯度引起的湍动能k的产生项，k是湍动能；
耗散率ε方程
Figure FDA0003647532930000018
Figure FDA0003647532930000019
Figure FDA00036475329300000110
Figure FDA0003647532930000021
Figure FDA0003647532930000022
式中，αε为耗散率ε对应的Prandtl数常数1.39；模型常数C1ε与C2ε分别取值1.42和1.68；η0为4.377；Eij为时均应变率，i、j为非同方向的x或y或z方向的张量指标符号，ui、uj为非同一坐标方向的时均速度；
Figure FDA0003647532930000023
为
Figure FDA0003647532930000024
中的一个张量，
Figure FDA0003647532930000025
为
Figure FDA0003647532930000026
中的一个张量；β是热膨胀系数，取值为
Figure FDA0003647532930000027
T为烟气平均温度；
按以下边界条件在Ansys软件中设置边界位置处的流动参数和传热参数：
Figure FDA0003647532930000028
qw为热流密度，λ为导热系数，vs为烟气速度，TC为竖井内最低温度，T1为外界空气温度；
步骤2：基于竖井通道高宽比Ar，改变喷嘴喷出的射流的速度v1，使雷诺数Re表征的单排喷射射流产生的惯性力改变，通过Ansys模拟软件进行数据模拟计算，获得单排喷嘴喷射射流时的速度分布、射流方向和竖直方向的流函数图及喷嘴局部速度矢量图，当喷嘴喷出的射流水雾喷到矩形竖井通道内壁上的位置不高于喷嘴的高度时，喷嘴水平界面处热浮力与射流水雾的惯性力处于平衡状态，提取各工况中喷嘴水平界面处的平均速度，得到所述平衡状态时Re的临界值，和对应的格拉晓夫数Gr；
Ar＝d1/W，
Figure FDA0003647532930000029
ρ为气体密度，d1为喷嘴高度，即最接近火源高度顶层的喷嘴上部至底部面的距离；μ为粘性系数；β为热膨胀系数；g为重力加速度；Δt＝TH-TC：TH和TC分别为竖井内最高温度和竖井内最低温度；Hi为喷嘴高度，喷嘴上部至底部面的距离Hi＝d1；
步骤3：以宽高比Ar与所述平衡状态时对应的浮力参数Gr/Re的关系进行拟合，得到水雾喷射抑烟降温模型中Ar与Gr/Re关系：Gr/Re＝A1+B1Ar：A1与B1为拟合常数；
步骤4：改变火源位置，按步骤1～步骤3重新获得Ar与Gr/Re的关系，根据Ar与Gr/Re的关系得到最接近火源高度顶层的喷嘴最小的喷射速度。
2.根据权利要求1所述的方法，其特征在于，所述步骤4之后还包括，改变火源功率，按步骤1～4获得火源在竖井内任一火源位置时，最接近火源高度顶层的喷嘴最小的喷射速度。
3.基于权利要求1～2任一项所述方法的一种高层建筑竖井内抑烟降温的水雾喷射装置，其特征在于，所述高层建筑竖井中每层楼安装有所述单排水雾喷射喷嘴，喷嘴通过管道、泵与水源连通，相邻楼层的单排水雾喷射喷嘴之间的竖井壁面设置有连通竖井内外的狭缝，所述竖井各楼层之间还安装有火灾传感器，用于检测火源顶部和底部的位置以及火源功率，并传输给高层建筑内的工控机，工控机内存储有按所述方法获得的不同火源功率和火源位置下的Ar与Gr/Re的关系，工控机根据当前火源功率和火源位置和对应的Ar与Gr/Re的关系查询到最接近火源高度顶层的喷嘴最小的喷射速度，然后控制泵将水以不小于最小喷射速度的速度，从最接近火源高度顶层的喷嘴喷出。
CN202210538749.6A 2022-05-17 2022-05-17 一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法 Pending CN114970255A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210538749.6A CN114970255A (zh) 2022-05-17 2022-05-17 一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210538749.6A CN114970255A (zh) 2022-05-17 2022-05-17 一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法

Publications (1)

Publication Number Publication Date
CN114970255A true CN114970255A (zh) 2022-08-30

Family

ID=82983655

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210538749.6A Pending CN114970255A (zh) 2022-05-17 2022-05-17 一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法

Country Status (1)

Country Link
CN (1) CN114970255A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116050289A (zh) * 2022-12-12 2023-05-02 宁波大学 一种计算城市排水竖井射流冲击力的方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116050289A (zh) * 2022-12-12 2023-05-02 宁波大学 一种计算城市排水竖井射流冲击力的方法
CN116050289B (zh) * 2022-12-12 2023-10-10 宁波大学 一种计算城市排水竖井射流冲击力的方法

Similar Documents

Publication Publication Date Title
Ji et al. Experimental investigation on the rising characteristics of the fire-induced buoyant plume in stairwells
Mousavi et al. Heat and fluid flow in entrance region of a channel with staggered baffles
Chochua et al. A computational and experimental investigation of turbulent jet and crossflow interaction
CN114970255A (zh) 一种高层建筑竖井内抑烟降温的水雾喷射速率确定方法
Chow et al. Numerical modeling for interaction of a water spray with smoke layer
CN104729824B (zh) 一种用于冷却高马赫数喷管喉道的换热装置及其构造方法
Zhou et al. The effect of forced ventilation by using two movable fans on thermal smoke movement in a tunnel fire
Ji et al. Ambient wind effect on combustion characteristics in compartment with simultaneous door and window opened
Lin et al. Buoyancy-driven flows by a heat source at different levels
Wang et al. Reduced-scale experimental and numerical study of fire in a hybrid ventilation system in a large underground subway depot with superstructures under fire scenario
Velayati et al. Numerical study of convective heat transfer from an array of parallel bluff plates
He et al. Study on the effect of tunnel dimensions on the smoke layer thickness in naturally ventilated short tunnel fires
Cao et al. Numerical study on the effect of buoyancy-driven pollution source on vortex ventilation performance
Kinell et al. Comparison of gas turbine vane pressure side and suction side film cooling performance and the applicability of superposition
Yang et al. Experimental and numerical investigations of vane endwall film cooling with different density ratios
Paul et al. A multiphase computational framework for deposit formation and growth
Sugawa Simple estimation model on ceiling temperature and velocity of fire induced flow under ceiling
Wan et al. Experimental and numerical studies on the spray cooling of deionized water in non-boiling regime
Cho et al. Comparison between wall confluent jets and displacement ventilation in aspect of the spreading ratio on the floor
Tkachenko et al. Computational study of natural convection flow in an open-ended channel coupled with a room: Application to building-integrated photovoltaic (BIPV) systems
Hara et al. Numerical Simulation of Fire Plume-Induced Ceiling Jets Using the Standard k ε Model
Hadjisophocleous et al. Evaluation of Atrium Smoke Exhaust Make-Up Air Velocity.
Hu et al. A Comprehensive CFD study of tile flow rate distribution in a compact data center laboratory
Mei et al. Energy conservation and heat transfer enhancement for mixed convection on the vertical galvanizing furnace
Phan et al. Air Flow Velocity Field Validation and Turbulence Studies on a Single Rack Model in Data Centers

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