一种基于OpenFOAM的模拟复杂地形地表粗糙度的方法技术

本发明专利技术属于计算风工程领域，并公开了一种基于OpenFOAM的模拟复杂地形地表粗糙度的方法，包含以下步骤：1)建立复杂地形CFD模型，生成精细化网格；2)提取复杂地形CFD模型的地表粗糙度长度数据，将地表粗糙度长度转化为地表覆盖植被高度，为模拟地表覆盖植被对复杂地形风场影响提供准确的数据基础；3)计算地表覆盖植被在不同高度处对风场的阻力系数，然后将地表覆盖植被高度转换为相应的阻力系数；4)使用OpenFOAM添加阻力系数和地表覆盖植被阻力项，从而获得更精确的复杂地形风场特性模拟结果。本发明专利技术使用计算机支持的程序语言，考虑了地表粗糙度对复杂地形风场的影响，有利于提升复杂地形风场数值模拟的精确度。

A method for simulating surface roughness of complex terrain based on OpenFOAM

The invention belongs to the field of wind engineering, and discloses a method for simulating the surface roughness of complex terrain based on OpenFOAM. The method comprises the following steps: 1) establishing CFD model of complex terrain to generate refined meshes; 2) extracting the surface roughness length data of CFD model of complex terrain to convert the surface roughness length into surface roughness length. The height of vegetation cover provides accurate data basis for simulating the influence of vegetation cover on wind field of complex terrain; 3) Calculate the resistance coefficient of vegetation cover to wind field at different heights, and then convert the height of vegetation cover to the corresponding resistance coefficient; 4) Use OpenFOAM to add resistance coefficient and vegetation cover the surface. A more accurate simulation result of wind field characteristics in complex terrain is obtained. The invention uses a computer-supported programming language and considers the influence of surface roughness on wind field of complex terrain, which is beneficial to improving the accuracy of wind field numerical simulation of complex terrain.

