pith. sign in

arxiv: 2606.31784 · v1 · pith:72EIR65Unew · submitted 2026-06-30 · 🧮 math-ph · math.MP

Kinetic derivation of thermal viscous models for nematic liquid crystal dynamics

P. E. Farrell , J. M\'alek , O. Sou\v{c}ek , U. Zerbinati This is my paper

Pith reviewed 2026-07-01 02:46 UTC · model grok-4.3

classification 🧮 math-ph math.MP
keywords nematic liquid crystalskinetic theoryBGK collision operatorChapman-Enskog expansionconstitutive equationsviscous modelsthermodynamic theorycouple stress
0
0 comments X

The pith

A kinetic theory with BGK collisions and time scale separation produces constitutive equations for viscous and thermal effects in nematic liquid crystals.

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

The paper starts from a kinetic model of ordered fluids that includes mean-field alignment and a BGK collision operator. It assumes orientational relaxation happens faster than momentum relaxation. Chapman-Enskog expansions then yield the free energy and the structure of entropy production. A constrained maximization procedure determines the explicit forms of the Cauchy stress, couple stress, and the energy and entropy fluxes. This creates both compressible and incompressible viscous thermal models that extend earlier inviscid kinetic theories.

Core claim

Starting from a kinetic theory of ordered fluids with a BGK-type collision operator and Vlasov potential, and relying on separation of time scales with faster orientational relaxation, the zeroth and first order Chapman-Enskog expansions establish the balance equations for mass, momentum, energy and entropy together with a constitutive equation for the Helmholtz free energy and the associated structural form of the entropy production rate. Additional information from the kinetic description determines a constitutive relation for the entropy production rate itself. Application of the constrained maximisation procedure of Rajagopal and Srinivasa then yields constitutive equations for the Cauch

What carries the argument

The constrained maximisation procedure of Rajagopal and Srinivasa applied to the entropy production rate identified via Chapman-Enskog expansion of the BGK kinetic model.

If this is right

  • Viscous dissipation, thermal conduction and spin diffusion are incorporated into the macroscopic balance laws.
  • The resulting models satisfy the second law through the entropy production obtained from the kinetic description.
  • Both compressible and incompressible formulations are available from the same derivation.
  • The approach supplies a systematic kinetic basis for dissipative effects in nematic dynamics.

Where Pith is reading between the lines

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

  • The explicit constitutive relations could be matched term-by-term to coefficients appearing in established phenomenological models of nematic flow.
  • Numerical solutions of the derived equations in simple shear or Poiseuille geometries would allow direct comparison with molecular dynamics trajectories.
  • Replacing the BGK operator with other collision models would test how sensitive the final stress expressions are to the choice of kinetic approximation.

Load-bearing premise

Orientational relaxation occurs on a much faster time scale than translational momentum relaxation.

What would settle it

Direct measurement of the Cauchy stress and couple stress components in a sheared nematic sample under a temperature gradient, compared to the explicit forms obtained from the maximization procedure.

Figures

Figures reproduced from arXiv: 2606.31784 by J. M\'alek, O. Sou\v{c}ek, P. E. Farrell, U. Zerbinati.

Figure 1
Figure 1. Figure 1: Graphical representation of the self-consistent equation for the order parameter s in the case of a Kuramoto-type mean-field potential. The dashed lines represent the function I(s) = coth  s kBθ  − kBθ s while the solid line represents the function s. Each plot corresponds to a different value of the temperature θ; from left to right, the temperature θ decreases. (i) The isotropic distribution λiso(ν) = … view at source ↗
read the original abstract

We develop a macroscopic thermodynamic theory of nematic liquid crystals starting from a kinetic theory of ordered fluids with a collision operator of Bhatnagar-Gross-Krook (BGK) type. The kinetic description incorporates mean-field alignment interactions through a Vlasov potential and relies on a separation of time scales, with orientational relaxation occurring on a faster time scale than translational momentum relaxation. At the continuum level, we establish the balance equations for mass, linear and angular momentum, energy, and entropy. Using the zeroth and first order Chapman-Enskog expansions, we derive a constitutive equation for the Helmholtz free energy and identify the associated structural form of the entropy production rate. We then exploit additional information from the kinetic description to determine a constitutive relation for the entropy production rate itself. Finally, by applying the constrained maximisation procedure of Rajagopal and Srinivasa, we obtain constitutive equations for the Cauchy stress and couple-stress tensors, as well as for the energy and entropy fluxes. In this way we generalise the recent inviscid kinetic theory of Farrell, Russo, and Zerbinati to account for viscous, thermal, and spin-diffusive effects, using the simplest BGK-type approximation of the collision operator. Both compressible and incompressible variants of the theory are presented.

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 paper starts from a kinetic model of ordered fluids with BGK-type collisions and Vlasov mean-field alignment, assumes a separation of orientational and translational relaxation times, performs zeroth- and first-order Chapman-Enskog expansions to obtain balance laws and a structural form for the entropy production, uses kinetic information to fix the entropy-production constitutive relation, and then applies the Rajagopal-Srinivasa constrained-maximization procedure to derive constitutive equations for the Cauchy stress, couple stress, energy flux, and entropy flux. Both compressible and incompressible versions are presented, generalizing an earlier inviscid kinetic theory to include viscous, thermal, and spin-diffusive effects.

