Results of FDA’s First Interlaboratory Computational Study of a Nozzle with a Sudden Contraction and Conical Diffuser

The U.S. Food and Drug Administration recently hosted an interlaboratory study to assess the suitability and methodology of computational fluid dynamics (CFD) for demonstrating medical device safety in regulatory submissions. The benchmark study was performed in a generic medical device consisting of a 0.012 m diameter cylindrical nozzle with a sudden contraction and 10° conical diffuser, on either side of a 0.04 m long, 0.004 m diameter throat. Results from 28 simulations from around the world were compared to planar particle image velocimetry (PIV) measurements performed at three laboratories. Five flow rates were chosen that produced laminar, transitional, and turbulent regimes. In general, the CFD results showed modest agreement in global and local flow behaviors. However, all CFD data sets contained wide degrees of velocity variation in comparison to experiment and with each other, much of which could be attributed to the turbulence models used. Some of the velocity discrepancies in the sixteen three dimensional (3D) simulations resulted from substantial flow asymmetries within and downstream of the conical diffuser. Large differences in velocity symmetry were found even when using the same code and many of the same simulation parameters. In contrast, no velocity asymmetries of the same magnitude were observed in any of the planar PIV experiments; relatively minor asymmetries in one experiment were linked to slight asymmetries in the velocity profile at the inlet to the nozzle. CFD predictions of peak wall shear stress at the sudden contraction (normalized to that 0.015 m downstream) varied over two orders of magnitude, with variations of greater than 10 times even when the same software and turbulence model were used. This degree of shear stress variation will necessarily propagate through calculations of blood damage in medical devices. We conclude that CFD simulations used in qualifying medical devices for regulatory purposes need to be conducted following the best practices available, and that targeted experimental validation is essential. The results of this interlaboratory study are freely available in an internet repository (https://fdacfd.nci.nih.gov). We encourage the use of this model in further studies, and support the development of additional benchmarks, better modeling techniques, and consensus standards and guidelines for using CFD in the evaluation of medical devices.

[1]  Pavlos Vlachos,et al.  In Vitro, Time-Resolved PIV Comparison of the Effect of Stent Design on Wall Shear Stress , 2009, Annals of Biomedical Engineering.

[2]  D. K. Walters,et al.  Laminar, Turbulent, and Transitional Simulations in Benchmark Cases with Cardiovascular Device Features , 2013 .

[3]  Xiao-Chen Yang,et al.  Computational fluid dynamics-based hydraulic and hemolytic analyses of a novel left ventricular assist blood pump. , 2011, Artificial organs.

[4]  Pavlos P. Vlachos,et al.  Robust wall gradient estimation using radial basis functions and proper orthogonal decomposition (POD) for particle image velocimetry (PIV) measured fields , 2009 .

[5]  David D. Apsley,et al.  Advanced Turbulence Modelling of Separated Flow in a Diffuser , 2000 .

[6]  T. Trucano,et al.  Verification, Validation, and Predictive Capability in Computational Engineering and Physics , 2004 .

[7]  Yubo Fan,et al.  A new way to reduce flow disturbance in endovascular stents: a numerical study. , 2011, Artificial organs.

[8]  Kerem Pekkan,et al.  Right ventricular outflow tract reconstruction with bicuspid valved polytetrafluoroethylene conduit. , 2011, The Annals of thoracic surgery.

[9]  Richard B. Medvitz,et al.  Computational fluid dynamics design and analysis of a passively suspended Tesla pump left ventricular assist device. , 2011, Artificial organs.

[10]  Gábor Janiga,et al.  Numerical Prediction of Hemolysis Based on Computational Fluid Dynamics , 2011 .

[11]  Christian J. Kähler,et al.  Wall-shear-stress and near-wall turbulence measurements up to single pixel resolution by means of long-distance micro-PIV , 2006 .

[12]  Gianluca Iaccarino,et al.  Predictions of a turbulent separated flow using commercial CFD codes , 2001 .

[13]  Kameswararao Anupindi,et al.  Large Eddy Simulation of FDA’s Idealized Medical Device , 2013, Cardiovascular engineering and technology.

[14]  Steven Deutsch,et al.  Assessment of CFD Performance in Simulations of an Idealized Medical Device: Results of FDA’s First Computational Interlaboratory Study , 2012 .

[15]  Ulrich Steinseifer,et al.  Transient, Three-Dimensional Flow Field Simulation through a Mechanical, Trileaflet Heart Valve Prosthesis , 2011, ASAIO journal.

[16]  F. Durst,et al.  Asymmetric flows and instabilities in symmetric ducts with sudden expansions , 1978, Journal of Fluid Mechanics.

[17]  Linden Alyssa Down Computational investigations of red blood cell mechanical trauma and of diseased renal artery hemodynamics , 2011 .

[18]  Mitsuo Umezu,et al.  Fluid Dynamic Investigation of Mechanical Blood Hemolysis , 1996 .

[19]  Majid Nabavi Three-dimensional asymmetric flow through a planar diffuser: Effects of divergence angle, Reynolds number and aspect ratio , 2010 .

[20]  Ali Etebari,et al.  Improvements on the accuracy of derivative estimation from DPIV velocity measurements , 2005 .

[21]  Bruce R. Munson,et al.  Viscous buckling of thin fluid layers , 1981 .

[22]  小柳 仁,et al.  Heart replacement : artificial heart 4 , 1993 .

[23]  L Goubergrits,et al.  Flow dynamics of a novel counterpulsation device characterized by CFD and PIV modeling. , 2011, Medical engineering & physics.

[24]  David H. Frakes,et al.  Physics-Driven CFD Modeling of Complex Anatomical Cardiovascular Flows—A TCPC Case Study , 2005, Annals of Biomedical Engineering.

[25]  Mitsuo Umezu,et al.  Effects of Turbulent Stresses upon Mechanical Hemolysis: Experimental and Computational Analysis , 2004, ASAIO journal.

[26]  Pavlos Vlachos,et al.  In Vitro Comparison of the Effect of Stent Configuration on Wall Shear Stress Using Time-resolved Particle Image Velocimetry , 2010, Annals of Biomedical Engineering.

[27]  C. K. Chong,et al.  Rotational invariance in the three-dimensional lattice Boltzmann method is dependent on the choice of lattice , 2011, J. Comput. Phys..

[28]  Mehrdad Massoudi,et al.  On the Representation of Turbulent Stresses for Computing Blood Damage. , 2010, International journal of engineering science.

[29]  Shmuel Einav,et al.  Device Thrombogenicity Emulator (DTE)--design optimization methodology for cardiovascular devices: a study in two bileaflet MHV designs. , 2010, Journal of biomechanics.

[30]  Alessandra Nigro,et al.  Numerical and Experimental Flow Analysis of the Wang-Zwische Double-Lumen Cannula , 2011, ASAIO journal.

[31]  F. Durst,et al.  Low Reynolds number flow over a plane symmetric sudden expansion , 1974, Journal of Fluid Mechanics.

[32]  J. Westerweel,et al.  Single-pixel resolution ensemble correlation for micro-PIV applications , 2004 .

[33]  D. K. Walters,et al.  A dynamic hybrid Reynolds-averaged Navier Stokes–Large eddy simulation modeling framework , 2012 .

[34]  D. Giddens,et al.  Local hemodynamics affect monocytic cell adhesion to a three-dimensional flow model coated with E-selectin. , 2001, Journal of biomechanics.

[35]  Lois Gold,et al.  A New Way , 1992 .

[36]  Matthew F. Barone,et al.  Measures of agreement between computation and experiment: Validation metrics , 2004, J. Comput. Phys..

[37]  Steven Deutsch,et al.  Wall shear-rate estimation within the 50cc Penn State artificial heart using particle image velocimetry. , 2004, Journal of biomechanical engineering.

[38]  Steven Deutsch,et al.  Multilaboratory particle image velocimetry analysis of the FDA benchmark nozzle model to support validation of computational fluid dynamics simulations. , 2011, Journal of biomechanical engineering.