pith. sign in

arxiv: 2604.11776 · v1 · submitted 2026-04-13 · ⚛️ physics.flu-dyn · physics.bio-ph· physics.comp-ph· physics.med-ph

Bicuspid Valve Closure and Backflow Prevention: Role of Leaflet Geometry

B. Kaoui , A. Bou Orm , P. Navet , J. Baish , L.L. Munn This is my paper

Pith reviewed 2026-05-10 15:47 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn physics.bio-phphysics.comp-phphysics.med-ph
keywords bicuspid valveleaflet lengthreflux preventionfluid-structure interactionbackward flowvalve closurelymphatic valvesvenous valves
0
0 comments X

The pith

Leaflet length controls whether bicuspid valves achieve complete closure and block reflux under backward flow.

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

Bicuspid valves with crescent-shaped leaflets prevent backward flow in veins and lymphatic vessels through passive closure driven by fluid forces. The study uses numerical simulations of fluid-structure interaction to track how leaflet length alters the flux through the valve when subjected to steady reverse pressure, similar to gravity. A sharp transition occurs from allowing reflux to achieving full blockage as leaflet length increases, and this critical length shifts with the overall valve shape and the rigidity of the leaflets. The results directly account for experimental observations of leakage in valves that have shorter leaflets, whether due to immaturity or abnormality.

Core claim

The flux through the valve orifice depends on valve length, leaflet length, and leaflet rigidity. Varying only the leaflet length produces a transition from reflux to complete flow blockage, with the threshold value depending strongly on valve shape and stiffness. These parameters were captured numerically to evaluate valve ability to close and prevent reflux, explaining why certain incompetent valves with shorter leaflets permit leakage.

What carries the argument

Numerical fluid-structure interaction model that imposes backward flow and computes the resulting leaflet deformation and orifice flux as leaflet length is varied.

If this is right

  • Increasing leaflet length alone can drive the valve from leaky to fully sealed without changes in material stiffness.
  • Shorter leaflets keep the valve in a reflux regime even when shape and stiffness would otherwise allow closure.
  • Valve shape sets the precise leaflet length at which the flow regime switches from reflux to blockage.
  • The same numerical mapping of length, shape, and stiffness can be used to predict performance of observed abnormal valves.

Where Pith is reading between the lines

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

  • Valve engineering in synthetic devices could prioritize leaflet proportion over material tuning to achieve reliable anti-reflux behavior.
  • Temporary reflux during early development may arise simply because leaflets have not yet reached the critical length relative to vessel size.
  • Testing the same length threshold under pulsatile rather than steady reverse flow would reveal whether the transition remains robust in physiological conditions.

Load-bearing premise

The numerical fluid-structure interaction model accurately captures real leaflet deformation and flow patterns under imposed backward flow conditions similar to gravity.

What would settle it

Measure the actual leaflet lengths in valves that show reflux versus those that achieve closure in controlled experiments with steady backward flow, and check whether the observed transition matches the simulated critical length for given shape and stiffness.

Figures

