Pith. sign in

REVIEW 2 major objections 1 minor 49 references

A pressure-gradient wall model combined with the fifth-order compact gas-kinetic scheme reproduces separated compressible flows on near-wall meshes twenty times coarser than wall-resolved requirements.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · grok-4.3

2026-06-30 03:54 UTC pith:5RMRUV7D

load-bearing objection The paper couples a pressure-gradient non-equilibrium wall model to CGKS-5th and demonstrates usable accuracy on 20x-coarsened near-wall grids for two separated-flow cases, but the model itself is validated only indirectly through overall flow agreement. the 2 major comments →

arxiv 2606.30061 v1 pith:5RMRUV7D submitted 2026-06-29 physics.flu-dyn cs.NAmath.NA

Efficient Wall-Modeled High-Order Compact Gas-Kinetic Scheme for Compressible Turbulent Flows

classification physics.flu-dyn cs.NAmath.NA
keywords wall-modeled simulationgas-kinetic schemecompressible turbulent flowsseparated flowshigh-order compact schemenear-wall modelingskin-friction prediction
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a wall-modeled version of the CGKS-5th scheme to reduce the high computational cost of resolving near-wall turbulence in compressible flows. It pairs the high-order outer solver with a non-equilibrium wall model that uses pressure gradients to reconstruct wall stress from coarse data. Tests on flow past a cylinder and the RAE 2822 airfoil show that the method captures separation and improves skin-friction predictions compared to the wall-model-free version. The coupling adds negligible overhead. This approach makes scale-resolving simulations of high-Reynolds-number wall-bounded flows more practical.

Core claim

The wall-modeled CGKS-5th framework, which uses the outer solver to provide data at an exchange location to a pressure-gradient-based non-equilibrium wall model that reconstructs viscous wall stress while the scheme supplies inviscid flux, reproduces separated flow structures and improves near-wall predictions such as skin-friction coefficient on meshes far coarser than those needed for wall-resolved simulations.

What carries the argument

The pressure-gradient-based non-equilibrium wall model that retains a pressure-gradient source term and a corrected damping function, coupled to CGKS-5th at the exchange location to form the wall momentum flux.

Load-bearing premise

The pressure-gradient-based non-equilibrium wall model accurately reconstructs the under-resolved viscous wall stress when given data from the outer solver, even in adverse pressure gradient and separated regions.

What would settle it

A direct comparison where the wall-modeled scheme's predicted skin-friction coefficient on the coarse mesh deviates substantially from experimental or high-resolution reference values in the RAE 2822 or cylinder cases.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Reproduces separated flow structures in bluff-body and shock-induced separation cases.
  • Markedly improves skin-friction coefficient predictions over wall-model-free CGKS-5th.
  • Maintains accuracy with twentyfold coarsening in wall-normal direction for the airfoil case.
  • Adds less than 1% runtime overhead in multi-GPU implementation.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the model generalizes, it could enable routine high-fidelity simulations of aircraft or turbomachinery at flight Reynolds numbers.
  • Further tests on other geometries with stronger separation might reveal limits of the pressure-gradient correction.
  • The lightweight coupling suggests easy integration into other high-order compressible solvers.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 1 minor

Summary. The paper develops a wall-modeled extension of the fifth-order compact gas-kinetic scheme (CGKS-5th) for compressible turbulent flows. A pressure-gradient-based non-equilibrium wall model supplies the viscous wall stress on coarse near-wall meshes while the outer CGKS-5th solver supplies the inviscid flux; the two are combined at the wall. The framework is tested on bluff-body separation past a circular cylinder and shock-induced separation on the RAE 2822 airfoil, using meshes coarsened by a factor of approximately 20 in the wall-normal direction. The authors report that the wall-modeled scheme reproduces separated flow structures and markedly improves skin-friction predictions relative to the wall-model-free CGKS-5th.

Significance. If the central claim holds, the work would offer a practical route to reducing the near-wall resolution cost of high-order scale-resolving simulations of compressible separated flows while retaining the accuracy and robustness of CGKS-5th in the outer region. The reported <1% runtime overhead of the coupling is a concrete practical advantage. The absence of direct verification of the wall-model reconstruction in non-equilibrium regions, however, limits the strength of the evidence for the claimed accuracy on coarsened meshes.

