CN113987779A - 一种基于热带海岛气候特征的城市冠层模型建立方法 - Google Patents
一种基于热带海岛气候特征的城市冠层模型建立方法 Download PDFInfo
- Publication number
- CN113987779A CN113987779A CN202111236079.4A CN202111236079A CN113987779A CN 113987779 A CN113987779 A CN 113987779A CN 202111236079 A CN202111236079 A CN 202111236079A CN 113987779 A CN113987779 A CN 113987779A
- Authority
- CN
- China
- Prior art keywords
- street
- valley
- flux
- heat
- air
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000004907 flux Effects 0.000 claims abstract description 106
- 230000005855 radiation Effects 0.000 claims abstract description 57
- 230000000694 effects Effects 0.000 claims abstract description 15
- 238000005338 heat storage Methods 0.000 claims abstract description 8
- 238000012545 processing Methods 0.000 claims abstract description 8
- 239000010410 layer Substances 0.000 claims description 58
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 9
- 239000000126 substance Substances 0.000 claims description 9
- 230000008859 change Effects 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 6
- 238000004519 manufacturing process Methods 0.000 claims description 5
- 238000009423 ventilation Methods 0.000 claims description 5
- 230000010354 integration Effects 0.000 claims description 4
- 239000002344 surface layer Substances 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000004378 air conditioning Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000004060 metabolic process Effects 0.000 claims description 3
- 238000002310 reflectometry Methods 0.000 claims description 3
- 238000009834 vaporization Methods 0.000 claims description 3
- 230000008016 vaporization Effects 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 230000020169 heat generation Effects 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 6
- 230000008901 benefit Effects 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 238000001816 cooling Methods 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 208000035473 Communicable disease Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 239000005431 greenhouse gas Substances 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
- G01W1/10—Devices for predicting weather conditions
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
- Y02A30/60—Planning or developing urban green infrastructure
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Data Mining & Analysis (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Mathematical Physics (AREA)
- Computer Hardware Design (AREA)
- Life Sciences & Earth Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Road Paving Structures (AREA)
Abstract
Description
技术领域
本发明涉及城市冠层领域,尤其涉及一种基于热带海岛气候特征的城市冠层模型建立方法。
背景技术
在目前对于城市热岛效应的研究过程中,在给城市街谷建模时一般采用经典的城市冠层模型。但是由于城市的经纬度、本地气候特征、地貌和用地特征等差别巨大,城市冠层模型中的参数化处理方式差异很大,城市冠层模型中水平通风特性对热岛效应的影响基本采用“很小而忽略”的处理方式,对研究“热带海岛”城市的热岛效应是不足的,因为热带海岛具有常年大风带来的强水平通风特性。
发明内容
鉴于上述问题,提出了本发明以便提供克服上述问题或者至少部分地解决上述问题的一种基于热带海岛气候特征的城市冠层模型建立方法。
根据本发明的一个方面,提供了一种基于热带海岛气候特征的城市冠层模型建立方法,所述建立方法包括:
将相邻区域联系在一起,区域内的有限长度的多个街道及相互影响;
街道内的能量平衡方程如下:
Q* s+QF,s=QH,s+QE,s+ΔQS,s+ΔQA,s
其中,Q*为净辐射热通量,W/m2;QF为人为产热,W/m2;QH为显热通量,W/m2;QE为潜热通量,W/m2;ΔQS为净储热通量,W/m2;ΔQA为净对流热通量,W/m2;s为街道;
所述净辐射热通量Q*=长波辐射+短波辐射。
可选的,所述长波辐射包括:
路面的长波净辐射
墙面的长波净辐射
其中,L↓为太阳辐射量,∈r、∈w为路面和墙壁的发射率,对路面来说,若Ψr为路面对天空的视角系数,则(1-Ψr)为路面对两边墙壁的视角系数;墙壁对天空视角系数为Ψw,对路面的视角系数为Ψw,则对对面墙的视角系数为(1-2Ψw),对屋顶来说,天空视角为1.0;
采用平面角度计算天空视角系数,路面w/2位置的天空视角系数为
可选的,所述短波辐射包括:
根据街道与太阳方向的垂直角度计算路面、西墙、东墙和屋顶的平均太阳直接辐射通量是:
根据所述街谷的朝向变化,将所述街谷的宽度w修正为w/sinθ;在得到墙体热通量后,需要再乘上sinθ来修正,θ0为路面完全没有得到阳光直射的街道峡谷的朝向
墙体,路面和屋顶的直射太阳平均通量为:
S↓为水平表面可获得的散射太阳辐射,街谷内表面接受到的散射太阳辐射量由天空视角系数直接得到,由于受到街谷的形状和高反照率的建筑表面材料的影响,短波辐射收支计算是求解一个具有无限多个反射面的几何系统,这些反射过程假设为等熵过程;
当每个表面的直射和散射反照率都是一样的,当第一次反射发生时,路面和墙所存储的能量为:
其中αr、αw分别表示路面和墙壁的反射率;
反射部分的能量Rr和Rw是:
在发生n次反射之后,
Ar(n+1)=Ar(n)+(1-αr)(1-Ψr)Rw(n)
Aw(n+1)=Aw(n)+(1-αr)ΨwRr(n)+(1-αw)(1-2Ψw)Rw(n)
Rr(n+1)=αr(1-Ψr)Rw(n)
Rw(n+1)=αwΨwRr(n)+αw(1-2Ψw)Rw(n)
递推公式就得到,
以及
对于这种无限次反射的情况,求解该几何系统可得,
设M作为路面和墙反射的总和,
其中,
每个表面类型吸收的总太阳辐射为:
可选的,所述人为产热具体包括:
当前街谷内的人为热通量QF=QFV+QFH+QFM;
其中,QFV、QFH和QFM分别为车辆、固定热源和生物代谢产生的热量。
可选的,所述显热通量QH包括:
QH,r,ww,we=ρCpCH1Ucan(Tr,ww,we-Tcan)
QH,R=ρCpCH2Utop(TR-Tair)
QH,can=ρCpCH2Uair(Tcan-Tair)
其中,r、ww、we和R分别指路面、西墙、东墙和屋顶;ρ为空气密度;cp为恒定压力下的比热;Tcan为街谷中心(w/2,h/2)的温度;Ucan和Utop为街谷中心(w/2,h/2)的风速和街谷上方的风速;
Uair和Tair为湍流模型参考高度的输入风速和温度,CH1和CH2为无量纲化的速度传递系数;两者的区别仅在于参考层的高度和粗糙度取值;采用相同的零平面层和粗糙度,两者值相等,采用如下方式计算,
其中,k为Von Karman常数,u*为参考层的摩擦速度,Ψh为通用积分函数,
其中,ζ'=(za-d)/L;ζT=zT/L,zT为热流的粗糙度长度;L为Obukhov稳定性长度,
其中,T为该层的平均温度,Ha是街谷与大气之间的空气通量,L为一个隐式函数,通过简化迭代的方式求解,在忽略城市冠层空气的比热时,Ha为街谷内墙壁通量和道路通量的加权平均,即有,
Ha=2(h/w)Qw+QR
在Masson提出的TEB模型中,CH*u*为空气动力学阻力的倒数,即1/RES*,通过对街谷内和顶部的风速情况来确定;
若不考虑绿地等植物覆盖的表面,平均街谷显热热流采用屋顶、墙壁与路面的面积加权平均值,
可选的,所述潜热通量QE包括:
在建筑屋顶与大气直接的潜热热流
QE,R=lvBRρCH2Utop(qR-qair)
其中,lv是蒸发潜热,BR是屋顶的湿度参数,0~1之间,0为完全干,1为完全湿,B的值取决于表面的植物和水体情况,qR为屋顶表面湿度;qair为参考高度处的湿度,
路面和墙面与街谷内的空气采用相似定律来计算潜热热流
QE,r=lvBrρCH1Ucan(qr-qcan)
QE,w=0
街谷内部与顶部大气之间的潜热热流
QE,can=lvρCH2Uair(qcan-qair)。
可选的,所述净储热通量ΔQS包括:
由于建筑或路面内部存在温度梯度,所以在假设屋顶、墙面和路面结构时都至少假定为3层结构,对最外部的一层结构来说,三种平面的传热方程写为,
相邻两层之间的平均导热率采用几何平均方法:
其中,λ*i为第i层的热导率;
对于表面的内部第一层,假设为非常薄的表面,把第一层的温度简化为外部的表面温度;对内部的其他的第i层,计算相邻层之间的导热量,对最内部的一层,如第n层,屋顶和墙壁表面可采用建筑的内部温度,而路面则可采用0通量;
Grn-n+1=0
假设建筑内部温度与外部街谷温度处于准稳态平衡状态,那么,若假设在热带海岛气候下,空调或自然通风下建筑内部的温度Tin基本恒定不变,将建筑内部中心的平均温度,
其中,b为建筑物的平均宽度。
可选的,所述风速包括:
在街谷内,风速分解为沿着墙壁的垂直速度Wcan,沿街道长度方向的水平速度Ucan;忽略沿街道宽度方向的水平速度;
根据观测,在接近街谷顶部的部分,不论街谷上方的空气稳定性及风向如何,垂直风速的标准方差σw等于摩擦速度u*;
接近屋顶的部分σw/u*为1.15,与观测结果是同一个数量级,对惯性边界层来说,u*的偏差不超过10%,所以,对任意的街谷高宽比,垂直速度都假定为,
其中,Uair是大气模型第一层的风速,Cd是曳力系数,通过街谷内及上方的温度/湿度、粗糙度Z0、及稳定性效应来计算获得;
街谷顶部的水平风速是通过Log近似曲线来获得的,Log曲线的处理范围是从屋顶下部h/3处开始到大气模型的第一层高度,当考虑到所有的街谷朝向,作360o的积分处理,那么街谷顶部的速度是,
其中,Δz是屋顶到大气模型第一层的高度;
水平风速Ucan是根据在街谷1/2高度处的风速确定的;
为了计算Ucan,需要假设一个合理的Ucan在垂直方向的变化规律;
根据风速的连续性假设,Ucan在垂直方向的变化曲线具有如下形式,
Ucan=Utopexp(-N/2)
其中,N的取值稍有差别;
根据街谷的高宽比(h/w=1~4),Ucan的值变化范围是0.75Utop [248]~0.4Utop;
N=0.5(h/w),街谷内的水平风速为,
简化街谷内路面和墙壁的空气动力学粗糙度计算,认为两者相等,且与街谷内外的稳定性无关,
其中,参数RESw和RESr,为CpCH1和CpCH2的倒数,用于计算显热和潜热热流。
可选的,所述净对流热通量为:
在准稳态条件下,街谷顶部沿街谷长度方向,沿街谷长度方向为x轴方向的风速流动已经达到稳定;若不考虑街道内部人和车辆走动的影响,在x方向上采用质量守恒和动量守恒定律,从而得到空气在街谷内部的水平运动情况;
如果街谷内的空气密度和水平速度参数都按照准稳态变量处理,那么在入口速度和出口速度满足出流条件的情况下,在x方向上,质量守恒和动量守恒定律按照一维流动方程写为:
其中,ρ是空气密度;是在x方向的平均速度;是进入或离开控制体的空气比体积质量流量,所述比体积为进出控制体的空气质量流量与控制体的体积之比;p为街谷横截面的平均压力;τw为墙壁表面和街道表面对空气流动的平均摩擦应力;和分别为街谷的入口和出口空气平均流速;x0为街谷的入口位置,街谷入口的流速通过仪器测量获得;
从区域规模的角度,街道通常都是连接起来,形成一个街道网络;城市热岛现象的水平通量可通过研究道路网来代替仅仅一条条独立存在的街道所形成的街谷;在计算十字路口的水平通量时,采用计算流体网络时的基尔霍夫原理;根据街道网络的拓扑结构和平面图论,可将街道网络用相应的邻接矩阵表示,求解每个支路的水平空气流量;
对于一个十字路口的场景,设街道m出口的水平通量为Qm,out,街道j出口的水平通量为Qj,out;街道i入口的水平通量为Qi,out,街道n入口的水平通量为Qn,out;那么,根据能量守恒定律,在节点处的水平通量为,最后根据能量守恒定律,在节点处的水平通量为:
混合后的水平通量Qmix为节点处流入街谷n和i的入口的水平通量,
本发明提供的热带海岛城市冠层模型,考虑了热带城市强烈的太阳辐射和高温高湿度气候特点和海岛城市受到常年季风的影响,在城市冠层模型的基础上改进了处理长波辐射通量、短波辐射通量、显热和潜热通量、街谷风速、热存储通量、人为热通量及水平热通量的方法,对于研究热带海岛型城市热岛效应的适配性更高。
上述说明仅是本发明技术方案的概述,为了能够更清楚了解本发明的技术手段,而可依照说明书的内容予以实施,并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂,以下特举本发明的具体实施方式。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图。
图1为本发明实施例提供的街谷冠层模型横截面图;
图2为本发明实施例提供的沿街道长度方向的空气流动分析图;
图3为本发明实施例提供的街谷网络在节点处的水平通量。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。
本发明的说明书实施例和权利要求书及附图中的术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元。
下面结合附图和实施例,对本发明的技术方案做进一步的详细描述。
热带海岛城市冠层模型的主要假设条件;
道路与其两边的建筑物相对面构成的空间称为街谷。
道路的长度远大于其宽度。街道的高宽比(aspect ratio)随道路的实际情况而改变。
道路的方向与太阳照射的夹角是任意的。
道路两边的建筑物的高度和宽度相同,屋顶高度为大气模型的表面层高度。
街谷内的空气流动和能量交换变化缓慢,以小时为时间分辨率,热岛效应的热力和动力过程可视为准稳态过程。
忽略城市热岛对外部大气效应影响的反馈。
采用2层模型来描述路面,墙壁和屋顶的热传导。
如图1所示,一种基于热带海岛气候特征的城市冠层模型建立方法,所述建立方法包括:
将相邻区域联系在一起,区域内的有限长度的多个街道及相互影响;
街道内的能量平衡方程如下:
Q* s+QF,s=QH,s+QE,s+ΔQS,s+ΔQA,s
其中,Q*为净辐射热通量,W/m2;QF为人为产热,W/m2;QH为显热通量,W/m2;QE为潜热通量,W/m2;ΔQS为净储热通量,W/m2;ΔQA为净对流热通量,W/m2;s为街道;
净辐射热通量Q*=长波辐射+短波辐射。
长波辐射包括:
路面的长波净辐射
墙面的长波净辐射
其中,L↓为太阳辐射量,∈r、∈w为路面和墙壁的发射率,对路面来说,若Ψr为路面对天空的视角系数,则(1-Ψr)为路面对两边墙壁的视角系数;墙壁对天空视角系数为Ψw,对路面的视角系数为Ψw,则对对面墙的视角系数为(1-2Ψw),对屋顶来说,天空视角为1.0;
采用平面角度计算天空视角系数,路面w/2位置的天空视角系数为
短波辐射包括:
根据街道与太阳方向的垂直角度计算路面、西墙、东墙和屋顶的平均太阳直接辐射通量是:
根据所述街谷的朝向变化,将所述街谷的宽度w修正为w/sinθ;在得到墙体热通量后,需要再乘上sinθ来修正,θ0为路面完全没有得到阳光直射的街道峡谷的朝向
墙体,路面和屋顶的直射太阳平均通量为:
S↓为水平表面可获得的散射太阳辐射,街谷内表面接受到的散射太阳辐射量由天空视角系数直接得到,由于受到街谷的形状和高反照率的建筑表面材料的影响,短波辐射收支计算是求解一个具有无限多个反射面的几何系统,这些反射过程假设为等熵过程;
当每个表面的直射和散射反照率都是一样的,当第一次反射发生时,路面和墙所存储的能量为:
其中αr、αw分别表示路面和墙壁的反射率;
反射部分的能量Rr和Rw是:
在发生n次反射之后,
Ar(n+1)=Ar(n)+(1-αr)(1-Ψr)Rw(n)
Aw(n+1)=Aw(n)+(1-αr)ΨwRr(n)+(1-αw)(1-2Ψw)Rw(n)
Rr(n+1)=αr(1-Ψr)Rw(n)
Rw(n+1)=αwΨwRr(n)+αw(1-2Ψw)Rw(n)
递推公式就得到,
以及
对于这种无限次反射的情况,求解该几何系统可得,
设M作为路面和墙反射的总和,
其中,
每个表面类型吸收的总太阳辐射为:
人为产热具体包括:
当前街谷内的人为热通量QF=QFV+QFH+QFM;
其中,QFV、QFH和QFM分别为车辆、固定热源和生物代谢产生的热量。
显热通量QH包括:
QH,r,ww,we=ρCpCH1Ucan(Tr,ww,we-Tcan)
QH,R=ρCpCH2Utop(TR-Tair)
QH,can=ρCpCH2Uair(Tcan-Tair)
其中,r、ww、we和R分别指路面、西墙、东墙和屋顶;ρ为空气密度;cp为恒定压力下的比热;Tcan为街谷中心(w/2,h/2)的温度;Ucan和Utop为街谷中心(w/2,h/2)的风速和街谷上方的风速;
Uair和Tair为湍流模型参考高度的输入风速和温度,CH1和CH2为无量纲化的速度传递系数;两者的区别仅在于参考层的高度和粗糙度取值;采用相同的零平面层和粗糙度,两者值相等,采用如下方式计算,
其中,k为Von Karman常数,u*为参考层的摩擦速度,Ψh为通用积分函数,
其中,ζ'=(za-d)/L;ζT=zT/L,zT为热流的粗糙度长度;L为Obukhov稳定性长度,
其中,T为该层的平均温度,Ha是街谷与大气之间的空气通量,L为一个隐式函数,通过简化迭代的方式求解,在忽略城市冠层空气的比热时,Ha为街谷内墙壁通量和道路通量的加权平均,即有,
Ha=2(h/w)Qw+QR
在Masson提出的TEB模型中,CH*u*为空气动力学阻力的倒数,即1/RES*,通过对街谷内和顶部的风速情况来确定;
若不考虑绿地等植物覆盖的表面,平均街谷显热热流采用屋顶、墙壁与路面的面积加权平均值,
可选的,所述潜热通量QE包括:
在建筑屋顶与大气直接的潜热热流
QE,R=lvBRρCH2Utop(qR-qair)
其中,lv是蒸发潜热,BR是屋顶的湿度参数,0~1之间,0为完全干,1为完全湿,B的值取决于表面的植物和水体情况,qR为屋顶表面湿度;qair为参考高度处的湿度,
路面和墙面与街谷内的空气采用相似定律来计算潜热热流
QE,r=lvBrρCH1Ucan(qr-qcan)
QE,w=0
街谷内部与顶部大气之间的潜热热流
QE,can=lvρCH2Uair(qcan-qair)。
净储热通量ΔQS包括:
由于建筑或路面内部存在温度梯度,所以在假设屋顶、墙面和路面结构时都至少假定为3层结构,对最外部的一层结构来说,三种平面的传热方程写为,
相邻两层之间的平均导热率采用几何平均方法:
其中,λ*i为第i层的热导率;
对于表面的内部第一层,假设为非常薄的表面,把第一层的温度简化为外部的表面温度;对内部的其他的第i层,计算相邻层之间的导热量,对最内部的一层,如第n层,屋顶和墙壁表面可采用建筑的内部温度,而路面则可采用0通量;
Grn-n+1=0
假设建筑内部温度与外部街谷温度处于准稳态平衡状态,那么,若假设在热带海岛气候下,空调或自然通风下建筑内部的温度Tin基本恒定不变,将建筑内部中心的平均温度,
其中,b为建筑物的平均宽度。
将建筑内部温度简化的好处是无需假设一个存在于建筑内部源项,由于空间加热或冷却系统带来的;而且,建筑内部的加热或冷却系统很难估计其功耗。这样,建筑内部的热通量存储都可统一采用温度边界条件处理。
风速包括:如图2所示,在街谷内,风速分解为沿着墙壁的垂直速度Wcan,沿街道长度方向的水平速度Ucan;忽略沿街道宽度方向的水平速度;
根据观测,在接近街谷顶部的部分,不论街谷上方的空气稳定性及风向如何,垂直风速的标准方差σw等于摩擦速度u*;
接近屋顶的部分σw/u*为1.15,与观测结果是同一个数量级,对惯性边界层来说,u*的偏差不超过10%,所以,对任意的街谷高宽比,垂直速度都假定为,
其中,Uair是大气模型第一层的风速,Cd是曳力系数,通过街谷内及上方的温度/湿度、粗糙度Z0、及稳定性效应来计算获得;
街谷顶部的水平风速是通过Log近似曲线来获得的,Log曲线的处理范围是从屋顶下部h/3处开始到大气模型的第一层高度,当考虑到所有的街谷朝向,作360°的积分处理,那么街谷顶部的速度是,
其中,Δz是屋顶到大气模型第一层的高度;
水平风速Ucan是根据在街谷1/2高度处的风速确定的;
为了计算Ucan,需要假设一个合理的Ucan在垂直方向的变化规律;
根据风速的连续性假设,Ucan在垂直方向的变化曲线具有如下形式,
Ucan=Utopexp(-N/2)
其中,N的取值稍有差别;
根据街谷的高宽比(h/w=1~4),Ucan的值变化范围是0.75Utop [248]~0.4Utop;N=0.5(h/w),街谷内的水平风速为,
简化街谷内路面和墙壁的空气动力学粗糙度计算,认为两者相等,且与街谷内外的稳定性无关,
其中,参数RESw和RESr,为CpCH1和CpCH2的倒数,用于计算显热和潜热热流。
净对流热通量为:
在准稳态条件下,街谷顶部沿街谷长度方向,沿街谷长度方向为x轴方向的风速流动已经达到稳定;若不考虑街道内部人和车辆走动的影响,在x方向上采用质量守恒和动量守恒定律,从而得到空气在街谷内部的水平运动情况;
如果街谷内的空气密度和水平速度参数都按照准稳态变量处理,那么在入口速度和出口速度满足出流条件的情况下,在x方向上,质量守恒和动量守恒定律按照一维流动方程写为:
其中,ρ是空气密度;是在x方向的平均速度;是进入或离开控制体的空气比体积质量流量,所述比体积为进出控制体的空气质量流量与控制体的体积之比;p为街谷横截面的平均压力;τw为墙壁表面和街道表面对空气流动的平均摩擦应力;和分别为街谷的入口和出口空气平均流速;x0为街谷的入口位置,街谷入口的流速通过仪器测量获得;
从区域规模的角度,街道通常都是连接起来,形成一个街道网络;城市热岛现象的水平通量可通过研究道路网来代替仅仅一条条独立存在的街道所形成的街谷;在计算十字路口的水平通量时,采用计算流体网络时的基尔霍夫原理;根据街道网络的拓扑结构和平面图论,可将街道网络用相应的邻接矩阵表示,求解每个支路的水平空气流量;
如图3所示,对于一个十字路口的场景,设街道m出口的水平通量为Qm,out,街道j出口的水平通量为Qj,out;街道i入口的水平通量为Qi,out,街道n入口的水平通量为Qn,out;那么,根据能量守恒定律,在节点处的水平通量为,最后根据能量守恒定律,在节点处的水平通量为:
混合后的水平通量Qmix为节点处流入街谷n和i的入口的水平通量,
有益效果:
随着城市不断高速发展,城市热岛效应日益严重。缓解城市热岛,有助于抑制传染病传播、减少温室气体排放、降低建筑能耗。为了研究热岛效应,需用城市冠层模型给区域建模,提高模型适配性是提高模型对热岛效应研究精确性的基础。
此处提出的新模型热带海岛城市冠层模型考虑了水平对流通量对热岛效应的影响,这弥补了经典的城市冠层模型对“热带海岛”城市的热岛效应研究的不足,对于“热带海岛”城市具有更高的适配性。另外,新模型基于有限长度的街谷进行建模,构造了一种三维空间的城市冠层模型,这种改进可大幅提高街谷内水平热通量的计算精度。
以上的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于热带海岛气候特征的城市冠层模型建立方法,其特征在于,所述建立方法包括:
将相邻区域联系在一起,区域内的有限长度的多个街道及相互影响;
街道内的能量平衡方程如下:
Q* s+QF,s=QH,s+QE,s+ΔQS,s+ΔQA,s 式1
其中,Q*为净辐射热通量,W/m2;QF为人为产热,W/m2;QH为显热通量,W/m2;QE为潜热通量,W/m2;ΔQS为净储热通量,W/m2;ΔQA为净对流热通量,W/m2;s为街道;
所述净辐射热通量Q*=长波辐射+短波辐射。
3.根据权利要求1所述的一种基于热带海岛气候特征的城市冠层模型建立方法,其特征在于,所述短波辐射包括:
根据街道与太阳方向的垂直角度计算路面、西墙、东墙和屋顶的平均太阳直接辐射通量是:
根据所述街谷的朝向变化,将所述街谷的宽度w修正为w/sinθ;在得到墙体热通量后,需要再乘上sinθ来修正,θ0为路面完全没有得到阳光直射的街道峡谷的朝向
墙体,路面和屋顶的直射太阳平均通量为:
S↓为水平表面可获得的散射太阳辐射,街谷内表面接受到的散射太阳辐射量由天空视角系数直接得到,由于受到街谷的形状和高反照率的建筑表面材料的影响,短波辐射收支计算是求解一个具有无限多个反射面的几何系统,这些反射过程假设为等熵过程;
当每个表面的直射和散射反照率都是一样的,当第一次反射发生时,路面和墙所存储的能量为:
其中αr、αw分别表示路面和墙壁的反射率;
反射部分的能量Rr和Rw是:
在发生n次反射之后,
Ar(n+1)=Ar(n)+(1-αr)(1-Ψr)Rw(n)
Aw(n+1)=Aw(n)+(1-αr)ΨwRr(n)+(1-αw)(1-2Ψw)Rw(n)
Rr(n+1)=αr(1-Ψr)Rw(n)
Rw(n+1)=αwΨwRr(n)+αw(1-2Ψw)Rw(n)
那么,递推公式就得到,
以及
对于这种无限次反射的情况,求解该几何系统可得,
设M作为路面和墙反射的总和,
其中,
每个表面类型吸收的总太阳辐射为:
4.根据权利要求1所述的一种基于热带海岛气候特征的城市冠层模型建立方法,其特征在于,所述人为产热具体包括:
当前街谷内的人为热通量QF=QFV+QFH+QFM;
其中,QFV、QFH和QFM分别为车辆、固定热源和生物代谢产生的热量。
5.根据权利要求1所述的一种基于热带海岛气候特征的城市冠层模型建立方法,其特征在于,所述显热通量QH包括:
QH,r,ww,we=ρCpCH1Ucan(Tr,ww,we-Tcan)
QH,R=ρCpCH2Utop(TR-Tair)
QH,can=ρCpCH2Uatr(Tcan-Tair)
其中,r、ww、we和R分别指路面、西墙、东墙和屋顶;ρ为空气密度;cp为恒定压力下的比热;Tcan为街谷中心(w/2,h/2)的温度;Ucan和Utop为街谷中心(w/2,h/2)的风速和街谷上方的风速;
Uair和Tair为湍流模型参考高度的输入风速和温度,CH1和CH2为无量纲化的速度传递系数;两者的区别仅在于参考层的高度和粗糙度取值;采用相同的零平面层和粗糙度,两者值相等,采用如下方式计算,
其中,k为Von Karman常数,u*为参考层的摩擦速度,Ψh为通用积分函数,
其中,ζ′=(za-d)/L;ζT=zT/L,zT为热流的粗糙度长度;L为Obukhov稳定性长度,
其中,T为该层的平均温度,Ha是街谷与大气之间的空气通量,L为一个隐式函数,通过简化迭代的方式求解,在忽略城市冠层空气的比热时,Ha为街谷内墙壁通量和道路通量的加权平均,即有,
Ha=2(h/w)Qw+QR
在Masson提出的TEB模型中,CH*u*为空气动力学阻力的倒数,即1/RES*,通过对街谷内和顶部的风速情况来确定;
若不考虑绿地等植物覆盖的表面,平均街谷显热热流采用屋顶、墙壁与路面的面积加权平均值,
6.根据权利要求1所述的一种基于热带海岛气候特征的城市冠层模型建立方法,其特征在于,所述潜热通量QE包括:
在建筑屋顶与大气直接的潜热热流
QE,R=lvBRρCH2Utop(qR-qair)
其中,lv是蒸发潜热,BR是屋顶的湿度参数,0~1之间,0为完全干,1为完全湿,B的值取决于表面的植物和水体情况,qR为屋顶表面湿度;qair为参考高度处的湿度,
路面和墙面与街谷内的空气采用相似定律来计算潜热热流
QE,r=lvBrρCH1Ucan(qr-qcan)
QE,w=0
街谷内部与顶部大气之间的潜热热流
QE,can=lvρCH2Uair(qcan-qair)。
7.根据权利要求1所述的一种基于热带海岛气候特征的城市冠层模型建立方法,其特征在于,所述净储热通量ΔQS包括:
由于建筑或路面内部存在温度梯度,所以在假设屋顶、墙面和路面结构时都至少假定为3层结构,对最外部的一层结构来说,三种平面的传热方程写为,
相邻两层之间的平均导热率采用几何平均方法:
其中,λ*i为第i层的热导率;
对于表面的内部第一层,假设为非常薄的表面,把第一层的温度简化为外部的表面温度;对内部的其他的第i层,计算相邻层之间的导热量,对最内部的一层,如第n层,屋顶和墙壁表面可采用建筑的内部温度,而路面则可采用0通量;
Grn-n+1=0
假设建筑内部温度与外部街谷温度处于准稳态平衡状态,那么,若假设在热带海岛气候下,空调或自然通风下建筑内部的温度Tin基本恒定不变,将建筑内部中心的平均温度,
其中,b为建筑物的平均宽度。
8.根据权利要求1所述的一种基于热带海岛气候特征的城市冠层模型建立方法,其特征在于,所述风速包括:
在街谷内,风速分解为沿着墙壁的垂直速度Wcan,沿街道长度方向的水平速度Ucan;忽略沿街道宽度方向的水平速度;
根据观测,在接近街谷顶部的部分,不论街谷上方的空气稳定性及风向如何,垂直风速的标准方差σw等于摩擦速度u*;
接近屋顶的部分σw/u*为1.15,与观测结果是同一个数量级,对惯性边界层来说,u*的偏差不超过10%,所以,对任意的街谷高宽比,垂直速度都假定为,
其中,Uair是大气模型第一层的风速,Cd是曳力系数,通过街谷内及上方的温度/湿度、粗糙度Z0、及稳定性效应来计算获得;
街谷顶部的水平风速是通过Log近似曲线来获得的,Log曲线的处理范围是从屋顶下部h/3处开始到大气模型的第一层高度,当考虑到所有的街谷朝向,作360°的积分处理,那么街谷顶部的速度是,
其中,Δz是屋顶到大气模型第一层的高度;
水平风速Ucan是根据在街谷1/2高度处的风速确定的;
为了计算Ucan,需要假设一个合理的Ucan在垂直方向的变化规律;
根据风速的连续性假设,Ucan在垂直方向的变化曲线具有如下形式,
Ucan=Utopexp(-N/2)
其中,N的取值稍有差别;
根据街谷的高宽比(h/w=1~4),Ucan的值变化范围是0.75Utop [248]~0.4Utop;
N=0.5(h/w),街谷内的水平风速为,
简化街谷内路面和墙壁的空气动力学粗糙度计算,认为两者相等,且与街谷内外的稳定性无关,
其中,参数RESw和RESr,为CpCH1和CpCH2的倒数,用于计算显热和潜热热流。
9.根据权利要求1所述的一种基于热带海岛气候特征的城市冠层模型建立方法,其特征在于,所述净对流热通量为:
在准稳态条件下,街谷顶部沿街谷长度方向,沿街谷长度方向为x轴方向的风速流动已经达到稳定;若不考虑街道内部人和车辆走动的影响,在x方向上采用质量守恒和动量守恒定律,从而得到空气在街谷内部的水平运动情况;
如果街谷内的空气密度和水平速度参数都按照准稳态变量处理,那么在入口速度和出口速度满足出流条件的情况下,在x方向上,质量守恒和动量守恒定律按照一维流动方程写为:
其中,ρ是空气密度;是在x方向的平均速度;是进入或离开控制体的空气比体积质量流量,所述比体积为进出控制体的空气质量流量与控制体的体积之比;p为街谷横截面的平均压力;τw为墙壁表面和街道表面对空气流动的平均摩擦应力;和分别为街谷的入口和出口空气平均流速;x0为街谷的入口位置,街谷入口的流速通过仪器测量获得;
从区域规模的角度,街道通常都是连接起来,形成一个街道网络;城市热岛现象的水平通量可通过研究道路网来代替仅仅一条条独立存在的街道所形成的街谷;在计算十字路口的水平通量时,采用计算流体网络时的基尔霍夫原理;根据街道网络的拓扑结构和平面图论,可将街道网络用相应的邻接矩阵表示,求解每个支路的水平空气流量;
对于一个十字路口的场景,设街道m出口的水平通量为Qm,out,街道j出口的水平通量为Qj,out;街道i入口的水平通量为Qi,out,街道n入口的水平通量为qn,out;那么,根据能量守恒定律,在节点处的水平通量为,最后根据能量守恒定律,在节点处的水平通量为:
混合后的水平通量Qmix为节点处流入街谷n和i的入口的水平通量,
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111236079.4A CN113987779A (zh) | 2021-10-22 | 2021-10-22 | 一种基于热带海岛气候特征的城市冠层模型建立方法 |
US17/873,185 US20230129704A1 (en) | 2021-10-22 | 2022-07-26 | Method for building urban canopy model based on tropical island climate characteristics |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111236079.4A CN113987779A (zh) | 2021-10-22 | 2021-10-22 | 一种基于热带海岛气候特征的城市冠层模型建立方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113987779A true CN113987779A (zh) | 2022-01-28 |
Family
ID=79740616
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111236079.4A Pending CN113987779A (zh) | 2021-10-22 | 2021-10-22 | 一种基于热带海岛气候特征的城市冠层模型建立方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20230129704A1 (zh) |
CN (1) | CN113987779A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115146564A (zh) * | 2022-08-02 | 2022-10-04 | 上海市气候中心(上海区域气候中心) | 基于垂直分层降尺度技术的城市地面风速精细化模拟方法 |
CN116720416A (zh) * | 2023-08-09 | 2023-09-08 | 武汉大学 | 基于三维模拟的城市通风性和污染物扩散的分析方法 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116843193B (zh) * | 2023-06-14 | 2024-02-20 | 华南理工大学 | 一种基于通量观测数据的人为热通量估算方法及系统 |
CN116882034B (zh) * | 2023-09-06 | 2023-11-17 | 武汉大学 | 基于三维模拟的城市立体绿化分布方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101320098A (zh) * | 2008-07-11 | 2008-12-10 | 重庆大学 | 基于数字图像分析的城镇热岛特性预测方法及系统 |
-
2021
- 2021-10-22 CN CN202111236079.4A patent/CN113987779A/zh active Pending
-
2022
- 2022-07-26 US US17/873,185 patent/US20230129704A1/en active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101320098A (zh) * | 2008-07-11 | 2008-12-10 | 重庆大学 | 基于数字图像分析的城镇热岛特性预测方法及系统 |
Non-Patent Citations (1)
Title |
---|
周海珠: "热带海岛城市热岛效应研究与城市绿地布局优化", 万方数据库, 15 July 2021 (2021-07-15), pages 44 - 66 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115146564A (zh) * | 2022-08-02 | 2022-10-04 | 上海市气候中心(上海区域气候中心) | 基于垂直分层降尺度技术的城市地面风速精细化模拟方法 |
CN116720416A (zh) * | 2023-08-09 | 2023-09-08 | 武汉大学 | 基于三维模拟的城市通风性和污染物扩散的分析方法 |
CN116720416B (zh) * | 2023-08-09 | 2023-10-17 | 武汉大学 | 基于三维模拟的城市通风性和污染物扩散的分析方法 |
Also Published As
Publication number | Publication date |
---|---|
US20230129704A1 (en) | 2023-04-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113987779A (zh) | 一种基于热带海岛气候特征的城市冠层模型建立方法 | |
Antoniou et al. | CFD simulation of urban microclimate: Validation using high-resolution field measurements | |
Kanda et al. | A simple energy balance model for regular building arrays | |
Antoniou et al. | CFD and wind-tunnel analysis of outdoor ventilation in a real compact heterogeneous urban area: Evaluation using “air delay” | |
He et al. | High-resolution dataset of urban canopy parameters for Beijing and its application to the integrated WRF/Urban modelling system | |
Zhao et al. | Isothermal and non-isothermal flow in street canyons: A review from theoretical, experimental and numerical perspectives | |
van Hooff et al. | CFD evaluation of natural ventilation of indoor environments by the concentration decay method: CO2 gas dispersion from a semi-enclosed stadium | |
Piroozmand et al. | Coupled CFD framework with mesoscale urban climate model: Application to microscale urban flows with weak synoptic forcing | |
Van Hooff et al. | Coupled urban wind flow and indoor natural ventilation modelling on a high-resolution grid: A case study for the Amsterdam ArenA stadium | |
Van Hooff et al. | On the effect of wind direction and urban surroundings on natural ventilation of a large semi-enclosed stadium | |
Martilli et al. | An urban surface exchange parameterisation for mesoscale models | |
Otte et al. | Implementation of an urban canopy parameterization in a mesoscale meteorological model | |
Asawa et al. | Thermal design tool for outdoor spaces based on heat balance simulation using a 3D-CAD system | |
Yamada et al. | Downscaling mesoscale meteorological models for computational wind engineering applications | |
Gülten et al. | Influence of trees on heat island potential in an urban canyon | |
He et al. | Measurement and evaluation of the summer microclimate in the semi-enclosed space under a membrane structure | |
Benzaama et al. | Study of the effect of sun patch on the transient thermal behaviour of a heating floor in Algeria | |
M’Saouri El Bat et al. | Integration of a practical model to assess the local urban interactions in building energy simulation with a street canyon | |
Sedighi et al. | Thermal study of domed roofs in a traditional bazaar (the case of old Ganj-Alikhan bazaar in Kerman, Iran) | |
Hadavand et al. | Thermal analysis of vaulted roofs | |
Mateus et al. | Natural ventilation of large air masses: experimental and numerical techniques review | |
Xiong et al. | Impacts of uneven surface heating of an ideal street canyon on airflows and indoor ventilation: Numerical study using OpenFOAM coupled with EnergyPlus | |
Piotrowski et al. | Mathematical model of the thermal-air regime of a ventilated attic | |
Kawai et al. | Validation of a numerical model for urban energy‐exchange using outdoor scale‐model measurements | |
Najafi et al. | Thermal study of a cistern's dome (the case of Motamed cistern in Lar, Iran) |
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 |