Flow field calculation method and device for large-scale wind power plant

Abstract

The invention relates to the technical field of flow field calculation, in particular to a flow field calculation method and device for a large-scale wind power plant. The method comprises the following steps: carrying out grid division on a target area of a wind power plant to obtain a grid target area; based on a preset segmentation direction, dividing the gridded target area into a plurality of sub-areas, each adjacent sub-area including an overlapping area; for each sub-region, an overlapping region between the sub-region and the adjacent sub-region is divided into a first shared region and a second shared region along the segmentation direction, the two shared regions are respectively subjected to flow field calculation in different sub-regions, and flow field data in the shared regions in the adjacent sub-regions are exchanged in the calculation process; obtaining a flow field calculation result of each sub-region; and superposing the flow field calculation results of the sub-regions to obtain the flow field parameters of the target region. The method is high in calculation efficiency and accurate in result.

Description

Technical Field

[0001] The present invention relates to the technical field of flow field calculation, and in particular to a flow field calculation method and device for a large-scale wind farm. Background Art

[0002] With the rapid development of wind power generation, the scale of wind turbines is increasing, and the terrain range of wind farms is becoming more and more complex. Before developing a wind farm, it is necessary to simulate the flow field of the entire field to evaluate the wind resources. However, for large wind farms, since their terrain range is dozens of times larger than that of conventional small wind farms, the overall calculation amount increases by thousands of times, and the amount of computer hardware resources required is very huge. Therefore, existing ordinary workstations cannot complete the flow field evaluation of wind farms.

[0003] Based on this, there is an urgent need for a flow field calculation method and device for large-scale wind farms to solve the above problems. Summary of the invention

[0004] The present invention provides a flow field calculation method and device for a large-scale wind farm, which has high calculation efficiency and accurate results. The technical solution is as follows: In a first aspect, an embodiment of the present invention provides a flow field calculation method for a large-scale wind farm, the method comprising: Gridding the target area of ​​the wind farm to obtain a gridded target area; Based on a preset segmentation direction, the gridded target area is divided into a plurality of sub-areas, and each adjacent sub-area includes an overlapping area; For each of the sub-regions, the overlapping area between the sub-region and its adjacent sub-region is divided into a first shared region and a second shared region along the segmentation direction, the flow field of the first shared region is calculated by a computing device corresponding to the sub-region, and the flow field of the second shared region is calculated by a computing device corresponding to the adjacent sub-region; Different computing devices are used to calculate the flow field of each sub-area, and during the calculation process, the flow field data in the shared area of ​​adjacent sub-areas are exchanged to obtain the flow field calculation results of each sub-area; The flow field calculation results of each sub-area are superimposed to obtain the flow field parameters of the target area.

[0005] In a second aspect, an embodiment of the present invention further provides a flow field calculation device for a large-scale wind farm, the device comprising: A first division unit is used to divide the target area of ​​the wind farm into grids to obtain a gridded target area; A second division unit, configured to divide the gridded target area into a plurality of sub-areas based on a preset division direction, wherein each adjacent sub-area includes an overlapping area; A third division unit is used to divide, for each of the sub-regions, an overlapping area between the sub-region and its adjacent sub-region into a first shared region and a second shared region along the segmentation direction, wherein the flow field of the first shared region is calculated by a computing device corresponding to the sub-region, and the flow field of the second shared region is calculated by a computing device corresponding to the adjacent sub-region; A calculation unit, used to respectively use different calculation devices to calculate the flow field of each sub-area, and exchange the flow field data in the shared area of ​​adjacent sub-areas during the calculation process to obtain the flow field calculation results of each sub-area; The superposition unit is used to superimpose the flow field calculation results of each sub-area to obtain the flow field parameters of the target area.

[0006] In a third aspect, an embodiment of the present invention further provides an electronic device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the method described in any embodiment of this specification is implemented.

[0007] In a fourth aspect, an embodiment of the present invention further provides a computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, enables the computer to execute the method described in any embodiment of this specification.

[0008] In a fifth aspect, an embodiment of the present invention further provides a computer program product, including a computer program, which implements the steps of the method described above when executed by a processor.

[0009] The embodiment of the present invention provides a method and device for calculating the flow field of a large-scale wind farm. On the one hand, by dividing a large-scale wind farm into multiple sub-areas, each sub-area can be calculated separately in parallel, thereby improving the calculation efficiency. On the other hand, when calculating the flow field of each sub-area, by exchanging the flow field data in the first shared area and the second shared area between adjacent sub-areas, the calculation results between adjacent sub-areas can be coupled to each other, preventing the calculation results between the sub-areas from diverging, and improving the accuracy of the calculation results. It can be seen that the present application has high calculation efficiency and accurate results. BRIEF DESCRIPTION OF THE DRAWINGS

