CN113657055A - Inflow condition generation method and system for numerical simulation of wind field in complex terrain - Google Patents
Inflow condition generation method and system for numerical simulation of wind field in complex terrain Download PDFInfo
- Publication number
- CN113657055A CN113657055A CN202110998569.1A CN202110998569A CN113657055A CN 113657055 A CN113657055 A CN 113657055A CN 202110998569 A CN202110998569 A CN 202110998569A CN 113657055 A CN113657055 A CN 113657055A
- 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.)
- Granted
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 183
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 36
- 230000001052 transient effect Effects 0.000 claims abstract description 27
- 238000013507 mapping Methods 0.000 claims abstract description 18
- 238000005070 sampling Methods 0.000 claims abstract description 15
- 238000004590 computer program Methods 0.000 claims description 16
- 230000000737 periodic effect Effects 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 230000003746 surface roughness Effects 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 5
- 238000005259 measurement Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 12
- 230000006872 improvement Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 239000012530 fluid Substances 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
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process 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
- 238000012876 topography Methods 0.000 description 1
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
- 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
The invention relates to an inflow condition generation method and a system for numerical simulation of a complex terrain wind field, wherein the method comprises the steps of sampling a flow statistical variable contour line under a steady state, mapping the flow statistical variable contour line to an inflow boundary surface of a main simulation area where the complex terrain wind field is located, and performing flow field simulation calculation in the main simulation area; comparing the simulated wind speed at the representative anemometer tower in the main simulation area with the actually measured average wind speed; feeding back and correcting the drive wind speed parameters of the leading simulation according to the difference value after the comparison between the simulated wind speed and the actually measured average wind speed, and further generating a steady-state inflow boundary condition of the main simulation area; and expanding the model to a large vortex simulation type atmospheric turbulence mode, establishing a time sequence database of unsteady transient wind speed, and generating a transient wind speed inflow boundary condition of a main simulation area according to the time sequence database. The invention overcomes the influence of the traditional simplified inflow boundary condition on the flow simulation precision and uncertainty of the whole wind field area, and establishes the inflow condition generation method suitable for different atmospheric turbulence modes.
Description
Technical Field
The invention belongs to the technical field of computational fluid mechanics, and particularly relates to an inflow condition generation method and system for numerical simulation of a wind field in a complex terrain.
Background
The wind resource reserves are abundant, and the potential large-scale land wind power plants under construction and active service mostly sit in wind energy gathering areas under complex terrains such as mountainous regions, hills and plateaus. Computational Fluid Dynamics (CFD) technology is used as a hydrodynamics numerical solving mode combined with modern computing technology and atmospheric Dynamics, and has become an important tool in a wind energy resource evaluation link in the development process of wind power projects in complicated terrain areas.
The geometric modeling range of wind resource simulation calculation of the area where the wind power plant is located is a hexahedral area which is usually formed by taking an underground bedding surface as a reference surface and vertically extending upwards to the height (the order of 1 to 2 kilometers) of an atmospheric boundary layer. In the numerical simulation of the wind field, boundary conditions of flow variables such as wind speed and air pressure need to be reasonably set at six boundary surfaces, so that the solvability of a flow control equation in a calculation region is ensured. Research shows that as a driving factor for simulating the whole atmospheric flow, the generation method of the wind speed boundary condition at the inflow interface directly influences the accuracy of wind condition simulation and wind resource prediction in the wind field area.
At present, in wind resource engineering calculation, a wind speed boundary condition at an inflow interface is numerically simulated in relation to a wind field flow CFD of a complex terrain, and the boundary condition is mainly used by solving a first type of boundary condition (Dirichlet boundary condition) in a hydrodynamics control equation set (namely a Navier-Stokes equation set) directly or in a modified mode: (1) the average wind speed and the standard deviation of the wind speed (or turbulence intensity) over the entire inflow boundary surface are given as fixed values; (2) the average wind speed on the inflow boundary surface is set as a space function of the ground clearance height, namely, the change of the average wind speed value along with the height conforms to the theoretical model of the log-log profile of the atmospheric wind speed (Logiathmic velocity profile), and the wind speed standard deviation (or turbulence intensity) with proper magnitude is given.
Comprehensively, the setting method of the inflow boundary condition is convenient and fast in calculation program realization, but when the terrain and the topography of the area where the wind field is located are complex and rugged, and the bottom edge of the inflow interface of the wind field is an irregular geometric boundary, the distribution of the actual wind speed of the inflow area of the wind power plant in the vertical direction is not a constant value or does not completely accord with the theoretical model of the atmospheric logarithmic wind profile, so that the problem of excessive simplification of the actual inflow wind condition exists by using the setting method of the inflow condition, the actual flow statistical state at the inflow boundary cannot be correctly reflected, and the accuracy of the wind speed solution in the calculation area is influenced; meanwhile, when a uniform boundary condition setting mode of average wind speed and standard deviation of wind speed is used, the statistical data of wind measured on the spot by a representative anemometer tower in a wind field area is difficult to reproduce, and the uncertainty of wind resource simulation prediction in a target area is increased to a certain extent.
Disclosure of Invention
The invention aims to provide an inflow condition generation method and system for numerical simulation of a complex terrain wind field.
In order to achieve the purpose, the invention adopts the technical scheme that:
an inflow condition generation method for numerical simulation of a wind field of a complex terrain comprises the following steps:
sampling a flow statistical variable contour line under a steady state obtained by performing leading simulation in a Reynolds average atmospheric turbulence mode, mapping the contour line to an inflow boundary surface of a main simulation area where a complex terrain wind field is located, and performing flow field simulation calculation of the main simulation area;
comparing the simulated wind speed at the representative anemometer tower in the main simulation area with the actually measured average wind speed; according to the comparison result of the simulated wind speed and the actually measured average wind speed, the driving wind speed parameters of the leading simulation are fed back and corrected, and then the steady-state inflow boundary condition of the main simulation area is generated;
and expanding the leading simulation driving wind speed parameters corresponding to the steady-state inflow boundary conditions to a large vortex simulation type atmospheric turbulence mode, establishing a time sequence database of unsteady transient wind speeds, and generating transient wind speed inflow boundary conditions of a main simulation area according to the time sequence database.
As a further improvement of the method, after sampling the flow statistical variable contour line under the steady state, mapping the flow statistical variable contour line to an inflow boundary surface of a main simulation area where a complex terrain wind field is located, and performing flow field simulation calculation of the main simulation area; the method specifically comprises the following steps:
determining a hexahedron main simulation area for numerical simulation of atmospheric flow in a complex terrain area where a wind field is located, and taking the complex terrain of a target area as a ground surface boundary surface;
counting the average wind speed and the standard deviation of the wind speed at each height measurement layer of a representative wind measuring tower in a wind power plant area within a target time period; setting the roughness of the earth surface according to the terrain of the area and the vegetation coverage condition of the earth surface; taking the average wind speed as the initial driving wind speed of the leading simulation, and combining the surface roughness to perform leading simulation calculation; stopping calculation when the leading simulation result reaches a stable state, sampling the contour line of the flow variable changing along with the height in the main calculation simulation area,
and mapping the contour line of the flow variable as an inflow boundary condition to an inflow boundary of a main simulation area, continuously adopting a Reynolds average turbulence model to perform numerical calculation in the main simulation area, and sampling and recording the simulated wind speed of the representative anemometer tower at the hub height layer after the simulation result reaches a statistical stable state.
As a further improvement of the invention, the geometric area of the preamble simulation is a standard hexahedron, the underlying surface is a plane, the length and the width are kilometer orders, and the height is the height of an atmospheric boundary layer; periodic boundary conditions were set at 4 circumferential boundary surfaces thereof perpendicular to the underlying surface.
As a further improvement of the present invention, the step of comparing the simulated wind speed at the representative anemometer tower in the main simulation area with the measured average wind speed, and feeding back and correcting the driving wind speed parameter of the leading simulation according to the comparison result of the simulated wind speed and the measured average wind speed, so as to generate the steady-state inflow boundary condition of the main simulation area specifically includes:
if the difference value between the simulated wind speed and the actually measured average wind speed is larger than the threshold value, returning to adjust the initial driving wind speed value of the leading simulation:
when the simulated wind speed is larger than the actually measured average wind speed, the initial value of the initial driving wind speed is reduced; when the simulated wind speed is smaller than the actually measured average wind speed, the initial driving wind speed is increased;
if the difference value between the simulated wind speed and the actually measured average wind speed is smaller than the threshold value, recording the flow variable contour line characteristics corresponding to the final leading simulated driving wind speed value, and finally mapping the flow variable contour line to the inflow boundary of the main simulation area to be used as the steady-state inflow boundary condition of the flow statistical variable in the Reynolds average turbulence-like mode.
As a further improvement of the present invention, the expanding the leading simulation driving wind speed parameter corresponding to the steady-state inflow boundary condition to a large vortex simulation type atmospheric turbulence mode and establishing a time series database of unsteady transient wind speeds specifically includes:
adopting a large vortex simulation type atmospheric turbulence mode, selecting the finally obtained driving wind speed under a Reynolds average type turbulence mode, and performing leading simulation by combining set surface roughness and adopting periodic boundary conditions on 4 circumferential boundary surfaces in a leading simulation calculation main simulation area;
and when the leading simulation result reaches a statistical stable state, selecting any section perpendicular to the main flow direction in the main simulation area of the leading simulation calculation area, storing the wind speed time course data of the section area, and establishing a time sequence database of the unsteady transient wind speed in the statistical stable state.
As a further improvement of the present invention, the statistically stable state refers to an unsteady state in which the average value of the flow field variable no longer changes in different time periods, and the instantaneous value thereof still fluctuates around the average value.
As a further improvement of the present invention, the generating the transient wind speed inflow boundary condition of the main simulation region includes:
and extracting and mapping the section wind speed time-course data in the unsteady transient wind speed time-series database on the inflow boundary surface of the main simulation area time-step by time-step based on the actual elevation of each space point of the inflow boundary surface in the main simulation area from the complex earth surface to form the transient wind speed inflow boundary condition in the large vortex simulation turbulence mode.
An inflow condition generating system for numerical simulation of a wind field in complex terrain, comprising:
the simulated wind speed calculation module is used for sampling a flow statistical variable contour line under a steady state obtained by performing leading simulation in a Reynolds average atmospheric turbulence mode, mapping the flow statistical variable contour line to an inflow boundary surface of a main simulation area where a wind field with complex terrain is located, and performing flow field simulation calculation of the main simulation area;
the comparison and correction module is used for comparing the simulated wind speed at the representative anemometer tower in the main simulation area with the actually measured average wind speed; according to the comparison result of the simulated wind speed and the actually measured average wind speed, the driving wind speed parameters of the leading simulation are fed back and corrected, and then the steady-state inflow boundary condition of the main simulation area is generated;
and the boundary condition generating module is used for expanding the leading simulation driving wind speed parameter corresponding to the steady-state inflow boundary condition to a large vortex simulation type atmospheric turbulence mode, establishing a time sequence database of unsteady transient wind speed, and generating the transient wind speed inflow boundary condition of the main simulation area according to the time sequence database.
An electronic device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of inflow condition generation for numerical simulation of the complex terrain wind field when executing the computer program.
A computer-readable storage medium, storing a computer program which, when executed by a processor, performs the steps of inflow condition generation for the complex terrain wind farm numerical simulation.
Compared with the prior art, the method has the beneficial effects that:
according to the inflow condition generation method for numerical simulation of the wind field in the complex terrain, the wind statistical characteristics measured at the representative wind measuring tower in the wind field area can be simulated and reproduced to the greatest extent in a lead simulation iterative correction mode, and the uncertainty of numerical simulation of wind resources is reduced; the inflow boundary condition generated by the leading simulation can consider the local change effect of the complex terrain on the flow of the near-stratum region, the generated flow field characteristics are more consistent with the actual atmospheric flow state near the inflow boundary surface, and the influence of the setting mode of the boundary condition of the uniform average wind speed and the standard deviation of the wind speed on the simulation accuracy of the whole wind field region is reduced to a certain extent. The invention better overcomes the influence of the traditional simplified inflow boundary condition on the flow simulation precision and uncertainty of the whole wind field area, and establishes an inflow condition generation method which is closer to the actual flow condition of the wind field.
Furthermore, the inflow condition generation method for wind field numerical Simulation provided by the invention is applicable to two common atmospheric turbulence modes, namely a Reynolds-average Navier-Stokes (RANS) mode and a Large vortex Simulation (LES) mode, so that the method has strong universality.
Drawings
FIG. 1 is a flow chart of a method for generating inflow conditions for numerical simulation of a wind field in a complex terrain in accordance with the present invention;
FIG. 2 is a schematic diagram of an inflow condition generation method of complex terrain wind field numerical simulation in an RANS type atmospheric turbulence mode;
FIG. 3 is a schematic diagram of an inflow condition generation method for numerical simulation of a wind field in a complex terrain in an LES-type atmospheric turbulence mode;
FIG. 4 is a schematic structural diagram of an inflow condition generation system for wind field numerical simulation of complex terrain in accordance with the present invention;
fig. 5 is a schematic structural diagram of an electronic device according to the present invention.
Detailed Description
The present invention will be described in detail below with reference to the embodiments with reference to the attached drawings. It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict.
The following detailed description is exemplary in nature and is intended to provide further details of the invention. Unless otherwise defined, all technical terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the invention.
The invention aims to provide a method for generating inflow conditions for numerical simulation of a wind field of a complex terrain, which comprises the following steps:
on the basis of a Reynolds Average (RANS) type atmospheric turbulence mode with high calculation simulation efficiency, sampling a flow statistic variable contour line under a steady state by utilizing a lead simulation iteration mode, mapping the flow statistic variable contour line to an inflow boundary surface of a main simulation area of a wind field with a complex terrain, performing numerical calculation, comparing a simulated wind speed at a representative anemometer tower in the field with an actually measured average wind speed, and performing feedback correction (generating) on a steady state boundary condition of the flow statistic variable at the inflow boundary surface; and further expanding the method to the establishment of a time sequence database of unsteady transient wind speed under a high-precision Large Eddy Simulation (LES) type atmospheric turbulence mode, and generating transient wind speed inflow boundary conditions of a main simulation area according to the time sequence database.
The invention better overcomes the influence of the traditional simplified inflow boundary condition on the flow simulation precision and uncertainty of the whole wind field area, and establishes the inflow condition generation method which is more approximate to the actual flow condition of the wind field and is suitable for the common atmospheric turbulence mode.
The present invention will be described in further detail with reference to the accompanying drawings.
Examples
As shown in fig. 1, the method for generating an inflow condition of a complex terrain wind field numerical simulation provided by the invention comprises the following steps:
step 1, determining a hexahedron calculation domain D of atmospheric flow numerical simulation in a complex terrain area where a wind field is locatedT(main simulation area), taking the complex terrain of the target area as a surface boundary surface, and setting the inflow boundary surface to be vertical to the main (prevailing) wind direction of the area;
step 2, counting the average wind speed and the standard deviation of the wind speed at each height measurement layer (including the height measurement layer of a hub of a wind driven generator) of a representative anemometer tower in a wind electric field area in a specific (target) time period; the actually measured average wind speed of the wind measuring tower at the hub height measuring layer (hub level) of the wind turbine is calculated in an important way
Step 3, setting the surface roughness z according to the terrain and the surface vegetation coverage condition of the area0；
The surface roughness is a general method for those skilled in the art, and will not be described in detail here.
Step 4, the product obtained in step 2 is processed
Initial driving wind speed as a leading simulation (precondition simulation)
(that is to say
Is directly assigned to
) In combination with z set in step 30Performing leading simulation; what is needed isThe atmospheric turbulence mode in the precursor simulation adopts a classical k-epsilon model in an RANS class, wherein k is turbulence energy, and epsilon is the dissipation rate of the turbulence energy;
preamble simulation is a well-known method for those skilled in the art, and will not be described in detail herein, and only the setting of key parameters thereof will be discussed).
Geometric calculation region D of the preamble simulationpThe cushion is a standard hexahedron, the lower cushion surface is a plane, the length and the width are kilometers (usually more than 3km), and the height is 1 to 2 kilometers (the height of an atmospheric boundary layer); periodic Boundary Conditions (PBS) are set at 4 circumferential Boundary surfaces perpendicular to the underlying surface to perform a leading simulation, as shown in fig. 2 (a);
step 5, stopping calculating when the leading simulation result in the step 4 reaches a steady state, and enabling contour lines of flow variables such as wind speed U, turbulence energy k and turbulence energy dissipation rate epsilon in the steady state to change along with height to be in a calculation domain DpInternally sampling, as shown in fig. 2 (b);
step 6, using the contour line of the flow variable (U, k, epsilon) in the step 5 as an inflow boundary condition, and mapping (mapping) the flow variable to the main simulation area DTOn the inflow boundary of the main simulation region DTContinuously adopting an RANS turbulence model for numerical calculation, and after the simulation result reaches a statistical stable state, simulating the wind speed of the representative anemometer tower on the hub height layer
Sampling and recording;
step 7, the method is to carry out the steps of step 6
And described in step 4
And (3) comparison:
(1) if there is a large error between the two
Then returning to step 4, will
Adjusting the value: when in use
Far greater than
Then will be
The initial value is adjusted to be smaller than that in the step 4; when in use
Much less than
Then will be
Is adjusted to be larger than the initial value in the step 4.
(2) If the error between the two is within the acceptable range
The final driving wind speed in step 4 is recorded
The values correspond to the flow variable (U, k, ε) contour features in step 6, and the flow variable (U, k, ε) contour is ultimately mapped to the main simulation region DTSteady-state (steady state) inflow boundary conditions on the inflow boundary, as flow statistical variables under RANS-like turbulent flow modes, as shown in FIG. 2 (c);
step 8, when a Large Eddy Simulation (LES) type atmospheric turbulence mode is required to be adopted for wind field numerical calculation, the driving wind speed finally obtained in the step 4 can be selected
In combination with z set in step 20In the domain D of the pilot simulation calculationpThe 4 circumferential boundary surfaces in (a) adopt periodic boundary conditions to carry out leading simulation, as shown in fig. 3 (a);
step 9, when the leading simulation result in the step 8 reaches a statistical steady state (statistical steady state), selecting a leading simulation calculation domain DpAnd (c) storing the wind speed time-course data of the cross-section area in any cross section perpendicular to the main flow direction, and establishing a time sequence database of unsteady transient wind speeds in a statistical stable state, as shown in fig. 3 (b).
The statistical steady state refers to a non-steady (unsteady) state that the system (time and space) average value of the flow field variable does not change any more in different time periods, and the instantaneous value of the flow field variable still fluctuates around the average value.
The above description is only a preferred guiding embodiment of the present invention, and the protection scope of the present invention is not limited thereto, and any person skilled in the art can easily think of the replacement of the specific model in the two major turbulence modes of the Reynolds Average (RANS) and the Large Eddy Simulation (LES) and other simulation parameters by using the inflow boundary condition generation method as the core within the technical scope of the present invention, and shall be covered by the protection scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
As shown in fig. 4, a second object of the present invention is to provide an inflow condition generating system for numerical simulation of a wind field in a complex terrain, comprising:
the simulated wind speed calculation module is used for sampling a steady-state flow statistical variable contour line obtained by performing leading simulation in a Reynolds Average (RANS) type atmospheric turbulence mode, mapping the contour line to an inflow boundary surface of a main simulation area where a wind field in a complex terrain is located, and performing flow field simulation calculation of the main simulation area;
the comparison and correction module is used for comparing the simulated wind speed at the representative anemometer tower in the main simulation area with the actually measured average wind speed; after the simulated wind speed is compared with the actually measured average wind speed, the driving wind speed parameters of the leading simulation are fed back and corrected, and the steady-state inflow boundary conditions of the main simulation area are generated;
and the boundary condition generating module is used for expanding the leading simulation driving wind speed parameter corresponding to the steady-state inflow boundary condition to a Large Eddy Simulation (LES) type atmospheric turbulence mode, establishing a time sequence database of unsteady transient wind speed, and generating the transient wind speed inflow boundary condition of the main simulation area.
A third object of the present invention is to provide an electronic device, as shown in fig. 5, including a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the steps of generating the inflow conditions for the numerical simulation of the complex terrain wind field when executing the computer program.
A fourth object of the present invention is to provide a computer readable storage medium, which stores a computer program that, when being executed by a processor, performs the step of generating the inflow condition of the complex terrain wind field numerical simulation.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solutions of the present invention and not for limiting the same, and although the present invention is described in detail with reference to the above embodiments, those of ordinary skill in the art should understand that: modifications and equivalents may be made to the embodiments of the invention without departing from the spirit and scope of the invention, which is to be covered by the claims.
Claims (10)
1. An inflow condition generation method for numerical simulation of a wind field in a complex terrain is characterized by comprising the following steps of:
sampling a flow statistical variable contour line under a steady state obtained by performing leading simulation in a Reynolds average atmospheric turbulence mode, mapping the contour line to an inflow boundary surface of a main simulation area where a complex terrain wind field is located, and performing flow field simulation calculation of the main simulation area;
comparing the simulated wind speed at the representative anemometer tower in the main simulation area with the actually measured average wind speed; according to the comparison result of the simulated wind speed and the actually measured average wind speed, the driving wind speed parameters of the leading simulation are fed back and corrected, and then the steady-state inflow boundary condition of the main simulation area is generated;
and expanding the leading simulation driving wind speed parameters corresponding to the steady-state inflow boundary conditions to a large vortex simulation type atmospheric turbulence mode, establishing a time sequence database of unsteady transient wind speeds, and generating transient wind speed inflow boundary conditions of a main simulation area according to the time sequence database.
2. The method of claim 1,
after sampling the flow statistical variable contour line under the steady state, mapping the flow statistical variable contour line to an inflow boundary surface of a main simulation area where a complex terrain wind field is located, and performing flow field simulation calculation on the main simulation area; the method specifically comprises the following steps:
determining a hexahedron main simulation area for numerical simulation of atmospheric flow in a complex terrain area where a wind field is located, and taking the complex terrain of a target area as a ground surface boundary surface;
counting the average wind speed and the standard deviation of the wind speed at each height measurement layer of a representative wind measuring tower in a wind power plant area within a target time period; setting the roughness of the earth surface according to the terrain of the area and the vegetation coverage condition of the earth surface; taking the average wind speed as the initial driving wind speed of the leading simulation, and combining the surface roughness to perform leading simulation calculation; stopping calculation when the leading simulation result reaches a stable state, sampling the contour line of the flow variable changing along with the height in the main calculation simulation area,
and mapping the contour line of the flow variable as an inflow boundary condition to an inflow boundary of a main simulation area, continuously adopting a Reynolds average turbulence model to perform numerical calculation in the main simulation area, and sampling and recording the simulated wind speed of the representative anemometer tower at the hub height layer after the simulation result reaches a statistical stable state.
3. The method of claim 2,
the geometric area of the preamble simulation is a standard hexahedron, the underlying surface is a plane, the length and the width are kilometers, and the height is the height of an atmospheric boundary layer; periodic boundary conditions were set at 4 circumferential boundary surfaces thereof perpendicular to the underlying surface.
4. The method of claim 1,
the step of comparing the simulated wind speed at the representative anemometer tower in the main simulation area with the actually measured average wind speed, and feeding back and correcting the drive wind speed parameter of the leading simulation according to the comparison result of the simulated wind speed and the actually measured average wind speed to generate the steady-state inflow boundary condition of the main simulation area specifically comprises the following steps:
if the difference value between the simulated wind speed and the actually measured average wind speed is larger than the threshold value, returning to adjust the initial driving wind speed value of the leading simulation:
when the simulated wind speed is larger than the actually measured average wind speed, the initial value of the initial driving wind speed is reduced; when the simulated wind speed is smaller than the actually measured average wind speed, the initial driving wind speed is increased;
if the difference value between the simulated wind speed and the actually measured average wind speed is smaller than the threshold value, recording the flow variable contour line characteristics corresponding to the final leading simulated driving wind speed value, and finally mapping the flow variable contour line to the inflow boundary of the main simulation area to be used as the steady-state inflow boundary condition of the flow statistical variable in the Reynolds average turbulence-like mode.
5. The method of claim 1,
the method for establishing the time series database of the unsteady transient wind speed by expanding the leading simulation driving wind speed parameter corresponding to the steady-state inflow boundary condition to a large vortex simulation type atmospheric turbulence mode specifically comprises the following steps:
adopting a large vortex simulation type atmospheric turbulence mode, selecting the finally obtained driving wind speed under a Reynolds average type turbulence mode, and performing leading simulation by combining set surface roughness and adopting periodic boundary conditions on 4 circumferential boundary surfaces in a leading simulation calculation main simulation area;
and when the leading simulation result reaches a statistical stable state, selecting any section perpendicular to the main flow direction in the main simulation area of the leading simulation calculation area, storing the wind speed time course data of the section area, and establishing a time sequence database of the unsteady transient wind speed in the statistical stable state.
6. The method of claim 5,
the statistical steady state refers to an unsteady state that the average value of the flow field variable in different time periods does not change any more, and the instantaneous value of the flow field variable still fluctuates around the average value.
7. The method of claim 1,
the generating of the transient wind speed inflow boundary condition of the main simulation area comprises:
and extracting and mapping the section wind speed time-course data in the unsteady transient wind speed time-series database on the inflow boundary surface of the main simulation area time-step by time-step based on the actual elevation of each space point of the inflow boundary surface in the main simulation area from the complex earth surface to form the transient wind speed inflow boundary condition in the large vortex simulation turbulence mode.
8. An inflow condition generating system for numerical simulation of a wind field in complex terrain, comprising:
the simulated wind speed calculation module is used for sampling a flow statistical variable contour line under a steady state obtained by performing leading simulation in a Reynolds average atmospheric turbulence mode, mapping the flow statistical variable contour line to an inflow boundary surface of a main simulation area where a wind field with complex terrain is located, and performing flow field simulation calculation of the main simulation area;
the comparison and correction module is used for comparing the simulated wind speed at the representative anemometer tower in the main simulation area with the actually measured average wind speed; according to the comparison result of the simulated wind speed and the actually measured average wind speed, the driving wind speed parameters of the leading simulation are fed back and corrected, and then the steady-state inflow boundary condition of the main simulation area is generated;
and the boundary condition generating module is used for expanding the leading simulation driving wind speed parameter corresponding to the steady-state inflow boundary condition to a large vortex simulation type atmospheric turbulence mode, establishing a time sequence database of unsteady transient wind speed, and generating the transient wind speed inflow boundary condition of the main simulation area according to the time sequence database.
9. An electronic device comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing the steps of generating a numerically simulated inflow condition of a complex terrain wind field according to any of claims 1 to 7 when executing the computer program.
10. A computer readable storage medium storing a computer program which, when executed by a processor, performs the steps of generating an inflow condition for a numerical simulation of a complex terrain wind field according to any of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110998569.1A CN113657055B (en) | 2021-08-27 | 2021-08-27 | Inflow condition generation method and system for complex terrain wind field numerical simulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110998569.1A CN113657055B (en) | 2021-08-27 | 2021-08-27 | Inflow condition generation method and system for complex terrain wind field numerical simulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113657055A true CN113657055A (en) | 2021-11-16 |
| CN113657055B CN113657055B (en) | 2024-05-07 |
Family
ID=78482310
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110998569.1A Active CN113657055B (en) | 2021-08-27 | 2021-08-27 | Inflow condition generation method and system for complex terrain wind field numerical simulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113657055B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114595643A (en) * | 2022-03-04 | 2022-06-07 | 重庆大学 | High-resolution mountain wind field measurement method coupling multi-point anemometer and micro-scale steady-state simulation |
| CN116070414A (en) * | 2022-12-12 | 2023-05-05 | 中广核风电有限公司 | Wind resource simulation method and device based on surface relief complexity |
| CN116341410A (en) * | 2023-02-07 | 2023-06-27 | 成都流体动力创新中心 | Typical terrain wind field simulation method and system considering different stabilities |
Citations (6)
| 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 (en) * | 2017-10-09 | 2018-02-23 | 中国民航大学 | Anisotropy In The Atmospheric Turbulent Field method for numerical simulation |
| US20200117842A1 (en) * | 2016-12-21 | 2020-04-16 | Northwestern University | Systems and methods for computational simulation of self-propelling vehicles for aerodynamic design |
| CN112163381A (en) * | 2020-09-25 | 2021-01-01 | 西安热工研究院有限公司 | Lateral boundary condition setting method suitable for complex terrain wind field flow numerical simulation |
| US20210124861A1 (en) * | 2019-10-28 | 2021-04-29 | Vorcat, Inc. | System and method for simulating turbulence |
| CN112949227A (en) * | 2021-03-30 | 2021-06-11 | 中国华能集团清洁能源技术研究院有限公司 | Method, system, equipment and storage medium suitable for determining turbulence intensity of wind field in complex terrain |
-
2021
- 2021-08-27 CN CN202110998569.1A patent/CN113657055B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2012105973A1 (en) * | 2011-02-02 | 2012-08-09 | Michigan Aerospace Corporation | Atmospheric measurement system and method |
| US20200117842A1 (en) * | 2016-12-21 | 2020-04-16 | Northwestern University | Systems and methods for computational simulation of self-propelling vehicles for aerodynamic design |
| CN107729638A (en) * | 2017-10-09 | 2018-02-23 | 中国民航大学 | Anisotropy In The Atmospheric Turbulent Field method for numerical simulation |
| US20210124861A1 (en) * | 2019-10-28 | 2021-04-29 | Vorcat, Inc. | System and method for simulating turbulence |
| CN112163381A (en) * | 2020-09-25 | 2021-01-01 | 西安热工研究院有限公司 | Lateral boundary condition setting method suitable for complex terrain wind field flow numerical simulation |
| CN112949227A (en) * | 2021-03-30 | 2021-06-11 | 中国华能集团清洁能源技术研究院有限公司 | Method, system, equipment and storage medium suitable for determining turbulence intensity of wind field in complex terrain |
Non-Patent Citations (1)
| Title |
|---|
| 梁思超;张晓东;康顺;康雅兰;赵永锋;: "基于数值模拟的复杂地形风场风资源评估方法", 空气动力学学报, no. 03, 15 June 2012 (2012-06-15) * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114595643A (en) * | 2022-03-04 | 2022-06-07 | 重庆大学 | High-resolution mountain wind field measurement method coupling multi-point anemometer and micro-scale steady-state simulation |
| CN114595643B (en) * | 2022-03-04 | 2024-06-07 | 重庆大学 | High-resolution mountain wind field measurement method coupling multi-point anemometer and microscale steady-state simulation |
| CN116070414A (en) * | 2022-12-12 | 2023-05-05 | 中广核风电有限公司 | Wind resource simulation method and device based on surface relief complexity |
| CN116070414B (en) * | 2022-12-12 | 2023-10-17 | 中广核风电有限公司 | Wind resource simulation method and device based on surface relief complexity |
| CN116341410A (en) * | 2023-02-07 | 2023-06-27 | 成都流体动力创新中心 | Typical terrain wind field simulation method and system considering different stabilities |
| CN116341410B (en) * | 2023-02-07 | 2024-01-30 | 成都流体动力创新中心 | Typical terrain wind field simulation method and system considering different stabilities |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113657055B (en) | 2024-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2017352549B2 (en) | Method and device for calculating power generation of wind farm | |
| CN113657055A (en) | Inflow condition generation method and system for numerical simulation of wind field in complex terrain | |
| Ti et al. | Artificial Neural Networks based wake model for power prediction of wind farm | |
| CN102663251B (en) | Physical prediction method for wind power station power based on computational fluid mechanics model | |
| CN104699936A (en) | Sector management method based on CFD short-term wind speed forecasting wind power plant | |
| US20190370418A1 (en) | Method and apparatus for arranging wind turbines based on rapid accessment fluid model and wake model | |
| CN107194097A (en) | Analysis method based on wind power plant pneumatic analog and wind speed and direction data | |
| CN108763825B (en) | Numerical simulation method for simulating wind field of complex terrain | |
| Syed et al. | Partial repowering analysis of a wind farm by turbine hub height variation to mitigate neighboring wind farm wake interference using mesoscale simulations | |
| CN112163381B (en) | Lateral boundary condition setting method suitable for complex terrain wind field flow numerical simulation | |
| CN108717593A (en) | A kind of microcosmic structure generated energy appraisal procedure based on wind wheel face equivalent wind speed | |
| CN108509718B (en) | Far-field wake two-dimensional analytic model based on mass conservation | |
| Diaz et al. | Full wind rose wind farm simulation including wake and terrain effects for energy yield assessment | |
| Branlard et al. | A digital twin solution for floating offshore wind turbines validated using a full-scale prototype | |
| CN114707437A (en) | Method for establishing three-dimensional full wake flow model of horizontal axis wind turbine | |
| Zhang et al. | Discussion on the spatial-temporal inhomogeneity characteristic of horizontal-axis wind turbine's wake and improvement of four typical wake models | |
| Lin et al. | Large-eddy simulation of a wind-turbine array subjected to active yaw control | |
| Wang et al. | Turbulence intensity identification and load reduction of wind turbine under extreme turbulence | |
| CN114186407A (en) | Wind power plant wake velocity field calculation method and system with parameters capable of being adjusted in self-adaptive mode | |
| CN115329690A (en) | Wake flow simulation method and system for wind power plant, control device and readable storage medium | |
| CN114154325A (en) | Method and system for evaluating steady and steady mixed wind energy resources | |
| Hwang et al. | Wind farm layout optimization using multidisciplinary model | |
| CN117313399B (en) | Method for establishing and applying three-dimensional anisotropic ultra-high-speed full wake model of horizontal axis wind turbine suitable for complex terrain | |
| Liu et al. | Data-driven parameter calibration in wake models | |
| 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 |