{"id":"531afdef-a4e2-44fb-912a-3ba7d8266b2b","arxiv_id":"1908.05704","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Global linear stability modes define the resonance of a heaving flexible plate in viscous flow; nonlinear simulations show this resonance controls small-amplitude performance but nonlinear flow effects weaken and complicate it at large amplitudes.","lead":"This paper computes the natural oscillation modes of a flexible plate flapping in a viscous fluid, then runs simulations at increasing motion amplitudes to show when resonance controls thrust and when nonlinear flow effects take over. It offers a way to define resonance for fluid-structure interaction and explains why flexible swimmers and flapping wings respond differently at large and small amplitudes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Resonance definition depends on an untested Er-based mode selection; a direct linear forced-response calculation is needed to confirm the Er-max mode governs the observed performance peaks.","rationale":"The paper's central contribution is an unambiguous resonance definition via global linear modes. The selection of the mode with the largest Er is physically motivated but not derived. Because the linearized FSI operator may be non-normal, eigenvalues alone need not determine the peak of a harmonically forced response; receptivity and projection onto the forcing also matter. The small-amplitude nonlinear simulations provide indirect evidence that the Er-max mode is the one that matters, but they are sampled at 0.1 frequency increments and only for four stiffnesses, so they do not establish that Er is the correct selector or that another selector would fail. The forced linear response test is the direct way to settle this: it uses the same operator and removes the Er heuristic entirely. The reader's conditional verdict already flags this assumption, so I do not change the verdict; the concrete test would either validate the definition or require the authors to replace the Er criterion with a response-based one. The concern is about the argument's load-bearing step, not about the authors' conduct, and no ad hominem is intended.","tokens_in":15151,"tokens_out":5700,"duration_ms":61411,"concrete_test":"Using the same discretized linearized operator A from Eq. (3.2) and the same grid/time-step family used in Appendix A, compute the steady harmonic response to the linearized leading-edge heave forcing for each stiffness S = 0.02, 0.2, 2, and 20: solve (iωM - A)y = f_heave(ω) over the frequency sweep f in [1, 3.2] in increments of 0.1, and evaluate the response amplitudes of trailing-edge displacement, thrust, and input power as in Table 2. Locate the peak frequency of each response curve. Then compare those peak frequencies with the imaginary part of the eigenvalue having the largest Er (Sec. 4). If any response peak differs from the Er-max eigenvalue frequency by more than one frequency increment (0.1) for any stiffness, the Er selection does not predict the forced linear resonance and the claim of an unambiguous resonance definition is undermined.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that the natural frequency of the fully coupled FSI system, identified as the imaginary part of the global mode with the largest energy ratio Er (Eq. 4.1), is the resonance frequency that organizes swimming-performance peaks. This selection is a convention: Er measures the fraction of modal energy in the plate, but nothing in the paper shows that this mode is the one with the largest response to the imposed leading-edge heave, nor that it is the mode dominating thrust and input power. Because the linearized FSI operator may be non-normal, the harmonic response amplitude can peak at frequencies different from the eigenvalues of individual modes; an eigenvalue with large Er but small projection onto the forcing or weak receptivity would not be the resonance seen in hTE, CT, and CP. The small-amplitude nonlinear simulations (h0 = 0.001, 0.01) do show peaks near the Er-max frequency for the four stiffnesses, so the convention is not contradicted, but those simulations are themselves the forced-response experiment, and they are sampled at 0.1 frequency increments for only four stiffness values. The paper does not test whether a different energy weight, a receptivity-weighted mode measure, or a resolvent-based response measure would change the predicted frequency. Since Section 4 and the abstract claim an 'unambiguous' definition, this untested selection rule is the load-bearing step: if the Er-max mode is not the most receptive or most forced mode, the resonance frequency could shift, and the alignment of performance peaks in Figs. 7-15 could be an artifact of the chosen metric rather than a property of the FSI system.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":15397,"tokens_out":6358,"duration_ms":62131,"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":[{"comment":"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).","section":"§4, Eq. (4.1)"},{"comment":"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.","section":"§3, Eq. (3.2)"},{"comment":"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.","section":"§5.2.4, Figs. 14–15"}],"minor_comments":[{"comment":"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.","section":"Table 1 and §5.2.4"},{"comment":"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.","section":"Eq. (2.3)"},{"comment":"The reference to 'Criesfield 1991' appears to be a typo for Crisfield; please correct the spelling.","section":"References"},{"comment":"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.","section":"Appendix A"},{"comment":"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.","section":"Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"This is a competent and well-executed study, and the proposed resonance definition is a worthwhile contribution. My main reservation is the Er-based mode selection: it is a convention, and the manuscript does not yet demonstrate that this mode governs the forced response. I believe this can be addressed within the scope of the paper by adding a linear forced-response analysis and sensitivity tests. I would not accept the paper in its current form, but the concern is fixable and the central claim remains plausible."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"At bottom: this is a careful, well-executed computational study that gives the first viscous-flow definition of resonance for a heaving flexible plate, by computing global linear modes of the fully coupled FSI system. The real contribution is the combination of that eigenvalue analysis with a systematic amplitude sweep, showing that at small amplitudes the resonant mode predicts peaks in trailing-edge amplitude, thrust, and input power, and that at large amplitudes these peaks broaden and weaken while non-resonant mechanisms (leading-edge separation, aperiodic shedding) take over. That separation of linear/resonant from nonlinear/non-resonant effects is genuinely useful. The grid convergence study in Appendix A is real evidence, the eigenvalue residuals are small, and the nonlinear solver was validated in prior published work. The numerical load-bearing is solid.\n\nThe soft spots are real but mostly minor. The Er selection rule is a heuristic: they pick the mode with the largest plate-energy fraction as \"the\" resonance frequency. Nothing in the linear analysis proves that mode is the most receptive to the imposed heave forcing; in principle non-normality could make the forced response peak elsewhere. That said, the paper itself supplies the forced-response experiment: at h0 = 0.001, the nonlinear simulations show peaks in hTE, CT, and CP at the Er-max frequency for all four stiffnesses. So the convention is empirically corroborated, at least for the parameters tested. Calling the definition \"unambiguous\" in the abstract and Section 4 is a bit strong, given the heuristic and the 0.1 frequency resolution, but the central claim holds.\n\nThe bigger limitations are scope and reproducibility: one Reynolds number (240), one mass ratio, two dimensions, no released code or data. That keeps the significance moderate rather than high. The connections to higher-Reynolds experimental work are drawn carefully, with appropriate caveats.\n\nWho this is for: anyone working on flexible-body propulsion, bioinspired swimming, or FSI stability. I'd bring it to a reading group and would cite it. It deserves a serious referee. I'd recommend sending it to JFM after minor revision, mainly softening \"unambiguous\" and adding a sentence or two on why Er is the appropriate energy weight, ideally with a quick sensitivity check.","headline":"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.","tokens_in":15942,"tokens_out":2050,"would_cite":true,"duration_ms":21515,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A viscous flexible plate has a true resonant frequency, and it predicts thrust peaks.","keywords":["flexible heaving plate","fluid-structure interaction","global linear stability","resonance","swimming performance","thrust","leading-edge separation","viscous flow"],"falsifier":"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.","tokens_in":14933,"feed_emoji":"🐟","tokens_out":5642,"duration_ms":50795,"temperature":0.7,"pith_summary":"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.","feed_headline":"Resonant mode predicts a heaving plate's thrust peaks","feed_subtitle":"Global linear modes fix a resonant frequency for a viscous heaving plate; small-amplitude peaks align, large-amplitude peaks broaden.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the inviscid definition of resonance via eigenvalues of the fully coupled fluid-structure system that this paper extends to viscous flow.","marker":"Michelin & Llewellyn Smith (2009)"},{"why":"Establishes in the inviscid setting that actuation at the resonant frequency gives peaks in thrust and input power at linear order, the baseline this paper generalizes to $Re=240$.","marker":"Floryan & Rowley (2018)"},{"why":"Provides the global-mode formulation for a flexible plate in viscous flow that is used here to compute the eigenvalues and eigenvectors.","marker":"Goza et al. (2018)"},{"why":"Supplies the strongly coupled immersed-boundary fluid-structure interaction solver used for the nonlinear simulations.","marker":"Goza & Colonius (2017)"},{"why":"Describes the iterative eigenvalue algorithm used to compute the global modes of the large generalized eigenproblem.","marker":"Lehoucq et al. (1998)"},{"why":"Provides higher-Reynolds-number experimental evidence of thrust and efficiency peaks near resonance, used as a point of comparison.","marker":"Quinn et al. (2014)"},{"why":"Offers an experimental account of nonlinear fluid damping and superharmonic resonance at large amplitudes, compared with the observed broadening and weakening of peaks.","marker":"Ramananarivo et al. (2011)"},{"why":"Documents wake asymmetry and leading-edge vortex effects in flexible foils, used to contextualize the separated-flow regimes at large amplitude.","marker":"Zhu et al. (2014)"},{"why":"Argues that optimal efficiency occurs near the wake's natural frequency, connected here to the alignment of thrust with maximal wake circulation.","marker":"Moored et al. (2014)"}],"fun_headline_variants":["Global linear modes define resonance for viscous heaving plates","Linear modes predict thrust peaks at low heave amplitude","Resonance from coupled modes predicts heaving plate thrust","Nonlinearity weakens and shifts resonance peaks at large amplitude"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Global linear modes define resonance for viscous heaving plates","Linear modes predict thrust peaks at low heave amplitude","Resonance from coupled modes predicts heaving plate thrust","Nonlinearity weakens and shifts resonance peaks at large amplitude"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001543,"raw_usage":{"total_tokens":6244,"prompt_tokens":1093,"completion_tokens":5151,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":709,"completion_tokens_details":{"reasoning_tokens":5085}},"tokens_in":709,"tokens_out":5151,"duration_ms":31630,"temperature":1.0,"reasoning_tokens":5085,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:06:11.187606+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"Physics of Fluids 21 (7), 071902","cited_arxiv_id":null,"evidence_quote":"Supplies the inviscid definition of resonance via eigenvalues of the fully coupled fluid-structure system that this paper extends to viscous flow."},{"cited_title":"Journal of Fluid Mechanics 857, 312–344","cited_arxiv_id":null,"evidence_quote":"Provides the global-mode formulation for a flexible plate in viscous flow that is used here to compute the eigenvalues and eigenvectors."},{"cited_title":"Journal of Computational Physics 336, 401–411","cited_arxiv_id":null,"evidence_quote":"Supplies the strongly coupled immersed-boundary fluid-structure interaction solver used for the nonlinear simulations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the iterative eigenvalue algorithm used to compute the global modes of the large generalized eigenproblem."},{"cited_title":"Journal of Fluid Mechanics 738, 250–267","cited_arxiv_id":null,"evidence_quote":"Provides higher-Reynolds-number experimental evidence of thrust and efficiency peaks near resonance, used as a point of comparison."},{"cited_title":"Proceedings of the National Academy of Sciences 108 (15), 5964–5969","cited_arxiv_id":null,"evidence_quote":"Offers an experimental account of nonlinear fluid damping and superharmonic resonance at large amplitudes, compared with the observed broadening and weakening of peaks."},{"cited_title":"Journal of Fluid Mechanics 751, 164–183","cited_arxiv_id":null,"evidence_quote":"Documents wake asymmetry and leading-edge vortex effects in flexible foils, used to contextualize the separated-flow regimes at large amplitude."},{"cited_title":"Physics of Fluids 26 (4), 041905","cited_arxiv_id":null,"evidence_quote":"Argues that optimal efficiency occurs near the wake's natural frequency, connected here to the alignment of thrust with maximal wake circulation."}],"review_version":1}