pith. sign in

arxiv: 2512.06433 · v2 · submitted 2025-12-06 · ⚛️ physics.flu-dyn

Model of incompressible turbulent flows via a kinetic theory

Ziyang Xin , Zhaoli Guo , Hudong Chen This is my paper

Pith reviewed 2026-05-17 00:48 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn
keywords kinetic theoryturbulence modelingChapman-Enskog expansionincompressible flowReynolds stresseddy viscositywall-bounded turbulenceCouette flow
0
0 comments X

The pith

A kinetic model for incompressible turbulent flows derives linear eddy-viscosity closures at first order in Chapman-Enskog expansion and nonlinear corrections plus a new turbulent kinetic energy flux expression at second order.

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

This paper extends a prior kinetic model for unbounded incompressible turbulent flows by reselecting the relaxation time to align turbulent transport coefficients more closely with established theory. Chapman-Enskog analysis of the revised model recovers the conventional linear eddy-viscosity and gradient-diffusion closures for Reynolds stress and turbulent kinetic energy flux at first order while producing nonlinear eddy-viscosity models and a previously unreported solution for the turbulent kinetic energy flux at second order. The work further adapts the model to wall-bounded flows through added low-Reynolds-number damping and viscous-diffusion terms that maintain required near-wall asymptotic behavior, enabling both resolved viscous-sublayer and wall-function treatments. Validation on turbulent plane Couette flow shows close agreement with experimental and DNS data for mean velocity, skin friction, and Reynolds stresses. The approach frames averaged turbulent flow as analogous to a rarefied gas at finite Knudsen number, thereby capturing non-Newtonian effects that linear eddy-viscosity models miss.

Core claim

The Chapman-Enskog analysis of the kinetic model reproduces the traditional linear eddy viscosity and gradient diffusion models for Reynolds stress and turbulent kinetic energy flux at the first order, and yields nonlinear eddy viscosity and closure models at the second order. Particularly, a previously unreported CE solution for turbulent kinetic energy flux is obtained. The second extension is to enable the model for wall-bounded turbulent flows with preserved near-wall asymptotic behaviours through a low-Reynolds number kinetic model incorporating wall damping effects and viscous diffusion.

What carries the argument

Chapman-Enskog expansion of a Boltzmann-like kinetic equation for turbulent flows after reselecting the relaxation time

If this is right

  • First-order Chapman-Enskog solution recovers the standard linear eddy-viscosity model for Reynolds stresses.
  • Second-order solution supplies nonlinear eddy-viscosity corrections and a new closure relation for turbulent kinetic energy flux.
  • Wall-bounded extension reproduces correct near-wall asymptotics for both viscous sublayer resolution and wall-function application.
  • Plane Couette flow simulations match measured mean-velocity profiles, skin-friction coefficients, and Reynolds-stress distributions.

Where Pith is reading between the lines

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

  • Treating averaged turbulence as a finite-Knudsen-number rarefied-gas problem suggests a route to derive non-Newtonian constitutive relations without additional empirical fitting.
  • Higher-order terms obtained from the same kinetic framework could be tested directly against DNS budgets for production, dissipation, and transport in channel or pipe flows.
  • The mesoscopic description offers a natural way to incorporate compressibility or heat-transfer effects by extending the underlying distribution function.

Load-bearing premise

Reselection of the relaxation time produces transport coefficients consistent with established turbulence theory and the added low-Reynolds-number damping plus viscous-diffusion terms preserve near-wall asymptotic behaviors without introducing inconsistencies.

What would settle it

Quantitative mismatch between the second-order nonlinear Reynolds-stress components or the newly derived turbulent kinetic energy flux and corresponding DNS or experimental measurements at moderate-to-high Reynolds numbers in simple shear flows.

Figures

Figures reproduced from arXiv: 2512.06433 by Hudong Chen, Zhaoli Guo, Ziyang Xin.

