[go: up one dir, main page]

CN107301276B - Method for calculating convection heat transfer load of large-space nozzle air supply layered air conditioner - Google Patents

Method for calculating convection heat transfer load of large-space nozzle air supply layered air conditioner Download PDF

Info

Publication number
CN107301276B
CN107301276B CN201710403639.8A CN201710403639A CN107301276B CN 107301276 B CN107301276 B CN 107301276B CN 201710403639 A CN201710403639 A CN 201710403639A CN 107301276 B CN107301276 B CN 107301276B
Authority
CN
China
Prior art keywords
air
heat transfer
area
model
indoor
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.)
Active
Application number
CN201710403639.8A
Other languages
Chinese (zh)
Other versions
CN107301276A (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.)
Shanghai Shinelong Air Conditioning Co ltd
Original Assignee
University of Shanghai for Science and Technology
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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN201710403639.8A priority Critical patent/CN107301276B/en
Publication of CN107301276A publication Critical patent/CN107301276A/en
Application granted granted Critical
Publication of CN107301276B publication Critical patent/CN107301276B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Air Conditioning Control Device (AREA)

Abstract

The invention relates to a method for calculating convection heat transfer load of a large-space-nozzle air supply layered air conditioner, which comprises the steps of establishing a Block-Gebhart model to calculate the convection heat transfer load of the large-space-nozzle air supply layered air conditioner, verifying the convection heat transfer load by using a reduced scale model experiment result, analyzing and discussing key factors influencing the convection heat transfer load, and making a line calculation chart. The invention can optimize the jet flow model, and the obtained optimized Block-Gebhart extended model can deeply analyze the change condition of the indoor thermal environment of the large-space layered air conditioner, and on the basis, the coupling calculation of the vertical air temperature and the wall surface temperature of the large space can be carried out, so that the calculation method of the convection heat transfer load of the layered air conditioner which is primarily suitable for air supply at the nozzle side is obtained by exploration.

Description

大空间喷口送风分层空调对流热转移负荷计算方法Calculation method of convective heat transfer load for large space nozzle air supply stratified air conditioner

技术领域technical field

本发明涉及一种空调负荷计算技术,特别涉及一种大空间喷口送风分层空调对流热转移负荷计算方法。The invention relates to an air-conditioning load calculation technology, in particular to a method for calculating the convective heat transfer load of a large-space spout air supply stratified air conditioner.

背景技术Background technique

大空间建筑分层空调下的热环境,室内空气垂直温度分层、梯度明显和上下温差大以及屋顶壁面温度分布的特点导致大空间分层空调负荷与普通房间计算方法不同。即大空间室内复杂的气流组织、室外环境变化、室内热源分布等的影响,都会导致室内热环境状态参数的变化和对流热转移量的变化,因此对流热转移负荷计算一直是设计人员的难题。In the thermal environment under stratified air conditioning in large space buildings, the vertical temperature stratification of indoor air, the obvious gradient and the large temperature difference between the upper and lower sides, as well as the characteristics of the temperature distribution on the roof and wall surface, lead to different calculation methods for large space stratified air conditioning load and ordinary rooms. That is, the influence of complex airflow organization, outdoor environment changes, and indoor heat source distribution in large spaces will lead to changes in indoor thermal environment state parameters and changes in convective heat transfer. Therefore, the calculation of convective heat transfer load has always been a problem for designers.

目前大空间建筑分层空调负荷计算方法采用由陆耀庆教授主编的《实用供热空调设计手册》中提及的方法,该方法指出分层空调负荷包括常规空调区负荷(如围护结构传热负荷、室内热源负荷、新风或渗透风负荷等)以及非空调区向空调区转移的辐射转移热和对流转移热形成的负荷,并通过模型实验和现场实测对照,得出影响热转移负荷的因素以及计算方法。具体计算方法为:(1)常规空调区负荷则采用传统全室空调负荷计算方法计算;(2)辐射转移热负荷采用冷负荷系数法:首先通过计算非空调区向空调区地板辐射转移热,再乘以得热系数C1(一般取1.3),得到非空调区向空调区的辐射转移热。继而非空调区向空调区的辐射转移热乘以冷负荷系数C2(取值范围为0.45~0.72,一般取0.5),得到辐射转移负荷;(3)对流热转移负荷计算是通过对某厂房多次试验与现场实测,将结果回归拟合成计算曲线并推广应用至大空间建筑中,将形成对流热转移负荷的关键因素——由喷口送风、内热源、排风三者作用推动的空调区与非空调区间的气流运动及其携带的能量转移简化为空调区热强度、非空调区热强度与排热量三者之间的对比关系,通过查图可以得到对流热转移负荷。At present, the calculation method of layered air-conditioning load in large-space buildings adopts the method mentioned in the "Practical Heating and Air-conditioning Design Manual" edited by Professor Lu Yaoqing. , indoor heat source load, fresh air or infiltration wind load, etc.) and the load formed by the radiative transfer heat and convective transfer heat transferred from the non-air-conditioned area to the air-conditioned area. calculation method. The specific calculation methods are: (1) The conventional air-conditioning area load is calculated by the traditional whole-room air-conditioning load calculation method; (2) The radiation transfer heat load is calculated by the cooling load coefficient method: First, by calculating the radiative transfer heat from the non-air-conditioned area to the air-conditioned area floor, Then multiply by the heat gain coefficient C 1 (usually 1.3) to obtain the radiation transfer heat from the non-air-conditioned area to the air-conditioned area. Then, the radiation transfer heat from the non-air-conditioned area to the air-conditioned area is multiplied by the cooling load coefficient C 2 (the value range is 0.45 to 0.72, generally 0.5) to obtain the radiation transfer load; (3) The convective heat transfer load is calculated by calculating the load of a certain factory building. Repeated tests and field measurements, fitting the results into calculation curves and applying them to large-space buildings will form the key factor of convective heat transfer load - driven by the functions of vent air supply, internal heat source and exhaust air. The airflow movement in the air-conditioned area and the non-air-conditioned area and the energy transfer carried by it are simplified as the comparison between the heat intensity of the air-conditioned area, the heat intensity of the non-air-conditioned area and the heat exhausted. The convective heat transfer load can be obtained by looking at the graph.

从上述可以看出对流热转移负荷作为分层空调冷负荷的重要组成部分,由于上述对流热转移负荷的计算方法研究程度受限于当时大空间建筑热环境,其中采用了许多经验取值与汽轮机高大厂房实测数据,且实验和实测条件单一。基于目前大空间建筑室内热环境研究程度,当时的计算方法存在一定的缺陷与不足。It can be seen from the above that the convective heat transfer load is an important part of the cooling load of the stratified air conditioner. Since the calculation method of the above-mentioned convective heat transfer load is limited by the thermal environment of the large space building at that time, many empirical values and steam turbines are used. The measured data of the tall factory building, and the experimental and measured conditions are single. Based on the current research level of indoor thermal environment of large-space buildings, the calculation method at that time has certain defects and deficiencies.

发明内容SUMMARY OF THE INVENTION

本发明是针对目前大空间建筑分层空调负荷计算方法受限当时的数据的问题,提出了一种大空间喷口送风分层空调对流热转移负荷计算方法,在实验的基础上,提出并建立Block-Gebhart模型,求出稳定传热下的对流热转移负荷理论计算模型。发明中利用缩尺模型实验的结果验证理论模型的正确性与可靠性,进而利用数学模型分析影响对流热转移负荷的主要因素。得到计算对流热转移负荷的线算图并进行实验验证并辅以现场实测,后期再结合辐射转移负荷理论计算模型研究喷口送风下的分层空调辐射转移负荷,探究出系统初步适用于喷口侧送风的分层空调计算方法。Aiming at the problem that the current load calculation method of layered air conditioning in large space buildings is limited by the data at that time, the invention proposes a calculation method of convective heat transfer load of layered air conditioners with large space nozzles. On the basis of experiments, it proposes and establishes The Block-Gebhart model is used to obtain the theoretical calculation model of convective heat transfer load under stable heat transfer. In the invention, the results of the scale model experiment are used to verify the correctness and reliability of the theoretical model, and then the main factors affecting the convective heat transfer load are analyzed by using the mathematical model. The line calculation diagram for calculating the convective heat transfer load was obtained and verified by experiments and supplemented by on-site measurements. Later, combined with the theoretical calculation model of radiation transfer load, the radiation transfer load of stratified air conditioners under the air supply from the nozzle was studied, and it was found that the system is initially suitable for the nozzle side. Calculation method of stratified air conditioning for supply air.

本发明的技术方案为:一种大空间喷口送风分层空调对流热转移负荷计算方法,具体包括如下步骤:The technical scheme of the present invention is: a method for calculating the convective heat transfer load of a large space spout air supply layered air conditioner, which specifically includes the following steps:

