Wind turbine airfoil surface drag reduction rib size optimization design method and device

Abstract

The embodiment of the invention provides a wind turbine airfoil surface drag reduction rib size optimization design method and device. The wind turbine airfoil surface drag reduction rib size optimization design method comprises the steps that the initial dimensioned size of a drag reduction rib and a first frictional resistance coefficient of a wind turbine airfoil surface are determined; on the basis of the initial dimensional size and the first frictional resistance coefficient, the dimensionless size of the drag reduction rib located on the airfoil surface of the wind turbine is determined; determining a first drag reduction rate according to the dimensionless size; determining the total drag reduction rate of the drag reduction rib on the airfoil surface of the wind turbine by combining the first frictional resistance coefficient and the first drag reduction rate; determining a corresponding target dimensional size when the total drag reduction rate is maximum; the dimension of the drag reduction rib is associated with the drag reduction effect generated when the drag reduction rib is applied to the wind turbine airfoil, the size of the drag reduction rib is optimally designed, rapid design of the size of the rib is achieved, the drag reduction effect of the airfoil is optimal, and the rib can be conveniently applied to actual engineering.

Description

Technical Field

[0001] The present invention relates to the technical field of wind power generation, and in particular to a method for optimizing the design of drag reduction ribs on the surface of a wind turbine airfoil, a device for optimizing the design of drag reduction ribs on the surface of a wind turbine airfoil, an electronic device and a storage medium. Background Art

[0002] As one of the core components of wind turbines, improving the aerodynamic performance of wind turbine blades can greatly improve the overall operating efficiency of the unit, and blade resistance is one of the main factors restricting the improvement of blade aerodynamic performance. The blades are composed of a series of airfoils. Researchers have proposed using a shark skin-like rib structure to reduce the wind turbine airfoil resistance and achieve the purpose of blade drag reduction. Studies have shown that it has a good drag reduction effect. The drag reduction performance of the ribs is affected by their structural shape and size, and is directly related to the incoming flow state. Therefore, when using the rib structure to reduce the wind turbine airfoil drag, it is necessary to optimize the rib size under the condition of considering the influence of multiple factors in order to achieve the optimal drag reduction effect.

[0003] Since the operating Reynolds number of wind turbine airfoils is relatively high, generally in the millions, the ribs used on wind turbine airfoils are very small, usually in the micron range, which increases the difficulty of experimental testing. In numerical calculations, the mesh volume of airfoils with ribs increases dramatically due to their small size, which consumes a lot of computing resources and restricts the actual engineering application of the ribs. Summary of the invention

[0004] In view of the above problems, embodiments of the present invention are proposed to provide a method for optimizing the size of drag reduction ribs on the surface of a wind turbine airfoil, a device for optimizing the size of drag reduction ribs on the surface of a wind turbine airfoil, an electronic device and a storage medium that overcome the above problems or at least partially solve the above problems.

[0005] In order to solve the above problems, in a first aspect of the present invention, an embodiment of the present invention discloses a method for optimizing the size of drag reduction ribs on the surface of a wind turbine airfoil, comprising:

[0006] Determining the initial dimensioned size of the drag reducing rib and the first friction drag coefficient of the wind turbine airfoil surface;

[0007] Determining a dimensionless size of the drag reduction rib located on the surface of the wind turbine airfoil based on the initial dimensioned size and the first friction drag coefficient;

[0008] determining a first drag reduction ratio based on the dimensionless size;

[0009] Determine the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface in combination with the first friction drag coefficient and the first drag reduction rate;

[0010] Determine the target dimensioned size corresponding to the maximum total drag reduction rate.

[0011] Optionally, the step of determining the initial dimensioned size of the drag reducing rib and the first friction drag coefficient of the airfoil surface comprises:

[0012] Calculating the square root of the cross-sectional area of ​​the drag reducing rib as the initial dimensioned size;

[0013] Modeling the wind turbine airfoil surface to generate a wind turbine wing model;

[0014] A numerical simulation calculation is performed on the wind turbine wing model to determine the first friction resistance coefficient.

[0015] Optionally, the step of determining the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface based on the initial dimensioned size and the first friction drag coefficient comprises:

[0016] Based on the dimensionless dimension formula, the initial dimensioned dimension and the first friction resistance coefficient are calculated to determine the dimensionless dimension; wherein the dimensionless dimension formula is:

[0017]

[0018] lg is the initial dimensioned size, lg + is a dimensionless dimension, u τ is the friction velocity of the wind turbine airfoil surface, ν is the air kinematic viscosity, U ∞ is the incoming wind speed, C f,s is the first friction resistance coefficient.

[0019] Optionally, the step of determining a first drag reduction rate according to the dimensionless size comprises:

[0020] A preset rib drag reduction curve is looked up in a table based on the dimensionless size to determine that the drag reduction rate corresponding to the dimensionless size is the first drag reduction rate.

[0021] Optionally, the step of determining the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction rate comprises:

[0022] determining a second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate;

[0023] Integrating the second friction resistance coefficient to determine the viscous resistance;

[0024] The total drag reduction rate is determined based on the viscous drag.

[0025] Optionally, the step of determining the second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate includes:

[0026] Based on the friction resistance coefficient calculation formula, the second friction resistance coefficient is determined in combination with the first friction resistance coefficient and the first drag reduction rate, wherein the friction resistance coefficient calculation formula is:

[0027] C f,r =C f,s (1+DR),

[0028] C f,r The second friction coefficient, C f,s is the first friction resistance coefficient, and DR is the first drag reduction rate.

