Flow field simulation test method and device for wind turbine

Abstract

The present invention relates to the field of new energy technologies, and provides a method and device for simulating and testing the flow field of a wind turbine. The method includes: based on a three-dimensional flow field model of the wind turbine, performing transient fluid dynamics simulation on the wind turbine at the current moment to obtain simulation data of the wind turbine and the three-dimensional flow field model at the current moment; based on the simulation data and the control strategy of the wind turbine, obtaining the simulation control parameters of the wind turbine at the next moment; based on the simulation control parameters, updating the rotational speed of the reference coordinate system in the revolution area and the rotational speed of the grid nodes in the rotation area in the three-dimensional flow field model. The method and device for simulating and testing the flow field of the wind turbine provided by the present invention can perform unsteady computational fluid dynamics analysis on the process of wind turbine variable speed and blade pitch change under the operating state of the wind turbine. The simulation takes into account the variable speed and pitch change actions of the wind turbine, and the obtained transient flow field calculation results are closer to the actual situation, and more accurate simulation test results can be obtained.

Description

Technical Field

[0001] The present invention relates to the field of new energy technology, and in particular to a wind turbine flow field simulation test method and device. Background Art

[0002] As a way to utilize clean and green energy, wind power generation is increasingly gaining attention and importance. The simulation test of wind turbines is of great significance to the efficiency, reliability, correctness and uniformity of wind turbine design and implementation.

[0003] In the prior art, wind turbines can be simulated and tested based on computational fluid dynamics (CFD) methods. However, the existing wind turbine flow field simulation test methods are difficult to accurately simulate the flow field changes of wind turbines during the dynamic changes of the wind rotor speed and blade pitch angle. Therefore, how to accurately simulate the flow field changes of wind turbines during the dynamic changes of the wind rotor speed and blade pitch angle is a technical problem that needs to be solved in this field. Summary of the invention

[0004] The present invention provides a wind turbine flow field simulation test method and device, which are used to solve the problem in the prior art that it is difficult to accurately simulate the flow field changes of the wind turbine during the dynamic changes of the wind rotor speed and the blade pitch angle, and to achieve more accurate simulation of the flow field changes of the wind turbine during the dynamic changes of the wind rotor speed and the blade pitch angle.

[0005] The present invention provides a wind turbine flow field simulation test method, comprising:

[0006] Constructing and based on a three-dimensional flow field model of a wind turbine, performing a transient fluid dynamics simulation on the wind turbine at the current moment, and obtaining simulation data of the wind turbine and the three-dimensional flow field model at the current moment;

[0007] Based on the simulation data and the control strategy of the wind turbine, obtaining simulation control parameters of the wind turbine at the next moment;

[0008] Based on the simulation control parameters, updating the reference coordinate system rotation speed of the revolution zone in the three-dimensional flow field model and the grid node rotation speed of the rotation zone in the three-dimensional flow field model;

[0009] Among them, the three-dimensional flow field model includes: a rotation area, a revolution area, and a stationary area; the rotation area extends circumferentially along the pitch axis direction of the blade in the wind turbine, and the radius of the rotation area is greater than the maximum distance between the blade surface and the pitch axis of the blade; the revolution area and the stationary area extend circumferentially along the direction of the wind turbine's rotor rotation axis, the radius of the revolution area is greater than the length of the blade, and the stationary area is nested outside the revolution area.

[0010] According to a wind turbine flow field simulation test method provided by the present invention, the revolution area in the three-dimensional flow field model is established based on a multiple reference frame model, the spatial coordinates of the grid nodes in the revolution area are fixed, and the multiple reference frame model is used to add a relative rotational speed to the revolution area.

[0011] According to a wind turbine flow field simulation test method provided by the present invention, the rotation area in the three-dimensional flow field model is established based on a sliding mesh model, the grid nodes can rotate around the central axis of the rotation area, and the sliding mesh model is used to change the positions of the grid nodes and the boundary of the three-dimensional flow field model to simulate the flow field change during blade pitching.

[0012] According to a wind turbine flow field simulation test method provided by the present invention, in the three-dimensional flow field model, boundary layer grids are arranged on the surface of the blade, and the thickness of the first layer of the boundary layer grid meets a preset condition; the revolution area and the stationary area adopt a contact method of sharing grid nodes; data transfer is realized between the rotation area and the revolution area through an interface of non-shared nodes; the grid size of the revolution area is greater than the grid size of the rotation area.

[0013] According to a wind turbine flow field simulation test method provided by the present invention, the rotation area is a cylinder; the central axis of the rotation area coincides with the pitch axis of the blade; the distance between one end of the rotation area and the rotation center of the wind is a preset value; the radius of the rotation area is determined based on the maximum distance between the blade surface and the pitch axis of the blade; the length of the rotation area is determined based on the length of the blade.

[0014] According to a wind turbine flow field simulation test method provided by the present invention, the revolution area is a cylinder, the central axis of the revolution area coincides with the rotation axis of the rotor, and the rotation center of the rotor coincides with the midpoint of the central axis of the revolution area; the radius of the revolution area is determined based on the length of the rotation area; the length of the revolution area is determined based on the radius of the rotation area.

[0015] A method for simulating and testing the flow field of a wind turbine provided by the present invention, wherein the static region is an annular cylinder; the central axis of the static region coincides with the rotation axis of the wind wheel; the outer diameter of the static region is determined based on the radius of the revolution region; the distance between the inlet boundary of the static region and the wind wheel group in the wind turbine, and the distance between the outlet boundary of the static region and the wind wheel group are determined based on the radius of the rotation region.

[0016] A method for simulating and testing the flow field of a wind turbine provided by the present invention, wherein the three-dimensional flow field model is composed of the three-dimensional flow field models of each blade; the boundary of the three-dimensional flow field model of any blade is a periodic symmetric boundary.

[0017] The present invention also provides a wind turbine simulation test device, including:

[0018] A data acquisition module, configured to construct and perform transient fluid dynamics simulation on the wind turbine at the current moment based on the three-dimensional flow field model of the wind turbine, and acquire the simulation data of the wind turbine and the three-dimensional flow field model at the current moment;

[0019] A data calculation module, configured to acquire the simulation control parameters of the wind turbine at the next moment based on the simulation data and the control strategy of the wind turbine;

[0020] A simulation control module, configured to update the rotational speed of the reference coordinate system in the revolution region and the rotational speed of the grid nodes in the rotation region of the three-dimensional flow field model based on the simulation control parameters;

[0021] Wherein, the three-dimensional flow field model includes: a rotation region, a revolution region and a static region; the rotation region extends circumferentially along the pitch axis direction of the blades in the wind turbine, and the radius of the rotation region is greater than the maximum distance between the blade surface and the pitch axis of the blade; the revolution region and the static region extend circumferentially along the rotation axis direction of the wind wheel in the wind turbine, the radius of the revolution region is greater than the length of the blade, and the static region is nested outside the revolution region.

[0022] The present invention also provides an electronic device, including a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein when the processor executes the program, the method for simulating and testing the flow field of the wind turbine as described in any one of the above is implemented.

[0023] The present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the method for simulating and testing the flow field of the wind turbine as described in any one of the above is implemented.

[0024] The wind turbine flow field simulation test method and device provided by the present invention perform transient fluid dynamics simulation on the wind turbine at the current moment based on the three-dimensional flow field model of the wind turbine. After obtaining the simulation data of the wind turbine and the three-dimensional flow field model at the current moment, based on the above simulation data and the control strategy of the wind turbine, the simulation control parameters of the wind turbine at the next moment are obtained. Based on the above simulation control parameters, the rotational speed of the reference coordinate system in the revolution area and the rotational speed of the grid nodes in the rotation area of the three-dimensional flow field model are updated. It can perform unsteady computational fluid dynamics analysis on the wind turbine speed change and blade pitch change processes during the operation of the wind turbine. The variable speed and variable pitch actions of the wind turbine are considered in the simulation, and the obtained transient flow field calculation results are closer to the actual situation, enabling more accurate simulation test results, which is of great significance for improving the efficiency, reliability, correctness, and unity of the design and implementation of wind turbines. When performing steady-state simulation calculations on the wind turbine flow field at different blade pitch angles, by controlling the rotation of the rotation area to make the blade at the target blade pitch angle, the time cost of the preprocessing of the flow field calculation can be effectively reduced. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] In order to more clearly illustrate the technical solutions in the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the following drawings are some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0026] Figure 1 is one of the flow schematic diagrams of the wind turbine flow field simulation test method provided by the present invention;

[0027] Figure 2 is the sectional view of the three-dimensional flow field model of the blade in the wind turbine flow field simulation test method provided by the present invention;

[0028] Figure 3 is the front view of the three-dimensional flow field model of the blade in the wind turbine flow field simulation test method provided by the present invention;

[0029] Figure 4 is the left view of the three-dimensional flow field model of the blade in the wind turbine flow field simulation test method provided by the present invention;

[0030] Figure 5 is the second flow schematic diagram of the wind turbine flow field simulation test method provided by the present invention;

[0031] Figure 6 is the structural schematic diagram of the wind turbine simulation test device provided by the present invention;

[0032] Figure 7 It is a schematic structural diagram of the electronic device provided by the present invention. Specific embodiments

[0033] To make the objectives, technical solutions and advantages of the present invention clearer, the technical solutions in the present invention will be clearly and completely described below with reference to the accompanying drawings in the present invention. Apparently, the described embodiments are some but not all of the embodiments of the present invention. All other embodiments obtained by those of ordinary skill in the art based on the embodiments in the present invention without making creative efforts shall fall within the protection scope of the present invention.

[0034] In the description of the present invention, it should be noted that unless otherwise clearly defined and limited, the terms "installed", "connected" and "connected" should be understood in a broad sense. For example, it may be a fixed connection, a detachable connection or an integral connection; it may be a mechanical connection or an electrical connection; it may be directly connected or indirectly connected through an intermediate medium, and it may be the internal connection of two components. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood according to specific situations.

[0035] It should be noted that when a wind turbine is simulated and tested based on the traditional wind turbine flow field simulation test method, the flow field where the wind turbine is located can be simulated through steady-state calculation under the condition that the wind wheel speed and the blade pitch angle in the wind turbine are constant. However, the traditional wind turbine flow field simulation test method is difficult to accurately simulate the flow field change of the wind turbine during the dynamic change process of the above wind wheel speed and blade pitch angle, and it is difficult to reasonably predict the aerodynamic performance of the wind wheel. In addition, in the steady-state analysis of the traditional wind turbine flow field simulation test method, if the flow field under multiple pitch angles needs to be calculated, it is necessary to readjust the blade pitch angle and then perform grid division on the new flow field, resulting in a relatively high time cost for preprocessing.

[0036] In view of this, the present invention provides a wind turbine flow field simulation test method and device. Based on the wind turbine flow field simulation test method provided by the present invention, the flow field change of the wind turbine can be simulated during the dynamic change process of the wind wheel speed and / or the blade pitch angle in the wind turbine, and thus the simulation test result of the wind turbine can be obtained more accurately, which is of great significance for improving the efficiency, reliability, correctness and unity of the aerodynamic design and implementation of the wind turbine.

[0037] Figure 1 It is one of the flow schematic diagrams of the wind turbine flow field simulation test method provided by the present invention. The following combines Figure 1 to describe the wind turbine flow field simulation test method of the present invention. As Figure 1As shown in the figure, the method includes: Step 101, constructing and based on a three-dimensional flow field model of a wind turbine, performing transient fluid dynamics simulation on the wind turbine at the current moment, and obtaining the simulation data of the wind turbine at the current moment.

[0038] Among them, the three-dimensional flow field model includes: a self-rotation area, a revolution area, and a stationary area; the self-rotation area extends circumferentially along the pitch axis direction of the blade in the wind turbine, and the radius of the self-rotation area is greater than the maximum distance between the blade surface and the pitch axis of the blade; the revolution area and the stationary area extend circumferentially along the direction of the wind turbine's rotor rotation axis, and the stationary area is nested outside the revolution area; the radius of the revolution area is greater than the length of the blade.

[0039] The simulation data includes: the wind speed distribution, the aerodynamic torque of the rotor, the rotor speed, the blade pitch angle, and the motor torque in the above three-dimensional flow field model.

[0040] The simulation control parameters include: the rotor speed, the pitch speed of the blade, the rotation speed of the reference coordinate system in the revolution area, and the rotation speed of the grid nodes in the self-rotation area. Among them, the positive and negative of the above rotation speed can represent the rotation direction of the above grid nodes.

[0041] Specifically, the computational fluid dynamics method is the product of the combination of fluid mechanics, numerical mathematics, and computer science. The computational fluid dynamics method approximately represents the integral and differential terms in the control equations of fluid mechanics as discrete algebraic forms, making them into algebraic equation systems, and then by solving these discrete algebraic equation systems through a computer, numerical solutions at discrete time / space points can be obtained.

[0042] Figure 2 It is the cross-sectional view of the three-dimensional flow field model of the blade in the wind turbine flow field simulation test method provided by the present invention. Figure 3 It is the front view of the three-dimensional flow field model of the blade in the wind turbine flow field simulation test method provided by the present invention. Figure 4 It is the left view of the three-dimensional flow field model of the blade in the wind turbine flow field simulation test method provided by the present invention.

[0043] Based on the content of the above embodiments, the three-dimensional flow field model is composed of the three-dimensional flow field models of each blade; the boundary of the three-dimensional flow field model of any blade is a periodic symmetric boundary.

[0044] It should be noted that Figures 2 to 4 As shown in the figure, it is the three-dimensional flow field model of any blade 1 in the wind turbine. The three-dimensional flow field model of the above blade 1 includes a self-rotation area 2, a revolution area 3, and a stationary area 4. Combining the three-dimensional flow field models of each blade 1 can obtain the three-dimensional flow field model of the wind turbine.

[0045] A three-dimensional flow field model of a wind turbine can be constructed based on the structural parameters of the wind turbine. Turbulent flow field modeling can be carried out based on the k-omega SST model or the Transition SST model, pressure-velocity coupling can be performed based on the Coupled algorithm, and a numerical format that satisfies second-order accuracy can be used for the spatial difference of the convection term. Among them, the k-omega SST model and the Transition SST model are widely used turbulent models.

[0046] It should be noted that due to the periodicity and symmetry of the blade distribution in the wind turbine, when creating a three-dimensional flow field model of the wind turbine, a three-dimensional flow field model of any one blade 1 in the wind turbine can be created only, and the periodicity and symmetry of the distribution of blade 1 can be utilized to obtain the three-dimensional flow field model of the wind turbine. Among them, in the three-dimensional flow field model of blade 1, the two sides of the flow field are periodic symmetric boundaries, and the periodic symmetric boundaries can reduce the number of grids in the above three-dimensional flow field model to one-third of the three-dimensional flow field model of the wind wheel in the entire wind turbine.

[0047] The inlet of the three-dimensional flow field model of the wind turbine can adopt a velocity boundary condition, and the outlet can adopt a pressure boundary condition.

[0048] Each blade 1 in the wind turbine corresponds to a self-rotation area 2. The self-rotation area 2 extends circumferentially along the pitch axis direction of the blade 1, so that the blade 1 is completely wrapped by the self-rotation area 2.

[0049] It should be noted that in the case where the self-rotation area is not set in the three-dimensional flow field model of the wind turbine, for different pitch angles, the pitch angle of the blade needs to be adjusted multiple times, and the flow field needs to be meshed multiple times. In the self-rotation area of the three-dimensional flow field model in the embodiment of the present invention, during the blade pitching process, the grids of the self-rotation area can be directly rotated, and there is no need to mesh the entire flow field model.

[0050] It should be noted that the grid nodes in the self-rotation area 2 of the embodiment of the present invention can rotate, and the rotation speed of the above grid nodes is related to the blade pitch angle of the wind turbine. Based on the blade pitch angle of the wind turbine, the rotation speed of the above grid nodes can be controlled, so that during the simulation test of the wind turbine, during the dynamic change of the above blade pitch angle, the change of the flow field of the wind turbine can be simulated.

[0051] It should be noted that the grid nodes in the self-rotation area 2 include the grid nodes on the blade wall surface.

[0052] The revolution area 3 extends circumferentially along the axis of rotation of the wind wheel in the wind turbine, so that the self-rotation area 2 is completely wrapped by the revolution area 3.

[0053] The static region 4 is nested outside the revolution region 3 and extends circumferentially along the direction of the wind turbine's rotor rotation axis, such that the revolution region 3 is completely wrapped by the static region 4.

[0054] Under normal circumstances, the number of blades in a wind turbine is 3. Taking the number of blades 1 in a wind turbine as 3 as an example, the wind turbine flow field simulation test method provided by the present invention will be described below.

[0055] As Figures 2 to 4 shown, in the three-dimensional flow field model of blade 1 in a wind turbine, the cross-sections of the static region 4 and the revolution region 3 corresponding to blade 1 are both fan-shaped. The center of the above-mentioned fan is the rotation center of the wind turbine's rotor, and the central angle of the above-mentioned fan is 360° / the number of blades. When the number of blades in a wind turbine is 3, the central angles of the above-mentioned fans are all 120°.

[0056] When performing fluid dynamics simulation on a wind turbine, the wind speed data of the wind turbine's oncoming flow can be set according to the simulation requirements. The above-mentioned wind speed data can be time-series data that is monotonically increasing, time-series data that is monotonically decreasing, pulsating time-series data that satisfies the turbulent wind spectrum, or time-series data that satisfies a specific function.

[0057] Optionally, in the embodiments of the present invention, Fluent software can be used to perform fluid dynamics simulation on a wind turbine.

[0058] When the wind speed data of the wind turbine's oncoming flow is pulsating time-series data that satisfies the turbulent wind spectrum, the profile file function of Fluent software can be used to simulate the change of wind speed over time, that is, the wind speed at the inlet boundary of the flow field at different times can be defined according to the data in the profile file.

[0059] When the wind speed data of the wind turbine's oncoming flow is time-series data that satisfies a specific function, the wind speed at the inlet boundary of the flow field at different times can be defined through the DEFINE_PROFILE macro in the secondary development module (UDF) of Fluent software.

[0060] During the process of performing fluid dynamics simulation on a wind turbine based on the three-dimensional flow field model of the wind turbine, the wind speed at the inlet boundary of the above-mentioned three-dimensional flow field model at the current moment can be obtained through the above method.

[0061] Step 102: Based on the simulation data and the control strategy of the wind turbine, obtain the simulation control parameters of the wind turbine at the next moment.

[0062] Specifically, based on the wind speed, the aerodynamic torque of the wind turbine, the rotational speed of the wind turbine (equivalent to the rotational speed of the reference coordinate system in the revolution area), and the blade pitch angle at the current moment in the three-dimensional flow field model of the wind turbine, the computational fluid dynamics method can be used to execute the corresponding wind turbine rotational speed control algorithm according to the control strategy of the wind turbine, perform transient numerical calculations, and obtain the wind turbine rotational speed and the pitch speed of the blade at the next moment of the wind turbine.

[0063] According to the aerodynamic torque and the rotational speed of the wind turbine at the current moment, the output power of the wind turbine at the current moment can be calculated; according to the wind speed, the output power of the wind turbine, the rotational speed of the wind turbine, and the blade pitch angle in the three-dimensional flow field model at the current moment, the motor torque at the current moment can be calculated based on the control algorithm; according to the aerodynamic torque of the wind turbine, the motor torque, and the moment of inertia of the blade (inherent property of the blade) at the current moment, the rotational acceleration of the wind turbine at the current moment can be calculated; according to the time interval between the current moment and the next moment, the rotational speed of the wind turbine at the current moment, and the acceleration of the wind turbine, the rotational speed of the wind turbine at the next moment of the wind turbine can be calculated.