1)分析对流热转移负荷,建立Block-Gebhart模型:大空间室内喷口送风系统下,喷口安装高度为分层高度,该高度为非空调区与空调区的分界面,上为非空调区,下为空调区,对流热转移负荷包括非空调区向空调区流动带来的热量以及分界面上的温差换热,Block模型将室内环境在竖直方向上划分为若干个控制区,预测计算大空间建筑室内垂直温度分布,采用Gebhart吸收系数方法,综合考虑室外环境传热、室内内表面辐射换热及室内热源辐射三者共同影响的导热、对流以及辐射耦合换热,求解室内壁面温度分布,并以此作为Block模型计算的边界条件,使用Block-Gebhart模型同步求解对流热转移负荷;1) Analyze the convective heat transfer load and establish the Block-Gebhart model: under the large-space indoor nozzle air supply system, the nozzle installation height is the layer height, which is the interface between the non-air-conditioned area and the air-conditioned area, and the upper is the non-air-conditioned area, The lower part is the air-conditioning area. The convective heat transfer load includes the heat from the non-air-conditioned area to the air-conditioned area and the temperature difference heat transfer on the interface. The Block model divides the indoor environment into several control areas in the vertical direction, and the prediction calculation is large. For the indoor vertical temperature distribution of the space building, the Gebhart absorption coefficient method is used to comprehensively consider the heat conduction, convection and radiation coupled heat transfer jointly affected by the outdoor environmental heat transfer, the indoor inner surface radiative heat transfer and the indoor heat source radiation to solve the indoor wall temperature distribution. Taking this as the boundary condition of the Block model calculation, the Block-Gebhart model is used to solve the convective heat transfer load simultaneously;

2)Block-Gebhart模型求解对流热转移负荷实验验证:在气态缩尺模型实验台喷口送风系统下,进行大空间室内温度场与对流转移热负荷的实验,对模型进行验证,分析影响对流热转移负荷的关键性因素,利用非空调区与空调区热强度比、排热比、无因次对流热转移负荷三者的关系制出线算图,简化计算对流热转移负荷,为实际所用。2) Experiment verification of block-Gebhart model solving convective heat transfer load: Under the gaseous scale model test bench nozzle air supply system, experiments on indoor temperature field and convective transfer heat load in a large space are carried out to verify the model and analyze the influence of convective heat. The key factor of the transfer load is to use the relationship between the heat intensity ratio, heat rejection ratio and dimensionless convective heat transfer load of the non-air-conditioned area and the air-conditioned area to make a line calculation diagram to simplify the calculation of the convective heat transfer load for practical use.

所述步骤1)中Block模型采用基于多区热质平衡Block模型,用于预测大空间建筑室内垂直温度分布;In the described step 1), the Block model adopts the Block model based on multi-zone heat and mass balance, for predicting the indoor vertical temperature distribution of large-space buildings;

假设除空调射流影响区域外,室内环境在水平方向上是均匀分布的,利用集总参数法的概念,针对每个控制主流区域分析质量流动和能量传递情况,分别建立质量和能量平衡方程,然后初设各控制区的主流区域空气温度及其壁面温度,经过一系列迭代运算,最后计算出允许误差内各控制区的温度,从而得到室内垂直温度分布(T1,T2,T3···Tn),n为竖直方向上划分的控制区域总个数;取喷口安装高度为分层高度,该高度为非空调区与空调区的分界面,即1-F层为空调区,F层以上为非空调区,对流热转移负荷计算公式如下:Assuming that the indoor environment is uniformly distributed in the horizontal direction except for the air-conditioning jet affected area, the concept of lumped parameter method is used to analyze the mass flow and energy transfer for each control mainstream area, and the mass and energy balance equations are established respectively, and then The air temperature and wall temperature of the mainstream area of each control area are initially set. After a series of iterative operations, the temperature of each control area within the allowable error is finally calculated, so as to obtain the indoor vertical temperature distribution (T 1 , T 2 , T 3 . . . . ·T n ), n is the total number of control areas divided in the vertical direction; take the installation height of the nozzle as the layer height, which is the interface between the non-air-conditioned area and the air-conditioned area, that is, floors 1-F are the air-conditioned area, Above the F floor is the non-air-conditioned area, and the calculation formula of convective heat transfer load is as follows:

qd=CpMcF(TF+1-TF)+CBFABF(TF+1-TF)q d =C p M cF (T F+1 -T F )+C BF A BF (T F+1 -T F )

式中:qd为非空调区对空调区的对流热转移负荷,单位W;CP为空气定压比热,单位kJ/(kg·℃),对于空气可取值1.01;MCF为非空调区对空调区的空气质量净流量,单位kg/s;TF+1,TF分别为非空调区最低层和空调区域最高层空气温度,单位℃;CBF为非空调区对空调区温差换热系数,单位W/(m2·℃);ABF为非空调区与空调区的分层界面面积,单位m2In the formula: q d is the convective heat transfer load from the non-air-conditioned area to the air-conditioned area, in W; C P is the specific heat of air at constant pressure, in kJ/(kg °C), which can be 1.01 for air; Air quality net flow from air-conditioned area to air-conditioned area, unit kg/s; TF +1 , TF are the air temperature of the lowest floor of non-air-conditioned area and the highest floor of air-conditioned area, unit ℃; C BF is non-air-conditioned area to air-conditioned area Temperature difference heat transfer coefficient, unit W/(m 2 ·°C); A BF is the layered interface area between non-air-conditioning area and air-conditioning area, unit m 2 .

所述步骤1)使用Block-Gebhart模型同步求解对流热转移负荷具体包括步骤如下:Described step 1) uses Block-Gebhart model to solve the convective heat transfer load synchronously and specifically comprises the steps as follows:

A:优化喷口多股射流计算模型,建立壁面流动模型和多区域间热质交换模型;A: Optimize the multi-jet calculation model of the nozzle, and establish the wall flow model and the heat and mass exchange model between multiple regions;

B:利用Block边界条件-壁面温度求解对流辐射耦合换热:根据建筑的几何条件,划分为壁面,并求解壁面之间的角系数和Gebhart吸收系数,通过假设壁面的初始温度分布,根据室内空气温度分布、室外环境参数、围护结构热工参数以及室内热源为边界条件,计算出导热量、辐射量,最后根据耦合热平衡方程,联立求解可得壁面温度的分布;B: Solve convective radiation coupled heat transfer using Block boundary condition - wall temperature: According to the geometric conditions of the building, it is divided into walls, and the angle coefficient between the walls and the Gebhart absorption coefficient are solved. By assuming the initial temperature distribution of the walls, according to the indoor air The temperature distribution, outdoor environmental parameters, thermal parameters of the envelope structure and indoor heat source are used as boundary conditions to calculate the thermal conductivity and radiation.

C:根据Block-Gebhart模型建立风量与能量平衡方程;完成各个区域子模型的计算后,建立每个Block区域的质量平衡和能量平衡方程,然后联立方程组迭代计算并求解出室内空气垂直温度分布、壁面温度分布以及对流热转移负荷,对于任意一块Block主流区域,建立方程;C: According to the Block-Gebhart model, establish the air volume and energy balance equation; after completing the calculation of each area sub-model, establish the mass balance and energy balance equation of each Block area, then iteratively calculate and solve the indoor air vertical temperature of the simultaneous equations Distribution, wall temperature distribution and convective heat transfer load, establish equations for any block mainstream area;

D:Block-Gebhart模型同步求解:D: Block-Gebhart model is solved simultaneously:

(1)假设初始温度,假设初始空气垂直温度分布和壁面温度分布,初始壁面温度分布输入Gebhart模型壁面导热、对流、辐射耦合换热方程中,以Block模型计算的空气垂直温度结果作为边界条件,初始空气垂直温度分布输入Block热质平衡方程中,该模型以Gebhart模型计算的壁面温度结果为计算的边界条件,两个模型参数相互输入;(1) Assuming the initial temperature, assuming the initial vertical air temperature distribution and wall temperature distribution, the initial wall temperature distribution is input into the Gebhart model wall heat conduction, convection and radiation coupled heat transfer equations, and the air vertical temperature results calculated by the Block model are used as boundary conditions, The initial air vertical temperature distribution is input into the Block heat-mass balance equation. The model uses the wall temperature calculated by the Gebhart model as the boundary condition for the calculation, and the two model parameters are input to each other;

(2)迭代计算:将第二步得到的室内空气温度、壁面温度分布与第一步的初设值比较,当两者误差不满足<10-6时,将其赋值于初设值,返回第一步开始重复计算;(2) Iterative calculation: compare the indoor air temperature and wall temperature distribution obtained in the second step with the initial value of the first step, when the error between the two does not satisfy <10 -6 , assign it to the initial value, and return The first step starts to repeat the calculation;