[0010] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0011] Figure 1This is a flow chart of a flow field calculation method for a large-scale wind farm provided by an embodiment of the present invention; Figure 2 This is a structural diagram of a flow field calculation device for a large-scale wind farm provided by an embodiment of the present invention; Figure 3 is a hardware architecture diagram of a computer device provided by an embodiment of the present invention; Figure 4 is a schematic diagram of dividing overlapping areas between adjacent sub-areas provided by an embodiment of the present invention; Figure 5 It is a schematic diagram of dividing overlapping areas between adjacent sub-areas provided by another embodiment of the present invention. DETAILED DESCRIPTION

[0012] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of the present invention.

[0013] The specific implementation of the above concept is described below.

[0014] Please refer to Figure 1 , an embodiment of the present invention provides a flow field calculation method for a large-scale wind farm, the method comprising: Step 100, gridding the target area of ​​the wind farm to obtain a gridded target area; Step 102, based on a preset segmentation direction, the gridded target area is divided into a plurality of sub-areas, and each adjacent sub-area includes an overlapping area; Step 104, for each sub-region, the overlapping area between the sub-region and its adjacent sub-region is divided into a first shared region and a second shared region along the segmentation direction, the flow field of the first shared region is calculated by the computing device corresponding to the sub-region, and the flow field of the second shared region is calculated by the computing device corresponding to the adjacent sub-region; Step 106, using different computing devices to perform flow field calculation on each sub-area, and in the calculation process, exchanging flow field data in a shared area between adjacent sub-areas to obtain flow field calculation results for each sub-area; Step 108, superimposing the flow field calculation results of each sub-area to obtain the flow field parameters of the target area.

[0015] In this embodiment, on the one hand, by dividing a large-scale wind farm into multiple sub-areas, each sub-area can be calculated separately in parallel, thereby improving the calculation efficiency. On the other hand, when calculating the flow field of each sub-area, by exchanging the flow field data in the first shared area and the second shared area between adjacent sub-areas, the calculation results between adjacent sub-areas can be coupled to each other, preventing the calculation results between the sub-areas from diverging, and improving the accuracy of the calculation results. It can be seen that the present application has high calculation efficiency and accurate results.

[0016] Described below Figure 1 How the various steps are performed.

[0017] First, with respect to step 100 , the target area of ​​the wind farm is gridded to obtain a gridded target area.

[0018] In this step, the target area is determined according to actual needs. In addition, the smaller the grid, the higher the calculation accuracy, but the amount of calculation also increases; conversely, the calculation accuracy decreases and the amount of calculation also decreases. The user can determine it independently according to needs, and this application does not make specific restrictions.

[0019] With respect to step 102 , the gridded target area is divided into a plurality of sub-areas based on a preset segmentation direction, and each adjacent sub-area includes an overlapping area.

[0020] In this step, the segmentation direction can be the wind direction or the terrain trend direction. The wind direction can be the main wind direction throughout the year. If the wind is northwest most of the year, the northwest to southeast direction will be used as the segmentation direction. This step adopts the one-way partitioning method, which can simplify the model and reduce the difficulty of calculation. The size of the sub-area and the overlapping area is determined according to the user's needs, and this application does not make specific restrictions.

[0021] In addition, the number of grids between each sub-area is as balanced as possible to ensure that the computing resources required by each sub-area are relatively balanced.

[0022] For step 104, for each sub-region, the overlapping area between the sub-region and its adjacent sub-region is divided into a first shared region and a second shared region along the segmentation direction, and the flow field of the first shared region is calculated by the computing device corresponding to the sub-region, and the second shared region is calculated by the computing device corresponding to its adjacent sub-region.

[0023] In this step, if Figure 4 and Figure 5 As shown, taking a two-dimensional grid as an example, the relationship between adjacent sub-areas and the method of dividing overlapping areas are illustrated.

[0024] In some embodiments, the grid dimension of the overlapping area is an even number; the number of grids in the first sharing area is equal to the number of grids in the second sharing area, and the first sharing area and the second sharing area are respectively allocated to the sub-area closest to them for calculation.