[0064] According to the time interval between the current moment and the next moment, the blade pitch angle at the current moment, and the pitch acceleration of the blade, the pitch speed of the blade at the next moment can be calculated.

[0065] It should be noted that the rotational speed of the wind turbine at the next moment of the wind turbine can be the same as or different from the rotational speed of the wind turbine at the current moment. The pitch speed of the wind turbine at the next moment can be the same as or different from the pitch speed of the wind turbine at the current moment.

[0066] It should be noted that the control strategy of the wind turbine can include, but is not limited to, the critical wind speed at which the wind turbine starts to pitch. The control strategy of the wind turbine can be determined in advance based on prior knowledge. In the embodiments of the present invention, the control strategy of the wind turbine is not specifically limited.

[0067] Step 103: Update the rotational speed of the reference coordinate system in the revolution area and the rotational speed of the grid nodes in the rotation area in the three-dimensional flow field model based on the simulation control parameters.

[0068] Specifically, based on the rotational speed of the wind turbine at the next moment of the wind turbine, the rotational speed of the reference coordinate system in the revolution area 3 at the next moment can be obtained, and the rotational speed of the reference coordinate system in the revolution area 3 at the next moment can be updated.

[0069] Based on the pitch speed of the blade at the next moment of the wind turbine, the rotational speed of the grid nodes in the rotation area 2 at the next moment can be obtained, and the rotational speed of the grid nodes in the rotation area 2 at the next moment can be updated.

[0070] The following uses an example to illustrate the wind turbine flow field simulation test method provided by the present invention.

[0071] Set the wind speed in the above three-dimensional flow field model at the initial moment of simulation to 5 m / s. During the simulation, the change in the wind speed in the above three-dimensional flow field model is determined by a randomly generated time series data. Set the wind turbine rotational speed at the initial moment of simulation to 3 rpm and the blade pitch angle to 0°. Among them, when the blade pitch angle is 0°, the blade is fully deployed.

[0072] Based on the three-dimensional flow field model of the wind turbine, perform transient fluid dynamics simulation of the wind turbine at the current moment, and the wind speed, the torque received by each blade, the wind turbine rotational speed, and the blade pitch angle in the above three-dimensional flow field model at the current moment can be obtained. Further obtain the sum of the torques received by each blade at the current moment, and the aerodynamic torque of the wind turbine at the current moment can be obtained.

[0073] According to the aerodynamic torque of the wind turbine and the wind turbine rotational speed at the current moment, the output power of the wind turbine at the current moment can be calculated.

[0074] According to the output power, the wind turbine rotational speed, and the blade pitch angle of the wind turbine at the current moment, based on the control algorithm, the motor torque and the pitch acceleration of the blade at the current moment can be calculated.

[0075] According to the aerodynamic torque of the wind turbine, the motor torque, and the moment of inertia of the blade at the current moment, the rotational acceleration of the wind turbine at the current moment can be calculated.

[0076] According to the time interval between the current moment and the next moment, the wind turbine rotational speed at the current moment, and the rotational acceleration of the wind turbine, calculate the wind turbine rotational speed at the next moment.

[0077] According to the time interval between the current moment and the next moment, the blade pitch angle at the current moment, and the pitch acceleration of the blade, calculate the pitch speed of the blade at the next moment.

[0078] According to the wind turbine rotational speed at the next moment, update the rotational speed of the reference coordinate system in the revolution area, and according to the pitch speed of the blade at the next moment, update the rotational speed of the grid nodes in the rotation area, so that the blade rotates to the blade pitch angle at the next moment.

[0079] In the embodiment of the present invention, by constructing and based on the three-dimensional flow field model of the wind turbine, transient fluid dynamics simulation of the wind turbine at the current moment is carried out. After obtaining the simulation data of the wind turbine and the three-dimensional flow field model at the current moment, based on the above simulation data and the control strategy of the wind turbine, the simulation control parameters of the wind turbine at the next moment are obtained. Based on the above simulation control parameters, the rotational speed of the reference coordinate system in the revolution area and the rotational speed of the grid nodes in the rotation area in the above three-dimensional flow field model are updated. It can perform unsteady computational fluid dynamics analysis on the wind turbine variable speed and blade pitch change processes during the operation of the wind turbine. The variable speed and pitch change actions of the wind turbine are considered in the simulation, and the obtained transient flow field calculation results are closer to the actual situation, and more accurate simulation test results can be obtained, which is of great significance for improving the efficiency, reliability, correctness and unity of the design and implementation of the wind turbine. When performing steady-state simulation calculations on the flow field of the wind turbine at different blade pitch angles, by controlling the rotation of the rotation area to make the blade in the target blade pitch angle, the time cost of the pre-processing of the flow field calculation can be effectively reduced.

[0080] Based on the content of the above embodiments, including: the revolution area in the three-dimensional flow field model is established based on the multiple reference frame model, the spatial coordinates of the grid nodes in the revolution area are fixed, and the multiple reference frame model is used to add a relative rotational speed to the revolution area.

[0081] Specifically, the three-dimensional flow field model can be constructed based on the Cartesian coordinate system.

[0082] The multiple coordinate system model (MRF model) is a steady-state calculation model. In the model, it is assumed that the grid cells move at a constant speed. This method is applicable to problems where the relative motions of points on the boundary of the grid region are basically the same. Most time-averaged flows can be calculated using the MRF model, especially when the interaction between the moving grid region and the stationary grid region is relatively weak, the MRF model can be used for calculation, such as the internal flow field calculation of stirrers, the internal flow field calculation of pumps and fans, etc. Another use of the MRF model is to provide an initial flow field for the sliding grid model calculation, that is, first roughly calculate the initial flow field using the MRF model, and then complete the entire calculation using the sliding grid model.

[0083] It should be noted that by adding a rotational speed component to the reference coordinate system of the revolution area 3 through the MRF model, the flow field changes caused by the rotation of the wind turbine can be equivalently simulated. By changing the rotational speed of the MRF model, the simulation of the flow field near the wind turbine during the variable speed process of the wind turbine can be realized. The rotational speed of the wind turbine is related to the aerodynamic torque of the wind turbine, the motor torque, and the moment of inertia of the blade (the inherent property of the blade).

[0084] In the embodiments of the present invention, the MRF model is used to add a rotational speed component to the reference coordinate system of the revolution area 3. In combination with the periodic symmetric boundary, it is possible to simulate the rotational flow field near the entire wind turbine rotor during the dynamic change of the wind turbine rotor speed.

[0085] Based on the content of the above embodiments, the rotation area in the three-dimensional flow field model is established based on the sliding mesh model. The grid nodes can rotate around the central axis of the rotation area. The sliding mesh model is used to change the positions of the grid nodes and the boundary of the three-dimensional flow field model to simulate the flow field change during blade pitch change.

