REVIEW 3 major objections 4 minor 44 references
Stress concentration between two adjacent rigid particles in Navier-Stokes flow
T0 review · 3 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read The paper establishes the exact blow-up rates for the velocity gradient and pressure in stationary Navier-Stokes flow between two nearly touching rigid particles, with matching upper and lower bounds in two and three dimensions.
desk verdict First optimal stress blow-up rates for stationary Navier-Stokes flow between nearly touching rigid particles, matching the Stokes rates, but the main upper-bound proof only works for a symmetric quadratic gap and a normalized coefficient; still worth refereeing. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The proof rests on a decomposition of the solution into components attached to the six (three in 2D) rigid displacements of each particle, plus a background Stokes flow. Only the component $(u_1^3,p_1^3)$ carries the Navier-Stokes nonlinearity; all other components solve Stokes equations and are controlled by the sharp estimates already available for the Stokes problem. The geometric engine is a Keller-type function $k(x)=x_3/\delta(x')$, where $\delta(x')=\varepsilon+h_1(x')+h_2(x')$ is the local gap between the two surfaces. Auxiliary fields $v_i^\alpha$ are built from $k$ so that $\mu\Delta v_i^\alpha-\nabla p_i^\alpha$ is small in the neck, and the remainder $w=u-v$ is controlled through a Caccioppoli-type inequality, a dyadic iteration, and rescaled $W^{2,\infty}$ estimates for the Stokes operator. The unknown coefficients $C_i^\alpha$ in the rigid-motion expansion are fixed by force and torque balance, which leads to a linear system whose diagonal dominance gives the key differences $|C_1^\alpha-C_2^\alpha|$. The lower bounds are obtained by comparing the Navier-Stokes solution with the Stokes solution at the narrowest point.
What would settle it
Numerically solve the stationary Navier-Stokes system for two equal unit spheres in $\mathbb{R}^3$ with $\phi=(0,0,x_3)$ and measure $|\nabla u|$ at the midpoint $(0,0)$ for $\varepsilon=10^{-2},10^{-3},10^{-4}$; if $\varepsilon|\log\varepsilon|\,|\nabla u|$ does not stay within two positive constants as $\varepsilon$ shrinks, the optimal rate in Theorem 1.4 is false.
Extended reading notes
Core claim
The central claim is that the optimal blow-up rates for stationary Navier-Stokes flow between two nearly touching convex rigid particles coincide with the Stokes rates. In three dimensions, Theorem 1.1 gives $|\nabla u(x)|\le C(1+|\log\varepsilon|\,|x'|)/(|\log\varepsilon|(\varepsilon+|x'|^2))\|\phi\|_{C^{2,\alpha}}$ in the neck region, together with matching bounds for $|\nabla^2 u|$ and $|\nabla p|$; Theorem 1.4 shows that for two unit balls with boundary data $\phi=(0',x_3)$, $|\nabla u(0',x_3)|\ge 1/(C\varepsilon|\log\varepsilon|)$ for $|x_3|\le\varepsilon$. In two dimensions, the upper bound is $C/\sqrt{\varepsilon}$ and the lower bound is $1/(C\sqrt{\varepsilon})$. The authors state that these are the first such estimates for the stationary Navier-Stokes system and that the Cauchy stress tensor inherits the same blow-up rates.
Load-bearing premise
The load-bearing premise is that the two facing surfaces have the same leading curvature at the closest point, reducing the local gap to exactly $\varepsilon+|x'|^2$; the paper does not supply a perturbation argument for two different convex bodies, even though the stated assumptions only bound their $C^3$ norms.
Editorial extensions
If this is right
- The Cauchy stress tensor $\sigma[u,p]=2\mu e(u)-pI$ blows up at the same rate as $|\nabla u|$ in the neck, giving a quantitative statement about force concentration between nearby particles.
- The matching lower bounds show the upper-bound rates are sharp for the canonical configuration of two equal balls (or disks) with axial boundary data.
- In dimensions $d\ge4$ the paper records the rate $C/\varepsilon^{3/2}$ with minor modifications, so the logarithmic correction in 3D is tied to the borderline dimension.
- The result extends the known optimal Stokes rates to the nonlinear stationary problem, implying that inertia enters only at lower order in the near-contact singularity.
Reading between the lines
- Editorial extension: the theorem's equal-curvature assumption (1.6) is the main restriction; for two convex bodies with different principal curvatures the same rate is plausible, but that case is not reduced here and would need a separate argument.
- Editorial extension: the lower bound is proved only for two equal balls and one special boundary field, so optimality for general shapes and general data remains open even if the upper bounds cover them.
- Editorial extension: a time-dependent version of this problem, where the gap closes under an applied force, would translate these stationary blow-up rates into a resistance law; the authors signal a forthcoming evolution result, but the stationary analysis is the first step.
- Editorial extension: one could test the rates numerically with a spectral or finite-element solver on two spheres; the predicted product $\varepsilon|\log\varepsilon|\,|\nabla u|$ should be order one and independent of $\varepsilon$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies stationary incompressible Navier-Stokes flow in the narrow region between two closely placed convex rigid particles. It claims pointwise upper bounds for |∇u|, |∇^2u| and |∇p| in the neck region in dimensions three and two, with blow-up rates (ε|log ε|)^{-1} and ε^{-1/2}, respectively, together with matching lower bounds for the gradient when the two particles are unit balls. The strategy decomposes the solution into linear Stokes components and one nonlinear Navier-Stokes component, constructs Keller-type auxiliary fields for the quadratic symmetric gap, proves energy and Caccioppoli estimates by iteration, solves a linear system for the rigid-body coefficients, and compares the Navier-Stokes solution with the corresponding Stokes solution in the lower-bound proof.
Significance. If the estimates were proved under the stated C^3 convex geometry, this would be a significant contribution: it would give the first pointwise stress-concentration estimates for stationary Navier-Stokes flow between adjacent rigid particles, would match the known optimal Stokes rates, and would improve the pressure estimate in [35]. The proof contains many careful energy and local regularity estimates, and the lower-bound comparison with the Stokes results [34,35] is a sensible strategy. However, as written, the main theorems are proved only for a special exact quadratic symmetric geometry and for a normalized rigid coefficient, and no reduction or perturbation argument is supplied for the announced C^3 convex domains. The current significance is therefore limited to the special geometry actually treated.
major comments (3)
- [Section 2.2, Eqs. (2.12)-(2.13), (4.4)-(4.5)] The proof of the 3D upper bound is carried out under the explicit reduction h1=h2=|x'|^2/2. This is not a harmless normalization: the Keller function k=x3/δ(x') takes the values ±1/2 on the two particle boundaries only when h1=h2, and the identities ∂_{x_j}δ=2x_j used throughout (4.1)-(4.5) require the exact quadratic form. For surfaces satisfying only assumption (1.6), with possibly different leading curvatures and O(|x'|^3) remainders, δ_{x_j} has an extra radius-dependent factor and k no longer has the required boundary values. Consequently the cancellations (4.4)-(4.5), the bound (4.10), Lemma 4.2, and Proposition 2.6 do not follow as written. Since these estimates feed directly into Proposition 2.7 and the final proof of Theorem 1.1, the main upper bound is not proved for the class of C^3 convex domains announced. This also affects Theorem 1.4, whose balls satisfy h=1-√(1-r^2) with nonzero higher-order terms, so the proof of the upper bound used in (5.7) is not available for that geometry without an additional perturbation argument.
- [Section 2.1, Eqs. (2.3)-(2.7)] The decomposition into Stokes components and one nonlinear component is introduced, but the text states: "Here we only consider the special case that C_1^3=1; for the other cases, a different decomposition may be required." Since C_1^3 is one of the unknown rigid-body coefficients determined by the solution and by the data, the condition C_1^3=1 cannot be imposed without loss of generality. In particular, if C_1^3=0, the nonlinear component (u_1^3,p_1^3) does not satisfy the stated equation (2.7). Thus Theorem 1.1 is not established for arbitrary boundary data φ. The same issue occurs in the two-dimensional proof in Section 6, where the normalization C_2^1=1 is used without justification.
- [Section 6, Propositions 6.4-6.11] The two-dimensional upper bound theorem is proved only after the same simplification h1(x1)=h2(x1)=x1^2/2, and with the same unjustified normalization C_2^1=1. The statement of Theorem 6.1, however, assumes general C^3 convex domains under (1.6). Since the two-dimensional auxiliary fields, the estimates for |C_1^α-C_2^α| in Proposition 6.11, and the final proof of Theorem 6.1 all rely on this exact quadratic geometry, Theorem 6.1 is not proved for the stated domain class. Theorem 6.13, the two-dimensional lower bound, is likewise affected because its proof relies on the unproved upper bound for the same ball geometry, whose surfaces are not exactly quadratic.
minor comments (4)
- [Section 4.4, proof of (2.10)] In the proof of Proposition 2.2 the reference "in view of (6.3)" appears to be a typographical error; the relevant equation is (2.4), since (6.3) is the two-dimensional Stokes decomposition introduced later.
- [Proposition 2.5] The function spaces in Proposition 2.5 are stated as C^{2,γ}(Ω;R^2) and C^{1,γ}(Ω), but the problem is three-dimensional; the target spaces should be R^3 and the ambient domain should be consistent with Section 2.
- [Theorems 1.4 and 6.13] The notation for the two unit balls, e.g. B_1(0,1+ε/2) and B_1(0',1+ε/2), is ambiguous: the center and radius should be written explicitly, such as B_1((0',1+ε/2)) and B_1((0',-1-ε/2)), to avoid confusion with balls centered at the origin.
- [Theorem 1.1 and Section 1.1] Theorem 1.1 assumes u∈W^{1,2}(D)∩C^2(Ωbar) and p∈L^2(D)∩C^1(Ωbar), but the only existence statement cited in Section 1.1 is Ladyzhenskaya's weak solution, which does not provide this regularity. The authors should clarify explicitly whether Theorem 1.1 is an a priori estimate for sufficiently regular solutions or a statement about any solution of (1.1)-(1.3), and, in the latter case, should comment on the regularity needed to justify the pointwise estimates.
Circularity Check
No significant circularity: the Navier-Stokes stress-concentration estimates are obtained by a new iteration argument over independent published Stokes estimates, not by assuming the target bounds.
full rationale
The claimed derivation chain is not circular. The proof decomposes the Navier-Stokes solution into linear Stokes components and one nonlinear component in (2.2)-(2.5). The estimates for the linear components are quoted from the first author's earlier Stokes papers [34,35] (Propositions 2.2-2.4 and 2.6 for i=2). These are published, parameter-free results about the linear Stokes system; their assumptions do not include the Navier-Stokes bound being proved, and their conclusions are the Stokes blow-up rates rather than the target theorem. Under the stated rules they therefore count as independent support and do not raise the circularity score. The genuinely new content is the treatment of the nonlinear component (u^3_1, p^3_1) in Sections 3-4: the energy boundedness lemma (Lemma 3.1), the Caccioppoli inequality (Lemma 3.2), the rescaling proposition (Proposition 3.5), and the local energy estimates (Lemmas 4.1-4.2). The nonlinear term u·∇u is treated as a forcing term and is bounded only after the linear estimates are in hand; no term is fitted to reproduce the predicted blow-up rates. The coefficient differences |C^α_1 - C^α_2| in Proposition 2.7 are solved from the force/torque algebraic system (4.22)-(4.23) via Cramer's rule, with the bound |C^3_1 - C^3_2| ≤ Cε obtained by absorbing a small multiple of itself (Cε ≤ 1/2); this bound is not presupposed. The lower bounds in Theorems 1.4 and 6.13 compare the Navier-Stokes solution with the Stokes solution, using the already-proved upper bound to show the difference is bounded and citing the external Stokes lower bounds [35, Theorem 1.4] and [34, Theorem 1.6]. This is a legitimate comparison argument, not a renaming of the Stokes result. The one caveat found in the derivation chain is not a circularity: Section 2.2 states "we assume, for simplicity, that h1 and h2 are quadratic and symmetric ... h1 = h2 = 1/2 |x'|^2", and the construction of k, v^3_1, pbar^3_1 and the cancellations (4.4)-(4.5) is performed for this exact symmetric gap. Theorem 1.1, however, assumes h1, h2 = κ/2 |x'|^2 + O(|x'|^3) with no exact symmetry or vanishing third-order remainder for general C^3 convex particles, and the paper supplies no perturbation or change-of-variables argument that reduces the stated geometry to the quadratic symmetric case. This is a completeness/generality gap in the proof as written, not an instance of an output being equal to an input by construction.
Assumptions & free parameters
assumptions (4)
- domain assumption Existence of a smooth solution u in W^{1,2} cap C^2(Omega bar), p in C^1(Omega bar) to (1.1)-(1.3)
- domain assumption The gap surfaces have equal leading curvature: h1=h2=kappa/2 |x'|^2 + O(|x'|^3) from (1.6), and the proof reduces further to h1=h2=|x'|^2/2 in Section 2.2
- standard math Prior Stokes estimates from [34,35] (Propositions 2.2-2.4, 2.6 for i=2, and Theorem 1.4 lower bounds) hold
- standard math The mobility/stiffness matrix a11 is positive definite, as used in the proof of Proposition 2.7 from [8]
Cite this review
Pith. "Pith review of Stress concentration between two adjacent rigid particles in Navier-Stokes flow." pith.science (2026). https://pith.science/paper/LRLPHG6X
@misc{pith2026241116413,
author = {Pith},
title = {Pith review of: Stress concentration between two adjacent rigid particles in Navier-Stokes flow},
year = {2026},
howpublished = {\url{https://pith.science/paper/LRLPHG6X}},
note = {Machine review of arXiv:2411.16413}
}
read the original abstract
In this paper we investigate the stress concentration problem that occurs when two convex rigid particles are closely immersed in a fluid flow. The governing equations for the fluid flow are the stationary incompressible Navier-Stokes equations. We establish precise upper bounds for the gradients and second-order derivatives of the fluid velocity as the distance between particles approaches zero, in dimensions two and three. The optimality of these blow-up rates of the gradients is demonstrated by deriving corresponding lower bounds. New difficulties arising from the nonlinear term in the Navier-Stokes equations is overcome. Consequently, the blow up rates of the Cauchy stress are studied as well.
Reference graph
Works this paper leans on
-
[35]
H.G. Li; L.J. Xu, Estimates for stress concentration be tween two adjacent rigid inclusions in Stokes flow. J. Funct. Anal. 286 (2024), no. 6, 110313
work page 2024
- [1]
- [2]
- [3]
- [4]
- [5]
- [6]
- [7]
Show all 44 references
-
[8]
Bao; H.G
J.G. Bao; H.G. Li; Y.Y. Li, Gradient estimates for soluti ons of the Lam´ e system with partially infinite coefficients in dimensions greater than two, Adv. Math. 305 (2017), 298–3 38
2017
-
[9]
Bonnetier, D
E. Bonnetier, D. Manceau, F. Triki, Asymptotic of the vel ocity of a dilute suspension of droplets with interfacial tension, Q. Appl. Math. 71 (1) (2013) 89–117
2013
-
[10]
Bonnetier; F
E. Bonnetier; F. Triki, On the spectrum of the Poincar´ e variational problem for two close-to-touching inclusions in 2D, Arch. Ration. Mech. Anal. 209 (2013), 541–567
2013
-
[11]
touching
E. Bonnetier; M. Vogelius, An elliptic regularity resu lt for a composite medium with “touching” fibers of circular cross-section, SIAM J. Math. Anal. 31 (2000), no. 3, 651–677
2000
-
[12]
Cheng, Calculation of elastohydrodynamic film thi ckness in high speed rolling and sliding contacts, in MTI Technical Report, MIT-67TR24, (1967)
H.S. Cheng, Calculation of elastohydrodynamic film thi ckness in high speed rolling and sliding contacts, in MTI Technical Report, MIT-67TR24, (1967)
1967
-
[13]
Amrouche; N
C. Amrouche; N. Seloula, Lp-theory for vector potentia ls and the Sobolev’s inequalities for vector fields: appli- cation to the Stokes equations with pressure boundary condi tions, Math. Models Methods Appl. Sci. 23 (2013) 37–92
2013
-
[14]
Dong; Y.Y
H.J. Dong; Y.Y. Li; Z.L. Yang, Optimal gradient estimat es of solutions to the insulated conductivity problem in dimension greater than two. arXiv:2110.11313v1 [math.AP] . To appear in J. Eur. Math. Soc, 2022
2022 arXiv
-
[15]
H.J. Dong; D. Kim, W eighted Lq-estimates for stationary Stokes system with partially BMO coefficients, J. Differ. Equ. 264 (7) (2018) 4603–4649
2018
-
[16]
Duerinckx, A
M. Duerinckx, A. Gloria, Effective viscosity of semi-di lute suspensions. arXiv:2301.00165
-
[17]
Galdi, An introduction to the mathematical theory of the Navier-Stokes equations: steady-state problems, Springer, Cham (2011)
G.P. Galdi, An introduction to the mathematical theory of the Navier-Stokes equations: steady-state problems, Springer, Cham (2011)
2011
-
[18]
Gunzburger, H.C
M. Gunzburger, H.C. Lee, G.A. Seregin, Global existenc e of weak solutions for viscous incompressible flows around a moving rigid body in three dimensions, J. Math. Flui d Mech. 2 (3) (2000), 219–266
2000
-
[19]
Galdi, Large-Time Behavior of a Rigid Body o f Arbitrary Shape in a Viscous Fluid Under the Action of Prescribed Forces and Torques, J
Giovanni P. Galdi, Large-Time Behavior of a Rigid Body o f Arbitrary Shape in a Viscous Fluid Under the Action of Prescribed Forces and Torques, J. Math. Fluid Mech. 25 (20 23), no. 3, Paper No. 43, 11 pp
-
[20]
G´ erard-Varet, M
D. G´ erard-Varet, M. Hillairet, Analysis of the viscos ity of dilute suspensions beyond Einstein’s formula, Arch. Ration. Mech. Anal. 238 (3) (2020) 1349–1411
2020
-
[21]
G´ erard-Varet; M
D. G´ erard-Varet; M. Hillairet, Regularity issues in t he problem of fluid structure interaction. Arch. Ration. Mec h. Anal. 195 (2010), 375–407
2010
-
[22]
G´ erard-Varet; M
D. G´ erard-Varet; M. Hillairet, Existence of weak solu tions up to collision for viscous fluid-solid systems with sl ip. Comm. Pure Appl. Math. 67 (2014), 2022–2075
2014
-
[23]
Jeffery, The motion of ellipsoidal particles immer sed in a viscous fluid, Proc
G.B. Jeffery, The motion of ellipsoidal particles immer sed in a viscous fluid, Proc. Roy. Soc. London Ser. A, 102 (1922), 161–179
1922
-
[24]
Kang, Quantitative analysis of field concentration i n presence of closely located inclusions of high contrast, Proceedings of the International Congress of Mathematicia ns 2022
H. Kang, Quantitative analysis of field concentration i n presence of closely located inclusions of high contrast, Proceedings of the International Congress of Mathematicia ns 2022
2022
-
[25]
H. Kang, H. Lee, K. Yun, Optimal estimates and asymptoti cs for the stress concentration between closely located stiff inclusions, Math. Ann. 363 (2015), 1281–1306. STRESS CONCENTRATION PROBLEM IN NA VIER-STOKES FLOW 33
2015
-
[26]
H. Kang; M. Lim; K. Yun, Asymptotics and computation of t he solution to the conductivity equation in the presence of adjacent inclusions with extreme conductiviti es, J. Math. Pures Appl. 99 (2013), 234–249
2013
-
[27]
H. Kang; S. Yu, Quantitative characterization of stres s concentration in the presence of closely spaced hard inclusions in two-dimensional linear elasticity, Arch. Ra tion. Mech. Anal. 232 (2019), 121–196
2019
-
[28]
Hillairet; L
M. Hillairet; L. Sabbagh, Global solutions to coupled ( Navier-)Stokes Newton systems in R3. Asymptotic Analysis 132(2) (2022), 1-27
2022
-
[29]
Hillairet; T
M. Hillairet; T. Takahashi, Collisions in three-dimen sional fluid structure interaction problems, SIAM J. Math. Anal. 40 (6) (2009) 2451–2477
2009
-
[30]
Kirchhoff, ¨Uber die Bewegung eines Rotationsk¨ orpers in einer Fl¨ ussigkeit, J
G. Kirchhoff, ¨Uber die Bewegung eines Rotationsk¨ orpers in einer Fl¨ ussigkeit, J. Reine Angew. Math. 71 (1869), 237-281
-
[31]
Ladyzhenskaya, The mathematical theory of viscou s incompressible flow, Mathematics and its applications, vol 2, 2nd edn
O.A. Ladyzhenskaya, The mathematical theory of viscou s incompressible flow, Mathematics and its applications, vol 2, 2nd edn. Gordon and Breach, New York (1969)
1969
-
[32]
Li, Lower bounds of gradient’s blow-up for the Lam´ e system with partially infinite coefficients, J
H.G. Li, Lower bounds of gradient’s blow-up for the Lam´ e system with partially infinite coefficients, J. Math. Pures Appl. 149 (2021), 98–134
2021
-
[33]
H.G. Li; Y.Y. Li; E.S. Bao; B. Yin, Derivative estimates of solutions of elliptic systems in narrow regions, Quart. Appl. Math. 72 (2014), no. 3, 589–596
2014
-
[34]
H.G. Li; L.J. Xu, Estimates for stress concentration be tween two adjacent rigid inclusions in two-dimensional Stokes flow. arXiv: 2204.00254v2. To appear in Sci. China Mat h
-
[36]
Y.Y. Li; Z.L. Yang, Gradient estimates of solutions to t he insulated conductivity problem in dimension greater than two, Math. Ann. 385 (2023), no. 3-4, 1775–1796
2023
-
[37]
Giaquinta
M. Giaquinta. Multiple integrals in the calculus of var iations and nonlinear elliptic systems. Princeton Univers ity Press: Princeton, NJ, 1983
1983
-
[38]
Oscar, The Reynolds centennial: a brief history of th e theory of hydrodynamic lubrication
P. Oscar, The Reynolds centennial: a brief history of th e theory of hydrodynamic lubrication. Trans. ASME 109, (1987), 2–20
1987
-
[39]
Reynolds, On the theory of lubrication and its applic ation to Mr
O. Reynolds, On the theory of lubrication and its applic ation to Mr. Beauchamp Towers experiments, including an experimental determination of the viscosity of olive oil . Philos. Trans. R. Soc. 177, (1886) 157-234
-
[40]
Sylvain; M
E. Sylvain; M. Debayan; T. Marius, Large time behaviour for the motion of a solid in a viscous incompressible fluid. Math. Ann. 385 (2023), no. 1-2, 631–691
2023
-
[41]
San Mart ´ ın, V
J.A. San Mart ´ ın, V. Starovoitov, M. Tucsnak, Global we ak solutions for the two-dimensional motion of several rigid bodies in an incompressible viscous fluid, Arch. Ratio n. Mech. Anal. 161 (2002), no. 2, 113–147
2002
-
[42]
Stokes, On the effect of internal friction of fluids on t he motion of pendulums, Trans
G. Stokes, On the effect of internal friction of fluids on t he motion of pendulums, Trans. Cambridge Philos. Soc. 9 (1851), 8-85
-
[43]
W einberger, On the steady fall of a body in a Navier- Stokes fluid, Proc
H.F. W einberger, On the steady fall of a body in a Navier- Stokes fluid, Proc. Symp. Pure Math. 23 (1973), 421–440
1973
-
[44]
W einkove, The insulated conductivity problem, effec tive gradient estimates and the maximum principle, Math
B. W einkove, The insulated conductivity problem, effec tive gradient estimates and the maximum principle, Math. Ann. 385 (2023), no. 1-2, 1–16. (H.G. Li) School of Mathematical Sciences, Beijing Normal University, Laboratory of Mathematics and Complex Systems, Ministry of Educ...
2023
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.