major comments (2)
  1. [Abstract and numerical-results sections] The central claim that the pressure-gradient-based wall model correctly reconstructs the viscous wall stress from exchange-location data supplied by CGKS-5th in adverse-pressure-gradient and separated regions rests on indirect evidence only. No direct comparison is shown between the wall-model-predicted τ_w and the corresponding value extracted from a wall-resolved reference simulation at identical exchange points inside the cylinder separation or the RAE 2822 shock bubble (see Abstract and the two test-case sections).
  2. [Abstract and numerical-results sections] Quantitative error norms (e.g., L2 or L∞ errors in skin friction or velocity profiles), explicit grid dimensions, and direct comparison baselines against wall-resolved CGKS-5th or other reference data are not reported. This makes it impossible to quantify the improvement or to confirm that the observed skin-friction gains are not artifacts of the particular coarsening chosen (Abstract).
minor comments (1)
  1. [Implementation and performance] The statement that the coupling adds less than 1% runtime overhead would benefit from a brief description of the measurement (number of GPUs, problem size, timing method).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive review and for recognizing the practical advantages of the proposed wall-modeled CGKS-5th framework, including the low runtime overhead. We address each major comment below and agree that additional direct evidence will strengthen the manuscript.

read point-by-point responses
  1. Referee: [Abstract and numerical-results sections] The central claim that the pressure-gradient-based wall model correctly reconstructs the viscous wall stress from exchange-location data supplied by CGKS-5th in adverse-pressure-gradient and separated regions rests on indirect evidence only. No direct comparison is shown between the wall-model-predicted τ_w and the corresponding value extracted from a wall-resolved reference simulation at identical exchange points inside the cylinder separation or the RAE 2822 shock bubble (see Abstract and the two test-case sections).

    Authors: We agree that direct comparison of the wall-model-predicted τ_w with reference values extracted at identical exchange locations from wall-resolved simulations would constitute stronger validation of the non-equilibrium wall model. The current evidence is indeed indirect, relying on integrated quantities such as skin-friction distributions and overall flow structures. In the revised manuscript we will post-process the wall-resolved reference fields to obtain τ_w at the precise exchange locations employed in the wall-modeled runs and include these direct comparisons for both the cylinder and RAE 2822 cases. revision: yes

  2. Referee: [Abstract and numerical-results sections] Quantitative error norms (e.g., L2 or L∞ errors in skin friction or velocity profiles), explicit grid dimensions, and direct comparison baselines against wall-resolved CGKS-5th or other reference data are not reported. This makes it impossible to quantify the improvement or to confirm that the observed skin-friction gains are not artifacts of the particular coarsening chosen (Abstract).

    Authors: We concur that the absence of quantitative error norms and explicit grid specifications limits the ability to assess the magnitude and robustness of the reported improvements. The manuscript currently presents results through figures and qualitative statements. We will revise the abstract and numerical-results sections to report L2 and L∞ error norms for skin-friction coefficient and selected velocity profiles (where reference data exist), provide the precise grid dimensions for all meshes, and include direct quantitative comparisons against wall-resolved CGKS-5th solutions. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected; derivation is self-contained.

full rationale

The manuscript presents a new wall-modeled CGKS-5th framework whose central elements—the pressure-gradient-based non-equilibrium wall model and its coupling to the outer solver—are introduced as physically motivated extensions rather than reductions of prior fitted quantities or self-citations. No equation chain shows a prediction (e.g., skin-friction or separation structures) being recovered by construction from parameters fitted to the same data. Self-citations to earlier CGKS work supply the base high-order scheme but do not carry the load-bearing claim about wall-stress reconstruction on coarse meshes. The reported improvements are empirical outcomes on cylinder and RAE 2822 cases, not tautological restatements of inputs. This satisfies the default expectation of an independent derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, axioms beyond standard fluid mechanics, or invented entities are stated.

axioms (1)
  • standard math Compressible Navier-Stokes equations govern the flow
    Implicit background for any gas-kinetic scheme.

pith-pipeline@v0.9.1-grok · 5872 in / 1093 out tokens · 42265 ms · 2026-06-30T03:54:04.704844+00:00 · methodology

0 comments
read the original abstract

