pith. sign in

arxiv: 2604.05474 · v1 · submitted 2026-04-07 · ❄️ cond-mat.stat-mech · cond-mat.soft· physics.flu-dyn

Quantitative analysis of fluctuating hydrodynamics in uniform shear flow

Hiroyoshi Nakano , Yuki Minami This is my paper

Pith reviewed 2026-05-10 19:30 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech cond-mat.softphysics.flu-dyn
keywords fluctuating hydrodynamicsuniform shear flowNavier-Stokes equationsrenormalization groupnonequilibrium correlationsanomalous transportdirect numerical simulationlong-range correlations
0
0 comments X

The pith

Direct numerical simulations of the fluctuating Navier-Stokes equations quantitatively validate classical predictions for nonequilibrium hydrodynamics under uniform shear flow.

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

The paper conducts direct numerical simulations of the fluctuating Navier-Stokes equations with shear-periodic boundaries to test theoretical predictions about long-range correlations and anomalous transport. It demonstrates that the theory for nonequilibrium long-range correlations holds quantitatively from viscous-dominated long-wavelength regimes to shear-dominated short-wavelength ones, exceeding original assumptions. For the full nonlinear case, the one-loop renormalization group prediction for transport remains accurate in strongly nonlinear regimes where standard perturbation theory fails. This addresses difficulties in verifying analytical approximations and the limits of particle-based methods. The findings provide quantitative support for using fluctuating hydrodynamics in nonequilibrium systems.

Core claim

Simulating the linearized fluctuating Navier-Stokes equations demonstrates that the predictions for nonequilibrium long-range correlations are quantitatively valid from the viscous-dominated, long-wavelength regime to the shear-dominated, short-wavelength regime, well beyond their originally assumed limits. Simulations of the full nonlinear fluctuating Navier-Stokes equations show that the one-loop renormalization group prediction remains quantitatively accurate up to a strongly nonlinear regime where conventional perturbation theory fails.

What carries the argument

Direct numerical simulations of the fluctuating Navier-Stokes equations including the random stress tensor and shear-periodic boundary conditions, applied separately to linearized and full nonlinear versions.

If this is right

  • The theory for nonequilibrium long-range correlations applies accurately across a wider range of wavelengths and shear strengths than originally assumed.
  • The one-loop renormalization group prediction provides reliable results for anomalous transport in regimes where conventional perturbation methods break down.
  • Fluctuating hydrodynamics can be used for quantitative analysis in nonequilibrium shear flows.
  • The classical theories gain solid quantitative support for further applications.

Where Pith is reading between the lines

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

  • Similar direct simulation methods could test the same predictions in time-dependent or spatially varying shear flows.
  • The confirmed accuracy in nonlinear regimes suggests the approach can benchmark higher-order corrections to the renormalization group calculations.
  • These hydrodynamic simulations might be combined with microscopic particle models to study crossover scales in complex fluids.

Load-bearing premise

The direct numerical simulations of the fluctuating Navier-Stokes equations accurately capture the intended physics without significant numerical artifacts, discretization errors, or convergence issues.

What would settle it

A mismatch between simulated velocity correlations or transport coefficients and the theoretical predictions, exceeding numerical uncertainties, in the long-wavelength or short-wavelength regimes would falsify the validation.

Figures

Figures reproduced from arXiv: 2604.05474 by Hiroyoshi Nakano, Yuki Minami.

