CN-122021397-A - High-efficiency high-fidelity coupling simulation method and simulation system for floating wind turbine

CN122021397ACN 122021397 ACN122021397 ACN 122021397ACN-122021397-A

Abstract

The invention discloses a high-efficiency high-fidelity coupling simulation method of a floating wind turbine, which comprises the steps of determining the current position of a hub reference point under global coordinates based on rigid body pose, determining a unit where the hub position is located and an owned process thereof in a parallel partition grid, reading the flow velocity in the unit in the owned process, constructing the relative wind velocity at the hub as the local flow velocity, subtracting the linear velocity obtained by cross multiplying the translational velocity of the rigid body and the angular velocity by the sagittal diameter, converting the reference axis into the current axis direction of the wind turbine according to quaternion pose to obtain the axial component of the relative wind velocity, performing quick binary search and linear interpolation on the axial wind velocity according to a preset wind speed-thrust discrete curve to obtain a thrust scalar and generating a thrust vector along the axis direction, cross multiplying the sagittal diameter of the hub relative rotation center by the thrust vector to obtain moment, combining the thrust and the flow volume component and gravity to serve as the external force and the external force moment of a 6DoF rigid body dynamics equation, and propelling the pose. The method avoids explicit modeling impeller geometry and body pasting grid division thereof, is different from an actuating line model, has higher robustness while reducing the calculation cost, ensures the gas-wave-mooring coupling precision, and is suitable for parallel calculation environments of engineering dimensions.

Inventors

LUO JIE
YANG XUELIANG
YU DENGMIN
LI CHUANGYE
Hong Chengxiao
JIANG TIAN

Assignees

中国电建集团贵州工程有限公司

Dates

Publication Date: 20260512
Application Date: 20251225

Claims (10)

