A method and device for optimizing the size of a wind turbine airfoil surface drag reduction rib

Abstract

The embodiment of the application provides a kind of wind turbine airfoil surface drag reduction rib size optimization design method and device, comprising: determining the initial dimensionless size of drag reduction rib and the first friction drag coefficient of wind turbine airfoil surface;Based on the initial dimensionless size and the first friction drag coefficient, determine the dimensionless size of the drag reduction rib located on the wind turbine airfoil surface;According to the dimensionless size, determine the first drag reduction rate;The total drag reduction rate of the drag reduction rib located on the wind turbine airfoil surface is determined by combining the first friction drag coefficient and the first drag reduction rate;Determine the target dimensionless size corresponding to the maximum total drag reduction rate;By associating the dimensionless size of drag reduction rib with the drag reduction effect generated when it is applied to wind turbine airfoil, the size of drag reduction rib is optimized and designed, the rib size is quickly designed, the airfoil drag reduction effect is optimal, and the rib is applied in actual engineering.

Description

Technical Field

[0001] This invention relates to the field of wind power generation technology, and in particular to a method for optimizing the design of drag-reducing rib dimensions on the surface of a wind turbine airfoil, a device for optimizing the design of drag-reducing rib dimensions on the surface of a wind turbine airfoil, an electronic device, and a storage medium. Background Technology

[0002] As one of the core components of wind turbine generators, improving the aerodynamic performance of wind turbine blades can greatly enhance the overall operating efficiency of the unit. Blade drag is one of the main factors restricting the improvement of blade aerodynamic performance. Blades consist of a series of airfoils. Researchers have proposed using a sharkskin-inspired rib structure to reduce airfoil drag in wind turbines, thereby achieving blade drag reduction. Studies have shown that this method 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 conditions. Therefore, when using rib structures to reduce airfoil drag in wind turbines, it is necessary to optimize the rib dimensions under the influence of multiple factors to achieve the best drag reduction effect.

[0003] Because wind turbine airfoils operate at high Reynolds numbers, typically in the millions, the size of ribs applied to them is extremely small, usually on the order of micrometers. This increases the difficulty of experimental testing. Furthermore, the small size of ribs leads to a sharp increase in the amount of mesh during numerical calculations, consuming a large amount of computational resources and hindering the practical engineering application of 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-reducing ribs on the surface of a wind turbine airfoil, a device for optimizing the size of drag-reducing ribs on the surface of a wind turbine airfoil, an electronic device, and a storage medium to overcome or at least partially solve the above problems.

[0005] To address the aforementioned problems, in a first aspect of this invention, an embodiment of the invention discloses a method for optimizing the size of drag-reducing ribs on the surface of a wind turbine airfoil, comprising:

[0006] Determine the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the wind turbine airfoil surface;

[0007] Based on the initial dimensional dimensions and the first friction drag coefficient, the dimensionless dimensions of the drag-reducing ribs located on the surface of the wind turbine airfoil are determined.

[0008] The first drag reduction ratio is determined based on the dimensionless dimensions.

[0009] The total drag reduction ratio of the drag-reducing ribs located on the surface of the wind turbine airfoil is determined by combining the first friction drag coefficient and the first drag reduction ratio.

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

[0011] Optionally, the step of determining the initial dimensional dimensions of the drag-reducing ribs and the first frictional drag coefficient of the airfoil surface includes:

[0012] The square root of the cross-sectional area of ​​the drag-reducing rib is calculated as the initial dimensional dimension;

[0013] Model the surface of the wind turbine airfoil to generate a wind turbine airfoil model;

[0014] Numerical simulation calculations were performed on the wind turbine wing model to determine the first friction drag coefficient.

[0015] Optionally, the step of determining the dimensionless dimension of the drag-reducing rib on the surface of the wind turbine airfoil based on the initial dimensional dimension and the first friction drag coefficient includes:

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

[0017]

[0018] lg is the initial dimensional dimension, lg + For dimensionless dimensions, u τ U is the frictional velocity of the airfoil surface, ν is the kinematic viscosity of air, and U ∞ For the incoming airflow velocity, C f,s This is the coefficient of first frictional resistance.

[0019] Optionally, the step of determining the first drag reduction ratio based on the dimensionless dimension includes:

[0020] Based on the dimensionless size, the drag reduction curve of the preset rib is looked up in a table to determine the drag reduction rate corresponding to the dimensionless size as the first drag reduction rate.

[0021] Optionally, the step of determining the total drag reduction ratio of the drag-reducing ribs located on the airfoil surface by combining the first friction drag coefficient and the first drag reduction ratio includes:

[0022] The second friction coefficient is determined based on the first friction coefficient and the first drag reduction ratio;

[0023] The viscous resistance is determined by integrating the second frictional resistance coefficient.

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

