{"id":"edfc696d-91d8-42ed-944c-69e09a85e328","arxiv_id":"2608.04161","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"In pulsatile particle-laden channel flow, low-frequency forcing destabilizes the flow while high-frequency forcing stabilizes it, with particle inertia and loading shifting the transition boundary.","lead":"Using computer simulations of the equations describing a fluid carrying tiny suspended particles, the authors map when a pulsating flow through a channel becomes unstable. They show that slow pulses destabilize the flow while fast pulses confine their effect near the walls and stabilize it, with particle inertia shifting the boundary between regimes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The omission of particle number-density perturbations from the momentum equations (after Eq. 3.15) is quantitatively unverified at the largest relaxation times and could shift the reported destabilization-stabilization boundary.","rationale":"The paper is generally careful: it validates against Klinkenberg et al. (2011) and Tsigklifis & Lucey (2017), reports spectral convergence, and its qualitative trend (low-Wo destabilization and high-Wo stabilization) is physically plausible. The Squire extension is explicitly marked as non-rigorous, but a Squire-type transformation for a time-periodic base flow with no spanwise dependence is likely valid, so I do not treat that as the main risk. The n' decoupling is the true soft spot: it is introduced as a simplifying assumption, the paper only says it holds for sufficiently small relaxation time, and the quantitative measure S Wo^2 used to justify strong coupling reaches about 0.156 at the largest S and Wo, where the base slip is roughly 15% of the oscillatory amplitude. Because the same prefactor f/(S Re) appears in both the retained drag terms and the dropped n' term, there is no demonstrated hierarchy showing the dropped term is negligible. Restoring n' in the eigenvalue problem is straightforward and would settle whether the phase diagram in Fig. 11 is robust; this concern does not overturn the qualitative picture but justifies the conditional verdict.","tokens_in":22852,"tokens_out":8428,"duration_ms":84407,"concrete_test":"Restore n' as an active unknown in the Chebyshev-Floquet eigenvalue problem, i.e., solve the full coupled system Eqs. (3.12)-(3.15) without setting U=U_p, and recompute the Delta Re_cr phase diagram for S=2.5e-4, f=0.1, delta=0.1 over Wo=6-25. Compare the Womersley number at which Delta Re_cr changes sign with Fig. 11b. If the zero crossing moves by more than 1-2 Womersley units, or if Re_cr changes by more than 5% at Wo=20, the quantitative boundary in the central claim is not robust to the discarded n' coupling; if the results match, the assumption is validated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing step is the decoupling of n' from the momentum equations, introduced just after Eq. (3.15) and restated in Appendix A. The authors assume U approximately equals U_p so that the term -(f/(S Re))(U-U_p)n' in Eqs. (3.4)/(3.12) vanishes. However, the Floquet system (3.17)-(3.19) retains distinct base-state harmonics q1 and q2 for the fluid and particle phases, and q2 differs from q1 by the factor (1+i S Wo^2)^{-1} (Eqs. 2.19 and 2.22). For S=2.5e-4 and Wo=20-25, S Wo^2 is about 0.1-0.156, so the base slip is not negligible. The dropped term has the same prefactor f/(S Re) as the retained drag terms, so there is no a priori hierarchy showing it is small. Since the paper's central mechanism is interphase momentum exchange, an unquantified contribution from n' coupling could shift the transition boundary in Fig. 11, which the conclusions are built around. The particle continuity equation (3.13) is simply dropped from the eigenvalue problem without a convergence test that restores it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the linear stability of a particle-laden channel flow driven by a sinusoidally varying pressure gradient. A two-phase dusty-gas model with Stokes drag is linearized about a time-periodic base flow, and Floquet theory is used to compute growth rates from a truncated harmonic expansion. The paper reports neutral curves and critical Reynolds numbers over ranges of Reynolds number, Womersley number, pulsation amplitude, particle relaxation time, and mass fraction. The central claim is that pulsation destabilizes the flow at low Womersley numbers and stabilizes it at high Womersley numbers, with particle relaxation time and mass loading shifting the transition through interphase momentum exchange; the small values of S Wo^2 along the transition are used to argue against a resonance-like mechanism.","tokens_in":23119,"tokens_out":8401,"duration_ms":83584,"significance":"The problem is timely, and mapping the destabilization-stabilization transition in (Wo, S) space is a potentially useful contribution to the pulsatile multiphase literature. The paper has clear strengths: no parameters are fitted to the results; the numerical implementation is validated against two independent benchmarks (steady particle-laden channel flow and single-phase pulsatile flow, Figure 4); and Table 1 provides a convergence study. If the two approximations discussed in the major comments are properly quantified, the proposed physical picture of oscillatory penetration combined with drag-mediated momentum exchange is compelling and falsifiable. At present, however, the central quantitative claims rest on assumptions that are not verified within the manuscript.","major_comments":[{"comment":"The decoupling of the particle number-density perturbation n' is the most load-bearing step of the analysis. The text justifies the decoupling by U≈U_p, but the retained base flow already contains distinct harmonics q1 and q2 in Eqs. (2.19) and (2.22), with q2 = q1/(1 + i S Wo^2). At S = 2.5e-4 and Wo = 20-25 this gives S Wo^2 = 0.1-0.156, so the base-state slip (U - U_p) is not negligible. The dropped terms in Eq. (3.12) are -(f/(S Re)) n' (∂U/∂y - ∂U_p/∂y) and -(f/(S Re)) ∂n'/∂y (U - U_p); both carry the same prefactor f/(S Re) as the retained drag terms, so there is no a priori hierarchy showing they are small. Since the paper attributes the transition-boundary shifts in Fig. 11 to interphase momentum exchange, the authors should either retain Eq. (3.13) together with the n' couplings in the eigenvalue problem and show convergence of ΔRe_cr, or provide a quantitative error estimate that demonstrates the neglected terms are negligible in the parameter ranges used.","section":"§3, after Eq. (3.15); Eq. (3.12); Appendix A"},{"comment":"The Squire-type reduction to two-dimensional disturbances is not established for this system. The appendix first assumes U = U_p to remove the n' coupling, which is inconsistent with the distinct q1 and q2 harmonics retained throughout the paper. The transformed particle momentum equations then use 1/(S Re) in Eq. (A13) but 1/(S \\tilde Re) in Eq. (A14), so the transformed system is not identical to the two-dimensional problem at the transformed Reynolds number. Consequently, the statement that \\tilde Re ≤ Re implies the least stable modes are two-dimensional does not follow from the given transformation. Because all neutral curves and phase diagrams are computed in two dimensions, the authors should either give a correct transformation with appropriately rescaled particle parameters or test representative three-dimensional modes numerically to confirm that the two-dimensional restriction does not change the reported transition boundaries.","section":"Appendix A"}],"minor_comments":[{"comment":"The expression for δ/Λ contains e^{±it} even though δ is defined as the constant amplitude of the flow-rate oscillation; the correct amplitude relation should be stated explicitly.","section":"Eq. (2.23)"},{"comment":"The imaginary part σ_i jumps between different branches (0.2338302, -0.766169, -0.7661589) as M increases; the authors should state that Floquet exponents are defined modulo integer shifts and report the branch convention used.","section":"Table 1"},{"comment":"The phrase 'a critical corresponding value remains small' is vague; it should specify S Wo^2, since that is the quantity used to rule out resonance-like particle dynamics.","section":"Abstract and §5"},{"comment":"The sentence 'increasing δ, t reduces the growth rate' contains a stray 't' and should be corrected.","section":"§4.2.3, p. 21"},{"comment":"'Womersely' is a typo for 'Womersley'.","section":"Figure 4 caption"},{"comment":"'systemically varying' should read 'systematically varying'.","section":"§4, first paragraph"},{"comment":"Several JFM template artifacts appear in the text (e.g., 'Rapids articles must not exceed this page length' and 'Focus on Fluids articles must not exceed this page length') and should be removed before resubmission.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for JFM and the numerical validation is a genuine strength. The decisive issue is the unquantified neglect of n' perturbations: if the authors can restore the coupling or provide a rigorous error bound, the central claims may be salvageable; if not, the quantitative transition-boundary claims in Fig. 11 and the resonance-exclusion interpretation would need to be substantially softened. The self-citation to Pier and Schmid (2017) is a relevant co-authored prior study and is not a concern."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague—\n\nBottom line: this is a serious, competently executed Floquet stability study that gives the first unified (Wo,S) phase diagram for pulsatile particle-laden channel flow. The central claim—low-Womersley pulsation destabilizes, high-Womersley stabilizes, and particle inertia/loading shifts the boundary via interphase momentum exchange rather than resonance—looks qualitatively right. The numerical work is solid: convergence table, validation against Klinkenberg et al. (2011) and Tsigklifis & Lucey (2017), no fitted parameters. I would send this to a referee.\n\nThe main soft spot is the one the stress test flags: the particle number-density perturbation n' is dropped from the momentum equations on the assumption U≈Up, made just after Eq. (3.15) and again in Appendix A. But the base state itself has q2 = q1/(1+i S Wo^2) (Eqs. 2.19, 2.22), so the slip is not zero. At S=2.5e-4 and Wo=20–25, |S Wo^2| is 0.1–0.156, which is not negligible, and the discarded term carries the same prefactor f/(S Re) as the retained drag terms. There is no estimate or convergence test that restores n' and shows the boundary in Fig. 11 is insensitive. This is a real gap, but it is addressable: include the n' equation in the eigenproblem, or at least do a few cases with it retained and show ΔRe_cr shifts by less than some tolerance. My guess is the qualitative transition survives, but the boundary location could move.\n\nSecond soft spot is more minor: the Squire extension to Floquet systems is explicitly non-rigorous. The transformation changes the effective Womersley number of the 2D problem (since Wo^2/Re transforms), so the usual “2D is worst” conclusion does not follow automatically. The authors admit this and proceed anyway. That is a defensible choice, but the phase diagram should be read as a 2D prediction until a 3D check is done.\n\nThe citation pattern is fine. The only self-citation is Pier & Schmid (2017), directly relevant and coauthored; no parameter fitting, no invented results. The paper is honest about its assumptions.\n\nWho it is for: people working on stability of time-periodic shear flows or particle-laden channel flows. It deserves a serious referee and would likely come back with a request for the n' sensitivity test and a sharper Squire discussion. I would not desk-reject it.","headline":"Solid Floquet phase-diagram paper with one under-tested assumption (dropped n' perturbations) that should be verified before publication.","tokens_in":23626,"tokens_out":3017,"would_cite":true,"duration_ms":31448,"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":"Pulsation destabilizes particle-laden channel flow at low frequency but stabilizes it at high frequency.","keywords":["hydrodynamic stability","particle-fluid interactions","pulsatile flow","time-periodic flows","Floquet analysis","dusty-gas model","Womersley number","suspension dynamics"],"falsifier":"Evaluate the dropped term $(f/(S\\,Re))(U-U_p)n'$ using the published base-state profiles at a case where $S Wo^2$ is not tiny (e.g., $S=2.5\\times10^{-4}$, $Wo=25$, where $S Wo^2=0.156$); if its magnitude is comparable to the retained drag terms, the decoupling assumption is not safe. A sharper experimental test is to measure the critical Reynolds number at $Wo \\approx 12$-$15$ for two mass fractions: the predicted shift of the crossover to higher $Wo$ with loading is a specific signature that would confirm or destroy the drag-mediated mechanism.","tokens_in":22662,"feed_emoji":"🌊","tokens_out":9888,"duration_ms":88351,"temperature":0.7,"pith_summary":"This paper asks how a periodically varying pressure gradient changes the first onset of instability in a channel carrying a dilute suspension of small spherical particles. The central claim is that pulsation reverses its effect as the Womersley number, the ratio of the forcing frequency to the rate at which viscous effects spread across the channel, increases: at low Womersley numbers the oscillation reaches across the channel and increasing pulsation amplitude lowers the critical Reynolds number, while at high Womersley numbers the oscillation is confined to thin wall layers and increasing amplitude raises it. Particles do not create this crossover; relaxation time and mass loading shift it through interphase drag, and the smallness of $S Wo^2$ along the boundary rules out a resonance-like particle response. A sympathetic reader would care because the same mechanism should appear wherever pulsing and suspended particles meet, from blood and aerosol transport to periodically forced multiphase devices.","feed_headline":"Pulsing flow: low frequency destabilizes, high frequency stabilizes","feed_subtitle":"The same penetration rule tells when periodic pumping promotes or suppresses instability in particle suspensions.","key_machinery":"The load-bearing object is the Floquet formulation of the linearized two-phase dusty-gas equations: perturbations are written as $e^{ikx}e^{\\sigma t}\\sum_m \\hat\\psi_m(y)e^{imt}$, with analogous expansions for the particle velocities, turning the time-periodic problem into a generalized eigenvalue problem $\\mathbf{A}X=\\sigma \\mathbf{B}X$, discretized with Chebyshev collocation in the wall-normal direction and a truncated Fourier series in time. The key derived quantity is the base-state wavenumber $A$, whose imaginary part controls the penetration depth of the oscillatory motion, and the coupling parameter $S Wo^2$, whose smallness indicates that particles stay strongly coupled to the fluid. The analysis also uses the assumption $U\\approx U_p$ to decouple the particle number-density perturbation and a Squire-type transformation to restrict the search to two-dimensional disturbances. This machinery allows the authors to map the stability boundary in the $(Wo,S)$ plane and to quantify the crossover through $\\Delta Re_{cr}=Re_{cr}(\\delta=0.1)-Re_{cr}(\\delta=0)$.","core_discovery":"The paper's central discovery is that pulsatile forcing has a frequency-dependent, non-monotonic effect on the linear stability of particle-laden channel flow, and that the reversal between destabilization and stabilization is controlled by the penetration of the unsteady shear layer together with particle-fluid drag. In the steady limit, very light tracer-like particles slightly destabilize the flow, while particles with finite relaxation time stabilize it through interphase slip and drag-mediated damping. With pulsation, low-Womersley-number forcing ($Wo \\approx 6$-$8$) modulates shear through the bulk and lowers the critical Reynolds number as the amplitude $\\delta$ increases; high-Womersley-number forcing ($Wo \\ge 15$-$20$) confines the oscillation to near-wall Stokes layers and raises the critical Reynolds number. Particle relaxation time $S$ and mass fraction $f$ do not change this mechanism but systematically shift the crossover, with larger mass loading moving the boundary toward higher $Wo$. Because $S Wo^2$ remains small along the transition, the authors conclude that the crossover emerges from the coupling of oscillatory penetration and interphase momentum exchange rather than from resonance-like particle dynamics.","pith_inferences":["Beyond the paper, the same penetration-depth logic should carry over to pulsatile pipe flow and other periodically forced shear layers, where the Womersley number plays the same role and the crossover frequency would scale with the viscous diffusion time across the geometry.","A direct experimental check is possible: in a controlled oscillatory channel, time-resolved velocity measurements of disturbance growth at $Wo \\approx 12$-$15$ with varying particle loading should show the crossover shifting to higher frequency as the mass fraction rises.","Since the analysis is linear and dilute, at larger pulsation amplitudes or volume fractions one might expect nonlinear disturbance interactions, concentration migration, and the neglected number-density coupling to modify the boundary; the present framework gives a baseline against which those effects could be measured."],"forward_implications":["Raising the pulsation amplitude at $Wo \\approx 6$-$8$ lowers the critical Reynolds number across all particle relaxation times studied, so pulsation acts as a destabilizer in that regime.","At $Wo \\ge 15$-$20$, raising the pulsation amplitude increases the critical Reynolds number for every particle relaxation time considered, so pulsation suppresses instability.","For tracer-like particles ($S=10^{-7}$), adding particle mass lowers the critical Reynolds number at all Womersley numbers; for finite-inertia particles ($S=2.5\\times10^{-4}$), adding mass raises it.","Increasing the mass fraction from $f=0.05$ to $f=0.1$ shifts the destabilization-to-stabilization boundary toward larger $Wo$, expanding the destabilizing region.","The smallness of $S Wo^2$ along the boundary rules out a resonance-like particle response and indicates strong fluid-particle coupling throughout the parameter space."],"supporting_citations":[{"why":"Supplies the small-Stokes-number limit in which fluid and particle base velocities are taken nearly equal, justifying the decoupling of the number-density perturbation.","marker":"Chen & Chung 1995"},{"why":"Provides the steady particle-laden channel stability results used as a validation baseline and as the reference for particle-induced destabilization and stabilization.","marker":"Klinkenberg et al. 2011"},{"why":"Provides the pulsatile single-phase channel result against which the particle-free limit of the present Floquet computation is validated.","marker":"Tsigklifis & Lucey 2017"},{"why":"Introduces the dusty-gas two-phase formulation used here and the tracer-like destabilization mechanism quoted for very short relaxation times.","marker":"Saffman 1962"},{"why":"Demonstrates in single-phase oscillatory channel flow that pulsation can stabilize or destabilize depending on forcing conditions, the effect the particles are shown to shift.","marker":"Von Kerczek 1982"},{"why":"Provides the modern linear and nonlinear treatment of pulsatile channel flow that the present work extends to a particle-laden suspension.","marker":"Pier & Schmid 2017"}],"fun_headline_variants":["Pulsatile forcing has dual effect on particle flow stability","Low-frequency pulsation destabilizes, high-frequency stabilizes suspension","Penetration depth governs pulsatile stability reversal in laden flow","Pulsing channel flow: particle drag shifts stability crossover","Frequency flips stability of particle-laden pulsing flows"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the particles almost keep up with the fluid in the base flow, so that particle number-density fluctuations can be ignored in the momentum equations; if particle-fluid slip is not small, the neglected coupling term could shift or erase the predicted transition.","fun_headline_variants_meta":{"raw":{"variants":["Pulsatile forcing has dual effect on particle flow stability","Low-frequency pulsation destabilizes, high-frequency stabilizes suspension","Penetration depth governs pulsatile stability reversal in laden flow","Pulsing channel flow: particle drag shifts stability crossover","Frequency flips stability of particle-laden pulsing flows"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000247,"raw_usage":{"total_tokens":1561,"prompt_tokens":981,"completion_tokens":580,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":597,"completion_tokens_details":{"reasoning_tokens":497}},"tokens_in":597,"tokens_out":580,"duration_ms":5505,"temperature":1.0,"reasoning_tokens":497,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T00:21:04.364753+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate the dropped term $(f/(S\\,Re))(U-U_p)n'$ using the published base-state profiles at a case where $S Wo^2$ is not tiny (e.g., $S=2.5\\times10^{-4}$, $Wo=25$, where $S Wo^2=0.156$); if its magnitude is comparable to the retained drag terms, the decoupling assumption is not safe. A sharper experimental test is to measure the critical Reynolds number at $Wo \\approx 12$-$15$ for two mass fractions: the predicted shift of the crossover to higher $Wo$ with loading is a specific signature that would confirm or destroy the drag-mediated mechanism.","supporting_citations":[],"review_version":1}