Pith. sign in

REVIEW 2 major objections 7 minor 2 references

Floquet stability analysis of pulsatile particle-laden channel flow

T0 review · 2 major / 7 minor · reviewed 2026-08-08 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Pulsation destabilizes particle-laden channel flow at low frequency but stabilizes it at high frequency.

desk verdict Solid Floquet phase-diagram paper with one under-tested assumption (dropped n' perturbations) that should be verified before publication. read the letter →

arxiv 2608.04161 v1 pith:FJOGSAOV submitted 2026-08-04 physics.flu-dyn

Ananthapadmanabhan Ramesh , Benoit Pier , Parisa Mirbod This is my paper · ORCID
classification physics.flu-dyn
keywords hydrodynamicstabilityparticle-fluidinteractionspulsatileflowtime-periodicflowsFloquetanalysisdusty-gasmodelWomersleynumbersuspensiondynamics
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

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.

What carries the argument

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)$.

What would settle it

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.

Watch

Extended reading notes

Core claim

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.

Load-bearing premise

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.

Editorial extensions

If this is right

  • 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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 7 minor

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.

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 (2)
  1. [§3, after Eq. (3.15); Eq. (3.12); Appendix A] 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.
  2. [Appendix A] 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.
minor comments (7)
  1. [Eq. (2.23)] 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.
  2. [Table 1] 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.
  3. [Abstract and §5] 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.
  4. [§4.2.3, p. 21] The sentence 'increasing δ, t reduces the growth rate' contains a stray 't' and should be corrected.
  5. [Figure 4 caption] 'Womersely' is a typo for 'Womersley'.
  6. [§4, first paragraph] 'systemically varying' should read 'systematically varying'.
  7. [General] 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.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: all stability results are computed directly from the stated linearized two-phase equations, and the literature self-citations are contextual rather than load-bearing.

full rationale

The paper's derivation chain is self-contained. The base-state profiles in Eqs. (2.17)-(2.22) and the parameter A in Eq. (2.20) are obtained by solving the same two-phase momentum equations, not fitted to the stability results. The linearized perturbation system in Eqs. (3.12)-(3.15) is then reduced by the explicit assumption U approximately equal to U_p introduced after Eq. (3.15), which drops the n' coupling terms; this is a transparent modeling approximation and not a circular reuse of the target conclusion. The Floquet eigenvalue problem in Eqs. (3.17)-(3.19) is solved numerically with Chebyshev collocation, and the claimed destabilization-to-stabilization transition is read off from the computed neutral curves and Delta Re_cr phase diagrams in Figs. 6-11, so the central claim is a consequence of the calculation rather than an input. Validation is against two independent benchmarks external to the author group: Klinkenberg et al. (2011) for steady dusty-gas channel flow and Tsigklifis and Lucey (2017) for single-phase pulsatile flow, as shown in Fig. 4. The only self-citations are Pier and Schmid (2017) and Lebbal et al. (2022), both used merely as literature context in the introduction and not as justification for any equation, uniqueness claim, or fitting procedure. The reviewer concern that the n' decoupling is quantitatively unverified is a completeness or correctness issue, not a circularity, because no fitted parameter or prior conclusion is being recycled to produce that step.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The paper relies on the standard dusty-gas model and standard linear stability methods. The load-bearing simplifications are the neglect of particle number density perturbations and the two-dimensional restriction; both are acknowledged but not rigorously quantified.

assumptions (4)
  • domain assumption The dusty-gas model with Stokes drag only is appropriate for heavy particles in a dilute suspension.
    Section 2 states the restriction to heavy particles, for which only Stokes drag is considered, and excludes added mass, Basset history force, and particle-particle interactions for dilute suspensions.
  • ad hoc to paper The base-state particle number density is uniform, n=1, and perturbations of number density do not feed back into the velocity field.
    Section 2.1 sets n=1 in the base state; Section 3 eliminates the n' coupling using the U approximately equal to U_p assumption. This is a simplification specific to this paper's analysis.
  • domain assumption Two-dimensional disturbances are most unstable, so the analysis can be restricted to 2D modes.
    Appendix A presents a Squire-type transformation, but the authors acknowledge that a rigorous extension to time-periodic Floquet systems is non-trivial. The transformation preserves the equation structure but does not strictly prove that two-dimensional modes dominate for the Floquet problem.
  • domain assumption The base flow is taken to be spatially parallel and driven by a single-frequency pressure gradient.
    Section 2 restricts to single-frequency forcing to isolate the fundamental interaction, and assumes a uniform pressure gradient in x, so the base flow depends only on y and t.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Floquet stability analysis of pulsatile particle-laden channel flow." pith.science (2026). https://pith.science/paper/FJOGSAOV

