the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Going Beyond BEM with BEM: an Insight into Dynamic Inflow Effects on Floating Wind Turbines
Abstract. Blade Element Momentum (BEM) theory is the backbone of many industry-standard wind turbine aerodynamic models. To be applied to a broader set of engineering problems, BEM models have been extended since their inception and now include several empirical corrections. These models have benefitted from decades of development and refinement and have been extensively used and validated, proving their adequacy in predicting aerodynamic forces of horizontal axis wind turbine rotors in most scenarios. However, the analysis of Floating Offshore Wind Turbines (FOWTs) introduces new sets of challenges, especially if new-generation large and flexible machines are considered. In fact, due to the combined action of wind and waves and their interaction with the turbine structure and control system, these machines are subject to unsteady motion, and thus unsteady inflow on the wind turbine’s blades, which could put BEM models to the test. Consensus is not present yet on the accuracy limits of BEM in these conditions. This study contributes to the ongoing research on the topic by systematically comparing four different aerodynamic models, ranging from BEM to Computational Fluid Dynamics (CFD), in an attempt to shed light on the unsteady aerodynamic phenomena that are at stake in FOWTs and whether BEM is able to model them appropriately. Simulations are performed on the UNAFLOW 1:75 scale rotor during imposed harmonic surge and pitch motion. Experimental results are available for these conditions and are used for baseline validation. The rotor is analysed both in rated operating conditions and in low wind speeds, where unsteady aerodynamic effects are expected to be more pronounced. Results show how BEM, despite its simplicity, if augmented with a dynamic inflow model, is able to adequately model the aerodynamics of FOWTs in most conditions.
- Preprint
(4300 KB) - Metadata XML
- BibTeX
- EndNote
Status: closed
-
CC1: 'Comment on wes-2023-109', Christian Schulz, 07 Sep 2023
Dear Francesco and colleagues,
thanks for sharing this extensive work! It is great to see that the presence of the returning wake effect at the UNAFLOW rotor can also be observed in OLAF and even in ALM. I feel we finally make some progress in the field of unsteady aerodynamics of FOWT :) Finally, the preprint containing the results from my presentation you cite is available since two days. I promised to upload it soon after the conference, but - for whatever reasons - it took nearly two month from handing in until first publication. It could be beneficial to the reader to also refer to the preprint, as things (especially the occurrence of the returning wake effect) are described in more detail there.
10.5194/wes-2023-81
Best regards,
Christian
Disclaimer: this community comment is written by an individual and does not necessarily reflect the opinion of their employer.Citation: https://doi.org/10.5194/wes-2023-109-CC1 -
AC1: 'Reply on CC1', Alessandro Bianchini, 07 Sep 2023
Dear Christian,
thank you for your comments.
It will be of course a pleasure to cite your prepint. As you can imagine, ours was submitted some weeks ago, when it was not available yet. We do believe that the studies together could really give a new and improved perspective on the problem!
Citation: https://doi.org/10.5194/wes-2023-109-AC1 -
CC2: 'Reply on AC1', Christian Schulz, 07 Sep 2023
Dear Alessandro,
good to see that we are aligned here:) Good luck with the further review process!
Best regards,
Christian
Disclaimer: this community comment is written by an individual and does not necessarily reflect the opinion of their employer.Citation: https://doi.org/10.5194/wes-2023-109-CC2
-
CC2: 'Reply on AC1', Christian Schulz, 07 Sep 2023
-
AC1: 'Reply on CC1', Alessandro Bianchini, 07 Sep 2023
-
RC1: 'Comment on wes-2023-109', Anonymous Referee #1, 21 Nov 2023
Dear authors,
thank you very much for the very well elaborated work. It is rare to have a submitted paper to be so well prepared on first submission. Thus, I only have a few very minor aspects to improve:
1. Please make sure, that all variables in the equations are explained at least once. Even if they seem very clear, in other papers, someone might use the exact same variable for something different, which also seems clear. To avoid this confusion, please define once. In most cases this is done, but not in all. So check.
2. On line 303 there seems to be a typo before “LC 2.12”
3. In the paragraph from line 317 on, please use one or two sentences to explain first what step tests are and what they are used for before you go into the discussion.
4. The sentence starting in line 320 is not really good and hard to read. Please check it.
5. On line 331 “and on the operating point in exam.” What is exam?
6. In line 340 the sentence “While not to the same extend, this consideration holds true for most of the blade.” is not very specific and clear. Either leave it away or get or precise, I would say.
7. In line 352 “This is indeed very different form a step test”, most likely it is “from a step test”, isn’t it?
8. On line 416: “BEM-based models are called to perform reliably in …”, I would always say: “BEM-based models are said to perform reliably in …”. What do you think?
9. On line 539 you write about windmill wakes – do you really mean windmills?
Otherwise, the paper opens a lot of room for specific aerodynamic discussion, which most likely cannot be finalized in one single paper. Thus, I find it good this way. But there are still open questions also coming from this paper, which will most likely need to be answered in the further discussions. I’m looking forward to it.
Citation: https://doi.org/10.5194/wes-2023-109-RC1 - AC2: 'Reply on RC1', Alessandro Bianchini, 11 Jan 2024
-
RC2: 'Comment on wes-2023-109', Martin Hansen, 08 Dec 2023
The topic of investigating the performance of BEM codes for FOWT is very interesting and of big practical importance, but not very new. A very similar study was made 10 years ago by
J.B. de Vaal et al., Effect of wind turbine surge motion on rotor thrust and induced velocity, Wind Energy (2012)
The conclusion in the submitted paper is reconfirmed, that a BEM code with a proper dynamic inflow model and an empirical Glauert model for high thrust coefficients used on FOWTs performs quite well for the exposed structural oscillations (amplitudes and frequencies). It would be really nice, if the expected range of frequencies and amplitudes for the DTU 10MW rotor (pitch and surge) was included in the paper for different foundation types. This will also show if going above possible frequencies in the PoliMi tunnel is of practical importance for a real FOWT.
The paper discusses different ways reported in the literature of how to treat the momentum equation in case of a dynamically oscillating wind turbine rotor. In basic fluid mechanics the conservation of momentum is used to determine an unknown force by keeping track of the total momentum deficit out of a control volume and including an inertia term in case of unsteadiness. The mean position of the turbine is not moving, so the control volume and velocities at the boundaries should be in a fixed frame of reference and not include the velocities of the rotor. These should in my opinion only be used when evaluating the angles of attack for the blade elements.
On page 5 and 6 the LLFVW is described as a dynamic vortex model where the vorticity is shed as vortex rings. Note that the coefficients in the Øye Dynamic Inflow model are actually calibrated from a similar dynamic ring vortex approach and the results using an unsteady BEM are therefore expected to be similar to the LLFVW output.
The axial induction velocities from the LLFVW shown in Figure 7 are very small and in the order of 0.03 m/s and compared to the inflow velocity of 4 m/s correspond to an axial induction factor of around a=0.03/4=0.008. If this is true then there is practically no induction for this case and what is the corresponding CT ?
On page 14 is reported a time constant for the dynamic wake LLFVW computations of around 3-4 seconds. In the Øye Dynamic Wake model the time constant is approximately the rotor diameter divided by the free wind speed and in the case of the UNAFLOW wind turbine should be around tau=2.4/4=0.6 seconds. That is the LLFVW model responds quite much slower to a dynamically changed force than will the Øye model. Since both the Øye dynamic wake model and LLFVW are based on a similar vortex ring model for the shed vorticity what is then the reason for this difference in time constant ? An Actuator Line simulation that through the N-S equations resolves the real physics and inertia of the wake response could be used to check these LLFVW results.
The result shown in figure 12 is interesting. Here the simulations show that the axial induction factor for a pitching motion of amplitude between 1 and 2 degrees at a frequency of 0.1Hz and at a wind speed of 5 m/s can be as high as a=4 near the blade tip, meaning that the blade will experience a velocity from behind of about 3 times the wind speed. This is estimated to occur at wave heights of 10-13 meters. Is it a realistic scenario to have a wind speed of only 5 m/s at wave heights of more than 10 meters ? And how should a Glauert correction be when the axial induction becomes so large corresponding to a thrust coefficient way above 2 ? And is it the free wind speed or the apparent wind speed taking the structural velocity into account one should use when computing the thrust coefficient CT in a BEM based model ?
It is well known that BEM becomes inaccurate for large blade deflections for a bottom fixed wind turbine. This is because the blade elements are moved away from the rotor plane where the induction is computed combined with a strong streamwise gradient of the induced velocity near the rotor. This effect can become very severe in the case of a floating wind turbine where the position of the rotor plane is continuously moving along the wind direction and how to possibly treat this in a dynamic BEM code should also be discussed in a paper like this. The question is whether the induced velocity field follows the oscillating FOWT or the rotor moves in a velocity gradient fixed in space. This effect is very important and depends on the time constants of the rotor oscillation compared to the inertia (time constants) of the flow as also discussed in the paper by J.B. de Vaal et al.
Citation: https://doi.org/10.5194/wes-2023-109-RC2 - AC3: 'Reply on RC2', Alessandro Bianchini, 11 Jan 2024
Status: closed
-
CC1: 'Comment on wes-2023-109', Christian Schulz, 07 Sep 2023
Dear Francesco and colleagues,
thanks for sharing this extensive work! It is great to see that the presence of the returning wake effect at the UNAFLOW rotor can also be observed in OLAF and even in ALM. I feel we finally make some progress in the field of unsteady aerodynamics of FOWT :) Finally, the preprint containing the results from my presentation you cite is available since two days. I promised to upload it soon after the conference, but - for whatever reasons - it took nearly two month from handing in until first publication. It could be beneficial to the reader to also refer to the preprint, as things (especially the occurrence of the returning wake effect) are described in more detail there.
10.5194/wes-2023-81
Best regards,
Christian
Disclaimer: this community comment is written by an individual and does not necessarily reflect the opinion of their employer.Citation: https://doi.org/10.5194/wes-2023-109-CC1 -
AC1: 'Reply on CC1', Alessandro Bianchini, 07 Sep 2023
Dear Christian,
thank you for your comments.
It will be of course a pleasure to cite your prepint. As you can imagine, ours was submitted some weeks ago, when it was not available yet. We do believe that the studies together could really give a new and improved perspective on the problem!
Citation: https://doi.org/10.5194/wes-2023-109-AC1 -
CC2: 'Reply on AC1', Christian Schulz, 07 Sep 2023
Dear Alessandro,
good to see that we are aligned here:) Good luck with the further review process!
Best regards,
Christian
Disclaimer: this community comment is written by an individual and does not necessarily reflect the opinion of their employer.Citation: https://doi.org/10.5194/wes-2023-109-CC2
-
CC2: 'Reply on AC1', Christian Schulz, 07 Sep 2023
-
AC1: 'Reply on CC1', Alessandro Bianchini, 07 Sep 2023
-
RC1: 'Comment on wes-2023-109', Anonymous Referee #1, 21 Nov 2023
Dear authors,
thank you very much for the very well elaborated work. It is rare to have a submitted paper to be so well prepared on first submission. Thus, I only have a few very minor aspects to improve:
1. Please make sure, that all variables in the equations are explained at least once. Even if they seem very clear, in other papers, someone might use the exact same variable for something different, which also seems clear. To avoid this confusion, please define once. In most cases this is done, but not in all. So check.
2. On line 303 there seems to be a typo before “LC 2.12”
3. In the paragraph from line 317 on, please use one or two sentences to explain first what step tests are and what they are used for before you go into the discussion.
4. The sentence starting in line 320 is not really good and hard to read. Please check it.
5. On line 331 “and on the operating point in exam.” What is exam?
6. In line 340 the sentence “While not to the same extend, this consideration holds true for most of the blade.” is not very specific and clear. Either leave it away or get or precise, I would say.
7. In line 352 “This is indeed very different form a step test”, most likely it is “from a step test”, isn’t it?
8. On line 416: “BEM-based models are called to perform reliably in …”, I would always say: “BEM-based models are said to perform reliably in …”. What do you think?
9. On line 539 you write about windmill wakes – do you really mean windmills?
Otherwise, the paper opens a lot of room for specific aerodynamic discussion, which most likely cannot be finalized in one single paper. Thus, I find it good this way. But there are still open questions also coming from this paper, which will most likely need to be answered in the further discussions. I’m looking forward to it.
Citation: https://doi.org/10.5194/wes-2023-109-RC1 - AC2: 'Reply on RC1', Alessandro Bianchini, 11 Jan 2024
-
RC2: 'Comment on wes-2023-109', Martin Hansen, 08 Dec 2023
The topic of investigating the performance of BEM codes for FOWT is very interesting and of big practical importance, but not very new. A very similar study was made 10 years ago by
J.B. de Vaal et al., Effect of wind turbine surge motion on rotor thrust and induced velocity, Wind Energy (2012)
The conclusion in the submitted paper is reconfirmed, that a BEM code with a proper dynamic inflow model and an empirical Glauert model for high thrust coefficients used on FOWTs performs quite well for the exposed structural oscillations (amplitudes and frequencies). It would be really nice, if the expected range of frequencies and amplitudes for the DTU 10MW rotor (pitch and surge) was included in the paper for different foundation types. This will also show if going above possible frequencies in the PoliMi tunnel is of practical importance for a real FOWT.
The paper discusses different ways reported in the literature of how to treat the momentum equation in case of a dynamically oscillating wind turbine rotor. In basic fluid mechanics the conservation of momentum is used to determine an unknown force by keeping track of the total momentum deficit out of a control volume and including an inertia term in case of unsteadiness. The mean position of the turbine is not moving, so the control volume and velocities at the boundaries should be in a fixed frame of reference and not include the velocities of the rotor. These should in my opinion only be used when evaluating the angles of attack for the blade elements.
On page 5 and 6 the LLFVW is described as a dynamic vortex model where the vorticity is shed as vortex rings. Note that the coefficients in the Øye Dynamic Inflow model are actually calibrated from a similar dynamic ring vortex approach and the results using an unsteady BEM are therefore expected to be similar to the LLFVW output.
The axial induction velocities from the LLFVW shown in Figure 7 are very small and in the order of 0.03 m/s and compared to the inflow velocity of 4 m/s correspond to an axial induction factor of around a=0.03/4=0.008. If this is true then there is practically no induction for this case and what is the corresponding CT ?
On page 14 is reported a time constant for the dynamic wake LLFVW computations of around 3-4 seconds. In the Øye Dynamic Wake model the time constant is approximately the rotor diameter divided by the free wind speed and in the case of the UNAFLOW wind turbine should be around tau=2.4/4=0.6 seconds. That is the LLFVW model responds quite much slower to a dynamically changed force than will the Øye model. Since both the Øye dynamic wake model and LLFVW are based on a similar vortex ring model for the shed vorticity what is then the reason for this difference in time constant ? An Actuator Line simulation that through the N-S equations resolves the real physics and inertia of the wake response could be used to check these LLFVW results.
The result shown in figure 12 is interesting. Here the simulations show that the axial induction factor for a pitching motion of amplitude between 1 and 2 degrees at a frequency of 0.1Hz and at a wind speed of 5 m/s can be as high as a=4 near the blade tip, meaning that the blade will experience a velocity from behind of about 3 times the wind speed. This is estimated to occur at wave heights of 10-13 meters. Is it a realistic scenario to have a wind speed of only 5 m/s at wave heights of more than 10 meters ? And how should a Glauert correction be when the axial induction becomes so large corresponding to a thrust coefficient way above 2 ? And is it the free wind speed or the apparent wind speed taking the structural velocity into account one should use when computing the thrust coefficient CT in a BEM based model ?
It is well known that BEM becomes inaccurate for large blade deflections for a bottom fixed wind turbine. This is because the blade elements are moved away from the rotor plane where the induction is computed combined with a strong streamwise gradient of the induced velocity near the rotor. This effect can become very severe in the case of a floating wind turbine where the position of the rotor plane is continuously moving along the wind direction and how to possibly treat this in a dynamic BEM code should also be discussed in a paper like this. The question is whether the induced velocity field follows the oscillating FOWT or the rotor moves in a velocity gradient fixed in space. This effect is very important and depends on the time constants of the rotor oscillation compared to the inertia (time constants) of the flow as also discussed in the paper by J.B. de Vaal et al.
Citation: https://doi.org/10.5194/wes-2023-109-RC2 - AC3: 'Reply on RC2', Alessandro Bianchini, 11 Jan 2024
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
587 | 190 | 23 | 800 | 14 | 15 |
- HTML: 587
- PDF: 190
- XML: 23
- Total: 800
- BibTeX: 14
- EndNote: 15
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1