[0029] Optionally, the step of integrating and calculating the second friction resistance coefficient to determine the viscous resistance comprises:

[0030] The second friction resistance coefficient is integrated based on an integral formula to determine the viscous resistance, wherein the integral formula is:

[0031]

[0032] C D,r is the viscous resistance.

[0033] Optionally, the step of determining the total drag reduction rate according to the viscous resistance comprises:

[0034] The total drag reduction rate is determined based on a difference between the viscous drag and the viscous drag of the wind turbine airfoil.

[0035] Optionally, the step of determining the total drag reduction rate based on the difference between the viscous drag and the viscous drag of the wind turbine airfoil comprises:

[0036] Based on the viscous drag formula, the difference between the viscous drag and the viscous drag of the wind turbine airfoil and the percentage of the difference are calculated to determine the total drag reduction rate, wherein the viscous drag formula is:

[0037]

[0038] DR tota l is the total drag reduction rate, C D is the viscous drag of the wind turbine airfoil.

[0039] Optionally, the step of determining a target dimensioned size corresponding to the maximum total drag reduction rate comprises:

[0040] A preset target optimization algorithm is adopted, the maximum total drag reduction rate is taken as the optimization target of the preset target optimization algorithm, and the target dimension size is determined.

[0041] In a second aspect of the present invention, an embodiment of the present invention further discloses a device for optimizing the size of drag reduction ribs on an airfoil surface of a wind turbine, comprising:

[0042] A first determination module, used to determine the initial dimensioned size of the drag reduction rib and the first friction drag coefficient of the wind turbine airfoil surface;

[0043] A first calculation module, configured to determine a dimensionless size of the drag reduction rib located on the surface of the wind turbine airfoil based on the initial dimensioned size and the first friction drag coefficient;

[0044] A second determination module, configured to determine a first drag reduction rate according to the dimensionless size;

[0045] A second calculation module, configured to determine a total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction rate;

[0046] The third calculation module is used to determine the target dimensioned size corresponding to the maximum total drag reduction rate.

[0047] Optionally, the first determining module includes:

[0048] A first calculation submodule, used for calculating the square root of the cross-sectional area of ​​the drag reduction rib as the initial dimensioned size;

[0049] A modeling submodule, used for modeling the wind turbine airfoil surface to generate a wind turbine wing model;

[0050] The second calculation submodule is used to perform numerical simulation calculation on the wind turbine wing model to determine the first friction resistance coefficient.

[0051] Optionally, the second calculation module includes:

[0052] a fourth calculation submodule, configured to determine a second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate;

[0053] a fifth calculation submodule, configured to perform integral calculation on the second friction resistance coefficient to determine viscous resistance;

[0054] A sixth calculation submodule is used to determine the total drag reduction rate according to the viscous resistance.

[0055] In the third aspect of the present invention, an embodiment of the present invention further discloses an electronic device, comprising a processor, a memory, and a computer program stored in the memory and capable of running on the processor, wherein the computer program, when executed by the processor, implements the steps of the method for optimizing the size of the drag reducing ribs on the surface of a wind turbine airfoil as described above.

[0056] In the fourth aspect of the present invention, an embodiment of the present invention further discloses a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the steps of the method for optimizing the size of the drag reducing ribs on the surface of a wind turbine airfoil are implemented as described above.

[0057] The embodiments of the present invention include the following advantages:

[0058] The embodiment of the present invention determines the initial dimensional size of the drag reduction rib and the first friction resistance coefficient of the wind turbine airfoil surface; determines the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface based on the initial dimensional size and the first friction resistance coefficient; determines the first drag reduction rate according to the dimensionless size; determines the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface in combination with the first friction resistance coefficient and the first drag reduction rate; determines the target dimensional size corresponding to the maximum total drag reduction rate; and optimizes the size of the drag reduction rib by associating the dimensional size of the drag reduction rib with the drag reduction effect generated when it is applied to the wind turbine airfoil, thereby achieving rapid design of the rib size, optimizing the airfoil drag reduction effect, and facilitating the application of the rib in actual engineering. BRIEF DESCRIPTION OF THE DRAWINGS

[0059] Figure 1 It is a flow chart of steps of an embodiment of a method for optimizing the size of drag-reducing ribs on the surface of a wind turbine airfoil of the present invention;

[0060] Figure 2 It is a flowchart of steps of an example of a method for optimizing the size of drag-reducing ribs on the surface of a wind turbine airfoil of the present invention;

[0061] Figure 3 It is a distribution diagram of friction drag coefficient on the airfoil surface of an example of a method for optimizing the size of drag reduction ribs on the airfoil surface of a wind turbine according to the present invention;

[0062] Figure 4 It is a schematic cross-sectional view of a rib according to an example of a method for optimizing the size of a drag-reducing rib on a wind turbine airfoil surface of the present invention;

[0063] Figure 5 It is a schematic diagram of a rib drag reduction curve of an example of a method for optimizing the size of a drag reduction rib on a wind turbine airfoil surface of the present invention;

[0064] Figure 6It is a structural block diagram of an embodiment of a device for optimizing the size of drag-reducing ribs on the surface of a wind turbine airfoil of the present invention;

[0065] Figure 7 is a structural block diagram of an electronic device provided by an embodiment of the present invention;

[0066] Figure 8 It is a structural block diagram of a storage medium provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0067] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0068] Reference Figure 1 , shows a flow chart of steps of an embodiment of a method for optimizing the size of drag reduction ribs on the surface of a wind turbine airfoil of the present invention, wherein the method for optimizing the size of drag reduction ribs on the surface of a wind turbine airfoil may specifically include the following steps:

[0069] Step 101, determining the initial dimensioned size of the drag reduction rib and the first friction drag coefficient of the wind turbine airfoil surface;

[0070] In the embodiment of the present invention, the original dimensioned size of the drag reduction rib, such as the height, width, cross-sectional area, etc. of the drag reduction rib, can be determined as the initial dimensioned size. At the same time, the frictional resistance coefficient of the wind turbine airfoil surface can be determined when the wind turbine airfoil surface is smooth, that is, the first frictional resistance coefficient.

[0071] Step 102, determining a dimensionless size of the drag reduction rib located on the surface of the wind turbine airfoil based on the initial dimensioned size and the first friction drag coefficient;

[0072] Since the initial dimensional size has limitations such as units, it cannot be applied to various scenarios. In order to accurately calculate, the initial dimensional size and the first friction resistance coefficient can be converted into the dimensionless size of the drag reduction rib located on the surface of the wind turbine airfoil, so as to use dimensionless parameters to describe the characteristics of the deceleration rib.

[0073] Step 103, determining a first drag reduction rate according to the dimensionless size;

[0074] Based on the characteristics of the deceleration rib corresponding to the dimensionless size, the drag reduction rate of the deceleration rib at various positions on the wind turbine airfoil surface, ie, the first drag reduction rate, is determined.

[0075] Step 104, determining a total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction rate;

[0076] The first friction drag coefficient and the first drag reduction rate are calculated together to determine the overall total drag reduction rate of the drag reduction ribs located on the wind turbine airfoil surface.

[0077] Step 105, determining the target dimensioned size corresponding to the maximum total drag reduction rate.

[0078] Based on the optimal range of rib size, the rib size on the wind turbine airfoil surface is optimized and designed, and the target dimensional size corresponding to the maximum total drag reduction rate is determined, so as to obtain the size of the deceleration rib with the best drag reduction effect.

[0079] In an optional embodiment of the present invention, the step of determining the initial dimensioned size of the drag reducing rib and the first friction drag coefficient of the airfoil surface comprises:

[0080] Sub-step S1011, calculating the square root of the cross-sectional area of ​​the drag reduction rib as the initial dimensioned size;

[0081] In practical applications, the square root of the cross-sectional area of ​​the drag reduction rib can be determined as the initial dimensional size.

[0082] Sub-step S1012, modeling the wind turbine airfoil surface to generate a wind turbine wing model;

[0083] The wind turbine airfoil surface can be modeled, and the wind turbine wing model can be generated through processes such as geometric shape modeling, unit mesh division, laminate material laying, finite element calculation and post-processing analysis.

[0084] Sub-step S1013, performing numerical simulation calculation on the wind turbine wing model to determine the first friction resistance coefficient.

[0085] Then, numerical simulation calculation is performed on the wind turbine wing model to determine the friction drag coefficient of the smooth wind turbine airfoil surface, that is, the first friction drag coefficient.

[0086] In addition, the first friction resistance coefficient may also be calculated according to an empirical formula or an experimental measurement method, which is not limited in the embodiment of the present invention.

[0087] In an optional embodiment of the present invention, the step of determining the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface based on the initial dimensioned size and the first friction drag coefficient comprises: calculating the initial dimensioned size and the first friction drag coefficient based on a dimensionless size formula to determine the dimensionless size; wherein the dimensionless size formula is:

[0088]

[0089] lg is the initial dimensioned size, lg +is a dimensionless dimension, u τ is the friction velocity of the wind turbine airfoil surface, ν is the air kinematic viscosity, U ∞ is the incoming wind speed, C f,s is the first friction resistance coefficient.

[0090] For the calculation of dimensionless dimensions, the initial dimensioned dimensions can be substituted into the dimensionless dimension formula to calculate the actual dimensionless dimension lg + .

[0091] In an optional embodiment of the present invention, the step of determining the first drag reduction rate based on the dimensionless size includes: looking up a preset rib drag reduction curve based on the dimensionless size, and determining that the drag reduction rate corresponding to the dimensionless size is the first drag reduction rate.

[0092] The first drag reduction rate of the speed reduction rib depends on the dimensionless size of the rib, lg + According to the dimensionless size of the deceleration rib at various positions on the wind turbine airfoil surface, the drag reduction curve (DR-lg + ) obtains the drag reduction rate corresponding to the dimensionless size as the first drag reduction rate. The preset rib drag reduction curve is determined by a preliminary experiment, and the process of determining the preset rib drag reduction curve can be determined according to actual conditions. The embodiment of the present invention does not limit the process of determining the preset rib drag reduction curve.

[0093] In an optional embodiment of the present invention, the step of determining the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction rate comprises:

[0094] Sub-step S1041, determining a second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate;

[0095] Firstly, the frictional resistance coefficient of the retrieval rib on the wind turbine airfoil surface, ie, the second frictional resistance coefficient, can be determined based on the first frictional resistance coefficient and the first drag reduction rate.

[0096] Specifically, the step of determining the second friction resistance coefficient based on the first friction resistance coefficient and the first drag reduction rate includes: determining the second friction resistance coefficient based on a friction resistance coefficient calculation formula in combination with the first friction resistance coefficient and the first drag reduction rate, wherein the friction resistance coefficient calculation formula is:

[0097] C f,r =C f,s (1+DR), (2)

[0098] C f,r The second friction coefficient, C f,sis the first friction resistance coefficient, and DR is the first drag reduction rate.

