arxiv: 2604.09934 · v2 · submitted 2026-04-10 · 🧮 math.AP

Recognition: no theorem link

A Carleman contraction method for inverse initial data recovery in the Navier-Stokes equations with unknown body force

Phuong M.Nguyen , Loc H.Nguyen

Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:39 UTC · model grok-4.3

classification 🧮 math.AP

keywords inverse initial data problemNavier-Stokes equationsCarleman estimatescontraction mappingLegendre polynomial reductionboundary observationsunknown body force

0 comments

The pith

A Carleman-weighted contraction mapping recovers initial velocity and pressure for Navier-Stokes equations from boundary observations even when the body force is unknown.

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

The paper establishes a method to recover the initial velocity and pressure in the incompressible Navier-Stokes system using only lateral boundary observations, without any knowledge of the time-independent body force. Time differentiation of the momentum equation removes the unknown force, after which a Legendre polynomial-exponential reduction converts the problem into a coupled system of elliptic equations. A map defined on an admissible set and equipped with a Carleman-weighted norm is shown to be contractive, so its fixed point approximates the reduced model and the associated Picard iteration converges globally. Numerical experiments confirm that this procedure reconstructs the initial data accurately from synthetic boundary measurements.

Core claim

Differentiating the Navier-Stokes momentum equation with respect to time eliminates the unknown body force. A Legendre polynomial-exponential time-dimensional reduction then produces a coupled elliptic system for the expansion coefficients. On a suitable admissible set equipped with a Carleman-weighted norm, a map is constructed that is contractive; its fixed point yields an approximate solution of the reduced model, and the contraction property guarantees global convergence of the Picard iteration.

What carries the argument

The contractive mapping on an admissible set equipped with a Carleman-weighted norm for the elliptic system obtained after time differentiation and Legendre polynomial-exponential reduction.

Load-bearing premise

The admissible set must be non-empty and the Carleman-weighted norm must make the constructed map contractive on the reduced elliptic system.

What would settle it

For a manufactured solution of the Navier-Stokes equations with known initial data and boundary traces, the Picard iteration fails to converge or produces a reconstruction that violates the original equations within the expected error tolerance.

Figures

Figures reproduced from arXiv: 2604.09934 by Loc H.Nguyen, Phuong M.Nguyen.

**Figure 2.** Figure 2: Test 1. Convergence history of the Picard iteration. The left image shows the relative changes [PITH_FULL_IMAGE:figures/full_fig_p023_2.png] view at source ↗

**Figure 3.** Figure 3: Test 2. First row: true initial velocity components and pressure. Second row: reconstructed initial [PITH_FULL_IMAGE:figures/full_fig_p024_3.png] view at source ↗

**Figure 4.** Figure 4: Test 2. Convergence history of the Picard iteration. The left image shows the relative changes [PITH_FULL_IMAGE:figures/full_fig_p025_4.png] view at source ↗

**Figure 5.** Figure 5: Test 3. First row: true initial velocity components and pressure. Second row: reconstructed initial [PITH_FULL_IMAGE:figures/full_fig_p026_5.png] view at source ↗

**Figure 6.** Figure 6: Test 3. Convergence history of the Picard iteration. The left image shows the relative changes [PITH_FULL_IMAGE:figures/full_fig_p027_6.png] view at source ↗

read the original abstract

We solve an inverse initial data problem for the incompressible Navier-Stokes system. The objective is to recover the initial velocity and pressure from lateral boundary observations, without assuming that the time-independent body force is known. To eliminate this unknown force, we differentiate the momentum equation with respect to time and then apply a Legendre polynomial-exponential time-dimensional reduction. This procedure yields a coupled system of elliptic equations for the expansion coefficients. We then construct a contractive map for this reduced system on a suitable admissible set equipped with a Carleman-weighted norm. Its fixed point yields an approximate solution of the time-dimensional reduction model, and the contraction property gives rise to a globally convergent Picard iteration. Finally, we present a numerical algorithm based on this framework and numerical experiments showing accurate reconstructions of the initial velocity and pressure from synthetic boundary data.

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.

Desk Editor's Note private letter to a colleague

The paper sets up a Carleman contraction iteration for recovering initial velocity and pressure in Navier-Stokes from boundary data with unknown force, but the global convergence hinges on an unverified contraction constant in the weighted norm.

read the letter

