CN113657055B - 一种复杂地形风场数值模拟的入流条件生成方法及*** - Google Patents
一种复杂地形风场数值模拟的入流条件生成方法及*** Download PDFInfo
- Publication number
- CN113657055B CN113657055B CN202110998569.1A CN202110998569A CN113657055B CN 113657055 B CN113657055 B CN 113657055B CN 202110998569 A CN202110998569 A CN 202110998569A CN 113657055 B CN113657055 B CN 113657055B
- Authority
- CN
- China
- Prior art keywords
- wind speed
- simulation
- area
- inflow
- main
- 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.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 189
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000004364 calculation method Methods 0.000 claims abstract description 36
- 230000001052 transient effect Effects 0.000 claims abstract description 28
- 238000013507 mapping Methods 0.000 claims abstract description 17
- 238000005070 sampling Methods 0.000 claims abstract description 15
- 238000004590 computer program Methods 0.000 claims description 16
- 230000003746 surface roughness Effects 0.000 claims description 8
- 230000000737 periodic effect Effects 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 239000000758 substrate Substances 0.000 claims 4
- 238000010586 diagram Methods 0.000 description 10
- 239000012530 fluid Substances 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 230000006872 improvement Effects 0.000 description 5
- 238000012876 topography Methods 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 1
- 239000002243 precursor Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- 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
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明一种复杂地形风场数值模拟的入流条件生成方法及***，该方法包括将稳态下的流动统计变量轮廓线进行采样后，映射至复杂地形风场所在主模拟区域的入流边界面处，在主模拟区域进行流场模拟计算；将主模拟区域内代表性测风塔处的模拟风速与实测平均风速进行对比；根据模拟风速与实测平均风速进行对比后的差值，反馈修正前导模拟的驱动风速参数，进而生成主模拟区域的稳态入流边界条件；拓展至大涡模拟类大气湍流模式，建立非定常瞬态风速的时间序列数据库，并以此生成主模拟区域的瞬态风速入流边界条件。本发明克服了现有的简化型入流边界条件对整个风场区域流动模拟精度与不确定度的影响，建立了适用于不同大气湍流模式下的入流条件生成方法。
Description
技术领域
本发明属于计算流体力学技术领域，具体涉及一种复杂地形风场数值模拟的入流条件生成方法及***。
背景技术
风资源储量丰富，潜在、在建与现役的大型陆上风电场大多坐落于山地、丘陵、高原等复杂地形下的风能聚集区域。计算流体力学(Computational Fluid Dynamics,CFD)技术作为现代计算技术贺大气动力学相结合的一种流体力学数值求解模式，已成为复杂地形区域风电项目开发过程中风能资源评估环节的重要工具。
风电场所在区域风资源模拟计算的几何建模范围，通常是以地表下垫面为基准面，竖直向上延伸至大气边界层高度(1至2公里量级)形成的六面体区域。在风场的数值模拟中，需要在六个边界面处对风速、气压等流动变量的边界条件进行合理的设置，从而确保流动控制方程在计算区域内的可解性。研究表明，作为整场大气流动模拟的驱动因子，入流界面处风速边界条件的生成方法直接影响着风场区域内风况模拟与风资源预测的精确度。
目前，风资源工程计算中有关复杂地形风场流动CFD数值模拟在入流界面处的风速边界条件，主要是基于求解流体力学控制方程组(即Navier-Stokes方程组)中的第一类边界条件(Dirichlet边界条件)进行直接或修正使用：(1)在整个入流边界面上的平均风速和风速标准差(或湍流强度)被赋予为固定值；(2)入流边界面上的平均风速设置为离地高度的空间函数，即平均风速值随高度的变化符合大气风速对数廓线(Logarithmicvelocity profile)理论模型，并给定合适量级的风速标准差(或湍流强度)。
综合而言，上述入流边界条件的设置方法在计算程序实现上较为便捷，但是当风场所在区域地形地貌较为复杂和崎岖，其入流界面底边为不规则几何边界时，风电场入流区域的实际风速在竖直方向的分布并非恒定值，或者并非完全符合大气对数风廓线理论模型，使用上述入流条件设置方法存在着对实际入流风况的过度简化问题，导致无法正确反映入流边界处实际的流动统计状态，对计算区域内风速求解的准确性造成影响；同时，使用均一化的平均风速与风速标准差边界条件设定方式时，难以对风场区域内代表性测风塔实地测得的风统计数据进行复现，在一定程度上增加了目标区域内的风资源模拟预测的不确定度。
发明内容
本发明的目的在于提供了一种复杂地形风场数值模拟的入流条件生成方法及***，该方法在求解山地、丘陵等复杂地形地貌上空大气流动时，通过利用前导模拟迭代方式与风场所在区域测风塔实测统计数据进行对比，反馈修正入流边界面处的流动变量边界条件，克服了现有的简化型入流边界条件对整个风场区域流动模拟精度与不确定度的影响。
为了达到上述目的，本发明采用的技术方案是：
一种复杂地形风场数值模拟的入流条件生成方法，包括以下步骤：
将以雷诺平均类大气湍流模式进行前导模拟得到的稳态下流动统计变量轮廓线进行采样后，映射至复杂地形风场所在主模拟区域的入流边界面处，进行主模拟区域的流场模拟计算；
将主模拟区域内代表性测风塔处的模拟风速与实测平均风速进行对比；根据模拟风速与实测平均风速的对比结果，反馈修正前导模拟的驱动风速参数，进而生成主模拟区域的稳态入流边界条件；
将稳态入流边界条件对应的前导模拟驱动风速参数拓展至大涡模拟类大气湍流模式，建立非定常瞬态风速的时间序列数据库，并以此生成主模拟区域的瞬态风速入流边界条件。
作为本发明的进一步改进，所述将稳态下的流动统计变量轮廓线进行采样后，映射至复杂地形风场所在主模拟区域的入流边界面处，主模拟区域的流场模拟计算；具体包括：
在风场所在复杂地形区域中确定大气流动数值模拟的六面体主模拟区域，将目标区域的复杂地形作为地表边界面；
统计目标时间段内，风电场区中代表性测风塔各高度测层处的平均风速与风速标准差；根据该区域的地形与地表植被覆盖情况，设定地表粗糙度；将平均风速作为前导模拟的初始驱动风速，结合地表粗糙度，进行前导模拟计算；当前导模拟结果达到稳定状态后停止计算，将稳定状态下流动变量随高度变化的轮廓线在计算主模拟区内进行采样，
将所述流动变量的轮廓线作为入流边界条件，映射到主模拟区域的入流边界上，在主模拟区域内继续采用雷诺平均类湍流模型进行数值计算，待模拟结果达到统计稳定状态后，对代表性测风塔所处位置在轮毂高度层的模拟风速进行采样记录。
作为本发明的进一步改进，所述前导模拟的几何区域为标准的六面体，下垫面为平面，长度和宽度为公里量级，高度为大气边界层高度；其与下垫面垂直的4个周向边界面处均设置为周期性边界条件。
作为本发明的进一步改进，所述将主模拟区域内代表性测风塔处的模拟风速与实测平均风速进行对比，根据模拟风速与实测平均风速的对比结果，反馈修正前导模拟的驱动风速参数，进而生成主模拟区域的稳态入流边界条件具体包括：
若模拟风速与实测平均风速差值大于阈值,则返回将前导模拟的初始驱动风速值进行调整：
当模拟风速大于实测平均风速，则将初始驱动风速的初始值调小；当模拟风速小于实测平均风速，则将初始驱动风速调大；
若模拟风速与实测平均风速差值小于阈值，则记录最终的前导模拟的驱动风速值对应流动变量轮廓线特征，将流动变量轮廓线最终映射到主模拟区域入流边界上，作为雷诺平均类湍流模式下流动统计变量的稳态入流边界条件。
作为本发明的进一步改进，所述将稳态入流边界条件对应的前导模拟驱动风速参数拓展至大涡模拟类大气湍流模式，建立非定常瞬态风速的时间序列数据库，具体包括：
采用大涡模拟类大气湍流模式，选取雷诺平均类湍流模式下最终得到的驱动风速，并结合设定的地表粗糙度，在前导模拟计算主模拟区域中的4个周向边界面采用周期性边界条件，进行前导模拟；
当前导模拟结果达到统计稳定状态后，选定前导模拟计算域主模拟区域中与主流动方向垂直的任一截面，将截面区域的风速时程数据进行存储，建立统计稳定状态下，非定常瞬态风速的时间序列数据库。
作为本发明的进一步改进，所述统计稳定状态是指流场变量在不同时间段内的平均值不再发生变化，而其瞬时值依旧是在平均值附近波动变化的非定常状态。
作为本发明的进一步改进，所述并以此生成主模拟区域的瞬态风速入流边界条件，具体包括：
将所述非定常瞬态风速时序数据库中的截面风速时程数据，基于主模拟区域中入流边界面各空间点位距离复杂地表的实际高程，逐时间步长地提取并映射到主模拟区域的入流边界面上，形成大涡模拟类湍流模式下瞬态的风速入流边界条件。
一种复杂地形风场数值模拟的入流条件生成***，包括：
模拟风速计算模块，用于将以雷诺平均类大气湍流模式进行前导模拟得到的稳态下流动统计变量轮廓线进行采样后，映射至复杂地形风场所在主模拟区域的入流边界面处，进行主模拟区域的流场模拟计算；
对比修正模块，用于将主模拟区域内代表性测风塔处的模拟风速与实测平均风速进行对比；根据模拟风速与实测平均风速的对比结果，反馈修正前导模拟的驱动风速参数，进而生成主模拟区域的稳态入流边界条件；
边界条件生成模块，用于将稳态入流边界条件对应的前导模拟驱动风速参数拓展至大涡模拟类大气湍流模式，建立非定常瞬态风速的时间序列数据库，并以此生成主模拟区域的瞬态风速入流边界条件。
一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述复杂地形风场数值模拟的入流条件生成的步骤。
一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现所述复杂地形风场数值模拟的入流条件生成的步骤。
与现有技术方法相比，本发明的有益效果是：
本发明提供的一种复杂地形风场数值模拟的入流条件生成方法，可以通过前导模拟迭代修正的方式，最大程度地模拟复现出风电场区域内代表性测风塔处测量得到的风统计特征，降低了风资源数值模拟的不确定度；通过前导模拟生成的入流边界条件，能够考虑到复杂地形对近地层区域流动的局部改变效应，生成的流场特征在入流边界面附近更加符合实际大气流动状态，在一定程度上减少了均一化的平均风速与风速标准差边界条件设定方式对整个风场区域模拟的精确度影响。本发明较好地克服了现有的简化型入流边界条件对整个风场区域流动模拟精度与不确定度的影响，建立了与风场流动实际情况更为接近的入流条件生成方法。
进一步地，本发明提供的风场数值模拟的入流条件生成方法，对于两种常用的大气湍流模式，即雷诺平均(Reynolds-Averaged Navier-Stokes，RANS)类模式和大涡模拟(Large-eddy Simulation，LES)类模式均可适用，故具有较强的通用性。
附图说明
图1为本发明所涉及的一种复杂地形风场数值模拟的入流条件生成方法的流程图；
图2为RANS类大气湍流模式下，复杂地形风场数值模拟的入流条件生成方法示意图；
图3为LES类大气湍流模式下，复杂地形风场数值模拟的入流条件生成方法示意图；
图4为本发明复杂地形风场数值模拟的入流条件生成***结构示意图；
图5为本发明电子设备结构示意图。
具体实施方式
下面将参考附图并结合实施例来详细说明本发明。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。
以下详细说明均是示例性的说明，旨在对本发明提供进一步的详细说明。除非另有指明，本发明所采用的所有技术术语与本申请所属领域的一般技术人员的通常理解的含义相同。本发明所使用的术语仅是为了描述具体实施方式，而并非意图限制根据本发明的示例性实施方式。
本发明第一个目的是提供一种复杂地形风场数值模拟的入流条件生成方法，该方法包括：
以计算模拟效率较高的雷诺平均(RANS)类大气湍流模式为基础，通过利用前导模拟迭代的方式，将稳态下的流动统计变量轮廓线进行采样后，映射至复杂地形风场所在主模拟区域的入流边界面处，进行数值计算，之后将场区域内代表性测风塔处的模拟风速与实测平均风速进行对比，反馈修正(并生成)入流边界面处流动统计变量的稳态边界条件；进而拓展至精度较高的大涡模拟(LES)类大气湍流模式下，非定常瞬态风速的时间序列数据库的建立，并以此生成主模拟区域的瞬态风速入流边界条件。
本发明的较好地克服了现有的简化型入流边界条件对整个风场区域流动模拟精度与不确定度的影响，建立了与风场流动实际情况更为接近、适用于常用大气湍流模式的入流条件生成方法。
下面结合附图，对本发明进一步详细说明。
实施例
如图1所示，本发明提供的一种复杂地形风场数值模拟的入流条件生成方法，包括如下步骤：
步骤1，在风场所在复杂地形区域中确定大气流动数值模拟的六面体计算域DT(主模拟区域)，将目标区域的复杂地形作为地表边界面，入流边界面的设置应垂直于该区域的主(盛行)风向；
步骤2，统计特定(目标)时间段内，风电场区中代表性测风塔各高度测层处(包括风力发电机轮毂高度测层)的平均风速与风速标准差；重点计算测风塔在风力机轮毂高度测层(hub level)处的实测平均风速
步骤3，根据该区域的地形与地表植被覆盖情况，设定地表粗糙度z0；
其中，设定地表粗糙度为本领域技术人员的通识性方法，在此不进行详细说明。
步骤4，将步骤2中的作为前导模拟(precursor simulation)的初始驱动风速/>(即将/>的值直接赋给/>)，结合步骤3中设定的z0，进行前导模拟；所述前导模拟中大气湍流模式采用RANS类中经典的k-ε模型，其中k为湍动能，ε为湍动能的耗散率；
前导模拟为业界为本领域技术人员的通识性方法，在此不进行详细说明，这里只论述其关键参数的设置)。
所述前导模拟的几何计算区域Dp为标准的六面体，下垫面为平面，长度和宽度为公里量级(通常大于3km)，高度为1至2公里量级(大气边界层高度)；其与下垫面垂直的4个周向边界面处均设置为周期性边界条件(Periodic Boundary Conditions，PBS)进行前导模拟，如图2(a)所示；
步骤5，当步骤4所述的前导模拟结果达到定常(稳定)状态(steady state)后停止计算，将稳定状态下的风速U、湍动能k，湍动能耗散率ε等流动变量随高度变化的轮廓线在计算域Dp内进行采样，如图2(b)所示；
步骤6，将步骤5中所述流动变量(U、k、ε)的轮廓线作为入流边界条件，映射(mapping)到主模拟区域DT的入流边界上，在主模拟区域DT内继续采用RANS类湍流模型进行数值计算，待模拟结果达到统计稳定状态后，对代表性测风塔所处位置在轮毂高度层的模拟风速进行采样记录；
步骤7，将步骤6中所述和步骤4中所述/>进行对比：
(1)若二者之间存在较大误差则返回步骤4，将/>值进行调整：当/>远大于/>则将/>较步骤4中的初始值调小；当/>远小于则将/>较步骤4中的初始值调大。
(2)若二者之间的误差在可接受范围之内则记录下步骤4中最终的驱动风速/>值对应步骤6中的流动变量(U、k、ε)轮廓线特征，将流动变量(U、k、ε)轮廓线最终映射到主模拟区域DT入流边界上，作为RANS类湍流模式下流动统计变量的稳态(定常)入流边界条件，如图2(c)所示；
步骤8，当风场数值计算需要采用大涡模拟(LES)类大气湍流模式时，可选取步骤4中最终得到的驱动风速并结合步骤2中设定的z0，在前导模拟计算域Dp中的4个周向边界面采用周期性边界条件，进行前导模拟，如图3(a)所示；
步骤9，当步骤8所述的前导模拟结果达到统计稳定状态(statistically steadystate)后，选定前导模拟计算域Dp中与主流动方向垂直的任一截面，将该截面区域的风速时程数据进行存储，建立统计稳定状态下，非定常瞬态风速的时间序列“数据库”，如图3(b)所示。
所述统计稳定状态是指流场变量在不同时间段内的***(时、空间)平均值不再发生变化，而其瞬时值依旧是在平均值附近波动变化的非定常(unsteady)状态。
步骤10，将步骤9所述非定常瞬态风速时序“数据库”中的截面(空间)风速时程数据，基于主模拟区域DT中入流边界面各空间点位距离复杂地表的实际高程，逐时间步长地提取并映射到DT的入流边界面上，形成LES湍流模式下，瞬态(非定常)的风速入流边界条件，如图3(c)所示。该入流边界条件为复杂地形上空流场模拟提供了可进行能谱、时间和空间相关统计特性分析的湍流漩涡结构，更加接近真实大气湍流的入流特性。
以上所述，仅为本发明较佳的指导性实施方式，本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明所述的技术范围内，可轻易想到的以该入流边界条件生成方法为核心，对雷诺平均(RANS)、大涡模拟(LES)两大类湍流模式中具体模型以及其它模拟参数的替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权力要求书的保护范围为准。
如图4所示，本发明第二个目的是提供一种复杂地形风场数值模拟的入流条件生成***，包括：
模拟风速计算模块，用于将以雷诺平均(RANS)类大气湍流模式进行前导模拟得到的稳态流动统计变量轮廓线进行采样后，映射至复杂地形风场所在主模拟区域的入流边界面处，进行主模拟区域的流场模拟计算；
对比修正模块，用于将主模拟区域内代表性测风塔处的模拟风速与实测平均风速进行对比；根据模拟风速与实测平均风速进行对比后，反馈修正前导模拟的驱动风速参数，并生成主模拟区域的稳态入流边界条件；
边界条件生成模块，用于将稳态入流边界条件对应的前导模拟驱动风速参数拓展至大涡模拟(LES)类大气湍流模式，建立非定常瞬态风速的时间序列数据库，并以此生成主模拟区域的瞬态风速入流边界条件。
如图5所示，本发明第三个目的是提供一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现所述复杂地形风场数值模拟的入流条件生成的步骤。
本发明第四个目的是提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现所述复杂地形风场数值模拟的入流条件生成的步骤。
本领域内的技术人员应明白，本发明的实施例可提供为方法、***、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。
Claims (8)
1.一种复杂地形风场数值模拟的入流条件生成方法，其特征在于，包括以下步骤：
将以雷诺平均类大气湍流模式进行前导模拟得到的稳态下流动统计变量轮廓线进行采样后，映射至复杂地形风场所在主模拟区域的入流边界面处，进行主模拟区域的流场模拟计算；
将主模拟区域内代表性测风塔处的模拟风速与实测平均风速进行对比；根据模拟风速与实测平均风速的对比结果，反馈修正前导模拟的驱动风速参数，进而生成主模拟区域的稳态入流边界条件；
将稳态入流边界条件对应的前导模拟驱动风速参数拓展至大涡模拟类大气湍流模式，建立非定常瞬态风速的时间序列数据库，并以此生成主模拟区域的瞬态风速入流边界条件；
所述将稳态下的流动统计变量轮廓线进行采样后，映射至复杂地形风场所在主模拟区域的入流边界面处，主模拟区域的流场模拟计算；具体包括：
在风场所在复杂地形区域中确定大气流动数值模拟的六面体主模拟区域，将目标区域的复杂地形作为地表边界面；
统计目标时间段内，风电场区中代表性测风塔各高度测层处的平均风速与风速标准差；根据该区域的地形与地表植被覆盖情况，设定地表粗糙度；将平均风速作为前导模拟的初始驱动风速，结合地表粗糙度，进行前导模拟计算；当前导模拟结果达到稳定状态后停止计算，将稳定状态下流动变量随高度变化的轮廓线在计算主模拟区内进行采样，
将所述流动变量的轮廓线作为入流边界条件，映射到主模拟区域的入流边界上，在主模拟区域内继续采用雷诺平均类湍流模型进行数值计算，待模拟结果达到统计稳定状态后，对代表性测风塔所处位置在轮毂高度层的模拟风速进行采样记录；
所述将主模拟区域内代表性测风塔处的模拟风速与实测平均风速进行对比，根据模拟风速与实测平均风速的对比结果，反馈修正前导模拟的驱动风速参数，进而生成主模拟区域的稳态入流边界条件具体包括：
若模拟风速与实测平均风速差值大于阈值,则返回将前导模拟的初始驱动风速值进行调整：
当模拟风速大于实测平均风速，则将初始驱动风速的初始值调小；当模拟风速小于实测平均风速，则将初始驱动风速调大；
若模拟风速与实测平均风速差值小于阈值，则记录最终的前导模拟的驱动风速值对应流动变量轮廓线特征，将流动变量轮廓线最终映射到主模拟区域入流边界上，作为雷诺平均类湍流模式下流动统计变量的稳态入流边界条件。
2.根据权利要求1所述的方法，其特征在于，
所述前导模拟的几何区域为标准的六面体，下垫面为平面，长度和宽度为公里量级，高度为大气边界层高度；其与下垫面垂直的4个周向边界面处均设置为周期性边界条件。
3.根据权利要求1所述的方法，其特征在于，
所述将稳态入流边界条件对应的前导模拟驱动风速参数拓展至大涡模拟类大气湍流模式，建立非定常瞬态风速的时间序列数据库，具体包括：
采用大涡模拟类大气湍流模式，选取雷诺平均类湍流模式下最终得到的驱动风速，并结合设定的地表粗糙度，在前导模拟计算主模拟区域中的4个周向边界面采用周期性边界条件，进行前导模拟；
当前导模拟结果达到统计稳定状态后，选定前导模拟计算域主模拟区域中与主流动方向垂直的任一截面，将截面区域的风速时程数据进行存储，建立统计稳定状态下，非定常瞬态风速的时间序列数据库。
4.根据权利要求1所述的方法，其特征在于，
所述统计稳定状态是指流场变量在不同时间段内的平均值不再发生变化，而其瞬时值依旧是在平均值附近波动变化的非定常状态。
5.根据权利要求1所述的方法，其特征在于，
所述并以此生成主模拟区域的瞬态风速入流边界条件，具体包括：
将所述非定常瞬态风速时序数据库中的截面风速时程数据，基于主模拟区域中入流边界面各空间点位距离复杂地表的实际高程，逐时间步长地提取并映射到主模拟区域的入流边界面上，形成大涡模拟类湍流模式下瞬态的风速入流边界条件。
6.一种复杂地形风场数值模拟的入流条件生成***，基于权利要求1至5任一项所述的方法，其特征在于，包括：
模拟风速计算模块，用于将以雷诺平均类大气湍流模式进行前导模拟得到的稳态下流动统计变量轮廓线进行采样后，映射至复杂地形风场所在主模拟区域的入流边界面处，进行主模拟区域的流场模拟计算；
对比修正模块，用于将主模拟区域内代表性测风塔处的模拟风速与实测平均风速进行对比；根据模拟风速与实测平均风速的对比结果，反馈修正前导模拟的驱动风速参数，进而生成主模拟区域的稳态入流边界条件；
边界条件生成模块，用于将稳态入流边界条件对应的前导模拟驱动风速参数拓展至大涡模拟类大气湍流模式，建立非定常瞬态风速的时间序列数据库，并以此生成主模拟区域的瞬态风速入流边界条件。
7.一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现权利要求1-5任一项所述复杂地形风场数值模拟的入流条件生成的步骤。
8.一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现权利要求1-5任一项所述复杂地形风场数值模拟的入流条件生成的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110998569.1A CN113657055B (zh) | 2021-08-27 | 2021-08-27 | 一种复杂地形风场数值模拟的入流条件生成方法及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110998569.1A CN113657055B (zh) | 2021-08-27 | 2021-08-27 | 一种复杂地形风场数值模拟的入流条件生成方法及*** |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113657055A CN113657055A (zh) | 2021-11-16 |
| CN113657055B true CN113657055B (zh) | 2024-05-07 |
Family
ID=78482310
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110998569.1A Active CN113657055B (zh) | 2021-08-27 | 2021-08-27 | 一种复杂地形风场数值模拟的入流条件生成方法及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113657055B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114595643B (zh) * | 2022-03-04 | 2024-06-07 | 重庆大学 | 耦合多点测风仪和微尺度稳态模拟的高分辨率山地风场测量方法 |
| CN116070414B (zh) * | 2022-12-12 | 2023-10-17 | 中广核风电有限公司 | 基于地表起伏复杂程度的风资源模拟方法及装置 |
| CN116341410B (zh) * | 2023-02-07 | 2024-01-30 | 成都流体动力创新中心 | 一种考虑不同稳定度的典型地形风场模拟方法及*** |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012105973A1 (en) * | 2011-02-02 | 2012-08-09 | Michigan Aerospace Corporation | Atmospheric measurement system and method |
| CN107729638A (zh) * | 2017-10-09 | 2018-02-23 | 中国民航大学 | 各向异性大气湍流场数值模拟方法 |
| CN112163381A (zh) * | 2020-09-25 | 2021-01-01 | 西安热工研究院有限公司 | 一种适用于复杂地形风场流动数值模拟的侧向边界条件设置方法 |
| CN112949227A (zh) * | 2021-03-30 | 2021-06-11 | 中国华能集团清洁能源技术研究院有限公司 | 适用于复杂地形风场湍流强度确定方法、***、设备及存储介质 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11281826B2 (en) * | 2016-12-21 | 2022-03-22 | Northwestern University | Systems and methods for computational simulation of self-propelling vehicles for aerodynamic design |
| WO2021087011A1 (en) * | 2019-10-28 | 2021-05-06 | Vorcat, Inc. | System and method for simulating turbulence |
-
2021
- 2021-08-27 CN CN202110998569.1A patent/CN113657055B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012105973A1 (en) * | 2011-02-02 | 2012-08-09 | Michigan Aerospace Corporation | Atmospheric measurement system and method |
| CN107729638A (zh) * | 2017-10-09 | 2018-02-23 | 中国民航大学 | 各向异性大气湍流场数值模拟方法 |
| CN112163381A (zh) * | 2020-09-25 | 2021-01-01 | 西安热工研究院有限公司 | 一种适用于复杂地形风场流动数值模拟的侧向边界条件设置方法 |
| CN112949227A (zh) * | 2021-03-30 | 2021-06-11 | 中国华能集团清洁能源技术研究院有限公司 | 适用于复杂地形风场湍流强度确定方法、***、设备及存储介质 |
Non-Patent Citations (1)
| Title |
|---|
| 基于数值模拟的复杂地形风场风资源评估方法;梁思超;张晓东;康顺;康雅兰;赵永锋;;空气动力学学报;20120615(第03期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113657055A (zh) | 2021-11-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113657055B (zh) | 一种复杂地形风场数值模拟的入流条件生成方法及*** | |
| AU2017352549B2 (en) | Method and device for calculating power generation of wind farm | |
| Porté-Agel et al. | Large-eddy simulation of atmospheric boundary layer flow through wind turbines and wind farms | |
| US11713747B2 (en) | Wind turbine control method | |
| CN109086534B (zh) | 一种基于cfd流体力学模型的风电场尾流订正方法及*** | |
| CN110088463B (zh) | 风力涡轮电场级载荷管理控制 | |
| CN112459965B (zh) | 考虑风电场尾流的偏航优化控制方法、装置、设备及介质 | |
| CN102663251B (zh) | 基于计算流体力学模型的风电场功率物理预测方法 | |
| CN107194097A (zh) | 基于风电场气动模拟和风速风向数据的分析方法 | |
| CN110414102A (zh) | 海上风力发电结构疲劳分析方法 | |
| Han et al. | Large eddy simulation for atmospheric boundary layer flow over flat and complex terrains | |
| CN106815773A (zh) | 一种风电场功率特性评估方法 | |
| CN108717593A (zh) | 一种基于风轮面等效风速的微观选址发电量评估方法 | |
| CN112163381B (zh) | 一种适用于复杂地形风场流动数值模拟的侧向边界条件设置方法 | |
| Diaz et al. | Full wind rose wind farm simulation including wake and terrain effects for energy yield assessment | |
| Agarwal | Structural reliability of offshore wind turbines | |
| Branlard et al. | A digital twin solution for floating offshore wind turbines validated using a full-scale prototype | |
| Creech et al. | Modeling wind turbine wakes for wind farms | |
| CN114154325A (zh) | 一种定常与非定常混合风能资源评估方法及*** | |
| Li et al. | Wind velocity prediction at wind turbine hub height based on CFD model | |
| CN117313399B (zh) | 一种适用于复杂地形的水平轴风力机三维各向异性超高斯全尾流模型的建立及应用方法 | |
| Ichenial et al. | Investigation of the impact of planetary boundary layer turbulence on the performance of onshore wind farm | |
| Ivanell et al. | The impact of wind direction in atmospheric BL on interacting wakes at Horns Rev wind farm | |
| Natarajan et al. | Demonstration of a basis for Tall Wind Turbine Design, EUDP project final report | |
| Li et al. | Simulation-based micro-site selection for wind farm in complex terrain |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |