Pith. sign in

REVIEW 3 major objections 5 minor 19 references

Connections between resonance and nonlinearity in swimming performance of a flexible heaving plate

T0 review · 3 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Pith's one-line read A viscous flexible plate has a true resonant frequency, and it predicts thrust peaks.

desk verdict A solid computational study giving the first viscous-flow resonance definition for a heaving flexible plate; the Er-based mode selection is a heuristic, but the paper's own small-amplitude simulations corroborate it. read the letter →

arxiv 1908.05704 v1 pith:ZDBRE422 submitted 2019-08-15 physics.flu-dyn

classification physics.flu-dyn
keywords flexibleheavingplatefluid-structureinteractiongloballinearstabilityresonanceswimmingperformancethrustleading-edgeseparationviscousflow
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper tries to settle whether performance peaks in finite-amplitude swimming of a flexible plate are caused by resonance, and to define resonance properly for a viscous fluid-structure system. It defines resonance through global linear stability modes of the fully coupled Navier-Stokes plus beam system, selecting the frequency of the mode that maximizes an energy ratio $E_r$; this gives an unambiguous resonant frequency that accounts for added mass and viscosity. High-fidelity nonlinear simulations at heave amplitudes from $0.001$ to $0.1$ show that at small amplitude, trailing-edge motion, thrust, and input power all peak at that frequency, while at large amplitude the peaks broaden and weaken and nonlinear effects such as leading-edge separation complicate but do not erase the resonant signature. A sympathetic reader would care because the result supplies a principled definition of resonance for viscous flexible swimmers and clarifies when linear resonance arguments can and cannot predict swimming performance.

What carries the argument

The carrying object is the global linear stability mode of the fully coupled fluid-structure system: eigenvalues $\lambda_i$ and modes $\varphi_i$ of the generalized problem $\lambda_i M \varphi_i = A \varphi_i$, where $M$ is the singular mass matrix and $A$ is the linearization of the spatially discrete Navier-Stokes and beam system about a stationary base state. Each mode is scored by the energy ratio $E_r = E_\text{plate}/(E_\text{plate}+E_\text{flow})$, computed from the eigenvector; the resonant frequency is the imaginary part of the eigenvalue with the largest $E_r$. This object does the work of defining resonance without ambiguous vacuum-beam scalings or finite-amplitude response criteria, and it provides a concrete frequency against which the nonlinear simulations are compared.

What would settle it

Compute the global modes linearized about a time-periodic base state that reproduces the $h_0 = 0.01$ heave motion; if the mode with the largest energy ratio has a frequency that differs from the stationary-base-state resonant frequency by more than the frequency resolution used here ($\Delta f = 0.1$), the claimed unambiguous resonance definition is not robust to the base state. Equivalently, re-define $E_r$ using only bending strain energy or only plate kinetic energy; if the predicted resonant frequency shifts materially, the resonance metric is an artifact of the chosen weighting.

Watch

Extended reading notes

Core claim

The central claim is that a fully coupled viscous fluid-plate system has a well-defined resonant frequency, obtained from the eigenvalues of the global linear stability problem linearized about the stationary state, and that this frequency organizes swimming performance. The resonant frequency is taken to be the imaginary part of the eigenvalue whose mode maximizes the energy ratio $E_r = E_\text{plate}/(E_\text{plate}+E_\text{flow})$, which weights modes by how much energy they put into the plate. For small heave amplitudes ($h_0 \le 0.01$) the trailing-edge amplitude, thrust, and input power all peak at or very near this frequency for every stiffness studied ($S=0.02$ to $20$), and the plate shape in the nonlinear simulations matches the corresponding eigenmode. At larger amplitude ($h_0=0.1$) the peaks broaden, weaken, and shift slightly, and nonlinear mechanisms—leading-edge separation, asymmetric wakes, aperiodic dynamics, qualitatively different input-power phase behavior—become visible, but for flexible plates the resonant frequency still marks the frequency of maximal wake circulation and near-maximal thrust. The paper concludes that resonance is a genuine, identifiable contributor to finite-amplitude swimming performance in viscous flow, not merely a linear artifact.