[0099] The first friction resistance coefficient and the first drag reduction rate can be substituted into the friction resistance coefficient calculation formula, and the calculated result is the second friction resistance coefficient.

[0100] Sub-step S1042, performing integral calculation on the second friction resistance coefficient to determine viscous resistance;

[0101] The second friction drag coefficient is integrated and calculated to obtain the viscous drag of the wind turbine airfoil with ribs.

[0102] Specifically, the step of integrating the second friction resistance coefficient to determine the viscous resistance includes: integrating the second friction resistance coefficient based on an integral formula to determine the viscous resistance, wherein the integral formula is:

[0103]

[0104] C D,r is the viscous resistance.

[0105] Substitute the second friction resistance coefficient into the integral formula, then solve the integral formula, and the result obtained is the viscous resistance.

[0106] Sub-step S1043, determining the total drag reduction rate according to the viscous resistance.

[0107] The viscous drag is calculated to determine the overall drag reduction rate of the speed reduction ribs on the wind turbine airfoil, that is, the total drag reduction rate.

[0108] Specifically, the step of determining the total drag reduction rate according to the viscous drag includes: determining the total drag reduction rate based on the difference between the viscous drag and the viscous drag of the wind turbine airfoil.

[0109] For the calculation of the total drag reduction rate, the viscous drag and the viscous drag corresponding to the wind turbine airfoil in a smooth state can be determined, the difference between the two can be calculated, and the total drag reduction rate can be determined based on the size of the difference.

[0110] Further, the step of determining the total drag reduction rate based on the difference between the viscous drag and the viscous drag of the wind turbine airfoil comprises: calculating the difference between the viscous drag and the viscous drag of the wind turbine airfoil and the percentage of the difference based on a viscous drag formula to determine the total drag reduction rate, wherein the viscous drag formula is:

[0111]

[0112] DR total is the total drag reduction rate, C Dis the viscous drag of the wind turbine airfoil.

[0113] In practical applications, the viscous drag can be substituted into the viscous drag formula for calculation, and the result of the calculation is the total drag reduction rate.

[0114] In an optional embodiment of the present invention, the step of determining the target dimensional size corresponding to the maximum total drag reduction rate includes: using a preset target optimization algorithm, taking the maximum total drag reduction rate as the optimization target of the preset target optimization algorithm, and determining the target dimensional size.

[0115] The target optimization algorithm can be preset, the optimization range of the rib size can be set, the rib size of the wind turbine airfoil surface can be optimized, and the rib size with the best drag reduction effect can be obtained. The input parameters of the preset target optimization algorithm are the first drag reduction rate, the incoming wind speed, the preset rib drag reduction curve, the independent variable is the dimensioned size, and the optimization target is the overall drag reduction rate DR of the airfoil total The maximum is solved based on the preset target optimization algorithm, and the corresponding dimensional size is the target dimensional size.

[0116] The embodiment of the present invention determines the initial dimensional size of the drag reduction rib and the first friction resistance coefficient of the wind turbine airfoil surface; determines the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface based on the initial dimensional size and the first friction resistance coefficient; determines the first drag reduction rate according to the dimensionless size; determines the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface in combination with the first friction resistance coefficient and the first drag reduction rate; determines the target dimensional size corresponding to the maximum total drag reduction rate; and optimizes the size of the drag reduction rib by associating the dimensional size of the drag reduction rib with the drag reduction effect generated when it is applied to the wind turbine airfoil, thereby achieving rapid design of the rib size, optimizing the airfoil drag reduction effect, and facilitating the application of the rib in actual engineering.

[0117] In order to make the implementation process of the embodiments of the present invention clear to those skilled in the art, refer to Figure 2 Let me explain this with an example:

[0118] The wind turbine airfoil is DU-96-W-18, and the airfoil is modeled. The friction coefficient C of the airfoil surface is obtained by numerical simulation. f,s The numerical calculation model adopts the Transition-SST model, the airfoil chord length is c = 1m, the airfoil angle of attack is α = 0°, and the incoming wind speed is U ∞ The airfoil surface friction coefficient is calculated as follows: Figure 3 shown.

[0119] Select a rib structure with a trapezoidal cross section, as shown in the attached Figure 4 As shown, an initial rib size lg is set, and the dimensionless size lg of the rib structure at each position on the airfoil surface is calculated using formula (1): + Then, using formula (4), we can get the overall drag reduction effect DR of the rib structure on the wind turbine airfoil: total .

[0120] Adaptive simulated annealing algorithm is used as the optimization design method, and the rib optimization range is set to lg = 10 ~ 200μm. The input parameters of the optimization design method are the friction coefficient C of the wind turbine airfoil surface f,s , Inflow wind speed U ∞ , rib drag reduction curve DR-lg + See attached Figure 5 , the independent variable is the dimensioned size lg, and the optimization target is the overall drag reduction rate DR of the airfoil total maximum.

[0121] Two optimization cases are calculated: the case of full rib coverage and the case of partial rib coverage of 0-0.5c airfoil. The specific settings are shown in Table 1:

[0122] Case Rib optimization range lg / μm Coverage area Case 1 10～200 Full coverage Case 2 10～200 0~0.5c Partial coverage

[0123] Table 1

[0124] Case Rib size lg / μm <![CDATA[DR total ]]> Case 1 78 -5.01％ Case 2 74 -3.69％

[0125] Table 2

[0126] It can be seen that the optimal rib size lg = 78μm for Case 1, and the total drag reduction rate is 5.01%. The optimal rib size lg = 74μm for Case 2, and the total drag reduction rate is 3.69%. Both have achieved significant drag reduction effects, proving that the patented method can achieve the optimal design of the drag reduction rib size on the wind turbine airfoil surface.

