pith. sign in

arxiv: 2605.23380 · v1 · pith:7Q7BU3MTnew · submitted 2026-05-22 · 🪐 quant-ph · physics.flu-dyn

Lowest order Carleman linearization for steady state fluid flow simulations

Pith reviewed 2026-05-25 04:22 UTC · model grok-4.3

classification 🪐 quant-ph physics.flu-dyn
keywords Carleman linearizationfluid dynamicssteady stateNavier-Stokestruncationquantum simulation
0
0 comments X

The pith

The second-order truncation of Carleman linearization for the fluid equations recovers the steady-state solution.

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

The paper shows that truncating the Carleman linearization of the Navier-Stokes equations at second order produces not only the early time evolution but also the correct long-time steady state. This property is proved analytically for a forced decaying logistic equation and confirmed numerically for two-dimensional flows at moderate Reynolds number. A reader would care because the result turns a nonlinear problem into a linear one that still reaches the physical equilibrium, which matters for methods that solve linear systems efficiently.

Core claim

It is shown that the lowest (second) order truncation of the Carleman linearization of the fluid equations (C2) recovers not only the initial transient of the time evolution but also its late stage, namely the steady-state solution. This asymptotic property is first proved analytically for the decaying logistic with external forcing and then shown to hold to a significant degree of accuracy also for the fairly more complex case of two-dimensional fluid flows at moderate Reynolds number.

What carries the argument

The second-order Carleman linearization (C2), which replaces the nonlinear fluid equations with a larger linear system by discarding all terms beyond quadratic order in the expansion.

If this is right

  • The steady state satisfies the truncated linear system exactly, so it can be obtained by solving a linear problem rather than integrating the nonlinear equations in time.
  • The property holds analytically for the forced logistic equation and numerically for two-dimensional flows at moderate Reynolds number.
  • Linear-system methods become applicable to steady fluid problems that were previously treated as nonlinear initial-value problems.

Where Pith is reading between the lines

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

  • The same truncation might preserve fixed points in other nonlinear systems whose Carleman expansion has a similar structure.
  • Quantum linear solvers could be applied directly to the C2 system to obtain steady fluid solutions without classical time marching.
  • Higher-dimensional or higher-Reynolds cases can be checked by repeating the numerical comparison performed for the two-dimensional moderate-Re flows.

Load-bearing premise

The error from discarding terms beyond second order does not move or destabilize the fixed point that the full fluid equations approach at late times.

What would settle it

Solve the steady Navier-Stokes equations and the C2-truncated linear system for the same two-dimensional domain, boundary conditions, and moderate Reynolds number; the velocity fields must differ by more than numerical tolerance if the claim is false.

Figures

Figures reproduced from arXiv: 2605.23380 by Luca Cappelli, Sauro Succi.

