[发明专利]一种街谷水平通风影响城市热岛效应的数值模拟方法

说明书

本发明公开的街谷水平通风影响城市热岛效应的数值模拟方法，涉及模拟计算技术领域，通过采用高斯‑赛戴尔Gauss–Seidel迭代算法求解各条街谷交点处的空气流量初值及空气压力初值，采用牛顿迭代Newton‑Raphson算法求解各条街谷交点处的空气流量值及空气压力值，随着街谷数量的增加，采用Gauss–Seidel迭代算法与Newton‑Raphson算法的组合比仅采用Newton‑Raphson算法迭代次数更少，提高了收敛速度及精度，既避免了采用大规模模拟计算城市热岛效应的负担，又可在街区内存在多条街谷的规模下，保证了CFD数学模型的计算精度。

技术领域

本发明涉及模拟计算技术领域，具体涉及一种街谷水平通风影响城市热岛效应的数值模拟方法。

背景技术

随着城市不断高速发展，城市热岛效应日益严重。缓解城市热岛，有助于抑制传染病传播、减少温室气体排放、降低建筑能耗。

目前对城市热岛效应的研究过程中，在给城市街谷建模时一般采用经典的城市冠层模型。而经典城市冠层模型的水平通风特性基本采用“对热岛效应的影响很小而忽略”的处理方式，这对许多处于多风环境下城市的热岛效应进行研究显然是不足的。比如许多沿海城市的大气运动会受到热岛环流和常年的海陆风环流的相互作用，此时城市的水平通风特性对其热岛效应的影响很大。

近年来，由于计算机流体力学CFD技术的大幅提高，CFD数学模型也逐渐被广泛采用。但是采用CFD数学模型模拟整个城市热岛现象在计算技术和理论建模上都将极具挑战性。在建模方面，CFD技术对计算机资源消耗巨大，而且需要详细的地貌资料，而提供城市冠层内各建筑外部结构的三维几何尺寸数据库是代价相当昂贵的一件事。因此，在采用CFD数学模型模拟时，也通常需要简化建筑的几何形状，如此一来，空间分辨率就有所下降。计算技术上，CFD数学模型将面对大量的控制体或节点来进行计算，边界条件的测量较为困难。另外，准确地对湍流的处理仍然是CFD数学模型的一个重要弱点，湍流建模的准确性会对CFD数学模型模拟空气流场的精度有很大影响。目前主流的两种对城市热岛效应的数值模拟方法都有一定缺陷。

发明内容

为解决现有技术的不足，本发明实施例提供了一种街谷水平通风影响城市热岛效应的数值模拟方法，该方法包括以下步骤：

S1，利用高斯-赛戴尔Gauss-Seidel迭代算法，计算各条街谷交点处的空气流量初值和空气压力初值；

S2，利用牛顿迭代Newton-Raphson算法，计算各条街谷交点处的空气流量值和空气压力值；

S3，根据公式计算处于上游各条街谷出口处的混合空气内能值e_mix，其中，V_in、e_in分别为流入各条街谷的空气流量及空气内能值，ρ为空气密度，V_out为流出各条街谷的空气流量；

S4，根据空气物性方程获得各条街谷交点处的温度值。

优选地，其特征在于，利用高斯-赛戴尔Gauss-Seidel迭代算法，计算各条街谷交点处的空气流量初值和空气压力初值包括：

S11分别对各条街谷及各条街谷的交点进行编号；

S12根据各条街谷的编号及各条街谷的交点的编号，得到相应的邻接矩阵A，其中，邻接矩阵A定义为A＝[a_ji]_(j+1)×i，

S13假设根据公式计算对角阵r₀，其中，ρ为空气密度，f为街谷横截面的水力粗糙度，L为街谷长度，为街谷湿周直径，S为摩阻系数，为街谷内初始平均空气流量；

S14根据公式计算各条街谷交点处的空气压力P₁，V_c是街谷建筑外泄空气流量或边界交点处的空气流量；