[0086] Specifically, in the sliding mesh model, during the calculation process, the moving unit area slides along the grid interface, and the internal grid of the moving grid area remains unchanged. This feature makes the sliding mesh model have great advantages in numerical simulation when there are problems related to the rotating area.

[0087] In the embodiments of the present invention, the sliding mesh model is used to model the rotation of the grid nodes in the rotation area 2. Thus, after obtaining the rotational speed of the grid nodes in the rotation area 2 at the next moment, the movement of the grid nodes in the rotation area 2 can be controlled based on the sliding mesh model, and rotate around the pitch axis of the blade wrapped by the rotation area 2 according to the above rotational speed. Therefore, during the dynamic change of the blade pitch angle in the wind turbine, the rotational flow field near each blade in the wind turbine can be simulated.

[0088] Figure 5 It is the second flow schematic diagram of the wind turbine flow field simulation test method provided by the present invention. The implementation process of the wind turbine flow field simulation test method provided by the present invention in the Fluent software is as Figure 5 shown.

[0089] Through multiple secondary development modules (UDF) of the Fluent software, the solver in the Fluent software is secondary developed. The wind speed function at the flow field inlet can be realized through the UDF module with the DEFINE_PROFILE macro as the main body.

[0090] After the numerical iteration at the current moment is completed, the instantaneous simulation data of the wind turbine at the current moment is obtained, and the instantaneous simulation data of the wind turbine at the current moment is transmitted to the UDF with the DEFINE_ADJUST macro as the main body. The aerodynamic load of the blade is extracted by using the Compute_Force_And_Moment macro, and the wind turbine torque and output power are calculated. By judging whether the wind turbine reaches the rated speed and whether it exceeds the rated power, according to the control strategy of the wind turbine, the corresponding wind turbine speed and blade pitch angle control algorithms can be executed to determine the wind turbine speed and blade pitch angle at the next moment.

[0091] Through the DEFINE_ZONE_MOTION macro, the rotational speed of the sliding mesh model in the self-rotation zone 2 and the rotational speed of the MRF multiple reference frame model in the revolution zone 3 are adjusted.

[0092] As the simulation time progresses, when the next moment becomes the current moment, the above calculation process is repeated.

[0093] Based on the content of the above embodiments, in the three-dimensional flow field model, boundary layer grids are set on the surface of the blade, and the thickness of the first layer of boundary layer grids meets the preset conditions; the revolution zone and the stationary zone adopt a contact method of sharing grid nodes; data transfer is realized between the self-rotation zone and the revolution zone through an interface with non-shared nodes; the grid size of the revolution zone is larger than that of the self-rotation zone.

[0094] Specifically, in the spatial discretization of the three-dimensional flow field model of the wind turbine, boundary layer grids need to be set on the surface of blade 1, and the thickness of the first layer of boundary layer grids meets the preset conditions.

[0095] Optionally, considering the influence of the transition effect, the preset conditions may include: determining the thickness of the first layer of boundary layer grids with the goal of Y+ being less than 1. Among them, Y+ is a basic index in the art for measuring whether the grid thickness is appropriate.

[0096] The revolution zone 3 and the stationary zone 4 can adopt a contact method of sharing grid nodes, and data transfer between the self-rotation zone 2 and the revolution zone 3 is realized through an interface with non-shared nodes.

[0097] When meshing the three-dimensional flow field model of the wind turbine, the distribution of the Reynolds number (Re number) along the blade span can be calculated first according to parameters such as the wind turbine rotational speed, chord length, wind speed, air density, and dynamic viscosity; secondly, with the goal of Y+ being less than 1, the thickness of the first layer of boundary layer grids is determined, and the thickness of the first layer of wall grids is calculated; thirdly, 15 to 20 layers of boundary layer grids are generated with an increment factor of 1.05.

[0098] The self-rotation zone 2 adopts a structured grid generation method (for example: the "O-block" in ICEM software can be used for meshing). The grid nodes in the self-rotation zone 2 do not correspond to those in the revolution zone 3, but it is necessary to ensure that the grid area ratio on both sides of the interface between the self-rotation zone 2 and the revolution zone 3 is less than 4.

[0099] The stationary zone 4 can share grid nodes with the revolution zone 3. The stationary zone 4 and the revolution zone 3 adopt unstructured grids.

[0100] The grid sizes of the stationary zone 4 and the revolution zone 3 are larger than those of the self-rotation zone 2.

[0101] In the three-dimensional flow field model according to the embodiments of the present invention, boundary layer grids are provided on the surface of the blade, and the thickness of the first layer of boundary layer grids meets the preset conditions. The rotation area and the stationary area adopt a contact method with shared grid nodes, and data transfer is achieved through an interface with non-shared nodes between the rotation area and the revolution area. The grid size of the revolution area is larger than that of the rotation area, and it can more accurately simulate the flow field change of the wind turbine during the dynamic change of the wind turbine speed and / or the blade pitch angle.

[0102] Based on the content of the above embodiments, the rotation area is a cylinder; the central axis of the rotation area coincides with the pitch axis of the blade; the distance between one end of the rotation area and the rotation center of the wind wheel is a preset value; the radius of the rotation area is determined based on the maximum distance between the blade surface and the pitch axis of the blade; the length of the rotation area is determined based on the length of the blade.

[0103] Specifically, as Figures 2 to 4 shown, the rotation area 2 is a cylinder, and the central axis of the rotation area 2 coincides with the pitch axis of the blade 1.

[0104] The rotation area 2 has two ends, and the distance between one of the ends and the rotation center of the wind wheel of the wind turbine is a preset value.

[0105] Optionally, the value range of the above preset value can be between 4 and 8 cm, for example: the value of the preset value can be 4 cm, 6 cm or 8 cm.

[0106] The radius r1 of the rotation area 2, that is, the bottom radius of the above cylinder, can be determined based on the maximum distance between the surface of the blade 1 and the pitch axis of the blade 1.

[0107] Optionally, the radius r1 of the rotation area 2 can be 1 to 2 times the above maximum distance. For example, the radius r1 of the rotation area 2 can be 1 time, 1.5 times or 2 times the above maximum distance.

[0108] Preferably, the radius r1 of the rotation area 2 can be 1.5 times the above maximum distance.

[0109] The length h1 of the rotation area 2, that is, the height of the above cylinder, can be determined based on the length h of the blade 1.

[0110] Optionally, the length h1 of the rotation area 2 can be 0.5 to 1 m longer than the length h of the blade 1. For example: the length h1 of the rotation area 2 can be 0.5 m, 0.75 m or 1 m longer than the length h of the blade 1.