Figures reproduced from arXiv: 2604.11776 by A. Bou Orm, B. Kaoui, J. Baish, L.L. Munn, P. Navet.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic illustration of the numerical set-up (left [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (top panel) Top view of the resting shape of valves of t [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Reflux Q as a function of leaflet length e for various val [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Effect of the valve length L on valve performance in pre [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. State diagrams summarizing the performance of one-w [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Experimental data points (∆ on the right [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Apparent fluid leakage through the leaflet, arising fr [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Comparison between the steady-state shapes of soft, [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Comparison of computed reflux as a function of valve le [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
read the original abstract

Bicuspid valves with crescent-shaped leaflets are found in lymphatic vessels and veins, where their primary function is to prevent reflux and ensure unidirectional flow toward the heart. These valves are passive, and their functionality emerges spontaneously from a complex interplay between the properties of the valve leaflets and the flow patterns developing within the vessel sinus region surrounding the valve. The main function of the valves is to limit retrograde flow, or reflux, but the optimal valve structure has not been well-characterized. Here we investigate numerically how the length of the leaflets affects the valve efficiency in preventing reflux. The valves are subjected to backward flow, akin to that imposed by gravity. We report the flux through the valve orifice as a function of key parameters: valve length, leaflet length, and leaflet rigidity. We monitor the transition in the flow regime - from reflux to complete flow blockage - by varying only the leaflet length. The transition threshold is found to depend strongly on the valve shape and stiffness. We captured these control parameters numerically to evaluate the ability of the valve to close and prevent reflux. This study allowed us to explain reflux observed experimentally in certain incompetent abnormal and immature valves, particularly those with shorter leaflets.

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

3 major / 2 minor

Summary. The manuscript presents a numerical fluid-structure interaction (FSI) study of bicuspid valves with crescent-shaped leaflets, focusing on how leaflet length affects reflux prevention under imposed backward flow akin to gravity. Simulations vary leaflet length (while holding other parameters fixed) and track flux through the orifice, identifying a transition from leakage to complete blockage; the threshold depends strongly on valve shape and leaflet stiffness. This parametric exploration is used to explain experimental observations of reflux in short-leaflet, abnormal, or immature valves.

Significance. If the numerical results are reliable, the work provides a direct mechanistic link between leaflet geometry/stiffness and passive valve competence in lymphatic and venous systems, derived from solution of the coupled fluid and structural equations without fitted parameters or circular definitions. The approach allows clear isolation of geometric control parameters and could guide interpretation of valve incompetence.

major comments (3)
  1. [Numerical Methods] Numerical Methods section: No mesh convergence or grid-independence study is reported for the FSI simulations. Since the central claim rests on identifying a sharp transition in reflux as a function of leaflet length, it is necessary to show that the reported threshold is insensitive to spatial and temporal discretization.
  2. [Results] Results section: Flux values and the transition threshold are presented without error quantification, uncertainty estimates, or sensitivity analysis to numerical parameters. This undermines confidence that the regime change is physical rather than numerical artifact, particularly given the imposed backward-flow boundary conditions.
  3. [Numerical Methods] Validation and model description: The FSI model is not validated against analytical benchmarks (e.g., simple beam deformation under pressure) or experimental valve-closure data. The assumption that the numerical leaflet deformation and flow patterns accurately represent real bicuspid valves under gravity-like reflux therefore remains untested and is load-bearing for the claimed explanatory power.
minor comments (2)
  1. [Abstract] Abstract: The text states that flux is reported 'as a function of key parameters: valve length, leaflet length, and leaflet rigidity' yet only leaflet length is varied; clarify the distinction and whether 'valve length' is a fixed geometric parameter.
  2. [Results] Figure captions and text: Parameter values (e.g., specific rigidity moduli, vessel dimensions) used in the simulations should be listed explicitly in a table or methods summary to enable reproducibility.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed comments, which have helped us improve the clarity and rigor of our numerical study. We address each major comment point by point below, agreeing where revisions are warranted and providing explanations where we maintain our original approach.

read point-by-point responses

  1. Referee: [Numerical Methods] Numerical Methods section: No mesh convergence or grid-independence study is reported for the FSI simulations. Since the central claim rests on identifying a sharp transition in reflux as a function of leaflet length, it is necessary to show that the reported threshold is insensitive to spatial and temporal discretization.

    Authors: We agree that a mesh convergence study is essential to support the reliability of the transition threshold. In the revised manuscript, we have added a new subsection to the Numerical Methods section that reports results from successively refined meshes (doubling and quadrupling the element count in the fluid and solid domains). The critical leaflet length at which reflux ceases and the associated flux values converge to within 4% between the two finest meshes, confirming that the reported regime transition is insensitive to spatial discretization. Temporal convergence with respect to time-step size is also demonstrated for the same threshold. revision: yes

  2. Referee: [Results] Results section: Flux values and the transition threshold are presented without error quantification, uncertainty estimates, or sensitivity analysis to numerical parameters. This undermines confidence that the regime change is physical rather than numerical artifact, particularly given the imposed backward-flow boundary conditions.

    Authors: We acknowledge the absence of explicit uncertainty quantification in the original submission. The revised Results section now includes a sensitivity analysis in which we vary mesh density, time-step size, and the stiffness parameter by ±10% around the reported values. The transition threshold in leaflet length remains unchanged to within one discrete increment, and time-averaged flux differs by less than 3% across these cases. We also report the standard deviation of flux over multiple simulated cycles to quantify variability under the imposed backward flow. revision: yes

  3. Referee: [Numerical Methods] Validation and model description: The FSI model is not validated against analytical benchmarks (e.g., simple beam deformation under pressure) or experimental valve-closure data. The assumption that the numerical leaflet deformation and flow patterns accurately represent real bicuspid valves under gravity-like reflux therefore remains untested and is load-bearing for the claimed explanatory power.

    Authors: We agree that explicit validation strengthens the manuscript. We have added a validation subsection that compares leaflet tip displacement under a uniform pressure load to the analytical solution for a cantilever beam (agreement within 2%). We also include qualitative comparison of the simulated closure dynamics and sinus flow patterns to published experimental observations of bicuspid valves in veins and lymphatics. Direct quantitative validation against high-resolution experimental data for the exact gravity-driven reflux scenario is not available in the literature; we have therefore noted this as a limitation while emphasizing that the model employs standard, parameter-free FSI formulations. revision: partial

Circularity Check

0 steps flagged

No significant circularity; direct numerical parametric study

full rationale

The paper conducts a fluid-structure interaction simulation of bicuspid valve leaflets under imposed backward flow. Leaflet length is varied as an independent input parameter while holding other quantities fixed; the resulting orifice flux and flow-regime transition are computed outputs. No parameter is fitted to a data subset and then re-used as a prediction of a closely related quantity, no self-citation supplies a load-bearing uniqueness theorem or ansatz, and the central claims (threshold dependence on shape and stiffness) follow directly from the numerical solution of the governing equations rather than from any definitional equivalence or renaming of known results.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard continuum mechanics assumptions for incompressible flow and elastic leaflets, with leaflet length and rigidity treated as free parameters varied across simulations.

free parameters (2)
  • leaflet length
    Varied to identify the transition threshold from reflux to complete blockage.
  • leaflet rigidity
    Varied as a control parameter affecting closure efficiency.
axioms (2)
  • standard math Incompressible Navier-Stokes equations govern the fluid flow
    Standard assumption for low-Reynolds-number biofluid flow in vessels.
  • domain assumption Leaflets behave as elastic structures under fluid loading
    Required for modeling passive deformation and closure.

pith-pipeline@v0.9.0 · 5530 in / 1191 out tokens · 76681 ms · 2026-05-10T15:47:05.946506+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

33 extracted references · 33 canonical work pages

  1. [1]

    Scallan, S.D

    J.P. Scallan, S.D. Zawieja, J.A. Castorena-Gonzalez, M .J. Davis, Lymphatic pumping: mechanics, mechanisms and ma l- function, The Journal of Physiology 594.20, pp 5749–5768 (2 016)

    work page

  2. [2]

    Padera, E.F.J

    T.P. Padera, E.F.J. Meijer, L.L. Munn, The lymphatic sys tem in disease processes and cancer progression, Annual Rev iew of Biomedical Engineering 18, 125–58 (2016)

    work page 2016

  3. [3]

    Moore Jr., C.D

    J.E. Moore Jr., C.D. Bertram, Lymphatic system flows, Ann ual Review of Fluid Mechanics 50, 459–82 (2018)

    work page 2018

  4. [4]

    Kunert, J.W

    C. Kunert, J.W. Baish, S. Liao, T.P. Padera, L.L. Munn, Me chanobiological oscillators control lymph flow, Proc. Natl . Acad. Sci. U.S.A. 112 (35) 10938-10943 (2015),

    work page 2015

  5. [5]

    Munn, Mechanobiology of lymphatic contractions, S eminars in Cell & Developmental Biology 38, 67-74 (2015)

    L.L. Munn, Mechanobiology of lymphatic contractions, S eminars in Cell & Developmental Biology 38, 67-74 (2015)

    work page 2015

  6. [6]

    Baish , C

    J.W. Baish , C. Kunert, T.P. Padera, L.L. Munn, Synchroni zation and random triggering of lymphatic vessel contracti ons. PLOS Computational Biology 12(12): e1005231 (2016)

    work page 2016

  7. [7]

    Munger, M.J

    S.J. Munger, M.J. Davis, A.M. Simon, Defective lymphati c valve development and chylothorax in mice with a lymphatic - specific deletion of Connexin43, Developmental Biology 421 , 204–218 (2017)

    work page 2017

  8. [8]

    Davis, S.D

    M.J. Davis, S.D. Zawieja, Y. Yang, Developmental progre ssion of lymphatic valve morphology and function, Frontier s in Cell and Developmental Biology 12:1331291 (2024)

    work page 2024

  9. [9]

    Margaris and R.A

    K.N. Margaris and R.A. Black, Modelling the Lymphatic Sy stem: Challenges and Opportunities, J. R. Soc. Interface 9, 601–612 (2012)

    work page 2012

  10. [10]

    D. J. Watson, I. Sazonov, D. C. Zawieja, J. E. Moore Jr., R . van Loon, Integrated Geometric and Mechanical Analysis of an Image-Based Lymphatic Valve, J. Biomech. 64, 172–179 (20 17)

    work page

  11. [11]

    J. T. Wilson, R. van Loon, W. Wang, D. C. Zawieja, et J. E. M oore, Determining the combined effect of the lymphatic valve leaflets and sinus on resistance to forward flow, Journa l of Biomechanics 48(13), 3584-3590 (2015)

    work page 2015

  12. [12]

    J. T. Wilson, L. T. Edgar, S. Prabhakar, M. Horner, R. van Loon, J. E. Moore Jr., A fully Coupled Fluid-Structure Interaction Model of the Secondary Lymphatic Valve, Comput . Methods Biomech. Biomed. Eng. 21(16), 813-823 (2018)

    work page 2018

  13. [13]

    C. D. Bertram, Modelling Secondary Lymphatic Valves wi th a Flexible Vessel Wall: How Geometry and Material Propert ies Combine to Provide Function, Biomech. Model. Mechanobiol 1 9, 2081–2098 (2020)

    work page 2081

  14. [14]

    C. D. Bertram, M. J. Davis, An Enhanced 3D Model of Intrav ascular Lymphatic Valves to Assess Leaflet Apposition and Transvalvular Differences in Wall Distensibility, Biology 12, 379 (2023)

    work page 2023

  15. [15]

    K. T. Wolf, J. B. Dixon, A. Alexeev, Fluid Pumping of Peri staltic Vessel Fitted with Elastic Valves, J. Fluid Mech. 91 8, A28 (2021)

    work page 2021

  16. [16]

    Ballard, K

    M. Ballard, K. T. Wolf, Z. Ne, P. Cheikh, J. B. Dixon, et A. Alexeev, Probing the effect of morphology on lymphatic valve dynamic function, Biomechanics and Modeling in Mechanobio logy 17(4), 1031-1044 (2018)

    work page 2018

  17. [17]

    Succi, The lattice Boltzmann equation for fluid dynam ics and beyond (Oxford University Press, Oxford, UK, 2001)

    S. Succi, The lattice Boltzmann equation for fluid dynam ics and beyond (Oxford University Press, Oxford, UK, 2001)

    work page 2001

  18. [18]

    Ostoja-Starzewski, Lattice models in micromechani cs, Applied Mechanics Reviews 55(1), 35-60 (2002)

    M. Ostoja-Starzewski, Lattice models in micromechani cs, Applied Mechanics Reviews 55(1), 35-60 (2002)

    work page 2002

  19. [19]

    Peskin, The Immersed Boundary Method, Acta Numeri ca 11, 479 (2002)

    C.S. Peskin, The Immersed Boundary Method, Acta Numeri ca 11, 479 (2002)

    work page 2002

  20. [20]

    Hofferberth, M.Y

    S.C. Hofferberth, M.Y. Saeed, L. Tomholt, M.C. Fernande s, C.J. Payne, K. Price, G.R. Marx, J. J. Esch, D.W. Brown, J. Brown, P.E. Hammer, R.W. Bianco, J.C. Weaver, E.R. Edelma n, P.J. del Nido, A geometrically adaptable heart valve replacement, Science Translational Medicine 12, eaay4006 (2020)

    work page 2020

  21. [21]

    Hammer, E.G

    P.E. Hammer, E.G. Roberts, S.M. Emani, P.J. del Nido, Su rgical reconstruction of semilunar valves in the growing ch ild: Should we mimic the venous valve? A simulation study, The Jou rnal of Thoracic and Cardiovascular Surgery 153(2), 389 (2017)

    work page 2017

  22. [22]

    H. Li, Y. Mei, N. Maimon, T.P. Padera, J.W. Baish, L.L. Mu nn, The effects of valve leaflet mechanics on lymphatic pumping assessed using numerical simulations, Scientific R eports 9, 10649 (2019)

    work page 2019

  23. [23]

    Lurie, Robert L

    F. Lurie, Robert L. Kistner, B. Eklof, The mechanism of v enous valve closure in normal physiologic conditions, Jour nal of Vascular Surgery 35(4), 713-717 (2002)

    work page 2002

  24. [24]

    Lurie, R.L

    F. Lurie, R.L. Kistner, B. Eklof, D. Kessler, Mechanism of venous valve closure and role of the valve in circulation: A new concept, Journal of Vascular Surgery 38 (5), 955-961 (2003)

    work page 2003

  25. [25]

    F. C. Wells, The Heart of Leonardo (Springer-Verlag, Lo ndon, 2013)

    work page 2013

  26. [26]

    X. Provot, Deformation constraints in a mass-spring mo del to describe rigid cloth behaviour, Proceedings of Graph ics Interface ’95: Qu´ ebec, Qu´ ebec, Canada, 17 - 19 May 1995, 147-154

    work page 1995

  27. [27]

    Kaoui, Multiphysics simulations in biomedical: flui d–structure interaction coupled with mass transfer, Compt es Rendus

    B. Kaoui, Multiphysics simulations in biomedical: flui d–structure interaction coupled with mass transfer, Compt es Rendus. M´ ecanique 353, 779-790 (2025)

    work page 2025

  28. [28]

    Kaoui, Valve 3D, Gitlab repository (2026), https://gitlab.utc.fr/kaouibad/valve_3d

    B. Kaoui, Valve 3D, Gitlab repository (2026), https://gitlab.utc.fr/kaouibad/valve_3d

    work page 2026

  29. [29]

    Sabine et al

    A. Sabine et al. , Mechanotransduction, PROX1, and FOXC2 Cooperate to Contr ol Connexin37 and Calcineurin during Lymphatic-Valve Formation, Developmental Cell 22, 430–44 5 (2012)

    work page 2012

  30. [30]

    Bielinski, L

    C. Bielinski, L. Xia, G. Helbecque, B. Kaoui, Mass trans fer from a sheared spherical rigid capsule, Physics of Fluid s 34, 031902 (2022)

    work page 2022

  31. [31]

    Kr¨ uger, B

    T. Kr¨ uger, B. Kaoui, J. Harting Interplay of inertia an d deformability on rheological properties of a suspension o f capsules Journal of Fluid Mechanics 751, 725-745 (2014)

    work page 2014

  32. [32]

    M. M. Dupin, I. Halliday, C. M. Care, L. Alboul, L. L. Munn , Modeling the flow of dense suspensions of deformable particles in three dimensions, Physical Review E 75, 066707 (2007)

    work page 2007

  33. [33]

    Kaoui, J

    B. Kaoui, J. Harting, C. Misbah, Two-dimensional vesic le dynamics under shear flow: effect of confinement, Physical Review E 83, 066319 (2011)

    work page 2011