CN114386215A - 一种建筑物可分辨尺度的城市三维风场快速分析方法 - Google Patents

一种建筑物可分辨尺度的城市三维风场快速分析方法 Download PDF

Info

Publication number
CN114386215A
CN114386215A CN202111367080.0A CN202111367080A CN114386215A CN 114386215 A CN114386215 A CN 114386215A CN 202111367080 A CN202111367080 A CN 202111367080A CN 114386215 A CN114386215 A CN 114386215A
Authority
CN
China
Prior art keywords
building
wind field
wind
area
urban
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
CN202111367080.0A
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.)
Shanghai Meteorological Information And Technical Support Center Shanghai Meteorological Archives
Shanghai Institute Of Typhoon China Meteorological Administration
Original Assignee
Shanghai Meteorological Information And Technical Support Center Shanghai Meteorological Archives
Shanghai Institute Of Typhoon China Meteorological Administration
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 Shanghai Meteorological Information And Technical Support Center Shanghai Meteorological Archives, Shanghai Institute Of Typhoon China Meteorological Administration filed Critical Shanghai Meteorological Information And Technical Support Center Shanghai Meteorological Archives
Priority to CN202111367080.0A priority Critical patent/CN114386215A/zh
Publication of CN114386215A publication Critical patent/CN114386215A/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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明公开了一种建筑物可分辨尺度的城市三维风场快速分析方法，其特征在于，包括：获取目标区域气象数据，气象数据包括气象数据时间和入流风信息；获取目标区域建筑物数据，并确认各个建筑物在模型网格内的尺寸、相对位置及每个建筑物的迎风区、背风区和顶部涡旋区的尺度；根据气象数据和建筑物数据构建初始风场；根据质量守恒方程调整初始风场，得到最终风场。本发明通过利用简单的气象输入以及基本的建筑物信息能够快速模拟出超大城市中超高层建筑物对周围风场环境的影响，计算量小，同时加入坡度因素，模拟结果更加准确。

Description