(3)如此循环迭代求解壁面温度分布和室内垂直温度分布,当前后两次两者计算结果的相对误差同时<10-6时,即认为最后两次的室内空气温度分布及壁面温度分布计算结果为问题的解,并取此时平衡方程组计算得到的分层界面上的净流量来计算非空调区对空调区的对流热转移负荷;根据得到的内壁面温度、室内垂直温度以及非空调区与空调区分层界面上的空气净流量,再根据对流热转移负荷计算公式从而得到分层空调的对流热转移负荷。(3) Iteratively solve the wall temperature distribution and the indoor vertical temperature distribution in this way. When the relative errors of the first and last two calculation results are <10 -6 at the same time, it is considered that the last two calculation results of the indoor air temperature distribution and the wall temperature distribution is the solution of the problem, and the net flow at the layered interface calculated by the equilibrium equations is used to calculate the convective heat transfer load from the non-air-conditioned area to the air-conditioned area; according to the obtained inner wall temperature, indoor vertical temperature and non-air-conditioned area The net air flow at the interface with the air conditioner is calculated, and then the convective heat transfer load of the stratified air conditioner is obtained according to the calculation formula of the convective heat transfer load.

本发明的有益效果在于:本发明大空间喷口送风分层空调对流热转移负荷计算方法,可对射流模型进行优化,所得到优化的Block-Gebhart扩展模型可以深入分析大空间分层空调室内热环境变化情况,在此基础上可以进行大空间垂直空气温度和壁面温度的耦合计算,探究得到初步适用于喷口侧送风的分层空调对流热转移负荷的计算方法。The beneficial effects of the present invention are as follows: the method for calculating the convection heat transfer load of the large space spout air supply stratified air conditioner of the present invention can optimize the jet model, and the optimized Block-Gebhart extended model can deeply analyze the indoor heat of the large space stratified air conditioner. Based on the environmental changes, the coupled calculation of the vertical air temperature and the wall temperature in the large space can be carried out, and the calculation method of the convective heat transfer load of the stratified air conditioner initially suitable for the air supply from the nozzle side can be obtained.

附图说明Description of drawings

图1为本发明大空间喷口送风分层空调对流热转移负荷计算方法技术路线图;Fig. 1 is the technical roadmap of the method for calculating the convective heat transfer load of the large-space spout air supply stratified air conditioner of the present invention;

图2为本发明壁面流模型图;Fig. 2 is a wall flow model diagram of the present invention;

图3为本发明室内Block模型示意图;3 is a schematic diagram of an indoor Block model of the present invention;

图4为本发明非空调区与空调区热强度比与对流热转移负荷之间关系图;4 is a graph showing the relationship between the heat intensity ratio of the non-air-conditioned area and the air-conditioned area and the convective heat transfer load of the present invention;

图5为本发明缩尺模型线算图验证图。FIG. 5 is a verification diagram of a line calculation diagram of a scaled model of the present invention.

具体实施方式Detailed ways

如图1所示大空间喷口送风分层空调对流热转移负荷计算方法技术路线图,方法具体包括以下步骤:As shown in Figure 1, the technical roadmap of the calculation method for the convection heat transfer load of the large space nozzle air supply stratified air conditioner, the method specifically includes the following steps:

步骤1:分析对流热转移负荷的形成原因;Step 1: Analyze the cause of the convective heat transfer load;

大空间室内喷口送风系统下,喷口安装高度为分层高度,该高度为非空调区与空调区的分界面,上为非空调区,在喷口送风射流的卷吸作用下,使得非空调区域部分热量转移到空调区域,立即全部成为空调区域的冷负荷。对流热转移发生的根本原因主要归结为非空调区与空调区的以下两点:一是两区域的温差,二是两区域的气流流动,两者缺一不可。根据上述分析,对流热转移负荷包括非空调区向空调区流动带来的热量以及分界面上的温差换热。因此,室内气流的流动和温差决定了对流热转移的发生以及对流热转移热量的大小,而通常影响气流流动和温差的因素主要有如下几点:空调区域送风、回风量,空调区域送风温度,非空调区的排风量,非空调区域与空调区域的得热比。另外有些建筑形式,设有非空调区进风口、以及其他区域设有排风口也均会影响室内气流流动,本发明只针对上述罗列的4个主要因素进行分析。Under the large-space indoor nozzle air supply system, the nozzle installation height is the layered height, which is the interface between the non-air-conditioned area and the air-conditioned area, and the upper part is the non-air-conditioned area. Part of the heat in the area is transferred to the air-conditioned area, which immediately becomes the cooling load of the air-conditioned area. The root cause of convective heat transfer is mainly attributed to the following two points between the non-air-conditioned area and the air-conditioned area: one is the temperature difference between the two areas, and the other is the airflow between the two areas, both of which are indispensable. According to the above analysis, the convective heat transfer load includes the heat flow from the non-air-conditioned area to the air-conditioned area and the temperature difference heat transfer on the interface. Therefore, the flow and temperature difference of indoor air flow determine the occurrence of convective heat transfer and the size of convective heat transfer heat, and the factors that usually affect air flow and temperature difference mainly include the following: air supply and return air volume in the air-conditioning area, and air supply in the air-conditioning area Temperature, exhaust air volume in non-air-conditioned area, heat gain ratio between non-air-conditioned area and air-conditioned area. In addition, in some architectural forms, air inlets in non-air-conditioned areas and air outlets in other areas will also affect the flow of indoor air flow. The present invention only analyzes the four main factors listed above.

步骤2:采用基于多区热质平衡Block模型,用于预测大空间建筑室内垂直温度分布;Step 2: Use the Block model based on multi-zone heat and mass balance to predict the indoor vertical temperature distribution of large-space buildings;

使用Block模型确定室内空气温度分布时,壁面温度分布会影响空气温度分布,故壁面温度的确定方法是非常重要的。本文采用综合考虑室外环境传热、室内内表面辐射换热及室内热源辐射三者共同影响的导热、对流以及辐射耦合换热方法,即Gebhart吸收系数方法,求解室内壁面温度分布,并以此作为Block模型计算的边界条件,使用Block-Gebhart模型同步求解方法来研究大空间建筑室内热环境。由于内表面温度和室内空气温度是相互影响的耦合关系,需分别迭代计算,同步求解计算出以下参数,以此来计算对流热转移负荷:When using the Block model to determine the indoor air temperature distribution, the wall temperature distribution will affect the air temperature distribution, so the method of determining the wall temperature is very important. In this paper, the thermal conduction, convection and radiation coupling heat transfer method, namely the Gebhart absorption coefficient method, which comprehensively considers the joint effects of outdoor environmental heat transfer, indoor inner surface radiative heat transfer and indoor heat source radiation, is used to solve the indoor wall temperature distribution, and use this as the The boundary conditions calculated by the Block model, and the Block-Gebhart model synchronous solution method is used to study the indoor thermal environment of large-space buildings. Since the inner surface temperature and the indoor air temperature are coupled with each other, iterative calculations are required, and the following parameters are calculated simultaneously to calculate the convective heat transfer load:

(1)室内各区域内壁面温度;(1) The temperature of the inner wall in each area of the room;

(2)室内各区域空气温度分布;(2) Air temperature distribution in each area of the room;

(3)各区域之间的质量流量。(3) Mass flow between regions.

Block-Gebhart模型分析大空间喷口送风下的室内热环境时,假设除空调射流影响区域外,室内环境在水平方向上是均匀分布的,利用集总参数法的概念,针对每个控制主流区域分析质量流动和能量传递情况,分别建立质量和能量平衡方程。然后初设各控制区的主流区域空气温度及其壁面温度,经过一系列迭代运算,最后计算出允许误差内各控制区的温度,从而得到室内垂直温度分布(T1,T2,T3···Tn)。When the Block-Gebhart model analyzes the indoor thermal environment under the air supply from the large space nozzle, it is assumed that the indoor environment is uniformly distributed in the horizontal direction except for the area affected by the air-conditioning jet, and the concept of the lumped parameter method is used to control each mainstream area. Analyze mass flow and energy transfer, and establish mass and energy balance equations, respectively. Then, the air temperature and wall temperature of the mainstream area of each control area are initially set, and after a series of iterative operations, the temperature of each control area within the allowable error is finally calculated, so as to obtain the indoor vertical temperature distribution (T 1 , T 2 , T 3 · · · T n ).

