{"id":"e3adf4a4-256a-445b-bf54-4d0c50eb956b","arxiv_id":"2608.04004","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Localizing polymers near the velocity maximum of a rectilinear viscoelastic flow strongly promotes the centre-mode elastic instability, reducing the critical Weissenberg number.","lead":"A linear stability analysis shows that concentrating polymers near the fastest part of a shear flow makes it go unstable at much lower elastic stress levels, while putting them near high shear suppresses the instability. This suggests microfluidic devices could trigger mixing by injecting polymer into just the central stream rather than premixing it everywhere.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Frozen-concentration assumption is load-bearing: the destabilization is attributed to P_∇φ (Eqs. 5.8-5.9), which vanishes if the gradient diffuses; the paper never quantifies the diffusion-vs-growth timescale at the predicted critical Wi.","rationale":"The paper's central quantitative result is the reduction of critical Weissenberg number by localizing polymers near the velocity maximum, demonstrated in both Kolmogorov flow (§3-4) and channel flow (§6). For this result to be physically meaningful, the nonuniform concentration field must persist long enough for the instability to develop. The model Eq. (2.5) has no diffusion, making the base state exact, but real polymer solutions have small but nonzero diffusivity. The energy budget (§5) identifies the concentration-gradient stress term P_∇φ as the sole positive work source at moderate loading, and this term is zero if the gradient erodes. Thus the entire mechanism is tied to the frozen-gradient assumption. The paper cites D ~ 10^-12 m^2/s and asserts diffusion is slow, but it never computes the ratio of the diffusion timescale to the instability growth timescale, which is the condition for the quasi-static approximation to be valid. At the neutral point the growth rate is zero, so the approximation is formally singular; one must consider a small supercritical Wi and compare the e-folding time to T_diff. This is the weakest link in the argument. The reader's verdict correctly identifies this as a conditional issue. I agree with the reader's conditional acceptance; the paper's benchmarks against Lewy & Kerswell (2025) and Khalid et al. (2021a,b) provide confidence in the numerics, and the energy analysis is internally consistent, but the frozen-profile assumption needs a quantitative sensitivity check before the strong claim of two-to-three order-of-magnitude reduction is fully established.","tokens_in":29231,"tokens_out":14388,"duration_ms":155201,"concrete_test":"Recompute the neutral stability curves of Fig. 10(b) (and Fig. 4(b)) with a small but finite polymer diffusivity, either by solving the base state with D∇^2φ in (2.5) at successive times and performing a quasi-static stability analysis as the interface broadens, or by adding a physical Pe^{-1}∇^2φ̂ term with Pe = V L / D ≈ 10^8 and measuring the shift in Wi_c. For the channel case (L ≈ 100 μm, V ≈ 1 mm/s, D ≈ 10^-12 m^2/s), the interface width δ ≈ 0.1L gives T_diff ≈ 10^3 s; compare this with the e-folding time of the mode at Wi = 1.1 Wi_c. If Wi_c rises by more than ~20% when δ doubles, or if the growth rate at the frozen-profile critical Wi is smaller than 1/T_diff, the frozen-profile claim overstates the promotion.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that localizing polymers near the velocity maximum lowers Wi_c by two to three orders of magnitude—relies on the base state φ̄(y) being steady. Equation (2.5) sets ∂_t φ + u·∇φ = 0, i.e., zero polymer diffusivity, so φ̄ is an exact steady state. The energy analysis (§5) shows that for moderate loading (β≲0.8) the sole positive work term is P_∇φ (Eqs. 5.8-5.9), which is proportional to φ̄' and to the concentration perturbation φ̂; from (2.17), φ̂ itself is forced by φ̄' v_y. Thus the predicted destabilization is entirely contingent on the persistence of the base concentration gradient. If physical polymer diffusion (D ~ 10^-12 m^2/s, cited in §2) erodes this gradient on a timescale T_diff ~ δ^2/D, where δ is the interface width, then the quasi-static base is not realized. The paper's only justification is the assertion that diffusion is 'very slow' (§2), but no comparison is made between T_diff and the e-folding time of the critical mode at the predicted Wi_c (or at a small supercritical Wi). At the linear threshold the growth rate is zero, so the quasi-static approximation is formally singular there; the relevant question is whether growth at finite supercritical Wi is fast enough. The artificial diffusion Pe^{-1}∇^2φ̂ added to the perturbation equations (§2.3, Fig. 2c) does not address base-state erosion: it only regularizes the continuous spectrum of (2.17). Consequently, the predicted two-to-three-order Wi_c reduction in channel flow (§6) may be an artifact of the infinite-T_diff limit. A quantitative test is needed to determine whether the prediction survives a physically reasonable finite diffusivity.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the linear stability of inertialess Oldroyd-B flow with a spatially nonuniform polymer concentration, focusing on the centre-mode elastic instability. In Kolmogorov flow, the authors solve the coupled eigenvalue problem using Fourier spectral collocation and show that localizing a fixed amount of polymer near the velocity maximum reduces the critical Weissenberg number, while localizing near the shear maximum increases it. As polymer loading increases, the neutral curve acquires a double-lobed shape, producing an abrupt drop in Wi_c. A kinetic-energy budget identifies a new work term P_∇φ, arising from base-state concentration gradients, as the dominant destabilizing contribution at moderate loading. The results are extended to channel flow, where centreline localization is shown to promote instability at Wi of order ten.","tokens_in":29658,"tokens_out":8013,"duration_ms":77919,"significance":"The paper is a well-executed numerical study of a novel and practically relevant question: how the spatial distribution of polymers affects the purely-elastic centre-mode instability. The main result—that localization near the velocity maximum destabilizes while localization near the shear maximum stabilizes—is clearly demonstrated and quantified. The identification of P_∇φ as the operative mechanism via an energy budget is insightful. The numerical methodology is careful, with convergence checks (N≈250), benchmarks against known results, and cross-validation of the neutral curves with and without the added perturbation diffusion. If the frozen-concentration assumption holds on the instability timescale, the predictions are experimentally testable and could inform microfluidic mixing strategies. The main caveat is the absence of a quantitative justification of the timescale separation underlying the frozen-concentration assumption.","major_comments":[{"comment":"The frozen-concentration assumption is load-bearing for the central claim. The destabilization at moderate polymer loading is attributed entirely to P_∇φ, which is proportional to the base gradient φ̄′ and to the concentration perturbation φ̂; from Eq. (2.17), φ̂ is forced by φ̄′ v_y. If physical polymer diffusion (D ≈ 10⁻¹² m² s⁻¹, as cited in §2) smears the base profile on a timescale T_diff ≈ δ²/D, and if this timescale is not long compared with the growth time of the critical mode at the predicted Wi_c (or at small supercritical Wi), the quasi-static base state is not realized and the predicted enhancement may be an artifact. The paper's assertion that diffusion is \"very slow\" (§2) is not quantified. At the linear threshold the growth rate is zero, so the quasi-static approximation is formally singular; the relevant check is the e-folding time at finite supercritical Wi. The added numerical diffusion Pe⁻¹∇²φ̂ in §2.3 only regularizes the continuous spectrum and does not address base-state erosion. Please provide an explicit timescale comparison using the physical parameters that set Wi, or include a calculation with a slowly diffusing base state to demonstrate that the qualitative result survives.","section":"Section 2 (Eq. 2.5) and Section 5 (Eqs. 5.8–5.9)"}],"minor_comments":[{"comment":"The phrase \"two to three orders of magnitude\" compares the nonuniform case at β=0.7 with the uniform case at β=0.994; the uniform case at β=0.7 is linearly stable, so please specify the baseline clearly to avoid overstatement.","section":"Section 6, Fig. 10(b)"},{"comment":"The pulse-profile notation is unclear: the parameters δ and ε appear in a garbled form in the equation, and the text should state explicitly which parameter controls the width and which controls the steepness, and which values are fixed.","section":"Eq. (2.9)"},{"comment":"The caption text for panel (c) erroneously refers to it as \"(d)\".","section":"Fig. 12 caption"},{"comment":"The statement that Squire's theorem remains valid for transverse concentration variations is not demonstrated; please include the short derivation or refer to supplementary material for the proof.","section":"Section 2.3"},{"comment":"The definition of the inner product ⟨·⟩ is given after the equation; move it before or immediately after the equation for readability.","section":"Section 5, Eq. (5.1)"},{"comment":"The notation \"±1.0265 + 0.03393 i\" is ambiguous; please write \"1.0265 ± 0.03393 i\" for the pair of eigenvalues.","section":"Fig. 2"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope and the numerical work appears thorough. The main concern is the frozen-concentration assumption, which is load-bearing for the central claim. If the authors provide a convincing timescale comparison or an additional calculation with a slowly diffusing base state, the paper would be suitable for publication. No concerns about citation or novelty."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this paper does a clean job on its own terms. It extends the centre-mode linear stability analysis of Kolmogorov and channel flows to a base state with a frozen, transversely varying polymer concentration, and reports large destabilization when the polymers sit around the velocity maximum. The numerics look solid: spectral collocation, convergence around N≈250, benchmarks against Lewy & Kerswell and Khalid et al. The energy budget is informative: at moderate loading, the dominant positive work term is P_∇φ, which arises from the base concentration gradient, and the eigenfunctions are indeed focused at the interface. The double-lobed neutral curves are a nice, non-obvious signature. Nothing here is circular: the destabilization is an eigenvalue trend, not a fit.\n\nThe soft spot is exactly where the stress-test put it. The base concentration φ̄ is frozen by dropping diffusion (Eq. 2.5). All the destabilizing physics at moderate β comes from P_∇φ, which is present only because φ̄' is nonzero. If the polymer diffuses on a timescale comparable to the instability growth timescale, the gradient erodes and the effect weakens. The paper says polymer diffusion is 'very slow' and then proceeds quasi-statically, but it never computes the interface erosion time against the e-folding time of the critical mode at a small supercritical Wi. At threshold the growth rate is zero, so that comparison is exactly what is needed; without it, the two-to-three order-of-magnitude Wi_c reduction in channel flow remains a prediction of the frozen-concentration model rather than something I can treat as a physical quantity.\n\nThere are two smaller points. The extension of Squire's theorem to nonuniform concentration is asserted in one sentence; the derivation should be shown or at least sketched, though I have no reason to think it fails. And there is no code or data archive, which makes independent spot-checks harder than they should be for a computational study.\n\nWho should read this: anyone working on elastic instabilities in rectilinear flows or on microfluidic mixing with spatially structured polymer concentration. It deserves a serious referee and, if the timescale question is addressed with a finite-diffusivity check or an order-of-magnitude estimate, it could be a good publication. Send it to review, with the frozen-diffusion timescale as the major comment.","headline":"A credible linear-stability study claiming polymer localization near the velocity maximum lowers the critical Wi by orders of magnitude, but its frozen-concentration base state may undermine the quantitative promise.","tokens_in":30142,"tokens_out":2507,"would_cite":true,"duration_ms":28046,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Localizing a fixed amount of polymer near the base-flow velocity maximum strongly promotes the inertialess centre-mode elastic instability, lowering the critical Weissenberg number by two to three orders of magnitude in channel flow.","keywords":["centre-mode instability","elastic instability","Oldroyd-B fluid","nonuniform polymer concentration","Weissenberg number","Kolmogorov flow","viscoelastic channel flow","concentration-gradient feedback"],"falsifier":"Run the same base state with a finite polymer diffusivity (nonzero inverse Péclet number) in the linearized equations, or in a microfluidic experiment with a polymer stream of known diffusivity, and measure the critical Weissenberg number and the unstable eigenfunctions as the Péclet number varies; the claimed mechanism predicts a rise in the threshold and a loss of interface-focused eigenfunctions once diffusion erodes the concentration gradients on the instability timescale.","tokens_in":29036,"feed_emoji":"🧪","tokens_out":8789,"duration_ms":91103,"temperature":0.7,"pith_summary":"The paper asks whether the centre-mode elastic instability—the purely elastic, inertia-free instability of rectilinear flows of dilute polymer solutions—can be controlled by where the polymers sit rather than by how much polymer is present. Using an Oldroyd-B fluid in the Stokes limit with a frozen transverse concentration profile, it compares a uniformly premixed solution with a polymer-laden stream sandwiched between solvent streams while holding total polymer mass fixed. It finds that localizing polymers near the maximum of the base velocity strongly promotes the instability, lowering the critical Weissenberg number, while localizing them near the maximum shear suppresses it. In channel flow the reduction is by two to three orders of magnitude. An energy budget traces the enhancement to elastic feedback forces produced by concentration gradients, not to the modified base velocity or polymer stretching.","feed_headline":"Polymers at the flow centre cut instability onset by up to 1000x","feed_subtitle":"Channel flow becomes unstable at Weissenberg numbers near 10 instead of 1000 when the polymer stream hugs the velocity maximum.","key_machinery":"The machinery is the Oldroyd-B linear stability problem with a concentration-dependent stress, $T = \\phi (1-\\beta)/Wi\\,(C-I)$, and a scalar advection equation for the normalized polymer concentration $\\phi$. The base state is a unidirectional shear flow with a frozen, prescribed pulse profile $\\bar\\phi(y)$, with total polymer loading fixed by the normalization $\\int \\bar\\phi\\,dy = 1$. The argument is carried by separating the perturbation polymer stress into a part proportional to the base concentration and a part proportional to the concentration perturbation; this separation generates a new term $P_{\\nabla\\phi}$ in the kinetic-energy budget that is nonzero only when $\\bar\\phi'$ is nonzero. Comparison of neutral curves in the $Wi$–$k$ plane, with the continuous spectrum stabilized by weak artificial diffusion, and the eigenfunction energy budget together identify the concentration-gradient feedback as the mechanism that lowers the critical Weissenberg number.","core_discovery":"The central claim is that a spatially nonuniform polymer distribution is not a minor perturbation of the uniform case: placing a fixed amount of polymer in a stream centred on the velocity maximum destabilizes the flow in a way that a uniform solution of the same total loading cannot match. In periodic Kolmogorov flow the critical Weissenberg number $Wi_c$ decreases as the concentration profile sharpens from uniform to a near-square pulse, and for moderate polymer loading the neutral curve develops a second lobe, producing a sudden drop in $Wi_c$. For channel flow, localizing polymers about the centreline makes the purely elastic centre-mode unstable at $Wi = O(10)$, whereas the uniform Oldroyd-B channel requires $Wi \\approx 10^3$ and is unstable only in the very dilute regime. The paper argues that this is not simply a higher effective concentration: the perturbation eigenfunctions become confined to the interfaces between polymer-rich and polymer-free streams, and the kinetic-energy budget shows that the sole positive work source at moderate loading is $P_{\\nabla\\phi}$, the polymer power arising from base-state concentration gradients.","pith_inferences":["A testable microfluidic strategy follows: inject polymer solution only into the central inlet of a three-stream channel rather than premixing it everywhere, and measure the onset of elastic instability at lower flow rates; the paper's channel results predict a threshold reduction of two to three orders of magnitude if the interface can be kept sharp.","The concentration-gradient feedback is likely to modify other viscoelastic instabilities whenever the base concentration gradient overlaps the perturbation stress field, such as the elasto-inertial centre-mode, polymer-diffusive wall modes, and hoop-stress instabilities in curved flows, so polymer localization could serve as a general suppression or promotion route.","Because the frozen-profile idealization ignores diffusion, the physically relevant control parameter should be the ratio of interface diffusion time to instability growth time; a finite-diffusivity model should show the enhancement degrading once the polymer Péclet number falls below a threshold.","The double-lobed neutral curves suggest an avoided-crossing-like interaction between a centre mode and an interface mode, which could be probed by tracking the two eigenvalue branches as the polymer loading is varied."],"forward_implications":["In channel flow of a sufficiently loaded dilute solution, feeding a central polymer stream can make the inertialess centre-mode unstable at $Wi \\sim 10$, roughly three orders of magnitude below the uniform premixed case.","For a fixed total amount of polymer, the location of the polymer stream becomes a design choice: place it near the velocity maximum to promote instability, or near the maximum-shear region to suppress it.","At moderate polymer loading the neutral curve is double-lobed, so the critical mode switches from a fast centre-travelling mode to a slower interface-focused mode, with a sudden drop in the critical Weissenberg number.","The polymer-work term from concentration gradients becomes the sole positive energy source for $\\beta \\lesssim 0.8$ in Kolmogorov flow, marking where uniform-polymer intuition fails.","Stream width can be tuned: narrower polymer streams are most destabilizing at low loading, while wider streams are most destabilizing at higher loading."],"supporting_citations":[{"why":"Shows that uniform Oldroyd-B channel flow is inertially unstable via the centre-mode only for very dilute solutions with β ≳ 0.99, the baseline that polymer localization overturns.","marker":"Khalid et al. (2021b)"},{"why":"Establishes the requirement of a base-flow velocity maximum for the centre-mode and supplies the asymptotic description used to motivate the Kolmogorov and channel models.","marker":"Kerswell & Page (2024)"},{"why":"Confirms the centre-mode as the relevant instability of two-dimensional viscoelastic Kolmogorov flow and provides uniform-concentration spectra and neutral curves used for benchmarking.","marker":"Lewy & Kerswell (2025)"},{"why":"Supplies the governing equations and the scenario of polymers localized in a layer within a shear flow that the paper adapts to the inertialess regime.","marker":"Vaithianathan et al. (2006, 2007)"},{"why":"Documents the very small diffusivity of long polymer molecules, supporting the frozen-concentration quasi-static base state on which the analysis rests.","marker":"El-Kareh & Leal (1989)"},{"why":"Provides the kinetic-energy budget formulation whose polymer-work terms are split to isolate the concentration-gradient contribution.","marker":"Zhang et al. (2013)"},{"why":"Proves Squire's theorem for Oldroyd-B with uniform concentration, a result the paper extends to transverse concentration variations to justify the two-dimensional analysis.","marker":"Bistagnino et al. (2007)"},{"why":"Shows that weak polymer diffusion only weakly stabilizes the centre-mode in Kolmogorov flow, supporting the artificial-diffusion technique used to track neutral curves.","marker":"Nichols et al. (2025)"}],"fun_headline_variants":["Polymer localization cuts elastic instability onset by 100x","Center-focused polymers drop Wi onset to ~10 in channels","Localizing polymers to flow core ignites elastic instability","Polymer stream at centerline destabilizes flow at Wi ~10"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the base-state polymer concentration is a frozen, non-diffusing field, so the sharp interfaces between polymer-rich and polymer-free streams persist while the instability grows; if polymer diffusion smears those interfaces on the instability timescale, the concentration-gradient feedback identified as the driver would be weakened or lost.","fun_headline_variants_meta":{"raw":{"variants":["Polymer localization cuts elastic instability onset by 100x","Center-focused polymers drop Wi onset to ~10 in channels","Localizing polymers to flow core ignites elastic instability","Polymer stream at centerline destabilizes flow at Wi ~10"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000732,"raw_usage":{"total_tokens":3307,"prompt_tokens":1011,"completion_tokens":2296,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":627,"completion_tokens_details":{"reasoning_tokens":2228}},"tokens_in":627,"tokens_out":2296,"duration_ms":20190,"temperature":1.0,"reasoning_tokens":2228,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T00:34:15.421364+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same base state with a finite polymer diffusivity (nonzero inverse Péclet number) in the linearized equations, or in a microfluidic experiment with a polymer stream of known diffusivity, and measure the critical Weissenberg number and the unstable eigenfunctions as the Péclet number varies; the claimed mechanism predicts a rise in the threshold and a loss of interface-focused eigenfunctions once diffusion erodes the concentration gradients on the instability timescale.","supporting_citations":[],"review_version":2}