Load-bearing premise

The argument rests on identifying resonance with the frequency of the single global mode that maximizes the energy ratio $E_r$, computed about a stationary base state, and assuming that this same mode still governs finite-amplitude heaving.

Editorial extensions

If this is right

  • At small heave amplitudes ($h_0 \le 0.01$), peaks in trailing-edge amplitude, thrust, and input power all occur near the resonant frequency defined by the global mode with the largest energy ratio.
  • As heave amplitude increases to $h_0 = 0.1$, the resonant peaks broaden, weaken, and shift slightly in frequency, consistent with an effective damping supplied by nonlinear fluid mechanisms.
  • For all stiffnesses considered, maximum mean thrust occurs at or near the frequency of maximum wake circulation, which for flexible plates sits close to the resonant frequency.
  • At large amplitude and higher stiffness, leading-edge separation can appear without necessarily destroying performance, while for $S=2$ at $f>2.4$ it drives aperiodic dynamics and thrust peaks well away from resonance.
  • Because the resonance definition is amplitude-independent, the same linear stability framework can be carried over to other parameter regimes and to flexible-wing flight, where the plate inertia relative to the flow differs.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: The energy ratio $E_r$ weights total plate energy against total system energy, but other plausible weights, such as bending strain energy alone or flow kinetic energy alone, might select a different eigenvalue; a sensitivity sweep over energy weights would test whether the claimed unambiguous resonance frequency is robust.
  • Editorial inference: The global modes here are computed about a stationary base state; repeating the analysis about a time-periodic base state that reproduces the finite-amplitude heave cycle would separate added-mass and damping shifts from nonlinear vortex effects, and would show whether the linear resonant frequency remains the right predictor at $h_0 = 0.1$.
  • Editorial inference: The alignment of maximum mean thrust with maximum wake circulation suggests that a wake-circulation model built from linear global modes could serve as a cheap thrust predictor in regimes where leading-edge separation is weak.
  • Editorial inference: The broadening and weakening of the resonant peaks resembles classical damping; if a single amplitude-dependent damping coefficient could collapse the peak height and width across all four stiffnesses, it would make the effective-damping interpretation quantitative.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 5 minor

Summary. The paper investigates whether resonance governs the swimming performance of a flexible heaving plate in a viscous fluid. It proposes a definition of resonance for the fully coupled fluid-structure system: the imaginary part of the eigenvalue of the global linear stability mode with the largest energy ratio Er = plate energy / (plate energy + flow kinetic energy), introduced in Eq. (4.1). Global eigenvalue computations are performed at Re = 240 and M = 0.01 for a range of stiffnesses, and high-fidelity nonlinear simulations are run at heave amplitudes h0 = 0.001, 0.01, and 0.1 for S = 0.02, 0.2, 2, 20, with frequency sweeps in increments of 0.1. At the smallest amplitude, peaks in trailing-edge amplitude, thrust, and input power occur near the Er-max natural frequency; at large amplitude the peaks broaden and weaken, the input power changes qualitatively, and leading-edge separation appears for stiff plates. The authors interpret these changes as the onset of non-resonant, nonlinear mechanisms and connect their observations to earlier high-Reynolds-number studies.

Significance. If the central claim holds, the paper provides a tractable and potentially unambiguous definition of resonance for viscous fluid-structure interaction and demonstrates that, at small amplitude, linear resonance organizes thrust and input-power peaks across a broad range of stiffnesses. The paper's strengths include machine-checked convergence of the eigenvalue computations (residual below 10^-8), a grid-convergence study in Appendix A, a nonlinear solver validated against prior work, and a systematic amplitude comparison using the scalings in Table 2. The falsifiable prediction that the Er-maximal mode frequency governs small-amplitude performance peaks is clearly stated. The main unresolved issue is whether the Er-max mode selection, rather than some other linear response measure, is the mechanism behind the observed peaks.

