pith. sign in

arxiv: 2511.07661 · v4 · submitted 2025-11-10 · ⚛️ physics.flu-dyn · physics.bio-ph

Resonant spectral cascade in Womersley flow triggered by arterial geometry

Khalid M. Saqr This is my paper

Pith reviewed 2026-05-17 23:03 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn physics.bio-ph
keywords arterial geometryWomersley flowspectral broadeningpulsatile flowresonant energy transfervascular remodelingflow stability
0
0 comments X

The pith

Arterial geometry triggers resonant transfer of energy to short-wavelength flow components in pulsatile Womersley flow.

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

The paper solves a mathematical model to show that arterial geometry actively triggers a resonant transfer of energy to short-wavelength components of pulsatile flow. Across physiological Womersley numbers the overall wave energy decays with a negative growth rate, yet the spectral broadening ratio peaks sharply at an intermediate value. This identifies a resonant frequency where geometry generates maximum spectral complexity. The result reframes geometry from a passive dissipator to an active modulator of flow spectra, with possible use in vascular health diagnostics.

Core claim

Arterial geometry can trigger a resonant transfer of energy to short-wavelength components of the flow. The global wave energy consistently decays, confirmed by a negative growth rate G less than zero, indicating the flow does not become exponentially unstable. However, the spectral broadening ratio R, which quantifies energy in high-wavenumber versus low-wavenumber modes, exhibits a sharp non-monotonic peak at an intermediate Womersley number. This identifies a resonant frequency at which geometry is maximally efficient at generating spectral complexity even as the overall flow attenuates.

What carries the argument

The spectral broadening ratio R that measures the redistribution of energy from low- to high-wavenumber modes under the influence of arterial geometry in the Womersley flow model.

If this is right

  • Geometry modifies pulsatile phase relationships, near-wall shear, and axial transport through spectral energy redistribution.
  • Age-related changes such as elongation and tortuosity actively influence flow dynamics beyond simply increasing viscous resistance.
  • Spectral diagnostics could serve as sensitive markers for vascular health and geometric remodeling.
  • The flow remains globally stable yet develops greater spectral complexity precisely at the resonant frequency.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Resonant frequencies may align with clinical observations of altered flow in elongated or tortuous arteries.
  • The same geometry-driven mechanism could be tested in other pulsatile systems such as respiratory airways or industrial pipes.
  • Interventions that change vessel geometry might shift or suppress the resonant peak and thereby alter local shear patterns.

Load-bearing premise

The mathematical model accurately represents the effects of realistic arterial geometry on pulsatile flow without simplifications that would eliminate the reported resonance.

What would settle it

A numerical simulation or laboratory measurement of pulsatile flow through a realistic arterial geometry that shows no peak in the spectral broadening ratio at any intermediate Womersley number.

Figures

Figures reproduced from arXiv: 2511.07661 by Khalid M. Saqr.

