pith. sign in

arxiv: 2606.06766 · v1 · pith:4TDCCMI4new · submitted 2026-06-04 · ⚛️ physics.flu-dyn

Vortex gust interactions with a freely-flying rigid airfoil

Bingfei Yan , Jennifer A. Franck This is my paper

Pith reviewed 2026-06-27 23:12 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords vortex gustfreely-flying airfoilheave trajectorylift coefficient modelingadded-mass effectinduced angle of attackvortex sheddinggust interaction
0
0 comments X

The pith

A model using lift from a stationary airfoil plus induced angle of attack and added-mass terms predicts the heave trajectory of a freely-flying airfoil hit by a vortex gust.

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

The paper studies the coupled lift and heave response when an isolated vortex gust strikes a rigid airfoil that can move freely in heave. It builds a model that takes the lift coefficient measured or computed for the same gust on a fixed airfoil and adds two motion-induced contributions: one from the changing angle of attack caused by the airfoil's own velocity and one from added-mass effects. Simulations show the airfoil reaches a peak heave displacement after the vortex passes and then rebounds. The added terms explain the motion well before the vortex arrives but only partly explain the rebound, which is also shaped by new vortex shedding triggered by the gust. The resulting model therefore offers a direct route to estimate free-flight heave from data already available for fixed airfoils.

Core claim

The lift coefficient on the freely-flying airfoil is obtained by augmenting the lift from the corresponding stationary-airfoil interaction with contributions from the induced angle of attack and from added mass. Direct comparison with the numerical simulations shows that these two contributions dominate the pre-impingement dynamics, whereas the post-impingement rebound is only partially recovered because gust-induced vortex shedding also acts. When the stationary-airfoil lift coefficient is supplied as input, the model reproduces the observed heave trajectory.

What carries the argument

Augmented lift model that adds induced-angle-of-attack and added-mass contributions to the stationary-airfoil lift coefficient, then integrates to obtain heave displacement.

If this is right

  • Pre-impingement heave depends mainly on the sense of vortex rotation.
  • Post-impingement rebound and the pattern of induced shedding change with angle of attack and with the vortex's transverse location.
  • The same stationary-airfoil lift data can be reused to forecast heave for a range of angles of attack and vortex positions.

Where Pith is reading between the lines

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

  • Wind-tunnel measurements on fixed airfoils could be post-processed to estimate free-flight gust responses without repeating the full dynamic simulation.
  • The decomposition may apply to problems with additional degrees of freedom such as pitch or to sequences of multiple gusts.
  • Similar splitting of lift into quasi-steady and motion-induced parts could be tested on flexible or cambered sections.

Load-bearing premise

Pre-impingement motion is governed only by induced angle of attack and added mass while post-impingement rebound is shaped by additional vortex shedding that the model does not fully include.

What would settle it

A numerical or experimental run in which the heave time history computed from the stationary lift input plus the two motion terms deviates measurably from the actual free-airfoil trajectory.

Figures

Figures reproduced from arXiv: 2606.06766 by Bingfei Yan, Jennifer A. Franck.