Scale-resolving simulations of wall-bounded turbulent flows remain prohibitively expensive at high Reynolds numbers, owing to the stringent near-wall resolution requirements. High-order compact gas-kinetic schemes (CGKS) are accurate, robust, and efficient for compressible flows, making them an attractive foundation for reducing this cost. Building on the fifth-order scheme CGKS-5th, we develop a wall-modeled CGKS framework that alleviates the near-wall resolution burden through a pressure-gradient-based non-equilibrium wall model while preserving the resolving power of the outer solver. CGKS-5th resolves the outer flow and supplies the wall model with data at the exchange location. On coarse near-wall meshes, the wall model reconstructs the under-resolved viscous wall stress, while CGKS-5th provides the inviscid wall flux directly; the two combine to form the wall momentum flux. To capture non-equilibrium effects in adverse-pressure-gradient and separated regions, the wall model retains a pressure-gradient source term together with a pressure-gradient-corrected near-wall damping function. We assess the framework on two distinct flows: bluff-body separation past a circular cylinder, and a shock-induced separation bubble on the transonic RAE 2822 airfoil, using near-wall meshes far coarser than wall-resolved simulations require. For the RAE 2822 case, this corresponds to a twentyfold coarsening in the wallnormal direction, with comparable coarsening in other directions. In both cases, the wall-modeled CGKS-5th reproduces the separated flow structures and markedly improves near-wall predictions over its wall-model-free counterpart, most notably the skin-friction coefficient. The framework thus delivers accurate predictions of these separated flows at substantially reduced near-wall cost, while its lightweight coupling adds less than 1% runtime overhead in a multi-GPU implementation.

Figures

Figures reproduced from arXiv: 2606.30061 by Fengxiang Zhao, Kun Xu, Yaqing Yang.

