{"id":"b7246f65-9a5a-49fb-91ab-1cc8d2a3c6ba","arxiv_id":"2411.12266","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Numerical simulations map the SHM1 airfoil from low-speed to shock-induced separation, but the claimed quantitative drag benchmark is not yet validated.","lead":"This paper simulates the SHM1 natural laminar airfoil at six flight conditions, from low-speed climb to transonic shock-induced separation, and compares computed pressure, skin friction, and lift and drag. It aims to provide a benchmark for optimizing shock and separation control on this business-jet airfoil.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Table 2's L/D values are internally inconsistent (Cases 1 and 6 off by ~10x) and the reported minimum drag at shock-induced separation is physically implausible, undermining the benchmark claim as printed.","rationale":"The reader's weakest assumption focuses on the credibility of 2D ILES at high Reynolds numbers without grid-convergence evidence. That is a legitimate concern about whether the simulation resolves the relevant physics. However, the more decisive and concrete problem is that the paper's central quantitative dataset, Table 2, contains self-contradictory L/D ratios and a physically implausible drag ranking. This is not a speculation about resolution; it is an internal check that can be settled immediately by recomputing the printed numbers. The corrected Case-6 L/D (102.04) would overturn the stated conclusion that the cruise condition is the most efficient. Because the central claim of establishing benchmark results depends on these numbers, and because the numbers as printed contradict both arithmetic and physical expectation, the current manuscript should not be accepted. A future version that corrects Table 2, verifies the drag integration, and addresses the resolution concern could be reconsidered, but the present central quantitative argument is not reliable.","tokens_in":13911,"tokens_out":5088,"duration_ms":58851,"concrete_test":"Recompute every Cl/Cd ratio in Table 2 from the printed Cl and Cd values. For Case-1 and Case-6, if the recomputed ratios (53.73 and 102.04) do not match the table entries, the table is internally inconsistent. Then recompute Cd for Cases 3 and 6 by direct integration of the time-averaged pressure and skin-friction distributions used for Figures 8 and 9. If Case-6's recomputed Cd is not lower than Case-3, the reported drag ranking and the associated benchmark conclusions are invalid.","verdict_should_be":"REJECT","load_bearing_attack":"To support the central claim of reliable benchmark aerodynamic performance, Table 2 must present mutually consistent Cl, Cd, and Cl/Cd values. It does not. For Case-1, Cl/Cd = 0.25089 / 4.66969e-3 = 53.73, but the table lists 5.37273. For Case-6, Cl/Cd = 0.26214 / 2.56901e-3 = 102.04, but the table lists 10.20380. These are not round-off effects; they are order-of-magnitude internal contradictions in the quantitative core of the paper. The stated conclusion that Case-6 has the least drag is also physically suspect: a shock-induced separation case at M=0.78 is reported with Cd = 2.57e-3, lower than the M=0.62 cruise value of 3.06e-3, despite strong normal shocks and a large separated region on the suction surface, which should substantially increase drag. If the corrected Case-6 Cl/Cd is 102.04, the paper's conclusion that 'the cruise condition emerged as the most aerodynamically efficient' is reversed. This is more direct than the grid-resolution concern: even granting the ILES methodology, the reported performance data contradict themselves and the paper's own conclusions.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents two-dimensional implicit large eddy simulations of the SHM1 natural laminar flow airfoil at six conditions spanning low-speed, climb, cruise, transonic shock appearance, drag divergence, and shock-induced separation. The paper compares instantaneous numerical Schlieren, vorticity time series and spectra, time-averaged pressure and skin-friction distributions, and reports integrated lift and drag coefficients in Table 2. The authors claim that the results provide benchmark aerodynamic performance data for the SHM1 airfoil and its future optimization. Validation is limited to a visual comparison of time-averaged pressure coefficient with Fujino et al.'s flight test data at Mach 0.62 and 0.72.","tokens_in":14161,"tokens_out":3873,"duration_ms":44232,"significance":"The qualitative flow-physics comparisons in the manuscript, particularly the identification of Kutta waves, lambda shocks, and shock-boundary-layer interaction regimes, are useful and align with prior experimental and computational observations. The use of high-order dispersion-relation-preserving compact schemes and validation of pressure distributions against independent flight-test data at two Mach numbers are genuine strengths, and the absence of fitted free parameters in the reported simulations is a positive feature. However, the central quantitative claim of the paper, namely that Table 2 provides reliable benchmark lift and drag coefficients, is not supported as printed. The integrated performance table contains internal arithmetic inconsistencies that change the paper's conclusions, and no grid-convergence study or independent validation of lift/drag is provided. The manuscript's significance would be substantially strengthened by correcting the table, reconciling the text with the corrected values, and adding a resolution/verification study.","major_comments":[{"comment":"Table 2 is internally inconsistent. For Case-1, Cl/Cd = 0.25089 / 4.66969e-3 = 53.73, but the table lists 5.37273; for Case-6, Cl/Cd = 0.26214 / 2.56901e-3 = 102.04, but the table lists 10.20380. These are factor-of-ten errors, not rounding effects. Because the paper's benchmark claim and its conclusion that 'the cruise condition emerged as the most aerodynamically efficient' rely directly on these values, the table must be corrected and the discussion in Sections IV.B and V revisited. If the corrected Case-6 value is approximately 102, then Case-6, not Case-3, has the highest Cl/Cd in the table, reversing the stated conclusion.","section":"Table 2 and Section V"},{"comment":"The text states that the shock-induced-separation case (Case-6) 'shows the largest streamwise extent of the separated region with largest skin friction drag,' but Table 2 reports the lowest total drag coefficient (2.56901e-3) for this case, even lower than the cruise value (3.05563e-3). This is physically implausible for a case with strong normal shocks and a large separated region, and it directly contradicts the paper's own statement. The authors should either correct the drag data or explain why the integrated drag is nonetheless smallest; otherwise the quantitative performance discussion is unreliable.","section":"Section IV.B and Table 2"},{"comment":"No grid-convergence study or resolution assessment is reported for the 1251x401 O-grid, and the integrated lift and drag coefficients are not validated against any experimental or high-fidelity reference. The runs are described as implicit large eddy simulations, yet the governing equations are two-dimensional and no subgrid-scale model or transition/turbulence model is described. For Reynolds numbers up to 13.6e6, the claim that these simulations resolve the relevant boundary-layer and shock dynamics well enough for quantitative Cl and Cd values in Table 2 requires at least a grid-refinement study, a time-step sensitivity check, and preferably a comparison of integrated forces with experimental data. Without this, the benchmark results for aerodynamic performance are not supported.","section":"Sections II, III, and V"}],"minor_comments":[{"comment":"The caption for Figure 9 lists frames (a) through (f), but the figure displays only two panels, labeled (a) and (b), each containing multiple curves. The caption should be revised to describe the two-panel layout accurately.","section":"Figure 9 caption"},{"comment":"The text contains a typo: 'Kelvin-Helmotz' should be 'Kelvin-Helmholtz'.","section":"Section IV.A"},{"comment":"The sentence 'This numerical framework has been validated in a prior work [21]' is misleading because reference [21] is primarily a parallel-computing paper, not a validation study of the SHM1 airfoil. The actual validation comparison is presented in Figure 2 and should be described as such.","section":"Section III"},{"comment":"The statement that the climb condition (Case-2) shows elevated drag 'due to the elevated thrust requirement' is not a physical explanation; thrust is not an input to the isolated-airfoil simulation, and the drag increase should be traced to the modified pressure and skin-friction distributions that are shown in Figures 8 and 9.","section":"Section IV.B"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's core quantitative deliverable is Table 2, and the factor-of-ten internal inconsistencies make the paper unsuitable for publication in its current form. I would ask the authors to recompute the table, correct the associated conclusions, and add a grid-convergence study or at least an explicit resolution assessment before resubmission. The paper also leans heavily on the authors' own prior work as authority; independent validation of the integrated loads would materially improve the persuasiveness of the revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a visually rich, well-written CFD study of the SHM1 airfoil, but the quantitative core—the Cl, Cd, Cl/Cd table—is not usable as printed. The lift-to-drag ratios for Cases 1 and 6 are off by a factor of ten relative to the listed Cl and Cd values, and the conclusion that Case-6 has the lowest drag contradicts both the text's own statement about largest skin friction drag and basic physics for a shock-separated case at M=0.78. The paper needs a straightforward revision; the errors are glaring but fixable.\n\nWhat's genuinely new: the systematic comparison of six operating conditions for this airfoil, from low-speed to shock-induced separation, in one place. The Cp validation against Fujino et al.'s flight data at M=0.62 and 0.72 is convincing and gives some confidence in the qualitative flow pictures. The numerical Schlieren and the discussion of Kutta waves, lambda shock formation, and delayed separation with Mach number are plausible and well presented. The spectral analysis of vorticity probes is a nice touch.\n\nThe soft spots are the load-bearing ones. First, Table 2: Case-1 lists Cl/Cd=5.37 but 0.25089/0.00466969 is 53.7; Case-6 lists 10.20 but 0.26214/0.00256901 is 102.0. The other four rows are consistent. So it looks like two decimal slips, but the paper's claim that the cruise condition has the best Cl/Cd and that Case-6 has the worst is exactly backwards if those corrected values hold. Second, the text says Case-6 shows the largest skin friction drag, yet Table 2 gives it the lowest total drag. A normal shock plus a large separated region should increase both friction and pressure drag; a total Cd below the cruise value is implausible. Third, the quantitative benchmark rests on a 2D ILES with no grid-convergence study, no error bars, and no validation of Cd or Cl. The Cp agreement is nice, but Cp is a weak constraint on drag. Also, the validation is largely inherited from earlier papers by the same group, so independent evidence is thin.\n\nFor whom is this paper? Readers working on natural laminar flow airfoils or transonic shock-boundary layer interactions will find the flow insights useful. As a benchmark for optimization, it is not ready. It deserves a serious referee, but the referee should send it back for major revision.","headline":"Table 2's L/D values are off by a factor of ten for two cases, and that contradiction flips the paper's efficiency conclusion; the qualitative flow physics is solid but the quantitative benchmark is not.","tokens_in":14679,"tokens_out":3527,"would_cite":false,"duration_ms":35299,"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":"Six implicit large eddy simulations of the SHM1 natural-laminar-flow airfoil map its aerodynamic performance from low-speed cruise to shock-induced separation, and the results are offered as a benchmark for separation and shock control.","keywords":["natural laminar flow airfoil","SHM1","transonic shock boundary layer interaction","implicit large eddy simulation","numerical Schlieren","lift and drag coefficients","Kutta waves","drag divergence"],"falsifier":"Run any one of the six cases on a systematically refined grid, or with an explicit subgrid-scale model, and check whether the time-averaged Cl and Cd in Table 2 change by more than a few percent; alternatively, compare the M=0.62 and M=0.72 integrated lift and drag against three-dimensional flight-test measurements.","tokens_in":13711,"feed_emoji":"✈️","tokens_out":4714,"duration_ms":43505,"temperature":0.7,"pith_summary":"The paper uses six two-dimensional implicit large eddy simulations to compare the SHM1 natural-laminar-flow airfoil at low speed, climb, cruise, and three transonic off-design points. It reports time-averaged pressure, skin friction, vorticity spectra, and integrated lift and drag, and it argues these results, validated against flight-test pressure data at two Mach numbers, provide a benchmark for flow-separation and shock-control optimization. The central payoff is a consistent picture of how the flow breaks down as Mach number rises: a shock-free cruise, a first normal shock, drag divergence, and finally a Lambda shock with separation.","feed_headline":"Simulations map SHM1 airfoil drag from cruise to shock separation","feed_subtitle":"Six implicit large-eddy runs match flight-test pressure and locate where lift and drag break down.","key_machinery":"The workhorse is a set of six test cases run with dispersion-relation-preserving compact schemes solving the two-dimensional compressible Navier-Stokes equations on a 1251-by-401 O-grid, with implicit large eddy simulation providing the turbulence resolution. Time-averaging from t=20 to 100 yields the pressure and skin-friction distributions, while numerical Schlieren (density-gradient magnitude) reveals the shock structures. The mechanism carries the argument by connecting flow visualization, spectral analysis of vorticity, and integrated force coefficients to the validated Cp curves.","core_discovery":"The central discovery is a self-consistent set of computed aerodynamic states for the SHM1 airfoil spanning design and off-design operation. At M=0.62 the flow is shock-free and produces the highest lift-to-drag ratio; at M=0.72 a curved near-normal shock appears at about x/c=0.45; at M=0.73 the shock strengthens and moves aft; and at M=0.78 a Lambda shock forms on the suction surface and separates the boundary layer, causing a drop in lift despite the lowest drag of the transonic cases. The paper also shows that upstream-propagating Kutta waves interact nonlinearly with the boundary layer and shocks, making the vorticity spectra increasingly multi-periodic as Mach number rises.","pith_inferences":["If the two-dimensional ILES resolution is adequate, the reported table of Cl, Cd, and Cl/Cd could serve as a quick screening tool for flow-control strategies, but the quantitative values should be checked against a finer grid before being used for design decisions.","The plateau-and-spike Cp signature the paper identifies on the suction surface is probably a general marker for natural-laminar-flow airfoils entering transonic conditions; testing other NLF geometries would show whether the trend extends.","The delayed-separation trend with increasing Mach number hints that the Lambda-shock regime may be controllable with modest thermal or vortical excitation, a possibility the paper does not test but its benchmark setup would support."],"forward_implications":["At the cruise condition M=0.62, the SHM1 achieves a computed Cl/Cd near 99, the best of the six cases studied.","The first transonic shock appears near M=0.72 at x/c approximately 0.45, and moves aft to 0.5 and 0.6 as Mach number rises to 0.73 and 0.78.","The drag-divergence case, M=0.73, generates the highest lift but also the highest drag, while the M=0.78 shock-induced separation case has the lowest drag but a large lift penalty.","Vorticity spectra show a single dominant shedding frequency at low speed, and become chaotic and multi-periodic once transonic shocks and Kutta-wave interactions appear."],"supporting_citations":[{"why":"Supplies the flight-test pressure data used for validating the simulated Cp distributions.","marker":"[6]"},{"why":"Provides prior transonic shock-boundary-layer interaction simulations on the same airfoil that this study extends.","marker":"[8]"},{"why":"Describes the dispersion-relation-preserving compact schemes used for spatial discretization.","marker":"[20]"},{"why":"Provides the non-overlapping parallel closure strategy that keeps the numerical solution accurate in parallel.","marker":"[21]"},{"why":"Offers the benchmark wind-tunnel data for NACA0012 used to validate the numerical method.","marker":"[23]"},{"why":"Previous work on shock-induced separation in a natural laminar airfoil that informs the off-design cases.","marker":"[24]"},{"why":"Validates the same numerical approach for two-dimensional transonic airfoil flows against the NACA0012 benchmark.","marker":"[27]"}],"fun_headline_variants":["SHM1 airfoil aerodynamics mapped from cruise to shock separation","Six LES runs reveal SHM1 drag breakdown from cruise to shock","SHM1 simulations show lift loss at shock-induced separation","Simulations trace SHM1 airfoil from shock-free cruise to Lambda shock"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative lift and drag numbers rest on the assumption that a single two-dimensional grid, with no explicit turbulence model and no grid-convergence study, resolves the unsteady boundary layer and shock motion well enough to make the time-averaged forces trustworthy.","fun_headline_variants_meta":{"raw":{"variants":["SHM1 airfoil aerodynamics mapped from cruise to shock separation","Six LES runs reveal SHM1 drag breakdown from cruise to shock","SHM1 simulations show lift loss at shock-induced separation","Simulations trace SHM1 airfoil from shock-free cruise to Lambda shock"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000258,"raw_usage":{"total_tokens":1584,"prompt_tokens":946,"completion_tokens":638,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":562,"completion_tokens_details":{"reasoning_tokens":561}},"tokens_in":562,"tokens_out":638,"duration_ms":6237,"temperature":1.0,"reasoning_tokens":561,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T17:43:59.721492+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run any one of the six cases on a systematically refined grid, or with an explicit subgrid-scale model, and check whether the time-averaged Cl and Cd in Table 2 change by more than a few percent; alternatively, compare the M=0.62 and M=0.72 integrated lift and drag against three-dimensional flight-test measurements.","supporting_citations":[{"cited_title":"Natural-laminar-flow airfoil development for a lightweight business jet,","cited_arxiv_id":null,"evidence_quote":"Supplies the flight-test pressure data used for validating the simulated Cp distributions."},{"cited_title":"Thermal control of transonic shock-boundary layer interaction over a natural laminar flow airfoil,","cited_arxiv_id":null,"evidence_quote":"Provides prior transonic shock-boundary-layer interaction simulations on the same airfoil that this study extends."},{"cited_title":"Global spectral analysis: Review of numerical methods,","cited_arxiv_id":null,"evidence_quote":"Describes the dispersion-relation-preserving compact schemes used for spatial discretization."},{"cited_title":"Non-overlapping High-accuracy Parallel Closure for Compact Schemes: Application in Multiphysics and Complex Geometry,","cited_arxiv_id":null,"evidence_quote":"Provides the non-overlapping parallel closure strategy that keeps the numerical solution accurate in parallel."},{"cited_title":"Two-dimensional aerodynamic characteristics of the NACA 0012 airfoil in the Langley 8 foot transonic pressure tunnel,","cited_arxiv_id":null,"evidence_quote":"Offers the benchmark wind-tunnel data for NACA0012 used to validate the numerical method."},{"cited_title":"Thermal optimization of shock-induced separation in a natural laminar airfoil operating at off-design conditions,","cited_arxiv_id":null,"evidence_quote":"Previous work on shock-induced separation in a natural laminar airfoil that informs the off-design cases."},{"cited_title":"Direct numerical simulation of 2D transonic flows around airfoils,","cited_arxiv_id":null,"evidence_quote":"Validates the same numerical approach for two-dimensional transonic airfoil flows against the NACA0012 benchmark."}],"review_version":1}