Figure 1
Figure 1. Figure 1: Schematic of the configuration. An upstream airfoil generates a vortex through a prescribed heaving and pitching motion, which subsequently interacts with a downstream airfoil free to heave in the transverse flow direction. where the asterisk subscript indicates a dimensional quantity, ℎ ∗ is the heave amplitude, 𝐹 ∗ 𝑦 is the vertical force, and 𝜌 𝑓 is the density of the fluid. The heave motion is then gov… view at source ↗
Figure 2
Figure 2. Figure 2: Flowchart of the tightly coupled fluid–structure interaction (FSI) algorithm. At each timestep, the fluid and structural solvers iterate until the lift force converges before advancing the simulation. (a) Overall domain (b) Intermediate mesh region (c) Near-airfoil mesh [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Meshing strategy for the two airfoil domain includes structured body-fitted mesh layers on the airfoils and a structured far-field mesh connected with an unstructured mesh between the airfoils where maximum relative mesh motion occurs. A mesh resolution study was performed to assess the sensitivity of the solution to the choice of mesh. Five meshes of varying refinement levels were considered, as summarize… view at source ↗
Figure 4
Figure 4. Figure 4: Spanwise vorticity contours showing the vortex gust and freely-flying airfoil prior to impingement. The dash-dot lines correspond to vortex trajectories [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Airfoil’s heave and lift coefficient in response to vortex interactions. Dashed lines correspond to approximate vortex impingement time. ℎ(𝑡) = ℎ0+ 2𝑀ℎ¤ 0 𝑎  1 − 𝑒 − 𝑎 2𝑀 (𝑡−𝑡0 )  + 1 𝑎 ∫ 𝑡 𝑡0  1 − 𝑒 − 𝑎 2𝑀 (𝑡−𝜏)  [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Heave motion contributions to 𝐶𝐿 for CW2.5. Dashed lines correspond to approximate impingement time. (a) Freely-flying 𝐶𝐿 vs. stationary 𝐶𝐿 (b) Freely-flying 𝐶𝐿 vs. modeled 𝐶e𝐿 (c) Heave velocity and acceleration (d) Motion induced contributions to 𝐶𝐿 [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Heave motion contributions to 𝐶𝐿 for CCW2.5. Dashed lines correspond to approximate impingement time. acceleration shown in panels (c) determine the induced angle of attack and added-mass contributions shown in panels (d). Around the first lift peak, these contributions act to reduce the magnitude of the vortex-induced lift, thereby limiting the initial acceleration of the airfoil. As the lift changes sign… view at source ↗
Figure 8
Figure 8. Figure 8: Spanwise-vorticity contours (top) and pressure contours (bottom) in the region around the airfoil after vortex impingement for CW2.5 and CW2.5-S. In the pressure panels, contour lines of the 𝑄-criterion are overlaid to identify vortical structures. 𝐶e𝐿. As a result, the duration over which the modeled 𝐶e𝐿 exceeds the freely-flying value is reduced to 𝛥𝑡 = 1. In summary, the dominant contribution arises fro… view at source ↗
Figure 9
Figure 9. Figure 9: Spanwise-vorticity contours (top) and pressure contours (bottom) in the region around the airfoil after vortex impingement for CCW2.5 and CCW2.5-S. In the pressure panels, contour lines of the 𝑄-criterion are overlaid to identify vortical structures. (a) CW0 (b) CCW0 (c) CW5 (d) CCW5 [PITH_FULL_IMAGE:figures/full_fig_p015_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Comparison between 𝐶e𝐿 and 𝐶𝐿 for various direct impingement cases. Dashed lines correspond to approximate impingement time. 0 X0-15 [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Spanwise vorticity contour of the region around the airfoil in the post-impingement stage for CW0, CCW0, CW5, and CCW5. is nearly symmetric for the 0◦ cases, whereas a much stronger asymmetry is observed for the 5◦ cases. From [PITH_FULL_IMAGE:figures/full_fig_p016_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Spanwise vorticity contour of the region around the airfoil in the post-impingement stage for CW2.5a, CW2.5b, CCW2.5a, and CCW2.5b. with upper-surface shedding becoming stronger as the angle of attack increases, whereas lower-surface shedding becomes weaker. However, at lower angles of attack this shedding is often so weak that the shed structure is barely identifiable as a vortex. It becomes most evident… view at source ↗
Figure 13
Figure 13. Figure 13: Predicted vs. actual heave trajectory for cases at different density ratio. Dashed lines correspond to instance of vortex impingement. for the CCW interactions, all remain below 0.005 across all density ratios tested. Because the relative error also decreases with density ratio, this trend is not only a consequence of the weaker heave amplitude, but also suggests that higher density ratio leads to weaker … view at source ↗
Figure 14
Figure 14. Figure 14: Schematic of the FSI validation demonstrating an elastically mounted heaving circular cylinder in incoming flow. shifts the transverse position of the vortex with minimal influence on its structure. In the present study, 𝛥𝜃 and 𝛥𝛼eff are held fixed within each group of cases corresponding to clockwise and counterclockwise vortices, and only the timing parameters are varied to control the vortex location. … view at source ↗
read the original abstract

This study numerically investigates the interaction between an isolated vortex gust and a freely-flying airfoil, introducing a theoretical framework for interpreting the coupled lift and heave response. This complex and coupled dynamics is important for modern light-weight aircraft where gusts may easily perturb the wing, generating transient changes in trajectory and attitude. Here, the freely-flying airfoil is modeled with a single degree-of-freedom in heave, and is impacted by an isolated vortex gust generated upstream. Computational results demonstrate that the freely-flying airfoil reaches a maximum heave displacement after vortex impingement and subsequently rebounds with a comparable magnitude. The lift coefficient is then modeled by augmenting the lift from a corresponding stationary airfoil interaction with motion induced contributions associated with the induced angle of attack and added-mass. A comparison of the modeled lift with the simulation data confirms that the dynamics of the airfoil before impingement is dominated by these two terms, however the rebound after impingement is only partially explained by the model since it is also influenced by the gust-induced vortex shedding. Comparisons across various parameters show that the pre-impingement motion depends primarily on vortex rotation direction, whereas the post-impingement and induced shedding patterns vary with respect to angle of attack and vortex transverse position. With the lift coefficient of the corresponding stationary airfoil interaction as an input, the model can successfully predict the heave trajectory, thus providing a mechanism to assess the dynamic motion of an airfoil from experimental/computational data of gusts interacting with fixed airfoils.

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 / 0 minor

Summary. The paper numerically studies the interaction of an isolated vortex gust with a rigid airfoil free to heave (single DOF), and proposes a reduced-order model for the lift coefficient obtained by augmenting the lift history measured on the corresponding stationary airfoil with standard induced-angle-of-attack and added-mass contributions. The model is integrated to recover the heave trajectory; the abstract states that pre-impingement dynamics are captured by these two terms while the post-impingement rebound is only partially recovered because of additional gust-induced vortex shedding.

Significance. If the quantitative accuracy of the heave prediction were demonstrated, the framework would supply a practical route to infer dynamic gust response from stationary-airfoil data, which is relevant for lightweight aircraft gust-load analysis.

major comments (2)
  1. [Abstract] Abstract: the central claim that the model 'can successfully predict the heave trajectory' rests on qualitative visual comparisons; no L2 error norms, peak-displacement errors, or time-integrated residuals between the integrated model heave and the simulated trajectory are reported, nor is any mesh-convergence or validation study against independent test cases provided.
  2. [Abstract] Abstract and §4 (model description): the post-impingement rebound is explicitly stated to be 'only partially explained' by the induced-AoA + added-mass augmentation because of unmodeled gust-induced vortex shedding; without a quantitative decomposition of the residual lift force attributable to shedding versus the two motion terms, it is impossible to determine whether the heave prediction remains within acceptable engineering tolerance after impingement.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments highlight opportunities to strengthen the quantitative support for the model's claims, which we will address in revision. We respond to each major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the central claim that the model 'can successfully predict the heave trajectory' rests on qualitative visual comparisons; no L2 error norms, peak-displacement errors, or time-integrated residuals between the integrated model heave and the simulated trajectory are reported, nor is any mesh-convergence or validation study against independent test cases provided.

    Authors: We agree that quantitative error metrics would provide a more rigorous basis for the claim of successful prediction. In the revised manuscript we will add L2 norms, peak-displacement errors, and time-integrated residuals for the heave trajectories in all presented cases. The computational mesh follows the validated resolution used in our earlier stationary-airfoil studies; we will include a brief mesh-convergence check (lift and heave sensitivity) in an appendix. The stationary-airfoil data serve as the primary validation baseline for the augmentation terms; we will clarify this linkage in the text. revision: yes

  2. Referee: [Abstract] Abstract and §4 (model description): the post-impingement rebound is explicitly stated to be 'only partially explained' by the induced-AoA + added-mass augmentation because of unmodeled gust-induced vortex shedding; without a quantitative decomposition of the residual lift force attributable to shedding versus the two motion terms, it is impossible to determine whether the heave prediction remains within acceptable engineering tolerance after impingement.

    Authors: We accept that a quantitative decomposition of the residual lift would allow a clearer assessment of engineering tolerance. The present work demonstrates partial recovery via direct overlay of modeled and simulated lift histories, showing systematic deviation only after impingement. Performing a full decomposition requires additional flow-field post-processing to isolate shedding-induced forces, which exceeds the scope of the current reduced-order modeling focus. In revision we will report the magnitude of the post-impingement residual lift relative to the motion terms, quantify the resulting heave error, and discuss its implications for tolerance; a supplementary figure illustrating the residual will be added if space permits. revision: partial

Circularity Check

0 steps flagged

No significant circularity; stationary lift coefficient is external input, augmented by standard terms

full rationale

The paper obtains the stationary-airfoil lift coefficient from separate simulations of fixed airfoils interacting with the gust and uses it as an independent external input. This is then augmented only with standard induced-angle-of-attack and added-mass contributions (standard in unsteady aerodynamics) to integrate the heave equation of motion. The resulting trajectory prediction does not reduce by the paper's own equations to any quantity fitted from the dynamic free-flight data itself. The abstract explicitly acknowledges that post-impingement rebound is only partially captured due to unmodeled gust-induced vortex shedding, confirming the model is not claimed to be complete by construction. No self-citation chain or self-definitional step is load-bearing.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard assumptions of incompressible unsteady aerodynamics and the numerical fidelity of the CFD solver for both stationary and moving cases; no new entities or fitted parameters are introduced in the abstract.

axioms (1)
  • domain assumption Lift on a moving airfoil can be decomposed into the lift from the corresponding stationary interaction plus contributions from induced angle of attack and added-mass effects.
    This decomposition is the explicit basis of the modeling approach described in the abstract.

pith-pipeline@v0.9.1-grok · 5796 in / 1467 out tokens · 30683 ms · 2026-06-27T23:12:55.749337+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

25 extracted references

  1. [1]

    Physical Review Fluids , volume=

    Gust encounters of rigid wings: Taming the parameter space , author=. Physical Review Fluids , volume=. 2020 , publisher=

    2020

  2. [2]

    AIAA Journal , volume=

    Counterclockwise vortical-gust/airfoil interactions at a transitional Reynolds number , author=. AIAA Journal , volume=. 2018 , publisher=

    2018

  3. [3]

    AIAA Journal , volume=

    Clockwise vortical-gust/airfoil interactions at a transitional Reynolds number , author=. AIAA Journal , volume=. 2018 , publisher=

    2018

  4. [4]

    Journal of Fluids and Structures , volume=

    Analysis of vortical gust impact on airfoils at low Reynolds number , author=. Journal of Fluids and Structures , volume=. 2020 , publisher=

    2020

  5. [5]

    Physics of Fluids , volume=

    Vortex dynamics during blade-vortex interactions , author=. Physics of Fluids , volume=. 2015 , publisher=

    2015

  6. [6]

    AIAA Journal , volume=

    Vortical gusts: experimental generation and interaction with wing , author=. AIAA Journal , volume=. 2019 , publisher=

    2019

  7. [7]

    Journal of Fluids and Structures , volume=

    Aeroelastic interactions and trajectory selection of vortex gusts impinging upon Joukowski airfoils , author=. Journal of Fluids and Structures , volume=. 2020 , publisher=

    2020

  8. [8]

    Aerospace Science and Technology , volume=

    Gust response of rigid and elastically mounted airfoils at a transitional Reynolds number , author=. Aerospace Science and Technology , volume=. 2018 , publisher=

    2018

  9. [9]

    Journal of Fluids and Structures , volume=

    Flow-induced vibration of a circular cylinder at limiting structural parameters , author=. Journal of Fluids and Structures , volume=. 2001 , publisher=

    2001

  10. [10]

    Measurement Science and technology , volume=

    Combining PIV, POD and vortex identification algorithms for the study of unsteady turbulent swirling flows , author=. Measurement Science and technology , volume=. 2001 , publisher=

    2001

  11. [11]

    General theory of aerodynamic instability and the mechanism of flutter , author=

  12. [12]

    arXiv preprint arXiv:2603.21336 , year=

    Generation of an isolated vortex gust through a heaving and pitching foil , author=. arXiv preprint arXiv:2603.21336 , year=

    arXiv

  13. [13]

    AIAA Journal , volume=

    Angle of attack and core size effects on transitional vortical-gust--airfoil interactions , author=. AIAA Journal , volume=. 2020 , publisher=

    2020

  14. [14]

    Physical Review Fluids , volume=

    Extreme vortex-gust airfoil interactions at Reynolds number 5000 , author=. Physical Review Fluids , volume=. 2025 , publisher=

    2025

  15. [15]

    Journal of Fluid Mechanics , volume=

    Data-driven transient lift attenuation for extreme vortex gust--airfoil interactions , author=. Journal of Fluid Mechanics , volume=. 2024 , publisher=

    2024

  16. [16]

    Journal of Fluids and Structures , volume=

    Asymmetric distributions in pressure/load fluctuation levels during blade-vortex interactions , author=. Journal of Fluids and Structures , volume=. 2017 , publisher=

    2017

  17. [17]

    Experiments in Fluids , doi=

    Interaction of quasi-two-dimensional vortical gusts with airfoils, unswept and swept wings , author=. Experiments in Fluids , doi=. 2022 , publisher=

    2022

  18. [18]

    The Aeronautical Journal , doi=

    Lift alleviation in travelling vortical gusts , author=. The Aeronautical Journal , doi=. 2023 , publisher=

    2023

  19. [19]

    AIAA Journal , doi=

    Experimental investigation of transverse and vortex gust encounters at low Reynolds numbers , author=. AIAA Journal , doi=. 2021 , publisher=

    2021

  20. [20]

    Journal of Fluid Mechanics , volume=

    Airfoil response to periodic vertical and longitudinal gusts , author=. Journal of Fluid Mechanics , volume=. 2024 , publisher=

    2024

  21. [21]

    Journal of Fluid Mechanics , volume=

    The unsteady lift of an oscillating airfoil encountering a sinusoidal streamwise gust , author=. Journal of Fluid Mechanics , volume=. 2021 , publisher=

    2021

  22. [22]

    Journal of Fluid Mechanics , volume=

    A complete second-order theory for the unsteady flow about an airfoil due to a periodic gust , author=. Journal of Fluid Mechanics , volume=. 1976 , publisher=

    1976

  23. [23]

    Experiments in Fluids , volume=

    Controlling upwards and downwards gust loads on aerofoils by pitching , author=. Experiments in Fluids , volume=. 2023 , publisher=

    2023

  24. [24]

    Airfoil in sinusoidal motion in a pulsating stream , author=

  25. [25]

    Physical Review Fluids , doi=

    Wake-foil interactions and energy harvesting efficiency in tandem oscillating foils , author=. Physical Review Fluids , doi=. 2021 , publisher=

    2021