the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Efficient derivative computation for unsteady fatigue-constrained nonlinear aero-structural wind turbine blade optimization
Abstract. Gradient-based optimization offers significant efficiency advantages for wind turbine blade design, but its application has often been limited by the cost and accuracy of finite difference derivative calculations, especially when fatigue constraints are considered. In this work, we systematically compare and evaluate four differentiation techniques, namely algorithmic differentiation, implicit differentiation, sparsity exploitation, and parallelization, to determine their effectiveness in computing accurate gradients through time-domain aero-structural simulations. By integrating these techniques with unsteady nonlinear aerodynamic and structural models, we develop software designed for accurate gradient computation. We show that combining these techniques addresses memory and runtime challenges associated with long simulations required by design load cases. Specifically, the most effective combination reduces derivative computation wall time by over an order of magnitude compared to finite differencing while maintaining superior accuracy. We demonstrate this approach in a proof-of-concept aero-structural optimization of a wind turbine blade that improves the cost of energy by 11.4 %. This comparative study establishes a viable approach for fatigue-aware blade design that balances computational efficiency with modeling accuracy.
- Preprint
(1087 KB) - Metadata XML
- BibTeX
- EndNote
Status: open (until 25 Feb 2026)
-
RC1: 'Comment on wes-2026-10', Anonymous Referee #1, 30 Jan 2026
reply
-
AC1: 'Reply on RC1', Adam Cardoza, 04 Feb 2026
reply
Thanks for your time in reading our work and feedback. We'll incorporate that fix and do another review of grammar in our revision.
Citation: https://doi.org/10.5194/wes-2026-10-AC1
-
AC1: 'Reply on RC1', Adam Cardoza, 04 Feb 2026
reply
-
RC2: 'Comment on wes-2026-10', Anonymous Referee #2, 04 Feb 2026
reply
The paper addresses a highly important problem in blade design by investigating how the efficiency advantages of gradient-based optimization can be preserved when fatigue constraints are explicitly included. Specifically, the work focuses on integrating fatigue damage calculations based on unsteady loads into a gradient-based optimization framework.
The work is well structured and demonstrates a high degree of reproducibility supported by the availability of the GitHub repository. In my opinion, no major revisions appear to be necessary. I only suggest that it may be useful to include references in which a multi-start approach is employed in optimization problems to demonstrate the robustness of the optimization results and also motivate the use of 50 starting designs (Section 3.3.2).Citation: https://doi.org/10.5194/wes-2026-10-RC2 -
AC2: 'Reply on RC2', Adam Cardoza, 04 Feb 2026
reply
Thank you for your time in reading our work and giving feedback. We agree, adding some references to support the multi-start approach will strengthen the work. We'll include some in the next revision.
Citation: https://doi.org/10.5194/wes-2026-10-AC2
-
AC2: 'Reply on RC2', Adam Cardoza, 04 Feb 2026
reply
-
RC3: 'Comment on wes-2026-10', Anonymous Referee #3, 09 Feb 2026
reply
General Remarks:
This work presents a comparison study of different methods to obtain derivatives that allow efficient gradient-based optimization in a multi-disciplinary design context. This work is well presented and its findings are very relevant for the community.The governing equations required for aeroelastic simulation of wind turbines are implemented in Julia for computational performance and the results match the state of the art code OpenFAST very well. The fresh implementation enables the use of algorithmic differentiation and enables the comparative study of FD with more efficient methods to obtain gradients.
The background, results and discussion are of high quality and logically structured. I recommend the following minor revisions to address a few small points for clarification.
Minor Comments:
Line 328: The introduction of partial and material safety factors gamma_f and gamma_m could be added here.Figure 6: The authors note that the interpretation of the optimization result is not the main focus and this is understandable. However, looking at Figure 6, an additional point comes to mind: the pitch curve indicates a considerably later initiation of the blade pitch controller in the optimized solution compared to the initial one. I would presume this contributes to the gain in AEP. Does the optimized blade design allow for this change in operation? Furthermore, it is mentioned that TSR is another DV. Perhaps the rotational speed could be added in this figure (and figure 8) on a second y axis.
Figure 7: Additional context might be helpful to interpret this figure. The induction factor shown for the initial NREL 5MW seems relatively low. Do both curves reflect the same inflow/operational condition? The wind speed could be added in the caption for reference and the authors could clarify whether the results reflect a steady (rigid) case or if they are the mean/instantaneous snapshot from an unsteady aeroelastic simulation at 10m/s (from the case described in section 2.3.3).
Line 571: “The sparse parallelized forward mode approach proves over an order of magnitude faster than traditional finite differences while delivering superior accuracy. This approach also exhibits favorable scaling characteristics as problem size increases.”
While the discussion makes this a logical conclusion, this improvement in accuracy by the sparse method is not explicitly shown in a figure. I assume that the statement in Line 537-538: “The inherent error in finite differenced derivatives exceeds the optimizer’s default tolerances for gradient noise, causing premature termination” means that the authors attempted a direct comparison of FD with the sparse method that yielded unsatisfying results in the case of FD. If so, the authors could consider to clarify this in the conclusion.
Technical Corrections:
AD introduced three times in Lines 25, 82 and 172.Line 121: ODE is not introduced.
Line 313: is beam element the correct term here or should this perhaps be blade section?
Line 412: “Table” is usually capitalized in WES. Also in Line 518.
Line 476: …”best scaling[blankspace].”
Citation: https://doi.org/10.5194/wes-2026-10-RC3
Model code and software
Cardoza2025_Efficient_aeroelastic_wind_gradients Adam Cardoza https://github.com/byuflowlab/Cardoza2025_Efficient_aeroelastic_wind_gradients
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 127 | 53 | 15 | 195 | 11 | 12 |
- HTML: 127
- PDF: 53
- XML: 15
- Total: 195
- BibTeX: 11
- EndNote: 12
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
As reviewer of the Manuscript wes-2026-10 entitled "Efficient derivative computation for unsteady fatigue-constrained nonlinear aero-structural wind turbine blade optimization", I have thoroughly reviewed the manuscript, and I would recommend addressing the below minor comment:
Except this minor comment rest of the manuscript is consolidated effectively.
I hope my critique helps the authors to improve their work and find useful in this review. Thank you!