Lowest order Carleman linearization for steady state fluid flow simulations
Pith reviewed 2026-05-25 04:22 UTC · model grok-4.3
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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)
- [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
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
-
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
-
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
- 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
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
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.
- domain assumption The same truncation order remains accurate for the steady state of the two-dimensional Navier-Stokes equations at moderate Reynolds number.
Reference graph
Works this paper leans on
-
[1]
Schr\"odinger-Navier-Stokes Equation for the Quantum Simulation of Navier-Stokes Flows , author=. 2026 , eprint=
work page 2026
-
[2]
Explicit Error Bounds for Carleman Linearization , author=. 2017 , eprint=
work page 2017
-
[3]
AVS Quantum Science , volume =
Sanavio, Claudio and Succi, Sauro , title =. AVS Quantum Science , volume =. 2024 , month =. doi:10.1116/5.0195549 , url =
-
[4]
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]
Wael Itani and Sauro Succi , title =. Fluids , year =. doi:10.3390/fluids7010024 , url =
-
[6]
Succi, Sauro and Sanavio, Claudio and Cappelli, Luca and Love, Peter , title =. 2026 , month =. doi:10.1209/0295-5075/ae333e , url =
- [7]
- [8]
-
[9]
Carleman, Torsten , title =. Acta Mathematica , volume =. 1932 , doi =
work page 1932
- [10]
-
[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 =
work page 2025
-
[12]
SIAM Journal on Scientific Computing , year =
Hsuan-Cheng Wu and Jingyao Wang and Xiantao Li , title =. SIAM Journal on Scientific Computing , year =
-
[13]
Quantum Carleman linearization efficiency in nonlinear fluid dynamics , author =. Phys. Rev. Res. , volume =. 2025 , month =. doi:10.1103/PhysRevResearch.7.023254 , url =
- [14]
-
[15]
Sachin S. Bharadwaj and Katepalli R. Sreenivasan , title =. S. 2025 , doi =
work page 2025
-
[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]
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]
Antonio David Bastida Zamora and Ljubomir Budinski and Valtteri Lahtinen and Pierre Sagaut , title =. Physical Review E , volume =. 2026 , doi =
work page 2026
-
[19]
Globalizing the Carleman linear embedding method for nonlinear dynamics , author=. 2025 , eprint=
work page 2025
-
[20]
Succi, Sauro and Itani, W. and Sreenivasan, K. and Steijl, R. , title =. 2023 , month =. doi:10.1209/0295-5075/acfdc7 , url =
-
[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]
Potential quantum advantage for simulation of fluid dynamics , author =. Phys. Rev. Res. , volume =. 2025 , month =. doi:10.1103/PhysRevResearch.7.013036 , url =
-
[23]
An end-to-end quantum algorithm for nonlinear fluid dynamics with bounded quantum advantage , author=. 2025 , eprint=
work page 2025
-
[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 =
work page 2022
-
[25]
Sachin S. Bharadwaj and Katepalli R. Sreenivasan , title =. Physical Review Research , volume =. 2025 , doi =
work page 2025
-
[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 =
work page 2026
-
[27]
Communications Physics , volume =
Simulating unsteady flows on a superconducting quantum processor , author =. Communications Physics , volume =. 2024 , doi =
work page 2024
-
[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 =
work page 2026
-
[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]
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 =
work page 2021
-
[31]
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]
F. Tennie and L. Magri , title =. arXiv preprint arXiv:2401.13500 , year =. doi:10.48550/arXiv.2401.13500 , archivePrefix =. 2401.13500 , url =
-
[33]
Xue, Cheng and Wu, Yu-Chun and Guo, Guo-Ping , title =. 2021 , month =. doi:10.1088/1367-2630/ac3eff , url =
-
[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]
Pedro C. S. Costa and Philipp Schleich and Mauro E. S. Morales and Dominic W. Berry , title =. npj Quantum Information , volume =. 2025 , doi =
work page 2025
-
[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]
Annual Review of Fluid Mechanics , volume=
Scale-Invariance and Turbulence Models for Large-Eddy Simulation , author=. Annual Review of Fluid Mechanics , volume=. 2000 , doi=
work page 2000
-
[38]
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 =
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.