Figure 1
Figure 1. Figure 1: Initial setup for the simulation with Wo = 10. Top: The multi-harmonic initial [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Primary instability diagnostics for Wo = 10. Top: Total wave energy [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Modal energy evolution for the first five harmonics ( [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Waveform shape evolution at selected time snapshots for Wo = 10. The plot [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Evolution of the classical KdV invariants for Wo = 10. Top: Momentum [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Spatiotemporal evolution of the wave and its spectrum for Wo = 10. Top: [PITH_FULL_IMAGE:figures/full_fig_p013_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Summary of instability metrics from the parameter sweep versus Womersley [PITH_FULL_IMAGE:figures/full_fig_p014_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Parametric comparison of final states from the sweep. Top: Overlay of the [PITH_FULL_IMAGE:figures/full_fig_p015_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Parametric structure of final states represented as heatmaps. The vertical axis [PITH_FULL_IMAGE:figures/full_fig_p016_9.png] view at source ↗
read the original abstract

Age-related arterial remodeling is dominated by progressive loss of elastic-fiber function and concomitant stiffening, and in many vascular beds it is also accompanied by measurable geometric remodeling (e.g., elongation and tortuosity). These changes are clinically relevant because they modify pulsatile phase relationships, near-wall shear, and axial transport, yet the precise physical mechanisms by which geometry modulates spectral energy redistribution remain insufficiently resolved. While complex geometry is known to increase viscous resistance, its active role in modulating flow dynamics is not fully understood. Here we solve a mathematical model to show that arterial geometry can trigger a resonant transfer of energy to short-wavelength components of the flow. The investigation, conducted over a physiological range of Womersley numbers (Wo, a dimensionless measure of pulsation frequency), reveal a dual dynamic. The global wave energy consistently decays, confirmed by a negative growth rate (G < 0), indicating that the flow does not become exponentially unstable. However, a spectral broadening ratio (R), which quantifies the energy in high-wavenumber versus low-wavenumber modes, exhibits a sharp, non-monotonic peak at an intermediate Wo. This result identifies a resonant frequency at which geometry is maximally efficient at generating spectral complexity, even as the overall flow attenuates. These findings reframe the role of arterial geometry from a passive dissipator to an active modulator of the flow's spectral content, suggesting that spectral diagnostics could provide a sensitive marker for vascular health.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript solves a mathematical model of pulsatile Womersley flow in a geometrically complex arterial domain. It reports that while the global growth rate G remains negative (indicating overall decay) across a physiological range of Womersley numbers, a spectral broadening ratio R that compares energy in high- versus low-wavenumber modes exhibits a sharp non-monotonic peak at an intermediate Wo. The authors interpret this as evidence that arterial geometry actively triggers resonant energy transfer to short-wavelength components, reframing geometry as a modulator of spectral complexity rather than a passive dissipator.

Significance. If the central claim holds, the work offers a mechanistic link between age-related geometric remodeling and flow spectral features that could inform vascular diagnostics. The distinction between global stability (G < 0) and local spectral redistribution is conceptually useful, and the identification of a resonant Wo would be a falsifiable prediction. However, the strength of these implications depends on whether the reported transfer arises from nonlinear triadic interactions or from linear modal projection induced by the geometry.

major comments (2)
  1. Abstract and model description: The title and abstract frame the result as a 'resonant spectral cascade,' which presupposes nonlinear energy transfer. The provided text does not state whether the governing equations retain the convective term of the Navier-Stokes equations or whether a linearised formulation is used with geometry entering only through boundary conditions or coordinate mapping. If the latter, the non-monotonic peak in R is consistent with passive redistribution among eigenmodes rather than an active cascade, directly undermining the central claim. Please supply the full set of equations, boundary conditions, and numerical scheme (including any linearisation steps) so that the nature of the transfer can be assessed.
  2. Definition and independence of R: The spectral broadening ratio R is introduced to quantify the result, yet no explicit formula, wavenumber cutoff, or normalisation is given in the abstract. If R is constructed from the same modal decomposition used to compute G, or if it incorporates fitted parameters, the reported peak may be tautological rather than an independent diagnostic of resonance. A concrete expression for R (e.g., integral of energy above a cutoff wavenumber divided by energy below it) is required to evaluate whether the non-monotonicity is robust.
minor comments (2)
  1. The abstract states that the investigation was 'conducted over a physiological range of Womersley numbers' but does not list the specific values or the arterial geometry parameters (e.g., curvature radius, tortuosity amplitude). Adding these would improve reproducibility.
  2. The claim that geometry is 'maximally efficient at generating spectral complexity' would benefit from a brief comparison to the straight-tube Womersley solution to isolate the geometric contribution.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. The points raised help strengthen the clarity of our model formulation and diagnostics. We respond to each major comment below and have made revisions to the manuscript to address them directly.

read point-by-point responses

  1. Referee: Abstract and model description: The title and abstract frame the result as a 'resonant spectral cascade,' which presupposes nonlinear energy transfer. The provided text does not state whether the governing equations retain the convective term of the Navier-Stokes equations or whether a linearised formulation is used with geometry entering only through boundary conditions or coordinate mapping. If the latter, the non-monotonic peak in R is consistent with passive redistribution among eigenmodes rather than an active cascade, directly undermining the central claim. Please supply the full set of equations, boundary conditions, and numerical scheme (including any linearisation steps) so that the nature of the transfer can be assessed.

    Authors: We agree that specifying the precise nature of the governing equations is essential for interpreting the mechanism. Our model solves the linearised incompressible Navier-Stokes equations for perturbations about a base Womersley flow, with arterial geometry entering through a coordinate mapping that renders the wall boundary conditions non-trivial in the computational domain. The convective term is linearised about the base flow, so nonlinear triadic interactions are absent; the observed spectral broadening instead arises from linear coupling among non-orthogonal eigenmodes induced by the geometry. We have added a dedicated Methods subsection that states the full linearised equations, no-slip boundary conditions on the mapped arterial wall, inlet/outlet conditions, and the numerical scheme (Fourier spectral in the axial direction combined with finite-element discretisation in the cross-section). We have also revised the abstract, title phrasing, and discussion to describe the phenomenon as a 'geometry-triggered resonant spectral redistribution' rather than a nonlinear cascade, while preserving the central claim that geometry actively modulates spectral complexity even while global energy decays. This clarification directly addresses the concern and strengthens the mechanistic interpretation. revision: yes

  2. Referee: Definition and independence of R: The spectral broadening ratio R is introduced to quantify the result, yet no explicit formula, wavenumber cutoff, or normalisation is given in the abstract. If R is constructed from the same modal decomposition used to compute G, or if it incorporates fitted parameters, the reported peak may be tautological rather than an independent diagnostic of resonance. A concrete expression for R (e.g., integral of energy above a cutoff wavenumber divided by energy below it) is required to evaluate whether the non-monotonicity is robust.

    Authors: R is constructed as an independent diagnostic from the axial Fourier decomposition of the perturbation velocity field. Explicitly, R = (∫_{|k|>k_c} E(k) dk) / (∫_{|k|≤k_c} E(k) dk), where E(k) is the modal kinetic energy spectrum, and the cutoff k_c = 10 (non-dimensional) is chosen a priori from the characteristic axial scale of the arterial geometry (corresponding to wavelengths shorter than ~0.1 of the domain length). This partitioning is performed after the modal decomposition but is separate from the computation of the global growth rate G, which is obtained from the time derivative of the total L2 energy norm integrated over the entire domain. No parameters in R are fitted to the data. We have inserted the explicit formula, the rationale for k_c, and a robustness check (showing the non-monotonic peak persists for k_c between 5 and 15) into the revised Results section and have added a brief statement in the abstract. These additions confirm that the resonance peak is a genuine feature of the geometry-induced linear mode coupling rather than an artifact of the diagnostic definition. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is model-driven output

full rationale

The paper solves a mathematical model of Womersley flow incorporating arterial geometry over a range of Wo values. It reports a negative global growth rate G < 0 as a direct indicator of overall decay and defines the spectral broadening ratio R explicitly as the ratio of energy in high-wavenumber versus low-wavenumber modes, then computes its non-monotonic peak from the model solutions. These quantities are post-processed diagnostics from the governing equations rather than inputs, fitted parameters, or self-referential definitions. No load-bearing self-citations, uniqueness theorems, or ansatzes are invoked in the abstract or description to force the result. The central claim follows from independent model integration without reducing to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Review is based solely on the abstract; no explicit free parameters, axioms, or invented entities are stated. The model presumably relies on standard incompressible Navier-Stokes assumptions and a geometric perturbation, but these are not detailed.

pith-pipeline@v0.9.0 · 5556 in / 1312 out tokens · 52802 ms · 2026-05-17T23:03:00.694501+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages

  1. [1]

    Physiologic blood flow is turbulent,

    K. M. Saqr, S. Tupin, S. Rashad, T. Endo, K. Niizuma, T. Tominaga, and M. Ohta, “Physiologic blood flow is turbulent,”Scientific Reports, vol. 10, no. 1, 15492, 2020

    work page 2020

  2. [2]

    Numerical study on the energy cascade of pulsatile Newtonian and power-law flow models in an ICA bifurcation,

    S. A. Mahrous, N. A. C. Sidik, and K. M. Saqr, “Numerical study on the energy cascade of pulsatile Newtonian and power-law flow models in an ICA bifurcation,”PLOS ONE, vol. 16, no. 1, e0245775, 2021

    work page 2021

  3. [3]

    Effect of artery curvature on the coronary fractional flow reserve: a computational study,

    N. Freidoonimehr, A. F. Fardin, and M. Molla, “Effect of artery curvature on the coronary fractional flow reserve: a computational study,”Phys. Fluids, 33:031906, 2021

    work page 2021

  4. [4]

    Yamaguchi, G

    R. Yamaguchi, G. Tanaka, N. S. Shafii, K. B. Osman, Y. Shimizu, K. M. Saqr, and M. Ohta, “Characteristic effect of wall elasticity on flow instability and wall shear stress of a full-scale, patient-specific aneurysm model in the middle cerebral artery: An experimental approach,”Journal of Applied Physics, vol. 131, no. 18, 184701, 2022

    work page 2022

  5. [5]

    Effects of wall compliance on multiharmonic pulsatile flow in idealized cerebral aneurysm models: comparative PIV experiments,

    S. Tupin, K. M. Saqr, and M. Ohta, “Effects of wall compliance on multiharmonic pulsatile flow in idealized cerebral aneurysm models: comparative PIV experiments,”Experiments in Fluids, vol. 61, no. 7, 172, 2020. 20

    work page 2020

  6. [6]

    Non-Kolmogorov turbulence in carotid artery stenosis and the impact of carotid stenting on near-wall turbulence,

    K. M. Saqr, K. Kano, S. Rashad, K. Niizuma, Y. Kaku, T. Iwama, and T. Tominaga, “Non-Kolmogorov turbulence in carotid artery stenosis and the impact of carotid stenting on near-wall turbulence,”AIP Advances, vol. 12, no. 1, 015118, 2022

    work page 2022

  7. [7]

    Pulsatile entrance flow in a curved pipe,

    L. Talbot and K. O. Gong, “Pulsatile entrance flow in a curved pipe,”J. Fluid Mech., 127:1–25, 1983

    work page 1983

  8. [8]

    Fully developed pulsatile flow in a curved pipe,

    C. C. Hamakiotes and S. A. Berger, “Fully developed pulsatile flow in a curved pipe,”J. Fluid Mech., 195:23–55, 1988

    work page 1988

  9. [9]

    Pulsating flow in curved pipes: 2nd report, exper- iments,

    M. Sumida, K. Sudou, and T. Takami, “Pulsating flow in curved pipes: 2nd report, exper- iments,”Bull. JSME, 27(234):2714–2721, 1984

    work page 1984

  10. [10]

    The effect of Dean, Reynolds, and Womersley number on the flow in a spherical cavity on a curved round pipe,

    F. Chassagne et al., “The effect of Dean, Reynolds, and Womersley number on the flow in a spherical cavity on a curved round pipe,”Sci. Rep., 11:10605, 2021

    work page 2021

  11. [11]

    A review of fluid–structure interaction: blood flow in arteries,

    Z. A. Saib, H. M. Albarody, and A. R. A. Aziz, “A review of fluid–structure interaction: blood flow in arteries,”Eng. Rep., 7(5):e13141, 2025

    work page 2025

  12. [12]

    Impact of modelling surface roughness in an arterial stenosis model on wall shear stress metrics,

    J. Yi, B. J. Li, and B. Guo, “Impact of modelling surface roughness in an arterial stenosis model on wall shear stress metrics,”Fluids, 7(5):179, 2022

    work page 2022

  13. [13]

    Assessment of surface roughness and blood rheology on local coronary haemodynamics—A multiscale CFD study,

    D. G. Owen et al., “Assessment of surface roughness and blood rheology on local coronary haemodynamics—A multiscale CFD study,”Med. Eng. Phys., 74:105–115, 2020

    work page 2020

  14. [14]

    The Physics and Manipulation of Dean Vortices in Curved Microchannels: A Review,

    Y. Saffar, S. Kashanj, D. S. Nobes, and R. Sabbagh, “The Physics and Manipulation of Dean Vortices in Curved Microchannels: A Review,”Micromachines, 14(12):2202, 2023

    work page 2023

  15. [15]

    Literature Survey for In-Vivo Reynolds and Womersley Numbers of Various Arteries and Implications for Compliant In-Vitro Modelling,

    P. N. Williamson, J. D. Humphrey, and S. Baek, “Literature Survey for In-Vivo Reynolds and Womersley Numbers of Various Arteries and Implications for Compliant In-Vitro Modelling,”Cardiovasc. Eng. Technol., 15:1007–1027, 2024

    work page 2024

  16. [16]

    Non-modal transient growth of disturbances in pulsatile and oscillatory pipe flows,

    D. Xu, J. Peixinho, and R. J. Lingwood, “Non-modal transient growth of disturbances in pulsatile and oscillatory pipe flows,”J. Fluid Mech., 919:A16, 2021

    work page 2021

  17. [17]

    Floquet stability analysis of pulsatile flow in toroidal pipes,

    J. S. Kern, V. Lupi, and D. S. Henningson, “Floquet stability analysis of pulsatile flow in toroidal pipes,”Phys. Rev. Fluids, 9:043906, 2024

    work page 2024

  18. [18]

    Fractional-Order Modeling of Arterial Compliance in Vascular As- sessment: A Comprehensive Review,

    M. A. Bahloul et al., “Fractional-Order Modeling of Arterial Compliance in Vascular As- sessment: A Comprehensive Review,”Front. Cardiovasc. Med., 10:1297830, 2023

    work page 2023

  19. [19]

    Comparative assessment of classical and fractional models for blood flow: A review,

    W. F. W. Azmi, M. Z. Mohamad, and A. Q. Mohamad, “Comparative assessment of classical and fractional models for blood flow: A review,”Comput. Methods Programs Biomed., 258:108342, 2025. 21

    work page 2025

  20. [20]

    Applications of neuro-computing and fractional calculus to blood flow problems: A review,

    A. Ali and S. Das, “Applications of neuro-computing and fractional calculus to blood flow problems: A review,”Comput. Biol. Med., 172:107607, 2024

    work page 2024

  21. [21]

    On non-Kolmogorov turbulence in blood flow and its possible role in mechanobiological stimulation,

    K. M. Saqr and I. F. Zidane, “On non-Kolmogorov turbulence in blood flow and its possible role in mechanobiological stimulation,”Scientific Reports, vol. 12, no. 1, 12213, 2022. 22

    work page 2022