Figure 1
Figure 1. Figure 1: Sketch of turbulent plane Couette flow between two infinite parallel plates separated [PITH_FULL_IMAGE:figures/full_fig_p014_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Mean velocity profiles scaled in (a) outer and (b) wall units of the turbulent Couette flow. [PITH_FULL_IMAGE:figures/full_fig_p016_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Mean velocity profiles scaled in (a) outer and (b) wall units of turbulent Couette flow. The [PITH_FULL_IMAGE:figures/full_fig_p017_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Variation of skin friction coefficient with Reynolds number. The solid and dash-dot [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Reynolds shear stress profiles scaled in wall units of turbulent Couette flow at (a) [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Distribution of the velocity fluctuation intensities (Reynolds normal stress) scaled [PITH_FULL_IMAGE:figures/full_fig_p019_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Distribution of the velocity fluctuation intensities (Reynolds normal stress) scaled [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The reduced velocity distribution functions [PITH_FULL_IMAGE:figures/full_fig_p022_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: The reduced velocity distribution functions [PITH_FULL_IMAGE:figures/full_fig_p023_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Deviatoric Reynolds shear stress component in wall units of turbulent Couette flow [PITH_FULL_IMAGE:figures/full_fig_p024_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Deviatoric Reynolds normal stress components in wall units of turbulent Couette flow [PITH_FULL_IMAGE:figures/full_fig_p024_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: The reduced velocity distribution function [PITH_FULL_IMAGE:figures/full_fig_p032_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: The reduced velocity distribution function [PITH_FULL_IMAGE:figures/full_fig_p033_13.png] view at source ↗
read the original abstract

Kinetic theory offers a promising alternative to conventional turbulence modelling by providing a mesoscopic perspective that naturally captures non-equilibrium physics such as non-Newtonian effects. In this work, we present an extension and theoretical analysis of the recent kinetic model for incompressible turbulent flows developed by Chen et al. (Atmos. 14(7), 1109, 2023), constructed for unbounded flows. The first extension is to reselect a relaxation time such that the turbulent transport coefficients are obtained more consistently and better align with well-established turbulence theory. The Chapman-Enskog (CE) analysis of the kinetic model reproduces the traditional linear eddy viscosity and gradient diffusion models for Reynolds stress and turbulent kinetic energy flux at the first order, and yields nonlinear eddy viscosity and closure models at the second order. Particularly, a previously unreported CE solution for turbulent kinetic energy flux is obtained. The second extension is to enable the model for wall-bounded turbulent flows with preserved near-wall asymptotic behaviours. This involves developing a low-Reynolds number kinetic model incorporating wall damping effects and viscous diffusion, with boundary conditions enabling both viscous sublayer resolution and wall function application. Comprehensive validation against experimental and DNS data for turbulent plane Couette flow demonstrates excellent agreement in predicting mean velocity profiles, skin friction coefficients, and Reynolds stress distributions. It reveals that an averaged turbulent flow behaves similarly to a rarefied gas flow at a finite Knudsen number, capturing non-Newtonian effects inaccessible to linear eddy viscosity models. This kinetic model provides a physics-based foundation for turbulence modelling with reduced empirical dependence.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript extends the kinetic model of Chen et al. (Atmos. 14(7), 1109, 2023) for incompressible turbulent flows. It reselects the relaxation time to obtain transport coefficients more consistent with established turbulence theory, performs a Chapman-Enskog expansion that recovers the linear eddy-viscosity and gradient-diffusion closures at first order while producing nonlinear closures and a new expression for turbulent kinetic energy flux at second order, and adds low-Reynolds-number wall-damping plus viscous-diffusion terms to enable wall-bounded simulations. Boundary conditions support both viscous-sublayer resolution and wall-function use. Validation on plane Couette flow is reported to show excellent agreement with experimental and DNS data for mean velocity, skin friction, and Reynolds stresses, with the model capturing non-Newtonian effects.

Significance. If the Chapman-Enskog derivations hold without hidden assumptions, the work supplies a mesoscopic, physics-based route to turbulence closures that naturally incorporates non-equilibrium behavior and reduces empirical content relative to conventional RANS models. The explicit recovery of both linear and nonlinear forms plus the wall extensions constitute a coherent framework whose practical utility is supported by the Couette-flow comparisons.

major comments (2)
  1. [§3] §3 (Chapman-Enskog analysis): the second-order solution for the turbulent kinetic energy flux is presented as previously unreported and central to the nonlinear closure claim. An independent re-derivation is required to confirm that this expression follows directly from the kinetic equation and the reselected relaxation time without unstated moment closures or ordering assumptions; this step is load-bearing for the novelty and consistency assertions.
  2. [§5] §5 (Validation): the reported 'excellent agreement' with experimental and DNS data for plane Couette flow lacks quantitative error metrics, error-bar reporting, or explicit data-exclusion criteria. Without these, it is impossible to judge whether the model reproduces the data within stated uncertainties or whether post-hoc adjustments were applied, undermining the practical validation of the central claim.
minor comments (2)
  1. The notation and explicit functional form of the reselected relaxation time should be compared side-by-side with the original Chen et al. choice, preferably in a table, to make the improvement transparent.
  2. Boundary-condition implementation for viscous-sublayer resolution versus wall-function application would benefit from a short dedicated subsection or appendix with the precise discrete equations used at the wall.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed review of our manuscript. The comments have helped us identify areas where additional rigor and clarity can strengthen the presentation. We address each major comment point by point below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [§3] §3 (Chapman-Enskog analysis): the second-order solution for the turbulent kinetic energy flux is presented as previously unreported and central to the nonlinear closure claim. An independent re-derivation is required to confirm that this expression follows directly from the kinetic equation and the reselected relaxation time without unstated moment closures or ordering assumptions; this step is load-bearing for the novelty and consistency assertions.

    Authors: We appreciate the referee's emphasis on verifying the Chapman-Enskog derivation. The second-order expression for the turbulent kinetic energy flux is obtained directly from the kinetic equation by applying the standard Chapman-Enskog expansion to the reselected relaxation time, without introducing additional moment closures or non-standard ordering assumptions. To fully address this point, the revised manuscript will include a complete, step-by-step re-derivation in a new appendix, explicitly tracing each term from the kinetic model through the expansion. This will confirm that the result follows rigorously and support the claims of novelty and consistency with established turbulence theory. revision: yes

  2. Referee: [§5] §5 (Validation): the reported 'excellent agreement' with experimental and DNS data for plane Couette flow lacks quantitative error metrics, error-bar reporting, or explicit data-exclusion criteria. Without these, it is impossible to judge whether the model reproduces the data within stated uncertainties or whether post-hoc adjustments were applied, undermining the practical validation of the central claim.

    Authors: We agree that quantitative metrics are necessary for a rigorous assessment of the validation results. In the revised manuscript, we will add explicit quantitative error measures, including relative errors, root-mean-square deviations, and comparisons against reported uncertainties in the experimental and DNS datasets. We will also state the data sources, any selection criteria applied, and confirm that no post-hoc parameter adjustments were made; all coefficients were determined from the theoretical analysis and established values in the literature. These additions will allow readers to evaluate the agreement objectively. revision: yes

Circularity Check

1 steps flagged

Reselection of relaxation time forces first-order reproduction of standard turbulence closures in CE analysis

specific steps
  1. fitted input called prediction [Abstract]
    "The first extension is to reselect a relaxation time such that the turbulent transport coefficients are obtained more consistently and better align with well-established turbulence theory. The Chapman-Enskog (CE) analysis of the kinetic model reproduces the traditional linear eddy viscosity and gradient diffusion models for Reynolds stress and turbulent kinetic energy flux at the first order, and yields nonlinear eddy viscosity and closure models at the second order. Particularly, a previously unreported CE solution for turbulent kinetic energy flux is obtained."

    The relaxation time is chosen explicitly to force alignment of transport coefficients with established turbulence theory. The subsequent claim that CE analysis 'reproduces' the linear eddy-viscosity and gradient-diffusion models is then a direct outcome of that choice, rendering the reproduction a fitted result renamed as a first-principles prediction rather than an independent derivation from the kinetic equation.

full rationale

The paper reselects the relaxation time parameter specifically so that the resulting turbulent transport coefficients align with well-established theory, after which the Chapman-Enskog expansion is shown to reproduce the linear eddy-viscosity and gradient-diffusion forms. This reproduction is therefore a direct consequence of the parameter choice rather than an independent derivation. The base kinetic model is taken from the self-cited Chen et al. (2023) work by overlapping authors, so the overall framework and the first-order results reduce to the adjusted inputs by construction. The second-order nonlinear terms and the new TKE-flux expression may contain independent content, but the load-bearing first-order claim does not.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on reselecting a single relaxation time to enforce consistency with established turbulence closures and on adding ad-hoc wall-damping and viscous-diffusion terms whose functional forms are chosen to recover known near-wall asymptotics.

free parameters (1)
  • relaxation time
    Reselected so that the resulting turbulent transport coefficients align with well-established turbulence theory.
axioms (1)
  • domain assumption Standard kinetic-theory assumptions for incompressible flow (BGK-type collision operator, Chapman-Enskog expansion validity)
    Invoked as the foundation for both the unbounded and wall-bounded extensions.

pith-pipeline@v0.9.0 · 5578 in / 1561 out tokens · 55172 ms · 2026-05-17T00:48:01.856363+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

  1. [1]

    Cercignani, Carlo1988 The Boltzmann equation

    Cercignani, Carlo1972 Scattering kernels for gas-surface interactions.Transport Theory and Statistical Physics2(1), 27–53. Cercignani, Carlo1988 The Boltzmann equation. In The Boltzmann Equation and Its Applications, pp. 40–103. Springer. Cercignani, Carlo & Lampis, Maria1971 Kinetic models for gas-surface interactions.Transport Theory and Statistical Phy...

    work page 2003

  2. [2]

    Launder, Brian Edward, Reece, G Jr & Rodi, W1975 Progress in the development of a Reynolds-stress turbulence closure.Journal of Fluid Mechanics68(3), 537–566

    Kosuge, Shingo, Aoki, Kazuo, Giovangigli, Vincent & Golse, Franc ¸ois2025 Applications of new boundary conditions for the Boltzmann equation derived from a kinetic model of gas-surface interaction.Physical Review Fluids10(5), 053401. Launder, Brian Edward, Reece, G Jr & Rodi, W1975 Progress in the development of a Reynolds-stress turbulence closure.Journa...

    work page 2020

  3. [3]

    Xu, Kun, Liu, Hongwei & Jiang, Jianzheng2007 Multiple-temperature kinetic model for continuum and near continuum flows.Physics of Fluids19(1), 016101

    Wu, Lei, Reese, Jason M & Zhang, Yonghao2014 Solving the Boltzmann equation deterministically by the fast spectral method: application to gas microflows.Journal of Fluid Mechanics746, 53–84. Xu, Kun, Liu, Hongwei & Jiang, Jianzheng2007 Multiple-temperature kinetic model for continuum and near continuum flows.Physics of Fluids19(1), 016101. Yakhot, VSASTBC...

    work page arXiv