The paper gives a constructive way to recover initial velocity and pressure for the incompressible Navier-Stokes system from lateral boundary observations when the time-independent body force is unknown. They differentiate the momentum equation in time to cancel the force, apply a Legendre polynomial-exponential reduction to turn the problem into a coupled elliptic system for the coefficients, and then build a map on an admissible set that they equip with a Carleman-weighted norm. The fixed point of that map solves the reduced system, and the contraction gives a Picard iteration that converges globally on paper. They close with a numerical algorithm and tests on synthetic data that recover the initial fields accurately in the examples shown. This specific sequence of steps looks new for the inverse initial-data problem in fluids. The reduction is clean and the numerical part adds practical value that many inverse-problem papers skip. The soft spot is the contractivity claim. Carleman estimates produce constants that grow with the weight parameters, the domain, and the a-priori bounds, so it is not immediate that the map stays contractive with constant strictly less than one on the full admissible set without some implicit smallness restriction on the data. The abstract asserts global convergence, but the details of how they control the constant are not visible here and need direct checking in the proofs. If that step holds without extra assumptions, the result stands; if it requires small data, the claim needs adjustment. This is for researchers who work on Carleman methods or numerical inverse problems for fluid equations. A reader who already knows the standard Carleman toolkit for parabolic systems will see the value in the time-differentiation trick and the Legendre reduction. I would send it to peer review. The framework is concrete, the numerics are present, and the central estimates are checkable, so referees can settle whether the contraction works as stated.

Referee Report

1 major / 2 minor

Summary. The paper addresses an inverse initial-data recovery problem for the incompressible Navier-Stokes equations with unknown time-independent body force. Time differentiation of the momentum equation followed by a Legendre-polynomial exponential time-dimensional reduction produces a coupled elliptic system for the expansion coefficients. A contractive map is constructed for this system on an admissible set equipped with a Carleman-weighted norm; the fixed point furnishes an approximate solution of the reduced model and the contraction yields a globally convergent Picard iteration. A numerical algorithm based on the framework is presented together with experiments on synthetic boundary data.

Significance. If the contraction constant is shown to be strictly less than one on a non-empty admissible set without implicit smallness restrictions, the method would supply a globally convergent, theoretically justified algorithm for a practically relevant inverse problem in fluid dynamics. The numerical experiments demonstrating accurate reconstructions constitute a concrete strength that supports applicability.

major comments (1)

[Construction of the contractive map and associated Carleman estimates] The central claim that the Picard iteration converges globally rests on the map being contractive with constant strictly less than one in the Carleman-weighted norm on the admissible set for the time-differentiated, Legendre-reduced elliptic system. Carleman estimates typically produce constants that grow with the weight parameters, domain geometry, and a-priori bounds on coefficients; explicit verification that the constant remains <1 independently of the unknown force (and that the admissible set is non-empty) is required, as this is load-bearing for the global-convergence statement.

minor comments (2)

Clarify the precise regularity assumptions on the initial data, force, and boundary observations that guarantee the reduced elliptic system is well-defined.
In the numerical section, report the specific values of the Carleman weight parameters and the observed contraction ratios to allow readers to assess the practical margin below 1.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. The major comment focuses on the need for explicit verification that the contraction constant is strictly less than one independently of the unknown force and that the admissible set is non-empty. We address this point directly below and have made revisions to strengthen the presentation.

read point-by-point responses

Referee: [Construction of the contractive map and associated Carleman estimates] The central claim that the Picard iteration converges globally rests on the map being contractive with constant strictly less than one in the Carleman-weighted norm on the admissible set for the time-differentiated, Legendre-reduced elliptic system. Carleman estimates typically produce constants that grow with the weight parameters, domain geometry, and a-priori bounds on coefficients; explicit verification that the constant remains <1 independently of the unknown force (and that the admissible set is non-empty) is required, as this is load-bearing for the global-convergence statement.

Authors: We agree that explicit control over the contraction constant is essential to substantiate the global convergence of the Picard iteration. In the manuscript, time differentiation of the momentum equation eliminates the unknown time-independent body force from the reduced elliptic system, so the Carleman estimates and the resulting contraction constant depend only on the domain geometry, the viscosity, the Legendre truncation order, and the a priori bounds that define the admissible set. These quantities are independent of the specific force. The weight parameters in the Carleman estimates are chosen sufficiently large (depending only on the above quantities) to ensure the contraction constant is strictly less than one; the admissible set is non-empty by construction because it consists of all functions satisfying the same a priori bounds that the true solution is assumed to obey. To make this verification more explicit, we have added a new remark following the statement of the contraction mapping theorem that outlines the parameter selection procedure and confirms independence from the body force. We believe this addresses the concern while preserving the original argument. revision: partial