major comments (3)
  1. [§4, Eq. (4.1)] The central definition of resonance is the eigenvalue of the mode with the largest energy ratio Er, but the paper does not establish that this mode is the one most strongly excited by the prescribed leading-edge heave or that it dominates the harmonic forced response of the linearized system (3.2). Because the linearized FSI operator may be non-normal, forced-response peaks need not coincide with the eigenvalues of any single mode; a mode with large Er could have small projection onto the forcing or weak receptivity. The small-amplitude simulations at h0 = 0.001 are effectively a forced-response experiment, but they are sampled at Δf = 0.1 and only for four stiffness values, and they do not test alternative energy weightings. I recommend adding a direct linear forced-response calculation over the same frequency range, and reporting the sensitivity of the predicted resonant frequency to the definition of the energy ratio (for example, plate kinetic energy only, bending energy only, or flow kinetic energy only).
  2. [§3, Eq. (3.2)] The base state yb about which the global modes are computed is never specified. If yb is the steady uniform flow over an undeformed straight plate, the eigenvalue problem describes free (unforced) dynamics, and the link to the harmonically forced heaving plate is not formal: one must identify how the leading-edge heave enters the linearized equations and whether the forcing projects onto the Er-max mode. Please state yb explicitly and, if the connection is to be made through the nonlinear simulations, explain how the amplitude ramp of Eq. (3.3) ensures that the small-amplitude response is a perturbation about that same base state.
  3. [§5.2.4, Figs. 14–15] For S = 20, the resonant frequency f ≈ 3.1 is close to the upper end of the frequency sweep, leaving only a few data points beyond resonance. The statements that peaks occur 'near' resonance and 'broaden and weaken' at h0 = 0.1 are therefore weakly constrained for the stiffest case, and the broad peaks at f = 2.1–2.4 in Fig. 4 lie well below the resonant frequency. Extending the frequency range above 3.5, or refining the sampling around f = 3.1, would provide a clearer test of the resonance hypothesis for this stiffness.
minor comments (5)
  1. [Table 1 and §5.2.4] Table 1 lists the frequency range as 1–3.2, but §5.2.4 states that the range is f ∈ [1, 3.5] for S = 20; please reconcile these values.
  2. [Eq. (2.3)] The sign convention in the definition CP = -CL vLE should be stated explicitly, since positive input power must correspond to heaving against the lift force.
  3. [References] The reference to 'Criesfield 1991' appears to be a typo for Crisfield; please correct the spelling.
  4. [Appendix A] The grid-convergence study reports the plate-position error for one parameter combination at a single time instant; a sentence explaining why this case is representative of the full parameter sweep would strengthen the presentation.
  5. [Fig. 2] It is not stated how many eigenvalues are shown or how the subset with 'largest growth rate' is chosen; clarifying the selection criterion would help the reader understand whether the Er-max mode is always included.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the linear stability eigenvalues are computed independently of the nonlinear performance data, and the Er-based mode selection is a stated convention rather than a fitted parameter.

full rationale

The paper's central comparison is between the natural frequencies obtained from the linearized unforced FSI system and the performance peaks observed in forced nonlinear simulations. The eigenvalues are computed by solving lambda_i M phi_i = A phi_i, where A = dr/dy at the base state (Section 3, Eq. 3.2), and no performance quantity (h_TE, C_T, C_P) is used in that computation. The nonlinear simulation peaks for h0 = 0.001, 0.01, and 0.1 are independent forward outputs that could in principle disagree with the eigenvalue frequencies; the fact that they agree near the Er-max mode is an empirical finding, not a construction. The energy ratio Er (Eq. 4.1) is a selection rule for identifying which eigenvalue is structure-dominated, but the paper does not fit Er or the chosen frequency to the thrust, input power, or trailing-edge amplitude data. Thus the claim that 'peaks in trailing-edge motion and thrust occur near the resonant frequency' is a genuine prediction from the linear analysis. Self-citations (Goza & Colonius 2017; Goza et al. 2018) are used for numerical methods and global-mode formulation, providing independent method-level support rather than importing the paper's central conclusion. The known limitations, such as possible non-normality/receptivity effects of the linearized operator and the convention of maximizing Er, are scientific modeling choices or correctness risks, not circular steps.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central comparison relies on the chosen dimensionless parameters, the linearization base state, and the energy-ratio selection rule; none of these are independently benchmarked against experiments at this parameter point, and no code or data are provided. No free parameters are fitted to the performance data in this paper.