[0127] It should be noted that, for the sake of simplicity, the method embodiments are described as a series of action combinations, but those skilled in the art should be aware that the embodiments of the present invention are not limited by the order of the actions described, because according to the embodiments of the present invention, certain steps can be performed in other orders or simultaneously. Secondly, those skilled in the art should also be aware that the embodiments described in the specification are all preferred embodiments, and the actions involved are not necessarily required by the embodiments of the present invention.

[0128] Reference Figure 6 , shows a structural block diagram of an embodiment of a wind turbine airfoil surface drag reduction rib size optimization design device of the present invention, the wind turbine airfoil surface drag reduction rib size optimization design device may specifically include the following modules:

[0129] A first determination module 601 is used to determine the initial dimensioned size of the drag reduction rib and the first friction drag coefficient of the wind turbine airfoil surface;

[0130] A first calculation module 602 is used to determine the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface based on the initial dimensioned size and the first friction drag coefficient;

[0131] A second determination module 603, configured to determine a first drag reduction rate according to the dimensionless size;

[0132] A second calculation module 604 is used to determine a total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction rate;

[0133] The third calculation module 605 is used to determine the target dimensioned size corresponding to the maximum total drag reduction rate.

[0134] In an optional embodiment of the present invention, the first determining module 601 includes:

[0135] A first calculation submodule, used for calculating the square root of the cross-sectional area of ​​the drag reduction rib as the initial dimensioned size;

[0136] A modeling submodule, used for modeling the wind turbine airfoil surface to generate a wind turbine wing model;

[0137] The second calculation submodule is used to perform numerical simulation calculation on the wind turbine wing model to determine the first friction resistance coefficient.

[0138] In an optional embodiment of the present invention, the first calculation module 602 includes:

[0139] The third calculation submodule is used to calculate the initial dimensioned size and the first friction resistance coefficient based on a dimensionless size formula to determine the dimensionless size; wherein the dimensionless size formula is:

[0140]

[0141] lg is the initial dimensioned size, lg + is a dimensionless dimension, u τ is the friction velocity of the wind turbine airfoil surface, ν is the air kinematic viscosity, U ∞ is the incoming wind speed, C f,s is the first friction resistance coefficient.

[0142] In an optional embodiment of the present invention, the second determining module 603 includes:

[0143] A table lookup submodule is used to look up a preset rib drag reduction curve based on the dimensionless size, and determine that the drag reduction rate corresponding to the dimensionless size is the first drag reduction rate

[0144] In an optional embodiment of the present invention, the second calculation module 604 includes:

[0145] a fourth calculation submodule, configured to determine a second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate;

[0146] a fifth calculation submodule, configured to perform integral calculation on the second friction resistance coefficient to determine viscous resistance;

[0147] A sixth calculation submodule is used to determine the total drag reduction rate according to the viscous resistance.

[0148] In an optional embodiment of the present invention, the fourth calculation submodule includes:

[0149] A friction calculation unit is used to determine a second friction resistance coefficient based on a friction resistance coefficient calculation formula in combination with the first friction resistance coefficient and the first drag reduction rate, wherein the friction resistance coefficient calculation formula is:

[0150] C f,r =C f,s (1+DR),

[0151] C f,r The second friction coefficient, C f,s is the first friction resistance coefficient, and DR is the first drag reduction rate.

[0152] In an optional embodiment of the present invention, the fifth calculation submodule includes:

[0153] An integration unit is used to integrate the second friction resistance coefficient based on an integral formula to determine the viscous resistance, wherein the integral formula is:

[0154]

[0155] C D,r is the viscous resistance.

[0156] In an optional embodiment of the present invention, the sixth calculation submodule includes:

[0157] The total drag reduction rate calculation unit is used to determine the total drag reduction rate based on the difference between the viscous drag and the viscous drag of the wind turbine airfoil.

[0158] In an optional embodiment of the present invention, the total drag reduction rate calculation unit includes:

[0159] The total drag reduction rate calculation subunit is used to calculate the difference between the viscous drag and the viscous drag of the wind turbine airfoil, and the percentage of the difference, based on the viscous drag formula, to determine the total drag reduction rate, wherein the viscous drag formula is:

[0160]

[0161] DR total is the total drag reduction rate, C D is the viscous drag of the wind turbine airfoil.

[0162] In an optional embodiment of the present invention, the third calculation module 605 includes:

[0163] The target optimization submodule is used to adopt a preset target optimization algorithm, take the maximum total drag reduction rate as the optimization target of the preset target optimization algorithm, and determine the dimensioned size of the target.

[0164] The embodiment of the present invention determines the initial dimensional size of the drag reduction rib and the first friction resistance coefficient of the wind turbine airfoil surface; determines the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface based on the initial dimensional size and the first friction resistance coefficient; determines the first drag reduction rate according to the dimensionless size; determines the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface in combination with the first friction resistance coefficient and the first drag reduction rate; determines the target dimensional size corresponding to the maximum total drag reduction rate; and optimizes the size of the drag reduction rib by associating the dimensional size of the drag reduction rib with the drag reduction effect generated when it is applied to the wind turbine airfoil, thereby achieving rapid design of the rib size, optimizing the airfoil drag reduction effect, and facilitating the application of the rib in actual engineering.

[0165] As for the device embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0166] Reference Figure 7 , an embodiment of the present invention further provides an electronic device, including:

[0167] A processor 701 and a storage medium 702, wherein the storage medium 702 stores a computer program executable by the processor 701, and when the electronic device is running, the processor 701 executes the computer program to perform the method for optimizing the size of the drag reduction ribs on the surface of a wind turbine airfoil as described in any one of the embodiments of the present invention.

[0168] The wind turbine airfoil surface drag reduction rib size optimization design method comprises:

[0169] Determining the initial dimensioned size of the drag reducing rib and the first friction drag coefficient of the wind turbine airfoil surface;

[0170] Determine the dimensionless size of the drag reduction rib located on the surface of the wind turbine airfoil based on the initial dimensioned size and the first friction drag coefficient;

[0171] determining a first drag reduction ratio based on the dimensionless size;

[0172] Determine the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface in combination with the first friction drag coefficient and the first drag reduction rate;

[0173] Determine the target dimensioned size corresponding to the maximum total drag reduction rate.

[0174] Optionally, the step of determining the initial dimensioned size of the drag reducing rib and the first friction drag coefficient of the airfoil surface comprises:

[0175] Calculating the square root of the cross-sectional area of ​​the drag reducing rib as the initial dimensioned size;

[0176] Modeling the wind turbine airfoil surface to generate a wind turbine wing model;

[0177] A numerical simulation calculation is performed on the wind turbine wing model to determine the first friction resistance coefficient.

[0178] Optionally, the step of determining the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface based on the initial dimensioned size and the first friction drag coefficient comprises:

[0179] Based on the dimensionless dimension formula, the initial dimensioned dimension and the first friction resistance coefficient are calculated to determine the dimensionless dimension; wherein the dimensionless dimension formula is:

[0180]

[0181] lg is the initial dimensioned size, lg + is a dimensionless dimension, u τ is the friction velocity of the wind turbine airfoil surface, ν is the air kinematic viscosity, U ∞ is the incoming wind speed, C f,s is the first friction coefficient.

[0182] Optionally, the step of determining a first drag reduction rate according to the dimensionless size comprises:

[0183] A preset rib drag reduction curve is looked up in a table based on the dimensionless size to determine that the drag reduction rate corresponding to the dimensionless size is the first drag reduction rate.

[0184] Optionally, the step of determining the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction rate comprises:

[0185] determining a second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate;

[0186] Integrating the second friction resistance coefficient to determine the viscous resistance;

[0187] The total drag reduction rate is determined based on the viscous drag.

[0188] Optionally, the step of determining the second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate includes:

[0189] Based on the friction resistance coefficient calculation formula, the second friction resistance coefficient is determined in combination with the first friction resistance coefficient and the first drag reduction rate, wherein the friction resistance coefficient calculation formula is:

[0190] C f,r =C f,s (1+DR),

[0191] C f,r The second friction coefficient, C f,s is the first friction resistance coefficient, and DR is the first drag reduction rate.

[0192] Optionally, the step of integrating and calculating the second friction resistance coefficient to determine the viscous resistance comprises:

[0193] The second friction resistance coefficient is integrated based on an integral formula to determine the viscous resistance, wherein the integral formula is:

[0194]

[0195] C D,r is the viscous resistance.

[0196] Optionally, the step of determining the total drag reduction rate according to the viscous resistance comprises:

[0197] The total drag reduction rate is determined based on a difference between the viscous drag and the viscous drag of the wind turbine airfoil.

[0198] Optionally, the step of determining the total drag reduction rate based on the difference between the viscous drag and the viscous drag of the wind turbine airfoil comprises:

[0199] Based on the viscous drag formula, the difference between the viscous drag and the viscous drag of the wind turbine airfoil and the percentage of the difference are calculated to determine the total drag reduction rate, wherein the viscous drag formula is:

[0200]

[0201] DR total is the total drag reduction rate, C D is the viscous drag of the wind turbine airfoil.

[0202] Optionally, the step of determining a target dimensioned size corresponding to the maximum total drag reduction rate comprises:

[0203] A preset target optimization algorithm is adopted, the maximum total drag reduction rate is taken as the optimization target of the preset target optimization algorithm, and the target dimension size is determined.

[0204] The embodiment of the present invention determines the initial dimensional size of the drag reduction rib and the first friction resistance coefficient of the wind turbine airfoil surface; determines the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface based on the initial dimensional size and the first friction resistance coefficient; determines the first drag reduction rate according to the dimensionless size; determines the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface in combination with the first friction resistance coefficient and the first drag reduction rate; determines the target dimensional size corresponding to the maximum total drag reduction rate; and optimizes the size of the drag reduction rib by associating the dimensional size of the drag reduction rib with the drag reduction effect generated when it is applied to the wind turbine airfoil, thereby achieving rapid design of the rib size, optimizing the airfoil drag reduction effect, and facilitating the application of the rib in actual engineering.

[0205] The memory may include a random access memory (RAM) or a non-volatile memory, such as at least one disk memory. Optionally, the memory may also be at least one storage device located away from the processor.

[0206] The above-mentioned processor can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it can also be a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field programmable gate array (FPGA) or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.

[0207] Reference Figure 8The embodiment of the present invention further provides a computer-readable storage medium 801, on which a computer program is stored. When the computer program is executed by a processor, the method for optimizing the size of the drag reduction ribs on the surface of a wind turbine airfoil as described in any one of the embodiments of the present invention is executed.

[0208] The wind turbine airfoil surface drag reduction rib size optimization design method comprises:

[0209] Determining the initial dimensioned size of the drag reducing rib and the first friction drag coefficient of the wind turbine airfoil surface;