[0025] In this embodiment, if Figure 4 and Figure 5 As shown, both sub-areas include the first shared area and the second shared area. However, during calculation, the first shared area is calculated by sub-area 1, and the second shared area is calculated by sub-area 2. This can avoid repeated calculations of the two areas and cause data confusion. In addition, by setting the grid dimension of the overlapping area to an even number, it is convenient to divide the overlapping area into equal parts. The two shared areas after equal division are respectively assigned to different sub-areas for calculation, and the calculation results are exchanged to ensure the convergence and accuracy of the calculation results.

[0026] For step 106, different computing devices are used to perform flow field calculation on each sub-region, and flow field data in the shared region of adjacent sub-regions are exchanged during the calculation process to obtain flow field calculation results for each sub-region, including: Perform iterative calculations on the grid in each sub-region separately, and execute the following for each iteration: For the grids other than the shared grid in the sub-region, the flow field calculation result of the previous iteration of the grid is used as the initial value of the current iteration, and the current flow field calculation is performed based on the initial value; For the first shared area in the sub-area, its initial value in the current iteration is determined based on the flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area; and so on, until the flow field in the sub-area converges, and the flow field calculation result of the sub-area is obtained.

[0027] It should be noted that if the overlapped area is too large, the amount of calculation will increase and the accuracy of the calculation results will be little improved. Therefore, the dimension of the overlapped area along the segmentation direction is preferably two-dimensional or four-dimensional. The following is a specific description of the data exchange method of the shared area under the two dimensions.

[0028] The first case, such as Figure 4 As shown, the grid dimensions of the first shared area and the second shared area along the segmentation direction are both one-dimensional; For the first shared area in the sub-area, determining its initial value in the current iteration based on the flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area includes: The flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area is used as the initial value of the first shared area in the current iteration; Assuming the x direction is the segmentation direction, the data exchange formula for adjacent sub-regions is: In the formula, Indicates the grid of the first shared area at the k+1th iteration of the flow field parameter s The initial value of Indicates the x-direction i A grid; Indicates that when the flow field parameter s converges at the kth iteration, the grid of the second shared area The flow field calculation results of Indicates the x-direction i +1 grid; i =1~X, y=1~Y, z=1~Z, where X, Y, and Z are the maximum number of grids in the x, y, and z directions respectively.

[0029] The second case, such as Figure 5 As shown, the grid dimensions of the first shared area and the second shared area along the segmentation direction are both two-dimensional; For the first shared area in the sub-area, determining its initial value in the current iteration based on the flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area includes: The calculation result of the first dimension grid of the second shared area in the sub-area in the previous iteration in the adjacent sub-area is used as the initial value of the first dimension grid of the first shared area in the current iteration; The calculation result of the second dimension grid of the second shared area in the sub-area in the previous iteration in the adjacent sub-area is used as the initial value of the second dimension grid of the first shared area in the current iteration; Assuming the x direction is the segmentation direction, the data exchange formula for adjacent sub-regions is: In the formula, The first dimension grid of the first shared area representing the flow field parameter s at the k+1th iteration The initial value of Indicates the x-direction i A grid; Indicates the first dimension grid of the second shared area when the flow field parameter s converges at the kth iteration The flow field calculation results of Indicates the x-direction i +2 grids; The second dimensional grid representing the first shared area of ​​the flow field parameter s at the k+1th iteration The initial value of Indicates the x-direction i +1 grid; Indicates the second dimension grid of the second shared area when the flow field parameter s converges at the kth iteration The flow field calculation results of Indicates the x-direction i +3 grids.

[0030] Through the above calculation method, before each iterative calculation of each sub-region, the flow field calculation results in the shared area are exchanged, so as to ensure the convergence and accuracy of the calculation results of each sub-region.

[0031] In some embodiments, the flow field parameters include: wind energy density and turbulence intensity. Each flow field parameter is calculated by CFD, wherein the calculation process of turbulence intensity is a common method in the art and will not be repeated here. The wind energy density of each area is calculated by the following formula: in, In the formula, For sub-region n The wind energy density, n =1~N, N is the total number of sub-areas; is the air density; is the free wind speed outside the wind farm; For sub-region n The speed loss of the mth wind turbine in the middle, m=1~M, M is the sub-area n The total number of blowers; hour, ;when hour, ; is the impeller diameter of the mth fan; and Sub-areas n The flow distance and span distance between adjacent fans; For sub-region n The wake diffusion coefficient is According to sub-area n Determined by the terrain characteristics and wind turbine distribution; is the thrust coefficient of the mth wind turbine, which is determined according to the blade shape and wind speed of the wind turbine; G is the boundary layer height of the atmosphere above the sub-region; For sub-region n Average roughness of the surface; is the hub height of the mth wind turbine; is the eddy viscosity coefficient.