[0111] In the embodiment of the present invention, the self-rotation area is a cylinder, the central axis of the self-rotation area coincides with the pitch axis of the blade, one end of the self-rotation area is located at the rotation center of the wind turbine, the radius of the self-rotation area is determined based on the maximum distance between the blade surface and the pitch axis of the blade, and the length of the self-rotation area is determined based on the length of the blade, which can more accurately simulate the rotational flow field near each blade of the wind turbine during the dynamic change of the blade pitch angle of the wind turbine.

[0112] Based on the content of the above embodiments, the revolution area is a cylinder, the central axis of the revolution area coincides with the rotation axis of the wind turbine, and the rotation center of the wind turbine coincides with the midpoint of the central axis of the revolution area; the radius of the revolution area is determined based on the length of the self-rotation area; the length of the revolution area is determined based on the radius of the self-rotation area.

[0113] Specifically, as Figures 2 to 4 shown, the revolution area 3 is a cylinder, and the central axis of the revolution area 3 coincides with the rotation axis of the wind turbine in the wind turbine, and the rotation center of the wind turbine coincides with the midpoint of the central axis of the revolution area 3.

[0114] The radius r2 of the revolution area 3, that is, the bottom radius of the above cylinder, can be determined based on the length h1 of the self-rotation area 2.

[0115] Optionally, the radius r2 of the revolution area 3 can be 1.05 to 1.1 times the length h1 of the self-rotation area 2. For example, the radius r2 of the revolution area 3 can be 1.05 times, 1.075 times or 1.1 times the length h1 of the self-rotation area 2.

[0116] The length h2 of the revolution area 3, that is, the height of the above cylinder, can be determined based on the radius r1 of the self-rotation area 2.

[0117] Optionally, the length h2 of the revolution area 3 can be 3 to 4 times the radius r1 of the self-rotation area 2. For example: the length h2 of the revolution area 3 can be 3 times, 2.5 times or 4 times the radius r1 of the self-rotation area 2.

[0118] In the embodiment of the present invention, the revolution area is a cylinder, the central axis of the revolution area coincides with the rotation axis of the wind turbine, the rotation center of the wind turbine coincides with the midpoint of the central axis of the revolution area, the radius of the revolution area is determined based on the length of the self-rotation area, and the length of the revolution area is determined based on the radius of the self-rotation area, which can more accurately simulate the rotational flow field of the entire wind turbine in the wind turbine during the dynamic change of the wind turbine speed.

[0119] Based on the content of the above embodiments, the stationary area is an annular cylinder; the central axis of the stationary area coincides with the rotation axis of the wind turbine; the outer diameter of the stationary area is determined based on the radius of the revolution area; the distance between the inlet boundary of the stationary area and the wind turbine group in the wind turbine, and the distance between the outlet boundary of the stationary area and the wind turbine group are determined based on the radius of the self-rotation area.

[0120] Specifically, as Figures 2 to 4 shown, the static region 4 is nested outside the revolution region 3. The static region 4 is a cylindrical ring, and the central axis of the static region 4 coincides with the rotation axis of the wind turbine.

[0121] The inner diameter r3 of the static region 4 is the radius r2 of the revolution region 3.

[0122] The outer diameter r4 of the static region 4 can be determined based on the radius r2 of the revolution region 3.

[0123] Optionally, the outer diameter r4 of the static region 4 can be 1.5 to 2 times the radius r2 of the revolution region 3. For example, the outer diameter r4 of the static region 4 can be 1.5 times, 1.75 times, or 2 times the radius r2 of the revolution region 3.

[0124] The inlet boundary of the static region 4 is one of the two circular surfaces of the above-mentioned cylinder that the oncoming flow first passes through. The outlet boundary of the static region 4 is one of the two circular surfaces of the above-mentioned cylinder that the oncoming flow passes through later.

[0125] The distance h3 between the inlet boundary of the static region 4 and the wind turbine group can be determined based on the radius r1 of the rotation region 2.

[0126] Optionally, the distance h3 between the inlet boundary of the static region 4 and the wind turbine group needs to be greater than 8 to 12 times the radius r1 of the rotation region 2. For example, the distance h3 between the inlet boundary of the static region 4 and the wind turbine group needs to be greater than 8 times, 10 times, or 12 times the radius r1 of the rotation region 2.

[0127] Preferably, the distance h3 between the inlet boundary of the static region 4 and the wind turbine group needs to be greater than 10 times the radius r1 of the rotation region 2.

[0128] The distance h4 between the outlet boundary of the static region 4 and the wind turbine group can also be determined based on the radius r1 of the rotation region 2.

[0129] Optionally, the distance h4 between the outlet boundary of the static region 4 and the wind turbine group needs to be greater than 18 to 22 times the radius r1 of the rotation region 2. For example, the distance h4 between the outlet boundary of the static region 4 and the wind turbine group needs to be greater than 18 times, 20 times, or 22 times the radius r1 of the rotation region 2.

[0130] Preferably, the distance h4 between the outlet boundary of the static region 4 and the wind turbine group needs to be greater than 20 times the radius r1 of the rotation region 2.

[0131] In the embodiment of the present invention, the static area is a cylindrical ring, the central axis of the static area coincides with the rotation axis of the wind turbine, the outer diameter of the static area is determined based on the radius of the revolution area, and the distance between the inlet boundary of the static area and the wind turbine group, as well as the distance between the outlet boundary of the static area and the wind turbine group, are determined based on the radius of the rotation area, which can more accurately simulate the flow field change of the wind turbine during the dynamic change of the wind turbine speed and the blade pitch angle in the wind turbine.

[0132] Figure 6 It is a schematic structural diagram of the wind turbine simulation test device provided by the present invention. The following combines Figure 6 to describe the wind turbine simulation test device provided by the present invention. The wind turbine simulation test device described below can be correspondingly referred to the wind turbine flow field simulation test method provided by the present invention described above. As Figure 6 shown, the device includes: a data acquisition module 601, a data calculation module 602, and a simulation control module 603.

[0133] The data acquisition module 601 is used to construct and perform transient fluid dynamics simulation on the wind turbine based on the three-dimensional flow field model of the wind turbine, and obtain the simulation data of the wind turbine and the three-dimensional flow field model at the current moment.

[0134] The data calculation module 602 is used to obtain the simulation control parameters of the wind turbine at the next moment based on the simulation data and the control strategy of the wind turbine.

[0135] The simulation control module 603 is used to update the rotational speed of the reference coordinate system in the revolution area and the rotational speed of the grid nodes in the rotation area in the three-dimensional flow field model based on the simulation control parameters;

[0136] Among them, the three-dimensional flow field model includes: a rotation area, a revolution area, and a static area; the rotation area extends circumferentially along the pitch axis direction of the blades in the wind turbine, and the radius of the rotation area is greater than the maximum distance between the blade surface and the pitch axis of the blade; the revolution area and the static area extend circumferentially along the rotation axis direction of the wind turbine in the wind turbine, the radius of the revolution area is greater than the length of the blade, and the static area is nested outside the revolution area.

