On well-posedness theory of very weak solutions to Navier-Stokes equations on irregular domains with nonhomogeneous Dirichlet boundary data

Siran Li; Xiangxiang Su; Xiaojin Bai

REVIEW 3 major objections 5 minor 24 references

On well-posedness theory of very weak solutions to Navier-Stokes equations on irregular domains with nonhomogeneous Dirichlet boundary data

T0 review · 3 major / 5 minor · reviewed 2026-08-01 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read On bounded Lipschitz domains whose boundary graphs have small Sobolev-multiplier norm, the three-dimensional Navier-Stokes equations with nonhomogeneous Dirichlet data admit a unique very weak solution in L^s_t L^q_x, locally in time, and g

desk verdict The paper's internal argument is clean, but Theorem 1.1 is a deduction from an unverified external preprint; it deserves review, not blind acceptance. read the letter →

arxiv 2607.23234 v1 pith:LEAU2T3U submitted 2026-07-25 math.AP

Xiaojin Bai , Siran Li , Xiangxiang Su This is my paper · ORCID

classification math.AP MSC 35Q3076D0576D03

keywords veryweaksolutionsNavier-StokesequationsSobolevmultipliersLipschitzdomainsnonhomogeneousDirichletboundarydataStokessemigroupwell-posednessmaximalregularity

topics Navier-Stokes Regularity

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper establishes a well-posedness theory for very weak solutions of the three-dimensional Navier-Stokes equations on a class of rough bounded domains: bounded Lipschitz domains whose boundary charts belong to a Sobolev-multiplier class with sufficiently small norm. The main theorem states that, under natural scaling assumptions on the forcing, boundary data, and initial data, a unique very weak solution exists in L^s(0,T*;L^q(Ω)); the solution exists on the whole half-line if the data are small. This matters because the standard theory, which relies on W^{2,q} elliptic regularity, is limited to C^{2,1} domains, and the authors' substitution — semigroup estimates and a duality-based weak formulation that avoids boundary integration by parts — extends existence and uniqueness to a strictly larger class of irregular domains.

What carries the argument

