REVIEW 3 major objections 4 minor 36 references
Comparing design and off-design aerodynamic performance of a natural laminar airfoil
T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read 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.
desk verdict 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. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
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.
What would settle it
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.
Extended reading notes
Core claim
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.
Load-bearing premise
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.
Editorial extensions
If this is right
- 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.
Reading between the lines
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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.
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 (3)
- [Table 2 and Section V] 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 IV.B and Table 2] 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.
- [Sections II, III, and V] 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.
minor comments (4)
- [Figure 9 caption] 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 IV.A] The text contains a typo: 'Kelvin-Helmotz' should be 'Kelvin-Helmholtz'.
- [Section III] 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 IV.B] 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.
Circularity Check
No circular construction: the simulations are checked against independent Fujino flight-test Cp data, and the reported lift/drag coefficients are integrals of the simulated flow, not fitted inputs.
full rationale
The derivation chain is self-contained in the relevant sense. The paper solves the 2D compressible Navier-Stokes equations with high-accuracy compact schemes, validates the pressure distribution against independent flight-test data from Fujino et al. (Fig. 2), and then computes time-averaged Cp, Cf, Cl, and Cd by integrating the simulated surface pressure and shear. No parameter is fitted to the target quantities: the lift and drag coefficients are derived from the flow fields, and the only external comparisons are the experimental Cp distributions, which are independent of the authors' prior work. The paper does cite several prior papers by the same group for the numerical framework, grid strategy, and related SHM1 test-case context (e.g., refs. [5, 8, 21, 22, 24, 26, 27]), but the present validation is not merely an appeal to those citations: Fig. 2 directly overlays simulated Cp against Fujino et al.'s flight-test data, and the text states 'the results show good agreement between the simulated and experimental Cp.' The numerical methods are also checked against the external NACA0012 wind-tunnel benchmark of Harris [23]. Thus the central claims do not reduce to their inputs by construction, and no fitted parameter is renamed as a prediction. The internal arithmetic inconsistency in Table 2 (e.g., Case-1 Cl/Cd = 0.25089 / 4.66969e-3 = 53.73, not the tabulated 5.37273) is a correctness and self-consistency defect, but it is not a circular-reasoning defect, so it does not change this circularity verdict.
Assumptions & free parameters
assumptions (4)
- domain assumption Two-dimensional compressible Navier-Stokes equations with implicit large eddy simulation are adequate at Reynolds numbers up to 13.6e6; no explicit turbulence or transition model is applied.
- domain assumption Viscosity-temperature law and bulk viscosity value needed to close the stress tensor are unspecified.
- domain assumption Outer boundary at 16c with characteristic non-reflective conditions is sufficiently far to avoid spurious reflections.
- domain assumption Time averaging from t=20 to 100 at Delta t=2.5e-6 yields statistically converged mean Cp, Cf, Cl, and Cd.
Cite this review
Pith. "Pith review of Comparing design and off-design aerodynamic performance of a natural laminar airfoil." pith.science (2026). https://pith.science/paper/HDK5STU2
@misc{pith2026241112266,
author = {Pith},
title = {Pith review of: Comparing design and off-design aerodynamic performance of a natural laminar airfoil},
year = {2026},
howpublished = {\url{https://pith.science/paper/HDK5STU2}},
note = {Machine review of arXiv:2411.12266}
}
read the original abstract
Natural laminar flow airfoils are essential technologies designed to reduce drag and significantly enhance aerodynamic performance. A notable example is the SHM1 airfoil, created to meet the requirements of the small-business Honda jet. This airfoil has undergone extensive testing across various operational conditions, including low-speed wind tunnel tests and flight tests across a range of Reynolds numbers and free-stream Mach numbers, as detailed in "Natural-laminar-flow airfoil development for a lightweight business jet" by Fujino et al., J. Aircraft, 40(4), 2003. Additionally, investigations into drag-divergence behavior have been conducted using a transonic wind tunnel, with subsequent studies focusing on transonic shock boundary layer interactions through both experimental and numerical approaches. This study employs a series of numerical simulations to analyze the flow physics and aerodynamic performance across different free-stream Mach numbers in the subsonic and transonic regimes. This is achieved by examining computed instantaneous numerical Schlieren for various design conditions (such as low speed, climb, and cruise) and off-design scenarios (including transonic shock emergence, drag-divergence, and shock-induced separation). The dominant time scales, the time-averaged load distributions and boundary layer parameters are compared to provide a comprehensive overview of the SHM1's aerodynamics, establishing benchmark results for optimization of various flow separation and shock control techniques.
Reference graph
Works this paper leans on
-
[1]
Coder, J. G., and Maughmer, M. D., “Computational fluid dynamics compatible transition modeling using an amplification factor transport equation,”AIAA journal, Vol. 52, No. 11, 2014, pp. 2506–2512
work page 2014
-
[2]
Review of hybrid laminar flow control systems,
Krishnan, K., Bertram, O., and Seibel, O., “Review of hybrid laminar flow control systems,”Progress in Aerospace Sciences, Vol. 93, 2017, pp. 24–52
work page 2017
-
[3]
Investigation of Hybrid Laminar Flow Control Capabilities from the Flight Envelope Perspective,
Karpuk, S., and Mosca, V., “Investigation of Hybrid Laminar Flow Control Capabilities from the Flight Envelope Perspective,” Journal of Aircraft, 2024, pp. 1–26
work page 2024
-
[4]
Adjoint-based aerodynamic shape optimization including transition to turbulence effects,
Halila, G. L., Martins, J. R., and Fidkowski, K. J., “Adjoint-based aerodynamic shape optimization including transition to turbulence effects,”Aerospace Science and Technology, Vol. 107, 2020, p. 106243
work page 2020
-
[5]
Chakraborty, A., Roy, A. G., Sundaram, P., Sengupta, A., and Sengupta, T. K., “Controlling transonic shock–boundary layer interactions over a natural laminar flow airfoil by vortical and thermal excitation,”Physics of Fluids, Vol. 34, No. 8, 2022
work page 2022
-
[6]
Natural-laminar-flow airfoil development for a lightweight business jet,
Fujino, M., Yoshizaki, Y., and Kawamura, Y., “Natural-laminar-flow airfoil development for a lightweight business jet,”Journal of Aircraft, Vol. 40, No. 4, 2003, pp. 609–615
work page 2003
-
[7]
Effects of boundary layer transition on the aerodynamic analysis of high-lift systems,
Halila, G. L. O., Antunes, A. P., Da Silva, R. G., and Azevedo, J. L. F., “Effects of boundary layer transition on the aerodynamic analysis of high-lift systems,”Aerospace Science and Technology, Vol. 90, 2019, pp. 233–245
work page 2019
-
[8]
Thermal control of transonic shock-boundary layer interaction over a natural laminar flow airfoil,
Sengupta, T. K., Roy, A. G., Chakraborty, A., Sengupta, A., and Sundaram, P., “Thermal control of transonic shock-boundary layer interaction over a natural laminar flow airfoil,”Physics of Fluids, Vol. 33, No. 12, 2021
work page 2021
Show all 36 references
-
[9]
Building a practical natural laminar flow design capability,
Campbell, R. L., and Lynde, M. N., “Building a practical natural laminar flow design capability,”35th AIAA Applied Aerodynamics Conference, 2017, p. 3059
2017
-
[10]
D., and Cook, L
Cole, J. D., and Cook, L. P.,Transonic aerodynamics, Elsevier, 2012
2012
-
[11]
Reynolds-andMach-numbereffectsincanonicalshock–turbulenceinteraction,
Larsson,J.,Bermejo-Moreno,I.,andLele,S.K.,“Reynolds-andMach-numbereffectsincanonicalshock–turbulenceinteraction,” Journal of Fluid Mechanics, Vol. 717, 2013, pp. 293–321
2013
-
[12]
Experimentalstudyofanormalshock/homogeneousturbulenceinteraction,
Barre,S.,Alem,D.,andBonnet,J.,“Experimentalstudyofanormalshock/homogeneousturbulenceinteraction,” AIAAjournal, Vol. 34, No. 5, 1996, pp. 968–974
1996
-
[13]
An experimental study of transonic shock/boundary layer interactions subject to downstream pressure perturbations,
Bruce, P., and Babinsky, H., “An experimental study of transonic shock/boundary layer interactions subject to downstream pressure perturbations,”Aerospace Science and Technology, Vol. 14, No. 2, 2010, pp. 134–142. 22
2010
-
[14]
Numerical and experimental study of nominal 2-D shock-wave/turbulent boundary layer interactions,
Touré, P., and Schuelein, E., “Numerical and experimental study of nominal 2-D shock-wave/turbulent boundary layer interactions,”2018 Fluid Dynamics Conference, 2018, p. 3395
2018
-
[15]
Numerical analysis of laminar and turbulent shock-wave boundary layer interactions,
Gross, A., and Lee, S., “Numerical analysis of laminar and turbulent shock-wave boundary layer interactions,”2018 Fluid Dynamics Conference, 2018, p. 4033
2018
-
[16]
Numerical investigation of transitional shock-wave/boundary-layer interaction in supersonic regime,
Quadros, R., and Bernardini, M., “Numerical investigation of transitional shock-wave/boundary-layer interaction in supersonic regime,”AIAA Journal, Vol. 56, No. 7, 2018, pp. 2712–2724
2018
-
[17]
Shock/boundary-layer interactions: Possible sources of unsteadiness,
Dussauge, J.-P., and Piponniau, S., “Shock/boundary-layer interactions: Possible sources of unsteadiness,”Journal of Fluids and Structures, Vol. 24, No. 8, 2008, pp. 1166–1175
2008
-
[18]
A review of recent developments in the understanding of transonic shock buffet,
Giannelis, N., Vio, G., and Levinski, O., “A review of recent developments in the understanding of transonic shock buffet,” Progress in Aerospace Sciences, Vol. 92, 2017, pp. 39–84
2017
-
[19]
Self-sustained shock oscillations on airfoils at transonic speeds,
Lee, B., “Self-sustained shock oscillations on airfoils at transonic speeds,”Progress in Aerospace Sciences, Vol. 37, No. 2, 2001, pp. 147–196
2001
-
[20]
Global spectral analysis: Review of numerical methods,
Sagaut, P., Suman, V. K., Sundaram, P., Rajpoot, M. K., Bhumkar, Y. G., Sengupta, S., Sengupta, A., and Sengupta, T. K., “Global spectral analysis: Review of numerical methods,”Computers & Fluids, 2023, p. 105915
2023
-
[21]
Non-overlapping High-accuracy Parallel Closure for Compact Schemes: Application in Multiphysics and Complex Geometry,
Sundaram, P., Sengupta, A., Suman, V. K., and Sengupta, T. K., “Non-overlapping High-accuracy Parallel Closure for Compact Schemes: Application in Multiphysics and Complex Geometry,”ACM Transactions on Parallel Computing, Vol. 10, No. 1, 2023, pp. 1–28
2023
-
[22]
Comparativestudyoftransonicshock–boundary layer interactions due to surface heating and cooling on an airfoil,
Sengupta,T.K.,Chakraborty,A.,Roy,A.G.,Sengupta,A.,andSundaram,P.,“Comparativestudyoftransonicshock–boundary layer interactions due to surface heating and cooling on an airfoil,”Physics of Fluids, Vol. 34, No. 4, 2022
2022
-
[23]
Two-dimensional aerodynamic characteristics of the NACA 0012 airfoil in the Langley 8 foot transonic pressure tunnel,
Harris, C. D., “Two-dimensional aerodynamic characteristics of the NACA 0012 airfoil in the Langley 8 foot transonic pressure tunnel,” Tech. rep., 1981
1981
-
[24]
Thermal optimization of shock-induced separation in a natural laminar airfoil operating at off-design conditions,
Sengupta, A., and Shandilya, N., “Thermal optimization of shock-induced separation in a natural laminar airfoil operating at off-design conditions,”Physics of Fluids, Vol. 36, No. 4, 2024
2024
-
[25]
Implicit solution methods in computational fluid dynamics,
Pulliam, T. H., “Implicit solution methods in computational fluid dynamics,”Applied numerical mathematics, Vol. 2, No. 6, 1986, pp. 441–474
1986
-
[26]
A non-overlapping high accuracy parallel subdomain closure for compact scheme: Onset of Rayleigh-Taylor instability by ultrasonic waves,
Sundaram, P., Sengupta, A., and Sengupta, T. K., “A non-overlapping high accuracy parallel subdomain closure for compact scheme: Onset of Rayleigh-Taylor instability by ultrasonic waves,”Journal of Computational Physics, Vol. 470, 2022, p. 111593
2022
-
[27]
Direct numerical simulation of 2D transonic flows around airfoils,
Sengupta, T. K., Bhole, A., and Sreejith, N., “Direct numerical simulation of 2D transonic flows around airfoils,”Computers & Fluids, Vol. 88, 2013, pp. 19–37. 23
2013
-
[28]
Numerical study of wind-tunnel walls effects on transonic airfoil flow,
Garbaruk, A., Shur, M., Strelets, M., and Spalart, P. R., “Numerical study of wind-tunnel walls effects on transonic airfoil flow,” AIAA journal, Vol. 41, No. 6, 2003, pp. 1046–1054
2003
-
[29]
Subsonic natural-laminar-flow airfoils,
Somers, D. M., “Subsonic natural-laminar-flow airfoils,”Natural laminar flow and laminar flow control, Springer, 1992, pp. 143–176
1992
-
[30]
Shock-wave boundary layer interactions,
Délery, J., Marvin, J. G., and Reshotko, E., “Shock-wave boundary layer interactions,” AGARD, 1986
1986
-
[31]
Visualization,extractionandquantification of discontinuities in compressible flows,
Samtaney,R.,Morris,R.,Cheeseman,P.,Sunelyansky,V.,Maluf,D.,andWolf,D.,“Visualization,extractionandquantification of discontinuities in compressible flows,”International conference on computer vision and pattern recognition, 2000
2000
-
[32]
L., and Carpenter, P
Houghton, E. L., and Carpenter, P. W.,Aerodynamics for engineering students, Elsevier, 2003
2003
-
[33]
Investigations of the transonic flow around oscillating airfoils,
Tijdeman, H., “Investigations of the transonic flow around oscillating airfoils,”NLR-TR 77090 U, 1977
1977
-
[34]
New results for the normal shock in inviscid flow at a curved surface,
Zierep, J., “New results for the normal shock in inviscid flow at a curved surface,”ZAMM-Journal of Applied Mathematics and Mechanics/Zeitschrift für Angewandte Mathematik und Mechanik: Applied Mathematics and Mechanics, Vol. 83, No. 9, 2003, pp. 603–610
2003
-
[35]
Separation-induced transition on a T106A blade under low and elevated free stream turbulence,
Sengupta, A., Gupta, N., and Ubald, B., “Separation-induced transition on a T106A blade under low and elevated free stream turbulence,”Physics of Fluids, Vol. 36, No. 2, 2024
2024
-
[36]
Natural-laminar-flow airfoil for general-aviation applications,
Selig, M. S., Maughmer, M. D., and Somers, D. M., “Natural-laminar-flow airfoil for general-aviation applications,”Journal of aircraft, Vol. 32, No. 4, 1995, pp. 710–715. 24
1995
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.