[0137] Specifically, the data acquisition module 601, the data calculation module 602, and the simulation control module 603 are electrically connected.

[0138] During the process of performing fluid dynamics simulation on the wind turbine based on the three-dimensional flow field model of the wind turbine, the data acquisition module 601 can be used to construct the above three-dimensional flow field model and obtain the wind speed distribution, the aerodynamic torque of the wind turbine, the wind turbine speed, the blade pitch angle, and the motor torque in the above three-dimensional flow field model.

[0139] The data calculation module 602 can be used to perform transient numerical calculations based on the simulation data at the current moment and the control strategy of the wind turbine, and obtain the motor torque of the wind turbine at the next moment based on the computational fluid dynamics method.

[0140] The simulation control module 603 can be used to obtain the rotational speed of the reference coordinate system in the revolution area 3 at the next moment based on the motor torque of the wind turbine at the next moment, and update the rotational speed of the reference coordinate system in the revolution area 3 at the next moment. The simulation control module 603 can be used to obtain the rotational speed of the grid nodes in the rotation area 2 based on the pitch speed of the wind turbine at the next moment, and can update the rotational speed of the grid nodes in the rotation area 2 at the next moment.

[0141] Optionally, the wind turbine simulation test device may further include a model construction module.

[0142] The model construction module can be used to construct a three-dimensional flow field model.

[0143] In the embodiment of the present invention, by constructing and based on the three-dimensional flow field model of the wind turbine, performing transient fluid dynamics simulation on the wind turbine at the current moment, obtaining the simulation data of the wind turbine and the three-dimensional flow field model at the current moment, then based on the above simulation data and the control strategy of the wind turbine, obtaining the simulation control parameters of the wind turbine at the next moment, and based on the above simulation control parameters, updating the rotational speed of the reference coordinate system in the revolution area and the rotational speed of the grid nodes in the rotation area in the above three-dimensional flow field model, it is possible to perform unsteady computational fluid dynamics analysis on the wind wheel speed change and blade pitch change processes during the operation of the wind turbine. The variable speed and variable pitch actions of the wind turbine are considered in the simulation, and the obtained transient flow field calculation results are closer to the actual situation, and more accurate simulation test results can be obtained, which is of great significance for improving the efficiency, reliability, correctness and unity of the design and implementation of the wind turbine. When performing steady-state simulation calculations on the flow field of the wind turbine at different blade pitch angles, by controlling the rotation of the rotation area to make the blade in the target blade pitch angle, the time cost of the preprocessing of the flow field calculation can be effectively reduced.

[0144] Figure 7 An example of the physical structure diagram of an electronic device is shown in Figure 7As shown in the figure, the electronic device may include: a processor 710, a communications interface 720, a memory 730, and a communication bus 740. Among them, the processor 710, the communications interface 720, and the memory 730 complete communication with each other through the communication bus 740. The processor 710 may call the logical instructions in the memory 730 to execute the wind turbine flow field simulation test method, which includes: constructing and based on a three-dimensional flow field model of the wind turbine, performing transient fluid dynamics simulation on the wind turbine at the current moment to obtain simulation data of the wind turbine and the three-dimensional flow field model at the current moment; obtaining the simulation control parameters of the wind turbine at the next moment based on the simulation data and the control strategy of the wind turbine; updating the rotational speed of the reference coordinate system in the revolution area and the rotational speed of the grid nodes in the rotation area of the three-dimensional flow field model based on the simulation control parameters; where the three-dimensional flow field model includes: a rotation area, a revolution area, and a stationary area; the rotation area extends circumferentially along the pitch axis direction of the blade in the wind turbine, and the radius of the rotation area is greater than the maximum distance between the blade surface and the pitch axis of the blade; the revolution area and the stationary area extend circumferentially along the axis direction of the wind turbine rotor rotation, the radius of the revolution area is greater than the length of the blade, and the stationary area is nested outside the revolution area.

[0145] In addition, when the logical instructions in the above-mentioned memory 730 are implemented in the form of software function units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or a part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present invention. The foregoing storage medium includes: various media such as USB flash drives, mobile hard disks, read-only memories (ROM, Read-Only Memory), random access memories (RAM, Random Access Memory), magnetic disks, or optical discs that can store program codes.

[0146] On the other hand, the present invention also provides a computer program product. The computer program product includes a computer program stored on a non-transitory computer-readable storage medium. The computer program includes program instructions. When the program instructions are executed by a computer, the computer can execute the wind turbine flow field simulation test method provided by each of the above methods. The method includes: constructing and based on a three-dimensional flow field model of a wind turbine, performing a transient fluid dynamics simulation of the wind turbine at the current moment, and obtaining simulation data of the wind turbine and the three-dimensional flow field model at the current moment; based on the simulation data and the control strategy of the wind turbine, obtaining simulation control parameters of the wind turbine at the next moment; based on the simulation control parameters, updating the rotational speed of the reference coordinate system in the revolution area and the rotational speed of the grid nodes in the rotation area in the three-dimensional flow field model; wherein, the three-dimensional flow field model includes: a rotation area, a revolution area, and a stationary area; the rotation area extends circumferentially along the pitch axis direction of the blades in the wind turbine, and the radius of the rotation area is greater than the maximum distance between the blade surface and the pitch axis of the blade; the revolution area and the stationary area extend circumferentially along the direction of the wind turbine rotor rotation axis, the radius of the revolution area is greater than the length of the blade, and the stationary area is nested outside the revolution area.

[0147] In another aspect, the present invention also provides a non-transitory computer-readable storage medium, on which a computer program is stored. When the computer program is executed by a processor, it is configured to execute the wind turbine flow field simulation test method provided by each of the above. The method includes: constructing and based on a three-dimensional flow field model of a wind turbine, performing a transient fluid dynamics simulation of the wind turbine at the current moment, and obtaining simulation data of the wind turbine and the three-dimensional flow field model at the current moment; based on the simulation data and the control strategy of the wind turbine, obtaining simulation control parameters of the wind turbine at the next moment; based on the simulation control parameters, updating the rotational speed of the reference coordinate system in the revolution area and the rotational speed of the grid nodes in the rotation area in the three-dimensional flow field model; wherein, the three-dimensional flow field model includes: a rotation area, a revolution area, and a stationary area; the rotation area extends circumferentially along the pitch axis direction of the blades in the wind turbine, and the radius of the rotation area is greater than the maximum distance between the blade surface and the pitch axis of the blade; the revolution area and the stationary area extend circumferentially along the direction of the wind turbine rotor rotation axis, the radius of the revolution area is greater than the length of the blade, and the stationary area is nested outside the revolution area.