assumptions (4)
  • domain assumption Two-dimensional Navier-Stokes flow and a geometrically nonlinear Euler-Bernoulli beam model the swimming plate at Re=240, M=0.01.
    Section 2 defines the equations; prior solver validation on flapping flags is cited, but no experiment at this exact parameter point is compared.
  • domain assumption Linearization about the stationary base state y_b yields the natural frequencies that govern heaving response.
    Section 3 linearizes Eq. (3.1) about y_b; for the finite amplitude h0=0.1, a forced nonlinear oscillation is far from this base state, so the relevance of the global modes is assumed rather than derived.
  • ad hoc to paper The eigenvalue of the mode with the largest energy ratio Er is the system resonance relevant to swimming performance.
    Section 4, Eq. (4.1), defines Er as plate energy divided by total system energy; the equal weighting of bending and kinetic energy and the choice to single out the largest-Er mode are conventions, with no sensitivity analysis provided.
  • domain assumption Heave amplitudes h0=0.001 and 0.01 are small enough that nonlinear effects are negligible, so their dynamics can be treated as linear.
    Section 5.2 compares scaled amplitudes across h0 values and infers linearity from similarity; there is no direct projection onto the linear modes to verify proportionality for all quantities.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Connections between resonance and nonlinearity in swimming performance of a flexible heaving plate." pith.science (2026). https://pith.science/paper/ZDBRE422

@misc{pith2026190805704,
  author       = {Pith},
  title        = {Pith review of: Connections between resonance and nonlinearity in swimming performance of a flexible heaving plate},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ZDBRE422}},
  note         = {Machine review of arXiv:1908.05704}
}
read the original abstract

We investigate the role of resonance in finite-amplitude swimming of a flexible flat plate in a viscous fluid. The role of resonance in performance remains unclear for two reasons: i) a lack of definition of resonance for the fully-coupled fluid-structure interaction system in a viscous flow, and ii) the presence of nonlinear effects, which makes it difficult to disentangle resonant and non-resonant mechanisms in finite-amplitude swimming. We address point i) and provide an unambiguous definition for system resonance by computing global linear stability modes of the fully-coupled fluid-structure interaction system that account for the viscous fluid, the plate, and the coupling between them. We then resolve point ii) by considering high-fidelity nonlinear simulations of systematically increased amplitude. By comparing the results for different amplitudes with one another and with the linear stability modes, we separate linear and/or resonant effects from nonlinear and/or non-resonant effects. Resonant behavior is observed over a wide range of plate stiffnesses, with peaks in trailing-edge motion and thrust occurring near the resonant frequency defined by the global linear analysis. The peaks broaden and weaken with increasing heave amplitude, consistent with an increased damping effect from the fluid. At the same time, non-resonant mechanisms are present at large heave amplitudes. The input power exhibits qualitatively different dynamics at large heave amplitudes compared to smaller heave amplitudes, where resonance dominates. Moreover, leading-edge separation is present for stiff plates at large heave amplitudes, which can drastically alter the performance characteristics from what one would expect through linear predictions.

Figures

Figures reproduced from arXiv: 1908.05704 by the authors.