Circularity Check

0 steps flagged

Standard Carleman estimates and contraction mapping applied to derived elliptic system

full rationale

The paper differentiates the Navier-Stokes momentum equation in time to remove the unknown body force, applies a Legendre polynomial-exponential reduction to obtain a coupled elliptic system for the coefficients, and then constructs a map on an admissible set equipped with a Carleman-weighted norm whose contraction (with constant <1) is established via standard Carleman estimates. The fixed point of this map is shown to solve the reduced system approximately, and the contraction property directly implies global convergence of the Picard iteration by the Banach fixed-point theorem. No step equates a derived quantity to its own inputs by definition, renames a fitted parameter as a prediction, or relies on a load-bearing self-citation whose content is unverified; the admissible set is defined from a priori bounds independent of the final result, and the estimates are standard tools applied to the explicitly derived system.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Only the abstract is available, so the ledger is necessarily incomplete; the method assumes standard properties of the incompressible Navier-Stokes system and the existence of suitable Carleman estimates for the reduced elliptic operators.

axioms (2)

domain assumption The incompressible Navier-Stokes equations hold with a time-independent body force.
Invoked at the start of the inverse problem formulation.
standard math Carleman estimates exist and can be used to equip the admissible set with a norm that makes the solution map contractive.
Central to constructing the fixed-point argument.

pith-pipeline@v0.9.0 · 5439 in / 1141 out tokens · 33240 ms · 2026-05-10T16:39:32.175891+00:00 · methodology

discussion (0)

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Inverse initial data for nonlinear Schr\"odinger equation via Carleman estimates and the contraction principle
math.NA 2026-05 unverdicted novelty 5.0

A method combining Legendre basis reduction in time with a Carleman-estimate-based contraction mapping reconstructs initial data for nonlinear Schrödinger equations from boundary observations, with proven stability un...

Reference graph

Works this paper leans on

50 extracted references · 1 canonical work pages · cited by 1 Pith paper

[1]

Abney, T

R. Abney, T. T. Le, L. H. Nguyen, and C. Peters. A Carleman-Picard approach for reconstructing zero-order coefficients in parabolic equations with limited data.Applied Mathematics and Computation, 494:129286, 2025

2025
[2]

S. N. Antontsev and Kh. Khompysh. Inverse problems for a Boussinesq system for incompressible viscoelastic fluids.Mathematical Methods in the Applied Sciences, 46(9):11130–11156, jun 2023. 26 (a) relU (k) and relP (k) (b) ResU (k) and ResP (k) Figure 6: Test 3. Convergence history of the Picard iteration. The left image shows the relative changes relU(k) ...

2023
[3]

Beilina and M

L. Beilina and M. V. Klibanov.Approximate Global Convergence and Adaptivity for Coefficient Inverse Problems. Springer, New York, 2012

2012
[4]

A. L. Bukhgeim and M. V. Klibanov. Uniqueness in the large of a class of multidimensional inverse problems.Soviet Math. Doklady, 17:244–247, 1981

1981
[5]

Carleman