[0148] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. Those of ordinary skill in the art can understand and implement it without creative labor.

[0149] Through the description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware. Based on such an understanding, the essence of the above technical solution, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product can be stored in a computer-readable storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions for causing a computer device (which can be a personal computer, a server, or a network device, etc.) to execute the methods described in each embodiment or some parts of the embodiments.

[0150] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions described in the foregoing embodiments, or perform equivalent replacements for some of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of each embodiment of the present invention.

Claims

1. A method for simulating and testing the flow field of a wind turbine, characterized in that, Including: Construct and, based on a three-dimensional flow field model of a wind turbine, perform transient fluid dynamics simulation on the wind turbine at the current moment to obtain simulation data of the wind turbine and the three-dimensional flow field model at the current moment; Based on the simulation data and the control strategy of the wind turbine, obtain the simulation control parameters of the wind turbine at the next moment; Based on the simulation control parameters, update the rotational speed of the reference coordinate system in the revolution area and the rotational speed of the grid nodes in the rotation area in the three-dimensional flow field model; Wherein, the three-dimensional flow field model includes a rotation area, a revolution area, and a stationary area; the rotation area extends circumferentially along the pitch axis direction of the blades in the wind turbine, and the radius of the rotation area is greater than the maximum distance between the blade surface and the pitch axis of the blade; the revolution area and the stationary area extend circumferentially along the direction of the wind turbine's rotor rotation axis, the radius of the revolution area is greater than the length of the blade, and the stationary area is nested outside the revolution area; the simulation data includes the wind speed distribution, the aerodynamic torque of the rotor, the rotor speed, the blade pitch angle, and the motor torque in the three-dimensional flow field model; the simulation control parameters include the rotor speed, the pitch speed of the blade, the rotational speed of the reference coordinate system in the revolution area, and the rotational speed of the grid nodes in the rotation area; The revolution area in the three-dimensional flow field model is established based on the multiple reference frame model, the spatial coordinates of the grid nodes in the revolution area are fixed, and the multiple reference frame model is used to add a relative rotational speed to the revolution area; The rotation area in the three-dimensional flow field model is established based on the sliding mesh model, the grid nodes can rotate around the central axis of the rotation area, and the sliding mesh model is used to change the positions of the grid nodes and the boundary of the three-dimensional flow field model to simulate the flow field change during blade pitching; In the three-dimensional flow field model, boundary layer grids are arranged on the surface of the blade, and the thickness of the first layer of the boundary layer grid meets a preset condition; the revolution area and the stationary area adopt a contact method of sharing grid nodes; data transfer between the rotation area and the revolution area is achieved through an interface of non-shared nodes; the grid size of the revolution area is greater than the grid size of the rotation area.

2. The method for simulating and testing the flow field of a wind turbine according to claim 1, wherein The rotation area is a cylinder; the central axis of the rotation area coincides with the pitch axis of the blade; the distance from one end of the rotation area to the rotation center of the rotor is a preset value; the radius of the rotation area is determined based on the maximum distance between the blade surface and the pitch axis of the blade; the length of the rotation area is determined based on the length of the blade.

3. The method for simulating and testing the flow field of a wind turbine according to claim 2, wherein The revolution area is a cylinder, the central axis of the revolution area coincides with the rotation axis of the rotor, and the rotation center of the rotor coincides with the midpoint of the central axis of the revolution area; the radius of the revolution area is determined based on the length of the rotation area; the length of the revolution area is determined based on the radius of the rotation area.

4. The method for simulating and testing the flow field of a wind turbine according to claim 3, characterized in that, The static region is an annular cylinder; the central axis of the static region coincides with the rotation axis of the wind turbine; the outer diameter of the static region is determined based on the radius of the revolution region; the distance between the inlet boundary of the static region and the wind turbine rotor group, and the distance between the outlet boundary of the static region and the rotor group are determined based on the radius of the rotation region.

5. The method for simulating and testing the flow field of a wind turbine according to any one of claims 1 to 4, characterized in that, The three-dimensional flow field model is composed of the three-dimensional flow field models of each blade; the boundary of the three-dimensional flow field model of any blade is a periodic symmetry boundary.

6. A wind turbine simulation test device, characterized in that, Including: A data acquisition module, configured to construct and perform transient fluid dynamics simulation on the wind turbine at the current moment based on the three-dimensional flow field model of the wind turbine, and obtain the simulation data of the wind turbine and the three-dimensional flow field model at the current moment; A data calculation module, configured to obtain the simulation control parameters of the wind turbine at the next moment based on the simulation data and the control strategy of the wind turbine; A simulation control module, configured to update the rotational speed of the reference coordinate system in the revolution region and the rotational speed of the grid nodes in the rotation region in the three-dimensional flow field model based on the simulation control parameters; Wherein, the three-dimensional flow field model includes: a rotation region, a revolution region, and a static region; the rotation region extends circumferentially along the pitch axis direction of the blades in the wind turbine, and the radius of the rotation region is greater than the maximum distance between the blade surface and the pitch axis of the blade; the revolution region and the static region extend circumferentially along the rotation axis direction of the wind turbine rotor, the radius of the revolution region is greater than the length of the blade, and the static region is nested outside the revolution region; The simulation data includes: the wind speed distribution in the three-dimensional flow field model, the aerodynamic torque of the wind turbine rotor, the wind turbine rotor speed, the blade pitch angle, and the motor torque; the simulation control parameters include: the wind turbine rotor speed, the pitch speed of the blade, the rotational speed of the reference coordinate system in the revolution region, and the rotational speed of the grid nodes in the rotation region; The revolution region in the three-dimensional flow field model is established based on the multiple reference frame model, the spatial coordinates of the grid nodes in the revolution region are fixed, and the multiple reference frame model is used to add relative rotational speed to the revolution region; The rotation region in the three-dimensional flow field model is established based on the sliding mesh model, the grid nodes can rotate around the central axis of the rotation region, and the sliding mesh model is used to change the positions of the grid nodes and the boundary of the three-dimensional flow field model to simulate the flow field change during blade pitching; In the three-dimensional flow field model, boundary layer grids are arranged on the surface of the blade, and the thickness of the first layer of the boundary layer grids meets a preset condition; the revolution region and the static region adopt a contact method of sharing grid nodes; data transfer is realized between the rotation region and the revolution region through an interface with non-shared nodes; the grid size of the revolution region is larger than the grid size of the rotation region.

7. An electronic device, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized in that, When the processor executes the program, it realizes the steps of the wind turbine flow field simulation test method according to any one of claims 1 to 5.

8. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it realizes the steps of the wind turbine flow field simulation test method according to any one of claims 1 to 5.

Images