一种建筑物可分辨尺度的城市三维风场快速分析方法
技术领域
本发明涉及风环境分析技术领域，尤其涉及一种建筑物可分辨尺度的城市三维风场快速分析方法。
背景技术
随着城市化水平的提高以及城市建筑密度、人口密度、机动车保有量等不断加大，城市区域局地风灾环境以及城市冠层内风场与污染扩散等逐渐成为人们日益关注的问题。研究城市冠层内风环境对城市建设、城市风灾、城市空气质量改善均具有重要作用。
过去对城市风的研究，一般忽略城市建筑物的影响，然而世界经济的发展催生大量的高层建筑，给高层建筑带来了一些列的风灾事故，这使得人们不得不重视建筑周围风环境会受到建筑影响的问题。建筑周围的风环境与建筑单体的外形、尺寸、建筑间距及周围的地形地貌等因素都有密切关系，由于建筑物钝体的阻挡作用，造成建筑物周围风场在时间和空间上都具有非常复杂的非定常流性状。如何更好地观测和模拟城市风环境，一直是城市微气候研究的难题。
当前研究建筑周围风环境的方法分为实地测试、物理模拟和计算机数值模拟。三种方法各有优缺点，在现今的研究或应用当中均有所体现。总体而言，当前研究方法共同存在的问题是研究成本(包括计算成本和实验成本)较高，无法长时间模拟或快速模拟建筑物周围风场。在一些特殊情形下，比如突发事故，三种研究方法无法满足实际研究需求。并且城市建筑物中还有大量为斜坡型屋顶结构，在目前的城市风模拟中未充分考虑坡度影响。因此，针对上述问题，有必要提出进一步地解决方案。
发明内容
本发明旨在提供一种建筑物可分辨尺度的城市三维风场快速分析方法，以克服现有技术中存在的不足。
为解决上述技术问题，本发明的技术方案是：
一种建筑物可分辨尺度的城市三维风场快速分析方法，包括：
获取目标区域气象数据，所述气象数据包括气象数据时间和入流风信息；
获取目标区域建筑物数据，并确认各个建筑物在模型网格内的尺寸、相对位置及每个建筑物的迎风区、背风区和顶部涡旋区的尺度；
根据气象数据和建筑物数据构建初始风场；
根据质量守恒方程调整所述初始风场，得到最终风场。
本发明的一个较佳实施例中，所述质量守恒方程中包括斜坡前方、后方、右侧及左侧边界条件系数，所述斜坡包括建筑物倾斜屋顶及地表斜坡。
本发明的一个较佳实施例中，在所述质量守恒方程中对于底面及建筑物采用第二类边界条件
Figure RE-RE-GDA0003464590630000021
其他边界采用第一类自由边界条件λ＝0。
本发明的一个较佳实施例中，所述根据质量守恒方程调整所述初始风场，得到最终风场，包括：
采用变分方法求解所述质量守恒方程，且α1＝α2＝0.5。
本发明的一个较佳实施例中，所述建筑物数据包括屋顶是否为坡度、建筑物高度、建筑物长度及建筑物最长边和正东方向夹角角度。
本发明的一个较佳实施例中，所述确认各个建筑物在模型网格内的尺寸，包括：
当入流风和建筑物的夹角θ不超过20°时，建筑物在模型网格内的尺寸Weff、Leff与其实际尺寸W、L相等；
当入流风和建筑物的夹角大于20°时，Weff＝Wsinθ+Lcosθ，Leff＝Wcosθ+Lsinθ；
其中，W为宽度，L为长度。
本发明的一个较佳实施例中，所述根据气象数据和建筑物数据构建初始风场，包括：
计算空旷区域和建筑物区域的入流风的垂直廓线；
对单体建筑物周围风场参数化；
对建筑街渠风场参数化。
本发明的一个较佳实施例中，对单体建筑物周围风场参数化，包括：
当入流风和建筑物的夹角θ大于20°时，采用迎风面的入流风速分解为平行迎风面分量和垂直迎风面分量进行插值。
本发明的一个较佳实施例中，对建筑街渠风场参数化，包括：
当S/H>2.5或1.4<S/H<2.4时，采用单体建筑物周围风场参数化方案；
当S/H<1.4时，采用SR参数化方案或MR参数化方案；
其中，H为建筑物高度，S为建筑间距离。
与现有技术相比，本发明的有益效果是：
本发明通过利用简单的气象输入以及基本的建筑物信息能够快速模拟出超大城市中超高层建筑物对周围风场环境的影响，引入的参数化方案能够简化计算过程、减少计算量，同时加入坡度因素，使模拟结果更加准确。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明流程图；
图2为建筑物顶部涡旋区示意图；
图3为另一方案中建筑物顶部涡旋区示意图；
图4为入流风倾斜入射夹角示意图；
图5为建筑物等效长宽示意图；
图6为本发明的模拟网格图；
图7为本发明实施例中模拟区域建筑物分布图；
图8为本发明实施例中2-140m高度水平风场；
图9a为不同平面CEDVAL风洞观测，图9b为本方法模拟风矢量图；
图10为不同垂直平面全风速廓线比较图，三角为CEDVAL风洞观测数据，实线为本方法模式结果；
图11为不同垂直平面u、w风速廓线比较图，三角为CEDVAL 风洞观测数据u，实线为本方法模式结果u，圆圈为CEDVAL风洞观测数据w，虚线为本方法模式结果w。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的其他方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。
如图1所示，一种建筑物可分辨尺度的城市三维风场快速分析方法，包括：
S100获取目标区域气象数据，所述气象数据包括气象数据时间和入流风信息。
S200获取目标区域建筑物数据，并确认各个建筑物在模型网格内的尺寸、相对位置及每个建筑物的迎风区、背风区和顶部涡旋区的尺度。
其中，所述建筑物数据包括屋顶是否为坡度、建筑物高度、建筑物长度及建筑物最长边和正东方向夹角角度，以在模拟网格中构建目标区域建筑物分布。
S300根据气象数据和建筑物数据构建初始风场。
S400根据质量守恒方程调整所述初始风场，得到最终风场。
具体而言，S300根据气象数据和建筑物数据构建初始风场，包括：
S310计算空旷区域和建筑物区域的入流风的垂直廓线，如下：
(1)空旷区域
Figure RE-RE-GDA0003464590630000061
(2)建筑物区域
Figure RE-RE-GDA0003464590630000062
Figure RE-RE-GDA0003464590630000063
在(1)-(3)中，Hcan表示冠层高度，Ucan表示冠层顶Hcan风速大小，z0表示粗糙系数，一般约为0.1～0.2Hcan，根据一般城市冠层的粗糙度特征，以及现有城市模式的模拟经验，本申请取 z0＝0.1Hcan，d表示位移高度，一般约为0.7Hcan，本申请取d＝0.67Hcan。α(z)为衰减系数，为了简化模型，本申请将α(z)设为常数，α的范围一般为1～3，本申请取α＝1。α(z)和建筑物密度λf成正比，α＝kλf，其中k＝9.6。
S320对单体建筑物周围风场参数化，如下：
单体建筑物周围风场可分以下三个关键区域：迎风区(包括迎风位移区和迎风涡旋区)，尾流区(包括空腔区)，顶部涡旋区。
S321建筑物前迎风区
当入流风垂直于建筑物迎风墙面时，会在建筑物前迎风区域内形成位移区和迎风涡旋区。本模式采用的建筑物迎风区插值方案如下
Figure RE-RE-GDA0003464590630000071
Figure RE-RE-GDA0003464590630000072
Figure RE-RE-GDA0003464590630000073
Figure RE-RE-GDA0003464590630000074
在(4)-(7)中，H表示建筑物高度，L表示建筑物迎风方向的长度，W表示建筑物垂直于迎风方向的长度。(4)(5)公式左项表示Lfx、Lfx1与H的比值。从水平方向上看，Lfx表示迎风区(Upwind Zone) 长度，Lfx1表示涡旋区长度，迎风区和涡旋区的空间大小近似看作 1/4椭球，椭球公式为(6)(7)，Lfx即迎风区的椭球长轴，Lfx1即迎风区内涡旋区的椭球长轴，椭球z方向的半轴为0.6H，y方向上半轴为0.5W。椭球的原点为建筑物迎风面底边中点。假设迎风区内涡旋区为A区域，迎风区为B区域。
在涡旋区内，即A区域，当位置坐标(x,y,z)满足
Figure RE-RE-GDA0003464590630000081
水平和垂直方向风速的插值方案如下：
Figure RE-RE-GDA0003464590630000082
Figure RE-RE-GDA0003464590630000083
其中xf、zf表示x正方向、z正方向的位置，xf＝x，zf＝z。
位移区指B-A区域，当位置坐标(x,y,z)满足
Figure RE-RE-GDA0003464590630000084
且不在A区域内，水平方向风速的插值方案如下：
u0(z)＝CdzU(z) (10)
其中U(z)为S310中插值风速，Cdz是一个衰减系数，根据一般城市冠层的粗糙度特征，本申请取Cdz＝0.4。
S322建筑物尾流区
当入流风垂直于建筑物迎风墙面时，会在建筑物背后形成尾流区。本研究提供了两种背风区插值方案以供选择，通过风洞实验资料与模拟结果的对比试验，挑选合适方案。
第一种背风区插值方案
Figure RE-RE-GDA0003464590630000091
LW＝3LC (12)
Figure RE-RE-GDA0003464590630000092
Figure RE-RE-GDA0003464590630000093
在(11)-(14)中，H表示建筑物高度，L表示建筑物迎风方向的长度，W表示建筑物垂直于迎风方向的长度。(11)公式左项表示Lc与H的比值。从水平方向上看，Lw表示尾流区(Wake Zone) 长度，Lc表示空腔区(Cavity Zone)长度，尾流区和空腔区的空间大小近似看作1/4椭球，椭球公式为(13)(14)，Lw即背风区的椭球长轴，Lc即背风区内空腔区的椭球长轴，椭球z方向的半轴为 H，y方向上半轴为0.5W。椭球的原点为建筑物背风面底边中点。假设尾流区内空腔区为A区域，尾流区为B区域。
在空腔区内，即A区域，当位置坐标(x,y,z)满足
Figure RE-RE-GDA0003464590630000094
水平方向风速的插值方案如下：
Figure RE-RE-GDA0003464590630000095
Figure RE-RE-GDA0003464590630000096
在B-A区域，当位置坐标(x,y,z)满足
Figure RE-RE-GDA0003464590630000101
且不在A区域内，水平方向风速的插值方案如下。
Figure RE-RE-GDA0003464590630000102
在(15)(17)中，U(H)和U(z)表示建筑物顶和z高度的边界层风速。
第二种背风区插值方案
将尾流边界分为垂直区域和水平区域。
垂直区域指原尾流区的正上方，位置坐标(x,y,z)处的风速插值方案如下。
u(x,y,z)＝(URef_vert-U(z))exp(-ξ1.5)+U(z) (18)
Figure RE-RE-GDA0003464590630000103
Figure RE-RE-GDA0003464590630000104
Figure RE-RE-GDA0003464590630000105
(18)中U(z)表示z高度边界层风速，ξ是一个无量纲参数。(21) 中Ht_ellipse(x,y)表示当前位置的椭球高度，本申请中 Ht_ellipse(x,y)＝z_vert，因此(21)可以转化为
Figure RE-RE-GDA0003464590630000106
水平区域指原尾流区的两侧，位置坐标(x,y,z)处的风速插值方案如下。
u(x,y,z)＝(URef_hor-U(z))exp(-η2)+U(z) (23)
Figure RE-RE-GDA0003464590630000111
Figure RE-RE-GDA0003464590630000112
Figure RE-RE-GDA0003464590630000113
其中，U(z)表示z高度边界层风速，η是一个无量纲参数。
S323建筑物顶部涡旋区
本研究提供了两种顶部涡旋区插值方案以供选择，通过风洞实验资料与模拟结果的对比试验，挑选合适方案。
第一种顶部涡旋区插值方案
将建筑物顶部的插值区域看作椭圆柱体，长度和高度由建筑物的几何性质确定，尺度用R来表示。
Figure RE-RE-GDA0003464590630000114
HCM(Height of vortex)＝0.22*R (28)
LC(Length of vortex)＝0.9*R (29)
BS＝Smaller of upwind building height or width (30)
BL＝Larger of upwind building height or width (31)
如图2所示，建筑物顶部涡旋区可分为两个区域，区域1也被成为逆流区，采用调整的负对数风速插值方案，区域2采用正对数风速插值方案。
区域1
u(z)＝-U(z)*C1 (32)
Figure RE-RE-GDA0003464590630000121
区域2
u(z)＝U(z) (34)
(27)-(34)中H是建筑物高度，HCM是插值区的高度，U(z) 表示z高度边界层风速。
第二种顶部涡旋区插值方案
如图3所示，逆流区的风廓线形状应为倒楔形。回流区内采用高度线性垂直变化的插值公式。
Figure RE-RE-GDA0003464590630000122
z′＝H+HCM-z (36)
其中，H是建筑物高度，HCM是插值区(回流区)的高度， U(z)表示z高度的来流边界层风速。
注意的是需要考虑夹角入流情形。如图4所示，当入流风和建筑物的夹角θ过大时(一般认为θ>20°)，即当θ>20°时，迎风面的入流风速分解为平行迎风面分量u0(z)//和垂直迎风面分量u0(z)⊥， u0(z)⊥代替入流风速根据S321、S322和S323有关参数化方案进行插值，u0(z)//保持不变。
如图5所示，背部空腔区尾流区插值时，用建筑等效长度Leff和等效宽度Weff替代原有长宽。
Weff＝Wsinθ+Lcosθ (37)
Leff＝Wcosθ+Lsinθ (38)
综上，当入流风和建筑物的夹角过大时，建筑物前后的插值方案如表1、表2所示。
表1夹角入流建筑物位移区插值
Figure RE-RE-GDA0003464590630000131
表2夹角入流建筑物空腔区、尾流区插值
Figure RE-RE-GDA0003464590630000132
S330对建筑街渠风场参数化
当S/H>2.5或1.4<S/H<2.4时，采用单体建筑物周围风场参数化方案；
当S/H<1.4时，采用SR参数化方案或MR参数化方案。
具体而言，SR参数化方案(标准Rockle插值方案，Standard Rockle streetcanyon model)采用空间距离参数S*和比例形态参数W/H来确定街渠流场类型，S*的计算公式为
Figure RE-RE-GDA0003464590630000141
当S<S*，建筑街渠内形成一个顺时针涡旋，水平方向和垂直方向的速度插值公式为
Figure RE-RE-GDA0003464590630000142
Figure RE-RE-GDA0003464590630000143
上式中S为两个建筑物之间距离，d为位置坐标到上风建筑物墙面的距离，U(H)表示上风向建筑物顶的边界层风速。
MR参数化方案(改进的Rockle插值方案，Modified Rockle street canyonmodel)重新计算空间距离参数S*。
Figure RE-RE-GDA0003464590630000144
并将街渠流场分为三个区域：(1)中心区域，风速插值方案和传统街渠插值方案相同；(2)垂直湍流扩算区域，需考虑街渠上方的动量输入和输出；(3)水平湍流扩散区域，需考虑区域(1)水平两侧的动量输入和输出。区域(2)(3)的空间形态是由入流方向建筑物延伸出来的三角柱体或楔形体。区域(1)(2)(3)的边长分别为
δSC(xcan)＝W-2δlw(xcan) (43)
δvw(xcan)＝0.2xcan (44)
δlw(xcan)＝0.2xcan (45)
区域(1)的风速参数化方案为
Figure RE-RE-GDA0003464590630000151
Figure RE-RE-GDA0003464590630000152
区域(2)的风速参数化方案为
Figure RE-RE-GDA0003464590630000153
区域(3)的风速参数化方案为
Figure RE-RE-GDA0003464590630000154
Ylw＝[|ycan|-δSC(xcan)/2] (50)
上式中，xcan表示位置坐标到上风建筑物墙面的距离，ycan表示位置坐标到上风建筑物墙面中心平面的距离，p为经验参数，用于调整风洞实验结果与模拟结果的差异，此处取0.25，Uroof表示建筑物顶风速，ubl(z)表示z高度下边界层风速，γ为经验衰减系数，可根据观测资料或者风洞实验结果决定，此处取0.3。
S400根据质量守恒方程调整所述初始风场，得到最终风场。
由于上文所用风场参数化方案不能反映建筑物墙壁对气流的作用，因此在插值区域(特别是建筑物周围)不满足质量守恒连续方程。质量守恒是指速度的散度为0，假设插值风场为U0(u0,v0,w0)，有
Figure RE-RE-GDA0003464590630000155
Figure RE-RE-GDA0003464590630000156
本研究采用变分方法求解满足(51)的最终风场U(u,v,w)。首先，以U(u,v,w)和U0(u0,v0,w0)构建的变分函数公式
Figure RE-RE-GDA0003464590630000161
(52)中，α1和α2定义为高斯精度模，用于控制变分公式对水平 u、v和垂直w风速分量的调整力度。α1和α2的值越大，说明变分公式对速度的调整幅度较小，即相对于初始风速的改变较小。
引入拉格朗日乘数λ，联立(51)(52)可得
Figure RE-RE-GDA0003464590630000162
根据欧拉-拉格朗日方程组，U(u,v,w)和U0(u0,v0,w0)满足以下条件时，E(u,v,w,λ)取极小值。
Figure RE-RE-GDA0003464590630000163
Figure RE-RE-GDA0003464590630000164
Figure RE-RE-GDA0003464590630000165
将(54)代入(51)，得到关于λ的泊松方程。
Figure RE-RE-GDA0003464590630000166
α1和α2取常数，α1＝α2＝0.5。根据分辨率dx、dy、dz，以及计算区域长宽高的格点数nx、ny、nz，生成nx*ny*nz大小的三维网格。计算效率采用1/2网格，λi,j,k表示(i,j,k)网格的值，速度则表示为网格边界值(图6)。对(55)按照图6网格离散化，对于非边界格点，有
Figure RE-RE-GDA0003464590630000171
等式两边简化后，有
Figure RE-RE-GDA0003464590630000172
令
Figure RE-RE-GDA0003464590630000173
并定义边界系数为e，f，g，h， m，n，o，p，q。(57)可简化为
Figure RE-RE-GDA0003464590630000174
整理后得到
Figure RE-RE-GDA0003464590630000175
其中
Figure RE-RE-GDA0003464590630000176
采用超松弛算法(SOR)，根据(60)，有
Figure RE-RE-GDA0003464590630000181
ω为SOR收敛速度参数，本申请ω取1.78。
当
Figure RE-RE-GDA0003464590630000182
时，根据(62)求得λ的最优解。ε是一个小量，可按照计算速度要求进行调整，为了提高计算效率，加快收敛速度，本申请ε取0.01。
对于边界格点，离散化时需要考虑边界条件。本研究对于底面及建筑物采用第二类边界条件
Figure RE-RE-GDA0003464590630000183
其他边界采用第一类自由边界条件λ＝0。
比如，若正i方向存在建筑物，即建筑物墙壁位于i+1/2处，有
Figure RE-RE-GDA0003464590630000184
根据第二类边界条件，此时
Figure RE-RE-GDA0003464590630000185
Figure RE-RE-GDA0003464590630000186
代入(55)并按照图6网格离散化，有
Figure RE-RE-GDA0003464590630000187
简化后得
Figure RE-RE-GDA0003464590630000188
此时，o＝0.5，p＝1，q＝1，e＝0，f＝1，g＝1，h＝1，m＝1，n＝1。同理可得其他边界条件系数。
边界条件系数在质量守恒方程中，代表当前格点与建筑物墙面边界的相对位置。本研究将相对位置分为14类，详见表3。
表3边界条件系数表
Figure RE-RE-GDA0003464590630000191
本实施例中，具体对国内城市某地部分街区进行了模拟，建筑物信息来自地理信息系统(GIS)软件，模拟后的建筑物分布图如图7。通过以上方法计算出对应高度下水平面的格点风场，k表示输出的高度为2/4/6/8/10/20/40/.../140m，通过matlab画图程序read_wind.m，得到图8。
本实施例中气象数据来自欧洲数值预报中心(ECWMF)，每日08 时、20时预报结果，每一次预报53个时次，预报时长为10天，在 0-72hr内每3hr一个预报结果，72hr-240hr每6hr一个预报结果。
将本方法模拟结果与风洞试验结果进行比较，这里用的风洞数据来自于CEDVAL数据集。通过风洞观测和本方法模拟结构比较，如图9所示，本方法较准确地模拟出街区峡谷风场的关键区域。选取五个不同垂直剖面(X/H＝-2.58，X/H＝-0.5，X/H＝0，X/H＝0.5和 X/H＝2.58)的风洞观测和方法模拟结果，比较全风速U
Figure RE-RE-GDA0003464590630000201
u分量、w分量的垂直变化廓线(图10和图11)。从图中看出，廓线街渠中心区域(X/H＝0.5)变化幅度较大，建筑物高度以下的模式模拟结果存在一定高估，U和w分量的模拟误差较大。分析原因可能是由于街渠中心区域受到气流汇流作用，风速迅速增大，加之墙体对风速的阻挡和拖曳使得此区域风速梯度很大，增加了模拟难度，从而导致模拟风速偏大较多。X/H＝-2.58，X/H＝-0.5，X/H＝0.5和X/H＝2.58四个垂直平面U的模拟效果较好。
综上所述，本发明通过利用简单的气象输入以及基本的建筑物信息能够快速模拟出超大城市中超高层建筑物对周围风场环境的影响，计算量小，同时加入坡度因素，模拟结果更加准确。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。

