pith. sign in

arxiv: 2605.21596 · v1 · pith:OITOHQD4new · submitted 2026-05-20 · ⚛️ physics.flu-dyn

Study of flutter instability using the actuator line method for wind energy harvesting devices

Vitor G. Kleine , Matias Herrera This is my paper

Pith reviewed 2026-05-22 08:35 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords actuator line methodflutter instabilityunsteady aerodynamicsaeroelasticitywind energy harvestingTheodorsen theorysmearing parameter
0
0 comments X

The pith

An actuator line method including pitch-rate and non-circulatory terms reproduces classical flutter predictions when the smearing parameter to chord ratio is chosen carefully.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper tests whether the actuator line method can reliably predict flutter instability for airfoils in wind energy harvesting devices. It builds a two-dimensional linear model of the ALM in harmonic motion and compares three versions against Theodorsen's theory. The basic ALM misses key unsteady effects, but adding pitch-rate and non-circulatory terms lets the method match reference flutter speeds and frequencies. This matters for designing stable energy harvesters, wind turbines, and aircraft because accurate low-cost simulations reduce reliance on expensive experiments or full CFD. Guidance emerges on tuning the smearing parameter relative to chord length to achieve agreement with classical results.

Core claim

The classical actuator line method does not accurately predict flutter, yet an ALM that includes the pitch-rate and non-circulatory terms reproduces the flutter velocity and frequency obtained from Theodorsen's theory. This holds for both a typical airfoil section and an energy-harvesting device when the ratio of the ALM smearing parameter to chord is selected appropriately. The approach replaces Theodorsen's function in the lift and moment equations with the complex function derived from the unsteady ALM response under harmonic motion.

What carries the argument

The complex function that maps the lift produced by an unsteady actuator line model in harmonic motion onto the quasi-steady lift, substituted for Theodorsen's function in the aeroelastic equations.

If this is right

  • The tuned ALM can be used to set up aeroelastic simulations of energy harvesting devices.
  • The same parameter choice supplies guidance for flutter analysis of large horizontal-axis wind turbines.
  • The method extends directly to aeroelastic studies of fixed-wing aircraft.
  • Proper selection of the smearing-to-chord ratio enables reliable flutter prediction without reverting to full CFD.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the 2D model holds in practice, the same tuning could be applied to calibrate 3D ALM simulations of finite-span devices.
  • This approach might allow faster preliminary design loops for new harvester geometries before committing to viscous or nonlinear computations.
  • Testing the method on motions that depart from pure harmonic oscillation would reveal whether the linear transfer function remains useful outside the paper's assumptions.

Load-bearing premise

The two-dimensional linear harmonic-motion model of the actuator line method captures the unsteady aerodynamic forces that actually drive flutter in three-dimensional viscous flow around a real device.

What would settle it

A side-by-side comparison of the flutter speed and frequency predicted by the tuned ALM against measurements from a wind-tunnel test of a physical energy-harvesting airfoil would confirm or refute the claimed reproduction of classical results.

Figures

Figures reproduced from arXiv: 2605.21596 by Matias Herrera, Vitor G. Kleine.

Figure 1
Figure 1. Figure 1: Notations used for the geometric parameters of the two-dimensional airfoil. Dotted line indicates [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Graphical representation of ratio of flutter frequencies (full lines) and flutter speeds (dashed [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: figure 3. These results confirm that this case is more challenging for an aerodynamic model. Neglecting [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 3
Figure 3. Figure 3: Graphical representation of ratio of flutter frequencies (full lines) and flutter speeds (dashed [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
read the original abstract

The suitability of the actuator line method (ALM) to predict flutter instability is theoretically studied by employing a two-dimensional linear model of the ALM undergoing harmonic motion. Three different analytical models of the ALM, including or not the non-circulatory and pitch-rate terms, are compared to Theodorsen's theory. First, classical methods using Theodorsen's function are employed to calculate reference values of flutter velocity and frequency. Then, the theoretical response of the ALM is predicted by replacing Theodorsen's function in the lift and aerodynamic pitching moment models with the corresponding complex function that relates the lift calculated by an unsteady ALM and the quasi-steady lift in harmonic motion. This method is applied to an airfoil typical section and to an energy harvesting device based on aeroelastic vibrations of an airfoil. From the results, it is possible to conclude that the classical ALM does not accurately predict flutter. However, we show that an ALM that considers the pitch-rate and non-circulatory terms has the capability to reproduce the results of classical methods if the ratio between ALM smearing parameter and chord is carefully chosen. These results can guide aeroelastic simulations of energy harvesting devices, large horizontal-axis wind turbines and fixed-wing aircraft.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a theoretical study of flutter instability using the actuator line method (ALM) in a two-dimensional linear model of an airfoil undergoing harmonic motion. It derives complex functions from unsteady ALM equations for different model variants (classical ALM, with pitch-rate, with non-circulatory terms) and substitutes them into classical flutter equations to compare flutter velocity and frequency against Theodorsen's theory for both a typical section and an energy harvesting device. The central result is that including pitch-rate and non-circulatory terms allows reproduction of classical results when the ALM smearing parameter to chord ratio is suitably selected.

Significance. If the result holds, this provides a parameter-selection guideline that lets an enhanced ALM recover classical Theodorsen flutter predictions inside a 2D linear framework. Such a bridge could improve the fidelity of actuator-line-based aeroelastic simulations for wind-energy harvesting devices and turbines, especially where unsteady lift and moment modeling is required.

major comments (2)
  1. [Abstract] Abstract: The claim that the results 'can guide aeroelastic simulations of energy harvesting devices, large horizontal-axis wind turbines and fixed-wing aircraft' is not supported by the presented analysis. All derivations and comparisons remain inside a 2D linear harmonic-motion model; no evidence is given that the tuned smearing-to-chord ratio remains representative once tip vortices, spanwise flow, or viscous boundary layers are introduced.
  2. [Theoretical setup] Theoretical setup (described in the comparison to Theodorsen's theory): The substitution of the ALM-derived complex function into the classical flutter equations assumes that the 2D harmonic-motion response fully captures the unsteady aerodynamic forces driving flutter. The manuscript supplies no sensitivity study or error estimate quantifying how deviations from this assumption (e.g., finite-span effects) would shift the predicted flutter speed.
minor comments (2)
  1. The notation used for the ALM-derived complex function (lift relative to quasi-steady lift) should be introduced with an explicit equation early in the text to improve readability.
  2. A brief discussion of how the chosen smearing-to-chord ratio was determined (selection rule versus trial values) would clarify the practical applicability of the guideline.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We agree that the abstract claim regarding guidance for large-scale three-dimensional applications exceeds the scope of the two-dimensional linear analysis presented. We have revised the abstract to qualify the conclusions accordingly. For the theoretical setup, we maintain that the substitution is appropriate within the stated two-dimensional framework but acknowledge the absence of a sensitivity study for deviations such as finite-span effects; we will add a clarifying discussion of this limitation.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The claim that the results 'can guide aeroelastic simulations of energy harvesting devices, large horizontal-axis wind turbines and fixed-wing aircraft' is not supported by the presented analysis. All derivations and comparisons remain inside a 2D linear harmonic-motion model; no evidence is given that the tuned smearing-to-chord ratio remains representative once tip vortices, spanwise flow, or viscous boundary layers are introduced.

    Authors: We agree that the analysis is confined to a two-dimensional linear harmonic-motion model and provides no direct evidence for the persistence of the tuned smearing-to-chord ratio under three-dimensional effects such as tip vortices, spanwise flow, or viscous boundary layers. Accordingly, we have revised the abstract to remove the references to large horizontal-axis wind turbines and fixed-wing aircraft. The statement now limits the potential guidance to the two-dimensional context and to energy-harvesting devices, where the model assumptions hold. The smearing-parameter guideline remains a useful result inside the studied framework and may serve as an initial reference for more complex simulations, but we do not claim it has been validated beyond two dimensions. revision: yes

  2. Referee: [Theoretical setup] Theoretical setup (described in the comparison to Theodorsen's theory): The substitution of the ALM-derived complex function into the classical flutter equations assumes that the 2D harmonic-motion response fully captures the unsteady aerodynamic forces driving flutter. The manuscript supplies no sensitivity study or error estimate quantifying how deviations from this assumption (e.g., finite-span effects) would shift the predicted flutter speed.

    Authors: The substitution is performed deliberately inside the same two-dimensional linear framework employed by Theodorsen's theory, thereby isolating the effect of replacing the circulatory function with the ALM-derived complex function. This permits a direct, apples-to-apples comparison of flutter velocity and frequency. We do not assert that the two-dimensional harmonic response fully represents aerodynamic forces once three-dimensional or viscous effects are present. No sensitivity study for finite-span effects is included because such effects lie outside the two-dimensional assumption of the model. In the revised manuscript we will add a short paragraph explicitly stating this limitation and noting that a three-dimensional extension would be required to quantify any shift in predicted flutter speed due to finite-span or spanwise-flow contributions. revision: partial

Circularity Check

0 steps flagged

Derivation is self-contained parametric comparison of ALM response to Theodorsen within shared 2D linear framework

full rationale

The paper first computes reference flutter velocity and frequency using classical Theodorsen theory on the aeroelastic equations. It then derives a complex function directly from the unsteady ALM equations (including or excluding pitch-rate and non-circulatory terms) that relates ALM lift to quasi-steady lift under harmonic motion. This function is substituted into the same lift and pitching-moment expressions to obtain ALM-based flutter predictions. The smearing-to-chord ratio is varied explicitly as a model parameter to identify values that make the ALM flutter results match the classical references. This is a direct model-to-model equivalence check inside the identical 2D linear harmonic setup, not a reduction of any output to its inputs by construction, not a fitted prediction presented as independent, and not reliant on self-citation chains. The derivation therefore remains self-contained.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of the 2D linear unsteady aerodynamic model, the assumption that Theodorsen's function is the correct reference, and the existence of a suitable smearing parameter that makes the ALM transfer function match it.

free parameters (1)
  • smearing parameter to chord ratio
    Chosen so that the ALM-derived complex lift function reproduces Theodorsen results; the paper states it must be 'carefully chosen' but does not derive its value from first principles.
axioms (2)
  • domain assumption The flow is two-dimensional, inviscid, and the airfoil motion is purely harmonic.
    Invoked throughout the linear model construction and comparison to Theodorsen theory.
  • standard math Theodorsen's function provides the exact reference lift and moment for the unsteady airfoil problem.
    Used as the benchmark against which all ALM variants are judged.

pith-pipeline@v0.9.0 · 5739 in / 1467 out tokens · 43818 ms · 2026-05-22T08:35:53.584657+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

23 extracted references · 23 canonical work pages

  1. [1]

    Fluids Eng.124 393–399

    Sorensen J N and Shen W Z 2002 Numerical modeling of wind turbine wakesJ. Fluids Eng.124 393–399

    work page 2002

  2. [2]

    Technical University of Denmark Ph.D

    Mikkelsen R 2003Actuator disc methods applied to wind turbines. Technical University of Denmark Ph.D. thesis PhD Thesis

    work page

  3. [3]

    Series A, Mathematical, physical, and engineering sciences37320140071

    Sørensen J N, Mikkelsen R F, Henningson D S, Ivanell S, Sarmast S and Andersen S J 2015 Simu- lation of wind turbine wakes using the actuator line techniquePhilosophical transactions. Series A, Mathematical, physical, and engineering sciences37320140071

    work page 2015

  4. [4]

    Taschner E, Deskos G, Kuhn M B, van Wingerden J W and Mart´ ınez-Tossas L A 2025 Unsteady aerodynamic loads on pitching aerofoils represented by gaussian body force distributionsJournal of Fluid Mechanics1024A15

    work page 2025

  5. [5]

    Alva E, Kleine V G and Cavalieri A V 2026 On the applicability of the actuator line method for unsteady aerodynamicsJournal of Fluid Mechanics1028A7

    work page 2026

  6. [6]

    Theodorsen T 1935 General theory of aerodynamic instability and the mechanism of flutter Tech. Rep. 496 NACA

    work page 1935

  7. [7]

    Trigaux F, Chatelain P and Winckelmans G 2024 Investigation of blade flexibility effects on the loads and wake of a 15MW wind turbine using a flexible actuator line methodWind Energy Science 91765–1789

    work page 2024

  8. [8]

    Hodgson E L, Grinderslev C, Meyer Forsting A R, Troldborg N, Sørensen N N, Sørensen J N and Andersen S J 2022 Validation of aeroelastic actuator line for wind turbine modelling in complex flowsFrontiers in Energy Research10864645

    work page 2022

  9. [9]

    Chetan M, Yao S and Griffith D T 2022 Flutter behavior of highly flexible blades for two-and three-bladed wind turbinesWind Energy Science71731–1751

    work page 2022

  10. [10]

    Hayat K, de Lecea A G M, Moriones C D and Ha S K 2016 Flutter performance of bend–twist coupled large-scale wind turbine bladesJournal of Sound and Vibration370149–162

    work page 2016

  11. [11]

    Zhuang Y and Yuan G 2024 Study on coupled mode flutter parameters of large wind turbine blades Scientific Reports1412804

    work page 2024

  12. [12]

    Erturk A, Vieira W G R, De Marqui C and Inman D J 2010 On the energy harvesting potential of piezoaeroelastic systemsApplied physics letters96

    work page 2010

  13. [13]

    Li Z, Zhou S and Yang Z 2022 Recent progress on flutter-based wind energy harvestingInternational Journal of Mechanical System Dynamics282–98

    work page 2022

  14. [14]

    Fung Y C 2008An introduction to the theory of aeroelasticity(Courier Dover Publications)

    work page

  15. [15]

    Kleine V G, Hanifi A and Henningson D S 2023 Simulating airplane aerodynamics with body forces: actuator line method for nonplanar wingsAIAA Journal612048–2059

    work page 2023

  16. [16]

    Bisplinghoff R L, Ashley H and Halfman R L 1996Aeroelasticity(Courier Dover Publications)

    work page

  17. [17]

    thesis Technical University of Denmark

    Troldborg N 2009Actuator line modeling of wind turbine wakesPh.D. thesis Technical University of Denmark

    work page

  18. [18]

    Shives M and Crawford C 2013 Mesh and load distribution requirements for actuator line cfd simu- lationsWind Energy161183–1196 9

    work page 2013

  19. [19]

    Kleine V G, Hanifi A and Henningson D S 2023 Non-iterative vortex-based smearing correction for the actuator line methodJournal of Fluid Mechanics961A29

    work page 2023

  20. [20]

    thesis KTH Royal Institute of Technology

    Kleine V G 2022On stability of vortices and vorticity generated by actuator linesPh.D. thesis KTH Royal Institute of Technology

    work page

  21. [21]

    Forsythe J R, Lynch E, Polsky S and Spalart P 2015 Coupled flight simulator and CFD calculations of ship airwake using kestrel53rd AIAA Aerospace Sciences Meetingp 0556

    work page 2015

  22. [22]

    Mart´ ınez-Tossas L A, Churchfield M J and Meneveau C 2017 Optimal smoothing length scale for actuator line models of wind turbine blades based on gaussian body force distributionWind Energy 201083–1096

    work page 2017

  23. [23]

    Katz J and Plotkin A 1991Low-speed aerodynamics: From wing theory to panel methods(McGraw- Hill Inc) 10

    work page