[0025] Optionally, the step of determining the second frictional resistance coefficient based on the first frictional resistance coefficient and the first drag reduction ratio includes:

[0026] Based on the formula for calculating the friction resistance coefficient, the second friction resistance coefficient is determined by combining the first friction resistance coefficient and the first drag reduction ratio. The formula for calculating the friction resistance coefficient is as follows:

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

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

[0029] Optionally, the step of integrating the second frictional resistance coefficient to determine the viscous resistance includes:

[0030] The viscous resistance is determined by integrating the second frictional resistance coefficient using an integral formula, wherein the integral formula is:

[0031]

[0032] C D,r It is a viscous resistance.

[0033] Optionally, the step of determining the total drag reduction rate based on the viscous resistance includes:

[0034] The total drag reduction rate is determined based on the 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 includes:

[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. The viscous drag formula is as follows:

[0037]

[0038] DR tota l represents the overall drag reduction ratio, and C... D The viscous drag of the wind turbine airfoil is given.

[0039] Optionally, the step of determining the target dimensional size corresponding to the maximum total drag reduction rate includes:

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

[0041] In a second aspect, embodiments of the present invention also disclose a device for optimizing the size of drag-reducing ribs on the surface of a wind turbine airfoil, comprising:

[0042] The first determining module is used to determine the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the wind turbine airfoil surface;

[0043] The first calculation module is used to determine the dimensionless size of the drag-reducing rib on the surface of the wind turbine airfoil based on the initial dimensional size and the first friction drag coefficient.

[0044] The second determining module is used to determine the first drag reduction ratio based on the dimensionless dimension.

[0045] The second calculation module is used to determine the total drag reduction rate of the drag reduction rib located on the surface of the wind turbine airfoil by combining the first friction drag coefficient and the first drag reduction rate.

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

[0047] Optionally, the first determining module includes:

[0048] The first calculation submodule is used to calculate the square root of the cross-sectional area of ​​the drag-reducing rib as the initial dimensional dimension;

[0049] The modeling submodule is used to model the surface of the wind turbine airfoil and generate a wind turbine airfoil model;

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

[0051] Optionally, the second computing module includes:

[0052] The fourth calculation submodule is used to determine the second friction resistance coefficient based on the first friction resistance coefficient and the first drag reduction ratio;

[0053] The fifth calculation submodule is used to perform integral calculations on the second frictional resistance coefficient to determine the viscous resistance;

[0054] The sixth calculation submodule is used to determine the total drag reduction rate based on the viscous resistance.

[0055] In a third aspect, embodiments of the present invention also disclose an electronic device, including 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 wind turbine airfoil surface drag reduction rib size optimization design method as described above.

[0056] In a fourth aspect, embodiments of the present invention also disclose a computer-readable storage medium storing a computer program that, when executed by a processor, implements the steps of the wind turbine airfoil surface drag reduction rib size optimization design method as described above.

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

[0058] This invention, through its embodiments, determines the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the wind turbine airfoil surface; determines the dimensionless dimensions of the drag-reducing ribs on the wind turbine airfoil surface based on the initial dimensional dimensions and the first friction drag coefficient; determines the first drag reduction ratio based on the dimensionless dimensions; determines the total drag reduction ratio of the drag-reducing ribs on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction ratio; determines the target dimensional dimension corresponding to the maximum total drag reduction ratio; and optimizes the drag-reducing rib dimensions by correlating the dimensional dimensions of the drag-reducing ribs with their drag reduction effect when applied to the wind turbine airfoil, thereby achieving rapid rib size design, optimizing the airfoil drag reduction effect, and facilitating the application of the ribs in practical engineering. Attached Figure Description

[0059] Figure 1 This is a flowchart illustrating the steps of an embodiment of the wind turbine airfoil surface drag reduction rib size optimization design method of the present invention;

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

[0061] Figure 3 This is an example of the airfoil surface friction drag coefficient distribution diagram of the wind turbine airfoil surface drag reduction rib size optimization design method of the present invention;

[0062] Figure 4 This is a schematic diagram of the rib cross-section as an example of a method for optimizing the size of drag-reducing ribs on the surface of a wind turbine airfoil according to the present invention;

[0063] Figure 5 This is a schematic diagram of the drag reduction curve of the ribs in an example of the wind turbine airfoil surface drag reduction rib size optimization design method of the present invention;

[0064] Figure 6This is a structural block diagram of an embodiment of the wind turbine airfoil surface drag reduction rib size optimization design device of the present invention;

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

[0066] Figure 8 This is a structural block diagram of a storage medium provided in an embodiment of the present invention. Detailed Implementation

[0067] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments.

[0068] Reference Figure 1 The diagram illustrates a step-by-step flowchart of an embodiment of a method for optimizing the size of drag-reducing ribs on the surface of a wind turbine airfoil according to the present invention. The method specifically includes the following steps:

[0069] Step 101: Determine the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the airfoil surface of the wind turbine;

[0070] In this embodiment of the invention, the original dimensional dimensions of the drag-reducing ribs, such as the height, width, and cross-sectional area of ​​the drag-reducing ribs, can be determined first, and these are defined as the initial dimensional dimensions. Simultaneously, the friction drag coefficient of the wind turbine airfoil surface, i.e., the first friction drag coefficient, can be determined under the condition that the wind turbine airfoil surface is smooth.

[0071] Step 102: Determine the dimensionless dimension of the drag-reducing rib located on the surface of the wind turbine airfoil based on the initial dimensional dimension and the first friction drag coefficient;

[0072] Since the initial dimensional dimensions have limitations such as units, they cannot be applied to all scenarios. In order to make accurate calculations, the initial dimensional dimensions and the first friction drag coefficient can be converted into dimensionless dimensions of the drag-reducing ribs located on the surface of the wind turbine airfoil, so that dimensionless parameters can be used to describe the characteristics of the deceleration ribs.

[0073] Step 103: Determine the first drag reduction ratio based on the dimensionless dimensions;

[0074] Based on the characteristics of the deceleration ribs corresponding to dimensionless dimensions, the drag reduction ratio of the deceleration ribs at various positions on the surface of the wind turbine airfoil is determined, i.e., the first drag reduction ratio.

[0075] Step 104: Determine the total drag reduction rate of the drag-reducing ribs located on the surface of the wind turbine airfoil by combining the first friction drag coefficient and the first drag reduction rate;

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

[0077] Step 105: Determine the target dimensional size corresponding to the maximum total drag reduction rate.

[0078] Based on the optimization range of rib size, the rib size of the wind turbine airfoil surface is optimized, and the target dimensional size corresponding to the maximum total drag reduction rate is determined, thereby obtaining 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 dimensional dimensions of the drag-reducing ribs and the first frictional drag coefficient of the airfoil surface includes:

[0080] Sub-step S1011: Calculate the square root of the cross-sectional area of ​​the drag-reducing rib as the initial dimensional dimension;

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

[0082] Sub-step S1012: Model the surface of the wind turbine airfoil to generate a wind turbine airfoil model;

[0083] It can model the surface of wind turbine airfoils, and generate wind turbine airfoil models through processes such as geometric shape modeling, unit mesh generation, laminate material layup, finite element calculation and post-processing analysis.

[0084] Sub-step S1013: Perform numerical simulation calculations on the wind turbine wing model to determine the first friction drag coefficient.

[0085] Then, numerical simulation calculations were performed on the wind turbine airfoil model to determine the friction drag coefficient of the smooth wind turbine airfoil surface, i.e., the first friction drag coefficient.

[0086] In addition, the first frictional resistance coefficient can also be calculated based on empirical formulas or experimental measurement methods, which are not limited in the embodiments of the present invention.

[0087] In an optional embodiment of the present invention, the step of determining the dimensionless size of the drag-reducing rib located on the surface of the wind turbine airfoil based on the initial dimensional size and the first friction drag coefficient includes: calculating the dimensionless size based on the initial dimensional size and the first friction drag coefficient using a dimensionless size formula; wherein, the dimensionless size formula is:

[0088]

[0089] lg is the initial dimensional dimension, lg +For dimensionless dimensions, u τ U is the frictional velocity of the airfoil surface, ν is the kinematic viscosity of air, and U ∞ For the incoming airflow velocity, C f,s This is the coefficient of first frictional resistance.

[0090] For the calculation of dimensionless dimensions, the initial dimensional dimension 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 table of the preset rib drag reduction curve based on the dimensionless size to determine the drag reduction rate corresponding to the dimensionless size as the first drag reduction rate.

[0092] The first drag reduction ratio of the deceleration rib depends on the dimensionless dimension lg of the rib. + Based on the dimensionless dimensions of the deceleration ribs at various locations on the airfoil surface, the drag reduction curve (DR-lg) of the preset ribs can be consulted. + The drag reduction rate corresponding to the dimensionless dimension is obtained as the first drag reduction rate. The preset rib drag reduction curve is determined through prior experiments. The process for determining the preset rib drag reduction curve can be determined according to actual conditions. This embodiment of the invention does not limit the process for determining the preset rib drag reduction curve.

[0093] In an optional embodiment of the present invention, the step of determining the total drag reduction ratio of the drag-reducing rib located on the surface of the wind turbine airfoil by combining the first frictional drag coefficient and the first drag reduction ratio includes:

[0094] Sub-step S1041: Determine the second friction resistance coefficient based on the first friction resistance coefficient and the first drag reduction ratio;

[0095] First, based on the first friction drag coefficient and the first drag reduction ratio, the friction drag coefficient of the search rib on the surface of the wind turbine airfoil, i.e., the second friction drag coefficient, can be determined.

[0096] Specifically, the step of determining the second frictional resistance coefficient based on the first frictional resistance coefficient and the first drag reduction ratio includes: determining the second frictional resistance coefficient based on the frictional resistance coefficient calculation formula, combined with the first frictional resistance coefficient and the first drag reduction ratio, wherein the frictional resistance coefficient calculation formula is:

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

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

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

[0100] Sub-step S1042: Integrate the second frictional resistance coefficient to determine the viscous resistance;

[0101] The viscous drag of the ribbed wind turbine airfoil is obtained by integrating the second friction drag coefficient.

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

[0103]

[0104] C D,r It is a viscous resistance.

[0105] Substitute the second frictional resistance coefficient into the integral formula, and then solve the integral formula to obtain the viscous resistance.

[0106] Sub-step S1043: Determine the total drag reduction rate based on the viscous resistance.

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

[0108] Specifically, the step of determining the total drag reduction rate based on 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] To calculate the total drag reduction rate, the viscous drag and the viscous drag corresponding to the 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 magnitude 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 includes: calculating 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:

[0111]

[0112] DR total For the total drag reduction ratio, C DThe viscous drag of the wind turbine airfoil is given.

[0113] In practical applications, the viscous resistance can be substituted into the viscous resistance formula for calculation, and the result 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] A preset target optimization algorithm can be used to set the optimization range for rib dimensions, thereby optimizing the rib dimensions on the airfoil surface of the wind turbine to obtain the rib dimensions with the best drag reduction effect. The input parameters of the preset target optimization algorithm are the first drag reduction ratio, the incoming wind speed, and the preset rib drag reduction curve. The independent variable is the dimensionless dimension, and the optimization target is the overall drag reduction ratio (DR) of the airfoil. total The maximum is obtained by solving the optimization algorithm based on the preset target, and the corresponding dimensional size is the target dimensional size.

[0116] This invention, through its embodiments, determines the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the wind turbine airfoil surface; determines the dimensionless dimensions of the drag-reducing ribs on the wind turbine airfoil surface based on the initial dimensional dimensions and the first friction drag coefficient; determines the first drag reduction ratio based on the dimensionless dimensions; determines the total drag reduction ratio of the drag-reducing ribs on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction ratio; determines the target dimensional dimension corresponding to the maximum total drag reduction ratio; and optimizes the drag-reducing rib dimensions by correlating the dimensional dimensions of the drag-reducing ribs with their drag reduction effect when applied to the wind turbine airfoil, thereby achieving rapid rib size design, optimizing the airfoil drag reduction effect, and facilitating the application of the ribs in practical engineering.