Claims (9)

1.一种建筑物可分辨尺度的城市三维风场快速分析方法，其特征在于，包括：
获取目标区域气象数据，所述气象数据包括气象数据时间和入流风信息；
获取目标区域建筑物数据，并确认各个建筑物在模型网格内的尺寸、相对位置及每个建筑物的迎风区、背风区和顶部涡旋区的尺度；
根据气象数据和建筑物数据构建初始风场；
根据质量守恒方程调整所述初始风场，得到最终风场。
2.根据权利要求1所述的建筑物可分辨尺度的城市三维风场快速分析方法，其特征在于，所述质量守恒方程中包括斜坡前方、后方、右侧及左侧边界条件系数，所述斜坡包括建筑物倾斜屋顶。
3.根据权利要求1所述的建筑物可分辨尺度的城市三维风场快速分析方法，其特征在于，在所述质量守恒方程中对于底面及建筑物采用第二类边界条件
Figure FDA0003361053600000011
其他边界采用第一类自由边界条件λ＝0。
4.根据权利要求1所述的建筑物可分辨尺度的城市三维风场快速分析方法，其特征在于，所述根据质量守恒方程调整所述初始风场，得到最终风场，包括：
采用变分方法求解所述质量守恒方程，且α1＝α2＝0.5。
5.根据权利要求1所述的建筑物可分辨尺度的城市三维风场快速分析方法，其特征在于，所述建筑物数据包括屋顶是否为坡度、建筑物高度、建筑物长度及建筑物最长边和正东方向夹角角度。
6.根据权利要求1所述的建筑物可分辨尺度的城市三维风场快速分析方法，其特征在于，所述确认各个建筑物在模型网格内的尺寸，包括：
当入流风和建筑物的夹角θ不超过20°时，建筑物在模型网格内的尺寸Weff、Leff与其实际尺寸W、L相等；
当入流风和建筑物的夹角大于20°时，Weff＝Wsinθ+Lcosθ，Leff＝Wcosθ+Lsinθ；
其中，W为宽度，L为长度。
7.根据权利要求1所述的建筑物可分辨尺度的城市三维风场快速分析方法，其特征在于，所述根据气象数据和建筑物数据构建初始风场，包括：
计算空旷区域和建筑物区域的入流风的垂直廓线；
对单体建筑物周围风场参数化；
对建筑街渠风场参数化。
8.根据权利要求7所述的建筑物可分辨尺度的城市三维风场快速分析方法，其特征在于，对单体建筑物周围风场参数化，包括：
当入流风和建筑物的夹角θ大于20°时，采用迎风面的入流风速分解为平行迎风面分量和垂直迎风面分量进行插值。
9.根据权利要求7所述的建筑物可分辨尺度的城市三维风场快速分析方法，其特征在于，对建筑街渠风场参数化，包括：
当S/H>2.5或1.4<S/H<2.4时，采用单体建筑物周围风场参数化方案；
当S/H<1.4时，采用SR参数化方案或MR参数化方案；
其中，H为建筑物高度，S为建筑间距离。
CN202111367080.0A 2021-11-18 2021-11-18 一种建筑物可分辨尺度的城市三维风场快速分析方法 Pending CN114386215A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111367080.0A CN114386215A (zh) 2021-11-18 2021-11-18 一种建筑物可分辨尺度的城市三维风场快速分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111367080.0A CN114386215A (zh) 2021-11-18 2021-11-18 一种建筑物可分辨尺度的城市三维风场快速分析方法