【技术实现步骤摘要】
一种基于OpenFOAM的模拟复杂地形地表粗糙度的方法
本专利技术属于计算风工程领域，更具体地，涉及一种模拟复杂地形地表粗糙度的方法。
技术介绍
地表粗糙度长度，是指因地表起伏不平或地物本身几何形状的影响，风速廓线上风速为零的位置并不在地表(高度为零处)，而在离地表的一定高度。这一高度也称为空气动力学粗糙度长度，是地表对经过的流体的型态和阻滞力产生影响的一个综合特征参数。地表粗糙度长度是由于地表的植被等产生的，表示了地表的空气动力学特征，反应地表对风速的减弱作用，取决于地面粗糙因素，如粗糙元素的几何形状、大小高度、分布类型和分布密度。研究地表粗糙度长度对分析近地层中气候与环境的变化、风场模拟和风资源评估的准确性具有重要意义。对于地表粗糙度长度的模拟，一种方法是基于湍流模型的壁面函数法。Blocken等从理论上讨论了传统CFD(ComputationalFluidDynamics,计算流体动力学)设置在模拟粗糙壁面条件风场方面的制约性，并提出了相应的改进方法；Knopp等提出了CFD模拟中适用于模型的粗糙度修正方法。但这些研究没有针对实际地形测试风场分布对地表粗糙度的敏感性，也没有对实际复杂地形中的地表粗糙度长度变化做出适当调整。邓院昌等对整个模拟区域的下垫面进行粗糙度精细化设置来提高模拟精度，但是在数值模拟进行区域分块时均设置了相同的粗糙度，并未最大程度地考虑复杂地形的实际情况。而且使用基于湍流模型的壁面函数法模拟复杂地表粗糙度长度时，无法捕捉粗糙区域内的湍流结构，也没有考虑地表覆盖植被相关参数的影响。另一种方法是通过在Navier-Stokes动量方程中加入阻力项来模拟，根据Enoki与Ishihara的研究采用如下表达式计算阻力项：式中，为阻力项，ρ为密度，Cd为阻力系数，为来流风速，为来流风速分量。采用CFD方法能够模拟大气边界层中的湍流在复杂地形中产生的撞击、分离、环绕等现象；评估复杂地形中的风能分布时具有具有费用低、周期短、便于模拟真实环境、描述风场细节和给出风场定量结果的优点。CFD商业软件简单易用，但是其购买费用昂贵，并且源程序不公开，所提供的二次开发接口有较大局限性，难以满足用户的全部需求。
技术实现思路
针对现有技术的以上缺陷或改进需求，本专利技术提供了一种基于OpenFOAM的模拟复杂地形地表粗糙度的方法，该方法使用计算机支持的程序语言，能够基于GIS数据建立复杂地形CFD模型，提取地表粗糙度长度数据，使用OpenFOAM添加阻力系数和地表覆盖植被阻力项，考虑了地表粗糙度对复杂地形风场影响，有利于提升复杂地形风场数值模拟的精确度。为实现上述目的，按照本专利技术，提供了一种基于OpenFOAM的模拟复杂地形地表粗糙度的方法，其特征在于，该方法包含以下步骤：1)在OpenFOAM中基于GIS数据建立复杂地形CFD模型，生成精细化网格；2)基于GIS数据提取复杂地形CFD模型的地表粗糙度长度数据，然后将一系列分布在等值线上的随机点的地表粗糙度长度转化为地表覆盖植被高度；3)计算地表覆盖植被在不同高度处对风场的阻力系数，然后将地表覆盖植被高度转换为相应的阻力系数；4)使用OpenFOAM添加阻力系数和地表覆盖植被阻力项，从而获得复杂地形风场特性模拟结果。优选地，步骤1)具体包括以下子步骤：1.1)确定目标的中心经纬度，从GIS地理信息系统中获取目标所在区域的ASTERGDEM30米分辨率大地坐标下的地形数据，使用ArcMap处理GIS数据，从而获得WGS_1984_UTM坐标系下的地形数据网格节点；1.2)参数化建立三维平坦地形CFD模型并生成精细化网格：复杂地形CFD模型的计算域设置为长方体区域，根据实际的地形地貌找到对计算域有影响的复杂地形，计算域需要包含这些复杂地形，而且采用非结构化三棱柱网格以充分拟合复杂几何边界，其中，在计算域的不同区域采用分辨率不同的网格，复杂地形CFD模型在竖直方向上采用相邻网格尺寸比值为定值的σ网格；1.3)实现计算域的地形数据网格节点与三维平坦地形CFD模型网格节点之间的数据转换：根据三维平坦地形CFD模型中每一个底面网格数据节点的坐标，找到其在地形数据网格节点中的对应位置，然后使用该位置周围的3个地形数据网格节点的高程，通过三角形共面线性插值得到底面网格节点的高程，同时，在三维平坦地形CFD模型中，竖直方向上其余网格节点的高程会随着底面网格节点高程的变化产生相应变化，最终得到修改高程后的复杂地形CFD模型；1.4)通过转化器gambitToFoam导入OpenFOAM，将复杂地形CFD模型转换为OpenFOAM可识别的格式。优选地，步骤2)具体包括以下子步骤：2.1)从GIS地理信息系统中提取复杂地形CFD模型的地表粗糙度长度数据，地表粗糙度长度数据的分布形式为一系列地表粗糙度长度等值线上的随机点；2.2)处理地表粗糙度长度数据：判断各随机点是否在复杂地形CFD模型计算域内，如果在计算域内，将地表粗糙度长度转换为地表覆盖植被高度，否则弃之不用；其中，地表粗糙度长度与地表覆盖植被高度的转换表达式如下：h＝a·z0式中，h为地表覆盖植被高度，a为换算系数，z0为地表粗糙度长度；2.3)建立矩形分布的网格节点，使用最近邻点插值法进行赋值，得到矩形分布网格节点的地表覆盖植被高度，其中，最近邻点插值法在搜索时需要设置一个搜索半径，搜索半径不小于矩形分布网格的尺寸，如果搜索半径内没有邻近点，则该节点的地表覆盖植被高度为0。优选地，步骤3)具体包括以下子步骤：3.1)根据复杂地形CFD模型底层的每一个三棱柱网格单元的中心坐标(x，y)，找到其在矩形分布的网格中的对应位置，然后使用该中心坐标(x，y)周围的3个地表覆盖植被高度网格节点，通过三角形共面线性插值得到该中心坐标(x，y)处的地表覆盖植被高度；3.2)从该中心坐标(x，y)处的底层网格单元开始，依次向上，判断每一个三棱柱网格单元距地面的竖直高度z和中心坐标(x，y)处的地表覆盖植被高度之间的关系，如果三棱柱网格单元距地面的竖直高度z不大于中心坐标(x，y)处的地表覆盖植被高度，计算三棱柱网格单元对应的阻力系数，否则阻力系数为0，阻力系数的表达式如下：Cd＝Ct·at式中，Cd为阻力系数，Ct为抗力系数，at为叶面积密度。优选地，步骤4)具体包括以下子步骤：4.1)对于复杂地形CFD模型的边界条件，底面采用壁面，入口采用速度入口，出口采用压力出口，其余侧面和顶面采用对称边界条件；4.2)使用OpenFOAM添加每一个三棱柱网格单元对应的阻力系数；4.3)在Navier-Stokes动量方程中加入阻力项，考虑了地表粗糙度对风场的影响，使用OpenFOAM计算，从而获得复杂地形风场特性模拟结果，阻力项的表达式如下：式中，为阻力项，ρ为密度，Cd为阻力系数，为来流风速，为来流风速的速度分量。总体而言，通过本专利技术所构思的以上技术方案与现有技术相比，能够取得下列有益效果：(1)本专利技术可以实现在OpenFOAM中基于GIS数据建立复杂地形CFD模型，提取地表粗糙度长度数据并转化为地表覆盖植被高度，计算地表覆盖植被在不同高度处对风场的阻力系数，然后将地表覆盖植被高度转换为相应的阻力系数，使用OpenFOAM添加阻力系数和地表覆盖本文档来自技高网...