[0210] Determine the dimensionless size of the drag reduction rib located on the surface of the wind turbine airfoil based on the initial dimensioned size and the first friction drag coefficient;

[0211] determining a first drag reduction ratio based on the dimensionless size;

[0212] Determine the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction rate;

[0213] Determine the target dimensioned size corresponding to the maximum total drag reduction rate.

[0214] Optionally, the step of determining the initial dimensioned size of the drag reducing rib and the first friction drag coefficient of the airfoil surface comprises:

[0215] Calculating the square root of the cross-sectional area of ​​the drag reducing rib as the initial dimensioned size;

[0216] Modeling the wind turbine airfoil surface to generate a wind turbine wing model;

[0217] A numerical simulation calculation is performed on the wind turbine wing model to determine the first friction resistance coefficient.

[0218] Optionally, the step of determining the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface based on the initial dimensioned size and the first friction drag coefficient comprises:

[0219] Based on the dimensionless dimension formula, the initial dimensioned dimension and the first friction resistance coefficient are calculated to determine the dimensionless dimension; wherein the dimensionless dimension formula is:

[0220]

[0221] lg is the initial dimensioned size, lg + is a dimensionless dimension, u τ is the friction velocity of the wind turbine airfoil surface, ν is the air kinematic viscosity, U ∞ is the incoming wind speed, C f,sis the first friction resistance coefficient.

[0222] Optionally, the step of determining a first drag reduction rate according to the dimensionless size comprises:

[0223] A preset rib drag reduction curve is looked up in a table based on the dimensionless size to determine that the drag reduction rate corresponding to the dimensionless size is the first drag reduction rate.

[0224] Optionally, the step of determining the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction rate comprises:

[0225] determining a second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate;

[0226] Integrating the second friction resistance coefficient to determine the viscous resistance;

[0227] The total drag reduction rate is determined based on the viscous drag.

[0228] Optionally, the step of determining the second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate includes:

[0229] Based on the friction resistance coefficient calculation formula, the second friction resistance coefficient is determined in combination with the first friction resistance coefficient and the first drag reduction rate, wherein the friction resistance coefficient calculation formula is:

[0230] C f,r =C f,s (1+DR),

[0231] C f,r The second friction coefficient, C f,s is the first friction resistance coefficient, and DR is the first drag reduction rate.

[0232] Optionally, the step of integrating and calculating the second friction resistance coefficient to determine the viscous resistance comprises:

[0233] The second friction resistance coefficient is integrated based on an integral formula to determine the viscous resistance, wherein the integral formula is:

[0234]

[0235] C D,r is the viscous resistance.

[0236] Optionally, the step of determining the total drag reduction rate according to the viscous resistance comprises:

[0237] The total drag reduction rate is determined based on a difference between the viscous drag and the viscous drag of the wind turbine airfoil.

[0238] Optionally, the step of determining the total drag reduction rate based on the difference between the viscous drag and the viscous drag of the wind turbine airfoil comprises:

[0239] Based on the viscous drag formula, the difference between the viscous drag and the viscous drag of the wind turbine airfoil and the percentage of the difference are calculated to determine the total drag reduction rate, wherein the viscous drag formula is:

[0240]

[0241] DR total is the total drag reduction rate, C D is the viscous drag of the wind turbine airfoil.

[0242] Optionally, the step of determining a target dimensioned size corresponding to the maximum total drag reduction rate comprises:

[0243] A preset target optimization algorithm is adopted, the maximum total drag reduction rate is taken as the optimization target of the preset target optimization algorithm, and the target dimension size is determined.

[0244] The embodiment of the present invention determines the initial dimensional size of the drag reduction rib and the first friction resistance coefficient of the wind turbine airfoil surface; determines the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface based on the initial dimensional size and the first friction resistance coefficient; determines the first drag reduction rate according to the dimensionless size; determines the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface in combination with the first friction resistance coefficient and the first drag reduction rate; determines the target dimensional size corresponding to the maximum total drag reduction rate; and optimizes the size of the drag reduction rib by associating the dimensional size of the drag reduction rib with the drag reduction effect generated when it is applied to the wind turbine airfoil, thereby achieving rapid design of the rib size, optimizing the airfoil drag reduction effect, and facilitating the application of the rib in actual engineering.

[0245] The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.

[0246] It will be appreciated by those skilled in the art that the embodiments of the present invention may be provided as methods, devices, or computer program products. Therefore, the embodiments of the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the embodiments of the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.

[0247] The embodiments of the present invention are described with reference to the flowcharts and / or block diagrams of the methods, terminal devices (systems), and computer program products according to the embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of the processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing terminal device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing terminal device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0248] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing terminal device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a manufactured product including an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0249] These computer program instructions can also be loaded onto a computer or other programmable data processing terminal device so that a series of operating steps are executed on the computer or other programmable terminal device to produce a computer-implemented process, thereby providing instructions for implementing the process in the computer or other programmable terminal device. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0250] Although the preferred embodiments of the present invention have been described, those skilled in the art may make additional changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the embodiments of the present invention.

[0251] Finally, it should be noted that, in this article, relational terms such as first and second, etc. 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 terminal 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 terminal device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the existence of other identical elements in the process, method, article or terminal device including the elements.

[0252] The above describes in detail a method for optimizing the size of drag reduction ribs on the surface of a wind turbine airfoil, a device for optimizing the size of drag reduction ribs on the surface of a wind turbine airfoil, an electronic device and a storage medium provided by the present invention. Specific examples are used herein to illustrate the principles and implementation methods of the present invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on the present invention.

Claims