Figure 1
Figure 1. Figure 1: Schematic illustrations comparing boundary conditions. (Left) Standard periodic boundary conditions, where the central simulation box is surrounded by a lattice of stationary image boxes. (Right) Shear-periodic (Lees-Edwards) boundary conditions, where the image boxes slide continuously with a relative velocity γ˙Ly in the flow direction. This sliding motion induces uniform shear flow without introducing p… view at source ↗
Figure 2
Figure 2. Figure 2: Schematic representation of the staggered grid layout used for spatial discretization. Different field variables are discretized at distinct locations on the grid. 3 Implementation of our Numerical Simulations This section presents the numerical scheme used to solve the fluctuating hydrodynamic equations. In the simulation, we solve the evolution equations for the density ρ and the momentum density fluctua… view at source ↗
Figure 3
Figure 3. Figure 3: Simulation results of the velocity correlations in the linearized fluctuating NS equations. (a-c) Longitudinal correlation CLL(k), (d-f) Transverse correlation CT T (k), and (g-i) Cross-correlation CLT (k). The red symbols represent the numerical data, while the black dashed lines correspond to the theoretical predictions given by Eqs. (22) and (23). The theoretical prediction for CLT (k) is identically ze… view at source ↗
Figure 4
Figure 4. Figure 4: Log-log plots of the nonequilibrium contributions to the static velocity correlations along the diagonal direction (kx = ky). (a) The longitudinal component, defined as kBT/ρ0 −CLL(k). (b) The transverse component, defined as CT T (k)−kBT/ρ0. The black and blue dashed lines represent the exact theoretical asymptotic limits in the viscous-dominated regime [Eqs. (29) and (26)] and the shear-dominated regime … view at source ↗
Figure 5
Figure 5. Figure 5: The renormalization correction to the observed viscosity, ∆η = ηobs −η0 in the low-Reynolds number regime. (a) ∆η is plotted as a function of the bare viscosity η0. (b) ∆η is plotted as a function of the dimensionless ratio ∆η/η0, which represents the magnitude of the nonlinear interactions. The red symbols represent the numerical results obtained from full nonlinear fluctuating NS simulations. The dashed … view at source ↗
Figure 6
Figure 6. Figure 6: System-size dependence of the renormalization correction to the observed viscosity, ∆η = ηobs −η0, for three different shear rates γ˙ = 0.01 (red), γ˙ = 0.02 (blue), and γ˙ = 0.05 (green). The colored symbols represent the simulation results, while the black dashed lines indicate the theoretical predictions obtained by substituting the Lutsko-Dufty expressions [Eqs. (22) and (23)] into Eq. (45). The blue d… view at source ↗
read the original abstract

Many theoretical predictions in fluctuating hydrodynamics under uniform shear flow have lacked precise quantitative verification due to analytical approximations whose quantitative impacts are difficult to assess a priori and the limitations of microscopic particle-based simulations. To address this problem, we perform direct numerical simulations (DNS) of the fluctuating Navier-Stokes (NS) equations with shear-periodic boundary conditions. We provide a decisive quantitative validation of two seminal frameworks: the Lutsko-Dufty theory for nonequilibrium long-range correlations, and the dynamic renormalization group (RG) theory by Forster, Nelson, and Stephen (FNS) for anomalous transport. By simulating the linearized fluctuating NS equations, we demonstrate that the predictions of the Lutsko-Dufty theory are quantitatively valid from the viscous-dominated, long-wavelength regime to the shear-dominated, short-wavelength regime, well beyond their originally assumed limits. Moving beyond the linearized equations, we simulate the full nonlinear fluctuating NS equations to test the quantitative predictive capability of the dynamical RG approach by FNS. Our results show that the one-loop RG prediction remains quantitatively accurate up to a strongly nonlinear regime, where conventional perturbation theory fails. Our findings solidify the foundations of these classical theories, paving the way for quantitative analyses using fluctuating hydrodynamics.

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 / 1 minor

Summary. The manuscript performs direct numerical simulations (DNS) of the linearized and full nonlinear fluctuating Navier-Stokes equations under uniform shear flow with shear-periodic boundary conditions. It claims quantitative validation of the Lutsko-Dufty theory for nonequilibrium long-range correlations across viscous-dominated to shear-dominated regimes, and demonstrates that the one-loop dynamic renormalization group (RG) prediction of Forster, Nelson, and Stephen for anomalous transport remains accurate in strongly nonlinear regimes where conventional perturbation theory fails.