T. Carleman. Sur les syst` emes lin´ eaires aux deriv´ ees partielles du premier ordre a deux variables.C. R. Acad. Sci. Paris, 197:471–474, 1933

1933
[6]

L. Fang, R. Nei, and Z. Guo. Global well-posedness of the nonhomogeneous incompressible Navier– Stokes–Cahn–Hilliard system with Landau potential.Journal of Differential Equations, 445:113585, 2025

2025
[7]

Feireisl.Dynamics of Viscous Compressible Fluids

E. Feireisl.Dynamics of Viscous Compressible Fluids. Oxford University Press, Oxford, UK, 2004

2004
[8]

Fourestey and M

G. Fourestey and M. Moubachir. Solving inverse problems involving the navier–stokes equations dis- cretized by a lagrange–galerkin method.Computer Methods in Applied Mechanics and Engineering, 194(6):877–906, feb 2005

2005
[9]

Frehse, S

J. Frehse, S. Goj, and M. Steinhauer.L p-estimates for the Navier–Stokes Equations for Steady Com- pressible Flow.Manuscripta Mathematica, 116:265–275, 2005

2005
[10]

M. Fujii. Stationary Navier–Stokes equations on the half spaces in the scaling critical framework.Journal of Differential Equations, 435:113298, 2025

2025
[11]

He and H

J. He and H. Wu. Global well-posedness of a Navier–Stokes–Cahn–Hilliard system with chemotaxis and singular potential in 2D.Journal of Differential Equations, 297:47–80, oct 2021

2021
[12]

O. Yu. Imanuvilov and M. Yamamoto. Global uniqueness in inverse boundary value problems for the navier–stokes equations and lam´ e system in two dimensions.Inverse Problems, 31(3):035004, mar 2015

2015
[13]

V. A. Khoa, G. W. Bidney, M. V. Klibanov, L. H. Nguyen, L. Nguyen, A. Sullivan, and V. N. Astratov. Convexification and experimental data for a 3D inverse scattering problem with the moving point source. Inverse Problems, 36:085007, 2020

2020
[14]

V. A. Khoa, M. V. Klibanov, and L. H. Nguyen. Convexification for a 3D inverse scattering problem with the moving point source.SIAM J. Imaging Sci., 13(2):871–904, 2020

2020
[15]

M. V. Klibanov. Convexification of restricted Dirichlet to Neumann map.J. Inverse and Ill-Posed Problems, 25(5):669–685, 2017. 27

2017
[16]

M. V. Klibanov and J. Li.Inverse Problems and Carleman Estimates: Global Uniqueness, Global Convergence and Experimental Data. De Gruyter, 2021

2021
[17]

M. V. Klibanov and L. H. Nguyen. PDE-based numerical method for a limited angle X-ray tomography. Inverse Problems, 35:045009, 2019

2019
[18]

M. V. Klibanov, L. H. Nguyen, and H. V. Tran. Numerical viscosity solutions to Hamilton-Jacobi equations via a Carleman estimate and the convexification method.Journal of Computational Physics, 451:110828, 2022

2022
[19]

O. A. Ladyzhenskaya.The Mathematical Theory of Viscous Incompressible Flow. Gordon and Breach, New York, 1969

1969
[20]

R.-Y. Lai, G. Uhlmann, and J.-N. Wang. Inverse boundary value problem for the stokes and the navier–stokes equations in the plane.Archive for Rational Mechanics and Analysis, 215(3):811–829, mar 2015

2015
[21]

M. M. Lavrent’ev, V. G. Romanov, and S. P. Shishat·ski ˘i.Ill-Posed Problems of Mathematical Physics and Analysis. Translations of Mathematical Monographs. AMS, Providence: RI, 1986

1986
[22]

T. T. Le. Global reconstruction of initial conditions of nonlinear parabolic equations via the Carleman- contraction method. In D-L. Nguyen, L. H. Nguyen, and T-P. Nguyen, editors,Advances in Inverse problems for Partial Differential Equations, volume 784 ofContemporary Mathematics, pages 23–42. American Mathematical Society, 2023

2023
[23]

T. T. Le, M. V. Klibanov, L. H. Nguyen, A. Sullivan, and L. Nguyen. Carleman contraction mapping for a 1D inverse scattering problem with experimental time-dependent data.Inverse Problems, 38:045002, 2022

2022
[24]

T. T. Le and L. H. Nguyen. A convergent numerical method to recover the initial condition of nonlinear parabolic equations from lateral Cauchy data.Journal of Inverse and Ill-posed Problems, 30(2):265–286, 2022

2022
[25]

T. T. Le, L. H. Nguyen, T-P. Nguyen, and W. Powell. The quasi-reversibility method to numerically solve an inverse source problem for hyperbolic equations.Journal of Scientific Computing, 87:90, 2021

2021
[26]

T. T. Le, L. H. Nguyen, and H. V. Tran. A Carleman-based numerical method for quasilinear elliptic equations with over-determined boundary data and applications.Computers and Mathematics with Applications, 125:13–24, 2022

2022
[27]

T. T. Le, L. V. Nguyen, L. H. Nguyen, and H. Park. The time dimensional reduction method to deter- mine the initial conditions without the knowledge of damping coefficients.Computers and Mathematics with Applications, 166:77–90, 2024

2024
[28]

T. T. Le, C. B. Van, D. T. Dang, and L. H. Nguyen. Inverse initial data reconstruction for Maxwell’s equations via time-dimensional reduction method.Journal of Computational Physics, 559:114896, 2026

2026
[29]

P. L. Lions.Mathematical Topics in Fluid Dynamics, Volume 2: Compressible Models. Clarendon Press, Oxford, UK, 1998

1998
[30]

G. Liu. Determining the viscosity function from the boundary measurements for the stokes and the navier–stokes equations.Inverse Problems, 40(12):125011, dec 2024

2024
[31]

H. M. Nguyen and L. H. Nguyen. Cloaking using complementary media for the Helmholtz equation and a three spheres inequality for second order elliptic equations.Transaction of the American Mathematical Society, 2:93–112, 2015. 28

2015
[32]

L. H. Nguyen. An inverse space-dependent source problem for hyperbolic equations and the Lipschitz- like convergence of the quasi-reversibility method.Inverse Problems, 35:035007, 2019

2019
[33]

L. H. Nguyen. The Carleman contraction mapping method for quasilinear elliptic equations with over- determined boundary data.Acta Mathematica Vietnamica, 48:401–422, 2023

2023
[34]

L. H. Nguyen and M. V. Klibanov. Carleman estimates and the contraction principle for an inverse source problem for nonlinear hyperbolic equations.Inverse Problems, 38:035009, 2022

2022
[35]

L. H. Nguyen, Q. Li, and M. V. Klibanov. A convergent numerical method for a multi-frequency inverse source problem in inhomogenous media.Inverse Problems and Imaging, 13:1067–1094, 2019

2019
[36]

P. M. Nguyen, L. H. Nguyen, and H. T. Vu. Solving the inverse scattering problem via Carleman-based contraction mapping.Computers and Mathematics with Applications, 209:129–143, 2026

2026
[37]

Novotn´ y and I

A. Novotn´ y and I. Straˇ skraba.Introduction to the Mathematical Theory of Compressible Flow, volume 27 ofOxford Lecture Series in Mathematics and its Applications. Oxford University Press, Oxford, UK, 2004

2004
[38]

O’Connor, D

L. O’Connor, D. Lecoanet, E. H. Anders, K. C. Augustson, K. J. Burns, G. M. Vasil, J. S. Oishi, and B. P. Brown. Iterative methods for navier–stokes inverse problems.Physical Review E, 109(4):045108, apr 2024

2024
[39]

P. I. Plotnikov, E. V. Ruban, and J. Sokolowski. Inhomogeneous boundary value problems for compress- ible Navier–Stokes equations: Well-posedness and sensitivity analysis.SIAM Journal on Mathematical Analysis, 40(3):1152–1200, jan 2008

2008
[40]

P. I. Plotnikov and J. Sokolowski. Stationary boundary value problems for Navier–Stokes equations with adiabatic indexγ < 3 2.Doklady Mathematics, 70:535–538, 2004. In English; translated from Doklady Akademii Nauk, 397 (2004), pp. 166–169 (in Russian)

2004
[41]

P. I. Plotnikov and J. Sokolowski. Concentrations of solutions to time-discretized compressible Navier– Stokes equations.Communications in Mathematical Physics, 258:567–608, 2005

2005
[42]

M. H. Protter. Unique continuation for elliptic equations.Trans. Amer. Math. Soc., 95(1):81–91, 1960

1960
[43]

G. G. Stokes. On the effect of the internal friction of fluids on the motion of pendulums.Transactions of the Cambridge Philosophical Society, 9:8–106, 1851. Published by the Cambridge Philosophical Society, Cambridge, UK
[44]

Temam.Navier-Stokes Equations

R. Temam.Navier-Stokes Equations. North-Holland, Amsterdam, 1979

1979
[45]

D. D. Trong, C. V. Le, K. D. Luu, and L. H. Nguyen. Recovery of initial displacement and velocity in anisotropic elastic systems by the time dimensional reduction method.Journal of Computational Physics, 542:114371, 2025

2025
[46]

D. D. Trong, L. H. Nguyen, and H. T. Vu. Determining initial conditions for nonlinear hyperbolic equations with time dimensional reduction and the Carleman contraction principle.Inverse Problems, 40:125021, 2024

2024
[47]

Tsai.Lectures on Navier-Stokes Equations, volume 192 ofGraduate Studies in Mathematics

T.-P. Tsai.Lectures on Navier-Stokes Equations, volume 192 ofGraduate Studies in Mathematics. American Mathematical Society, Providence, RI, USA, 2018. MSC: Primary 35; 76

2018
[48]

C. B. Van, T. T. Le, and L. H. Nguyen. The inverse initial data problem for anisotropic Navier–Stokes equations via Legendre time reduction method.prepint arXiv:2507.16810, 2025

work page arXiv 2025
[49]

F. M. White.Viscous Fluid Flow. McGraw-Hill Mechanical Engineering. McGraw-Hill, 3 edition, 1991. 29

1991
[50]

W. Yan, Y. He, and Y. Ma. Shape reconstruction of an inverse boundary value problem of two- dimensional Navier–Stokes equations.International Journal for Numerical Methods in Fluids, 62(6):632– 646, feb 2010. 30

2010