根据上述分析,对流热转移负荷包括非空调区向空调区流动带来的热量以及分界面上的温差换热,如图2所示,取喷口安装高度为分层高度,该高度为非空调区与空调区的分界面,即1-4层为空调区,5层及以上为非空调区。因此对流热转移负荷计算公式如下:According to the above analysis, the convective heat transfer load includes the heat brought by the flow from the non-air-conditioned area to the air-conditioned area and the temperature difference heat transfer on the interface. As shown in Figure 2, the installation height of the nozzle is taken as the layer height, which is the non-air-conditioned area. The interface with the air-conditioned area, that is, floors 1-4 are air-conditioned areas, and floors 5 and above are non-air-conditioned areas. Therefore, the calculation formula of convective heat transfer load is as follows:

qd=CpMc5-4(T5-T4)+CB5-4AB5-4(T5-T4) (1)q d =C p M c5-4 (T 5 -T 4 )+C B5-4 A B5-4 (T 5 -T 4 ) (1)

式中:qd——非空调区对空调区的对流热转移负荷,单位W;In the formula: q d — the convective heat transfer load from the non-air-conditioned area to the air-conditioned area, in W;

CP——空气定压比热,单位kJ/(kg·℃),对于空气可取值1.01;C P ——The specific heat of air at constant pressure, in kJ/(kg·℃), it can be taken as 1.01 for air;

MC5-4——非空调区对空调区的空气质量净流量,单位kg/s;M C5-4 —— the air mass net flow from the non-air-conditioned area to the air-conditioned area, in kg/s;

T5,T4——非空调区最低层和空调区域最高层空气温度,单位℃;T 5 , T 4 —— the air temperature of the lowest floor in the non-air-conditioned area and the highest floor in the air-conditioned area, in °C;

CB5-4——非空调区对空调区温差换热系数,单位W/(m2·℃);C B5-4 ——The heat transfer coefficient of the temperature difference between the non-air-conditioned area and the air-conditioned area, in W/(m 2 ·℃);

AB5-4——非空调区与空调区的分层界面面积,单位m2A B5-4 ——The layered interface area between the non-air-conditioned area and the air-conditioned area, in m 2 .

为计算大空间建筑室内垂直温度分布,Block模型将室内环境在竖直方向上划分为若干个控制区。若不考虑自然通风、室内热源的影响,室内气流组织可归纳为3股气流:壁面流、空调射流以及因垂直方向上的温差形成的室内主体气流,假设室内热环境是由此3股射流之间进行热质交换,共同作用的结果。对应这3股气流,在Block模型中纳入3个子模型分别进行热质运动描述:壁面流模型、空调射流模型、区域间热质交换模型。根据区域热质平衡分析方法,Block模型将室内垂直方向上划分为若干区域,如图2所示。通过描述各Block区域之间空气热质传递过程,建立质量平衡和能量平衡方程,求解各Block区域温度,以此得到室内垂直温度分布。In order to calculate the indoor vertical temperature distribution of large-space buildings, the Block model divides the indoor environment into several control areas in the vertical direction. If the influence of natural ventilation and indoor heat source is not considered, the indoor airflow organization can be classified into three airflows: wall flow, air-conditioning jet and indoor main airflow formed by the temperature difference in the vertical direction. The exchange of heat and mass between them is the result of joint action. Corresponding to these three airflows, three sub-models are included in the Block model to describe the motion of heat and mass: the wall flow model, the air-conditioning jet model, and the inter-regional heat and mass exchange model. According to the regional heat and mass balance analysis method, the Block model divides the indoor vertical direction into several regions, as shown in Figure 2. By describing the air heat and mass transfer process between each block area, the mass balance and energy balance equations are established, and the temperature of each block area is solved to obtain the indoor vertical temperature distribution.

步骤3:优化喷口多股射流计算模型;包括单股射流运动轨迹,非等温射流轴心速度、温度衰减公式,多股射流的叠加,卷吸特性等。Step 3: Optimize the calculation model of multiple jets at the nozzle; including the trajectory of a single jet, the non-isothermal jet axial velocity, the temperature decay formula, the superposition of multiple jets, and the entrainment characteristics.

步骤4:建立壁面流动模型和多区域间热质交换模型;围护结构内壁面受室外环境温度、室内环境综合影响,在内表面附近区域会产生沿壁面向上或向下的气流流动。当壁面相对于室内空气温度是热壁面,则该气流流动方向为向上,计算壁面空气流动时须以室内最底层的Block为起点开始计算;当壁面相对于室内空气温度是冷壁面,则该气流流动方向为向下,计算时须从最高层的Block开始逐层向下计算。围护结构根据自身四周的壁面的朝向、材质等热工参数的不同,会造成内壁面温度不同,因此将壁面按对应的Block从最底层到最高层分成K个不同的壁面进行处理。如图3所示壁面流模型图,以夏季为例。壁面为热壁面,壁面温度高于附近空气温度,壁面与室内空气的对流换热引起每层(Block I)空气的出流量MOUT(I,K)汇合成沿壁面的气流,空气沿壁面流动方向向上,壁面流流动以最低Block 1为起点,开始逐层向上计算。中间任意一层Block I产生的质量出流空气(流量MOUT(I,K),温度为T(I))与从下层Block(I-1)流来的混合气流的空气(流量为MMD(I-1,K),温度为TM(I-1,K-1))相遇合成后,形成质量流量为MM(I,K),温度为TM(I,K)。该上升合成流继续沿热壁面向上流动再根据合成流温度TM(I,K)与Block I及Block(I+1)的空气温度三者的大小关系,判断部分或全部流回Block I,其余流量则流入上层Block(I+1),直到全部流入最顶层BlockN。Step 4: Establish a wall flow model and a multi-region heat and mass exchange model; the inner wall of the enclosure is affected by the outdoor ambient temperature and the indoor environment, and the area near the inner surface will generate upward or downward airflow along the wall. When the wall surface is a hot wall surface relative to the indoor air temperature, the airflow direction is upward, and the calculation of the air flow on the wall surface must start from the lowest block in the room; when the wall surface is a cold wall surface relative to the indoor air temperature, the airflow direction The flow direction is downward, and the calculation must start from the block at the highest level and calculate downwards layer by layer. According to the different thermal parameters such as the orientation and material of the surrounding walls of the enclosure structure, the temperature of the inner wall surface will be different. Therefore, the wall surface is divided into K different wall surfaces from the bottom layer to the top layer according to the corresponding Block for processing. As shown in Fig. 3, the wall flow model is taken as an example in summer. The wall is a hot wall, and the temperature of the wall is higher than the temperature of the nearby air. The convective heat exchange between the wall and the indoor air causes the outflow M OUT (I, K) of each layer of (Block I) air to merge into an airflow along the wall, and the air flows along the wall. In the upward direction, the wall flow starts from the lowest Block 1 and starts to calculate upward layer by layer. The mass outflow air (flow rate M OUT (I, K), temperature T(I)) generated by any block I in the middle layer and the mixed airflow from the lower block (I-1) (flow rate M MD) (I-1,K), the temperature is T M (I-1, K-1)) after meeting and synthesis, the mass flow rate is M M (I, K), and the temperature is T M (I, K). The rising synthetic flow continues to flow upward along the hot wall, and then according to the relationship between the synthetic flow temperature T M (I, K) and the air temperature of Block I and Block (I+1), it is judged that part or all of it flows back to Block I, The rest of the traffic flows into the upper block (I+1) until all flows into the topmost BlockN.

壁面附近空气温度可以根据边界层理论分析,Block I沿K壁面处边界层的空气温度按下式计算。The air temperature near the wall can be analyzed according to the boundary layer theory, and the air temperature of the boundary layer along the K wall of Block I is calculated as follows.

TD(I,K)=0.75T(I)+0.25TW(I,K) (2)T D (I,K)=0.75T(I)+0.25T W (I,K) (2)

式中,K为壁面序号,TD(I,K)为Block I沿K壁面处边界层内平均空气温度,T(I)为Block I的空气温度℃,TW(I,K)为壁面K的温度℃。In the formula, K is the wall serial number, T D (I, K) is the average air temperature in the boundary layer along the K wall of Block I, T(I) is the air temperature of Block I in °C, and T W (I, K) is the wall surface temperature in K °C.

建立壁面附近质量流动换热与自然对流换热的平衡式如下式:The equilibrium formula of mass flow heat transfer and natural convection heat transfer near the wall is established as follows:

CPMOUT(I,K)[TD(I,K)-T(I)]=αC(I,K)AW(I,K)[TW(I,K)-T(I)] (3)C P M OUT (I,K)[T D (I,K)-T(I)]=α C (I,K)A W (I,K)[T W (I,K)-T(I )] (3)

式中,MOUT(I,K)为Block I沿K壁面的质量出流量,kg/s,αC(I,K)为壁面K的对流放热系数,W/(m2·℃),AW(I,K)为Block I的壁面K的面积,m2,CP为空气定压比热,W/(m2·℃)。In the formula, M OUT (I, K) is the mass outflow of Block I along the K wall, kg/s, α C (I, K) is the convection heat release coefficient of the wall K, W/(m 2 ·℃), A W (I, K) is the area of the wall surface K of Block I, m 2 , and C P is the specific heat of air at constant pressure, W/(m 2 ·°C).

