Analysis of the effect of flow rate on the Doppler continuity equation for stenotic orifice area calculations: a numerical study.

BACKGROUND Flow-rate dependencies of the Doppler continuity equation are addressed in this study. METHODS AND RESULTS By use of computational fluid dynamic (CFD) software with turbulence modeling, three-dimensional axisymmetric models of round stenotic orifices were created. Flow simulations were run for various orifice area sizes (0.785, 1.13, 1.76, and 3.14 cm2) and flow rates (0.37 to 25.0 L/min). Reynolds numbers ranged from 100 to 8000. Once adequate convergence was obtained with each simulation, the location of the vena contracta was determined. For each run, maximum and average velocities across the cross section of the vena contracta were tabulated and vena contracta cross-sectional area (effective orifice area) determined. The difference between the maximum velocity and the average velocity at the vena contracta was smallest at high-flow states, with more of a difference at low-flow states. At lower-flow states, the velocity vector profile at the vena contracta was parabolic, whereas at high-flow states, the profile became more flattened. Also, the effective orifice area (vena contracta cross-sectional area) varied with flow rate. At moderate-flow states, the effective orifice area reached a minimum and expanded at low- and high-flow states, remaining relatively constant at high-flow states. CONCLUSIONS We have shown that significant differences exist between the maximum velocity and the average velocity at the vena contracta at low flow rates. A likely explanation for this is that viscous effects cause lower velocities at the edges of the vena contracta at low flow rates, resulting in a parabolic profile. At higher-flow states, inertial forces overcome viscous drag, causing a flatter profile. Effective orifice area itself varies with flow rate as well, with the smallest areas seen at moderate-flow states. These flow-dependent factors lead to flow rate-dependent errors in the Doppler continuity equation. Our results have strong relevance to clinical measurements of stenotic valve areas by use of the Doppler continuity equation under varying cardiac output conditions.

[1]  L Hatle,et al.  DIAGNOSTIC METHODS , 2005 .

[2]  Michael Jones,et al.  Hysteresis behavior of stenotic aortic bioprosthesis as a function of flow rate: In vitrostudies , 1996 .

[3]  Weeratunge Malalasekera,et al.  An introduction to computational fluid dynamics - the finite volume method , 2007 .

[4]  H. Q. Yang,et al.  An experimental-computational analysis of MHV cavitation: effects of leaflet squeezing and rebound. , 1994, The Journal of heart valve disease.

[5]  M Jones,et al.  In vitro hemodynamic characteristics of tissue bioprostheses in the aortic position. , 1986, The Journal of thoracic and cardiovascular surgery.

[6]  H Reul,et al.  Comparison of valvular resistance, stroke work loss, and Gorlin valve area for quantification of aortic stenosis. An in vitro study in a pulsatile aortic flow model. , 1995, Circulation.

[7]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[8]  O. Thual Turbulence and Random Processes in Fluid Mechanics , 1988 .

[9]  R. Shandas,et al.  Insights into catheter/Doppler discrepancies in congenital aortic stenosis. , 1999, The American journal of cardiology.

[10]  E D Verrier,et al.  Dependence of Gorlin Formula and Continuit Equation Valve Areas on Transvalvular Volume Flow Rate in Valvular Aortic Stenosis , 1994, Circulation.

[11]  D. Wilcox Turbulence modeling for CFD , 1993 .

[12]  W. Janna,et al.  Introduction to Fluid Mechanics , 2012 .

[13]  J. P. V. Doormaal,et al.  ENHANCEMENTS OF THE SIMPLE METHOD FOR PREDICTING INCOMPRESSIBLE FLUID FLOWS , 1984 .

[14]  A. J. Przekwas,et al.  General Strong Conservation Formulation of Navier-Stokes Equations in Nonorthogonal Curvilinear Coordinates , 1994 .

[15]  A. Yoganathan,et al.  Theoretical and practical differences between the Gorlin formula and the continuity equation for calculating aortic and mitral valve areas. , 1991, The American journal of cardiology.

[16]  J. Dumesnil,et al.  Validation and applications of mitral prosthetic valvular areas calculated by Doppler echocardiography. , 1990, The American journal of cardiology.

[17]  M. Crawford,et al.  Hydraulic estimation of stenotic orifice area: a correction of the Gorlin formula. , 1985, Circulation.

[18]  Physiologic Changes with Maximal Exercise in Asymptomatic Valvular Aortic Stenosis Assessed by Doppler Echocardiography , 1993 .

[19]  R. S. Mitchell,et al.  When should Doppler-determined valve area be better than the Gorlin formula?: Variation in hydraulic constants in low flow states. , 1987, Journal of the American College of Cardiology.

[20]  V. Makhijani,et al.  Application of computational fluid dynamics techniques to blood pumps. , 1996, Artificial organs.

[21]  Y. Fung,et al.  Mechanics of the Circulation , 2011, Developments in Cardiovascular Medicine.

[22]  C. Otto,et al.  Experimental Validation of Doppler Echocardiographic Measurement of Volume Flow Through the Stenotic Aortic Valve , 1988, Circulation.

[23]  N H Hwang,et al.  Numerical study of squeeze-flow in tilting disc mechanical heart valves. , 1996, The Journal of heart valve disease.

[24]  J. M. Siegel,et al.  Computational simulation of magnetic resonance angiograms in stenotic vessels: Effect of stenosis severity , 1996 .

[25]  R Gorlin,et al.  Calculations of cardiac valve stenosis: restoring an old concept for advanced applications. , 1987, Journal of the American College of Cardiology.

[26]  Hwang Nh,et al.  Numerical study of squeeze-flow in tilting disc mechanical heart valves. , 1996 .

[27]  A. R. Snider,et al.  Echocardiography in Pediatric Heart Disease , 1990 .

[28]  S. Goldberg,et al.  Noninvasive quantification of stenotic semilunar valve areas by Doppler echocardiography. , 1984, Journal of the American College of Cardiology.

[29]  M. Vanauker,et al.  Turbulent/viscous interactions control Doppler/catheter pressure discrepancies in aortic stenosis. The role of the Reynolds number. , 1996, Circulation.

[30]  A. D. Young,et al.  An Introduction to Fluid Mechanics , 1968 .

[31]  P J Dionne,et al.  Three‐Dimensional Coupled Fluid‐Structure Simulation of Pericardial Bioprosthetic Aortic Valve Function , 1997, ASAIO journal.

[32]  I. Palacios,et al.  Effects of dobutamine on Gorlin and continuity equation valve areas and valve resistance in valvular aortic stenosis. , 1992, The American journal of cardiology.

[33]  R GORLIN,et al.  Hydraulic formula for calculation of the area of the stenotic mitral valve, other cardiac valves, and central circulatory shunts. I. , 1951, American heart journal.