[0032] This embodiment adopts the above calculation method, comprehensively considering the influence of factors such as the atmospheric height, surface parameters, topography, wind turbine parameters (such as blade shape, diameter and hub height, etc.) and wind turbine distribution (the mutual influence of the installation position and distance of each wind turbine) in the area where the wind farm is located, so as to improve the calculation accuracy of the wind energy density in the sub-area and facilitate the flow field analysis of the wind farm.

[0033] Finally, for step 108, after the flow field of each sub-region is calculated, for adjacent sub-regions, since there are overlapping areas in the two sub-regions, only the overlapping area in one of the sub-regions is retained, and the overlapping area in the other sub-region is removed, and then splicing is performed. Similarly, the flow field parameters of the target area are obtained.

[0034] like Figure 2 , Figure 3 As shown, an embodiment of the present invention provides a flow field calculation device for a large-scale wind farm. The device embodiment can be implemented by software, or by hardware or a combination of software and hardware. From the hardware level, Figure 2 As shown in FIG. 1 , a hardware architecture diagram of a computing device where a flow field computing device for a large-scale wind farm is provided in an embodiment of the present invention is provided. Figure 2 In addition to the processor, memory, network interface, and non-volatile memory shown in the figure, the computing device in which the device is located in the embodiment may also generally include other hardware, such as a forwarding chip responsible for processing messages, etc. Taking software implementation as an example, Figure 3 As shown, as a device in a logical sense, the CPU of the computing device in which it is located reads the corresponding computer program in the non-volatile memory into the internal memory and runs it.

[0035] Please refer to Figure 3 The embodiment of the present invention provides a flow field calculation device for a large-scale wind farm, the device comprising: The first division unit 300 is used to divide the target area of ​​the wind farm into grids to obtain a gridded target area; The second division unit 302 is used to divide the gridded target area into a plurality of sub-areas based on a preset segmentation direction, and each adjacent sub-area includes an overlapped area; The third division unit 304 is used to divide, for each sub-region, an overlapping area between the sub-region and its adjacent sub-region into a first shared region and a second shared region along the segmentation direction, wherein the flow field of the first shared region is calculated by a computing device corresponding to the sub-region, and the flow field of the second shared region is calculated by a computing device corresponding to the adjacent sub-region; The calculation unit 306 is used to use different calculation devices to perform flow field calculation on each sub-area, and exchange the flow field data in the shared area of ​​adjacent sub-areas during the calculation process to obtain the flow field calculation results of each sub-area; The superposition unit 308 is used to superimpose the flow field calculation results of each sub-area to obtain the flow field parameters of the target area.

[0036] In some embodiments, the segmentation direction is: wind direction or terrain trend direction.

[0037] In some embodiments, along the segmentation direction, the grid dimension of the overlapping area is an even number; the number of grids in the first shared area is equal to the number of grids in the second shared area, and the first shared area and the second shared area are respectively allocated to the sub-area closest to them for calculation.

[0038] In some implementations, the computing unit 306 is configured to perform the following operations: Perform iterative calculations on the grid in each sub-region separately, and execute the following for each iteration: For the grids other than the shared grid in the sub-region, the flow field calculation result of the previous iteration of the grid is used as the initial value of the current iteration, and the current flow field calculation is performed based on the initial value; For the first shared area in the sub-area, its initial value in the current iteration is determined based on the flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area; and so on, until the flow field in the sub-area converges, and the flow field calculation result of the sub-area is obtained.

[0039] In some embodiments, the grid dimensions of the first shared area and the second shared area along the segmentation direction are both one-dimensional; For the first shared area in the sub-area, determining its initial value in the current iteration based on the flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area includes: The flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area is used as the initial value of the first shared area in the current iteration; Assuming the x direction is the segmentation direction, the data exchange formula for adjacent sub-regions is: In the formula, Indicates the grid of the first shared area at the k+1th iteration of the flow field parameter s The initial value of Indicates the x-direction i A grid; Indicates that when the flow field parameter s converges at the kth iteration, the grid of the second shared area The flow field calculation results of Indicates the x-direction i +1 grid; i =1~X, y=1~Y, z=1~Z, where X, Y, and Z are the maximum number of grids in the x, y, and z directions respectively.