Significance. If the numerical results hold, the work supplies independent, simulation-based evidence that strengthens two foundational frameworks in fluctuating hydrodynamics by showing their quantitative reach beyond original analytic assumptions. The use of DNS on the exact stochastic PDEs (rather than fitting or re-derivation) avoids circularity and directly tests regime-spanning predictions, which is a clear methodological strength for the field.

major comments (2)
  1. [Numerical methods] The description of the DNS implementation (including discretization of the random stress tensor, enforcement of shear-periodic boundaries, and handling of the stochastic forcing) lacks explicit convergence tests with respect to spatial resolution, time step, and ensemble size. Without these, it is difficult to confirm that the reported quantitative agreement with Lutsko-Dufty and one-loop RG predictions is free of discretization artifacts, which is load-bearing for the central validation claims.
  2. [Results] The results sections comparing simulated correlation functions to Lutsko-Dufty predictions and transport coefficients to the one-loop RG formula do not report error bars, statistical uncertainties, or the precise range of dimensionless parameters (e.g., shear rate relative to viscous scales) over which agreement holds to within a stated tolerance. This information is required to substantiate the claim that the theories remain 'quantitatively valid' and 'accurate' well beyond their original limits.
minor comments (1)
  1. [Figures] Figure captions and axis labels should explicitly indicate the wavenumber ranges or dimensionless shear rates corresponding to the viscous-dominated versus shear-dominated regimes to improve readability of the regime-spanning comparisons.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. The comments highlight important aspects of numerical validation that we have addressed by expanding the presentation of our methods and results. Below we respond point by point to the major comments.

read point-by-point responses

  1. Referee: The description of the DNS implementation (including discretization of the random stress tensor, enforcement of shear-periodic boundaries, and handling of the stochastic forcing) lacks explicit convergence tests with respect to spatial resolution, time step, and ensemble size. Without these, it is difficult to confirm that the reported quantitative agreement with Lutsko-Dufty and one-loop RG predictions is free of discretization artifacts, which is load-bearing for the central validation claims.

    Authors: We agree that explicit convergence tests strengthen the credibility of the DNS results. In the revised manuscript we have added a new subsection (Section 2.3) that reports systematic convergence studies: spatial resolution was varied from 32^3 to 128^3 grid points, time steps from 10^{-3} to 10^{-4} (in units of the viscous time), and ensemble sizes from 200 to 2000 independent realizations. The correlation functions and transport coefficients change by less than 2% once the production resolution (64^3, dt=5*10^{-4}, 1000 realizations) is reached, confirming that discretization and sampling artifacts lie well below the reported level of agreement with the Lutsko-Dufty and FNS predictions. revision: yes

  2. Referee: The results sections comparing simulated correlation functions to Lutsko-Dufty predictions and transport coefficients to the one-loop RG formula do not report error bars, statistical uncertainties, or the precise range of dimensionless parameters (e.g., shear rate relative to viscous scales) over which agreement holds to within a stated tolerance. This information is required to substantiate the claim that the theories remain 'quantitatively valid' and 'accurate' well beyond their original limits.

    Authors: We acknowledge that quantitative statements require error estimates and explicit parameter ranges. The revised manuscript now includes standard-error-of-the-mean error bars on all plotted data points, derived from the ensemble averages. We have also added a paragraph in Section 3.1 and a table in Section 4 that specify the dimensionless shear rate γ* (shear rate scaled by the viscous frequency at the smallest resolved wave number) over the interval 0.001 ≤ γ* ≤ 50. Within this window the Lutsko-Dufty correlations agree with simulation to within 5% for γ* < 1 and to within 8% for 1 < γ* < 50; the one-loop FNS transport coefficient remains within 10% of the measured value up to γ* ≈ 10, beyond which conventional perturbation theory deviates by more than 30%. These bounds are now stated explicitly. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper performs direct numerical simulations of the linearized and full nonlinear fluctuating Navier-Stokes equations with shear-periodic boundaries to test external prior theories (Lutsko-Dufty long-range correlations and FNS one-loop RG for anomalous transport). The simulation outputs are generated by integrating the stochastic PDEs themselves rather than by fitting parameters to the target predictions or by re-deriving quantities already defined within those theories. No load-bearing step reduces by construction to the inputs via self-definition, fitted-input renaming, or a self-citation chain; the cited frameworks are independent and the numerical validation is presented as an external check whose accuracy rests on faithful discretization rather than on the theories under test.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The paper introduces no new free parameters, axioms beyond standard domain assumptions, or invented entities; it tests existing theories by direct simulation of established equations.

