pith. sign in

arxiv: 1907.09191 · v1 · pith:EWG3MALRnew · submitted 2019-07-22 · 🧮 math.AP

Turbulent flows as generalized Kelvin-Voigt materials: modeling and analysis

Cherif Amrouche , Luigi C. Berselli , Roger Lewandowski , Dinh Duong Nguyen This is my paper

Pith reviewed 2026-05-24 18:22 UTC · model grok-4.3

classification 🧮 math.AP
keywords turbulent flowsKelvin-Voigt modelPrandtl mixing lengthReynolds averaged Navier-Stokesweak solutionsregularity estimatesstructural compactness
0
0 comments X

The pith

Turbulent flows modeled with a position-dependent Kelvin-Voigt term admit unique regular-weak solutions for the mean velocity and weak solutions for the full Reynolds-averaged system.

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

The paper augments the Navier-Stokes equations for the mean velocity and pressure with a Kelvin-Voigt-type term that involves the time derivative of the deformation tensor multiplied by a spatially varying Prandtl mixing length vanishing at the walls. This modeling choice produces a priori estimates placing the velocity in the intersection of L^∞ in time of H¹ and W^{1,2} in time of H^{1/2}. The estimates suffice to prove existence and uniqueness of regular-weak solutions when the eddy viscosity is held fixed. A structural compactness property of the model then permits passage to the limit inside the quadratic source term of the turbulent kinetic energy equation, establishing existence of a weak solution to the coupled Reynolds-averaged Navier-Stokes system.

Core claim

By treating the turbulent fluid as a generalized Kelvin-Voigt material with a wall-vanishing mixing length, the authors obtain estimates that guarantee unique regular-weak solutions to the augmented mean-field equations and, via compactness, weak solutions to the full RANS system including the turbulent kinetic energy.

What carries the argument

The term −α ∇·(ℓ(x) D v_t) where ℓ(x) is the Prandtl mixing length; it provides the additional viscous dissipation needed for the regularity estimates and compactness.

Load-bearing premise

The premise that the added Kelvin-Voigt term with the given mixing length accurately represents the statistical equilibrium reached by the turbulent flow.

What would settle it

A counter-example in which the quadratic source term for k fails to converge weakly despite the velocity estimates, or a physical flow in which the predicted damping near walls does not match measured strain-rate statistics.

read the original abstract

We model a 3D turbulent fluid, evolving toward a statistical equilibrium, by adding to the equations for the mean field $(v, p)$ a term like $-\alpha \nabla\cdot(\ell(x) D v_t)$. This is of the Kelvin-Voigt form, where the Prandtl mixing length $\ell$ is not constant and vanishes at the solid walls. We get estimates for velocity $v$ in $L^\infty_t H^1_x \cap W^{1,2}_t H^{1/2}_x$, that allow us to prove the existence and uniqueness of a regular-weak solutions $(v, p)$ to the resulting system, for a given fixed eddy viscosity. We then prove a structural compactness result that highlights the robustness of the model. This allows us to pass to the limit in the quadratic source term in the equation for the turbulent kinetic energy $k$, which yields the existence of a weak solution to the corresponding Reynolds Averaged Navier-Stokes system satisfied by $(v, p, k)$.

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

3 major / 1 minor

Summary. The paper models 3D turbulent fluids evolving toward statistical equilibrium by augmenting the mean-field Navier-Stokes equations for (v, p) with the Kelvin-Voigt term −α ∇·(ℓ(x) D v_t), where ℓ(x) is the Prandtl mixing length vanishing at walls. It derives a priori estimates placing v in L^∞_t H^1_x ∩ W^{1,2}_t H^{1/2}_x, proves existence and uniqueness of regular-weak solutions (v, p) for fixed eddy viscosity, establishes a structural compactness result, and passes to the limit in the quadratic source term of the k-equation to obtain a weak solution of the coupled RANS system for (v, p, k).

Significance. If the derivations hold, the work supplies rigorous existence/uniqueness for the augmented system and a compactness argument that justifies passage to the limit in the nonlinear k-source term. The structural compactness result is a clear strength, as it demonstrates robustness independent of specific approximations. The physical interpretation as a model for turbulent statistical equilibrium, however, rests entirely on the unverified modeling choice for the extra term.