[0040] In some embodiments, the grid dimensions of the first shared area and the second shared area along the segmentation direction are both two-dimensional; For the first shared area in the sub-area, determining its initial value in the current iteration based on the flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area includes: The calculation result of the first dimension grid of the second shared area in the sub-area in the previous iteration in the adjacent sub-area is used as the initial value of the first dimension grid of the first shared area in the current iteration; The calculation result of the second dimension grid of the second shared area in the sub-area in the previous iteration in the adjacent sub-area is used as the initial value of the second dimension grid of the first shared area in the current iteration; Assuming the x direction is the segmentation direction, the data exchange formula for adjacent sub-regions is: In the formula, The first dimension grid of the first shared area representing the flow field parameter s at the k+1th iteration The initial value of Indicates the x-direction i A grid; Indicates the first dimension grid of the second shared area when the flow field parameter s converges at the kth iteration The flow field calculation results of Indicates the x-direction i +2 grids; The second dimensional grid representing the first shared area of ​​the flow field parameter s at the k+1th iteration The initial value of Indicates the x-direction i +1 grid; Indicates the second dimension grid of the second shared area when the flow field parameter s converges at the kth iteration The flow field calculation results of Indicates the x-direction i +3 grids.

[0041] In some embodiments, the flow field parameters include: wind energy density and turbulence intensity.

[0042] It should be noted that the flow field calculation device for a large-scale wind farm provided in the above embodiment is only illustrated by the division of the above functional modules. In actual applications, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the device is divided into different functional modules to complete all or part of the functions described above. In addition, the flow field calculation device for a large-scale wind farm provided in the above embodiment and the flow field calculation method embodiment of a large-scale wind farm belong to the same concept. The specific implementation process is detailed in the method embodiment and will not be repeated here.

[0043] The embodiment of the present application also provides a computer device, please refer to Figure 3 The computer device includes a processor and a memory, in which at least one instruction, at least one program, a code set or an instruction set is stored, and the at least one instruction, at least one program, a code set or an instruction set is loaded and executed by the processor to implement the flow field calculation method for a large-scale wind farm provided by the above-mentioned method embodiments.

[0044] An embodiment of the present application also provides a computer-readable storage medium, on which is stored at least one instruction, at least one program, code set or instruction set, and the at least one instruction, at least one program, code set or instruction set is loaded and executed by a processor to implement the flow field calculation method for a large-scale wind farm provided by the above-mentioned method embodiments.

[0045] An embodiment of the present application also provides a computer program product, which includes a computer program. A processor of a computer device reads the computer program from a computer-readable storage medium, and the processor executes the computer program, so that the computer device executes the flow field calculation method for a large-scale wind farm described in any of the above embodiments.

[0046] For the convenience of description, the above system or device is described by dividing it into various modules or units according to its functions. Of course, when implementing the present application, the functions of each unit can be implemented in the same or multiple software and / or hardware.

[0047] It can be known from the description of the above implementation methods that those skilled in the art can clearly understand that the present application can be implemented by means of software plus a necessary general hardware platform. Based on such an understanding, the technical solution of the present application can be essentially or partly contributed to the prior art in the form of a software product, which can be stored in a storage medium such as ROM / RAM, a magnetic disk, an optical disk, etc., and includes several instructions for enabling a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in the various embodiments of the present application or certain parts of the embodiments.

[0048] Finally, it should be noted that, in this article, relational terms such as first, second, third and fourth are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, the elements defined by the statement "comprise a ..." do not exclude the presence of other identical elements in the process, method, article or device including the elements.

[0049] The above is only a preferred implementation of the present application. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principles of the present application. These improvements and modifications should also be regarded as the scope of protection of the present application.

Claims

1. A flow field calculation method for a large-scale wind farm, characterized in that: The method comprises: Gridding the target area of ​​the wind farm to obtain a gridded target area; Based on a preset segmentation direction, the gridded target area is divided into a plurality of sub-areas, and each adjacent sub-area includes an overlapping area; For each of the sub-regions, the overlapping area between the sub-region and its adjacent sub-region is divided into a first shared region and a second shared region along the segmentation direction, the flow field of the first shared region is calculated by a computing device corresponding to the sub-region, and the flow field of the second shared region is calculated by a computing device corresponding to the adjacent sub-region; Different computing devices are used to calculate the flow field of each sub-area, and during the calculation process, the flow field data in the shared area of ​​adjacent sub-areas are exchanged to obtain the flow field calculation results of each sub-area; The flow field calculation results of each sub-area are superimposed to obtain the flow field parameters of the target area.