Figure 1
Figure 1. Figure 1: FIG. 1: (a) Solution of the forced logistic equation ( [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Evolution of the velocity field in three different spatial locations: [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Evolution of the velocity field in three spatial locations of the previous figures. The initial condition in Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Amplitude and streamlines of the velocity field at [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Evolution of the velocity field along the [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Amplitude and streamlines of the velocity field obtained by evolving the initial condition in Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
read the original abstract

It is shown that the lowest (second) order truncation of the Carleman linearization of the fluid equations (C2) recovers not only the initial transient of the time evolution but also its late stage, namely the steady-state solution. This asymptotic property is first proved analytically for the decaying logistic with external forcing and then shown to hold to a significant degree of accuracy also for the fairly more complex case of two-dimensional fluid flows at moderate Reynolds number. This time-asymptotic property opens interesting prospects for the simulation of steady-state solutions of the fluid equations on quantum computers.

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 claims that the lowest (second) order truncation of the Carleman linearization of the fluid equations (denoted C2) recovers not only the initial transient but also the late-time steady-state solution. An analytic proof is provided for the forced logistic equation; for two-dimensional Navier-Stokes flows at moderate Reynolds number the property is asserted to hold to a significant degree of accuracy on the basis of numerical experiments. The result is presented as enabling quantum-computer simulation of steady fluid states.

Significance. If the fixed-point preservation property extends rigorously beyond the logistic case, the work would supply a concrete route for quantum linear-system solvers to reach steady fluid solutions without requiring the full nonlinear Carleman embedding or long-time integration. The analytic demonstration for the logistic equation is a clear strength; the numerical fluid evidence, however, remains preliminary.

major comments (2)
  1. [fluid-flow numerical section (implicit in abstract and introduction)] The central claim for the Navier-Stokes equations—that C2 leaves the steady-state fixed point unchanged (up to discretization error)—rests on numerical observation rather than derivation. No section demonstrates that the quadratic truncation terms vanish identically when the time derivative is set to zero, nor that the resulting algebraic system for the fixed point coincides with the original steady NS equations.
  2. [numerical results for 2D flows] The fluid-flow evidence is described only as holding 'to a significant degree of accuracy' with no error metrics, mesh details, Reynolds-number values, or comparison baselines supplied. This is insufficient to confirm that the truncation error does not alter the location or stability of the asymptotic fixed point.
minor comments (1)
  1. [abstract] The abstract and introduction should explicitly state the Reynolds-number range and the quantitative tolerance used to declare agreement with the reference steady solution.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the detailed and constructive report. We address the two major comments point by point below, indicating planned revisions where appropriate. The manuscript's central claim for the Navier-Stokes equations is presented as a numerical observation rather than a general theorem, consistent with the analytic proof being limited to the logistic case.

read point-by-point responses
  1. Referee: The central claim for the Navier-Stokes equations—that C2 leaves the steady-state fixed point unchanged (up to discretization error)—rests on numerical observation rather than derivation. No section demonstrates that the quadratic truncation terms vanish identically when the time derivative is set to zero, nor that the resulting algebraic system for the fixed point coincides with the original steady NS equations.

    Authors: We agree that no general derivation is provided for the Navier-Stokes equations; the analytic fixed-point preservation is proved only for the forced logistic equation. For the fluid case the manuscript relies on numerical evidence that the C2 truncation recovers the steady state to a significant degree of accuracy. We will revise the text to state this distinction more explicitly in the abstract, introduction, and conclusions, and we will add a short paragraph discussing the structural reason why quadratic truncation can be expected to preserve the fixed point (the time-derivative term vanishes and the remaining algebraic system is closed at second order), while acknowledging that a rigorous proof for NS remains open. revision: partial

  2. Referee: The fluid-flow evidence is described only as holding 'to a significant degree of accuracy' with no error metrics, mesh details, Reynolds-number values, or comparison baselines supplied. This is insufficient to confirm that the truncation error does not alter the location or stability of the asymptotic fixed point.

    Authors: We acknowledge the lack of quantitative details in the numerical section. In the revised manuscript we will expand the fluid-flow results to report: (i) specific Reynolds numbers and mesh resolutions used, (ii) quantitative error norms (L2 and L-infinity) between the C2 steady-state solution and the reference steady Navier-Stokes solution obtained by direct nonlinear iteration, (iii) comparison against a standard baseline solver, and (iv) a brief assessment of stability of the recovered fixed point. These additions will allow readers to evaluate the truncation error more precisely. revision: yes

standing simulated objections not resolved
  • A rigorous analytic demonstration that the second-order Carleman truncation preserves the steady-state fixed point for the full Navier-Stokes equations (beyond the logistic case) is not currently available and would require substantial additional theoretical development.

Circularity Check

0 steps flagged

No circularity; analytic logistic proof independent of fluid numerics

full rationale

The paper derives the time-asymptotic fixed-point property analytically for the forced logistic equation and reports separate numerical checks for 2D Navier-Stokes at moderate Re. No equation or claim reduces the fluid result to the logistic result by construction, and no self-citation is invoked as load-bearing justification for the fluid case. The fluid evidence is presented explicitly as numerical agreement rather than an identity or forced prediction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the analytic property shown for the logistic equation and on the assumption that the same truncation order suffices for the Navier-Stokes steady state; no free parameters or invented entities are mentioned.

axioms (2)
  • domain assumption The second-order truncation of the Carleman embedding preserves the location of the asymptotic fixed point for the forced logistic equation.
    This is the analytically proved base case on which the fluid claim is built.
  • domain assumption The same truncation order remains accurate for the steady state of the two-dimensional Navier-Stokes equations at moderate Reynolds number.
    This is the step that extends the logistic result to fluids and is supported only numerically in the abstract.

pith-pipeline@v0.9.0 · 5614 in / 1157 out tokens · 23919 ms · 2026-05-25T04:22:34.960756+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

38 extracted references · 38 canonical work pages

  1. [1]

    2026 , eprint=

    Schr\"odinger-Navier-Stokes Equation for the Quantum Simulation of Navier-Stokes Flows , author=. 2026 , eprint=

  2. [2]

    2017 , eprint=

    Explicit Error Bounds for Carleman Linearization , author=. 2017 , eprint=

  3. [3]

    AVS Quantum Science , volume =

    Sanavio, Claudio and Succi, Sauro , title =. AVS Quantum Science , volume =. 2024 , month =. doi:10.1116/5.0195549 , url =

  4. [4]

    and Scatamacchia, R

    Sanavio, C. and Scatamacchia, R. and de Falco, C. and Succi, S. , title =. Physics of Fluids , volume =. 2024 , month =. doi:10.1063/5.0204955 , url =

  5. [5]

    Fluids , year =

    Wael Itani and Sauro Succi , title =. Fluids , year =. doi:10.3390/fluids7010024 , url =

  6. [6]

    2026 , month =

    Succi, Sauro and Sanavio, Claudio and Cappelli, Luca and Love, Peter , title =. 2026 , month =. doi:10.1209/0295-5075/ae333e , url =

  7. [7]

    Dietrich and D

    K. Dietrich and D. Vautherin , title =. Journal de Physique , year =

  8. [8]

    1995 , doi =

    Uriel Frisch , title =. 1995 , doi =

  9. [9]

    Acta Mathematica , volume =

    Carleman, Torsten , title =. Acta Mathematica , volume =. 1932 , doi =

  10. [10]

    2006 , isbn =

    Pierre Sagaut , title =. 2006 , isbn =

  11. [11]

    Discrete and Continuous Dynamical Systems B , volume =

    Arash Amini and Cong Zheng and Qiyu Sun and Nader Motee , title =. Discrete and Continuous Dynamical Systems B , volume =. 2025 , doi =

  12. [12]

    SIAM Journal on Scientific Computing , year =

    Hsuan-Cheng Wu and Jingyao Wang and Xiantao Li , title =. SIAM Journal on Scientific Computing , year =

  13. [13]

    Quantum Carleman linearization efficiency in nonlinear fluid dynamics , author =. Phys. Rev. Res. , volume =. 2025 , month =. doi:10.1103/PhysRevResearch.7.023254 , url =

  14. [14]

    1991 , doi =

    Krzysztof Kowalski and Wojciech-Hubert Steeb , title =. 1991 , doi =

  15. [15]

    Bharadwaj and Katepalli R

    Sachin S. Bharadwaj and Katepalli R. Sreenivasan , title =. S. 2025 , doi =

  16. [16]

    npj Quantum Information , year =

    Wang, Boyuan and Meng, Zhaoyuan and Zhao, Yaomin and Yang, Yue , title =. npj Quantum Information , year =. doi:10.1038/s41534-025-01142-6 , publisher =

  17. [17]

    Nature Reviews Physics , year =

    Tennie, Felix and Laizet, Sylvain and Lloyd, Seth and Magri, Luca , title =. Nature Reviews Physics , year =. doi:10.1038/s42254-024-00799-w , url =

  18. [18]

    Physical Review E , volume =

    Antonio David Bastida Zamora and Ljubomir Budinski and Valtteri Lahtinen and Pierre Sagaut , title =. Physical Review E , volume =. 2026 , doi =

  19. [19]

    2025 , eprint=

    Globalizing the Carleman linear embedding method for nonlinear dynamics , author=. 2025 , eprint=

  20. [20]

    and Sreenivasan, K

    Succi, Sauro and Itani, W. and Sreenivasan, K. and Steijl, R. , title =. 2023 , month =. doi:10.1209/0295-5075/acfdc7 , url =

  21. [21]

    and Ralli, Alexis and Love, Peter and Succi, Sauro , title =

    Sanavio, Claudio and Simon, William A. and Ralli, Alexis and Love, Peter and Succi, Sauro , title =. Physics of Fluids , volume =. 2025 , month =. doi:10.1063/5.0254588 , url =

  22. [22]

    Potential quantum advantage for simulation of fluid dynamics , author =. Phys. Rev. Res. , volume =. 2025 , month =. doi:10.1103/PhysRevResearch.7.013036 , url =

  23. [23]

    2025 , eprint=

    An end-to-end quantum algorithm for nonlinear fluid dynamics with bounded quantum advantage , author=. 2025 , eprint=

  24. [24]

    and Debbasch, Fabrice , title =

    Zylberman, Julien and Di Molfetta, Giuseppe and Brachet, Marc and Loureiro, Nuno F. and Debbasch, Fabrice , title =. Physical Review A , volume =. 2022 , doi =

  25. [25]

    Bharadwaj and Katepalli R

    Sachin S. Bharadwaj and Katepalli R. Sreenivasan , title =. Physical Review Research , volume =. 2025 , doi =

  26. [26]

    Surrogate Quantum Circuit Design for the Lattice Boltzmann Collision Operator , journal =

    L. Surrogate Quantum Circuit Design for the Lattice Boltzmann Collision Operator , journal =. 2026 , doi =

  27. [27]

    Communications Physics , volume =

    Simulating unsteady flows on a superconducting quantum processor , author =. Communications Physics , volume =. 2024 , doi =

  28. [28]

    Nature Communications , volume =

    Ziteng Wang and Jiarun Zhong and Ke Wang and Zitian Zhu and Zehang Bao and Chenjia Zhu and Wenwen Zhao and Yaomin Zhao and Yue Yang and Chao Song and Shiying Xiong , title =. Nature Communications , volume =. 2026 , doi =

  29. [29]

    Bharadwaj and Philipp Pfeffer and Katepalli R

    Julia Ingelmann and Sachin S. Bharadwaj and Philipp Pfeffer and Katepalli R. Sreenivasan and J\"org Schumacher , keywords =. Two quantum algorithms for solving the one-dimensional advection–diffusion equation , journal =. 2024 , issn =. doi:https://doi.org/10.1016/j.compfluid.2024.106369 , url =

  30. [30]

    Efficient quantum algorithm for dissipative nonlinear differential equations , journal =

    Jin-Peng Liu and Hermund. Efficient quantum algorithm for dissipative nonlinear differential equations , journal =. 2021 , doi =

  31. [31]

    Katz and G

    J. Katz and G. Muraleedharan and A. Alase , title =. arXiv preprint arXiv:2512.06488 , year =. doi:10.48550/arXiv.2512.06488 , archivePrefix =. 2512.06488 , primaryClass =

  32. [32]

    Tennie and L

    F. Tennie and L. Magri , title =. arXiv preprint arXiv:2401.13500 , year =. doi:10.48550/arXiv.2401.13500 , archivePrefix =. 2401.13500 , url =

  33. [33]

    2021 , month =

    Xue, Cheng and Wu, Yu-Chun and Guo, Guo-Ping , title =. 2021 , month =. doi:10.1088/1367-2630/ac3eff , url =

  34. [34]

    Quantum algorithm for solving a quadratic nonlinear system of equations , author =. Phys. Rev. A , volume =. 2022 , month =. doi:10.1103/PhysRevA.106.032427 , url =

  35. [35]

    Pedro C. S. Costa and Philipp Schleich and Mauro E. S. Morales and Dominic W. Berry , title =. npj Quantum Information , volume =. 2025 , doi =

  36. [36]

    Diffusion and hydrodynamic dispersion with the lattice Boltzmann method , author =. Phys. Rev. A , volume =. 1992 , month =. doi:10.1103/PhysRevA.45.5771 , url =

  37. [37]

    Annual Review of Fluid Mechanics , volume=

    Scale-Invariance and Turbulence Models for Large-Eddy Simulation , author=. Annual Review of Fluid Mechanics , volume=. 2000 , doi=

  38. [38]

    2019 , issn =

    A review of current progress in multiscale simulations for fluid flow and heat transfer problems: The frameworks, coupling techniques and future perspectives , journal =. 2019 , issn =. doi:https://doi.org/10.1016/j.ijheatmasstransfer.2019.04.004 , url =