Significance. If the derivations are free of post-hoc choices and the time-scale separation is consistently implemented, the work supplies an explicit kinetic-to-macroscopic route to thermodynamically consistent viscous thermal nematic models. The use of an explicit BGK operator and the Rajagopal-Srinivasa step are strengths that could make the resulting constitutive relations falsifiable and reproducible once the expansions are written out.

major comments (2)
  1. [§3] §3 (Chapman-Enskog procedure): the claim that the entropy-production rate is fully determined by kinetic information appears to rest on the specific form chosen for the BGK operator and the Vlasov potential; it is not shown whether this determination remains independent of the relaxation-time ratio once the first-order correction is inserted into the entropy balance.
  2. [§4] §4 (Rajagopal-Srinivasa maximization): the constrained-maximization step yields constitutive relations for stress and fluxes, but the manuscript does not verify that the resulting dissipation inequality is satisfied identically for the derived expressions rather than only up to higher-order terms neglected in the expansion.
minor comments (2)
  1. Notation for the orientation distribution function and the Vlasov potential should be introduced once and used consistently; several symbols are redefined between the kinetic and macroscopic sections.
  2. The incompressible limit is stated to follow by a standard projection, but the precise constraint on the velocity field and the resulting simplification of the couple-stress equation are not written explicitly.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed reading and the two substantive comments on the Chapman-Enskog step and the Rajagopal-Srinivasa procedure. We respond to each below and indicate the revisions we are prepared to make.

read point-by-point responses

  1. Referee: §3 (Chapman-Enskog procedure): the claim that the entropy-production rate is fully determined by kinetic information appears to rest on the specific form chosen for the BGK operator and the Vlasov potential; it is not shown whether this determination remains independent of the relaxation-time ratio once the first-order correction is inserted into the entropy balance.

    Authors: The entropy-production expression is obtained directly from the chosen BGK collision operator and Vlasov potential after the first-order correction is substituted into the entropy balance that follows from the kinetic equation. The time-scale separation is an explicit modeling assumption of the paper; the ratio of relaxation times therefore enters the first-order terms by construction. We do not claim, nor does the derivation show, that the same entropy-production form would be recovered for an arbitrary collision operator or for arbitrary ratios. We will revise §3 to state this model dependence explicitly and to note that independence from the ratio would require a more general collision model outside the present scope. revision: partial

  2. Referee: §4 (Rajagopal-Srinivasa maximization): the constrained-maximization step yields constitutive relations for stress and fluxes, but the manuscript does not verify that the resulting dissipation inequality is satisfied identically for the derived expressions rather than only up to higher-order terms neglected in the expansion.

    Authors: The Rajagopal-Srinivasa procedure is applied to the entropy-production functional obtained from the kinetic expansion; the resulting constitutive relations therefore satisfy the dissipation inequality by construction at the order retained. Because the expansion is truncated at first order, the inequality holds identically only up to the neglected higher-order terms. We will add, in §4, an explicit substitution of the derived stress, couple-stress and flux expressions back into the dissipation inequality to confirm that it is satisfied identically within the first-order approximation. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation grounded in external kinetic model and Rajagopal-Srinivasa procedure

full rationale

The paper begins from an independent kinetic description (BGK collision operator, Vlasov mean-field, explicit time-scale separation) and applies the external constrained-maximisation framework of Rajagopal and Srinivasa. Chapman-Enskog expansions are used to obtain constitutive relations from the kinetic starting point rather than fitting parameters to the target macroscopic quantities. No load-bearing step reduces by construction to a self-citation, a fitted input renamed as prediction, or an ansatz smuggled via prior work by the same authors. The derivation chain remains self-contained against the stated external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Review limited to abstract; the derivation rests on standard kinetic-theory assumptions plus two explicit modeling choices stated in the abstract.

