Large-eddy simulation of the FDA benchmark blood pump: validation against experiments and implications for turbulent flow mechanisms
Pith reviewed 2026-05-10 08:14 UTC · model grok-4.3
The pith
Large-eddy simulation with transient rotor-stator coupling matches experimental velocity fields in blood pumps more closely than Reynolds-averaged methods, especially in turbulent diffuser regions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
LES with transient sliding-interface coupling produces velocity predictions that align more closely with particle image velocimetry data than either multiple-reference-frame or sliding-interface RANS formulations, particularly in the diffuser region. A mesh of approximately 80 million cells satisfies three complementary resolution metrics for a well-resolved LES. The validated simulations then reveal the spatial distribution of vortical structures, turbulent kinetic energy, and velocity spectra that characterize the unsteady internal flow.
What carries the argument
Large-eddy simulation (LES) with transient sliding-interface coupling between rotor and stator regions.
If this is right
- Scale-resolving transient methods are required to capture the highly unsteady turbulence that dominates flow inside ventricular assist devices.
- RANS approaches remain unreliable in diffuser sections where intermittency and wall interactions are strong.
- Validated LES supplies a practical route to detailed analysis of vortical structures and turbulent kinetic energy for hemocompatibility assessment.
- Mesh resolutions near 80 million cells suffice to reach a well-resolved regime for this class of pump geometry.
Where Pith is reading between the lines
- Design iterations for blood pumps could incorporate LES early to identify and reduce regions of high shear or recirculation before physical testing.
- The same transient coupling approach may transfer to other rotating fluid machines that combine blades and stationary passages.
- Accurate energy spectra from these simulations could support Lagrangian tracking of blood cells to estimate exposure to damaging flow conditions.
Load-bearing premise
An 80-million-cell mesh together with the selected subgrid-scale model and transient interface treatment reproduces the physical unsteady turbulence without dominant numerical artifacts or boundary-condition errors.
What would settle it
A simulation on a mesh at least twice as fine, or a direct numerical simulation of the same geometry and conditions, that produces velocity fields or spectra differing substantially from the reported LES results in the diffuser would indicate the current mesh and modeling choices are inadequate.
Figures
read the original abstract
This study presents a systematic validation and comparative assessment of computational fluid dynamics (CFD) strategies for centrifugal blood pump simulations using the U.S. Food and Drug Administration benchmark model. A scale-resolving large eddy simulation (LES) with transient sliding-interface (SI) coupling is evaluated and compared against Reynolds-averaged Navier-Stokes (RANS) approaches employing both multiple reference frame and SI formulations. Numerical predictions are validated through direct comparison with particle image velocimetry measurements under two representative operating conditions. The results indicate that LES with transient rotor-stator coupling achieves consistently improved agreement with experimental velocity fields compared with RANS-based methods, particularly in the diffuser region where strong intermittency and wall-bounded turbulence are present. In contrast, RANS-based approaches exhibit noticeable discrepancies in these regions. A mesh sensitivity study and an assessment of temporal averaging effects are conducted for LES. The quality of the LES results is further quantified using three complementary metrics, demonstrating that a mesh resolution of approximately 80 million cells achieves a well-resolved LES regime. Building on the validated scale-resolving simulations, detailed analyses of vortical structures, turbulent kinetic energy distributions, and velocity energy spectra are performed to characterize the internal flow physics of the pump. This study demonstrates that scale-resolving, transient simulation approaches are essential for accurately capturing the highly unsteady, turbulence-dominated flow features in ventricular assist devices and provides practical guidance for future high-fidelity hemodynamic and hemocompatibility studies.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents large-eddy simulations (LES) with transient sliding-interface (SI) coupling for the FDA benchmark centrifugal blood pump, comparing them to RANS simulations using both multiple reference frame (MRF) and SI formulations. Numerical results are validated directly against particle image velocimetry (PIV) measurements for two operating conditions. The authors claim that LES achieves consistently improved agreement with experimental velocity fields relative to RANS, especially in the diffuser region with strong intermittency and wall-bounded turbulence. A mesh sensitivity study and temporal averaging assessment are performed, and three complementary metrics are used to demonstrate that a mesh of approximately 80 million cells reaches a well-resolved LES regime. The work further analyzes vortical structures, turbulent kinetic energy distributions, and velocity energy spectra to characterize the pump's internal flow physics.
Significance. If the validation holds, the work is significant for biomedical fluid dynamics because it supplies direct experimental evidence that scale-resolving, transient LES is required to capture the unsteady turbulence that governs hemocompatibility in ventricular assist devices. The direct PIV comparison, mesh sensitivity study, and turbulence-mechanism analysis constitute concrete strengths that can guide future high-fidelity hemodynamic modeling.
major comments (1)
- The central claim that the ~80-million-cell mesh achieves a well-resolved LES regime (and that observed improvements over RANS therefore reflect faithful turbulence resolution) rests on three complementary quality metrics whose numerical values, thresholds, and verification against a finer grid or DNS are not reported. Without these quantities (e.g., resolved TKE fraction, Pope criterion, or spectral cutoff) or a demonstration that the transient SI does not inject spurious energy at the rotor frequency, it remains possible that the reported gains arise from reduced numerical dissipation rather than physical fidelity. This issue is load-bearing for the paper's recommendation of LES over RANS.
minor comments (2)
- The abstract would be strengthened by inclusion of quantitative error statistics (mean or maximum velocity discrepancies) and explicit locations of remaining discrepancies between LES and PIV data.
- The methods section should expand the description of the subgrid-scale model formulation and the precise boundary-condition treatment at the sliding interface to support reproducibility.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed review. The major comment identifies a key area where additional quantification and clarification are warranted to strengthen the LES resolution claims. We address the point below and have revised the manuscript to incorporate explicit numerical values, thresholds, and expanded discussion.
read point-by-point responses
-
Referee: The central claim that the ~80-million-cell mesh achieves a well-resolved LES regime (and that observed improvements over RANS therefore reflect faithful turbulence resolution) rests on three complementary quality metrics whose numerical values, thresholds, and verification against a finer grid or DNS are not reported. Without these quantities (e.g., resolved TKE fraction, Pope criterion, or spectral cutoff) or a demonstration that the transient SI does not inject spurious energy at the rotor frequency, it remains possible that the reported gains arise from reduced numerical dissipation rather than physical fidelity. This issue is load-bearing for the paper's recommendation of LES over RANS.
Authors: We agree that explicit numerical values and thresholds for the three quality metrics (resolved TKE fraction, Pope criterion, and spectral cutoff) should be reported to substantiate the well-resolved LES claim. In the revised manuscript we will add a dedicated table listing the exact values obtained on the 80-million-cell mesh together with the acceptance thresholds applied (e.g., resolved TKE fraction > 80 %, Pope criterion < 0.1, and inertial-range spectral cutoff consistent with -5/3 scaling). A mesh-sensitivity study comparing 40 M, 80 M, and 120 M cell resolutions already demonstrates convergence of mean velocity and TKE between the two finest meshes; we will expand this section to include the quantitative metric values at each resolution. A full DNS verification remains computationally prohibitive for the present Reynolds number and geometry (estimated cell count > 10^12), but the observed mesh convergence provides supporting evidence. Regarding the transient sliding-interface coupling, we have re-examined the velocity energy spectra in the diffuser and found no spurious peaks at the rotor blade-passing frequency; the spectra exhibit the expected inertial-range decay without artificial energy injection. We will add an explicit paragraph and supplementary spectra comparing SI and frozen-rotor results to demonstrate that the observed improvements over RANS arise from physical turbulence resolution rather than reduced numerical dissipation. These additions directly address the load-bearing concern for the LES-over-RANS recommendation. revision: yes
- A complete DNS verification is not feasible owing to prohibitive computational cost for this high-Re, complex-geometry flow.
Circularity Check
No circularity: central validation rests on external PIV experiments and mesh-sensitivity assessment
full rationale
The manuscript's primary claims are established by direct numerical-experimental comparison of velocity fields under two operating conditions, using independent particle-image-velocimetry data as the benchmark. Mesh resolution at ~80 million cells is supported by a sensitivity study plus three (unspecified) complementary quality metrics, none of which are shown to be fitted to or defined by the target flow quantities. No equations, parameter-fitting steps, or self-citations are invoked to derive the reported agreement; the simulation outputs remain falsifiable against the external measurements. Consequently the derivation chain contains no self-definitional, fitted-input, or self-citation reductions.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Incompressible Navier-Stokes equations govern the flow.
- domain assumption The chosen subgrid-scale model in LES sufficiently represents unresolved scales in the pump flow.
Reference graph
Works this paper leans on
-
[1]
G. Loor and G. Gonzalez-Stawinski. Pulsatile vs. continuous flow in ventricular assist device therapy.Best Practice & Research Clinical Anaesthesiology, 26(2):105–115, 2012
work page 2012
-
[2]
T. Miyamoto, K. Fukamachi, and J.H. Karimov. Chapter 6 - Continuous-Flow Ventricular Assist Devices. In J.H. Karimov, K.Fukamachi, and M.Gillinov, editors,Advances in Cardiovascular Technology, pages 79–119. Academic Press, 2022
work page 2022
- [3]
-
[4]
Z. Chen, A. Sun, H. Wang, Y. Fan, and X. Deng. Non-physiological shear stress-induced blood damage in ventricular assist device.Medicine in Novel Technology and Devices, 3:100024, 2019
work page 2019
-
[5]
C. O’Brien, J. Monteagudo, C. Schad, E. Cheung, and W. Middlesworth. Centrifugal pumps and hemolysis in pediatric extracorporeal membrane oxygenation (ECMO) patients: An analysis of Extracorporeal Life Support Organization (ELSO) registry data.Journal of Pediatric Surgery, 52(6):975–978, 2017
work page 2017
-
[6]
K.H. Fraser, M.E. Taskin, B.P. Griffith, and Z.J. Wu. The use of computational fluid dynamics in the development of ventricular assist devices.Medical engineering & physics, 33(3):263–280, 2011
work page 2011
-
[7]
Collected nondimensional performance of rotary dynamic blood pump
J. Wu. Letter to the editor: A possible major mistake in the paper entitled “Collected nondimensional performance of rotary dynamic blood pump”: Smith WA, Allaire P, Antaki J, Butler kc, Kerkhoffs W, Kink T, Loree H, Reul H. ASAIO Journal 50 : 25-32, 2004.ASAIO Journal, 53(2):255–256, 2007
work page 2004
-
[8]
B. Torner, L. Konnigk, S. Hallier, J. Kumar, M. Witte, and F.-H. Wurm. Large eddy simulation in a rotary blood pump: Viscous shear stress computation and comparison with unsteady Reynolds-averaged Navier- Stokes simulation.The International Journal of Artificial Organs, 41:752–763, 2018
work page 2018
-
[9]
US Food and Drug Administration. Benchmark dataset for validating computational fluid dynamic (CFD) simulation of blood flow through generalized medical device geometries.https://cdrh-rst.fda.gov/ benchmark-dataset-validating-computational-fluid-dynamic-cfd-simulation-blood-flow-through. Accessed: 2025-07-19
work page 2025
-
[10]
P. Hariharan, K.I. Aycock, M. Buesen, S.W. Day, B.C. Good, L.H. Herbertson, U. Steinseifer, K.B. Manning, B.A. Craven, and R.A. Malinauskas. Inter-Laboratory Characterization of the Velocity Field in the FDA Blood Pump Model Using Particle Image Velocimetry (PIV).Cardiovascular Engineering and Technology, 9:623–640, 2018
work page 2018
-
[11]
R. Malinauskas, P. Hariharan, S. Day, L. Herbertson, M. B¨ usen, K. Aycock, B. Good, S. Deutsch, K. Manning, and B. Craven. FDA Benchmark Medical Device Flow Models for CFD Validation.ASAIO Journal, 63(2):317– 327, 2017
work page 2017
-
[12]
S.V. Ponnaluri, P.Hariharan, L.H. Herbertson, K.B. Manning, R.A. Malinauskas, and B.A. Craven. Results of the Interlaboratory Computational Fluid Dynamics Study of the FDA Benchmark Blood Pump.Annals of Biomedical Engineering, 51:253–269, 2023
work page 2023
-
[13]
B. Good and K. Manning. Computational Modeling of the Food and Drug Administration’s Benchmark Centrifugal Blood Pump.Artificial Organs, 44:E263–E276, 2020
work page 2020
-
[14]
S. Gross-Hardt, S. Sonntag, F. Boehning, T. Schmitz-Rode, and T. Kaufmann. Crucial Aspects for Using Computational Fluid Dynamics as a Predictive Evaluation Tool for Blood Pumps.ASAIO Journal, 65(8):864– 873, 2019
work page 2019
-
[15]
C.S. Semenzin, B. Simpson, S.D. Gregory, and G. Tansley. Validated Guidelines for Simulating Centrifugal Blood Pumps.Cardiovascular Engineering and Technology, 12:273–285, 2021. 18
work page 2021
-
[16]
C. Miccoli, B. Collins, A. Scardigli, and F. Gallizio. Robust shape optimization of the FDA blood pump. Meccanica, 2024
work page 2024
-
[17]
J. Huo, P. Wu, L. Zhang, and W. Wu. Large eddy simulation as a fast and accurate engineering approach for the simulation of rotary blood pumps.The International Journal of Artificial Organs, 44(11):887–899, 2021
work page 2021
-
[18]
A. Gil, R. Navarro, P. Quintero, and A. Mares. Hemocompatibility and hemodynamic comparison of two centrifugal LVADs: HVAD and HeartMate3.Biomechanics and Modeling in Mechanobiology, 22:871–883, 2023
work page 2023
-
[19]
D. Han, J.L. Leibowitz, L. Han, S. Wang, G. He, B.P. Griffith, and Z.J. Wu. Computational fluid dynam- ics analysis and experimental hemolytic performance of three clinical centrifugal blood pumps: Revolution, Rotaflow and Centrimag.Medicine in Novel Technology and Devices, 15:100153, 2022
work page 2022
- [20]
-
[21]
I.B. Celik, Z.N. Cehreli, and I. Yavuz. Index of Resolution Quality for Large Eddy Simulations.Journal of Fluids Engineering, 127:949–958, 2005
work page 2005
-
[22]
P. Gousseau, B. Blocken, and G.J.F. van Heijst. Quality assessment of Large-Eddy Simulation of wind flow around a high-rise building: Validation and solution verification.Computers & Fluids, 79:120–133, 2013
work page 2013
-
[23]
M. Klein. An Attempt to Assess the Quality of Large Eddy Simulations in the Context of Implicit Filtering. Flow, Turbulence and Combustion, 75:131–147, 2005
work page 2005
-
[24]
F. Nicoud and F. Ducros. Subgrid-Scale Stress Modelling Based on the Square of the Velocity Gradient Tensor. Flow, Turbulence and Combustion, 62:183–200, 1999
work page 1999
- [25]
-
[26]
P. Wu, J. Huo, Z. Zhang, and C. Wang. The influence of non-conformal grid interfaces on the results of large eddy simulation of centrifugal blood pumps.Artificial Organs, 46(9):1804–1816, 2022
work page 2022
-
[27]
L. Pauli, J.W. Both, and M. Behr. Stabilized finite element method for flows with multiple reference frames. International Journal for Numerical Methods in Fluids, 78(11):657–669, 2015
work page 2015
-
[28]
M. Zadravecand S. Basic and M. Hriberˇ sek. The influence of rotating domain size in a rotating frame of reference approach for simulation of rotating impeller in a mixing vessel.Journal of Engineering Science and Technology, 2:126–138, 2007
work page 2007
-
[29]
P. Wu, K.J. Zhang, W.J. Xiang, and G.T. Du. Turbulent flow field in maglev centrifugal blood pumps of CH- VAD and HeartMate III: secondary flow and its effects on pump performance.Biomechanics and Modeling in Mechanobiology, 23(5):1571–1589, 2024
work page 2024
-
[30]
P. Wu, L. Zhang, Q. Gao, and W. Dai. Effect of turbulent inlet conditions on the prediction of flow field and hemolysis in the FDA ideal medical device.Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, 235:391–401, 2019
work page 2019
- [31]
-
[32]
S.B. Pope. Ten questions concerning the large-eddy simulation of turbulent flows.New Journal of Physics, 6(1):35, 2004
work page 2004
-
[33]
G. Matheou and D. Chung. Large-eddy simulation of stratified turbulence. Part II: Application of the stretched- vortex model to the atmospheric boundary layer.Journal of the Atmospheric Sciences, 71(12):4439–4460, 2014
work page 2014
-
[34]
C. Catellani.Development and Assessment of Large Eddy Simulation Methodology for Internal Combustion Engines. PhD thesis, Universit` a di Bologna, 2016. 19
work page 2016
-
[35]
B.J. Geurts and J. Fr¨ ohlich. A framework for predicting accuracy limitations in large-eddy simulation.Physics of Fluids, 14(6):L41–L44, 2002
work page 2002
-
[36]
R. Broglia, A. Pascarelli, and U. Piomelli. Large-eddy simulations of ducts with a free surface.Journal of Fluid Mechanics, 484:223–253, 2003
work page 2003
-
[37]
A.D. Scillitoe, P.G. Tucker, and P. Adami. Evaluation of RANS and ZDES Methods for the Prediction of Three-Dimensional Separation in Axial Flow Compressors. InASME Turbo Expo 2015: Turbine Technical Conference and Exposition, 2015
work page 2015
-
[38]
Y. Liu, N. Xie, Y. Tang, and Y. Zhang. Investigation of hemocompatibility and vortical structures for a centrifugal blood pump based on large-eddy simulation.Physics of Fluids, 34(11):115111, 2022
work page 2022
- [39]
-
[40]
J. C. R. Hunt, A. A. Wray, and P. Moin. Eddies, streams, and convergence zones in turbulent flows. InPro- ceedings of the 1988 Summer Program, pages 193–208. Center for Turbulence Research, NASA Ames/Stanford University, 1988
work page 1988
-
[41]
R. Mohammadi, M.S. Karimi, M. Raisee, and M. Sharbatdar. Probabilistic CFD analysis on the flow field and performance of the FDA centrifugal blood pump.Applied Mathematical Modelling, 109:555–577, 2022
work page 2022
-
[42]
R. Sahebi-Kuzeh Kanan, H. Niroomand-Oscuii, H. Badri Ghavifekr, and F. Ghalichi. Design and prototyping of a magnetically and hydrodynamically suspended blood pump via multiobjective optimization.Scientific Reports, 15(1):26171, 2025
work page 2025
-
[43]
M. Kameneva, K. Kono, B. Repko, J. Antaki, and M. Umezu. Effects of Turbulent Stresses upon Mechanical Hemolysis: Experimental and Computational Analysis.ASAIO journal, 50:418–423, 2004
work page 2004
- [44]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.