壁面与空气间的自然对流换热量等于边界层处理流量携带的能量。The natural convection heat transfer between the wall and the air is equal to the energy carried by the boundary layer to process the flow.

Figure BDA0001310376670000081
Figure BDA0001310376670000081

气流合成时,壁面与空气间合成上升流的风量、热量平衡方程如下:When the airflow is synthesized, the air volume and heat balance equation of the synthesized upflow between the wall and the air are as follows:

MOUT(I,K)T(I)+MMD(I-1,K-1)TM(I-1,K-1)=MM(I,K)TM(I,K) (5)M OUT (I,K)T(I)+M MD (I-1,K-1)T M (I-1,K-1)=M M (I,K)T M (I,K) ( 5)

MM(I,K)=MOUT(I,K)+MMD(I-1,K-1) (6)M M (I,K)=M OUT (I,K)+M MD (I-1,K-1) (6)

混合气流的温度TM(I,K):Temperature TM (I,K) of the mixed gas stream:

Figure BDA0001310376670000082
Figure BDA0001310376670000082

且热壁面上升流中从最上层Block(N)流来的混合气流的空气流量MMD(N,K)=0,从底层Block(0)流来的混合气流的空气流量MMD(0,K)=0。(当I=1时,没有前一层I-1向I层的流动量,因此MMD(0,K)=0;当I=N时,没有I层向后一层I+1的流动量,因此MMD(N,K)=0)And the air flow rate M MD (N, K) of the mixed air flow from the uppermost block (N) in the upflow of the hot wall surface is 0, and the air flow rate M MD (0, K)=0. (When I=1, there is no flow from the previous layer I-1 to the I layer, so MMD (0,K)=0; when I=N, there is no flow from the I layer to the next layer I+1 quantity, so M MD (N,K)=0)

每层合成上升流流向下一层具体分配时,混合气流分配过程:When each layer of synthetic upflow flows to the next layer for specific distribution, the mixed airflow distribution process:

MM(I,K)=MIN(I,K)+MMD(I,K) (8)M M (I,K)=M IN (I,K)+M MD (I,K) (8)

混合气流分配时,具体流动方向和流动量由混合气流平均温度TM(I,K)、本层Block I空气层的温度T(I)及下层Block(I-1)空气层的温度T(I+1)三者比较来判断,按比例全部或部分的流入Block I与Block(I+1),即流向本层的流动量MIN(I,K)和流向下一层的流动量MMD(I,K)其比例及判据如所表3示。When the mixed airflow is distributed, the specific flow direction and flow amount are determined by the average temperature of the mixed airflow T M (I, K), the temperature T (I) of the air layer of this layer of Block I and the temperature of the air layer of the lower Block (I-1) T ( I+1) to judge by comparing the three, all or part of the flow into Block I and Block(I+1) in proportion, that is, the flow volume M IN (I, K) flowing to the current layer and the flow volume M flowing to the next layer. The ratio and criterion of MD (I, K) are shown in Table 3.

表3table 3

Figure BDA0001310376670000083
Figure BDA0001310376670000083

Figure BDA0001310376670000091
Figure BDA0001310376670000091

大空间垂直方向上分成若干区域,假设每个Block的主流区域内温度分布是均匀的,由于相邻主流区域的温度存在差异性,会导致区域间的空气流动,同时在相邻Block区域的分界面上也产生热量交换。当上部区域比下部区域温度高时,此时空气温度分布呈稳态,热量会由上部区域传至下部区域;当上部区域比下部区域温度低时,此时空气温度分布不稳定,空气在浮升力作用下流动,直至上下空气温度均匀,此时,上下温差换热可以不予考虑。The large space is vertically divided into several areas. It is assumed that the temperature distribution in the mainstream area of each block is uniform. Due to the temperature difference between adjacent mainstream areas, air flow between areas will be caused. Heat exchange also occurs at the interface. When the temperature of the upper area is higher than that of the lower area, the air temperature distribution is stable, and the heat will be transferred from the upper area to the lower area; when the temperature of the upper area is lower than that of the lower area, the air temperature distribution is unstable at this time, and the air is floating Under the action of lift, it flows until the temperature of the upper and lower air is uniform. At this time, the heat exchange between the upper and lower temperature difference can be ignored.

当Block(I+1)比Block I主流区域空气温度高时,两者换热量如下式(9)所示。When Block(I+1) is higher than the air temperature in the mainstream area of Block I, the heat exchange between the two is shown in the following formula (9).

QB(I+1,I)=CPMC(I+1,I)[T(I+1)-T(I)]+CB(I)AB(I)[T(I+1)-T(I)] (9)Q B (I+1,I)=C P M C (I+1,I)[T(I+1)-T(I)]+C B (I)A B (I)[T(I+ 1)-T(I)] (9)

式中,AB(I)为Block I与Block(I+1)相邻分界面面积m2,CB为温差换热系数W/(m2·℃),T(I)为Block I的空气温度℃,MC(I+1,I)为Block I与Block(I+1)区域间的空气流动kg/s,温差换热系数CB值大小表示Block(I+1)向Block I传递热量的大小程度,随空间分层数量等因素变化。由实验测得,当T(I+1)>T(I)时,CB=2.3W/(m2·℃);当T(I+1)<T(I)时,由密度差引起的热传递流强度增大,CB=116W/(m2·℃)。In the formula, A B (I) is the adjacent interface area m 2 between Block I and Block (I+1), C B is the temperature difference heat transfer coefficient W/(m 2 ·℃), and T(I) is the Air temperature ℃, M C (I+1, I) is the air flow kg/s between Block I and Block (I+1), and the value of temperature difference heat transfer coefficient C B means Block (I+1) to Block I The degree of heat transfer varies with factors such as the number of spatial layers. Measured by experiment, when T(I+1)>T(I), C B =2.3W/(m 2 ·℃); when T(I+1)<T(I), it is caused by density difference The intensity of the heat transfer flow increases, C B =116W/(m 2 ·°C).

步骤5:利用Block边界条件-壁面温度求解对流辐射耦合换热;根据建筑的几何条件,划分为壁面,并求解壁面之间的角系数和Gebhart吸收系数。通过假设壁面的初始温度分布,根据室内空气温度分布、室外环境参数、围护结构热工参数以及室内热源等为边界条件,计算出导热量、辐射量等。最后根据耦合热平衡方程,联立求解可得壁面温度的分布。Step 5: Use the block boundary condition - wall temperature to solve the convection radiation coupled heat transfer; according to the geometric conditions of the building, divide it into walls, and solve the angle coefficient between the walls and the Gebhart absorption coefficient. By assuming the initial temperature distribution of the wall, according to the indoor air temperature distribution, outdoor environmental parameters, thermal parameters of the enclosure structure, and indoor heat sources as boundary conditions, the thermal conductivity, radiation, etc. are calculated. Finally, according to the coupled heat balance equation, the temperature distribution of the wall can be obtained by solving simultaneously.

步骤6:根据Block-Gebhart模型建立风量与能量平衡方程;完成上述各个子模型的计算后,可以建立每个Block区域的质量平衡和能量平衡方程,然后联立方程组迭代计算并求解出室内空气垂直温度分布、壁面温度分布以及对流热转移负荷。对于任意一块Block主流区域,建立方程。Step 6: Establish air volume and energy balance equations according to the Block-Gebhart model; after completing the calculations of the above sub-models, the mass balance and energy balance equations of each block area can be established, and then the simultaneous equations are iteratively calculated and solved to obtain the indoor air Vertical temperature distribution, wall temperature distribution, and convective heat transfer loads. For any block mainstream area, establish an equation.

步骤7:Block-Gebhart模型同步求解;此模型同步求解的计算步骤如下:Step 7: Simultaneous solution of the Block-Gebhart model; the calculation steps for the simultaneous solution of this model are as follows:

其计算步骤包括:(1)假设初始温度。假设初始空气垂直温度分布和壁面温度分布,初始壁面温度分布输入Gebhart模型壁面导热、对流、辐射耦合换热方程中,以Block模型计算的空气垂直温度结果作为边界条件,初始空气垂直温度分布输入Block热质平衡方程中,该模型以Gebhart模型计算的壁面温度结果为计算的边界条件,两个模型参数相互输入;(2)迭代计算。将第二步得到的室内空气温度、壁面温度分布与第一步的初设值比较,当两者误差不满足<10-6时,将其赋值于初设值,返回第一步开始重复计算。(3)如此循环迭代求解壁面温度分布和室内垂直温度分布,当前后两次两者计算结果的相对误差同时<10-6时,即认为最后两次的室内空气温度分布及壁面温度分布计算结果为问题的解,并取此时平衡方程组计算得到的分层界面上的净流量来计算非空调区对空调区的对流热转移负荷。根据得到的内壁面温度、室内垂直温度以及非空调区与空调区分层界面上的空气净流量,再根据qd=CpMc5-4(T5-T4)+CB5-4AB5-4(T5-T4)计算方法从而得到分层空调的对流热转移负荷,其中qd非空调区对空调区的对流热转移负荷,CP空气定压比热,MC5-4非空调区对空调区的空气质量净流量,T5,T4非空调区和空调区域空气温度,CB5-4温差换热系数,AB5-4非空调区与空调区的分层界面面积。The calculation steps include: (1) Assume the initial temperature. Assuming the initial vertical air temperature distribution and wall temperature distribution, the initial wall temperature distribution is input into the Gebhart model wall heat conduction, convection, and radiation coupled heat transfer equations, the air vertical temperature results calculated by the Block model are used as boundary conditions, and the initial air vertical temperature distribution is input into Block In the heat-mass balance equation, the model takes the wall temperature calculated by the Gebhart model as the boundary condition for the calculation, and the two model parameters are mutually input; (2) Iterative calculation. Compare the indoor air temperature and wall temperature distribution obtained in the second step with the initial value of the first step. When the error between the two does not satisfy <10 -6 , assign it to the initial value, and return to the first step to repeat the calculation. . (3) Iteratively solve the wall temperature distribution and the indoor vertical temperature distribution in this way. When the relative errors of the first and last two calculation results are <10 -6 at the same time, it is considered that the last two calculation results of the indoor air temperature distribution and the wall temperature distribution is the solution of the problem, and the net flow on the layered interface calculated by the equilibrium equations at this time is used to calculate the convective heat transfer load from the non-air-conditioned area to the air-conditioned area. According to the obtained inner wall surface temperature, indoor vertical temperature and air net flow at the interface between the non-air-conditioned area and the air-conditioned area, according to q d =C p M c5-4 (T 5 -T 4 )+C B5-4 A B5-4 (T 5 -T 4 ) calculation method to obtain the convective heat transfer load of stratified air conditioning, where q d is the convective heat transfer load of the non-air-conditioned area to the air-conditioned area, C P air constant pressure specific heat, M C5-4 Air mass net flow from non-air-conditioned area to air-conditioned area, T 5 , T 4 air temperature in non-air-conditioned area and air-conditioned area, C B5-4 temperature difference heat transfer coefficient, A B5-4 stratified interface area between non-air-conditioned area and air-conditioned area .

步骤8:Block-Gebhart模型求解对流热转移负荷实验验证;在气态缩尺模型实验台喷口送风系统下,进行大空间室内温度场与对流转移热负荷的实验研究,根据实验室的高度和喷口所在位置,在垂直方向上将缩尺模型分为4个区域,对这4个Block区域建立质量与能量平衡方程。将Block-Gebhart模型计算结果与实验数据结果进行对比,可以发现对流热转移负荷的实验值与理论值之间的最大误差为11.50%,六个工况的平均误差为6.60%,由此可以说明利用Block-Gebhart模型计算对流热转移负荷的可行性。Step 8: The block-Gebhart model solves the convective heat transfer load experimental verification; under the gas scale model test bench nozzle air supply system, the experimental research on the indoor temperature field and convective transfer heat load in a large space is carried out. According to the height of the laboratory and the nozzle At the location, the scaled model is divided into 4 regions in the vertical direction, and the mass and energy balance equations are established for these 4 block regions. Comparing the calculation results of the Block-Gebhart model with the experimental data results, it can be found that the maximum error between the experimental value and the theoretical value of the convective heat transfer load is 11.50%, and the average error of the six operating conditions is 6.60%, which can be explained Feasibility of calculating convective heat transfer loads using the Block-Gebhart model.

步骤9:将对流热转移负荷与其关键性的影响因素之间的相互关系转化为无因次对流热转移负荷(相对非空调区得热量下的对流热转移负荷)、非空调区与空调区热强度比、排热比(相对非空调区得热量下的排热量)三者之间的关系并制成线算图以便于实际应用。Step 9: Convert the relationship between the convective heat transfer load and its key influencing factors into dimensionless convective heat transfer load (convective heat transfer load relative to the heat gain in the non-air-conditioned area), heat in the non-air-conditioned area and the air-conditioned area. The relationship between the intensity ratio and the heat rejection ratio (relative to the heat rejection in the non-air-conditioned area) is made into a line graph for practical application.

构建先算图的具体思路为:(1)确定研究的热强度比分别是0.32,0.4,0.45,0.52,0.58;(2)分别计算出非空调区域和空调区各个外围护结构占总得热量的比例,将非空调区域热强度按照这个比例再分配到非空调区域的内热源上分别为:0,100W,200W,300W,400W;(3)按照不同热强度比下模型计算的对流转移热负荷,分别计算无因次对流热转移负荷;(4)改变排热比,分别为0%,10%,20%,按照上述方式计算,得到最终的线算图曲线,如图4所示。由图可知,在同一的排热比下,随着热强度比的增大,从非空调区域流向空调区域的对流热转移量也逐渐增大。这是由于得热量直接作用于非空调区域,使得非空调区域的温度上升,就会产生非空调区域向空调区域的热扩散。这样,就加大了非空调区域向空调区域的热转移量。随着排热比增大,非空调区热量被渐渐排除,减小了非空调区与空调区的温差,减弱了非空调区域向空调区域的热扩散,从而减小了非空调区域向空调区域的热转移量。The specific idea of constructing the first calculation map is as follows: (1) Determine the heat intensity ratios of the study to be 0.32, 0.4, 0.45, 0.52, 0.58; (2) Calculate the total heat gain of each outer protection structure in the non-air-conditioned area and the air-conditioned area. The heat intensity of the non-air-conditioned area is redistributed to the internal heat sources of the non-air-conditioned area according to this ratio: 0, 100W, 200W, 300W, 400W; (3) According to the convective transfer heat calculated by the model under different heat intensity ratios (4) Change the heat rejection ratio to 0%, 10%, and 20%, respectively, and calculate according to the above method to obtain the final line calculation curve, as shown in Figure 4. It can be seen from the figure that under the same heat rejection ratio, as the heat intensity ratio increases, the amount of convective heat transfer from the non-air-conditioned area to the air-conditioned area also increases gradually. This is because the heat gain directly acts on the non-air-conditioned area, causing the temperature of the non-air-conditioned area to rise, resulting in heat diffusion from the non-air-conditioned area to the air-conditioned area. In this way, the amount of heat transfer from the non-air-conditioned area to the air-conditioned area is increased. As the heat rejection ratio increases, the heat from the non-air-conditioned area is gradually removed, reducing the temperature difference between the non-air-conditioned area and the air-conditioned area, reducing the heat diffusion from the non-air-conditioned area to the air-conditioned area, thereby reducing the non-air-conditioned area to the air-conditioned area. amount of heat transfer.

步骤10:对线算图的缩尺模型验证,根据图5可以看出,0%排热比和10%排热比的实验数据与原手册计算曲线值差异较大,与本次研究的对流热转移负荷线算图基本吻合,可以基本说明线算图应用的准确性以及适用性。Step 10: Verify the scaled model of the line calculation diagram. According to Figure 5, it can be seen that the experimental data of 0% heat rejection ratio and 10% heat rejection ratio are quite different from the calculated curve values in the original manual, which are different from the convection in this study. The heat transfer load line calculation diagram is basically consistent, which can basically explain the accuracy and applicability of the line calculation diagram.

最后,本发明所制线算图从气流流动入手,使用理论和实验相结合并辅以现场实测的方式得到简易的计算对流热转移负荷的方法。本发明基于Block-Gebhart理论模型来分析大空间建筑的温度分布、气流流动等室内热环境,根据非空调区与空调区分界面区域间的温差和流动量来计算对流热转移负荷,继而分析出影响对流热转移负荷的关键因素喷口送风参数、室外传到室内的热量分布、排热等,利用非空调区与空调区热强度比、排热比、无因次对流热转移负荷三者的关系制出线算图,以便于简便计算对流热转移负荷。Finally, the line calculation diagram of the present invention starts from the airflow flow, and uses a combination of theory and experiment, supplemented by field measurement, to obtain a simple method for calculating the convective heat transfer load. Based on the Block-Gebhart theoretical model, the present invention analyzes the indoor thermal environment such as temperature distribution and airflow flow of large-space buildings, calculates the convective heat transfer load according to the temperature difference and flow between the non-air-conditioned area and the air-conditioned interface area, and then analyzes the influence The key factors of the convective heat transfer load are the air supply parameters of the nozzle, the heat distribution from the outside to the room, and the heat removal, etc., using the relationship between the heat intensity ratio, heat rejection ratio, and dimensionless convective heat transfer load between the non-air-conditioned area and the air-conditioned area Draw a line diagram for easy calculation of convective heat transfer loads.