axioms (2)
  • domain assumption Separation of time scales with orientational relaxation faster than translational momentum relaxation
    Explicitly invoked as the basis for the kinetic description and Chapman-Enskog ordering.
  • domain assumption BGK-type collision operator as the simplest approximation
    Chosen to close the kinetic model while retaining mean-field Vlasov alignment.

pith-pipeline@v0.9.1-grok · 5769 in / 1346 out tokens · 61332 ms · 2026-07-01T02:46:19.976795+00:00 · methodology

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

25 extracted references · 20 canonical work pages

  1. [1]

    2024 Kinetic derivation of an inviscid compressible Leslie– Ericksen equation for rarified calamitic gases.Multiscale Modeling & Simulation22, 1585–1607

    Farrell PE, Russo G, Zerbinati U. 2024 Kinetic derivation of an inviscid compressible Leslie– Ericksen equation for rarified calamitic gases.Multiscale Modeling & Simulation22, 1585–1607. (10.1137/24M1630529) 42royalsocietypublishing.org/journal/rspa Proc R Soc A 0000000

    work page doi:10.1137/24m1630529 2024

  2. [2]

    1956 Kinetic theory of nonspherical molecules.The Journal of Chemical Physics24, 225–241

    Curtiss CF. 1956 Kinetic theory of nonspherical molecules.The Journal of Chemical Physics24, 225–241. (10.1063/1.1742459)

    work page doi:10.1063/1.1742459 1956

  3. [3]

    1957 Kinetic theory of nonspherical molecules

    Curtiss CF, Muckenfuss C. 1957 Kinetic theory of nonspherical molecules. II.The Journal of Chemical Physics26, 1619–1636. (10.1063/1.1743595)

    work page doi:10.1063/1.1743595 1957

  4. [4]

    1958 Kinetic theory of nonspherical molecules

    Muckenfuss C, Curtiss CF. 1958 Kinetic theory of nonspherical molecules. III.The Journal of Chemical Physics29, 1257–1272. (10.1063/1.1744708)

    work page doi:10.1063/1.1744708 1958

  5. [5]

    1959 Kinetic theory of nonspherical molecules

    Livingston PM, Curtiss CF. 1959 Kinetic theory of nonspherical molecules. IV. Angular momentum transport coefficient.The Journal of Chemical Physics31, 1643–1645. (10.1063/1.1730668)

    work page doi:10.1063/1.1730668 1959

  6. [6]

    1963 Kinetic theory of nonspherical molecules

    Curtiss CF, Dahler JS. 1963 Kinetic theory of nonspherical molecules. V.The Journal of Chemical Physics38, 2352–2363. (10.1063/1.1733510)

    work page doi:10.1063/1.1733510 1963

  7. [7]

    2025 A kinetic theory approach to ordered fluids

    Carrillo JA, Farrell PE, Medaglia A, Zerbinati U. 2025 A kinetic theory approach to ordered fluids. arXiv:2508.10744v2 (10.48550/arXiv.2508.10744)

    work page doi:10.48550/arxiv.2508.10744 2025

  8. [8]

    2026 A kinetic interpretation of thermomechanical restrictions of continua.International Journal of Engineering Science225, 104557

    Farrell PE, Málek J, Souˇ cek O, Zerbinati U. 2026 A kinetic interpretation of thermomechanical restrictions of continua.International Journal of Engineering Science225, 104557. (10.1016/j.ijengsci.2026.104557)

    work page doi:10.1016/j.ijengsci.2026.104557 2026

  9. [9]

    1980Fundamentals of Maxwell’s kinetic theory of a simple monatomic gas: Treated as a branch of rational mechanicsvol

    Truesdell C, Muncaster RG. 1980Fundamentals of Maxwell’s kinetic theory of a simple monatomic gas: Treated as a branch of rational mechanicsvol. 83,Pure and Applied Mathematics. New York: Academic Press

  10. [10]

    2004 On thermomechanical restrictions of continua.Proceedings of the Royal Society of London

    Rajagopal KR, Srinivasa AR. 2004 On thermomechanical restrictions of continua.Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences460, 631–651. (10.1098/rspa.2002.1111)

    work page doi:10.1098/rspa.2002.1111 2004

  11. [11]

    1949 The effects of shape on the interaction of colloidal particles.Annals of the New York Academy of Sciences51, 627–659

    Onsager L. 1949 The effects of shape on the interaction of colloidal particles.Annals of the New York Academy of Sciences51, 627–659. (10.1111/j.1749-6632.1949.tb27296.x)

    work page doi:10.1111/j.1749-6632.1949.tb27296.x 1949

  12. [12]

    2017 Onsager’s missing steps retraced.Journal of Physics: Condensed Matter29, 475102

    Palffy-Muhoray P , Virga EG, Zheng X. 2017 Onsager’s missing steps retraced.Journal of Physics: Condensed Matter29, 475102. (10.1088/1361-648X/aa942b)

    work page doi:10.1088/1361-648x/aa942b 2017

  13. [13]

    1989Continua with microstructurevol

    Capriz G. 1989Continua with microstructurevol. 35,Springer Tracts in Natural Philosophy. Springer

  14. [14]

    2011 On Landau damping.Acta Mathematica207, 29–201

    Mouhot C, Villani C. 2011 On Landau damping.Acta Mathematica207, 29–201. (10.1007/s11511-011-0068-9)

    work page doi:10.1007/s11511-011-0068-9 2011

  15. [15]

    1946 On the Vibrations of the Electronic Plasma.Journal of Physics (USSR)10, 25–34

    Landau LD. 1946 On the Vibrations of the Electronic Plasma.Journal of Physics (USSR)10, 25–34. Translated from Zh. Eksp. Teor. Fiz., 16 (1946), 574–586

    1946

  16. [16]

    2016 Asymptotics of mean-field O(N) models.Journal of Statistical Physics165, 1114–1140

    Kirkpatrick K, Nawaz T. 2016 Asymptotics of mean-field O(N) models.Journal of Statistical Physics165, 1114–1140. (10.1007/s10955-016-1667-9)

    work page doi:10.1007/s10955-016-1667-9 2016

  17. [17]

    2019Variational theories for liquid crystalsvol

    Virga EG. 2019Variational theories for liquid crystalsvol. 8,Applied Mathematics and Mathematical Computation. Taylor & Francis Group

  18. [18]

    1960 Anisotropic fluids.Archive for Rational Mechanics and Analysis4, 231–237

    Ericksen JL. 1960 Anisotropic fluids.Archive for Rational Mechanics and Analysis4, 231–237. (10.1007/BF00281389)

    work page doi:10.1007/bf00281389 1960

  19. [19]

    1961 Conservation laws for liquid crystals.Transactions of the Society of Rheology5, 23–34

    Ericksen JL. 1961 Conservation laws for liquid crystals.Transactions of the Society of Rheology5, 23–34. (10.1122/1.548883)

    work page doi:10.1122/1.548883 1961

  20. [20]

    1962 Hydrostatic theory of liquid crystals.Archive for Rational Mechanics and Analysis9, 371–378

    Ericksen JL. 1962 Hydrostatic theory of liquid crystals.Archive for Rational Mechanics and Analysis9, 371–378. (10.1007/bf00253358)

    work page doi:10.1007/bf00253358 1962

  21. [21]

    1968 Some constitutive equations for liquid crystals.Archive for Rational Mechanics and Analysis28, 265–283

    Leslie FM. 1968 Some constitutive equations for liquid crystals.Archive for Rational Mechanics and Analysis28, 265–283. (10.1007/bf00251810)

    work page doi:10.1007/bf00251810 1968

  22. [22]

    1993The Physics of Liquid Crystals

    de Gennes PG, Prost J. 1993The Physics of Liquid Crystals. Number 83 in International Series of Monographs on Physics. Oxford University Press 2nd edition

  23. [23]

    1974 Physics of liquid crystals.Reviews of Modern Physics46, 617–704

    Stephen MJ, Straley JP . 1974 Physics of liquid crystals.Reviews of Modern Physics46, 617–704. (10.1103/RevModPhys.46.617)

    work page doi:10.1103/revmodphys.46.617 1974

  24. [24]

    2018 Derivation of Equations for Continuum Mechanics and Thermodynamics of Fluids

    Málek J, Pr ˚ uša V . 2018 Derivation of Equations for Continuum Mechanics and Thermodynamics of Fluids. In Giga Y, Novotný A, editors,Handbook of Mathematical Analysis in Mechanics of Viscous Fluids, pp. 3–72. Cham: Springer International Publishing. (10.1007/978-3-319-13344-7_1)

    work page doi:10.1007/978-3-319-13344-7_1 2018

  25. [25]

    1981 On the H-theorem for polyatomic gases.Journal of Statistical Physics26, 795–801

    Cercignani C, Lampis M. 1981 On the H-theorem for polyatomic gases.Journal of Statistical Physics26, 795–801. (10.1007/BF01010940)

    work page doi:10.1007/bf01010940 1981