the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
| 31 Jul 2024
Characterization of vortex shedding regimes and lock-in response of a wind turbine airfoil with two high-fidelity simulation approaches
Abstract. In this work, the Vortex-Induced Vibrations (VIV) phenomenon affecting a wind turbine airfoil section is analysed with two numerical approaches, a two-dimensional (2D) setup of the airfoil, simulated using the Unsteady Reynolds-averaged Navier–Stokes equations, and a three-dimensional (3D) setup with a span-to-chord aspect ratio of 1, employing the Delayed Detached Eddy Simulation model. A constant inflow velocity normal to the airfoil chord is considered, for a Reynolds number around 2 × 106. The only structural degree of freedom is the airfoil chordwise displacement. As a reference, simulations of the static airfoil are also performed. The 3D static simulation lift coefficient is shown to have intermittent periods of very different characteristics, including different Strouhal numbers. The VIV simulations are performed at different inflow velocities to cover the lock-in range. To make the lock-in range non-dimensional, a single Strouhal number is chosen for the 3D case, such that the non-dimensional lock-in ranges predicted by both approaches coincide. This Strouhal number is 5 % higher than the 2D Strouhal number and 14 % lower than the one previously reported for the same 3D airfoil setup. Inside the lock-in range, the 2D and 3D approaches predict a similar VIV development, characterized by the lift coefficient standard deviation, the mean drag coefficient and the airfoil vibration amplitude growth rate. These results are supported by the common hypothesis that the three-dimensional vortex shedding coherence increases when the body undergoes large and growing motions, becoming similar to a 2D case.
- Preprint
(18419 KB) - Metadata XML
- BibTeX
- EndNote
Status: closed
- RC1: 'Comment on wes-2024-92', Anonymous Referee #1, 22 Aug 2024
-
RC2: 'Comment on wes-2024-92', Anonymous Referee #2, 04 Sep 2024
The topic and the activity are very interesting and worthy of investigation. The authors performed a complete and in-depth analysis and the results are very well presented with many details.
The reviewer has few major concerns regarding the numerical methodology and the analysis of the results:
- The radius of the domain seems to be small, being just 30 times the airfoil chord. It was demonstrated that a 2D analysis needs more than 100 chords (preferably, 500 chords) to avoid any blockage effects (https://turbmodels.larc.nasa.gov/naca0012_val.html). How do the authors justify this choice?
- How was the 1c spanwise height of the domain chosen?
- URANS-2D FSI simulations are run for 1300 non-dimensional time units each, while the DDES-3D simulations are run for 300 non-dimensional time units. How do the authors can demonstrate that the length simulated intervals are sufficient for capturing all the relevant features and the collected data are statistically significant?
- It is not clear which is the physical phenomena associated to the St=0.132 that leads to comparable range of lock-in between the two approaches
- Finally, the reviewer suggests to specify in the abstract that the airfoil is at 90° incidence
Citation: https://doi.org/10.5194/wes-2024-92-RC2 -
RC3: 'Comment on wes-2024-92', Anonymous Referee #3, 24 Sep 2024
The authors present a very interesting, detailed and well designed numerical analysis to investigate the vortex shedding regimes and lock-in response of a wind turbine airfoil. To this aim two simulation approaches with different levels of fidelity are proposed (URANS 2D and DDES 3D).
Although the considered configuration is very simple, in my opinion the paper proposes interesting findings, among which the novel metrics for the definition of the lock-in region.
The paper is well written, although there are some minor technical corrections and suggestions for improvements that I have proposed in the following.
- Abstract: lines from 7 to the end are not really suited for an abstract. I suggest moving them elsewhere and put more focus on the most important scientific novelties and findings of the work.
- Introduction: I appreciated the variety and number of cited works. In order to better position their work within the proposed literature survey, I suggest to the authors to better (and explicitly) mention the novelties and findings of this work which had not been touched (or touched differently) by previous works.
- Section 2.3: to better compare computational costs of URANSE and DDES solvers, a table would be more suitable than many lines of text. Moreover, the choice of the quantities of interest for the analysis (max(x), std dev of Cl and mean Cd) should be better explained and motivated. Finally, lines from 202 to 214 (but in general the whole Section 2.3) would be better positioned at the beginning of the numerical results section.
- Section 3.1.1: one of the reason for the different frequency content observed for Cl in URANSE and DDES simulations should be the number of vortical structures modelled in the different solvers. I think that this aspect should be better investigated by the authors. In this regard, the comparative analysis of flowfield snapshots from both simulations is needed.
As a further general comment, the authors use an in-house solver which is mentioned to be well validated. Nevertheless, often in the paper they cite results from other authors which deal with a similar configuration (e.g. lines 268 onwards for a flat plate) but the authors don't investigate those cases with their solver. Another example is related to the oscillating cylinder case which is mentioned in the introduction. The authors should put some effort in finding in the literature a test case (either numerical or experimental) which is relevant (or simply similar, e.g. a pitching airfoil undergoing dynamic stall) to the VIV problem and which they can address with the proposed solver. This validation task is, in my opinion, fundamental to improve the value of this work. Finally, when mentioning results which can support their present findings, I suggest the authors not simply mentioning them but also including them in the paper by acknowledging the authors. For instance, the work by Lian which is cited in lines 489 and further. Moreover, this specific work would be suitable as a validation test for the present solvers even if the authors mention that Reynolds numbers of Lian case are much smaller than those typical of wind turbine airfoil flows.
My indication is to ACCEPT the paper only after MINOR REVISIONS following the comments listed above.Citation: https://doi.org/10.5194/wes-2024-92-RC3 -
AC1: 'Comment on wes-2024-92', Ricardo Fernandez-Aldama, 23 Oct 2024
Dear reviewers,
Many thanks for your feedback and valuable comments.
Our reply to all of your comments is in the supplement. This document shows the reviewers' comments, followed by our answers in dark blue font color.
Where page numbers and line numbers are given, they correspond to the revised manuscript with tracked changes, to make them easy to find. The Figure and Table numbers also correspond to the revised manuscript.
Best regards,
The authors
Status: closed
- RC1: 'Comment on wes-2024-92', Anonymous Referee #1, 22 Aug 2024
-
RC2: 'Comment on wes-2024-92', Anonymous Referee #2, 04 Sep 2024
The topic and the activity are very interesting and worthy of investigation. The authors performed a complete and in-depth analysis and the results are very well presented with many details.
The reviewer has few major concerns regarding the numerical methodology and the analysis of the results:
- The radius of the domain seems to be small, being just 30 times the airfoil chord. It was demonstrated that a 2D analysis needs more than 100 chords (preferably, 500 chords) to avoid any blockage effects (https://turbmodels.larc.nasa.gov/naca0012_val.html). How do the authors justify this choice?
- How was the 1c spanwise height of the domain chosen?
- URANS-2D FSI simulations are run for 1300 non-dimensional time units each, while the DDES-3D simulations are run for 300 non-dimensional time units. How do the authors can demonstrate that the length simulated intervals are sufficient for capturing all the relevant features and the collected data are statistically significant?
- It is not clear which is the physical phenomena associated to the St=0.132 that leads to comparable range of lock-in between the two approaches
- Finally, the reviewer suggests to specify in the abstract that the airfoil is at 90° incidence
Citation: https://doi.org/10.5194/wes-2024-92-RC2 -
RC3: 'Comment on wes-2024-92', Anonymous Referee #3, 24 Sep 2024
The authors present a very interesting, detailed and well designed numerical analysis to investigate the vortex shedding regimes and lock-in response of a wind turbine airfoil. To this aim two simulation approaches with different levels of fidelity are proposed (URANS 2D and DDES 3D).
Although the considered configuration is very simple, in my opinion the paper proposes interesting findings, among which the novel metrics for the definition of the lock-in region.
The paper is well written, although there are some minor technical corrections and suggestions for improvements that I have proposed in the following.
- Abstract: lines from 7 to the end are not really suited for an abstract. I suggest moving them elsewhere and put more focus on the most important scientific novelties and findings of the work.
- Introduction: I appreciated the variety and number of cited works. In order to better position their work within the proposed literature survey, I suggest to the authors to better (and explicitly) mention the novelties and findings of this work which had not been touched (or touched differently) by previous works.
- Section 2.3: to better compare computational costs of URANSE and DDES solvers, a table would be more suitable than many lines of text. Moreover, the choice of the quantities of interest for the analysis (max(x), std dev of Cl and mean Cd) should be better explained and motivated. Finally, lines from 202 to 214 (but in general the whole Section 2.3) would be better positioned at the beginning of the numerical results section.
- Section 3.1.1: one of the reason for the different frequency content observed for Cl in URANSE and DDES simulations should be the number of vortical structures modelled in the different solvers. I think that this aspect should be better investigated by the authors. In this regard, the comparative analysis of flowfield snapshots from both simulations is needed.
As a further general comment, the authors use an in-house solver which is mentioned to be well validated. Nevertheless, often in the paper they cite results from other authors which deal with a similar configuration (e.g. lines 268 onwards for a flat plate) but the authors don't investigate those cases with their solver. Another example is related to the oscillating cylinder case which is mentioned in the introduction. The authors should put some effort in finding in the literature a test case (either numerical or experimental) which is relevant (or simply similar, e.g. a pitching airfoil undergoing dynamic stall) to the VIV problem and which they can address with the proposed solver. This validation task is, in my opinion, fundamental to improve the value of this work. Finally, when mentioning results which can support their present findings, I suggest the authors not simply mentioning them but also including them in the paper by acknowledging the authors. For instance, the work by Lian which is cited in lines 489 and further. Moreover, this specific work would be suitable as a validation test for the present solvers even if the authors mention that Reynolds numbers of Lian case are much smaller than those typical of wind turbine airfoil flows.
My indication is to ACCEPT the paper only after MINOR REVISIONS following the comments listed above.Citation: https://doi.org/10.5194/wes-2024-92-RC3 -
AC1: 'Comment on wes-2024-92', Ricardo Fernandez-Aldama, 23 Oct 2024
Dear reviewers,
Many thanks for your feedback and valuable comments.
Our reply to all of your comments is in the supplement. This document shows the reviewers' comments, followed by our answers in dark blue font color.
Where page numbers and line numbers are given, they correspond to the revised manuscript with tracked changes, to make them easy to find. The Figure and Table numbers also correspond to the revised manuscript.
Best regards,
The authors
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 349 | 103 | 257 | 709 | 18 | 15 |
- HTML: 349
- PDF: 103
- XML: 257
- Total: 709
- BibTeX: 18
- EndNote: 15
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1