1. A high-efficiency high-fidelity coupling simulation method for a floating wind turbine comprises the following steps: a1, determining the current position of a preset hub reference point under global coordinates based on the pose of a rigid body; A2, determining grid cells and owned processes thereof where the current position is located in the parallel partition grid to form a unique corresponding relationship between the owners and the cells; A3, reading the flow velocity vector in the unit in the possession process, constructing rigid body motion speed at the hub, and calculating the relative wind speed; A4, converting a reference axis into the current axis direction of the fan based on the quaternion gesture, and projecting the relative wind speed into an axial component; a5, carrying out rapid interpolation on the axial wind speed according to a wind speed-thrust discrete curve to obtain a thrust scalar and generating a thrust vector along the current axis direction of the fan; a6, obtaining moment according to the radial vector of the hub relative to the rotation center and the cross multiplication of the thrust vector; A7, combining the thrust, the moment, the flow volume component and the gravity to serve as external force and external moment of a six-degree-of-freedom rigid body dynamics equation to push the rigid body pose; a8, carrying out global consistent synchronization on the thrust and the moment in a parallel environment.
2. The method for simulating high-efficiency and high-fidelity coupling of a floating wind turbine according to claim 1, wherein the relative wind speed is calculated as follows: , Wherein, the For the flow velocity vector in the unit, v is the rigid body translation speed, ω is the rigid body angular speed, r is the vector diameter of the hub position relative to the rotation center, For the current position of the hub reference point in global coordinates, For the location of the center of rotation in global coordinates, Is the relative wind speed.
3. The high-efficiency high-fidelity coupling simulation method of the floating wind turbine of claim 1, wherein the current axis direction of the wind turbine is obtained by performing quaternion attitude transformation on a reference axis of a machine body coordinate system, the reference axis is the x-axis direction of the machine body, and the axial component is calculated according to the following formula: , Wherein, the Is an axis unit vector, Q is a quaternion of the current posture of the rigid body and is used for transforming a reference axis from an engine system to a global system, For the reference axis of the machine system, For an axial component scalar in the current axis direction relative to wind speed, Is the inflow wind speed for interpolation.
4. The method for simulating high-efficiency and high-fidelity coupling of a floating wind turbine according to claim 1, wherein the fast interpolation in step A5 comprises: A51, firstly, judging the threshold of the axial wind speed based on the cut-in wind speed and the cut-out wind speed, and if the axial wind speed is lower than the cut-in wind speed or higher than the cut-out wind speed, the thrust scalar is 0; A52, if the axial wind speed is between the cut-in wind speed and the cut-out wind speed, locating adjacent discrete points on the wind speed-thrust discrete curve through binary search, and generating a thrust scalar through linear interpolation.
5. The method for simulating high-efficiency and high-fidelity coupling of a floating wind turbine according to claim 1, wherein the axial wind speed is threshold determined based on the cut-in wind speed and the cut-out wind speed before interpolation to suppress thrust in a non-working region.
6. The method for simulating high-efficiency and high-fidelity coupling of a floating wind turbine according to claim 1, wherein the determining of the correspondence between owners and units in step A2 comprises: a21, each parallel process executes local unit searching on the current position of the hub reference point; a22, summarizing the search results of all the processes, and selecting an owner process through reduction; a23, broadcasting owner process information and corresponding grid cell information to all parallel processes by taking the owner process as a root.
7. The method for simulating high-efficiency and high-fidelity coupling of a floating wind turbine according to claim 1, wherein the step A7 comprises the steps of: c1, fluid integral force and integral moment; C2, equivalent pneumatic load formed by thrust vector and moment thereof; C3, gravity and moment caused by gravity. The thrust vector is along the current axis direction of the fan, and the moment is the cross multiplication of the vector diameter and the thrust vector.
8. A floating wind turbine section coupling simulation system for implementing the method of any of claims 1 to 8, comprising: The hub position determining module is used for calculating the global position of the hub based on the pose of the rigid body; The parallel point positioning module is used for determining owner processes and corresponding grid cells in the parallel partition grid to form a unique corresponding relationship between the owners and the cells; the field value reading and relative wind speed construction module is used for reading flow velocity vectors in the corresponding grid units in the owner process, constructing rigid body motion speed at the hub and calculating relative wind speed; The axial determining and decomposing module is used for obtaining the current axial direction of the fan through quaternion attitude transformation and projecting the relative wind speed into an axial component; The thrust interpolation module is used for carrying out rapid interpolation on the axial wind speed corresponding to the axial component based on a wind speed-thrust discrete curve to obtain a thrust scalar and generate a thrust vector along the current axial direction of the fan; The external load generation and synchronization module is used for obtaining moment according to the vector diameter of the hub relative rotation center and the cross multiplication of the thrust vector, and carrying out global consistent synchronization on the thrust and the moment in a parallel environment; And the six-degree-of-freedom dynamics module is used for combining the thrust force, the moment, the flow volume component, the fluid integral moment and the moment caused by gravity as external force and external moment of a six-degree-of-freedom rigid body dynamics equation to push the rigid body pose.
9. The floating wind turbine section coupling simulation system of claim 8, wherein the parallel point location module comprises a local lookup unit, a global reduction unit, and a broadcast unit; the local searching unit is used for executing local unit searching by each parallel process; The global reduction unit is used for summarizing the search results of all the processes and selecting an owner process through a reduction algorithm; the broadcasting unit is used for broadcasting the owner and unit information to all parallel processes by taking the owner process as a root.
10. The floating wind turbine section coupling simulation system of claim 8, wherein the thrust interpolation module comprises a threshold determination unit, a binary search unit, and a linear interpolation unit; the threshold judging unit is used for judging the threshold of the axial wind speed based on the cut-in wind speed and the cut-out wind speed; The binary search unit is used for positioning adjacent discrete points on a wind speed-thrust discrete curve for axial wind speed between cut-in wind speed and cut-out wind speed; The linear interpolation unit is used for generating a thrust scalar through linear interpolation based on the adjacent discrete points and generating a thrust vector along the axis direction.