Claims (2)

1.一种大空间喷口送风分层空调对流热转移负荷计算方法,其特征在于,具体包括如下步骤:1. a large space spout air supply layered air conditioning convection heat transfer load calculation method, is characterized in that, specifically comprises the steps: 1)分析对流热转移负荷,建立Block-Gebhart模型:大空间室内喷口送风系统下,喷口安装高度为分层高度,该高度为非空调区与空调区的分界面,上为非空调区,下为空调区,对流热转移负荷包括非空调区向空调区流动带来的热量以及分界面上的温差换热,Block模型将室内环境在竖直方向上划分为若干个控制区,预测计算大空间建筑室内垂直温度分布,采用Gebhart吸收系数方法,综合考虑室外环境传热、室内内表面辐射换热及室内热源辐射三者共同影响的导热、对流以及辐射耦合换热,求解室内壁面温度分布,并以此作为Block模型计算的边界条件,使用Block-Gebhart模型同步求解对流热转移负荷;使用Block-Gebhart模型同步求解对流热转移负荷具体包括步骤如下:1) Analyze the convective heat transfer load and establish the Block-Gebhart model: under the large-space indoor nozzle air supply system, the nozzle installation height is the layer height, which is the interface between the non-air-conditioned area and the air-conditioned area, and the upper is the non-air-conditioned area, The lower part is the air-conditioning area. The convective heat transfer load includes the heat from the non-air-conditioned area to the air-conditioned area and the temperature difference heat transfer on the interface. The Block model divides the indoor environment into several control areas in the vertical direction, and the prediction calculation is large. For the indoor vertical temperature distribution of the space building, the Gebhart absorption coefficient method is used to comprehensively consider the heat conduction, convection and radiation coupled heat transfer jointly affected by the outdoor environmental heat transfer, the indoor inner surface radiative heat transfer and the indoor heat source radiation to solve the indoor wall temperature distribution. Taking this as the boundary condition for the calculation of the Block model, the Block-Gebhart model is used to solve the convective heat transfer load simultaneously; the steps to solve the convective heat transfer load simultaneously using the Block-Gebhart model are as follows: A:优化喷口多股射流计算模型,建立壁面流动模型和多区域间热质交换模型;A: Optimize the multi-jet calculation model of the nozzle, and establish the wall flow model and the heat and mass exchange model between multiple regions; B:利用Block边界条件-壁面温度求解对流辐射耦合换热:根据建筑的几何条件,划分为壁面,并求解壁面之间的角系数和Gebhart吸收系数,通过假设壁面的初始温度分布,根据室内空气温度分布、室外环境参数、围护结构热工参数以及室内热源为边界条件,计算出导热量、辐射量,最后根据耦合热平衡方程,联立求解可得壁面温度的分布;B: Solve convective radiation coupled heat transfer using Block boundary condition - wall temperature: According to the geometric conditions of the building, it is divided into walls, and the angle coefficient between the walls and the Gebhart absorption coefficient are solved. By assuming the initial temperature distribution of the walls, according to the indoor air The temperature distribution, outdoor environmental parameters, thermal parameters of the envelope structure and indoor heat source are used as boundary conditions to calculate the thermal conductivity and radiation. C:根据Block-Gebhart模型建立风量与能量平衡方程;完成各个区域子模型的计算后,建立每个Block区域的质量平衡和能量平衡方程,然后联立方程组迭代计算并求解出室内空气垂直温度分布、壁面温度分布以及对流热转移负荷,对于任意一块Block主流区域,建立方程;C: According to the Block-Gebhart model, establish the air volume and energy balance equation; after completing the calculation of each area sub-model, establish the mass balance and energy balance equation of each Block area, then iteratively calculate and solve the indoor air vertical temperature of the simultaneous equations Distribution, wall temperature distribution and convective heat transfer load, establish equations for any block mainstream area; D:Block-Gebhart模型同步求解:D: Block-Gebhart model is solved simultaneously: (1)假设初始温度,假设初始空气垂直温度分布和壁面温度分布,初始壁面温度分布输入Gebhart模型壁面导热、对流、辐射耦合换热方程中,以Block模型计算的空气垂直温度结果作为边界条件,初始空气垂直温度分布输入Block热质平衡方程中,该模型以Gebhart模型计算的壁面温度结果为计算的边界条件,两个模型参数相互输入;(1) Assuming the initial temperature, assuming the initial vertical air temperature distribution and wall temperature distribution, the initial wall temperature distribution is input into the Gebhart model wall heat conduction, convection and radiation coupled heat transfer equations, and the air vertical temperature results calculated by the Block model are used as boundary conditions, The initial air vertical temperature distribution is input into the Block heat-mass balance equation. The model uses the wall temperature calculated by the Gebhart model as the boundary condition for the calculation, and the two model parameters are input to each other; (2)迭代计算:将第二步得到的室内空气温度、壁面温度分布与第一步的初设值比较,当两者误差不满足<10-6时,将其赋值于初设值,返回第一步开始重复计算;(2) Iterative calculation: compare the indoor air temperature and wall temperature distribution obtained in the second step with the initial value of the first step, when the error between the two does not satisfy <10 -6 , assign it to the initial value, and return The first step starts to repeat the calculation; (3)如此循环迭代求解壁面温度分布和室内垂直温度分布,当前后两次两者计算结果的相对误差同时<10-6时,即认为最后两次的室内空气温度分布及壁面温度分布计算结果为问题的解,并取此时平衡方程组计算得到的分层界面上的净流量来计算非空调区对空调区的对流热转移负荷;根据得到的内壁面温度、室内垂直温度以及非空调区与空调区分层界面上的空气净流量,再根据对流热转移负荷计算公式从而得到分层空调的对流热转移负荷;(3) Iteratively solve the wall temperature distribution and the indoor vertical temperature distribution in this way. When the relative errors of the first and last two calculation results are <10 -6 at the same time, it is considered that the last two calculation results of the indoor air temperature distribution and the wall temperature distribution is the solution of the problem, and the net flow at the layered interface calculated by the equilibrium equations is used to calculate the convective heat transfer load from the non-air-conditioned area to the air-conditioned area; according to the obtained inner wall temperature, indoor vertical temperature and non-air-conditioned area The net air flow at the interface with the air conditioner, and then according to the convective heat transfer load calculation formula, the convective heat transfer load of the stratified air conditioner is obtained; 2)Block-Gebhart模型求解对流热转移负荷实验验证:在气态缩尺模型实验台喷口送风系统下,进行大空间室内温度场与对流转移热负荷的实验,对模型进行验证,分析影响对流热转移负荷的关键性因素,利用非空调区与空调区热强度比、排热比、无因次对流热转移负荷三者的关系制出线算图,简化计算对流热转移负荷,为实际所用。2) Experiment verification of block-Gebhart model solving convective heat transfer load: Under the gaseous scale model test bench nozzle air supply system, experiments on indoor temperature field and convective transfer heat load in a large space are carried out to verify the model and analyze the influence of convective heat. The key factor of the transfer load is to use the relationship between the heat intensity ratio, heat rejection ratio and dimensionless convective heat transfer load of the non-air-conditioned area and the air-conditioned area to make a line calculation diagram to simplify the calculation of the convective heat transfer load for practical use. 2.根据权利要求1所述大空间喷口送风分层空调对流热转移负荷计算方法,其特征在于,所述步骤1)中Block模型采用基于多区热质平衡Block模型,用于预测大空间建筑室内垂直温度分布;2. according to the described large space spout air supply layered air conditioner convection heat transfer load calculation method of claim 1, it is characterized in that, in described step 1), Block model adopts Block model based on multi-zone thermal mass balance, for predicting large space Building indoor vertical temperature distribution; 假设除空调射流影响区域外,室内环境在水平方向上是均匀分布的,利用集总参数法的概念,针对每个控制主流区域分析质量流动和能量传递情况,分别建立质量和能量平衡方程,然后初设各控制区的主流区域空气温度及其壁面温度,经过一系列迭代运算,最后计算出允许误差内各控制区的温度,从而得到室内垂直温度分布(T1,T2,T3···Tn),n为竖直方向上划分的控制区域总个数;取喷口安装高度为分层高度,该高度为非空调区与空调区的分界面,即1-F层为空调区,F层以上为非空调区,对流热转移负荷计算公式如下:Assuming that the indoor environment is uniformly distributed in the horizontal direction except for the air-conditioning jet affected area, the concept of lumped parameter method is used to analyze the mass flow and energy transfer for each control mainstream area, and the mass and energy balance equations are established respectively, and then The air temperature and wall temperature of the mainstream area of each control area are initially set. After a series of iterative operations, the temperature of each control area within the allowable error is finally calculated, so as to obtain the indoor vertical temperature distribution (T 1 , T 2 , T 3 . . . . ·T n ), n is the total number of control areas divided in the vertical direction; take the installation height of the nozzle as the layer height, which is the interface between the non-air-conditioned area and the air-conditioned area, that is, floors 1-F are the air-conditioned area, Above the F floor is the non-air-conditioned area, and the calculation formula of convective heat transfer load is as follows: qd=CpMcF(TF+1-TF)+CBFABF(TF+1-TF)q d =C p M cF (T F+1 -T F )+C BF A BF (T F+1 -T F ) 式中:qd为非空调区对空调区的对流热转移负荷,单位W;CP为空气定压比热,单位kJ/(kg·℃),对于空气可取值1.01;MCF为非空调区对空调区的空气质量净流量,单位kg/s;TF+1,TF分别为非空调区最低层和空调区域最高层空气温度,单位℃;CBF为非空调区对空调区温差换热系数,单位W/(m2·℃);ABF为非空调区与空调区的分层界面面积,单位m2In the formula: q d is the convective heat transfer load from the non-air-conditioned area to the air-conditioned area, in W; C P is the specific heat of air at constant pressure, in kJ/(kg °C), which can be 1.01 for air; Air quality net flow from air-conditioned area to air-conditioned area, unit kg/s; TF +1 , TF are the air temperature of the lowest floor of non-air-conditioned area and the highest floor of air-conditioned area, unit ℃; C BF is non-air-conditioned area to air-conditioned area Temperature difference heat transfer coefficient, unit W/(m 2 ·°C); A BF is the layered interface area between non-air-conditioning area and air-conditioning area, unit m 2 .
CN201710403639.8A 2017-06-01 2017-06-01 Method for calculating convection heat transfer load of large-space nozzle air supply layered air conditioner Active CN107301276B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710403639.8A CN107301276B (en) 2017-06-01 2017-06-01 Method for calculating convection heat transfer load of large-space nozzle air supply layered air conditioner

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710403639.8A CN107301276B (en) 2017-06-01 2017-06-01 Method for calculating convection heat transfer load of large-space nozzle air supply layered air conditioner