[0117] To enable those skilled in the art to clearly understand the implementation process of the embodiments of the present invention, please refer to... Figure 2 Let's illustrate with an example:

[0118] The DU-96-W-18 airfoil was selected for the wind turbine. The airfoil was modeled, and the surface friction drag coefficient C was obtained through numerical simulation. f,s The numerical model used for distribution calculation is the Transition-SST model, with an airfoil chord length of c = 1m, an airfoil angle of attack of α = 0°, and an incoming wind speed of U. ∞ The wind speed was 32.5 m / s (wind tunnel test wind speed). The calculated airfoil surface friction drag coefficient is shown in the attached figure. Figure 3 As shown.

[0119] Choose a ribbed structure with a trapezoidal cross-section, as shown in the attached diagram. Figure 4 As shown, an initial rib size lg is set, and the dimensionless size lg of the rib structure at various positions on the airfoil surface is calculated using formula (1). + Then, using formula (4), the overall drag reduction effect DR of the rib structure on the wind turbine airfoil is obtained. total .

[0120] An adaptive simulated annealing algorithm was used as the optimization design method, with the rib optimization range set to lg = 10–200 μm. The input parameter for the optimization design method was the surface friction drag coefficient C of the wind turbine airfoil. f,s Incoming wind speed U ∞ DR-lg of rib drag reduction curve + See attached Figure 5 The independent variable is the dimensionless size represented by lg, and the optimization objective is the overall drag reduction ratio DR of the airfoil. total maximum.

[0121] Two optimization cases were calculated: a case with full rib coverage and a case with rib coverage of the 0-0.5c airfoil portion. 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] As can be seen, for Case 1, the optimal rib size lg = 78 μm has a total drag reduction rate of 5.01%, and for Case 2, the optimal rib size lg = 74 μm has a total drag reduction rate of 3.69%. Both have achieved significant drag reduction effects, proving that the method of this patent can realize the optimized design of drag reduction rib size for wind turbine airfoil surfaces.

