{"id":"390e4d0f-3f53-4e5d-8a58-c91500f3ed34","arxiv_id":"2505.02412","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Pitch freeplay can make transonic buffet forces lock onto superharmonics of the heave frequency, driving 2:1 and 3:1 resonance limit cycles.","lead":"This paper simulates a NACA0012 wing section in transonic buffet coupled to a structural model with pitch freeplay, and finds the flow can lock onto twice or three times the heave frequency, creating large heave oscillations. The result points to a new buffet-driven limit-cycle risk for aircraft with loose hinges or linkages.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The §4.1 flutter exclusion uses a first-order spatial discretization that alters the mean flow and is not cross-checked against the second-order buffet flow; if the true flutter boundary is lower, the observed LCOs could be post-flutter rather than buffet-freeplay lock-in.","rationale":"I reviewed the paper in good faith. The central claim, that pitch freeplay induces aerodynamic lock-in to superharmonics of the heave frequency, is interesting and the paper provides supportive time histories, spectra, shock-location synchronization, and parameter sweeps. The rigid-airfoil buffet validation is reasonable for URANS (6.1% CL error, 1.29% frequency error). The strongest alternative explanation for the large LCOs is that they are a post-flutter instability, and the paper's own flutter screen (§4.1) is the only evidence against this. That screen relies on a first-order spatial discretization that removes the buffet but also changes the mean aerodynamic state; it is not validated against the second-order flow. The time-marching check in Fig. 6 uses the same approximation, so it is circular for this purpose. Thus the flutter exclusion is the load-bearing weak point. A second-order flutter computation via flow stabilization, or a growth test in the no-lock-in regime, would settle whether the LCOs are buffet-freeplay lock-in or flutter. The reader's weakest_assumption identified this same issue, and I agree. I therefore do not change the reader's conditional verdict; the concern strengthens the need for the proposed check before acceptance.","tokens_in":13979,"tokens_out":6130,"duration_ms":74217,"concrete_test":"Recompute the linear flutter speed ratio at M∞=0.72, α0=6◦ using the second-order spatial discretization with the buffet instability removed by a flow-stabilization technique (e.g., selective frequency damping), and compare V*/Vf* for kh=0.4-0.555 to the first-order results. If any case yields V*/Vf* ≥ 1, the observed LCOs are likely post-flutter. As a complementary check, run a time-marching second-order simulation at kh=0.4 with αs=0.1◦ (no-lock-in regime) and confirm that the heave response component at the heave frequency does not grow beyond the forced-buffet level.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 4.1 suppresses buffet by reducing the momentum-equation spatial accuracy to first order, then computes a linear flutter boundary via CFD-based indicial aerodynamics. Fig. 4 shows this first-order flow has a different mean shock location (at the maximum downstream position of the buffeting flow), implying substantially different aerodynamic loads. The resulting V*/Vf* values approach 0.96 at kh=0.4; the margin is thin, and the time-marching check in Fig. 6 uses the same first-order flow, so it cannot validate the flutter boundary for the actual second-order (buffet) flow. Consequently, the claim that the system operates below flutter at all tested conditions is not established. If the true flutter boundary is lower, the large-amplitude heave LCOs reported in §§4.2-4.4 would be freeplay-induced flutter LCOs rather than buffet lock-in to structural superharmonics, and the superharmonic spectral peaks would be a symptom of the nonlinear LCO rather than evidence of a buffet lock-in mechanism. This is the weakest load-bearing assumption in the paper.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents URANS simulations of a NACA0012 airfoil at M∞=0.72 and α0=6°, coupled to a two-degree-of-freedom heave-pitch structural model with pitch freeplay. The authors first validate the rigid-airfoil buffet prediction against the McDevitt-Okuno experiment, then sweep the heave frequency ratio (0.4–0.555) and freeplay angle (0°, 0.25°, 0.5°). They report that with sufficiently large freeplay, the aerodynamic force spectrum locks onto superharmonics 2kh or 3kh of the heave natural frequency while the structure oscillates at kh, producing 2:1 and 3:1 subharmonic resonances with heave limit-cycle amplitudes 16–37 times the linear baseline. The sensitivity of these phenomena to structural-to-fluid mass ratio and structural damping is also examined. The paper claims this is the first demonstration of buffet-freeplay lock-in.","tokens_in":14174,"tokens_out":6378,"duration_ms":76616,"significance":"If the reported phenomenon is confirmed, this is a novel aeroelastic interaction: freeplay-induced superharmonic excitation of the shock-buffet flow at frequency ratios well below the previously reported lock-in range, with potentially significant implications for fatigue and flight safety. The paper's strengths include the careful rigid-airfoil validation (buffet frequency within 1.29% of experiment), mesh refinement for the rigid configuration, and the use of time-frequency analysis, shock-trajectory tracking, and Lissajous curves to characterize the mechanisms. The kinetic-energy interpretation of the mass-ratio effect is a plausible physical explanation. The principal caveat is that the flutter exclusion, which is essential to the interpretation of the LCOs as buffet-driven rather than post-flutter, rests on an ad-hoc first-order spatial discretization that is not cross-validated against the buffeting flow.","major_comments":[{"comment":"The flutter exclusion relies entirely on reducing the momentum equation to first-order spatial accuracy to suppress the buffet instability, and the time-marching check in Fig. 6 uses the same first-order flow. Figure 4 shows that the first-order solution has a different mean shock location (at the maximum downstream position of the buffeting flow) than the second-order buffet flow, so the linear flutter boundary computed from this altered flow may not represent the actual buffet flow. Because the system is predicted to be at 96% of the flutter speed at the lowest heave frequency (kh=0.4) and the large-amplitude LCOs occur at or near this condition, the possibility that these LCOs are freeplay-triggered post-flutter motions rather than buffet lock-in is not eliminated. Please provide a cross-check using a buffet-stable flow that retains second-order accuracy (e.g., a flow-stabilization method) or a time-marching aeroelastic simulation with the second-order scheme at a lower velocity index to confirm that all reported LCOs are sub-flutter.","section":"§4.1, Figs. 4–6"},{"comment":"The coupled aeroelastic simulations are not subject to mesh or time-step sensitivity studies. The validation in §3 is for the rigid airfoil only; the dynamic-mesh, moving-airfoil, and freeplay-impact processes may have different resolution requirements. Since the reported amplification factors (16×, 24×, ~37×) and the abrupt transition between 3:1 and 2:1 lock-in near kh≈0.49 are central quantitative results, a grid and Δτ convergence study for at least one 2:1 case (e.g., kh=0.555, αs=0.5°) and one 3:1 case (e.g., kh=0.45, αs=0.5°) is needed to demonstrate that these observations are not numerical artifacts.","section":"§§4.2–4.5"}],"minor_comments":[{"comment":"There are several typographical errors: in §4.2, 'the structure oscillating oscillating at kh' should read 'oscillating at kh'; in the Fig. 13 caption, 'aligns wirg' should read 'aligns with'; in §5, 'signification interest' should read 'significant interest'; in §4.2, 'the the increase' should read 'the increase'.","section":"§4.2 and Fig. 13"},{"comment":"The text uses αs to denote the half freeplay gap, but the phrase 'freeplay αs=0.25°' may be ambiguous; please state explicitly that αs is the half-gap rather than the full gap, and confirm that the baseline aeroelastic results are computed with zero structural damping in both modes.","section":"Eq. (2) and §4"},{"comment":"The Short-Time Fourier Transform in Fig. 11 is used to support the beating transient discussion, but the window function, window length, and overlap are not specified; a brief note on these parameters would improve reproducibility.","section":"Fig. 11"},{"comment":"The claim of 'for the first time' is strong, especially in light of the authors' own related work reported in reference [38]. Please calibrate the novelty statement to distinguish the specific new mechanisms reported here from the earlier conference paper.","section":"Abstract and Introduction"},{"comment":"The computational cost of the 1×10^6 time-step simulations is not reported; a statement of typical wall-clock time and core count would help readers assess the practical feasibility of the study.","section":"§2.1"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a genuinely new and interesting problem, and the rigid-airfoil validation is solid. The main risk is the §4.1 flutter exclusion: the ad-hoc first-order momentum stabilization is not cross-checked against the actual buffet flow, and the reported 96% flutter margin leaves little room for error. This concern is fixable, but it is load-bearing for the interpretation of the LCOs. I would recommend major revision and suggest the authors add a corroborating flutter check (e.g., a time-marching aeroelastic simulation with second-order flow at a reduced velocity index, or a flutter boundary estimate from a second-order flow-stabilization method)."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know two things about this paper. First, it is genuinely new: the coupling of transonic shock buffet with a pitch-freeplay structural nonlinearity, and the finding that large freeplay can induce 2:1 and 3:1 superharmonic lock-in at frequency ratios below 1, is not in the cited literature. Prior work combined buffet with linear structure, or freeplay without the global buffet instability. Second, the paper's central claim is plausible but not fully established, because the authors rule out flutter using a flow with artificially increased numerical diffusion that has a different mean shock location, and they never check that approximation against the actual buffet flow.\n\nThe paper does several things well. The rigid-airfoil validation is solid: buffet frequency within 1.3% of experiment, grid convergence for ΔCL demonstrated. The parameter sweeps are informative, and the distinction between αs=0.25° (2:1 only) and αs=0.5° (3:1 also) is clean. The kinetic-energy explanation in Figure 24 is post hoc, but it is not circular and it matches the observed mass-ratio sensitivity. The authors do not hide their assumption; they explicitly call the flutter estimate \"approximate\" and recommend experimental verification.\n\nThe soft spots are proportionate to the evidence. The biggest one is Section 4.1. The first-order momentum equation changes the steady shock location to the maximum downstream position of the buffeting flow, as their own Figure 4 shows. That means the linear flutter boundary is computed about a mean flow that is not representative of the second-order buffet condition. At the lowest heave frequency ratio the system operates at 96% of this approximate flutter speed—a thin margin. The time-marching check in Figure 6 uses the same first-order flow, so it cannot validate the boundary for the actual buffet flow. If the true flutter speed is lower, the large heave LCOs could be freeplay-induced flutter rather than buffet lock-in, and the superharmonic peaks would be a symptom rather than a mechanism. I think the lock-in interpretation is quite plausible given the spectra and the distinct freeplay threshold, but the authors have not closed this door.\n\nMinor weaknesses: no aeroelastic validation (the rigid-airfoil validation does not cover the moving-mesh coupling), no mesh or time-step sensitivity for the aeroelastic cases, and no uncertainty quantification. Code and data are not released, so reproduction would require contacting the authors. These are normal for a URANS study but they do limit how definitive the claims can be.\n\nWho is this for? Researchers in transonic aeroelasticity, especially those working on buffet and freeplay, will find it valuable as a new data point and a clear hypothesis to test experimentally. It deserves a serious referee. I would send it to peer review with a request for major revision, and I would insist that the authors either cross-check the first-order flutter boundary against the second-order flow (for instance, by damping buffet through a different mechanism, or by running a few time-marching cases at the boundary with the full model) or substantially soften the claim that the system is below flutter.\n\nFor what it's worth, I would cite this paper if I were writing about buffet-freeplay interactions, but I would not rely on its flutter exclusion without more evidence.","headline":"Plausible new result on buffet-freeplay superharmonic lock-in, but the flutter exclusion rests on an unvalidated first-order flow approximation that needs a cross-check before the LCOs can be attributed to buffet rather than flutter.","tokens_in":14713,"tokens_out":1312,"would_cite":true,"duration_ms":18384,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Pitch freeplay can make transonic shock buffet lock onto heave superharmonics, producing 2:1 and 3:1 subharmonic resonances with heave limit-cycle amplitudes up to about 37 times the no-freeplay baseline.","keywords":["transonic shock buffet","freeplay nonlinearity","limit cycle oscillation","aerodynamic lock-in","superharmonic resonance","subharmonic resonance","URANS","aeroelasticity"],"falsifier":"A wind-tunnel test of a NACA 0012 section with pitch freeplay $\\alpha_s=0.5^\\circ$ at $M_\\infty=0.72$, $\\alpha_0=6^\\circ$, and $\\hat{k}_h=0.45$ should show a stable heave LCO about 37 times the no-freeplay amplitude with the lift spectrum locked onto $3\\hat{k}_h$; observing only small forced oscillations, or finding the large LCO only when the system is close to the flutter boundary, would settle against the claim.","tokens_in":1960,"feed_emoji":"✈️","tokens_out":2277,"duration_ms":107660,"temperature":0.7,"pith_summary":"This paper asks what happens when two nonlinearities meet: the self-sustained shock oscillations of transonic buffet and the bilinear stiffness of a hinge with freeplay. It finds that, for a two-degree-of-freedom heave-pitch airfoil in URANS flow, a sufficiently large pitch freeplay redistributes vibrational energy into superharmonics of the heave natural frequency, and the unsteady aerodynamics lock onto those superharmonics. The result is 2:1 and 3:1 subharmonic resonances at heave-to-buffet frequency ratios below 1, with heave limit-cycle amplitudes roughly 16–24 times (2:1) and about 37 times (3:1) the linear baseline. These are exactly the load conditions relevant to fatigue on control surfaces and empennages, and the regime identified is one that linear-structure studies would miss.","feed_headline":"Freeplay turns transonic buffet into 3:1 heave resonance","feed_subtitle":"Simulations show pitch freeplay makes shock buffet lock onto heave superharmonics, amplifying limit cycles up to 37 times.","key_machinery":"The central object is a two-degree-of-freedom heave-pitch NACA 0012 section with pitch freeplay modeled as a bilinear spring (zero restoring moment inside $\\pm\\alpha_s$), coupled to URANS aerodynamics through an embedded equation of motion. The load-bearing mechanism is superharmonic lock-in: each impact with the freeplay deadzone boundary injects energy into harmonics $p\\omega_h$; when the buffet frequency $k_{sb}$ sits close to $p\\hat{k}_h$, the aerodynamic force locks onto that superharmonic, while the structure responds at the subharmonic $\\omega_h$. Shock-location tracking and the equivalent angle of attack $\\alpha_e=\\alpha_0+\\alpha-\\dot h/U_\\infty$ are used to show how the shock and structure synchronize and why the shock dwells downstream in the 3:1 case.","core_discovery":"The paper claims that structural nonlinearity alone can move a buffet-driven airfoil from forced-harmonic-oscillator behavior into aerodynamic lock-in at superharmonics of the heave mode. With $\\alpha_s=0.25^\\circ$ freeplay and $\\hat{k}_h\\gtrsim0.515$, the lift locks onto $2\\hat{k}_h$ while the heave responds at $\\hat{k}_h$ (2:1 resonance), giving heave LCO amplification of 16 times, and 24 times with $\\alpha_s=0.5^\\circ$; with the larger freeplay and $\\hat{k}_h<0.49$, a 3:1 resonance locks onto $3\\hat{k}_h$, with heave amplification of about 37 times and pitch amplification of about 2 times. The lock-in band sits at $\\hat{k}_h=0.4$–$0.555$ with $\\hat{k}_\\alpha=0.755$, well below the previously known lock-in range near unity. The paper also reports that shock motion synchronizes with structural motion (phase about $8.6^\\circ$ with heave rate in the 2:1 case) and that during the 3:1 cycle the large heave rate pushes the equivalent angle of attack below the buffet onset angle, pinning the shock downstream. Sensitivity runs show that 2% structural damping or a doubled mass ratio suppresses the resonance, while 0.5% damping already reduces it substantially.","pith_inferences":["The same superharmonic lock-in mechanism should arise with other stiffness nonlinearities, such as cubic hardening or bilinear damping, because the energy-redistribution step does not depend on the deadzone specifically; substituting the force law in the structural equation would test this directly.","The equivalent-angle-of-attack pinning mechanism predicts that moving the freestream angle of attack closer to the buffet onset angle should narrow or shift the 3:1 lock-in band, since the large heave rate would no longer push $\\alpha_e$ below onset; this is testable with the same solver.","If these subharmonic resonances survive on a finite wing, certification fatigue spectra for tail surfaces and control-surface hinges with freeplay should include buffet-induced subharmonic loads even when the structural mode frequency lies well below the buffet frequency.","The 2-D URANS setting leaves open whether spanwise shock cells on a finite wing broaden or suppress the lock-in band; a half-span wing simulation would indicate whether the mechanism survives three-dimensional buffet."],"forward_implications":["Heave LCO amplitudes of 16–24 times (2:1) and about 37 times (3:1) the linear baseline occur with pitch freeplay at heave frequency ratios $\\hat{k}_h=0.4$–$0.555$, far below the canonical lock-in range near $\\hat{k}_h\\approx1$.","The onset of 2:1 lock-in requires only small freeplay ($0.1^\\circ<\\alpha_s<0.25^\\circ$), while 3:1 lock-in needs larger freeplay ($0.25^\\circ<\\alpha_s<0.5^\\circ$), so modest hinge wear can abruptly change the response.","Structural damping of 2% completely suppresses the 3:1 lock-in, and doubling the mass ratio also suppresses it, making the resonance far more sensitive to these parameters than linear-model lock-in.","During 3:1 lock-in the aerodynamic forces respond at $3\\omega_h$ while the structure responds at $\\omega_h$, a subharmonic resonance, with the shock pinned downstream when the equivalent angle of attack drops below the buffet onset angle.","Flutter screening indicates the test points lie below the linear flutter boundary, with the system at 96% of flutter speed at $\\hat{k}_h=0.4$, so the large LCOs are attributed to buffet-freeplay interaction rather than post-flutter motion."],"supporting_citations":[{"why":"Supplies the experimental rigid-airfoil buffet benchmark (pressure data, onset angle, buffet reduced frequency) used for validation.","marker":"[7]"},{"why":"Provides the mesh-refinement strategy and URANS buffet simulation approach that the grids here follow.","marker":"[11]"},{"why":"Furnishes the linear lock-in framework and the uncoupled 2-DOF heave-pitch equations of motion that this study extends with freeplay.","marker":"[14]"},{"why":"Gives the linear aeroelastic response baseline and the flutter-speed approximation approach used to rule out post-flutter LCO.","marker":"[31]"},{"why":"Supplies the frequency-ratio sweep and the sensitivity of lock-in to mass ratio and damping against which the nonlinear results are contrasted.","marker":"[6]"},{"why":"Identifies lock-in as coupled-mode flutter and provides the flow-stabilization idea used to obtain a linear flutter solution in buffeting conditions.","marker":"[16]"},{"why":"Records the companion result that freeplay of $0.1^\\circ$ produces no lock-in, fixing the onset threshold used to interpret the larger freeplay cases.","marker":"[38]"}],"fun_headline_variants":["Pitch freeplay induces buffet lock-in at 2:1 and 3:1 heave resonances","Freeplay drives shock buffet to lock onto heave superharmonics","Structural nonlinearity unlocks new buffet lock-in band near heave harmonics","Freeplay makes buffet lock onto heave overtones, amplifying heave 16-37x","Pitch freeplay triggers 3:1 resonance with 37x heave amplification in buffet"],"cache_read_input_tokens":16896,"weakest_assumption_plain":"The flutter check assumes that artificially increasing numerical diffusion in the flow solver removes the shock-buffet oscillations while leaving the linear flutter boundary believable; if that is wrong, the large limit cycles could be flutter rather than buffet-freeplay lock-in.","fun_headline_variants_meta":{"raw":{"variants":["Pitch freeplay induces buffet lock-in at 2:1 and 3:1 heave resonances","Freeplay drives shock buffet to lock onto heave superharmonics","Structural nonlinearity unlocks new buffet lock-in band near heave harmonics","Freeplay makes buffet lock onto heave overtones, amplifying heave 16-37x","Pitch freeplay triggers 3:1 resonance with 37x heave amplification in buffet"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000566,"raw_usage":{"total_tokens":2741,"prompt_tokens":1065,"completion_tokens":1676,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":681,"completion_tokens_details":{"reasoning_tokens":1564}},"tokens_in":681,"tokens_out":1676,"duration_ms":14242,"temperature":1.0,"reasoning_tokens":1564,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T00:52:35.208287+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A wind-tunnel test of a NACA 0012 section with pitch freeplay $\\alpha_s=0.5^\\circ$ at $M_\\infty=0.72$, $\\alpha_0=6^\\circ$, and $\\hat{k}_h=0.45$ should show a stable heave LCO about 37 times the no-freeplay amplitude with the lift spectrum locked onto $3\\hat{k}_h$; observing only small forced oscillations, or finding the large LCO only when the system is close to the flutter boundary, would settle against the claim.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the experimental rigid-airfoil buffet benchmark (pressure data, onset angle, buffet reduced frequency) used for validation."},{"cited_title":"Iovnovich, D","cited_arxiv_id":null,"evidence_quote":"Provides the mesh-refinement strategy and URANS buffet simulation approach that the grids here follow."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Furnishes the linear lock-in framework and the uncoupled 2-DOF heave-pitch equations of motion that this study extends with freeplay."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the linear aeroelastic response baseline and the flutter-speed approximation approach used to rule out post-flutter LCO."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the frequency-ratio sweep and the sensitivity of lock-in to mass ratio and damping against which the nonlinear results are contrasted."},{"cited_title":"Candon, V","cited_arxiv_id":null,"evidence_quote":"Records the companion result that freeplay of $0.1^\\circ$ produces no lock-in, fixing the onset threshold used to interpret the larger freeplay cases."}],"review_version":1}