Figure 1
Figure 1. Figure 1: Turbulent flow past a cylinder: visualization of the computational mesh in the x-y plane with a [PITH_FULL_IMAGE:figures/full_fig_p015_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Turbulent flow past a cylinder: the iso-surface of the Q-criterion ( [PITH_FULL_IMAGE:figures/full_fig_p015_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Turbulent flow past a cylinder: the pressure coefficient (left) and skin-friction coefficient (right) [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Turbulent flow past a cylinder: the streamwise velocity profiles along the centerline in the cylinder [PITH_FULL_IMAGE:figures/full_fig_p016_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Turbulent flow past a cylinder: the streamwise velocity profiles at different cylinder wake locations [PITH_FULL_IMAGE:figures/full_fig_p017_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Turbulent flow past a cylinder: the Reynolds stress components profiles at different cylinder wake [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Turbulent flow past an RAE 2822 airfoil: visualization of the computational mesh in the x-y plane. [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Turbulent flow past an RAE 2822 airfoil: the iso-surface of the Q-criterion ( [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Turbulent flow past an RAE 2822 airfoil: the iso-surface of the Q-criterion ( [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Turbulent flow past an RAE 2822 airfoil: the instantaneous density contour and Mach number [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Turbulent flow past an RAE 2822 airfoil: the instantaneous density gradient magnitude contour [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Turbulent flow past an RAE 2822 airfoil: the time- and spanwise-averaged [PITH_FULL_IMAGE:figures/full_fig_p023_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Turbulent flow past an RAE 2822 airfoil: the time- and spanwise-averaged pressure coefficient [PITH_FULL_IMAGE:figures/full_fig_p023_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Turbulent flow past an RAE 2822 airfoil: the time- and spanwise-averaged skin-friction coefficient [PITH_FULL_IMAGE:figures/full_fig_p024_14.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

49 extracted references

  1. [1]

    Lele, S. K. (1992). Compact finite difference schemes with spectral-like resolution. Journal of compu- tational physics, 103(1), 16-42

  2. [2]

    Mahesh, K. (1998). A family of high order finite difference schemes with good spectral resolution. Journal of Computational Physics, 145(1), 332-358

  3. [3]

    Harten, A., Engquist, B., Osher, S., & Chakravarthy, S. R. (1987). Uniformly high order accurate essentially non-oscillatory schemes, III. Journal of computational physics, 131(1), 3-47

  4. [4]

    W., & Osher, S

    Shu, C. W., & Osher, S. (1988). Efficient implementation of essentially non-oscillatory shock-capturing schemes. Journal of computational physics, 77(2), 439-471

  5. [5]

    D., Osher, S., & Chan, T

    Liu, X. D., Osher, S., & Chan, T. (1994). Weighted essentially non-oscillatory schemes. Journal of computational physics, 115(1), 200-212

  6. [6]

    S., & Shu, C

    Jiang, G. S., & Shu, C. W. (1996). Efficient implementation of weighted ENO schemes. Journal of computational physics, 126(1), 202-228

  7. [7]

    Qiu, J., & Shu, C. W. (2004). Hermite WENO schemes and their application as limiters for Runge–Kutta discontinuous Galerkin method: one-dimensional case. Journal of Computational Physics, 193(1), 115-135

  8. [8]

    H., & Hill, T

    Reed, W. H., & Hill, T. R. (1973). Triangular mesh methods for the neutron transport equation (No. LA-UR–73-479; CONF-730414–2). Los Alamos Scientific Lab., N. Mex. (USA)

  9. [9]

    Cockburn, B., & Shu, C. W. (1989). TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. II. General framework. Mathematics of computation, 52(186), 411-435

  10. [10]

    Huynh, H. T. (2007). A flux reconstruction approach to high-order schemes including discontinuous Galerkin methods. In 18th AIAA computational fluid dynamics conference (p. 4079)

  11. [11]

    J., & Gao, H

    Wang, Z. J., & Gao, H. (2009). A unifying lifting collocation penalty formulation including the discon- tinuous Galerkin, spectral volume/difference methods for conservation laws on mixed grids. Journal of Computational Physics, 228(21), 8161-8186

  12. [12]

    Pan, L., Xu, K., Li, Q., & Li, J. (2016). An efficient and accurate two-stage fourth-order gas-kinetic scheme for the Euler and Navier-Stokes equations. Journal of Computational Physics, 326, 197-221

  13. [13]

    Zhao, F., Ji, X., Shyy, W., & Xu, K. (2019). Compact higher-order gas-kinetic schemes with spectral-like resolution for compressible flow simulations. Advances in Aerodynamics, 1(1), 13

  14. [14]

    Ji, X., Shyy, W., & Xu, K. (2021). A gradient compression-based compact high-order gas-kinetic scheme on 3D hybrid unstructured meshes. International Journal of Computational Fluid Dynamics, 35(7), 485-509

  15. [15]

    Zhao, F., Ji, X., Shyy, W., & Xu, K. (2023). High-order compact gas-kinetic schemes for three- dimensional flow simulations on tetrahedral mesh. Advances in Aerodynamics, 5(1), 1

  16. [16]

    Yang, Y., Zhao, F., & Xu, K. (2026). An effective implementation of high-order compact gas-kinetic scheme on structured meshes for compressible flows. Journal of Computational Physics, 114729

  17. [17]

    Xu, K. (2001). A gas-kinetic BGK scheme for the Navier-Stokes equations and its connection with artificial dissipation and Godunov method. Journal of Computational Physics, 171(1), 289-335

  18. [18]

    Xu, K. (2015). Direct modeling for computational fluid dynamics: construction and application of unified gas-kinetic schemes. World Scientific

  19. [19]

    Xu, K. (2021). A unified computational fluid dynamics framework from rarefied to continuum regimes. Cambridge University Press

  20. [20]

    E., & Spalding, D

    Launder, B. E., & Spalding, D. B. (1983). The numerical computation of turbulent flows. In Numerical prediction of flow, heat transfer, turbulence and combustion (pp. 96-116). Pergamon. 27

  21. [21]

    Spalart, P. R. (2009). Detached-eddy simulation. Annual review of fluid mechanics, 41(1), 181-202

  22. [22]

    Piomelli, U., & Balaras, E. (2002). Wall-layer models for large-eddy simulations. Annual review of fluid mechanics, 34(1), 349-374

  23. [23]

    T., & Park, G

    Bose, S. T., & Park, G. I. (2018). Wall-modeled large-eddy simulation for complex turbulent flows. Annual review of fluid mechanics, 50, 535-561

  24. [24]

    Schumann, U. (1975). Subgrid scale model for finite difference simulations of turbulent flows in plane channels and annuli. Journal of computational physics, 18(4), 376-404

  25. [25]

    Shi, J., Yan, H., & Wang, Z. J. (2020). Flux reconstruction implementation of an algebraic wall model for large-eddy simulation. AIAA Journal, 58(7), 3051-3062

  26. [26]

    Balaras, E., Benocci, C., & Piomelli, U. (1996). Two-layer approximate boundary conditions for large- eddy simulations. AIAA journal, 34(6), 1111-1119

  27. [27]

    Cabot, W., & Moin, P. (2000). Approximate wall boundary conditions in the large-eddy simulation of high Reynolds number flow. Flow, Turbulence and Combustion, 63(1), 269-291

  28. [28]

    Wang, M., & Moin, P. (2002). Dynamic wall modeling for large-eddy simulation of complex turbulent flows. Physics of Fluids, 14(7), 2043-2051

  29. [29]

    R., & Subbareddy, P

    Mettu, B. R., & Subbareddy, P. K. (2022). Wall-modeled large eddy simulation of high speed flows. AIAA journal, 60(7), 4302-4324

  30. [30]

    M., & Brugi` ere, O

    Duprat, C., Balarac, G., M´ etais, O., Congedo, P. M., & Brugi` ere, O. (2011). A wall-layer model for large-eddy simulations of turbulent flows with/out pressure gradient. Physics of fluids, 23(1)

  31. [31]

    L., Gross, E

    Bhatnagar, P. L., Gross, E. P., & Krook, M. (1954). A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component systems. Physical review, 94(3), 511

  32. [32]

    Chapman, S., & Cowling, T. G. (1990). The mathematical theory of non-uniform gases: an account of the kinetic theory of viscosity, thermal conduction and diffusion in gases. Cambridge university press

  33. [33]

    Zhao, F., Ji, X., Shyy, W., & Xu, K. (2022). A compact high-order gas-kinetic scheme on unstructured mesh for acoustic and shock wave computations. Journal of Computational Physics, 449, 110812

  34. [34]

    Zhao, F., & Xu, K. (2025). A generalized ENO reconstruction in compact GKS for compressible flow simulations. Journal of Computational Physics, 114612

  35. [35]

    C., G¨ u¸ cl¨ u, Y., & Christlieb, A

    Seal, D. C., G¨ u¸ cl¨ u, Y., & Christlieb, A. J. (2014). High-order multiderivative time integrators for hyperbolic conservation laws. Journal of Scientific Computing, 60(1), 101-140

  36. [36]

    Li, J., & Du, Z. (2016). A two-stage fourth order time-accurate discretization for Lax–Wendroff type flow solvers I. Hyperbolic conservation laws. SIAM Journal on Scientific Computing, 38(5), A3046- A3069

  37. [37]

    Li, J. (2019). Two-stage fourth order: temporal-spatial coupling in computational fluid dynamics (CFD). Advances in Aerodynamics, 1(1), 3

  38. [38]

    Zhao, F., Ji, X., Shyy, W., & Xu, K. (2023). Direct modeling for computational fluid dynamics and the construction of high-order compact scheme for compressible flow simulations. Journal of Computational Physics, 477, 111921

  39. [39]

    Norberg, C. (1993). Pressure forces on a circular cylinder in cross flow. In Bluff-Body Wakes, Dynamics and Instabilities: IUTAM Symposium, Springer Berlin Heidelberg

  40. [40]

    S., & Karniadakis, G

    Ma, X., Karamanos, G. S., & Karniadakis, G. E. (2000). Dynamics and low-dimensionality of a turbu- lent near wake. Journal of fluid mechanics, 410, 29-65

  41. [41]

    I., Samtaney, R., Zhang, W., & Gao, W

    Cheng, W., Pullin, D. I., Samtaney, R., Zhang, W., & Gao, W. (2017). Large-eddy simulation of flow over a cylinder withRe D from 3.9×10 3 to 8.5×10 5: A skin-friction perspective. Journal of Fluid Mechanics, 820, 121-158

  42. [42]

    S., Matsuno, K

    Song, H., Ghate, A. S., Matsuno, K. V., West, J. R., Subramaniam, A., & Lele, S. K. (2024). A robust compact finite difference framework for simulations of compressible turbulent flows. Journal of Computational Physics, 519, 113419

  43. [43]

    M., Mikheev, N

    Molochnikov, V. M., Mikheev, N. I., Mikheev, A. N., Paereliy, A. A., Dushin, N. S., & Dushina, O. A. (2019). SIV measurements of flow structure in the near wake of a circular cylinder at Re= 3900. Fluid Dynamics Research, 51(5), 055505

  44. [44]

    Norberg, C. (1994). An experimental investigation of the flow around a circular cylinder: influence of 28 aspect ratio. Journal of Fluid Mechanics, 258, 287-316

  45. [45]

    Parnaudeau, P., Carlier, J., Heitz, D., & Lamballais, E. (2008). Experimental and numerical studies of the flow over a circular cylinder at Reynolds number 3900. Physics of fluids, 20(8)

  46. [46]

    H., Firmin, M

    Cook, P. H., Firmin, M. C. P., & McDonald, M. A. (1977). Aerofoil RAE 2822: pressure distributions, and boundary layer and wake measurements. RAE

  47. [47]

    Haase, W., Bradsma, F., Elsholz, E., Leschziner, M., & Schwamborn, D. (1993). EUROVAL-An Eu- ropean Initiative on Validation of CFD Codes. Notes on Numerical Fluid Mechanics, Vol. 42. Vieweg, Braunschweig/Wiebaden

  48. [48]

    J., & She, Z

    Xiao, M. J., & She, Z. S. (2020). Precise drag prediction of airfoil flows by a new algebraic model. Acta Mechanica Sinica, 36(1), 35

  49. [49]

    L., Fu, D

    Liang, X., Li, X. L., Fu, D. X., & Ma, Y. W. (2011). Ten thousand order cores expandable CFD software and its application. Journal of Huazhong University of Science and Technology (Natural Science Edition), 39(Sup. I), 67-70. (in Chinese) 29