[0127] It should be noted that, for the sake of simplicity, the method embodiments are all described as a series of actions. However, those skilled in the art should understand that the embodiments of the present invention are not limited to the described order of actions, because according to the embodiments of the present invention, some steps can be performed in other orders or simultaneously. Furthermore, those skilled in the art should also understand that the embodiments described in the specification are preferred embodiments, and the actions involved are not necessarily essential to the embodiments of the present invention.

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

[0129] The first determining module 601 is used to determine the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the airfoil surface of the wind turbine.

[0130] The first calculation module 602 is used to determine the dimensionless size of the drag-reducing rib on the surface of the wind turbine airfoil based on the initial dimensional size and the first friction drag coefficient.

[0131] The second determining module 603 is used to determine the first drag reduction ratio based on the dimensionless dimension.

[0132] The second calculation module 604 is used to determine the total drag reduction rate of the drag reduction rib located on the surface of the wind turbine airfoil 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 dimensional 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] The first calculation submodule is used to calculate the square root of the cross-sectional area of ​​the drag-reducing rib as the initial dimensional dimension;

[0136] The modeling submodule is used to model the surface of the wind turbine airfoil and generate a wind turbine airfoil model;

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

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

[0139] The third calculation submodule is used to calculate the initial dimensional dimension and the first frictional resistance coefficient based on the dimensionless dimension formula to determine the dimensionless dimension; wherein, the dimensionless dimension formula is:

[0140]

[0141] lg is the initial dimensional dimension, lg + For dimensionless dimensions, u τ U is the frictional velocity of the airfoil surface, ν is the kinematic viscosity of air, and U ∞ For the incoming airflow velocity, C f,s This is the coefficient of first frictional resistance.

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

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

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

[0145] The fourth calculation submodule is used to determine the second friction resistance coefficient based on the first friction resistance coefficient and the first drag reduction ratio;

[0146] The fifth calculation submodule is used to perform integral calculations on the second frictional resistance coefficient to determine the viscous resistance;

[0147] The sixth calculation submodule is used to determine the total drag reduction rate based on the viscous resistance.

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

[0149] The friction calculation unit is used to determine the second friction resistance coefficient based on the friction resistance coefficient calculation formula, combined with the first friction resistance coefficient and the first drag reduction ratio. The friction resistance coefficient calculation formula is as follows:

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

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

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

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

[0154]

[0155] C D,r It is a viscous resistance.

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

[0157] The total drag reduction 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 calculation unit includes:

[0159] The total drag reduction 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. The viscous drag formula is as follows:

[0160]

[0161] DR total For the total drag reduction ratio, C D The viscous drag of the wind turbine airfoil is given.

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

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

[0164] This invention, through its embodiments, determines the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the wind turbine airfoil surface; determines the dimensionless dimensions of the drag-reducing ribs on the wind turbine airfoil surface based on the initial dimensional dimensions and the first friction drag coefficient; determines the first drag reduction ratio based on the dimensionless dimensions; determines the total drag reduction ratio of the drag-reducing ribs on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction ratio; determines the target dimensional dimension corresponding to the maximum total drag reduction ratio; and optimizes the drag-reducing rib dimensions by correlating the dimensional dimensions of the drag-reducing ribs with their drag reduction effect when applied to the wind turbine airfoil, thereby achieving rapid rib size design, optimizing the airfoil drag reduction effect, and facilitating the application of the ribs in practical engineering.

[0165] As the device embodiment is basically similar to the method embodiment, the description is relatively simple, and relevant parts can be found in the description of the method embodiment.

[0166] Reference Figure 7 This invention also provides an electronic device, comprising:

[0167] The device includes a processor 701 and a storage medium 702, wherein the storage medium 702 stores a computer program executable by the processor 701. When the electronic device is running, the processor 701 executes the computer program to perform the wind turbine airfoil surface drag reduction rib size optimization design method as described in any one of the embodiments of the present invention.

[0168] The method for optimizing the size of drag-reducing ribs on the surface of the wind turbine airfoil includes:

[0169] Determine the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the wind turbine airfoil surface;

[0170] The dimensionless dimensions of the drag-reducing ribs located on the surface of the wind turbine airfoil are determined based on the initial dimensional dimensions and the first friction drag coefficient.

[0171] The first drag reduction ratio is determined based on the dimensionless dimensions.

[0172] The total drag reduction ratio of the drag-reducing ribs located on the surface of the wind turbine airfoil is determined by combining the first friction drag coefficient and the first drag reduction ratio.

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

[0174] Optionally, the step of determining the initial dimensional dimensions of the drag-reducing ribs and the first frictional drag coefficient of the airfoil surface includes:

[0175] The square root of the cross-sectional area of ​​the drag-reducing rib is calculated as the initial dimensional dimension;

[0176] Model the surface of the wind turbine airfoil to generate a wind turbine airfoil model;

[0177] Numerical simulation calculations were performed on the wind turbine wing model to determine the first friction drag coefficient.