major comments (3)
  1. [Abstract / Introduction] Modeling premise (abstract and introduction): the claim that −α ∇·(ℓ(x) D v_t) encodes the statistical equilibrium of unresolved fluctuations is load-bearing for applying the existence results to physical RANS, yet no derivation from the Navier-Stokes equations, filtered equations, or statistical closure is supplied and no concrete test (e.g., recovery of a standard eddy-viscosity closure or comparison with benchmark RANS solutions) is performed. Without such support the mathematical statements remain valid for the artificial system but do not automatically transfer to turbulent flows.
  2. [Existence section (energy estimates)] Energy estimates and variable coefficient (section containing the a priori bounds): the claimed regularity L^∞_t H^1_x ∩ W^{1,2}_t H^{1/2}_x is stated to follow from standard Sobolev estimates, but the degeneracy of ℓ(x) at the walls must be handled explicitly; if the estimates treat ℓ as bounded below in the interior only, the boundary-layer behavior and possible loss of coercivity require a separate argument that is not visible from the abstract sketch.
  3. [Compactness and limit passage section] Compactness and passage to the limit (section on structural compactness): while the compactness result is invoked to justify the limit inside the quadratic source term of the k-equation, the precise mode of convergence (strong in which space?) and the integrability needed to pass the limit must be verified against the variable-coefficient Kelvin-Voigt term; any gap here would block the final existence statement for the RANS system.
minor comments (1)
  1. [Notation / Preliminaries] Clarify the precise definition of the deformation tensor D and the notation v_t throughout; inconsistent use could obscure the time-derivative structure of the Kelvin-Voigt term.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight points where the manuscript can be clarified and strengthened. We respond to each major comment below and indicate the revisions we will make.

read point-by-point responses

  1. Referee: [Abstract / Introduction] Modeling premise (abstract and introduction): the claim that −α ∇·(ℓ(x) D v_t) encodes the statistical equilibrium of unresolved fluctuations is load-bearing for applying the existence results to physical RANS, yet no derivation from the Navier-Stokes equations, filtered equations, or statistical closure is supplied and no concrete test (e.g., recovery of a standard eddy-viscosity closure or comparison with benchmark RANS solutions) is performed. Without such support the mathematical statements remain valid for the artificial system but do not automatically transfer to turbulent flows.

    Authors: We agree that the manuscript presents the Kelvin-Voigt term as a modeling choice rather than a derived closure. The paper is devoted to the mathematical analysis of the resulting system and does not contain a first-principles derivation from the Navier-Stokes equations or numerical validation against standard RANS closures. We will revise the abstract and introduction to state explicitly that the extra term is a proposed regularization motivated by the Prandtl mixing length and Kelvin-Voigt viscoelasticity, without claiming a direct statistical derivation. The existence and compactness results therefore apply rigorously to the augmented equations; their relevance to physical turbulence remains a modeling question outside the scope of the present work. revision: partial

  2. Referee: [Existence section (energy estimates)] Energy estimates and variable coefficient (section containing the a priori bounds): the claimed regularity L^∞_t H^1_x ∩ W^{1,2}_t H^{1/2}_x is stated to follow from standard Sobolev estimates, but the degeneracy of ℓ(x) at the walls must be handled explicitly; if the estimates treat ℓ as bounded below in the interior only, the boundary-layer behavior and possible loss of coercivity require a separate argument that is not visible from the abstract sketch.

    Authors: The referee correctly notes that the degeneracy of ℓ(x) at the walls must be treated with care to preserve coercivity. While the full proof in the manuscript uses the positivity of ℓ in the interior together with the no-slip boundary condition, the argument is not written out in sufficient detail. We will add an explicit subsection deriving the a priori estimates, showing how the variable coefficient is handled near the boundary and confirming that the claimed regularity holds without loss of coercivity. revision: yes

  3. Referee: [Compactness and limit passage section] Compactness and passage to the limit (section on structural compactness): while the compactness result is invoked to justify the limit inside the quadratic source term of the k-equation, the precise mode of convergence (strong in which space?) and the integrability needed to pass the limit must be verified against the variable-coefficient Kelvin-Voigt term; any gap here would block the final existence statement for the RANS system.

    Authors: We will expand the structural-compactness section to state the precise convergence: the velocity sequence converges strongly in L^2(0,T;H^1(Ω)) and weakly in the spaces given by the a priori bounds. Because the Kelvin-Voigt term is linear in the time derivative and the coefficient ℓ(x) is fixed and bounded, the compactness carries over directly. We will verify the required integrability of the quadratic source term explicitly, confirming that the limit passage is justified under the obtained regularity. revision: yes

Circularity Check

0 steps flagged

No circularity; existence/uniqueness and compactness follow from standard Sobolev estimates on the augmented system

full rationale

The paper introduces the Kelvin-Voigt term −α ∇·(ℓ(x) D v_t) as an explicit modeling choice to represent statistical equilibrium, then derives a priori estimates in L^∞_t H¹_x ∩ W^{1,2}_t H^{1/2}_x, proves existence/uniqueness of regular-weak solutions for fixed eddy viscosity, and obtains structural compactness to pass to the limit in the k-equation. All steps rely on classical functional-analysis arguments applied to the chosen PDE; no quantity is fitted to data and then relabeled as a prediction, no self-citation supplies a load-bearing uniqueness theorem, and no ansatz is smuggled via prior work. The modeling premise itself is external to the mathematical derivation and does not create a definitional loop inside the proofs.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 1 invented entities