2. The method according to claim 1, characterized in that The segmentation direction is: wind direction or terrain trend direction.

3. The method according to claim 1, characterized in that Along the segmentation direction, the grid dimension of the overlapping area is an even number; the number of grids in the first shared area is equal to the number of grids in the second shared area, and the first shared area and the second shared area are respectively allocated to the sub-area closest to them for calculation.

4. The method according to claim 3, characterized in that The flow field calculation is performed on each sub-area using different computing devices, and the flow field data in the shared area of ​​adjacent sub-areas are exchanged during the calculation process to obtain the flow field calculation results of each sub-area, including: Perform iterative calculations on the grid in each sub-region separately, and execute the following for each iteration: For the grids other than the shared grid in the sub-region, the flow field calculation result of the previous iteration of the grid is used as the initial value of the current iteration, and the current flow field calculation is performed based on the initial value; For the first shared area in the sub-area, its initial value in the current iteration is determined based on the flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area; and so on, until the flow field in the sub-area converges, and the flow field calculation result of the sub-area is obtained.

5. The method according to claim 4, characterized in that The grid dimensions of the first shared area and the second shared area along the segmentation direction are both one-dimensional; The determining of the initial value of the first shared area in the sub-area in the current iteration based on the flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area includes: The flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area is used as the initial value of the first shared area in the current iteration; Assuming the x direction is the segmentation direction, the data exchange formula for adjacent sub-regions is: In the formula, Indicates the grid of the first shared area at the k+1th iteration of the flow field parameter s The initial value of Indicates the x-direction i A grid; Indicates that when the flow field parameter s converges at the kth iteration, the grid of the second shared area The flow field calculation results of Indicates the x-direction i +1 grid; i =1~X, y=1~Y, z=1~Z, where X, Y, and Z are the maximum number of grids in the x, y, and z directions respectively.

6. The method according to claim 4, characterized in that The grid dimensions of the first shared area and the second shared area along the segmentation direction are both two-dimensional; The determining of the initial value of the first shared area in the sub-area in the current iteration based on the flow field calculation result of the second shared area in the sub-area in the previous iteration in the adjacent sub-area includes: The calculation result of the first dimension grid of the second shared area in the sub-area in the previous iteration in the adjacent sub-area is used as the initial value of the first dimension grid of the first shared area in the current iteration; The calculation result of the second dimension grid of the second shared area in the sub-area in the previous iteration in the adjacent sub-area is used as the initial value of the second dimension grid of the first shared area in the current iteration; Assuming the x direction is the segmentation direction, the data exchange formula for adjacent sub-regions is: In the formula, The first dimension grid of the first shared area representing the flow field parameter s at the k+1th iteration The initial value of Indicates the x-direction i A grid; Indicates the first dimension grid of the second shared area when the flow field parameter s converges at the kth iteration The flow field calculation results of Indicates the x-direction i +2 grids; The second dimension grid representing the first shared area of ​​the flow field parameter s at the k+1th iteration The initial value of Indicates the x-direction i +1 grid; Indicates the second dimension grid of the second shared area when the flow field parameter s converges at the kth iteration The flow field calculation results of Indicates the x-direction i +3 grids.

7. The method according to claim 6, characterized in that The flow field parameters include: wind energy density and turbulence intensity.

8. A flow field calculation device for a large-scale wind farm, characterized in that: The device comprises: A first division unit is used to divide the target area of ​​the wind farm into grids to obtain a gridded target area; A second division unit, configured to divide the gridded target area into a plurality of sub-areas based on a preset division direction, wherein each adjacent sub-area includes an overlapping area; A third division unit is used to divide, for each of the sub-regions, an overlapping area between the sub-region and its adjacent sub-region into a first shared region and a second shared region along the segmentation direction, wherein the flow field of the first shared region is calculated by a computing device corresponding to the sub-region, and the flow field of the second shared region is calculated by a computing device corresponding to the adjacent sub-region; A calculation unit, used to respectively use different calculation devices to calculate the flow field of each sub-area, and exchange the flow field data in the shared area of ​​adjacent sub-areas during the calculation process to obtain the flow field calculation results of each sub-area; The superposition unit is used to superimpose the flow field calculation results of each sub-area to obtain the flow field parameters of the target area.

9. A computing device, comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium having a computer program stored thereon, which, when executed in a computer, causes the computer to execute the method according to any one of claims 1 to 7.

Images