[0178] Optionally, the step of determining the dimensionless size of the drag-reducing rib on the surface of the wind turbine airfoil based on the initial dimensional size and the first friction drag coefficient includes:

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

[0180]

[0181] lg is the initial dimensional dimension, lg + For dimensionless dimensions, u τ U is the frictional velocity of the airfoil surface, ν is the kinematic viscosity of air, and U ∞ For the incoming airflow velocity, C f,s This is the coefficient of first frictional resistance.

[0182] Optionally, the step of determining the first drag reduction ratio based on the dimensionless dimension includes:

[0183] Based on the dimensionless size, the drag reduction curve of the preset rib is looked up in a table to determine the drag reduction rate corresponding to the dimensionless size as the first drag reduction rate.

[0184] Optionally, the step of determining the total drag reduction ratio of the drag-reducing ribs located on the airfoil surface by combining the first friction drag coefficient and the first drag reduction ratio includes:

[0185] The second friction coefficient is determined based on the first friction coefficient and the first drag reduction ratio;

[0186] The viscous resistance is determined by integrating the second frictional resistance coefficient.

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

[0188] Optionally, the step of determining the second frictional resistance coefficient based on the first frictional resistance coefficient and the first drag reduction ratio includes:

[0189] Based on the formula for calculating the friction resistance coefficient, the second friction resistance coefficient is determined by combining the first friction resistance coefficient and the first drag reduction ratio. The formula for calculating the friction resistance coefficient is as follows:

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

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

[0192] Optionally, the step of integrating the second frictional resistance coefficient to determine the viscous resistance includes:

[0193] The viscous resistance is determined by integrating the second frictional resistance coefficient using an integral formula, wherein the integral formula is:

[0194]

[0195] C D,r It is a viscous resistance.

[0196] Optionally, the step of determining the total drag reduction rate based on the viscous resistance includes:

[0197] The total drag reduction rate is determined based on the 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 includes:

[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. The viscous drag formula is as follows:

[0200]

[0201] DR total For the total drag reduction ratio, C D The viscous drag of the wind turbine airfoil is given.

[0202] Optionally, the step of determining the target dimensional size corresponding to the maximum total drag reduction rate includes:

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

[0204] This invention, through its embodiments, determines the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the wind turbine airfoil surface; determines the dimensionless dimensions of the drag-reducing ribs on the wind turbine airfoil surface based on the initial dimensional dimensions and the first friction drag coefficient; determines the first drag reduction ratio based on the dimensionless dimensions; determines the total drag reduction ratio of the drag-reducing ribs on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction ratio; determines the target dimensional dimension corresponding to the maximum total drag reduction ratio; and optimizes the drag-reducing rib dimensions by correlating the dimensional dimensions of the drag-reducing ribs with their drag reduction effect when applied to the wind turbine airfoil, thereby achieving rapid rib size design, optimizing the airfoil drag reduction effect, and facilitating the application of the ribs in practical engineering.

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

[0206] The processors mentioned above can be general-purpose processors, including central processing units (CPUs), network processors (NPs), etc.; they can also be digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.

[0207] Reference Figure 8The present invention also provides a computer-readable storage medium 801, on which a computer program is stored. When the computer program is run by a processor, it executes the wind turbine airfoil surface drag reduction rib size optimization design method as described in any one of the embodiments of the present invention.

[0208] The method for optimizing the size of drag-reducing ribs on the surface of the wind turbine airfoil includes:

[0209] Determine the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the wind turbine airfoil surface;

[0210] The dimensionless dimensions of the drag-reducing ribs located on the surface of the wind turbine airfoil are determined based on the initial dimensional dimensions and the first friction drag coefficient.

[0211] The first drag reduction ratio is determined based on the dimensionless dimensions.

[0212] The total drag reduction ratio of the drag-reducing ribs located on the surface of the wind turbine airfoil is determined by combining the first friction drag coefficient and the first drag reduction ratio.

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

[0214] Optionally, the step of determining the initial dimensional dimensions of the drag-reducing ribs and the first frictional drag coefficient of the airfoil surface includes:

[0215] The square root of the cross-sectional area of ​​the drag-reducing rib is calculated as the initial dimensional dimension;

[0216] Model the surface of the wind turbine airfoil to generate a wind turbine airfoil model;

[0217] Numerical simulation calculations were performed on the wind turbine wing model to determine the first friction drag coefficient.

[0218] Optionally, the step of determining the dimensionless size of the drag-reducing rib on the surface of the wind turbine airfoil based on the initial dimensional size and the first friction drag coefficient includes:

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

[0220]

[0221] lg is the initial dimensional dimension, lg + For dimensionless dimensions, u τ U is the frictional velocity of the airfoil surface, ν is the kinematic viscosity of air, and U ∞ For the incoming airflow velocity, C f,sThis is the coefficient of first frictional resistance.

[0222] Optionally, the step of determining the first drag reduction ratio based on the dimensionless dimension includes:

[0223] Based on the dimensionless size, the drag reduction curve of the preset rib is looked up in a table to determine the drag reduction rate corresponding to the dimensionless size as the first drag reduction rate.

[0224] Optionally, the step of determining the total drag reduction ratio of the drag-reducing ribs located on the airfoil surface by combining the first friction drag coefficient and the first drag reduction ratio includes:

[0225] The second friction coefficient is determined based on the first friction coefficient and the first drag reduction ratio;

[0226] The viscous resistance is determined by integrating the second frictional resistance coefficient.

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

[0228] Optionally, the step of determining the second frictional resistance coefficient based on the first frictional resistance coefficient and the first drag reduction ratio includes:

[0229] Based on the formula for calculating the friction resistance coefficient, the second friction resistance coefficient is determined by combining the first friction resistance coefficient and the first drag reduction ratio. The formula for calculating the friction resistance coefficient is as follows:

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

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

[0232] Optionally, the step of integrating the second frictional resistance coefficient to determine the viscous resistance includes:

[0233] The viscous resistance is determined by integrating the second frictional resistance coefficient using an integral formula, wherein the integral formula is:

[0234]

[0235] C D,r It is a viscous resistance.

[0236] Optionally, the step of determining the total drag reduction rate based on the viscous resistance includes:

[0237] The total drag reduction rate is determined based on the 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 includes:

[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. The viscous drag formula is as follows:

[0240]

[0241] DR total For the total drag reduction ratio, C D The viscous drag of the wind turbine airfoil is given.

[0242] Optionally, the step of determining the target dimensional size corresponding to the maximum total drag reduction rate includes:

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

[0244] This invention, through its embodiments, determines the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the wind turbine airfoil surface; determines the dimensionless dimensions of the drag-reducing ribs on the wind turbine airfoil surface based on the initial dimensional dimensions and the first friction drag coefficient; determines the first drag reduction ratio based on the dimensionless dimensions; determines the total drag reduction ratio of the drag-reducing ribs on the wind turbine airfoil surface by combining the first friction drag coefficient and the first drag reduction ratio; determines the target dimensional dimension corresponding to the maximum total drag reduction ratio; and optimizes the drag-reducing rib dimensions by correlating the dimensional dimensions of the drag-reducing ribs with their drag reduction effect when applied to the wind turbine airfoil, thereby achieving rapid rib size design, optimizing the airfoil drag reduction effect, and facilitating the application of the ribs in practical engineering.

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

[0246] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, apparatus, or computer program products. Therefore, embodiments of the present invention can take the form of entirely hardware embodiments, entirely software embodiments, or embodiments combining software and hardware aspects. Furthermore, embodiments of the present invention can take the form of computer program products 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 code.

[0247] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, terminal devices (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing terminal device to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing terminal device, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0248] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing terminal device to operate in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The 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 equipment, causing a series of operational steps to be performed on the computer or other programmable terminal equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable terminal equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0250] Although preferred embodiments of the present invention have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments as well as all changes and modifications falling within the scope of the embodiments of the present invention.

[0251] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or terminal device that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or terminal device. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or terminal device that includes said element.

[0252] The foregoing has provided a detailed description of the wind turbine airfoil surface drag reduction rib size optimization design method, the wind turbine airfoil surface drag reduction rib size optimization design device, the electronic device, and the storage medium provided by the present invention. Specific examples have been used to illustrate the principles and implementation methods of the present invention. The above description of the embodiments is only for the purpose of helping to understand the method and core ideas of the present invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of the present invention. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for optimizing the design of drag-reducing rib dimensions on the surface of a wind turbine airfoil, characterized in that, include: Determine the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the wind turbine airfoil surface; Based on the initial dimensional dimensions and the first friction drag coefficient, the dimensionless dimensions of the drag-reducing ribs located on the surface of the wind turbine airfoil are determined. The first drag reduction ratio is determined based on the dimensionless dimensions. The total drag reduction ratio of the drag-reducing ribs located on the surface of the wind turbine airfoil is determined by combining the first friction drag coefficient and the first drag reduction ratio. Determine the target dimensionless size corresponding to the maximum total drag reduction rate; The step of determining the initial dimensional dimensions of the drag-reducing ribs and the first frictional drag coefficient of the airfoil surface includes: The square root of the cross-sectional area of ​​the drag-reducing rib is calculated as the initial dimensional dimension; Model the surface of the wind turbine airfoil to generate a wind turbine airfoil model; Numerical simulation calculations were performed on the wind turbine wing model to determine the first friction drag coefficient; The step of determining the dimensionless dimension of the drag-reducing rib on the surface of the wind turbine airfoil based on the initial dimensional dimension and the first friction drag coefficient includes: Based on the dimensionless dimension formula, the initial dimensional dimension and the first frictional resistance coefficient are calculated to determine the dimensionless dimension; wherein, the dimensionless dimension formula is: ， For the initial dimensional dimensions, Dimensionless dimensions The friction speed of the airfoil surface. The viscosity of air motion. For the incoming wind speed, The coefficient of friction is the first frictional resistance. The step of determining the total drag reduction ratio of the drag-reducing rib located on the surface of the wind turbine airfoil by combining the first friction drag coefficient and the first drag reduction ratio includes: The second friction coefficient is determined based on the first friction coefficient and the first drag reduction ratio; The viscous resistance is determined by integrating the second frictional resistance coefficient. The total drag reduction rate is determined based on the viscous resistance.

2. The method of claim 1, wherein, The step of determining the first drag reduction ratio based on the dimensionless dimension includes: Based on the dimensionless size, the drag reduction curve of the preset rib is looked up in a table to determine the drag reduction rate corresponding to the dimensionless size as the first drag reduction rate.

3. The method of claim 1, wherein, The step of determining the second friction resistance coefficient based on the first friction resistance coefficient and the first drag reduction ratio includes: Based on the formula for calculating the friction resistance coefficient, the second friction resistance coefficient is determined by combining the first friction resistance coefficient and the first drag reduction ratio. The formula for calculating the friction resistance coefficient is as follows: ， a second frictional resistance coefficient, is a first frictional resistance coefficient, is a first drag reduction rate.

4. The method of claim 3, wherein, The step of integrating the second frictional resistance coefficient to determine the viscous resistance includes: The viscous resistance is determined by integrating the second frictional resistance coefficient using an integral formula, wherein the integral formula is: ， is the viscous drag.

5. The method of claim 4, wherein, The step of determining the total drag reduction rate based on the viscous resistance includes: The total drag reduction rate is determined based on the difference between the viscous drag and the viscous drag of the wind turbine airfoil.

6. The method of claim 4, wherein, 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 includes: 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. The viscous drag formula is as follows: ， is the total drag reduction rate, is the viscous drag of the wind turbine airfoil.

7. The method of claim 1, wherein, The step of determining the target dimensionless size corresponding to the maximum total drag reduction rate includes: A preset target optimization algorithm is adopted, with the maximum total drag reduction rate as the optimization target of the preset target optimization algorithm, and the target dimension is determined.

8. A device for optimizing the size of drag-reducing ribs on the surface of a wind turbine airfoil, characterized in that, include: The first determining module is used to determine the initial dimensional dimensions of the drag-reducing ribs and the first friction drag coefficient of the wind turbine airfoil surface; The first calculation module is used to determine the dimensionless dimension of the drag-reducing rib located on the surface of the wind turbine airfoil based on the initial dimensional dimension and the first friction drag coefficient. The second determining module is used to determine the first drag reduction ratio based on the dimensionless dimension. The second calculation module is used to determine the total drag reduction rate of the drag reduction rib located on the surface of the wind turbine airfoil by combining the first friction drag coefficient and the first drag reduction rate. The third calculation module is used to determine the target dimensional size corresponding to the maximum total drag reduction rate. The first determining module includes: The first calculation submodule is used to calculate the square root of the cross-sectional area of ​​the drag-reducing rib as the initial dimensional dimension; The modeling submodule is used to model the surface of the wind turbine airfoil and generate a wind turbine airfoil model; The second calculation submodule is used to perform numerical simulation calculations on the wind turbine wing model to determine the first friction drag coefficient; The first computing module includes: The third calculation submodule is used to calculate the initial dimensional dimension and the first frictional resistance coefficient based on the dimensionless dimension formula to determine the dimensionless dimension; wherein, the dimensionless dimension formula is: ， For the initial dimensional dimensions, Dimensionless dimensions The friction speed of the airfoil surface. The viscosity of air motion. For the incoming wind speed, The first frictional resistance coefficient The second calculation module includes: The fourth calculation submodule is used to determine the second friction resistance coefficient based on the first friction resistance coefficient and the first drag reduction ratio; The fifth calculation submodule is used to perform integral calculations on the second frictional resistance coefficient to determine the viscous resistance; The sixth calculation submodule is used to determine the total drag reduction rate based on the viscous resistance.

9. An electronic device, comprising: The device includes a processor, a memory, and a computer program stored in the memory and capable of running on the processor. When executed by the processor, the computer program implements the steps of the wind turbine airfoil surface drag reduction rib size optimization design method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, A computer program is stored on the computer-readable storage medium, which, when executed by a processor, implements the steps of the wind turbine airfoil surface drag reduction rib size optimization design method as described in any one of claims 1 to 7.

Images