Publications (1)

Publication Number Publication Date
CN114386215A true CN114386215A (zh) 2022-04-22

Family

ID=81196342

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111367080.0A Pending CN114386215A (zh) 2021-11-18 2021-11-18 一种建筑物可分辨尺度的城市三维风场快速分析方法

Country Status (1)

Country Link
CN (1) CN114386215A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115146564A (zh) * 2022-08-02 2022-10-04 上海市气候中心(上海区域气候中心) 基于垂直分层降尺度技术的城市地面风速精细化模拟方法
CN116911214A (zh) * 2023-08-31 2023-10-20 中国科学院大气物理研究所 一种街区尺度气溶胶扩散模拟方法及装置

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115146564A (zh) * 2022-08-02 2022-10-04 上海市气候中心(上海区域气候中心) 基于垂直分层降尺度技术的城市地面风速精细化模拟方法
CN116911214A (zh) * 2023-08-31 2023-10-20 中国科学院大气物理研究所 一种街区尺度气溶胶扩散模拟方法及装置
CN116911214B (zh) * 2023-08-31 2023-12-01 中国科学院大气物理研究所 一种街区尺度气溶胶扩散模拟方法及装置

Similar Documents

Publication Publication Date Title
Palusci et al. Impact of morphological parameters on urban ventilation in compact cities: The case of the Tuscolano-Don Bosco district in Rome
Ricci et al. Local-scale forcing effects on wind flows in an urban environment: Impact of geometrical simplifications
Antoniou et al. CFD and wind-tunnel analysis of outdoor ventilation in a real compact heterogeneous urban area: Evaluation using “air delay”
Panagiotou et al. City breathability as quantified by the exchange velocity and its spatial variation in real inhomogeneous urban geometries: An example from central London urban area
Ng et al. Improving the wind environment in high-density cities by understanding urban morphology and surface roughness: A study in Hong Kong
Blocken et al. CFD simulation for pedestrian wind comfort and wind safety in urban areas: General decision framework and case study for the Eindhoven University campus
Wang et al. Air ventilation assessment under unstable atmospheric stratification—A comparative study for Hong Kong
Piroozmand et al. Coupled CFD framework with mesoscale urban climate model: Application to microscale urban flows with weak synoptic forcing
Shao et al. PIGNN-CFD: A physics-informed graph neural network for rapid predicting urban wind field defined on unstructured mesh
CN114386215A (zh) 一种建筑物可分辨尺度的城市三维风场快速分析方法
CN111767597B (zh) 一种城市模型验证方法、装置、设备及存储介质
Yuan et al. Practical application of CFD on environmentally sensitive architectural design at high density cities: A case study in Hong Kong
He et al. Hybrid framework for rapid evaluation of wind environment around buildings through parametric design, CFD simulation, image processing and machine learning
Ricci et al. Impact of surrounding environments and vegetation on wind comfort assessment of a new tower with vertical green park
Qin et al. Influence of site and tower types on urban natural ventilation performance in high-rise high-density urban environment
Li et al. Wind tunnel study on the morphological parameterization of building non-uniformity
CN114357571A (zh) 建成建筑环境下大气边界层风场特性的反演方法和系统
Węgrzyński et al. Wind and fire coupled modelling—Part II: Good practice guidelines
Li et al. Exploring urban space quantitative indicators associated with outdoor ventilation potential
Hågbo et al. Pedestrian wind comfort assessment using computational fluid dynamics simulations with varying number of wind directions
Xu et al. A review of simplified numerical approaches for fast urban airflow simulation
Wang et al. Evaluating a combined WRF and CityFFD method for calculating urban wind distributions
He et al. A novel methodology for architectural wind environment study by integrating CFD simulation, multiple parametric tools and evaluation criteria
Ding et al. Coupled natural ventilation modeling for contextual parametric design decision support
Lim et al. Multicollinearity issue for the parameterization of urban ventilation potential with urban morphology

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