{"id":"19429743-d935-4401-b04a-c003b9f9ee26","arxiv_id":"2607.16470","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A cluster-based method improves the allowable number of tiny rigid balls in a viscous fluid before their collective effect disappears, and shows the balls follow the fluid flow in the limit, even under gravity.","lead":"Researchers study a cloud of many tiny rigid balls in a viscous fluid as the balls shrink and their number grows. They introduce 'clusters' of nearby balls to prove that, under explicit size-count conditions, the cloud leaves the limiting Navier–Stokes equations unchanged and individual balls travel with the fluid even under gravity.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 2.3's trajectory convergence rests on the unproved estimate (7.13) and the imported rigid-velocity bound (6.6); if their r-power or cluster-dependence differs, the convergence h_n,ε → h_n in W^{1,2}(0,T) is unsupported.","rationale":"The paper is coherent and the cluster construction is genuinely new: Lemma 5.1 correctly handles nearby centers, and the estimates in Proposition 5.2 are plausible and internally consistent. No internal contradiction is visible. The weakest point is exactly the imported rigid-velocity bound (6.6), which the reader identified; I add that the final, most novel trajectory claim (2.19) depends on the unstated descendant (7.13), asserted only via ‘similarly to (6.6)’. Since Theorem 2.3 is the advertised first rate-of-convergence result and since the exponents in (2.17) are finely tuned to the r-powers in (6.6), a wrong r-power or a hidden dependence on cluster density would invalidate the convergence h_n,ε → h_n in W^{1,2}(0,T). This is a load-bearing concern that deserves a concrete proof, not just a citation, because collisions are allowed and the constants for perforated domains with touching holes are subtle. I therefore recommend conditional acceptance: the central ideas are sound, but Theorem 2.3 should be accepted only after (6.6) and (7.13) are fully derived with explicit constants, or the hypotheses adjusted accordingly.","tokens_in":25267,"tokens_out":49334,"duration_ms":429083,"concrete_test":"Derive (7.13) from (6.6) by applying the latter to v = u_ε − u_app, which is rigid on each ball, and combine with Lemma 3.1 and the relative-energy inequality (7.6). Compute the resulting power of r and track every constant through the local trace inequality on the shell B(βr,h_n) \\ B(r,h_n). In particular, check whether the constant remains bounded as the minimum inter-ball distance is exactly 2r (touching balls) and as the number of balls inside a ball of radius 4r grows without bound. A decisive analytical check: in d = 3, take two touching balls, fix v = a + ω×x on ball 1 and v = 0 on ball 2 and outside, and compute the optimal ratio |a|^2 r / ∫|∇v|^2 dx; if this ratio is unbounded as r → 0, the uniform r-power in (6.6)/(7.13) fails in the collision regime and Theorem 2.3 is not established.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central new claim, Theorem 2.3, is the first rate/convergence result for trajectories of the rigid balls in the dynamic homogenization limit. Its final step is the estimate (7.13):\n\nr^α |dh_{n,ε}/dt − u(t, X_{B,4r_ε}(h_{n,ε}))|^2 ≲ ∫_T^d (S(Du_ε) − S(Du_app)) : (D u_ε − D u_app) dx.\n\nThis is asserted with the one-line justification “similarly to (6.6)”, but it is not proved in the paper. The bound (6.6) itself is imported from the authors’ earlier work [12] and is the exact place where the r-powers in the hypotheses (2.17) are calibrated: the exponents in (2.17) are chosen so that the r^{-1/p} factor in (6.6) is compensated by the factors N^a r^b that appear in the estimates (5.18), (5.20), (7.7), (7.11). If (6.6) has a different r-power, e.g. r^{-β} with β ≠ 1/p, or if its constant depends on the local cluster geometry (for instance on the number of touching balls in a 4r-neighborhood), then the terms in Step 1, Step 6, and the final trace estimate (7.13) would not converge under the stated conditions (2.17). In particular, because the paper explicitly allows collisions (Section 1.1), the trace/Poincaré constant for the fluid domain with many touching small holes must be shown to be uniform in the number of bodies in a cluster; this is exactly the situation the cluster method is intended to control, but the proof does not supply the needed constant. The reader’s weakest-assumption flag is therefore on target, and the specific gap is that the rigorous derivation of (7.13) from (6.6) and the relative-energy inequality is omitted.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops a new 'cluster' approach to the dynamic homogenization of a cloud of N small rigid balls of common radius r moving in a viscous incompressible fluid. The balls may collide, and their number grows as r→0. The main results are Theorem 2.2, showing that under the growth conditions (2.12) the cloud has no effect on the limit, which is the incompressible Navier–Stokes system, and Theorem 2.3, showing that if the limit solution is smooth, the ball trajectories converge, h_{n,ε}→h_n in W^{1,2}(0,T), with dh_n/dt = u(t,h_n(t)). The central technical novelty is a cluster projection X_{B,δ}, introduced in §4.3.1, and the associated approximate test functions constructed in §5. The proof of Theorem 2.2 is completed via the estimates of §5–6; Theorem 2.3 is proved in §7 by a relative-energy inequality.","tokens_in":25768,"tokens_out":26960,"duration_ms":222183,"significance":"If correct, the paper substantially improves the previously known logarithmic critical number N≈log(1/r) from [12] to a polynomial-in-1/r bound, and it provides, to the authors' knowledge, the first convergence/rate result for the trajectories of the rigid balls in the dynamic homogenization limit. The cluster-projection construction is an elegant and potentially widely applicable idea, and Lemma 5.1, showing that nearby centers share the same cluster projection, is a clean device. The results are unconditional and global in time, and the paper explicitly allows collisions. The proof is mostly coherent, with explicit rate conditions and no free parameters in the main statements. However, two load-bearing trace estimates are not adequately justified, and the paper is not yet acceptable in its present form.","major_comments":[{"comment":"The final and essential estimate for Theorem 2.3 is asserted with the one-line justification 'similarly to (6.6)', but no proof is given. Equation (7.13) is the only step that transfers the relative-energy dissipation to the individual ball velocities; it must hold with a constant independent of N and, importantly, independent of the local cluster geometry, including touching balls. This is precisely the setting the cluster method is designed to control, so the missing proof is not a cosmetic gap. A trace/capacity lemma should be stated and proved: for each n, r^α |dh_{n,ε}/dt − u(t,X_{B,4r_ε}(h_{n,ε}))|^2 ≲ ∫ (S(Du_ε)−S(Du_app)):(Du_ε−Du_app) dx, with the r-power α matching (2.17) and no N-dependence. Without such a lemma, the convergence h_{n,ε}→h_n in W^{1,2}(0,T) in Theorem 2.3 is unsupported.","section":"§7, Eq. (7.13)"},{"comment":"The uniform rigid-velocity bound |dh_{n,ε}(t)/dt| ≤ c(p)||u_ε(t)||_{W^{1,2}} r^{-1/p} is imported from [12, Sections 3.1–3.2]. This bound is load-bearing: the exponent r^{-1/p} calibrates the rate conditions in (2.17) through the estimates (5.18), (5.20), (7.7), (7.11), and ultimately (7.13). Since collisions are allowed, one must be certain that the constant in (6.6) is uniform in the number of balls in a cluster and in the contact geometry. The paper should either state the precise theorem from [12] with its hypotheses or provide a self-contained proof. As written, a reader cannot verify that the exponents in (2.17) are the correct ones.","section":"§6.1, Eq. (6.6)"}],"minor_comments":[{"comment":"The convergence of the density-error term involving ϱ_ε^2 ||∂_t u_app||_{L^q(∪B)}^2 is not explicitly checked. It follows from (2.17) together with the uniform L^{3/2} bound on density and N r^3→0, but the paper should state this verification, since the displayed estimate (7.11) alone does not make the convergence obvious.","section":"§7, Step 6"},{"comment":"The proof of the time-differentiability estimate (4.27) is only sketched. Since this estimate feeds directly into Proposition 5.2, a short derivation using (4.23)–(4.26) and (4.10) would improve readability.","section":"§4.3.1, Lemma 4.5"},{"comment":"There are several typographical issues: in the Abstract, 'cluster- a collection' lacks a space; in Theorem 2.2, (2.13) states W^{1,2}(Ω;R^3) although d=2,3; and in §6.5, equation (6.16) contains 'dxdtdt'.","section":"Throughout"}],"recommendation":"major_revision","confidential_remarks":"The paper is potentially a strong contribution to dynamic homogenization, and the cluster idea is fresh. My main concern, as stated in the major comments, is that the proof of the critical trace estimate (7.13) is missing, while the imported bound (6.6) is too important to be taken on faith. These are fixable in principle, but the authors must supply a rigorous, self-contained derivation of (7.13) with explicit uniformity in the cluster geometry and verify that the r-powers in (2.17) are consistent with all error terms. If those gaps are filled, the paper would be a solid accept."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the cluster machinery: cluster distance, cluster projection, and the approximate test functions built on them. That is a real technical idea, and it earns the main improvement — moving the critical number of balls from logarithmic in [12] to polynomial in 1/r, and producing the first trajectory-level convergence rate in this dynamic homogenization problem. The proof of Theorem 2.2 is coherent, and the relative-energy argument in Section 7 is structurally sound. I found no circularity. The cluster projection estimates, especially (4.22), and Lemma 5.1 are the right load-bearing pieces, and the authors use them honestly.\n\nThe soft spots are concentrated in Theorem 2.3. The final step, estimate (7.13), is asserted with the one-line justification \"similarly to (6.6)\" and is not proved in the paper. That matters: (7.13) is what converts the relative-energy dissipation into a statement about the ball trajectories, and its r-power is matched against the hypotheses in (2.17). If the constants or exponents differ — for instance if the trace/Poincaré constant is not uniform when balls touch inside a cluster — the trajectory convergence does not follow. The imported bound (6.6) from [12] is the other load-bearing assumption. It may well be correct, and citing a prior paper for it is legitimate, but here it is used in several places where the cluster geometry is supposed to be controlled, and the paper does not state the precise conditions under which (6.6) holds with touching balls. This is fixable: the authors should either prove (7.13) directly or state (6.6) as a lemma with a proof sketch and explicit dependence on the cluster geometry.\n\nMinor point: the sedimentation discussion in Section 2.4 is more heuristic than the rest; it is not needed for the main theorems, so I would not treat it as a serious flaw.\n\nThis paper is worth a serious referee. It is not a desk reject. The right outcome is probably major revision: the authors need to close the gap at (7.13) and make the imported estimates precise, but the core cluster idea is strong and the main claims are likely true.","headline":"The cluster method is real and the polynomial improvement is credible, but Theorem 2.3 leans on an unproved estimate at (7.13) and an imported r-power bound that the paper should re-derive or pin down precisely.","tokens_in":26232,"tokens_out":1494,"would_cite":true,"duration_ms":17868,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q30","76D05","74F10","76M50"],"pacs":[],"model":"deepseek-v4-flash","headline":"Clouds of up to polynomial-in-1/r many tiny rigid balls leave the Navier–Stokes limit unchanged.","keywords":["fluid–structure interaction","dynamic homogenization","rigid balls in viscous fluid","Navier–Stokes system","cluster decomposition","relative energy","small rigid bodies","weak solutions"],"falsifier":"Run a two-dimensional periodic Navier–Stokes simulation with N ≈ r^{-2/5} small balls at the critical scaling and compare the time-averaged velocity with the unperforated solution; if the difference does not vanish as r→0 at the predicted rate, the cluster bound is violated. Alternatively, check the imported rigid-velocity bound directly for a pair of nearby balls: if the maximum speed grows faster than r^{-1/p}, the rate estimates break.","tokens_in":25159,"feed_emoji":"🌊","tokens_out":4354,"duration_ms":37854,"temperature":0.7,"pith_summary":"The paper studies a viscous fluid in a periodic domain containing N identical small rigid balls of radius r, with N growing and r shrinking. It aims to show that, under a growth condition on N of the form N ≈ r^{-β}, the collective effect of the balls disappears and the limiting velocity obeys the incompressible Navier–Stokes equations. It further aims to show that when the limit solution is smooth, each ball's center converges to a trajectory following the fluid velocity, even under gravity; this gives, the authors say, the first rate of convergence in this dynamic homogenization limit. The new idea is to group balls into 'clusters' whose center of mass acts like a single rigid body, so nearby balls do not force separate boundary conditions.","feed_headline":"Tiny ball clouds vanish from the Navier–Stokes limit","feed_subtitle":"A new cluster method pushes the allowed number of small rigid bodies far beyond previous logarithmic bounds, and the balls follow the flow.","key_machinery":"The conceptual engine is the cluster: a partition of the set of ball centers into connected components at scale δ, with a cluster distance d_K[x,y] measuring how far a curve must stay from the set K to join x and y. The cluster projection X_{B,δ}(h_i) is a weighted barycenter of centers whose cluster distance from h_i is less than δ, with exponentially decaying weights. Because close centers have the same projection, replacing each center by its cluster projection in the construction of approximate test functions makes the test function rigid on whole clusters rather than on each ball separately. These approximate test functions, with the exponential cutoff of the projection, supply the erro","core_discovery":"The central claim is that the cloud of balls is asymptotically invisible: the fluid velocity converges to the incompressible Navier–Stokes solution, and the ball trajectories converge to integral curves of that velocity field, d/dt h_n = u(t,h_n). The critical number N(r) for which this holds is improved from a logarithmic bound to a polynomial bound in 1/r. The proof constructs approximate divergence-free test functions that are rigid on the balls by replacing each ball's center by its cluster projection; nearby centers merge into one projection, so the boundary condition acts as a single rigid body on each cluster.","pith_inferences":["The proof's only r-power input from prior work is the rigid-velocity bound, so the same cluster construction should transfer to other fluid models (e.g., slightly compressible or non-Newtonian) with the critical N changed only by the available energy estimates — a testable cross-model prediction.","One could test numerically whether the invisible-cloud threshold is sharp: for N slightly above the polynomial bound, a nonzero correction to the Navier–Stokes limit should appear at order N r^{d/5}|log r|^{4/5}.","The cluster distance is defined on centers only; extending the argument to ellipsoidal or non-identical rigid bodies would require a cluster distance on orientations, but the merge-on-contact mechanism should survive.","The result suggests a coarse-graining principle for dilute suspensions: when small particles are free to move, their local velocity defects cancel in the bulk, a principle that the relative-energy method may quantify."],"forward_implications":["The Navier–Stokes equations are the correct effective equations for the fluid even with N up to about r^{-d/5} (up to logarithmic factors) small balls, and the same holds in two dimensions.","Individual ball trajectories are asymptotically material: each ball follows the fluid flow, with h_n converging in W^{1,2}(0,T) to solutions of d/dt h_n = u(t,h_n), provided the limit solution is smooth.","The result holds for global-in-time weak solutions and permits collisions; balls may cluster, touch, or stay together without changing the limit.","The balls' densities may grow moderately without altering the limit, so gravity does not prevent the balls from following the flow in this regime.","For a periodic sedimenting cloud, the vertical position of every ball tends to rest at the bottom (or top) of the slab, with kinetic energy decaying to zero."],"fun_headline_variants":["Cluster method makes tiny balls vanish from fluid equations","Invisible ball clouds: polynomial bound achieved","Small rigid balls don't disturb Navier-Stokes limit","Clusters act as single body in viscous fluid flow","Tiny ball swarm follows fluid in asymptotic limit"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The proof leans on the imported bound |dh_{n,ε}/dt| ≤ c(p) ||u_ε||_{W^{1,2}} r^{-1/p} from a preceding paper; if that bound is wrong or has a different power of r, the derived convergence rates and hence the main theorems collapse.","fun_headline_variants_meta":{"raw":{"variants":["Cluster method makes tiny balls vanish from fluid equations","Invisible ball clouds: polynomial bound achieved","Small rigid balls don't disturb Navier-Stokes limit","Clusters act as single body in viscous fluid flow","Tiny ball swarm follows fluid in asymptotic limit"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000251,"raw_usage":{"total_tokens":1337,"prompt_tokens":631,"completion_tokens":706,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":375,"completion_tokens_details":{"reasoning_tokens":633}},"tokens_in":375,"tokens_out":706,"duration_ms":7145,"temperature":1.0,"reasoning_tokens":633,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T20:52:07.429029+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a two-dimensional periodic Navier–Stokes simulation with N ≈ r^{-2/5} small balls at the critical scaling and compare the time-averaged velocity with the unperforated solution; if the difference does not vanish as r→0 at the predicted rate, the cluster bound is violated. Alternatively, check the imported rigid-velocity bound directly for a pair of nearby balls: if the maximum speed grows faster than r^{-1/p}, the rate estimates break.","supporting_citations":[],"review_version":1}