1. A method for optimizing the size of drag reduction ribs on a wind turbine airfoil surface, characterized in that: include: Determining the initial dimensioned size of the drag reducing rib and the first friction drag coefficient of the wind turbine airfoil surface; Determining a dimensionless size of the drag reduction rib located on the surface of the wind turbine airfoil based on the initial dimensioned size and the first friction drag coefficient; determining a first drag reduction ratio based on the dimensionless size; Determine the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface in combination with the first friction drag coefficient and the first drag reduction rate; Determine the target dimensioned size corresponding to the maximum total drag reduction rate.

2. The method according to claim 1, characterized in that The step of determining the initial dimensioned size of the drag reducing rib and the first friction drag coefficient of the airfoil surface comprises: Calculating the square root of the cross-sectional area of ​​the drag reducing rib as the initial dimensioned size; Modeling the wind turbine airfoil surface to generate a wind turbine wing model; A numerical simulation calculation is performed on the wind turbine wing model to determine the first friction resistance coefficient.

3. The method according to claim 2, characterized in that The step of determining the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface based on the initial dimensioned size and the first friction drag coefficient comprises: Based on the dimensionless dimension formula, the initial dimensioned dimension and the first friction resistance coefficient are calculated to determine the dimensionless dimension; wherein the dimensionless dimension formula is: lg is the initial dimensioned size, lg + is a dimensionless dimension, u τ is the friction velocity of the wind turbine airfoil surface, ν is the air kinematic viscosity, U ∞ is the incoming wind speed, C f,s is the first friction resistance coefficient.

4. The method according to any one of claims 1 to 3, characterized in that: The step of determining a first drag reduction rate according to the dimensionless size comprises: A preset rib drag reduction curve is looked up in a table based on the dimensionless size to determine that the drag reduction rate corresponding to the dimensionless size is the first drag reduction rate.

5. The method according to any one of claims 1 to 3, characterized in that: The step of determining the total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction rate comprises: determining a second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate; Integrating the second friction resistance coefficient to determine the viscous resistance; The total drag reduction rate is determined based on the viscous drag.

6. The method according to claim 5, characterized in that The step of determining the second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate comprises: Based on the friction resistance coefficient calculation formula, the second friction resistance coefficient is determined in combination with the first friction resistance coefficient and the first drag reduction rate, wherein the friction resistance coefficient calculation formula is: C f,r =C f,s (1+DR), C f,r The second friction coefficient, C f,s is the first friction resistance coefficient, and DR is the first drag reduction rate.

7. The method according to claim 6, characterized in that The step of integrating and calculating the second friction resistance coefficient to determine the viscous resistance comprises: The second friction resistance coefficient is integrated based on an integral formula to determine the viscous resistance, wherein the integral formula is: C D,r is the viscous resistance.

8. The method according to claim 7, characterized in that The step of determining the total drag reduction rate according to the viscous resistance comprises: The total drag reduction rate is determined based on a difference between the viscous drag and the viscous drag of the wind turbine airfoil.

9. The method according to claim 7, characterized in that: The step of determining the total drag reduction rate based on the difference between the viscous drag and the viscous drag of the wind turbine airfoil comprises: Based on the viscous drag formula, the difference between the viscous drag and the viscous drag of the wind turbine airfoil and the percentage of the difference are calculated to determine the total drag reduction rate, wherein the viscous drag formula is: DR tota l is the total drag reduction rate, C D is the viscous drag of the wind turbine airfoil.

10. The method according to claim 1, characterized in that The step of determining the target dimensioned size corresponding to the maximum total drag reduction rate comprises: A preset target optimization algorithm is adopted, the maximum total drag reduction rate is taken as the optimization target of the preset target optimization algorithm, and the target dimension size is determined.

11. A wind turbine airfoil surface drag reduction rib size optimization design device, characterized in that: include: A first determination module, used to determine the initial dimensioned size of the drag reduction rib and the first friction drag coefficient of the wind turbine airfoil surface; A first calculation module, configured to determine a dimensionless size of the drag reduction rib located on the surface of the wind turbine airfoil based on the initial dimensioned size and the first friction drag coefficient; A second determination module, configured to determine a first drag reduction rate according to the dimensionless size; A second calculation module, configured to determine a total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction rate; The third calculation module is used to determine the target dimensioned size corresponding to the maximum total drag reduction rate.

12. The device according to claim 11, characterized in that The first determining module comprises: A first calculation submodule, used for calculating the square root of the cross-sectional area of ​​the drag reduction rib as the initial dimensioned size; A modeling submodule, used for modeling the wind turbine airfoil surface to generate a wind turbine wing model; The second calculation submodule is used to perform numerical simulation calculation on the wind turbine wing model to determine the first friction resistance coefficient.

13. The device according to claim 11, characterized in that The second calculation module includes: a fourth calculation submodule, configured to determine a second friction resistance coefficient according to the first friction resistance coefficient and the first drag reduction rate; a fifth calculation submodule, configured to perform integral calculation on the second friction resistance coefficient to determine viscous resistance; A sixth calculation submodule is used to determine the total drag reduction rate according to the viscous resistance.

14. An electronic device, characterized in that: The invention comprises a processor, a memory and a computer program stored in the memory and capable of being run on the processor, wherein when the computer program is executed by the processor, the steps of the method for optimizing the size of the drag reducing ribs on the surface of a wind turbine airfoil are implemented as claimed in any one of claims 1 to 10.

15. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method for optimizing the size of the drag reducing ribs on the surface of a wind turbine airfoil according to any one of claims 1 to 10 are implemented.

Images