【技术保护点】
1.一种基于OpenFOAM的模拟复杂地形地表粗糙度的方法，其特征在于，该方法包含以下步骤：1)在OpenFOAM中基于GIS数据建立复杂地形CFD模型，生成精细化网格；2)基于GIS数据提取复杂地形CFD模型的地表粗糙度长度数据，然后将一系列分布在等值线上的随机点的地表粗糙度长度转化为地表覆盖植被高度；3)计算地表覆盖植被在不同高度处对风场的阻力系数，然后将地表覆盖植被高度转换为相应的阻力系数；4)使用OpenFOAM添加阻力系数和地表覆盖植被阻力项，从而获得复杂地形风场特性模拟结果。

【技术特征摘要】
1.一种基于OpenFOAM的模拟复杂地形地表粗糙度的方法，其特征在于，该方法包含以下步骤：1)在OpenFOAM中基于GIS数据建立复杂地形CFD模型，生成精细化网格；2)基于GIS数据提取复杂地形CFD模型的地表粗糙度长度数据，然后将一系列分布在等值线上的随机点的地表粗糙度长度转化为地表覆盖植被高度；3)计算地表覆盖植被在不同高度处对风场的阻力系数，然后将地表覆盖植被高度转换为相应的阻力系数；4)使用OpenFOAM添加阻力系数和地表覆盖植被阻力项，从而获得复杂地形风场特性模拟结果。2.根据权利要求1所述的一种基于OpenFOAM的模拟复杂地形地表粗糙度的方法，其特征在于，步骤1)具体包括以下子步骤：1.1)确定目标的中心经纬度，从GIS地理信息系统中获取目标所在区域的ASTERGDEM30米分辨率大地坐标下的地形数据，使用ArcMap处理GIS数据，从而获得WGS_1984_UTM坐标系下的地形数据网格节点；1.2)参数化建立三维平坦地形CFD模型并生成精细化网格：复杂地形CFD模型的计算域设置为长方体区域，根据实际的地形地貌找到对计算域有影响的复杂地形，计算域需要包含这些复杂地形，而且采用非结构化三棱柱网格以充分拟合复杂几何边界，其中，在计算域的不同区域采用分辨率不同的网格，复杂地形CFD模型在竖直方向上采用相邻网格尺寸比值为定值的σ网格；1.3)实现计算域的地形数据网格节点与三维平坦地形CFD模型网格节点之间的数据转换：根据三维平坦地形CFD模型中每一个底面网格数据节点的坐标，找到其在地形数据网格节点中的对应位置，然后使用该位置周围的3个地形数据网格节点的高程，通过三角形共面线性插值得到底面网格节点的高程，同时，在三维平坦地形CFD模型中，竖直方向上其余网格节点的高程会随着底面网格节点高程的变化产生相应变化，最终得到修改高程后的复杂地形CFD模型；1.4)通过转化器gambitToFoam导入OpenFOAM，将复杂地形CFD模型转换为OpenFOAM可识别的格式。3.根据权利要求1所述的一种基于OpenFOAM的模拟复杂地形地表粗糙度的方法，其特征在于，步骤2)具体包括以下子步骤：2.1)从GIS地理信息系统中提取复杂地形CFD模型的地表粗糙度长度数据...

【专利技术属性】
技术研发人员：刘震卿，吴晓波，熊世树，
申请(专利权)人：华中科技大学，
类型：发明
国别省市：湖北,42