{"id":"7c6916be-8c09-4524-b0ab-dbeea990926b","arxiv_id":"2608.03672","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In the critical dilute regime where particle diameter scales as ε^3 and mutual distance as ε, weak limits of steady Navier-Stokes solutions in a perforated pipe satisfy the Brinkman equation with a friction term, for both prescribed flux and prescribed pressure drop, without small-data assumptions.","lead":"This paper proves that fluid moving through a pipe filled with very small particles is described, in the limit of vanishing particle size, by the Brinkman equations with an extra drag force, even for strong flows. The result covers both prescribed total flux and prescribed pressure drop, removing the usual small-data restriction from this homogenization problem.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Prescribed-flux theorem is stated under (2.11)–(2.13) but the proof silently uses the strip condition (4.3) to identify the limit Bernoulli boundary values; Remark 4.4 leaves the case without (4.3) open.","rationale":"Read in good faith, the paper's central claim is that in the critical polydisperse regime, both the prescribed-flux and prescribed-pressure-drop Navier–Stokes problems homogenize to a Brinkman system with no smallness condition. For that to be true, the contradiction argument producing ε-uniform bounds must identify the boundary values of the normalized limit Bernoulli pressure. The most insecure step is the use of assumption (4.3) in Theorem 4.3: without particle-free strips at Γ_I and Γ_O, the W^{1,3/2} estimates (4.33)–(4.34) are not uniform, and the traces used to force bp=0 cannot be extracted from weak L^2 convergence. The authors themselves flag this in Remark 4.4, so my concern is not a speculative counterexample but an admitted unsupported case. Because the theorem statement lists only (2.11)–(2.13), the claim as written exceeds what is proved. This does not invalidate the substantial parts: Theorem 3.6 (restriction operator) and the proof of Theorem 4.3 are detailed, and the pressure-drop uniform bounds in Theorem 4.5 are given. The omitted proof of Theorem 2.5 is a second completeness gap but is less load-bearing because no hidden hypothesis is involved there. A conditional verdict is therefore appropriate, with the condition being to either add (4.3) to Theorem 2.4's hypotheses or supply the missing boundary argument.","tokens_in":46186,"tokens_out":6652,"duration_ms":59911,"concrete_test":"Test the boundary-identification step in a minimal one-particle geometry: let Ω be a straight pipe with inlet Γ_I, place one particle of radius ε^3 at distance ε from Γ_I (so (4.3) fails for any δ3>ε), and solve the normalized Stokes system used in the contraction argument with right-hand side f ≡ 0 and no flux. Compute the trace of the normalized Bernoulli pressure Φ^ε/Z_ε^2 on Γ_I and the W^{1,3/2}(Ω_I)-norm of Φ^ε/Z_ε^2 as ε→0. If the trace does not tend to 0, or the norm grows like ε^{-1/2} as Remark 4.4 suggests, then the inference bp=0 in Theorem 4.3 fails and Theorem 2.4 genuinely requires (4.3); if the trace still vanishes, the assumption may be an artifact of the proof and a sharpened argument could remove it.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 2.4 is the paper's first main homogenization result and is advertised, like the abstract, as holding under only (2.11)–(2.13). In the proof, however, the uniform bounds come from Theorem 4.3, and Theorem 4.3 has an additional hypothesis: no particles in the inlet/outlet strips (4.3). The contradiction argument in Theorem 4.3 needs (4.3) at exactly one point: after normalizing by Z_ε, the estimates (4.33)–(4.34) give uniform W^{1,3/2} bounds for u^ε/Z_ε^2 and Φ^ε/Z_ε^2 in Ω_I and Ω_O only because these strips contain no particles. These bounds produce the traces that force the limit Bernoulli pressure to vanish on Γ_I and Γ_O, yielding bp=0 and contradicting (4.21). If particles are allowed in the strips, the local-regularity constant in (4.33)–(4.34) becomes ε-dependent; Remark 4.4 records only a whole-domain W^{1,3/2} estimate with a 1/√ε factor. Thus the boundary identification step is unsupported exactly in the configuration the theorem statement claims to cover. Since Remark 4.4 explicitly leaves open the removal of (4.3), this is not a cosmetic gap: either (4.3) must be added to Theorem 2.4's hypotheses, or a new argument must identify the traces without uniform strip regularity. (Separately, Theorem 2.5 is stated without proof, but the flux-case gap is the load-bearing one.)","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies homogenization of the steady-state Navier-Stokes equations in a three-dimensional pipe perforated by many small particles of diameter ε^3 at mutual distance ε, in the critical regime where a Brinkman term is expected. Two boundary-value problems are considered: one with prescribed transversal flux rate and one with prescribed pressure drop, both with mixed boundary conditions involving the tangential velocity and the Bernoulli pressure on the inlet and outlet. The main results, Theorems 2.4 and 2.5, assert that any weak accumulation point of the ε-level solutions satisfies a Brinkman-type system with new terms Gu and J, without smallness assumptions on the data. The uniform bounds for the flux problem are obtained by a contradiction argument using Bernoulli's law for stationary Euler solutions, while the pressure-drop problem is handled by a direct energy estimate. The paper also constructs a flux carrier, a divergence-free extension of particle velocities, Stokes correctors, and a restriction operator with sharp estimates, with a detailed appendix on exterior Stokes problems and capacity densities.","tokens_in":46448,"tokens_out":10728,"duration_ms":98614,"significance":"If the main theorems hold as stated, this is a significant extension of the subcritical homogenization result in [59] to the critical particle-size regime, with non-periodic polydisperse particle configurations, non-zero constant particle velocities, and no smallness condition on the data. The paper contains genuinely useful technical contributions: explicit estimates for Stokes correctors, a genuine restriction operator that does not modify test functions vanishing on the particles, and a self-contained treatment of the exterior Stokes resistance matrix. The proof of Theorem 2.4 is largely detailed and the auxiliary estimates are explicit. However, the flux-case result currently depends on an additional strip condition (4.3) that is not stated in Theorem 2.4, and the second central claim, Theorem 2.5, is asserted without proof. These issues are fixable, but they are load-bearing for the advertised conclusions.","major_comments":[{"comment":"Theorem 2.4 is stated under assumptions (2.11)-(2.13), but its proof relies on the uniform bounds of Theorem 4.3, which are proved only under the additional strip assumption (4.3). The trace identification of the limit Bernoulli pressure on Γ_I and Γ_O in the contradiction argument uses the W^{1,3/2} regularity estimates (4.33)-(4.34), and those estimates are valid only because assumption (4.3) keeps the strips Ω_I and Ω_O free of particles. Remark 4.4 explicitly leaves open the possibility of recovering Theorem 4.3 without (4.3), so this is not a cosmetic gap: either (4.3) must be added to the hypotheses of Theorem 2.4 and to the abstract's phrasing, or a new argument must identify the boundary traces without uniform strip regularity.","section":"Theorem 2.4 and Remark 4.4"},{"comment":"Theorem 2.5 is one of the two central homogenization results, but its proof is omitted with the sentence that it resembles the proof of Theorem 2.4 with minor modifications. Given that Theorem 2.5 involves a different functional space, a different uniform-bounds theorem (Theorem 4.5), and an additional pressure boundary term, the proof is not literally identical. Please provide a complete proof or a detailed outline that identifies the test space, the treatment of the pressure-boundary term, and the changes in the passage to the limit.","section":"Section 4.2, Theorem 2.5"},{"comment":"Definitions 4.1 and 4.2 are internally inconsistent as written. In Definition 4.1, the condition u ∈ V(Ωε) forces u = 0 on each ∂K_n^ε, while the condition u − Ψε ∈ U(Ωε) forces u − Ψε = 0 on ∂K_n^ε and hence u = Ψε = μ_n^ε on ∂K_n^ε. Unless all μ_n^ε vanish, no such u exists. The same inconsistency appears in Definition 4.2 with Aε. The intended weak formulation should be stated for the affine space containing the flux carrier, e.g. u ∈ Ψε + U(Ωε), or the space V(Ωε) should be defined without the condition u = 0 on ∂K_n^ε. Since the subsequent proofs use the strong form (2.8)-(2.9), the intended meaning is recoverable, but the formal definitions need to be corrected.","section":"Section 4.1, Definitions 4.1-4.2"}],"minor_comments":[{"comment":"The sign of the nonlinear term in (4.19) appears to be wrong: testing with uε − Aε gives −∫((∇×buε)×buε)·Aε, not +∫((∇×buε)×buε)·Aε. The term vanishes in the limit by (4.20), so the final conclusion is unaffected, but the displayed identity should be corrected.","section":"Equation (4.19)"},{"comment":"There are small typographical errors: in Remark 3.5, 'this han been proved' should be 'this has been proved', and in Proposition A.4 the name 'Cauchy-Schwartz' should be 'Cauchy-Schwarz'.","section":"Remark 3.5 and Proposition A.4"},{"comment":"The notation p^∓_ε in (2.5) is confusing because the normalization p^-_ε = 0 is introduced only later. Please define the constants p^-_ε and p^+_ε at the point of first use and state the normalization explicitly in the displayed system.","section":"Section 2.3, notation p^∓_ε"},{"comment":"The sequence bvε = euε/Z^2_ε is introduced in (4.35) without a fresh notation paragraph, after buε has already been defined as euε/Zε. This makes the proof harder to follow; a short note that bvε is a further rescaling of the same extended velocity would improve readability.","section":"Section 4.1, around (4.35)"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope and the core inhomogeneous-homogenization ideas appear sound, but the flux theorem's dependence on the unstated strip condition (4.3), the missing proof of Theorem 2.5, and the inconsistent weak formulations in Definitions 4.1-4.2 need to be resolved. I do not see evidence of any fatal flaw that would require rejection; the issues are substantial but addressable within the manuscript's framework."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this paper genuinely extends the earlier subcritical result [59] to the critical scaling α=3, deriving a Brinkman term for steady Navier-Stokes in perforated pipes with prescribed flux or pressure drop, polydisperse particles, nonzero particle velocities, and no small-data restriction. That is a meaningful step, and the main flux-case proof is mostly convincing.\n\nWhat is actually new: the restriction operator built from Stokes correctors (Theorem 3.6) is a solid piece of work, and the contradiction argument for ε-uniform bounds via Bernoulli's law is clever and written out in real detail. I saw no fitting or circularity. The Brinkman force Gu and source term J come from exterior Stokes resistance matrices under an explicit convergence assumption (2.13), and the self-citations [59,60] are independent published results. Credit is earned there.\n\nThe soft spots are real, though. Theorem 2.4 is advertised under assumptions (2.11)–(2.13), but its proof relies on Theorem 4.3, which adds condition (4.3): no particles in the inlet and outlet strips. The stress-test note is right about where this bites. The uniform W^{1,3/2} estimates in (4.33)–(4.34) hold only because those strips are particle-free, and those estimates are what identify the boundary traces forcing bp=0. Without (4.3), the local-regularity constant becomes ε-dependent, and Remark 4.4 explicitly leaves removal of (4.3) open. So the theorem statement overclaims. The fix is either to add (4.3) to Theorem 2.4's hypotheses or to find a new boundary-identification argument. The latter may be possible, but it is not in the paper.\n\nSeparately, Theorem 2.5 is stated as a main result and then left unproved, “for the sake of brevity.” Given that the flux case is the delicate one and that Theorem 4.5 suggests the pressure-drop case avoids (4.3), a referee should ask for the proof rather than take it on faith. This is not a cosmetic omission; it is half the advertised contribution.\n\nWho gets value from this paper: people working on homogenization of fluid equations in perforated domains, especially those interested in Brinkman law, Leray-type flux problems, and mixed pressure boundary conditions. The detailed parts deserve a serious referee, and the gaps are addressable in revision. My recommendation: send it to peer review, but ask the authors to fix the mismatch between Theorem 2.4 and condition (4.3), and to either provide the proof of Theorem 2.5 or state clearly why the pressure-drop case avoids the strip assumption.","headline":"Solid homogenization result with a real incompleteness: the flux-case theorem as stated uses an extra strip condition that is not listed in its hypotheses, and the pressure-drop theorem is stated without proof.","tokens_in":47037,"tokens_out":1666,"would_cite":true,"duration_ms":17496,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76M50","76D05","35B27","35M32","35Q31"],"pacs":[],"model":"deepseek-v4-flash","headline":"In a pipe full of tiny particles, every weak limit of steady Navier-Stokes solutions obeys a Brinkman-type equation with an added friction term, with no smallness condition on the data.","keywords":["incompressible fluids","mixed boundary conditions","homogenization","perforated domain","Brinkman equation","Navier-Stokes equations","prescribed flux","Bernoulli pressure"],"falsifier":"Take the same pipe, same separations, and the same bound (2.11), and admit one particle whose $\\varepsilon$-scale ball intersects the $\\delta_3$-strip adjacent to $\\Gamma_I$ while all other assumptions hold; then check whether $\\|\\nabla u^\\varepsilon\\|_{L^2}$ stays bounded and whether the limits still satisfy (2.15). A bounded limit satisfying (2.15) would show (4.3) is removable; a blow-up, or a limit in which the inlet Bernoulli-pressure condition fails, would show (4.3) is genuinely load-bearing. The paper's Remark 4.4 already identifies this exact case as open.","tokens_in":45918,"feed_emoji":"🌊","tokens_out":7748,"duration_ms":69616,"temperature":0.7,"pith_summary":"The paper proves that a viscous incompressible fluid flowing through a three-dimensional pipe riddled with many tiny particles behaves, in the macroscopic limit, like a fluid governed by a Brinkman-type equation. The particles have diameter of order $\\varepsilon^3$ and mutual distance of order $\\varepsilon$, the critical homogenization scale: small enough to be dilute, but numerous enough to exert a collective friction. The setting is the steady Navier-Stokes system with mixed boundary conditions, involving tangential velocity and Bernoulli pressure on the inlet and outlet, with either the total flux or the pressure drop prescribed. The authors show that every weak accumulation point of the $\\varepsilon$-level solutions satisfies the homogenized system with an added drag term $Gu-J$, for arbitrarily large forces, flux rates, and pressure drops, with no smallness assumption. The main technical achievement is the corresponding $\\varepsilon$-uniform bound, obtained for the flux problem through a contradiction argument based on Bernoulli's law for stationary Euler flows.","feed_headline":"Tiny obstacles turn steady pipe flow into Brinkman flow","feed_subtitle":"For prescribed flux or pressure drop, the homogenized limit gains a drag term with no smallness condition.","key_machinery":"The argument is carried by a family of Stokes correctors and a restriction operator $R_\\varepsilon(\\varphi)$ that modifies a test vector field in thin shells around each particle so that it becomes admissible in the perforated domain and divergence-free, without changing fields that already vanish on the particles. The key identity (3.40) shows that the limit of $\\int_\\Omega \\nabla v^\\varepsilon : \\nabla R_\\varepsilon(\\varphi)$ contains the extra term $\\int_\\Omega (G v - J)\\cdot\\varphi$, which is exactly the Brinkman force. Supporting objects are the relative capacity potential (uniformly bounded in the critical regime), a uniform inverse divergence operator on perforated domains, a flux carrier built from Hagen-Poiseuille-type flows in the cylindrical outlets, and, for the prescribed-flux uniform bound, a contradiction argument that rescales a hypothetically unbounded solution and invokes Bernoulli's law for the limiting stationary Euler equation to show the scaled pressure must vanish on the walls.","core_discovery":"The central discovery is that in the critical regime (particle diameter $\\varepsilon^3$, spacing $\\varepsilon$), the steady Navier-Stokes equations with prescribed flux rate or prescribed pressure drop homogenize to the Brinkman-type system $$-\\$\\Delta$ u + (\\nabla\\times u)\\times u + \\nabla\\Phi + G u - J = f,\\qquad \\operatorname{div} u=0,$$ in the original unperforated admissible pipe domain, with the same mixed boundary conditions on the inlet, outlet, and walls, and with the same flux or pressure-drop prescription (Theorems 2.4 and 2.5). Here $G\\in L^2(\\Omega;\\mathbb{R}^{3\\times 3})$ is the limiting density of Stokes resistance matrices of the particles and $J\\in L^2(\\Omega;\\mathbb{R}^3)$ encodes their prescribed velocities; both arise as limits of $\\varepsilon$-scaled sums over the particles, as in (2.12)-(2.14). The result holds without any smallness condition on the data, and the pressure-drop version needs only the kinetic-energy bound (2.11), while the flux version additionally assumes the inlet and outlet strips stay particle-free (condition (4.3)).","pith_inferences":["If condition (4.3) turns out to be removable, the natural route suggested by the paper's structure is to extend the Bernoulli pressure inside the holes using the particle velocities rather than by zero; the pressure-drop proof already shows that boundary-adapted extensions can work.","For periodic or random stationary particle configurations the limit matrix $G$ should reduce to a multiple of the identity proportional to particle number density times mean resistance, recovering the classical scalar Brinkman permeability law; this is a testable specialization the paper does not spell out.","Because the uniform bounds hold without smallness, the same restriction operator is a plausible tool for an evolutionary Navier-Stokes version with prescribed flux, where the particle velocities would be coupled to their motion; the present paper treats only the stationary case.","The pressure-drop result suggests that prescribing the drop is the more robust formulation: it avoids the delicate boundary-layer condition (4.3), so applied settings should prefer the drop formulation when inlet and outlet obstructions cannot be excluded."],"forward_implications":["In the critical regime the particles neither disappear nor clog the pipe: their collective response survives as the zero-order drag term $Gu-J$, for arbitrary sizes of the data.","The prescribed-pressure-drop limit needs no strip condition near the inlet or outlet; only the uniform kinetic-energy bound (2.11) on the particle velocities is required.","The limiting flux problem keeps the prescribed flux condition $\\int_{\\Gamma_O} u\\cdot\\nu = F$ and the same Bernoulli-pressure boundary conditions, with the unknown outlet constant $p^+$ emerging as part of the limit solution.","Uniform bounds produce at least one weak accumulation point in both problems, so the effective Brinkman systems (2.15)-(2.16) admit weak solutions in this setting without smallness of the data.","The same scheme covers distorted pipes with cylindrical outlets and non-zero constant particle velocities, and it extends to several inlets and outlets."],"supporting_citations":[{"why":"Establishes the critical-regime homogenization framework for Navier-Stokes in perforated domains, including the Stokes corrector approach.","marker":"[2]"},{"why":"Introduced the Brinkman equation, the target effective model with an added friction force.","marker":"[11]"},{"why":"Provides the uniform inverse of the divergence operator on perforated domains, used for the pressure estimates and the Stokes problems at the $\\varepsilon$-level.","marker":"[19]"},{"why":"Develops the restriction-operator technique for Stokes homogenization, which the paper adapts and extends to polydisperse settings with nonzero particle velocities.","marker":"[28]"},{"why":"Supplies the kinetic-energy assumption on particle velocities and the framework for homogenization around arrays of moving spheres.","marker":"[30]"},{"why":"Provides the polydisperse Stokes-Brinkman force with non-zero particle velocities, plus the exterior Stokes estimates and resistance-matrix bounds used here.","marker":"[32]"},{"why":"Gives the weak formulation with Bernoulli pressure and pressure boundary conditions, which underpins the definitions of both problems.","marker":"[44]"},{"why":"Is the subcritical predecessor for homogenization of the flux or pressure-drop problem in a perforated pipe, which this paper extends to the critical regime.","marker":"[59]"},{"why":"Supplies existence and regularity for the steady Navier-Stokes problem in an obstructed pipe under mixed boundary conditions with prescribed flux.","marker":"[60]"}],"fun_headline_variants":["Pipe flow with tiny particles gains a Brinkman drag term","No smallness needed: Navier-Stokes homogenizes to Brinkman in pipes","Critical particle size yields Brinkman correction in steady pipe flow","From Navier-Stokes to Brinkman: drag emerges from sparse particles","Prescribed flux or pressure drop: Brinkman term emerges without smallness"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"For the prescribed-flux theorem, the load-bearing premise is assumption (4.3): a strip of fixed width $\\delta_3$ next to the inlet and outlet must never contain any particle; without it the boundary values of the limiting Bernoulli pressure cannot be identified and the uniform-bound proof stops.","fun_headline_variants_meta":{"raw":{"variants":["Pipe flow with tiny particles gains a Brinkman drag term","No smallness needed: Navier-Stokes homogenizes to Brinkman in pipes","Critical particle size yields Brinkman correction in steady pipe flow","From Navier-Stokes to Brinkman: drag emerges from sparse particles","Prescribed flux or pressure drop: Brinkman term emerges without smallness"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000779,"raw_usage":{"total_tokens":3453,"prompt_tokens":962,"completion_tokens":2491,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":578,"completion_tokens_details":{"reasoning_tokens":2399}},"tokens_in":578,"tokens_out":2491,"duration_ms":14685,"temperature":1.0,"reasoning_tokens":2399,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:47:18.964275+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same pipe, same separations, and the same bound (2.11), and admit one particle whose $\\varepsilon$-scale ball intersects the $\\delta_3$-strip adjacent to $\\Gamma_I$ while all other assumptions hold; then check whether $\\|\\nabla u^\\varepsilon\\|_{L^2}$ stays bounded and whether the limits still satisfy (2.15). A bounded limit satisfying (2.15) would show (4.3) is removable; a blow-up, or a limit in which the inlet Bernoulli-pressure condition fails, would show (4.3) is genuinely load-bearing. The paper's Remark 4.4 already identifies this exact case as open.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the critical-regime homogenization framework for Navier-Stokes in perforated domains, including the Stokes corrector approach."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduced the Brinkman equation, the target effective model with an added friction force."},{"cited_title":"Diening, E","cited_arxiv_id":null,"evidence_quote":"Provides the uniform inverse of the divergence operator on perforated domains, used for the pressure estimates and the Stokes problems at the $\\varepsilon$-level."},{"cited_title":"Giunti and R","cited_arxiv_id":null,"evidence_quote":"Develops the restriction-operator technique for Stokes homogenization, which the paper adapts and extends to polydisperse settings with nonzero particle velocities."},{"cited_title":"Hillairet","cited_arxiv_id":null,"evidence_quote":"Supplies the kinetic-energy assumption on particle velocities and the framework for homogenization around arrays of moving spheres."},{"cited_title":"Hillairet, A","cited_arxiv_id":null,"evidence_quote":"Provides the polydisperse Stokes-Brinkman force with non-zero particle velocities, plus the exterior Stokes estimates and resistance-matrix bounds used here."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Is the subcritical predecessor for homogenization of the flux or pressure-drop problem in a perforated pipe, which this paper extends to the critical regime."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies existence and regularity for the steady Navier-Stokes problem in an obstructed pipe under mixed boundary conditions with prescribed flux."}],"review_version":2}