@misc{pith2026260804161,
  author       = {Pith},
  title        = {Pith review of: Floquet stability analysis of pulsatile particle-laden channel flow},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FJOGSAOV}},
  note         = {Machine review of arXiv:2608.04161}
}
read the original abstract

We investigate the linear stability of particle-laden pulsatile channel flow using Floquet analysis within a two-phase dusty-gas framework. Uniformly distributed spherical particles are coupled to an incompressible Newtonian fluid through Stokes drag, and the governing equations are linearized about a time-periodic base flow driven by a sinusoidally varying pressure gradient. The effects of Reynolds number, Womersley number, pulsation amplitude, particle relaxation time, and particle mass fraction on temporal instability are examined. In the steady limit, particles with very short relaxation times destabilize the flow, whereas finite relaxation times introduce interphase slip and drag-mediated damping that stabilize disturbances. Under pulsatile forcing, increasing pulsation amplitude destabilizes the flow at low Womersley numbers but stabilizes it at sufficiently high Womersley numbers. This transition is governed by the penetration depth of oscillatory motion and is systematically shifted by particle relaxation time and mass loading through interphase momentum exchange. A critical corresponding value remains small throughout, indicating strong particle-fluid coupling and ruling out resonance-like particle dynamics. These findings provide a unified physical framework for the stability of pulsatile particle-laden flows relevant to physiological and periodically forced multiphase systems.

Figures

Figures reproduced from arXiv: 2608.04161 by the authors.

Figure 1
Figure 1. Schematic of pulsatile particle laden channel flow considered in the present [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Pulsating base state velocity profiles of fluid and particles presented in terms of [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Pulsating base state velocity profiles of fluid and particles shown in terms of the [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figures from the paper (8 more)
Figure 4
Figure 4. Figure 4: Validation of the model developed in the form of neutral stability curves. (a) [PITH_FULL_IMAGE:figures/full_fig_p013_4.png]
Figure 5
Figure 5. Figure 5: (a) Neutral stability curves for steady particle-laden Poiseuille flow at (f=0.1) for [PITH_FULL_IMAGE:figures/full_fig_p015_5.png]
Figure 6
Figure 6. Figure 6: Neutral stability curves for particle-laden pulsatile Poiseuille flow for different [PITH_FULL_IMAGE:figures/full_fig_p017_6.png]
Figure 7
Figure 7. Figure 7: Temporal evolution of the wall normal perturbation velocities of the fluid, ( [PITH_FULL_IMAGE:figures/full_fig_p018_7.png]
Figure 8
Figure 8. Figure 8: Real part of the Floquet growth rate, (𝜎𝑟 ), as a function of pulsation amplitude (𝛿) for different particle mass fractions ( 𝑓 ) and Womersley numbers (Wo), evaluated at fixed (𝑅𝑒=5772.22) and (𝑘=1.0206). All quantities are dimensionless. (a) 𝑆 = 1 × 10−7 , Wo=6,8,12;…
Figure 9
Figure 9. Figure 9: Variation of the critical Reynolds number, [PITH_FULL_IMAGE:figures/full_fig_p023_9.png]
Figure 10
Figure 10. Figure 10: Variation of critical Reynolds number, 𝑅𝑒𝑐𝑟 with pulsation amplitude 𝛿 for different particle mass fraction 𝑓 . (a) 𝑆 = 1 × 10−7 ; (b) 𝑆 = 2.5 × 10−4 . (a) (b) [PITH_FULL_IMAGE:figures/full_fig_p024_10.png]
Figure 11
Figure 11. Figure 11: Phase diagrams showing the variation of Δ𝑅𝑒𝑐𝑟 with Wo (6 -25) and 𝑆 (10−7 to 2.5 × 10−4 ) for 𝛿 = 0.1 and different 𝑓 : (a) 𝑓 = 0.05, and (b) 𝑓 = 0.1. tively. Figures 11a and 11b present phase diagrams of Δ𝑅𝑒𝑐𝑟 in the (Wo,𝑆) parameter space for 𝑓 = 0.05 and 𝑓 = 0.1. A…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

  1. [1]

    Asgharzadeh, Hafez & Borazjani, Iman2016 Effects of reynolds and womersley numbers on the hemodynamics of intracranial aneurysms.Computational and mathematical methods in medicine 2016(1), 7412926

    work page 2016

  2. [2]

    Straatman, AG, Khayat, RE, Haj-Qasem, E & Steinman, DA2002 On the hydrodynamic stability of pulsatile flow in a plane channel.Physics of Fluids14(6), 1938–1944

    work page 1938

Pith tools

Reviewed August 8, 2026 · model on record in the stance chip above.