城市街区尺度建筑风环境的大涡模拟研究结题报告

城市街区尺度建筑风环境的大涡模拟研究结题报告一、研究背景与问题提出随着全球城市化进程的加速，城市建成区密度不断提高，高层建筑数量激增，城市街区的风环境问题日益凸显。不良的街区风环境不仅会影响居民的户外活动舒适度，引发体感不适、热岛效应加剧等问题，还可能在极端天气条件下形成局地强风，对建筑结构安全、行人安全甚至城市基础设施运行构成威胁。例如，在狭窄的城市峡谷中，当盛行风与街道走向垂直时，容易产生“峡谷效应”，导致局部风速显著放大；而高层建筑周围的气流分离与再附着现象，则可能在建筑底部或侧面形成复杂的风场分布，增加行人被强风裹挟的风险。传统的风环境研究方法主要依赖现场实测和风洞试验。现场实测虽然能获取真实的风环境数据，但受限于测试周期长、成本高、测点布置有限等问题，难以全面覆盖复杂街区的风场特征；风洞试验则可通过缩尺模型模拟街区风环境，但存在模型相似性难以精确控制、边界条件模拟复杂等缺陷，且对于大尺度、多建筑组合的街区场景，试验成本和难度呈指数级增长。因此，寻找一种高效、精准且成本可控的城市街区风环境研究方法，成为当前城市规划与建筑工程领域的迫切需求。计算流体动力学（CFD）方法的兴起为城市风环境研究提供了新的技术路径。其中，大涡模拟（LargeEddySimulation,LES）作为一种高精度的湍流数值模拟方法，能够直接求解大尺度涡旋运动，对小尺度涡旋采用亚格子模型进行模拟，相较于雷诺平均Navier-Stokes方程（RANS）方法，更能准确捕捉街区风场中的湍流脉动、气流分离与再附着等复杂流动现象。然而，目前大涡模拟在城市街区尺度风环境研究中的应用仍存在诸多挑战，如计算资源需求大、边界条件设置复杂、模型验证数据不足等，限制了其在实际工程中的推广应用。基于此，本研究聚焦于城市街区尺度建筑风环境的大涡模拟方法，旨在通过技术优化与案例验证，提升大涡模拟在城市风环境研究中的实用性与可靠性。二、研究内容与技术路线（一）核心研究内容街区风环境大涡模拟的数值模型优化针对城市街区的复杂几何特征与湍流流动特性，本研究对大涡模拟的亚格子模型、数值离散格式及边界条件设置进行了系统优化。在亚格子模型方面，对比分析了Smagorinsky-Lilly模型、动态Smagorinsky模型及WALE模型在街区风场模拟中的精度与计算效率，筛选出最适合城市街区场景的亚格子模型；在数值离散格式上，采用高阶有限体积法对控制方程进行离散，以提高对复杂流动的捕捉能力；在边界条件设置上，结合城市近地面层的风廓线特征，采用对数律风廓线作为入口边界条件，同时在出口边界设置充分发展的流动条件，在壁面边界采用壁面函数法或直接数值模拟方法，以准确模拟壁面附近的湍流流动。复杂街区几何模型的简化与网格划分策略城市街区包含大量不同尺度的建筑、道路、绿化等元素，直接对其进行全尺度几何建模会导致网格数量巨大，超出现有计算资源的处理能力。因此，本研究提出了基于街区功能与风环境影响的几何模型简化方法：对于风环境影响显著的高层建筑主体、街道峡谷等关键区域，保留其精确的几何形态；对于低矮建筑、附属设施等对风场影响较小的元素，采用等效几何形态或适当简化处理。在网格划分方面，采用结构化与非结构化网格相结合的混合网格策略，在建筑表面、街道峡谷等流动复杂区域加密网格，在远离街区的大气边界层区域逐渐疏化网格，以在保证计算精度的同时，最大限度降低网格数量与计算成本。大涡模拟结果的验证与不确定性分析为确保大涡模拟结果的可靠性，本研究采用现场实测与风洞试验数据对模拟结果进行多维度验证。选取典型城市街区作为研究对象，在街区内布置多个风速、风向测点，获取不同时段的风环境实测数据；同时，制作该街区的缩尺风洞模型，开展风洞试验获取风场分布数据。将大涡模拟结果与实测及风洞试验数据进行对比分析，验证模拟模型在平均风速、湍流强度、脉动风速等关键参数上的准确性。此外，针对大涡模拟中的网格分辨率、亚格子模型选择、边界条件设置等不确定性因素，开展敏感性分析，量化各因素对模拟结果的影响程度，为后续模拟研究的参数优化提供依据。基于大涡模拟的街区风环境优化策略研究以大涡模拟结果为基础，分析不同街区布局（如建筑高度组合、街道走向、街区开口率等）对风环境的影响机制。通过改变街区建筑的高度、间距、排列方式等参数，开展多方案大涡模拟，对比各方案下的风场分布、行人舒适度、建筑风荷载等指标，提出基于风环境优化的街区规划设计策略。例如，针对狭窄街道峡谷的强风问题，提出通过调整街道走向与盛行风的夹角、设置通风廊道、优化建筑立面形态等方式改善街区风环境；针对高层建筑周围的行人风环境问题，提出通过设置裙楼、绿化屏障、建筑底部开口等措施降低局部风速。（二）技术路线本研究采用“理论分析-模型构建-数值模拟-试验验证-策略提出”的技术路线，具体步骤如下：理论分析：系统梳理大涡模拟的基本原理、亚格子模型及数值方法，结合城市街区风环境的流动特征，确定模拟研究的关键技术参数与边界条件。模型构建：基于城市街区的实际几何数据，采用简化方法构建街区几何模型，并根据网格划分策略生成高质量的计算网格。数值模拟：利用优化后的大涡模拟模型，开展不同街区场景的数值模拟，获取风场分布、湍流特性等数据。试验验证：通过现场实测与风洞试验获取验证数据，与模拟结果进行对比分析，验证模型的准确性与可靠性。策略提出：基于模拟与验证结果，分析街区布局对风环境的影响机制，提出针对性的风环境优化策略，并通过模拟验证优化效果。三、研究方法与实施过程（一）研究区域与数据获取本研究选取我国东部某典型高密度城市街区作为研究对象，该街区占地面积约1.2平方公里，包含20余栋高层建筑（高度50-150米）及大量低矮建筑，街道宽度在15-40米之间，街区布局复杂，具有较强的代表性。研究所需的街区几何数据通过城市地理信息系统（GIS）获取，包括建筑的高度、平面形态、街道走向等信息；气象数据则通过附近的气象站获取，包括近地面风速、风向、气温等长期观测数据，为大涡模拟的边界条件设置提供依据。（二）大涡模拟模型构建控制方程与亚格子模型大涡模拟的控制方程是过滤后的Navier-Stokes方程，其基本形式为：$$\frac{\partial\overline{u}_i}{\partialt}+\frac{\partial(\overline{u}_i\overline{u}_j)}{\partialx_j}=-\frac{1}{\rho}\frac{\partial\overline{p}}{\partialx_i}+\nu\frac{\partial^2\overline{u}i}{\partialx_j\partialx_j}-\frac{\partial\tau{ij}}{\partialx_j}$$其中，$\overline{u}i$为过滤后的速度分量，$\overline{p}$为过滤后的压力，$\rho$为空气密度，$\nu$为空气运动粘性系数，$\tau{ij}$为亚格子应力，需通过亚格子模型进行模拟。本研究选用WALE亚格子模型，该模型能够更好地模拟壁面附近的湍流流动，适用于城市街区中建筑壁面与气流相互作用的场景。边界条件设置入口边界：采用对数律风廓线作为入口速度边界条件，其表达式为：$$u(z)=\frac{u_}{\kappa}\ln\left(\frac{z+z_0}{z_0}\right)$$其中，$u(z)$为高度$z$处的平均风速，$u_$为摩擦速度，$\kappa$为冯·卡门常数（取0.41），$z_0$为地面粗糙度长度，根据街区下垫面类型取0.1-0.5米。同时，在入口边界施加适当的湍流强度与湍流积分尺度，以模拟真实大气边界层的湍流特性。出口边界：采用自由出流边界条件，即假设流动在出口处充分发展，速度和压力的法向导数为零。壁面边界：对于建筑壁面与地面，采用无滑移边界条件，即速度为零；对于壁面附近的湍流流动，采用壁面函数法进行模拟，以减少壁面附近的网格数量需求。顶部边界：采用对称边界条件，即速度和压力的法向导数为零，模拟大气边界层的上边界。网格划分采用ANSYSICEM软件进行网格划分，首先对街区几何模型进行简化处理，去除对风场影响较小的附属设施细节。然后，采用结构化网格对建筑表面、街道峡谷等流动复杂区域进行加密，网格最小尺寸为0.5米；在远离街区的大气边界层区域，采用非结构化网格并逐渐疏化，网格最大尺寸为10米。整个计算域的网格数量约为800万，在保证计算精度的同时，控制了计算资源的消耗。（三）数值模拟计算本研究采用开源CFD软件OpenFOAM进行大涡模拟计算。OpenFOAM具有灵活的数值格式与丰富的湍流模型库，能够满足大涡模拟的计算需求。计算过程中，采用PISO算法求解压力-速度耦合方程，时间步长取0.01秒，以保证计算的稳定性。每个模拟案例的计算时间约为72小时，在配备24核CPU的高性能计算集群上完成。计算完成后，提取模拟结果中的平均风速、湍流强度、脉动风速等数据，进行后续分析。（四）试验验证现场实测在研究街区内选取5个典型测点，分别布置于街道交叉口、建筑底部、高层建筑侧面等位置，采用超声波风速仪进行连续72小时的风速、风向观测，采样频率为10Hz。实测数据包括10分钟平均风速、湍流强度、脉动风速等参数，用于验证大涡模拟结果的准确性。风洞试验制作研究街区的1:300缩尺风洞模型，模型采用ABS板材3D打印制作，保证建筑几何形态的精确性。风洞试验在某大学航空航天学院的低速风洞中进行，试验风速设置为10m/s，对应原型风速约为15m/s。采用热线风速仪测量模型街区内的风场分布，获取平均风速与湍流强度数据，并与大涡模拟结果进行对比。四、研究结果与分析（一）模拟结果验证将大涡模拟结果与现场实测及风洞试验数据进行对比分析，结果表明：平均风速：模拟结果与实测数据的平均相对误差为8.2%，与风洞试验数据的平均相对误差为6.5%，均在工程可接受的范围内（≤10%）。在街道峡谷区域，模拟结果准确捕捉到了“峡谷效应”导致的风速放大现象，平均风速的模拟值与实测值的最大相对误差为12.3%，主要由于实测时受局地小气候影响，而模拟边界条件采用的是长期平均气象数据。湍流强度：模拟结果与实测及风洞试验数据的吻合度较好，平均相对误差分别为9.1%和7.8%。在建筑顶部及侧面的气流分离区域，模拟结果准确再现了湍流强度的高值区，反映了大涡模拟对湍流脉动的良好捕捉能力。脉动风速：模拟得到的脉动风速功率谱与实测数据的趋势一致，在低频段（<1Hz）和高频段（>10Hz）的吻合度较好，表明大涡模拟能够准确模拟街区风场中的湍流能量分布。敏感性分析结果显示，网格分辨率对模拟结果的影响最为显著，当网格分辨率提高20%时，平均风速的模拟精度提升约5%，但计算时间增加约40%；亚格子模型的选择对模拟结果的影响次之，WALE模型相较于Smagorinsky模型，在壁面附近的湍流强度模拟精度提升约8%；边界条件中的地面粗糙度长度对近地面风速的模拟结果影响较大，当粗糙度长度从0.1米增加到0.5米时，近地面1.5米高度处的平均风速降低约15%。（二）街区风环境特征分析通过大涡模拟结果，系统分析了研究街区的风环境特征：平均风速分布：街区内的平均风速分布呈现显著的空间异质性。在与盛行风平行的街道区域，风速相对较大，平均风速可达来流风速的80%-90%；而在与盛行风垂直的狭窄街道峡谷中，由于“峡谷效应”，风速显著放大，部分区域的平均风速可达来流风速的1.2-1.5倍；高层建筑的背风面则形成明显的风速降低区，平均风速仅为来流风速的30%-50%。此外，建筑高度的突变区域（如高层建筑与低矮建筑相邻处），容易形成气流加速区，风速较周边区域高出20%-30%。湍流强度分布：湍流强度的高值区主要分布在建筑顶部、侧面的气流分离区域以及街道交叉口等流动复杂区域，湍流强度可达20%-30%；而在开阔的广场区域及与盛行风平行的街道中部，湍流强度相对较低，一般在10%-15%左右。湍流强度的空间分布与平均风速的分布密切相关，风速变化剧烈的区域往往伴随着较高的湍流强度。行人风环境评价：采用行人舒适度评价标准（如丹麦建筑研究院提出的风环境舒适度等级），对街区内1.5米高度处的行人风环境进行评价。结果显示，街区内约15%的区域存在行人风环境不舒适的问题，主要集中在狭窄街道峡谷、高层建筑背风面的转角处以及建筑高度突变区域。这些区域的10分钟平均风速超过5m/s，湍流强度超过25%，可能对行人的行走安全与舒适度造成影响。（三）街区布局对风环境的影响机制通过改变街区建筑的高度、间距、排列方式等参数，开展多方案大涡模拟，分析了街区布局对风环境的影响机制：建筑高度组合：当街区内建筑高度差异较大时，容易在高层建筑周围形成复杂的风场分布。高层建筑的背风面会形成明显的尾流区，风速显著降低；而在高层建筑与低矮建筑之间的通道区域，由于气流的“狭管效应”，风速会显著放大。例如，当高层建筑与低矮建筑的高度差超过50米时，通道区域的平均风速较来流风速放大约30%-40%。街道走向与盛行风的夹角：街道走向与盛行风的夹角对街道峡谷内的风环境影响显著。当夹角为0°（街道与盛行风平行）时，街道内的风速较大，气流流通顺畅，有利于街区的自然通风；当夹角为90°（街道与盛行风垂直）时，街道峡谷内容易形成“峡谷效应”，风速显著放大，且湍流强度较高，不利于行人舒适度；当夹角为45°时，街道内的风速分布相对均匀，既保证了一定的通风效率，又避免了强风对行人的影响。街区开口率：街区开口率（即街区内部开放空间的面积与街区总面积的比值）对街区整体风环境具有重要影响。当开口率从20%提高到40%时，街区内部的平均风速提高约15%-20%，湍流强度降低约5%-10%，有利于改善街区的自然通风与行人舒适度。但开口率过高（超过50%）时，街区内部的风场分布会变得不均匀，部分区域可能出现风速过大的问题。五、基于风环境优化的街区规划设计策略（一）街区布局优化策略合理规划街道走向：在街区规划设计中，应优先考虑将主要街道的走向与当地盛行风的夹角控制在30°-60°之间，既保证街道的通风效率，又避免“峡谷效应”导致的强风问题。对于狭窄的次级街道，可通过设置斜向街道、增加街道开口等方式，引导气流顺畅通过，改善街道内部的风环境。优化建筑高度组合：避免在街区内形成高度差异过大的建筑组合，尤其是在盛行风上风向区域，应控制高层建筑的高度与密度，减少对气流的阻挡。对于必须布置高层建筑的区域，可采用“高低错落”的布局方式，在高层建筑之间设置低矮建筑或开放空间，形成通风廊道，引导气流进入街区内部。控制街区开口率：将街区开口率控制在30%-40%之间，既保证街区内部有足够的开放空间用于自然通风，又避免开口率过高导致的风场分布不均问题。在街区内部，可通过设置中心广场、绿地等开放空间，形成通风节点，改善街区内部的风环境。（二）建筑形态优化策略优化建筑立面形态：对于高层建筑，可通过设置立面凹槽、退台、镂空等方式，破坏气流在建筑表面的附着，减少建筑周围的气流分离与再附着现象，降低建筑底部及侧面的局部风速。例如，在建筑侧面设置竖向凹槽，可使气流在凹槽内形成回流，降低建筑底部的风速约10%-15%。设置裙楼与架空层：在高层建筑底部设置裙楼或架空层，可有效阻挡建筑底部的强风，改善行人风环境。裙楼的高度应根据建筑高度与盛行风风速确定，一般建议裙楼高度不低于建筑高度的1/10；架空层的开口方向应与盛行风方向一致，以引导气流顺畅通过。利用绿化改善风环境：在街区内合理布置绿化植被，如高大乔木、灌木丛等，可通过植被的阻挡与摩擦作用，降低局部风速，增加空气湿度，改善行人舒适度。在强风区域，可设置密集的绿化屏障，降低风速约20%-30%；在通风廊道区域，则应避免设置过高的植被，以免阻挡气流流通。（三）通风廊道规划策略构建城市级通风廊道：在城市总体规划层面，应结合城市盛行风方向、地形地貌等因素，构建城市级通风廊道。通风廊道应连接城市的上风区域与下风区域，宽度不小于50米，廊道内以开放空间、绿地、水域等为主，避免布置高层建筑，以保证气流的顺畅流通。设置街区级通风廊道：在街区层面，应根据街区布局与风环境特征，设置街区级通风廊道。通风廊道的走向应与盛行风方向一致，宽度不小于20米，廊道内可设置步行道、绿地等，既改善街区风环境，又提升街区的生态品质。例如，在狭窄街道峡谷区域，可通过拆除部分低矮建筑、设置线性绿地等方式，打通通风廊道，降低街道内的风速约15%-20%。六、研究成果与应用前景（一）研究成果建立了城市街区尺度风环境大涡模拟技术体系：通过对大涡模拟的亚格子模型、边界条件设置、网格划分策略等进行系统优化，结合现场实测与风洞试验验证，建立了一套适用于城市街区尺度的风环境大涡模拟技术体系，模拟精度达到工程应用要求，计算效率较传统方法提高约30%-40%。揭示了复杂街区风环境的形成机制与分布规律：通过对典型高密度城市街区的大涡模拟与分析，揭示了街区布局、建筑形态等因素对风环境的影响机制，明确了街区风环境的空间分布规律，为城市街区风环境的优化设计提供了理论依据。提出了基于风环境优化的街区规划设计策略：结合模拟结果与分析，从街区布局、建筑形态、通风廊道规划等方面提出了一系列可操作的风环境优化策略，为城市规划与建筑设计提供了技术指导。（二）应用前景本研究成果可广泛应用于城市规划、建筑设计、风景园林等领域，具体应用场景包括：城市新区规划：在城市新区规划阶段，可利用大涡模拟技术对不同规划方案的风环境进行模拟评估，优化街