The load-bearing mechanism is the Sobolev-multiplier boundary class M_W^{1+α,ρ}(ε), a class of bounded Lipschitz domains strictly between Lipschitz and C^{1,α}; on this class the Stokes-Dirichlet operator A_p is assumed to have the full set of semigroup properties: bounded Helmholtz projection, square-root domain identification D(A_p^{1/2}) = W^{1,p}_{0,σ}, maximal L^p regularity, exponential decay, and a gradient estimate. These imported estimates allow the construction of stationary fields absorbing the force and boundary data, and define the very weak formulation via duality with D(A_{q'}) test functions. The key operator is the bilinear map Q(w,v) = K(w⊗v), built from the Stokes semigrou

What would settle it

One concrete test: take a bounded Lipschitz domain in the multiplier class M_W^{1+α,ρ}(ε) with ε arbitrarily small, and solve the Stokes problem with a forcing term in L^p; if the Stokes semigroup fails any of the estimates of Proposition 2.3 (for instance, if the square-root domain D(A_p^{1/2}) ≠ W^{1,p}_{0,σ} for some p in the exponent set, or if maximal L^p regularity fails), then the theorem is false. Conversely, verifying these estimates for a nontrivial family of such domains would confirm the theory.

Watch

Extended reading notes

Core claim

The paper's central claim is that the very weak formulation of the Navier-Stokes system with nonhomogeneous Dirichlet boundary data is well-posed on domains of Sobolev-multiplier class M_W^{1+α,ρ}(ε) for sufficiently small ε, with a unique solution in L^s_t L^q_x. The proof reduces the problem to an equivalent integral equation u = U - Q(u,u), built from the Stokes semigroup, and defines the inverse Stokes operator on the rough nonlinear term by duality; the argument closes by a contraction. Local-in-time existence follows from the decay of the semigroup term as t→0, and global-in-time existence from an explicit smallness condition on the data.

Load-bearing premise

The whole proof assumes that the Stokes-Dirichlet operator on Sobolev-multiplier domains—with the specific smallness of the multiplier norm—satisfies the full set of semigroup estimates (bounded Helmholtz projection, square-root domain identification, maximal regularity, exponential decay); if those imported estimates fail for any exponent in the set used, the construction of the solution collapses.

Editorial extensions

If this is right

The class of domains on which very weak solutions to the Navier-Stokes equations are known to be well-posed is enlarged from C^{2,1} domains to a class of bounded Lipschitz domains with small Sobolev-multiplier boundary norm.
For initial data, forcing, and boundary data satisfying the natural scaling regularity, there is a unique very weak solution on a time interval (0,T*) whose length depends on the data; the estimate (1.4) holds with a constant independent of T.
If, in addition, the data are small in the sense of (1.5), the unique very weak solution exists for all positive times.
The solution satisfies the divergence-free condition and the boundary condition in a weak sense, and a pressure distribution is recovered from the equation.
The same proof, with the semigroup estimates replaced by Besov-space analogues, yields well-posedness for data in Besov spaces B^s_{p,q}(Ω) in the indicated range.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the same Stokes estimates hold on wider classes of rough domains than the Sobolev-multiplier class, the argument would transfer directly, since the construction only uses semigroup properties and duality.
The partition argument used for uniqueness gives a quantitative relation between the L^s_t L^q_x norm of the solution and the length of the existence interval; turning it around could yield a blow-up criterion for very weak solutions on irregular domains.
The duality-based treatment of the nonlinear term is not tied to the specific form of the quadratic term; analogues for other incompressible models (e.g., MHD or Boussinesq) would follow if their semigroups satisfy the same set of estimates.
A concrete check of the theorem's hypothesis would be to numerically or analytically verify the Stokes semigroup estimates (particularly the square-root domain property) on an explicit family of multiplier-class domains; non-uniformity across the exponent set would narrow the admissible q-range.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague,

The headline: the paper gives a coherent extension of very weak Navier–Stokes well-posedness to bounded Lipschitz domains whose boundary charts have sufficiently small Sobolev multiplier norm, and the internal argument from the auxiliary fields through the fixed point is clean. But Theorem 1.1 is not self-contained: every load-bearing Stokes estimate is imported from the unpublished Breit–Gaudin preprint [5]. If that toolkit holds, the theorem plausibly holds; if it does not, the construction collapses.

What is genuinely new is the duality-based formulation (Definition 3.3) that bypasses the W^{2,q} trace theory used on C^{2,1} domains, and the use of the R_q operator to interpret A_q^{-1} P_q div(u⊗u). The proof of the equivalence between this very weak notion and the integral equation (Prop 5.3) is careful, and the index relations in Lemma 5.1 check out (γ = 1/2 + 3/(2q) < 1 since q > 3). The contraction argument is standard but executed correctly, including the partition-of-unity uniqueness.

The soft spot is exactly where the stress-test points. Proposition 2.3 supplies the Helmholtz decomposition, square-root identification, maximal regularity, and the gradient estimate (2.10). All of it is cited from [5]. Appendix A shows how [5]'s results would be combined, but it does not prove them or even state the precise smallness assumptions in the multiplier norm under which they hold uniformly in p ∈ {q, q', r, r'} and in α, ρ. That uniformity matters, since Theorem 1.1 advertises a universal ε0. Without it, the statement is ambiguous. The paper is honest about the dependence — the authors flag it in the introduction and appendix — but honesty does not close the gap.

One more small thing: the abstract promises that the multiplier class contains bounded Lipschitz domains with small Lipschitz constants. I don't see that implication proved or even stated as a proposition in the text. It is plausible from Sobolev multiplier theory, but it should be shown.

Who should read this: people working on very weak solutions or Stokes theory on rough domains will want to know this framework, assuming [5] appears. It deserves a serious referee: the internal logic is sound, and the conditional result is a step beyond Farwig–Galdi–Sohr. My recommendation is to send it out, but instruct the referee to treat Proposition 2.3 as the crux and require the authors to either prove the needed estimates or explicitly rephrase Theorem 1.1 as a conditional statement with precise hypotheses on ε and the external a priori bounds.

Referee Report

3 major / 5 minor

Summary. The paper develops a well-posedness theory for very weak solutions of the Navier–Stokes equations on bounded Lipschitz domains whose boundary charts belong to a Sobolev multiplier class M_W^{1+α,ρ}(ε) with sufficiently small multiplier norm. The authors define very weak solutions via duality with the Stokes–Dirichlet operator, construct stationary auxiliary fields for the force and boundary data, establish the linear theory, and then reformulate the nonlinear problem as a fixed-point equation in L^s(0,T;L^q(Ω)). Local well-posedness for arbitrary data and global well-posedness under a smallness condition are claimed. The proof builds on a package of Stokes estimates — Helmholtz decomposition, square-root domain identification, maximal regularity, semigroup decay, and gradient estimates — collected in Proposition 2.3, which is only sketched in Appendix A and is said to follow from the unpublished preprint [5] by Breit and Gaudin.

Significance. If the underlying Stokes estimates hold, the paper makes a substantial advance by extending the very weak solution theory from C^{2,1} domains to a class of irregular Lipschitz domains that is strictly between Lipschitz and C^{1,α}. The internal argument is largely coherent: the duality-based very weak formulation, the construction of the stationary auxiliary fields, and the reduction to the integral equation are carefully presented. The paper does not fit parameters or rely on circular assumptions. However, the entire edifice rests on Proposition 2.3, whose proof is not included and whose uniformity in the relevant exponents is not demonstrated. The contribution is therefore conditional on the external toolkit in [5].

major comments (3)

[Section 2, Proposition 2.3 and Appendix A] Proposition 2.3 is load-bearing: every subsequent result uses it, including the Helmholtz decomposition in (3.8), the auxiliary Stokes solves in Propositions 3.1–3.2, the linear theory in Proposition 4.2, the definition of R_q in (5.1)–(5.2), and the K-estimate in Lemma 5.1. Appendix A does not prove these estimates; it only explains how they follow from [5, Theorems 4.38, 6.5; Propositions 4.40, 6.19; Meta-Theorem 6.20], and [5] is an unpublished preprint. As submitted, Theorem 1.1 is contingent on the correctness and completeness of [5]. The authors should either supply a full proof of Proposition 2.3 or state precisely the external theorem and verify all of its hypotheses, including the exact smallness condition on ε.
[Theorem 1.1 and Proposition 2.3 (uniformity of ε0)] Theorem 1.1 asserts a universal ε0 > 0 that works for all α∈(0,1), ρ∈[1,∞], and for the whole set of exponents p∈{q,q′,r,r′}. Proposition 2.3, however, only says that for each p there exists ε>0 sufficiently small. The proof requires one smallness threshold that simultaneously gives (2.8), (2.10), and the square-root identification for all four exponents. The uniformity of the smallness condition in p, α, and ρ is not shown in Appendix A. If the threshold in [5] is not uniform, then ε0 and the contraction threshold μ0 may depend on p in a way that invalidates the fixed-point argument in Section 6.
[Section 5, Eqs. (5.1)–(5.5) and Lemma 5.1] The boundedness of R_q and K is the only mechanism that makes the nonlinear term meaningful. It depends on the gradient estimate (2.10), i.e., on the bound ||∇e^{-tA_{q'}}ψ||_{L^{(q/2)'}} ≤ C t^{-1/2-3/(2q)} ||ψ||_{L^{q'}}. The exponent computation is correct for q>3, but (2.10) is cited from [5] and is not proved. If the constant in this estimate depends on the target exponent (q/2)' = q/(q-2) in an uncontrolled way, then (5.3) and (5.6) would fail and the bilinear estimate (5.7) would break. The authors should verify that (2.10) holds with a constant depending only on q and the domain class, uniformly in α and ρ, for the stated range.

minor comments (5)

[Theorem 1.1, display (1.4)] The norm in the estimate is written as ∥u∥_{L^s(0,T;L^q(Ω))}, but the theorem has already introduced T*∈(0,T]; it should be L^s(0,T*;L^q(Ω)).
[Definition 2.2] The phrase 'ε>0 sufficiently small' appears before the proposition that quantifies ε. This is ambiguous; state explicitly that the smallness condition is the one appearing in Proposition 2.3.
[Remark 1.2] The claimed adaptation to Besov spaces is stated without proof. Since it is not used in the main theorem, it should be either removed or clearly labeled as a conjecture, or the necessary Besov analogues should be stated.
[Throughout] There are numerous formatting issues: missing spaces in 'C 1 c', nonstandard 'L s t Lq x', and inconsistent use of L^s(0,T;L^q) versus L^s(0,T*;L^q). These should be corrected.
[References] Reference [5] is an arXiv preprint with DOI. If accepted for publication, the authors should update the reference to the published version and, more importantly, include the precise theorem statements they rely on.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the result is a fixed-point consequence of externally imported Stokes estimates, not an input renamed as a conclusion.

full rationale

The derivation chain in Sections 3–6 is: construct auxiliary fields (Propositions 3.1–3.2) using the Helmholtz decomposition and square-root identification imported from Proposition 2.3; solve the linear problem (Proposition 4.2) using maximal regularity of A_q; define the nonlinear operator K via the semigroup gradient estimate (2.10); prove equivalence to the integral equation (Proposition 5.3); and close by contraction (Proposition 6.1). None of these steps postulates the existence of a very weak NSE solution or the estimate (1.4) as an input. The only heavy analytic input is Proposition 2.3, and it is imported from the external preprint [5] (Breit–Gaudin) and Maz'ya–Shaposhnikova [22]; the present authors are not the authors of [5], so no self-citation chain is load-bearing. The external Stokes estimates are an independent benchmark: if they are valid, the proof yields a new Navier–Stokes result; if they fail, the theorem collapses, but that is a correctness or verification risk, not circularity. The smallness constants ε_0 and μ are chosen by the proof and are not part of the data. No fitted parameter is renamed as a prediction, and no uniqueness theorem from the same authors is used to forbid alternatives. I therefore find no step in which a conclusion is equivalent by construction to its inputs.

Assumptions & free parameters 1 free parameters · 4 assumptions · 0 invented entities

The paper introduces no new particles or physical entities. Its free parameter and domain axioms are analytical: the smallness threshold ε0 and the imported Stokes-operator toolkit from Breit–Gaudin, plus standard background theorems. No fitted numerical data appear.

free parameters (1)

ε0 (smallness threshold for boundary multiplier norm) = not explicit; 'sufficiently small'
Theorem 1.1 requires ε < ε0, and Proposition 2.3 only states the existence of some small ε for each p ∈ E_q. The uniformity of ε0 across q, q′, r, r′ is asserted without an explicit value or derivation; this quantitative premise makes the domain class work.

assumptions (4)

domain assumption On M_W^{1+α,ρ}(ε) domains, the Stokes–Dirichlet operator A_p satisfies the full set of estimates in Proposition 2.3: bounded Helmholtz projection, duality, invertibility, square-root domain identification, maximal regularity, exponential decay, and gradient bound (2.10).
Invoked throughout Sections 2–6; Appendix A only sketches the proof and delegates to the unpublished Breit–Gaudin preprint [5]. If these estimates fail, the auxiliary fields, duality formulation, and fixed-point argument collapse.
domain assumption Bounded Lipschitz domains with sufficiently small Lipschitz constants are contained in the Sobolev-multiplier class M_W^{1+α,ρ}(ε) for ε < ε0.
Asserted in the abstract but not proved in the body; needed to substantiate the advertised special case.
standard math Bogovskii operator gives a bounded right-inverse of the divergence on bounded Lipschitz domains; standard trace and Helmholtz decomposition facts.
Used in Proposition 3.2 to construct the solenoidal boundary extension U_g.
standard math Weis' maximal-regularity theorem, H∞-calculus for sectorial operators, Hardy–Littlewood–Sobolev inequality, and de Rham's theorem.
Used in Appendix A, Lemma 5.1, and the pressure reconstruction at the end of the proof of Theorem 1.1.

how reviews work

0 comments

Cite this review

Pith. "Pith review of On well-posedness theory of very weak solutions to Navier-Stokes equations on irregular domains with nonhomogeneous Dirichlet boundary data." pith.science (2026). https://pith.science/paper/LEAU2T3U

@misc{pith2026260723234,
  author       = {Pith},
  title        = {Pith review of: On well-posedness theory of very weak solutions to Navier-Stokes equations on irregular domains with nonhomogeneous Dirichlet boundary data},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LEAU2T3U}},
  note         = {Machine review of arXiv:2607.23234}
}

read the original abstract

The well-posedness theory of very weak solutions is a central topic in mathematical hydrodynamics, especially in the regularity theory for Navier-Stokes equations. It has been fully developed for incompressible fluid flows on bounded domains in R^3 of C^{2,1}-regularity. In this paper, based on the analytic theories in [D. Breit and A. Gaudin, ArXiv Preprint: 2511.19091 (2025)] and [V.G. Maz'ya and T.O. Shaposhnikova, Vol.337, Grundlehren der mathematischen Wissenschaften (2009)], we establish the well-posedness theory of very weak solutions to the Navier-Stokes equations on bounded Lipschitz domains whose boundary has local graphing functions with sufficiently small Sobolev multiplier norm, which contain the bounded Lipschitz domains with sufficiently small Lipschitz constants as a special case.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references

[5]

Breit and A

D. Breit and A. Gaudin. Optimal regularity results for the Stokes–Dirichlet problem. arXiv preprint, 2025, 10.48550/arXiv.2511.19091

doi:10.48550/arxiv.2511.19091 2025
[1]

H. Amann. Nonhomogeneous Navier-Stokes equations with integrable low-regularity data. InNonlinear problems in mathematical physics and related topics, Int. Math. Ser., pages 1–28. Kluwer/Plenum, New York, 2002

2002
[2]

H. Amann. Navier-Stokes equations with nonhomogeneous Dirichlet data.J. Non- linear Math. Phys., 10:1–11, 2003

2003
[3]

D. Breit. Regularity results in 2D fluid-structure interaction.Math. Ann., 388:1495– 1538, 2024

2024
[4]

D. Breit. Partial boundary regularity for the Navier-Stokes equations in irregular domains.J. Funct. Anal., 289:111188, 2025

2025
[6]

V. Coscia. Existence and uniqueness of very weak solutions to the steady-state Navier-Stokes problem in Lipschitz domains.J. Math. Fluid Mech., 19:819–829, 2017

2017
[7]

Farwig, G

R. Farwig, G. P. Galdi, and H. Sohr. Very weak solutions of stationary and insta- tionary Navier-Stokes equations with nonhomogeneous data. InNonlinear elliptic and parabolic problems, volume 64 ofProgr. Nonlinear Differential Equations Appl., pages 113–136. Birkh¨ auser, Basel, 2005

2005
[8]

Farwig, G

R. Farwig, G. P. Galdi, and H. Sohr. A new class of weak solutions of the Navier– Stokes equations with nonhomogeneous data.J. Math. Fluid Mech., 8:423–444, 2006

2006

Show all 24 references

[9]

Farwig, G

R. Farwig, G. P. Galdi, and H. Sohr. Very weak solutions and large uniqueness classes of stationary Navier-Stokes equations in bounded domains ofR 2.J. Differential Equations, 227:564–580, 2006

2006
[10]

Farwig, H

R. Farwig, H. Kozono, and H. Sohr. Very weak solutions of the Navier-Stokes equa- tions in exterior domains with nonhomogeneous data.J. Math. Soc. Japan, 59:127– 150, 2007

2007
[11]

Gabel and P

F. Gabel and P. Tolksdorf. The stokes operator in two-dimensional bounded Lipschitz domains.J. Differential Equations, 340:227–272, 2022

2022
[12]

G. P. Galdi.An Introduction to the Mathematical Theory of the Navier–Stokes Equa- tions: Steady-State Problems. Springer Monographs in Mathematics. Springer, New York, second edition, 2011

2011
[13]

Geng and Z

J. Geng and Z. Shen. Resolvent estimates for the Stokes operator in bounded and exteriorC 1 domains.Math. Ann., 391:1467–1503, 2025

2025
[14]

Geng and Z

J. Geng and Z. Shen. Resolvent estimates inL ∞ for the Stokes operator in nonsmooth domains.Invent. Math., 243:657–701, 2026. Very weak solutions to the Navier-Stokes equations22

2026
[15]

Y. Giga. Domains of fractional powers of the Stokes operator inL r spaces.Arch. Ration. Mech. Anal., 89:251–265, 1985

1985
[16]

Y. Giga. Solutions for semilinear parabolic equations inL p and regularity of weak solutions of the Navier–Stokes system.J. Differential Equations, 62:186–212, 1986

1986
[17]

Giga and T

Y. Giga and T. Miyakawa. Solutions inL r of the Navier–Stokes initial value problem. Arch. Ration. Mech. Anal., 89:267–281, 1985

1985
[18]

Giga and H

Y. Giga and H. Sohr. AbstractL p estimates for the Cauchy problem with applications to the Navier–Stokes equations in exterior domains.J. Funct. Anal., 102:72–94, 1991

1991
[19]

Haase.The Functional Calculus for Sectorial Operators, volume 169 ofOperator Theory: Advances and Applications

M. Haase.The Functional Calculus for Sectorial Operators, volume 169 ofOperator Theory: Advances and Applications. Birkh¨ auser, Basel, 2006

2006
[20]

Kunstmann and L

P.C. Kunstmann and L. Weis. New criteria for theH ∞-calculus and the Stokes operator on bounded Lipschitz domains.J. Evol. Equ, 17:387–409, 2017

2017
[21]

E. H. Lieb and M. Loss.Analysis, volume 14 ofGraduate Studies in Mathematics. American Mathematical Society, Providence, RI, second edition, 2001

2001
[22]

V. G. Maz’ya and T. O. Shaposhnikova.Theory of Sobolev multipliers: With appli- cations to differential and integral operators, volume 337 ofGrundlehren der mathe- matischen Wissenschaften. Springer-Verlag, Berlin, 2009

2009
[23]

Z. Shen. Resolvent estimates inL p for the Stokes operator in Lipschitz domains. Arch. Ration. Mech. Anal., 205:395–424, 2012

2012
[24]

Tolksdorf

P. Tolksdorf. On theL p-theory of the Navier-Stokes equations on three-dimensional bounded Lipschitz domains.Math. Ann., 371:445–460, 2018

2018

Pith tools

Reviewed August 1, 2026 · model on record in the stance chip above.

[1] [5]

Breit and A

D. Breit and A. Gaudin. Optimal regularity results for the Stokes–Dirichlet problem. arXiv preprint, 2025, 10.48550/arXiv.2511.19091

doi:10.48550/arxiv.2511.19091 2025

[2] [1]

H. Amann. Nonhomogeneous Navier-Stokes equations with integrable low-regularity data. InNonlinear problems in mathematical physics and related topics, Int. Math. Ser., pages 1–28. Kluwer/Plenum, New York, 2002

2002

[3] [2]

H. Amann. Navier-Stokes equations with nonhomogeneous Dirichlet data.J. Non- linear Math. Phys., 10:1–11, 2003

2003

[4] [3]

D. Breit. Regularity results in 2D fluid-structure interaction.Math. Ann., 388:1495– 1538, 2024

2024

[5] [4]

D. Breit. Partial boundary regularity for the Navier-Stokes equations in irregular domains.J. Funct. Anal., 289:111188, 2025

2025

[6] [6]

V. Coscia. Existence and uniqueness of very weak solutions to the steady-state Navier-Stokes problem in Lipschitz domains.J. Math. Fluid Mech., 19:819–829, 2017

2017

[7] [7]

Farwig, G

R. Farwig, G. P. Galdi, and H. Sohr. Very weak solutions of stationary and insta- tionary Navier-Stokes equations with nonhomogeneous data. InNonlinear elliptic and parabolic problems, volume 64 ofProgr. Nonlinear Differential Equations Appl., pages 113–136. Birkh¨ auser, Basel, 2005

2005

[8] [8]

Farwig, G

R. Farwig, G. P. Galdi, and H. Sohr. A new class of weak solutions of the Navier– Stokes equations with nonhomogeneous data.J. Math. Fluid Mech., 8:423–444, 2006

2006