Description

High-efficiency high-fidelity coupling simulation method and simulation system for floating wind turbine Technical Field The invention relates to a high-efficiency high-fidelity coupling simulation method and a simulation system for a floating wind turbine. Background The floating wind turbine generates complex coupling response under the action of wind, wave and current, and relates to interaction of multiple physical fields such as pneumatic, hydrodynamic, rigid dynamics, mooring dynamics and the like. With the development of offshore wind power to deep sea, the design, analysis and optimization of the floating wind turbine bring higher requirements to the high-fidelity numerical simulation technology. In the prior art, full-scale simulation of blade and tower geometry by adopting direct CFD analysis is a main method for obtaining high-precision results. According to the method, through the fine grid discretized blades, the tower and the surrounding flow field, the aerodynamic characteristics, wake effect and coupling effect with platform motion of the impeller can be accurately captured. The method has the technical problems that (1) the calculation cost is huge, the blade length can reach 60-120 m, the tower height is 80-150 m for a typical 5-15MW floating wind turbine, the full-scale CFD simulation usually needs tens of millions to hundreds of millions of grid units, the single simulation time can reach days to weeks, the calculation resource consumption is huge, (2) the parallel expansibility is limited, the parallel efficiency is obviously reduced after 96-256 cores due to the fact that the blade geometry is complex and boundary layer grids are dense, and the large-scale parallel calculation resource is difficult to fully utilize, (3) the grid generation is complex, the grid generation is difficult due to the complex geometry of the blade and the tower, particularly the quality control of the boundary layer grids, a large amount of manual intervention is needed, and (4) the time step is limited, the time step is usually limited to be about 10-3 seconds, so that the calculation cost is further increased. To reduce the computational cost, the prior art uses simplified methods such as an actuator disk model (ActuatorDiskModel, ADM) or an actuator wire model (ActuatorLineModel, ALM) to simulate the aerodynamic effects of the impeller by introducing distributed forces in the flow field. The method also has the technical problems that (1) grid sensitivity is that ADM/ALM needs to apply distribution force in a flow field, the grid sensitivity is high to grid resolution, grid quality and force distribution mode, grid change can cause significant difference of results, (2) numerical stability is that source items are introduced into the flow field, particularly, when large time step or grid quality is poor, divergence or oscillation is easy to occur, boundary condition processing is complex, coupling processing of an actuating disc/actuating line and the boundary condition of the flow field is complex, calculation convergence is possibly affected, and (4) force distribution parameterization is difficult, a force distribution function in the ADM/ALM needs to be parameterized according to specific fan model, and universality is lacking; in the integrated design process of the floating wind turbine, a key technical window period exists, namely a conceptual design stage and a preliminary design stage. At this stage, the designer is faced with the following typical challenges: Design iteration requirements are (1) the influence of different impeller specifications (power grade, impeller diameter, hub height) on the performance of the floating body platform, 2) the compatibility of the type of the floating body platform (semi-submersible type, tension leg, barge) and the impeller specifications, and (3) the geometric parameters of the platform (pontoon size, column spacing and draft) are required to be optimized rapidly to adapt to different impellers. Parameter uncertainty (1) impeller detailed aerodynamic parameters (blade airfoil, twist angle distribution, chord length distribution) are not completely determined, (2) impeller control strategies (pitch control, yaw control, power control) are still in optimization, and (3) impeller-platform coupling effects (aerodynamic-hydrodynamic-structural coupling) need to be rapidly evaluated. The time pressure is (1) the project is tense in early time, multiple schemes are required to be compared in limited time, (2) the design is frequently changed, quick response parameter adjustment is required, and (3) the cost control requirement is required, so that the economy of different configurations is required to be quickly evaluated. The prior art has limitations in the mid-early design stage. For full-scale CFD simulation, on the one hand, geometric modeling is complex, and each time impeller parameter adjustment needs to re-model blade geometry. On the other hand, grid gene