Publications (2)

Publication Number Publication Date
CN107301276A CN107301276A (en) 2017-10-27
CN107301276B true CN107301276B (en) 2020-08-04

Family

ID=60137283

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710403639.8A Active CN107301276B (en) 2017-06-01 2017-06-01 Method for calculating convection heat transfer load of large-space nozzle air supply layered air conditioner

Country Status (1)

Country Link
CN (1) CN107301276B (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109522662B (en) * 2018-11-23 2022-12-02 北京展坤博宇节能科技有限公司 Method for calculating mixed air temperature of roof type air conditioner
CN110704926B (en) * 2019-09-27 2021-07-20 上海理工大学 Calculation method of unsteady radiative heat transfer load for layered air conditioning in large space buildings
CN111473463B (en) * 2020-03-27 2021-06-25 中国铁建电气化局集团有限公司 A measuring point arrangement method for temperature and humidity monitoring in tall spaces
CN111706951B (en) * 2020-07-03 2021-08-24 湖南大学 A Thermal Comfort Ventilation and Pollutant Control Method for Air Stability in Limited Space
CN114413536A (en) * 2021-12-15 2022-04-29 合肥通用机械研究院有限公司 Heat load calculation method, detector, refrigeration control method and system for ice rink
CN114234384B (en) * 2021-12-27 2023-08-29 中铁第四勘察设计院集团有限公司 Air conditioning optimization control method and system for railway passenger station

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101049069A (en) * 2006-04-07 2007-10-10 上海科技馆 Technical method for planting rain forest plant in large indoor showing region
CN102819643A (en) * 2012-08-14 2012-12-12 中冶南方工程技术有限公司 Modeling method and application of heat transfer and heat storage simulation model of hot blast stove
CN106019070A (en) * 2016-07-05 2016-10-12 贵州电网有限责任公司电力科学研究院 On-line monitoring device and monitoring method for working state of copper-aluminum composite connection terminal
CN106202757A (en) * 2016-07-15 2016-12-07 东北大学 The numerical value of a kind of Con casting ladle space radiation heat flux distribution determines method

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6644244B2 (en) * 2001-12-05 2003-11-11 University Of Florida Apparatus and methods for testing pain sensitivity
US10208978B2 (en) * 2012-11-08 2019-02-19 Lennox Industries Inc. System for generating electrical energy from waste energy

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101049069A (en) * 2006-04-07 2007-10-10 上海科技馆 Technical method for planting rain forest plant in large indoor showing region
CN102819643A (en) * 2012-08-14 2012-12-12 中冶南方工程技术有限公司 Modeling method and application of heat transfer and heat storage simulation model of hot blast stove
CN106019070A (en) * 2016-07-05 2016-10-12 贵州电网有限责任公司电力科学研究院 On-line monitoring device and monitoring method for working state of copper-aluminum composite connection terminal
CN106202757A (en) * 2016-07-15 2016-12-07 东北大学 The numerical value of a kind of Con casting ladle space radiation heat flux distribution determines method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Application of the Gebhart-Block Model for Predicting Vertical Temperature Distribution in a Large Space Building with Natural Ventilation;chen huang 等;《Maximize Comfort: Temperature, Humidity and IAQ》;20061109;第1-4页 *
应用Gebhart-Block模型预测大空间自然通风下室内温度分布;宋岩 等;《暖通空调》;20081231;第38卷(第12期);第22-25页 *

Also Published As

Publication number Publication date
CN107301276A (en) 2017-10-27

Similar Documents

Publication Publication Date Title
CN107301276B (en) Method for calculating convection heat transfer load of large-space nozzle air supply layered air conditioner
Ye et al. Heating energy consumption of impinging jet ventilation and mixing ventilation in large-height spaces: A comparison study
Megri et al. Zonal modeling for simulating indoor environment of buildings: Review, recent developments, and applications
Su et al. Double-skin façade optimization design for different climate zones in China
Liu et al. Influence of air-conditioning systems on buoyancy driven air infiltration in large space buildings: A case study of a railway station
CN103995548B (en) Indoor thermal environment control method based on model reduction and multiple model predictive control
Zhang et al. A simplified approach to describe complex diffusers in displacement ventilation for CFD simulations
CN108830932B (en) An energy consumption prediction method for large-space buildings based on the coupling of EnergyPlus and CFD
Wang et al. Design configuration for a higher efficiency air conditioning system in large space building
Zhang et al. Numerical analysis of diffuse ceiling ventilation and its integration with a radiant ceiling system
Huang et al. Study of indoor thermal environment and stratified air-conditioning load with low-sidewall air supply for large space based on Block-Gebhart model
Yang et al. Heat transfer between occupied and unoccupied zone in large space building with floor-level side wall air-supply system
Liu et al. Theoretical model of buoyancy-driven air infiltration during heating/cooling seasons in large space buildings
Lu et al. Evaluation of stack ventilation in a large space using zonal simulation and a reduced-scale model experiment with particle image velocimetry
Xu et al. Study of convective heat transfer load induced by nozzle air supply in large spaces with thermal stratification based on Block-Gebhart model
CN111666723A (en) Large-space all-air conditioning system optimization design method based on numerical simulation
Qi et al. Froude-Stanton modeling of heat and mass transfer in large vertical spaces of high-rise buildings
KR102368921B1 (en) Method for predicting a vertical air temperature distribution in large space structure and mehtod for calculating heating load and cooling load in the large space sutructure using the said predicted vertical air temperature distribution
Fang et al. Study of underfloor air distribution using zonal model-based simulation and experimental measurements
Ayoub et al. Simplified thermal model of spaces cooled with combined positive displacement ventilation and chilled ceiling system
Vidhyashankar et al. Modelling spatial variations in thermal comfort in indoor open-plan spaces using a whole-building simulation tool
Novoselac et al. New convection correlations for cooled ceiling panels in room with mixed and stratified airflow
Wang et al. A simplified cooling load calculation method based on equivalent heat transfer coefficient for large space buildings with a stratified air-conditioning system
Ji et al. Numerical studies of displacement natural ventilation in multi-storey buildings connected to an atrium
Lu et al. Zoning strategy of zonal modeling for thermally stratified large spaces

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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20220106

Address after: 200241 room 1502, building 2, Lane 39, Shangyi Road, Minhang District, Shanghai

Patentee after: SHANGHAI SHINELONG AIR CONDITIONING Co.,Ltd.

Address before: 200093 No. 516, military road, Shanghai, Yangpu District

Patentee before: University of Shanghai for Science and Technology