The central claims rest on the unverified modeling hypothesis that the added Kelvin-Voigt term captures turbulent statistics and on standard but unlisted functional-analytic background results such as compactness embeddings in the relevant Sobolev spaces.

free parameters (1)
  • alpha
    Scaling coefficient in the added Kelvin-Voigt term; its value is not derived from first principles.
axioms (2)
  • domain assumption The Prandtl mixing length ℓ(x) vanishes at solid walls and the eddy viscosity is fixed and given.
    This is invoked to obtain the velocity estimates and is stated in the abstract as part of the model.
  • standard math Standard Sobolev embeddings and compactness theorems apply to the variable-coefficient term.
    Required for the existence and limit passage arguments but not proved in the abstract.
invented entities (1)
  • The term −α ∇·(ℓ(x) D v_t) as a model for turbulent statistical equilibrium no independent evidence
    purpose: Regularization that enables a priori estimates and compactness for the mean-field equations
    Introduced as a modeling choice without independent experimental or first-principles derivation in the abstract.

pith-pipeline@v0.9.0 · 5721 in / 1651 out tokens · 29842 ms · 2026-05-24T18:22:52.151533+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

30 extracted references · 30 canonical work pages

  1. [1]

    Amrouche, P

    C. Amrouche, P. G. Ciarlet, L. Gratie, and S. Kesavan. On t he characterizations of matrix fields as linearized strain tensor fields. J. Math. Pures Appl. (9) , 86(2):116– 132, 2006

    work page 2006

  2. [2]

    Berselli and R

    L.C. Berselli and R. Lewandowski. On the Reynolds time-a veraged equations and the long-time behavior of Leray-Hopf weak solutions, with a pplications to ensemble averages. Nonlinearity, to appear. Tech. Report on arXiv:1801.08721, 2019

    work page arXiv 2019

  3. [3]

    L. C. Berselli, T.-Y. Kim, and L. G. Rebholz. Analysis of a reduced-order approximate deconvolution model and its interpretation as a Navier-Sto kes-Voigt regularization. Discrete Contin. Dyn. Syst. Ser. B , 21(4):1027–1050, 2016

    work page 2016

  4. [4]

    Boussinesq

    J. Boussinesq. Theorie de l’´ ecoulement tourbillant. M´ em. pr´ es par div. savants ´ a la Acad. Sci., 23:46–50, 1877

    work page

  5. [5]

    Chac´ on-Rebollo and R

    T. Chac´ on-Rebollo and R. Lewandowski.Mathematical and Numerical Foundations of Turbulence Models and Applications. Modeling and Simulation in Science, Engineering and Technology. Springer New York, 2014

    work page 2014

  6. [6]

    P. G. Ciarlet, M. Malin, and C. Mardare. On a vector versio n of a fundamental lemma of J. L. Lions. Chin. Ann. Math. Ser. B , 39(1):33–46, 2018

    work page 2018

  7. [7]

    E. Emmrich. Discrete versions of Gronwall’s lemma and th eir application to the numerical analysis of parabolic problems. Technical repor t, 1999 available at

    work page 1999

  8. [8]

    P. Germain. M´ ecanique des milieux continus. Masson et cie, ´Editeurs, Paris, 1962, 1962

    work page 1962

  9. [9]

    Girault and P.-A

    V. Girault and P.-A. Raviart. Finite element methods for Navier-Stokes equations , volume 5 of Springer Series in Computational Mathematics . Springer-Verlag, Berlin,

    work page

  10. [10]

    Theory and algorithms

    work page

  11. [11]

    M. E. Gurtin. An introduction to continuum mechanics , volume 158 of Mathemat- ics in Science and Engineering . Academic Press, Inc. [Harcourt Brace Jovanovich, Publishers], New York-London, 1981

    work page 1981

  12. [12]

    V. K. Kalantarov and E. S. Titi. Global attractors and de termining modes for the 3D Navier-Stokes-Voigt equations. Chin. Ann. Math. Ser. B , 30(6):697–714, 2009

    work page 2009

  13. [13]

    Larios and E

    A. Larios and E. S. Titi. On the higher-order global regu larity of the inviscid Voigt- regularization of three-dimensional hydrodynamic models . Discrete Contin. Dyn. Syst. Ser. B , 14(2):603–627, 2010. 23

    work page 2010

  14. [14]

    Layton and R

    W. Layton and R. Lewandowski. On a well-posed turbulenc e model. Discrete Contin. Dyn. Syst. Ser. B , 6(1):111–128, 2006

    work page 2006

  15. [15]

    J. Leray. Sur le mouvement d’un liquide visqueux emplis sant l’espace. Acta Math. , 63(1):193–248, 1934

    work page 1934

  16. [16]

    Levant, F

    B. Levant, F. Ramos, and E. S. Titi. On the statistical pr operties of the 3D incom- pressible Navier-Stokes-Voigt model. Commun. Math. Sci. , 8(1):277–293, 2010

    work page 2010

  17. [17]

    Lewandowski

    R. Lewandowski. The mathematical analysis of the coupl ing of a turbulent kinetic energy equation to the Navier-Stokes equation with an eddy v iscosity. Nonlinear Anal., 28(2):393–417, 1997

    work page 1997

  18. [18]

    Lewandowski

    R. Lewandowski. Navier-stokes equations in the whole s pace with an eddy viscosity. J. Math. Anal. Appl. , 2019 Online first 10.1016/j.jmaa.2019.05.051

    work page doi:10.1016/j.jmaa.2019.05.051 2019

  19. [19]

    Lewandowski, B

    R. Lewandowski, B. Pinier, E. Memin, and Chandramouli P . Testing a one-closure equation turbulence model in neutral boundary layers. Subm itted. Paper avaible at hal-01875464, 2018

    work page 2018

  20. [20]

    Lions and E

    J.-L. Lions and E. Magenes. Non-homogeneous boundary value problems and appli- cations. Vol. I . Springer-Verlag, New York-Heidelberg, 1972. Translated from the French by P. Kenneth, Die Grundlehren der mathematischen Wi ssenschaften, Band 181

    work page 1972

  21. [21]

    Neˇ cas and I

    J. Neˇ cas and I. Hlav´ aˇ cek.Mathematical theory of elastic and elasto-plastic bodies: an introduction, volume 3 of Studies in Applied Mechanics . Elsevier Scientific Publishing Co., Amsterdam-New York, 1980

    work page 1980

  22. [22]

    A. M. Obuhov. Turbulence in an atmosphere with inhomoge neous temperature. Akad. Nauk SSSR. Trudy Inst. Teoret. Geofiz. , 1:95–115, 1946

    work page 1946

  23. [23]

    A. P. Oskolkov. The uniqueness and solvability in the la rge of boundary value prob- lems for the equations of motion of aqueous solutions of poly mers. Zap. Nauˇ cn. Sem. Leningrad. Otdel. Mat. Inst. Steklov. (LOMI) , 38:98–136, 1973. Boundary value prob- lems of mathematical physics and related questions in the th eory of functions, 7

    work page 1973

  24. [24]

    A. P. Oskolkov. On the theory of Voight fluids. Zap. Nauchn. Sem. Leningrad. Otdel. Mat. Inst. Steklov. (LOMI) , 96:233–236, 310, 1980. Boundary value problems of mathematical physics and related questions in the theory of functions, 12

    work page 1980

  25. [25]

    L. Prandtl. Prandtl—Essentials of fluid mechanics , volume 158 of Applied Mathe- matical Sciences. Springer, New York, third edition, 2010. Translated from t he 12th German edition by Katherine Asfaw and edited by Herbert Oert el, With contributions by P. Erhard, D. Etling, U. M¨ uller, U. Riedel, K. R. Sreeniva san and J. Warnatz

    work page 2010

  26. [26]

    Ramos and E

    F. Ramos and E. S. Titi. Invariant measures for the 3D Nav ier-Stokes-Voigt equations and their Navier-Stokes limit. Discrete Contin. Dyn. Syst. , 28(1):375–403, 2010

    work page 2010

  27. [27]

    Rong, W.J

    Y. Rong, W.J. Layton, and H. Zhao. Extension of a simplifi ed Baldwin-Lomax model to nonequilibrium turbulence: Model, analysis and algorit hms. Numer. Methods Par- tial Differential Equations , 5 2019

    work page 2019

  28. [28]

    R. Temam. Navier-Stokes equations. AMS Chelsea Publishing, Providence, RI, 2001. Theory and numerical analysis, Reprint of the 1984 edition. 24

    work page 2001

  29. [29]

    E. R. Van Driest. On turbulent flow near a wall. J. Aeronaut. Sci. , 23(11):1007–1011, 1956

    work page 1956

  30. [30]

    W. Walter. Ordinary differential equations , volume 182 of Graduate Texts in Mathe- matics. Springer-Verlag, New York, 1998. Translated from the sixt h German (1996) edition by Russell Thompson, Readings in Mathematics. 25

    work page 1998