Figure 1
Figure 1. A schematic of the problem setup and relevant dimensional quantities. [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. The main figures provide the first several eigenvalues with largest growth [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Left: eigenvalues associated with flow past a rigid stationary plate; [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (12 more)
Figure 4
Figure 4. Figure 4: Maximum transverse displacement of the trailing edge (left), mean thrust [PITH_FULL_IMAGE:figures/full_fig_p009_4.png]
Figure 5
Figure 5. Figure 5: Snapshots during a flapping period, with each row containing a different stiffness [PITH_FULL_IMAGE:figures/full_fig_p010_5.png]
Figure 6
Figure 6. Figure 6: Time traces of various quantities for S = 0.2, h0 = 0.1. Only the f = 3.2 curve is visible in the plot of hLE, as the prescribed leading-edge kinematics are identical. The limit-cycle behavior observed here is indicative of the plate-flow behavior except for f > 2.4 at…
Figure 7
Figure 7. Figure 7: The amplitude (top row), a, and phase shift (bottom row), ϕ, of the trailing￾edge displacement (left column), thrust (middle column), and input power (right column) versus frequency for S = 0.02. Each plot contains three sets of markers corresponding to different heave…
Figure 8
Figure 8. Figure 8: Maximum positive circulation in the plate’s wake versus frequency for [PITH_FULL_IMAGE:figures/full_fig_p013_8.png]
Figure 9
Figure 9. Figure 9: Analog of figure 7 for stiffness S = 0.2. observed for S = 0.02, the figure suggests the presence of both resonant and non-resonant mechanisms. For small heave amplitudes (h0 6 0.01), the dynamics appear to be primarily linear: the behavior of a and φ is similar for al…
Figure 10
Figure 10. Figure 10: Analog of figure 8 for stiffness S = 0.2. 5.2.3. Moderate stiffness: S = 2 [PITH_FULL_IMAGE:figures/full_fig_p015_10.png]
Figure 11
Figure 11. Figure 11: Analog of figure 7 for stiffness S = 2. The open markers for h0 = 0.1, f > 2.4 correspond to aperiodic dynamics for which this amplitude-phase analysis is not valid [PITH_FULL_IMAGE:figures/full_fig_p016_11.png]
Figure 12
Figure 12. Figure 12: Time traces of various quantities for S = 2, h0 = 0.1. Limit-cycle behavior is observed for f = 1.4, but not for f = 2.6, 3.2. This illustrates the presence of non-linear, non-resonant mechanisms for f > 2.2. Together, the results contained in figures 11–13 indicate t…
Figure 13
Figure 13. Figure 13: Analog of figure 8 for stiffness S = 2. with significantly larger thrust peaks and wake circulation than what is created by resonant mechanisms at lower frequencies. 5.2.4. Maximal stiffness: S = 20 We now turn to the stiffest case considered, S = 20 [PITH_FULL_IMAGE…
Figure 14
Figure 14. Figure 14: Analog of figure 7 for stiffness S = 20 [PITH_FULL_IMAGE:figures/full_fig_p018_14.png]
Figure 15
Figure 15. Figure 15: Analog of figure 8 for stiffness S = 20. analysis, do not necessarily indicate significant changes in performance from what one would predict through resonance-based arguments. 6. Conclusions and connections to other work In this article, we investigated the role of r…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

19 extracted references · 19 canonical work pages

  1. [1]

    Journal of Fluid Mechanics 614, 355–380

    Alben, Silas 2008 Optimal flexibility of a flapping appendage in an inviscid fluid. Journal of Fluid Mechanics 614, 355–380

  2. [2]

    Physics of Fluids 24 (5), 051901

    Alben, Silas, Witt, Charles, Baker, T Vernon, Anderson, Erik & Lauder, George V 2012 Dynamics of freely swimming flexible foils. Physics of Fluids 24 (5), 051901

  3. [3]

    Computer Methods in Applied Mechanics and Engineering 197 (25), 2131–2146

    Colonius, Tim & Taira, Kunihiko 2008 A fast immersed boundary method using a nullspace approach and multi-domain far-field boundary conditions. Computer Methods in Applied Mechanics and Engineering 197 (25), 2131–2146. Connecting resonance to swimming performance 21 Grid ∆x ∆t Smallest sub-domain size Total domain size ||χ−χf||∞ ||χf||∞ 1 0.0083 0.0004 [ −...

  4. [4]

    Criesfield, MA 1991 Non-linear finite element analysis of solids and structures, vol. 1 . Wiley, New York

  5. [5]

    Journal of Fluid Mechanics 732, 29–46

    Smits, Alexander J 2013 Scaling laws for the thrust production of flexible pitching panels. Journal of Fluid Mechanics 732, 29–46

  6. [6]

    Morgan Kaufmann

    Ebert, David S, Musgrave, F Kenton, Peachey, Darwyn, Perlin, Ken & Worley, Steven 2003 Texturing & modeling: a procedural approach. Morgan Kaufmann

  7. [7]

    Journal of Computational Physics 336, 401–411

    Goza, Andres & Colonius, Tim 2017 A strongly-coupled immersed-boundary formulation for thin elastic structures. Journal of Computational Physics 336, 401–411

  8. [8]

    Journal of Fluid Mechanics 857, 312–344

    Goza, Andres, Colonius, Tim & Sader, John Elie 2018 Global modes and nonlinear analysis of inverted-flag flapping. Journal of Fluid Mechanics 857, 312–344

Show all 19 references
  1. [9]

    Journal of Computational Physics 321, 860–873

    Goza, Andres, Liska, Sebastian, Morley, Benjamin & Colonius, Tim 2016 Accurate computation of surface stresses and forces with immersed boundary methods. Journal of Computational Physics 321, 860–873

  2. [10]

    Physics of Fluids 25 (12), 121901

    Hua, Ru-Nan, Zhu, Luoding & Lu, Xi-Yun 2013 Locomotion of a flapping flexible plate. Physics of Fluids 25 (12), 121901

  3. [11]

    Lehoucq, Richard B, Sorensen, Danny C & Yang, Chao 1998 ARPACK users’ guide: solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods . SIAM

  4. [12]

    Physics of Fluids 21 (7), 071902

    Michelin, S ´ebastien & Llewellyn Smith, Stefan G 2009 Resonance and propulsion performance of a heaving flexible wing. Physics of Fluids 21 (7), 071902

  5. [13]

    Physics of Fluids 26 (4), 041905

    Moored, KW, Dewey, PA, Boschitsch, BM, Smits, AJ & Haj-Hariri, H 2014 Linear instability mechanisms leading to optimally efficient locomotion with flexible propulsors. Physics of Fluids 26 (4), 041905

  6. [14]

    Journal of Fluid Mechanics 738, 250–267

    Quinn, Daniel B, Lauder, George V & Smits, Alexander J 2014 Scaling the propulsive performance of heaving flexible panels. Journal of Fluid Mechanics 738, 250–267

  7. [15]

    Proceedings of the National Academy of Sciences 108 (15), 5964–5969

    Ramananarivo, Sophie, Godoy-Diana, Ramiro & Thiria, Benjamin 2011 Rather than resonance, flapping wing flyers may play on aerodynamics to improve performance. Proceedings of the National Academy of Sciences 108 (15), 5964–5969

  8. [16]

    Journal of Experimental Biology 212 (1), 95–105

    Balachandran, Balakumar 2009 Influence of flexibility on the aerodynamic performance of a hovering wing. Journal of Experimental Biology 212 (1), 95–105

  9. [17]

    American Zoologist 23 (2), 709–725

    Webb, Paul W 1988 Simple physical principles and vertebrate aquatic locomotion. American Zoologist 23 (2), 709–725

  10. [18]

    Journal of Fluids and Structures 74, 385–400

    Zhang, Yang, Zhou, Chunhua & Luo, Haoxiang 2017 Effect of mass ratio on thrust production of an elastic panel pitching or heaving near resonance. Journal of Fluids and Structures 74, 385–400

  11. [19]

    Journal of Fluid Mechanics 751, 164–183

    Zhu, Xiaojue, He, Guowei & Zhang, Xing 2014 How flexibility affects the wake symmetry properties of a self-propelled plunging foil. Journal of Fluid Mechanics 751, 164–183

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.