{"id":"67fef4e8-5838-487a-94d2-790e1bf3464e","arxiv_id":"2607.21047","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"In developing airfoil boundary layers, free-stream turbulence suppresses the adverse-pressure-gradient wake, raises skin friction toward zero-pressure-gradient values, and this influence is amplified as the pressure gradient strengthens.","lead":"Wind-tunnel experiments on an airfoil show that stronger free-stream turbulence partly cancels the expected effects of an adverse pressure gradient: a slimmer wake in the velocity profile, a lower shape factor, and skin friction closer to the zero-pressure-gradient level. The effect grows as the pressure gradient strengthens, indicating that inflow turbulence should be treated as a governing variable when interpreting airfoil boundary layers in realistic conditions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The key quantitative claims rely on uτ and δ from one profile fit; if FST biases that fit, wake suppression and Cf→ZPG could be partly artifacts. Direct wall-shear measurement would settle it.","rationale":"The reader's weakest assumption already identified the same joint uτ/δ estimation issue, and I agree it is the most load-bearing concern. The paper's central claim—that FST systematically modifies APG TBL development and that APG amplifies this influence—rests on quantitative comparisons of H, Cf, wake strength, and λx/δ. All of these are sensitive to the fitted uτ and δ. Direct wall-shear validation exists for ZPG TBLs under FST (Esteban et al., 2017), but not for the combined APG+FST case, so the concern is real rather than hypothetical. However, I do not think it overturns the verdict: the qualitative trends are large, consistent across several independent observables (mean profiles, variance, spectra, scale decomposition), and the paper already acknowledges uncertainty in the low-Tu range. The appropriate disposition remains CONDITIONAL, pending the proposed direct validation and data release. Therefore the reader's verdict should be unchanged.","tokens_in":23801,"tokens_out":3939,"duration_ms":49463,"concrete_test":"Measure wall shear stress directly at x/c=0.625 for α=4° at Tu≈0.2% and Tu≈6% using oil-film interferometry (or a flush-mounted hot-film/Preston tube calibrated in ZPG), and compare with the uτ values from the Rodríguez-López fit. If the direct uτ differs by more than a few percent only for the high-Tu case, the 'Cf toward ZPG' and wake-suppression trends are partly fitting artifacts. As a complementary check, recompute H and λx/δ after perturbing the fitted δ by ±10% to verify that the Tu=6% vs Tu=0.2% differences and the ≈13δ peak survive the perturbation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central quantitative claims—wake suppression, H decreasing to ≈1.4, Cf rising toward the ZPG curve, the λx/δ≈13 energetic scale, and the 8δ scale-split—are all computed from or normalized by uτ and δ returned by a single Rodríguez-López et al. (2015) mean-profile fit (§2.1.3, Appendix A). This is the load-bearing joint because the fit is not independently validated for the combined APG+FST condition. The one direct cross-check attempted, the Vinuesa et al. (2016) δ99 diagnostic, is explicitly unusable under FST (Fig. A.1), and the hot-wire cannot resolve the near-wall gradient. If the fitted profile form absorbs FST-induced outer-region changes by increasing uτ and/or shifting δ, then the observed 'wake suppression' in U+ and the increase of Cf toward ZPG could be, at least partly, artifacts of the fitting parameterization rather than measured flow changes. The same concern propagates into the spectral normalization and the λx=8δ decomposition: because δ is a free parameter, the threshold moves with the fit. This does not disprove the paper's qualitative message—the Tu=6% contrasts are large and internally consistent—but it makes the headline mean-flow claims less secure than the text implies.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Jaroslawski and Scarano report wind-tunnel measurements of a spatially developing turbulent boundary layer on the suction side of a NACA 0015 airfoil at chord Reynolds number ≈2.5×10^5, two angles of attack (2° and 4°), and four freestream turbulence levels (Tu≈0.2%, 1.4%, 2.5%, 6%), at chordwise positions x/c=0.400–0.625. The central claim is that FST is a governing parameter in APG airfoil boundary layers: increasing Tu thickens the boundary layer, lowers the shape factor from ≈1.5 toward ≈1.4, partially suppresses the APG wake, and raises the skin-friction coefficient toward the ZPG correlation. The paper further claims that energetic freestream large-scale motions with λx/δ≈13 penetrate into the boundary layer, amplify inner- and outer-region streamwise variance, and modulate near-wall small scales, with the effect amplified by the adverse pressure gradient. Skin friction and boundary-layer thickness are not measured directly but are obtained together from a single Rodríguez-López et al. (2015) mean-velocity-profile fit; all integral quantities, Reynolds numbers, spectral normalizations, and the λx=8δ scale-split threshold depend on those two fitted parameters.","tokens_in":23932,"tokens_out":7305,"duration_ms":84566,"significance":"If substantiated, this is a valuable contribution: the parametric coverage of Tu and β on a spatially developing airfoil boundary layer is rare, and the combination of mean-flow, variance, spectral, scale-decomposition, and skewness diagnostics provides a rich dataset. The qualitative trends are large, internally consistent, and broadly consistent with earlier ZPG-FST studies and APG-TBL studies. The main weakness is that the headline quantitative results are expressed in inner/outer scaling built on uτ and δ from a single profile fit, with no direct wall-shear measurement and no independent usable check of δ. The qualitative conclusion that FST matters is not in doubt; however, the quantitative strength of the wake suppression, the Cf increase, and the APG-amplification effect is less secure than the text implies. This is fixable with uncertainty quantification and/or an independent skin-friction validation.","major_comments":[{"comment":"The two parameters on which nearly every quantitative statement rests — uτ and δ — are both free parameters of the same Rodríguez-López et al. (2015) profile fit; Appendix A explicitly states that δ is treated as an open fitting parameter together with the friction velocity. Consequently, H (Eq. 11), β (Eq. 6), Cf, Reτ, λx/δ≈13, and the λx=8δ filter threshold are all derived from those fits. Because the wake strength in U+ and Cf are inner-scaled quantities, an FST-dependent overestimate of uτ would by itself produce exactly the reported 'wake suppression' and 'Cf toward ZPG' signatures; a shifted δ would similarly move the spectral peak and the 8δ split. The one independent δ diagnostic attempted, the Vinuesa et al. (2016) method, is explicitly unusable under FST (Fig. A.1). I am not claiming the trends are artifacts, but the load-bearing mean-flow and spectral claims need an independen","section":"§2.1.3, Appendix A"},{"comment":"The headline integral quantities are presented without uncertainty bars or repeatability estimates. The text itself states that for low Tu the shape-factor variation is 'within the uncertainty of the measurements' (§3.1), while the high-Tu reduction from ≈1.5 to ≈1.4 is used to support the central mean-flow claim. Similarly, the Cf increase toward the ZPG curve in Fig. 8 has no error bars, even though Cf is obtained from a fitted uτ. Please provide quantified uncertainties for uτ, δ, H, and Cf for each condition — for example, bootstrap over profile points and fit parameters, or repeated traverses — and state explicitly whether the Tu=1.4% and Tu=2.5% cases are distinguishable from the clean case. Without this, the reader cannot separate genuine FST-driven shifts from fit uncertainty.","section":"§3.1, §3.2, Figs. 6 and 8"},{"comment":"The spectral and scale-decomposition evidence for the proposed mechanism is sensitive to the fitted δ. The wavelength normalization λx/δ and the filter threshold λx=8δ both use the fitted δ; Fig. 5 shows that δ itself increases with Tu. Thus the same physical wavelength is reported as a larger λx/δ, and more energy is assigned to the 'large-scale' band as Tu increases. The claims that the freestream peak is at λx/δ≈13 and that the λx<8δ variance collapses across Tu should therefore be tested for sensitivity to δ. Please show the premultiplied spectra and the filtered variance profiles for δ±10% (or using an independently estimated δ99 where feasible) and report the resulting uncertainty in the peak wavelength and in the large-scale variance amplification. This is a necessary check before the scale-penetration mechanism can be regarded as established.","section":"§3.4, §3.5, Figs. 12–15"},{"comment":"The conclusion that 'the influence of FST is amplified by the adverse pressure gradient' is presented using β as the abscissa, but β is computed from fitted uτ and δ via Eq. (6). Since increasing Tu changes both δ* and τw, part of the observed β shift with Tu is a restatement of the fit outputs rather than an independent measure of the pressure-gradient strength. To make the amplification claim robust, please additionally report the inner- and outer-peak trends against directly measured quantities such as x/c, local edge velocity Ue, or dPe/dx, and note explicitly that β is not an independent control parameter in this experiment. The physical interpretation may be correct, but the current presentation conflates the fitted quantities with the physical pressure-gradient forcing.","section":"§3.1, §3.5, Figs. 10 and 15; Eq. (6)"}],"minor_comments":[{"comment":"The sentence describing the boundary-layer trip is duplicated ('The boundary layer on the suction side of the airfoil was tripped upstream of the measurement region...'). Remove the repetition.","section":"§2.1.1"},{"comment":"Please proofread for typos and grammatical errors: 'usefull', 'despide', 'pronunced', 'thet', 'highligths', 'simular with', 'Figure figure 11', and similar. Several sentences are unfinished or ungrammatical, especially in §3.1.","section":"Throughout"},{"comment":"The figure appears to contain repeated panels with duplicated axes; check the final figure assembly so that each panel is distinct and labeled consistently.","section":"§3.1, Fig. 6"},{"comment":"The notation 'where the double overbar denotes normalization by u+23/2' is obscure. Define the normalization explicitly and state how each term is computed from the filtered signals.","section":"§3.5, Eq. (9)"},{"comment":"The text refers to x/c=0.6235, while all other locations are given as x/c=0.625. Make the reported coordinate consistent.","section":"§3.3"},{"comment":"Table 1 reports Tu=0.3% at x/c=0.400 for α=2°, while §2.1.2 quotes the clean tunnel condition as Tu≈0.2%. Clarify whether these are local values and indicate the precision of the Tu values.","section":"Table 1, §2.1.2"},{"comment":"The caption says the spectral slices are taken 'along the dashed white lines', but the dashed lines are not visible in the rendered figure. Label the row/column positions directly on the figure.","section":"Fig. 13"},{"comment":"There appears to be a mismatch between the panel labels in the figure and the descriptions in the text: the text calls panel (a) the total skewness, but the figure panel (a) appears to show the modulation term. Reorder the panels/caption so they match.","section":"Fig. 16"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about uτ and δ coming from the same profile fit is real and should be addressed before publication. The qualitative message of the paper is likely correct, but the quantitative strength of the wake-suppression, Cf, and APG-amplification claims is currently coupled to the fitting procedure. If direct wall-shear measurements are not available, the authors should at minimum add a thorough sensitivity analysis and soften the quantitative conclusions. The paper is within the scope of this journal and would be a useful contribution after the validation/uncertainty analysis is added."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is the first systematic database of FST×APG on a developing airfoil boundary layer—4 turbulence levels, 2 angles of attack, 4 chord stations, matched chord Reynolds number. The central trends at Tu=6% versus 0.2% are large, internally consistent, and corroborated by independent measurements (mean profiles, variance, spectra, scale decomposition, skewness). The λx/δ≈13 freestream scale is measured outside the boundary layer and independently identified inside it. That combination makes the paper a genuine step beyond the prior FST-in-ZPG and APG-with-clean-inflow literature.\n\nWhere it earns credit: the paper addresses a real gap, and the citation pattern is honest—Blair, Hancock & Bradshaw, Dogan, You & Zaki, Jooss, Monty, Bobke, Tanarro are all there. It also self-reports limitations: §3.1 admits the low-Tu shape-factor variations are within measurement uncertainty, and Appendix A explicitly shows why the Vinuesa δ99 diagnostic fails under FST (Fig. A.1). That transparency is good.\n\nThe soft spot is real but not fatal: uτ and δ come from the same Rodríguez-López profile fit (Appendix A states δ is an open fitting parameter together with uτ). Shape factor H, Clauser β, Cf, Reτ, the λx/δ normalization, and the 8δ scale-split threshold all inherit those two fitted numbers. If the profile form misattributes the FST-modified outer region, then wake suppression and the shift of Cf toward the ZPG curve could be partly fitting artifacts rather than measured changes. The one direct cross-check attempted is explicitly unusable under FST, and the hot-wire (ℓ+≈48–68) cannot resolve the near-wall gradient. So the qualitative message—FST matters and penetrates into the boundary layer, with APG enhancing the effect—is likely right, but the quantitative mean-flow claims are less secure than the text implies. Direct wall-shear measurement, or at least a sensitivity analysis over fit parameter choices, would settle it.\n\nMinor issues: no figure carries error bars despite the admitted low-Tu uncertainty; the \"amplified by APG\" claim rests on figures where β, Reθ, and streamwise position co-vary, with only two matched pairs in Fig. 11; grid geometry is deferred to the authors' 2023b paper; and there is no data-availability statement. These are addressable in revision.\n\nBottom line: this deserves a serious referee. For anyone working on FST or pressure-gradient boundary layers, it is a citable database and a useful statement of the fitting sensitivity. Recommend peer review, conditional on publishing data with uncertainties, tempering the amplification claim, and addressing the uτ/δ fit dependence—ideally with an independent wall-shear cross-check.","headline":"First controlled FST×APG airfoil boundary-layer database; the real risk is that uτ and δ both come from the same profile fit, so the headline mean-flow numbers need an independent cross-check.","tokens_in":24672,"tokens_out":1844,"would_cite":true,"duration_ms":20119,"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":"Free-stream turbulence is a governing parameter in adverse-pressure-gradient airfoil boundary layers: stronger free-stream turbulence reduces shape factor from ~1.5 to ~1.4, suppresses the APG wake, and raises skin friction toward ZPG value","keywords":["free-stream turbulence","adverse pressure gradient","turbulent boundary layer","airfoil boundary layer","shape factor","skin friction","large-scale motions","amplitude modulation"],"falsifier":"Measure skin friction on the same airfoil, at the same Tu and angle of attack, using an independent technique such as oil-film interferometry or a micro force balance, and measure boundary-layer thickness independently with PIV or a dual-probe total-pressure method; if the independently obtained Cf and δ no longer show the reported increase in Cf and decrease in H with Tu, the central claim collapses. Alternatively, a DNS of the same airfoil flow with clean versus FST inflow conditions at matching Rec and β would settle whether the λx/δ ≈ 13 mode is the actual agent.","tokens_in":23481,"feed_emoji":"🌪️","tokens_out":3828,"duration_ms":46950,"temperature":0.7,"pith_summary":"This paper tries to establish that free-stream turbulence (FST) is not a minor disturbance but a governing parameter in a spatially developing adverse-pressure-gradient (APG) turbulent boundary layer over an airfoil. Using wind-tunnel measurements on a NACA 0015 airfoil at fixed chord Reynolds number, the authors show that raising FST from about 0.2% to 6% systematically thickens the boundary layer, lowers the shape factor from roughly 1.5 toward 1.4, partially suppresses the APG-induced wake, and raises skin friction toward zero-pressure-gradient values. The mechanism they identify is the penetration of energetic large-scale free-stream motions, with wavelength around 13 times the boundary-layer thickness, which add energy to the outer layer, leave a near-wall footprint, and modulate near-wall small scales rather than simply superposing. A sympathetic reader would care because real airfoils in wind-turbine wakes, turbomachinery, urban and atmospheric flows experience exactly this combination, and the paper argues that APG makes boundary layers more receptive to external turbulence, so FST must be accounted for when interpreting mean-flow evolution, turbulence statistics, and scale interactions.","feed_headline":"Free-stream turbulence governs APG airfoil boundary layers","feed_subtitle":"At 6% turbulence the wake shrinks, H falls toward 1.4, and Cf climbs toward zero-pressure-gradient values.","key_machinery":"The central object is the set of large-scale free-stream motions with characteristic wavelength λx/δ ≈ 13 that penetrate the turbulent boundary layer and modulate near-wall small scales. The analysis machinery consists of three linked tools: (1) a mean-velocity-profile fit (Rodríguez-López et al.) that simultaneously returns friction velocity uτ and boundary-layer thickness δ, since direct edge detection was unusable under FST; (2) premultiplied spectra and a scale decomposition with a deliberately chosen threshold at λx = 8δ to isolate the free-stream-dominated large scales from the boundary-layer small scales; and (3) a skewness decomposition whose modulation term 3u_L^+ u_S^{+2} quantifie","core_discovery":"The authors claim that in a developing APG turbulent boundary layer over an airfoil, increasing free-stream turbulence systematically alters the mean flow and turbulence structure in ways that counteract canonical APG effects. At the highest FST level (Tu ≈ 6%), the shape factor drops to about 1.4, a value typical of zero-pressure-gradient boundary layers, despite the imposed adverse pressure gradient; the wake region of the mean velocity profile is partially suppressed; and the skin-friction coefficient rises toward the Coles–Fernholz ZPG correlation. Spectral analysis shows that the extra energy comes from large-scale free-stream motions with characteristic wavelength λx/δ ≈ 13, which pene","pith_inferences":["The paper leaves implicit a testable extension: varying the free-stream integral length scale Lu/δ independently of Tu would show whether the reported penetration and near-wall effect scale with Lu/δ and whether the APG amplification systematically increases with that ratio.","A direct consequence the authors do not pursue is that FST, by raising skin friction and lowering shape factor, may delay APG separation on airfoils at otherwise identical pressure distributions; this could be tested by measuring separation location or using surface pressure signatures.","The λx ≈ 13δ scale is close to the grid-generated turbulence's energy-containing scale; in numerical simulations of airfoil flows, inflow boundary conditions would need to reproduce such large-scale modes to capture the physics described here.","The saturation of the large-scale near-wall peak at high Tu and β, which the paper flags as needing further work, suggests there may be a limit to how much additional external energy can penetrate the near-wall region; identifying that limit would sharpen the claimed governing-parameter role."],"forward_implications":["If FST is a governing parameter in APG airfoil boundary layers, then wind-tunnel measurements and numerical simulations of airfoil flows must report and prescribe the free-stream turbulence level and its length scale, not just the pressure gradient and Reynolds number.","Because the APG amplifies FST effects, stronger adverse pressure gradients will make airfoil boundary layers more sensitive to incoming turbulence, meaning separation and skin-friction predictions in turbulent environments cannot rely on clean-inflow data alone.","The large-scale free-stream motions do not merely add energy; they modulate near-wall small scales, so simple additive corrections to turbulence statistics will be insufficient and scale-interaction models are needed.","Matching only Reτ or only β is not enough to characterize the boundary-layer state when FST varies, since profiles with nominally matched parameters differ substantially.","Combined APG and FST can produce features usually associated with higher-Reynolds-number wall turbulence at moderate friction Reynolds numbers, which may affect how flow-control strategies are designed and tested."],"fun_headline_variants":["Freestream turbulence resets APG boundary layer to ZPG-like state","High turbulence shrinks airfoil wake, boosts skin friction","Turbulent inflow counters adverse pressure gradient on airfoil","Large-scale turbulence reshapes airfoil boundary layer","APG boundary layer over airfoil tamed by freestream turbulence"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that the two fitted quantities from one mean-velocity-profile fit, friction velocity uτ and boundary-layer thickness δ, are accurate enough that every headline result—shape factor, Clauser parameter, skin friction, λx/δ ≈ 13, and the λx = 8δ split—faithfully reflects the flow; if the assumed profile form mislocates the edge or misattributes the FST-modified outer region, the reported wake suppression and skin-friction increase could be partly artif","fun_headline_variants_meta":{"raw":{"variants":["Freestream turbulence resets APG boundary layer to ZPG-like state","High turbulence shrinks airfoil wake, boosts skin friction","Turbulent inflow counters adverse pressure gradient on airfoil","Large-scale turbulence reshapes airfoil boundary layer","APG boundary layer over airfoil tamed by freestream turbulence"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000202,"raw_usage":{"total_tokens":1297,"prompt_tokens":897,"completion_tokens":400,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":641,"completion_tokens_details":{"reasoning_tokens":324}},"tokens_in":641,"tokens_out":400,"duration_ms":5899,"temperature":1.0,"reasoning_tokens":324,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T08:37:22.170195+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure skin friction on the same airfoil, at the same Tu and angle of attack, using an independent technique such as oil-film interferometry or a micro force balance, and measure boundary-layer thickness independently with PIV or a dual-probe total-pressure method; if the independently obtained Cf and δ no longer show the reported increase in Cf and decrease in H with Tu, the central claim collapses. Alternatively, a DNS of the same airfoil flow with clean versus FST inflow conditions at matching Rec and β would settle whether the λx/δ ≈ 13 mode is the actual agent.","supporting_citations":[],"review_version":1}