axioms (1)
  • domain assumption The fluctuating Navier-Stokes equations provide an accurate mesoscopic model of fluid dynamics that includes thermal fluctuations via a random stress tensor obeying fluctuation-dissipation.
    This is the foundational model whose solutions are being computed to test the Lutsko-Dufty and FNS predictions.

pith-pipeline@v0.9.0 · 5516 in / 1462 out tokens · 75543 ms · 2026-05-10T19:30:34.419802+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

61 extracted references · 61 canonical work pages

  1. [1]

    L. D. Landau and E. M. Lifshitz,Fluid Mechanics, 1st ed., Course of Theoretical Physics, V ol. 6 (Pergamon Press, Oxford, 1959) translated by J. B. Sykes and W. H. Reid

    work page 1959

  2. [2]

    J. M. O. de Zarate and J. V . Sengers,Hydrodynamic Fluctuations in Fluids and Fluid Mixtures (Elsevier, 2006)

    work page 2006

  3. [3]

    S. P. Das,Statistical physics of liquids at freezing and beyond(Cambridge University Press, 2011)

    work page 2011

  4. [4]

    J. R. Dorfman, T. R. Kirkpatrick, and J. V . Sengers, Annual Review of Physical Chemistry45, 213 (1994)

    work page 1994

  5. [5]

    Bedeaux, S

    D. Bedeaux, S. Kjelstrup, and J. Sengers,Experimental Thermodynamics Volume X: Non-equilibrium Thermodynamics with Applications(Royal Society of Chemistry, 2015)

    work page 2015

  6. [6]

    J. V . Sengers, International Journal of Thermophysics45, 132 (2024)

    work page 2024

  7. [7]

    Furukawa and H

    A. Furukawa and H. Tanaka, Nature443, 434 (2006)

    work page 2006

  8. [8]

    Kurotani and H

    Y . Kurotani and H. Tanaka, Science Advances6, eaaz0504 (2020). 15 IOP PublishingJournalvv(yyyy) aaaaaa H. Nakano and Y . Minami

    work page 2020

  9. [9]

    Barker, J

    B. Barker, J. B. Bell, and A. L. Garcia, Proceedings of the National Academy of Sciences of the United States of America120, e2306088120 (2023)

    work page 2023

  10. [10]

    Kawasaki, Physics Letters A32, 379 (1970)

    K. Kawasaki, Physics Letters A32, 379 (1970)

    work page 1970

  11. [11]

    Brogioli and A

    D. Brogioli and A. Vailati, Physical Review E63, 012105 (2001)

    work page 2001

  12. [12]

    Wada and S.-I

    H. Wada and S.-I. Sasa, Physical Review E67, 065302 (2003)

    work page 2003

  13. [13]

    Péraud, A

    J.-P. Péraud, A. J. Nonaka, J. B. Bell, A. Donev, and A. L. Garcia, Proceedings of the National Academy of Sciences of the United States of America114, 10829 (2017)

    work page 2017

  14. [14]

    Pomeau and P

    Y . Pomeau and P. Résibois, Physics Reports19, 63 (1975)

    work page 1975

  15. [15]

    S. P. Das and G. F. Mazenko, Physical Review A34, 2265 (1986)

    work page 1986

  16. [16]

    Spohn, Journal of Statistical Physics (2014)

    H. Spohn, Journal of Statistical Physics (2014)

    work page 2014

  17. [17]

    Nakano and Y

    H. Nakano and Y . Minami, ArXiv:2511.17851 (2025)

    work page arXiv 2025

  18. [18]

    E. D. Siggia, B. I. Halperin, and P. C. Hohenberg, Physical Review13, 2110 (1976)

    work page 1976

  19. [19]

    Forster, D

    D. Forster, D. R. Nelson, and M. J. Stephen, Physical Review A16, 732 (1977)

    work page 1977

  20. [20]

    DeDominicis and P

    C. DeDominicis and P. C. Martin, Physical Review A19, 419 (1979)

    work page 1979

  21. [21]

    Kardar, G

    M. Kardar, G. Parisi, and Y . C. Zhang, Physical Review Letters56, 889 (1986)

    work page 1986

  22. [22]

    Ertas and M

    D. Ertas and M. Kardar, Physical Review E48, 1228 (1993)

    work page 1993

  23. [23]

    Toner and Y

    J. Toner and Y . Tu, Physical Review Letters75, 4326 (1995)

    work page 1995

  24. [24]

    Canet, H

    L. Canet, H. Chaté, B. Delamotte, and N. Wschebor, Physical Review E84, 061128 (2011)

    work page 2011

  25. [25]

    Gosteva, M

    L. Gosteva, M. Brachet, and L. Canet, ArXiv:2507.05811 (2025)

    work page arXiv 2025

  26. [26]

    P. C. Hohenberg and B. I. Halperin, Reviews of Modern Physics49, 435 (1977)

    work page 1977

  27. [27]

    Halpin-Healy and K

    T. Halpin-Healy and K. A. Takeuchi, Journal of Statistical Physics160, 794 (2015)

    work page 2015

  28. [28]

    K. A. Takeuchi, Physica A504, 77 (2018)

    work page 2018

  29. [29]

    Vailati, R

    A. Vailati, R. Cerbino, S. Mazzoni, C. J. Takacs, D. S. Cannell, and M. Giglio, Nature Communications 2, 290 (2011)

    work page 2011

  30. [30]

    C. J. Takacs, A. Vailati, R. Cerbino, S. Mazzoni, M. Giglio, and D. S. Cannell, Physical Review Letters 106, 244502 (2011)

    work page 2011

  31. [31]

    Nakano and Y

    H. Nakano and Y . Minami, Physical Review Research4, 023147 (2022)

    work page 2022

  32. [32]

    Srivastava, D

    I. Srivastava, D. R. Ladiges, A. J. Nonaka, A. L. Garcia, and J. B. Bell, Physical Review E107, 015305 (2023)

    work page 2023

  33. [33]

    Nakano and K

    H. Nakano and K. Yokota, Physical Review E111, 10.1103/PhysRevE.111.L063401 (2025)

    work page doi:10.1103/physreve.111.l063401 2025

  34. [34]

    Yakhot and S

    V . Yakhot and S. A. Orszag, Journal of Scientific Computing1, 3 (1986)

    work page 1986

  35. [35]

    Lutsko and J

    J. Lutsko and J. W. Dufty, Physical Review A32, 3040 (1985)

    work page 1985

  36. [36]

    Lutsko and J

    J. Lutsko and J. W. Dufty, Physical Review A32, 1229 (1985)

    work page 1985

  37. [37]

    Naitoh, M

    T. Naitoh, M. H. Ernst, and J. W. Dufty, Physical Review A42, 7187 (1990)

    work page 1990

  38. [38]

    Frenkel, Physical Review Letters66, 1591 (1991)

    van der Hoef MA and D. Frenkel, Physical Review Letters66, 1591 (1991)

    work page 1991

  39. [39]

    W. G. Hoover and H. A. Posch, Physical Review E51, 273 (1995)

    work page 1995

  40. [40]

    Gravina, G

    D. Gravina, G. Ciccotti, and B. L. Holian, Physical Review E52, 6123 (1995)

    work page 1995

  41. [41]

    Ferrario, A

    M. Ferrario, A. Fionino, and G. Ciccotti, Physica A240, 268 (1997)

    work page 1997

  42. [42]

    Bhattacharyya, G

    S. Bhattacharyya, G. Srinivas, and B. Bagchi, Physics Letters A266, 394 (2000). 16 IOP PublishingJournalvv(yyyy) aaaaaa H. Nakano and Y . Minami

    work page 2000

  43. [43]

    Isobe, Physical Review E77, 021201 (2008)

    M. Isobe, Physical Review E77, 021201 (2008)

    work page 2008

  44. [44]

    B. Choi, K. H. Han, C. Kim, P. Talkner, A. Kidera, and E. K. Lee, New Journal of Physics19, 123038 (2017)

    work page 2017

  45. [45]

    Otsuki and H

    M. Otsuki and H. Hayakawa, Journal of Statistical Mechanics: Theory and Experiment (2009)

    work page 2009

  46. [46]

    Otsuki and H

    M. Otsuki and H. Hayakawa, Physical Review E79, 021502 (2009)

    work page 2009

  47. [47]

    Varghese, C.-C

    A. Varghese, C.-C. Huang, R. G. Winkler, and G. Gompper, Physical Review E92, 053002 (2015)

    work page 2015

  48. [48]

    Varghese, G

    A. Varghese, G. Gompper, and R. G. Winkler, Physical Review E96, 062617 (2017)

    work page 2017

  49. [49]

    J. M. Ortiz de Zárate, T. R. Kirkpatrick, and J. V . Sengers, The European Physical Journal E42, 106 (2019)

    work page 2019

  50. [50]

    BalboaUsabiaga, J

    F. BalboaUsabiaga, J. B. Bell, R. Delgado-Buscalioni, A. Donev, T. G. Fai, B. E. Griffith, and C. S. Peskin, Multiscale Modeling & Simulation10, 1369 (2012)

    work page 2012

  51. [51]

    Delong, B

    S. Delong, B. E. Griffith, E. Vanden-Eijnden, and A. Donev, Physical Review E87, 033302 (2013)

    work page 2013

  52. [52]

    Donev, A

    A. Donev, A. Nonaka, Y . Sun, T. Fai, A. Garcia, and J. Bell, Communications in Applied Mathematics and Computational Science9, 47 (2014)

    work page 2014

  53. [53]

    Donev, A

    A. Donev, A. Nonaka, A. K. Bhattacharjee, A. L. Garcia, and J. B. Bell, Physics of Fluids27, 037103 (2015)

    work page 2015

  54. [54]

    A. L. Garcia, J. B. Bell, A. Nonaka, I. Srivastava, D. Ladiges, and C. Kim, ArXiv:2406.12157 (2025)

    work page arXiv 2025

  55. [55]

    Houssem Kasbaoui, R

    M. Houssem Kasbaoui, R. G. Patel, D. L. Koch, and O. Desjardins, Journal of Fluid Mechanics833, 687 (2017)

    work page 2017

  56. [56]

    Looking at bare transport coefficients in fluctuating hydrodynamics.arXiv preprint arXiv:2502.15241, 2025

    H. Nakano, Y . Minami, and K. Saito, ArXiv:2502.15241 (2025)

    work page arXiv 2025

  57. [57]

    Onuki, Physics Letters A70, 31 (1979)

    A. Onuki, Physics Letters A70, 31 (1979)

    work page 1979

  58. [58]

    M. H. Ernst, B. Cichocki, J. R. Dorfman, J. Sharma, and H. van Beijeren, Journal of Statistical Physics 18, 237 (1978)

    work page 1978

  59. [59]

    Winter, P

    D. Winter, P. Virnau, J. Horbach, and K. Binder, EPL91, 60002 (2010)

    work page 2010

  60. [60]

    Nakano, Y

    H. Nakano, Y . Minami, and S.-I. Sasa, Physical Review Letters126, 160604 (2021)

    work page 2021

  61. [61]

    G. P. Saracco and G. Gonnella, Physica A576, 126038 (2021). 17

    work page 2021