Computational Analysis of Virtual Echocardiographic Assessment of Functional Mitral Regurgitation for Validation of Proximal Isovelocity Surface Area Methods

Patient-specific analysis of bicuspid aortic valve hemodynamics using a fully coupled fluid-structure interaction model

Bicuspid aortic valve (BAV), the most common congenital heart disease, is prone to develop significant valvular dysfunction and aortic wall abnormalities such as ascending aortic aneurysm. Growing evidence has suggested that abnormal BAV hemodynamics could contribute to disease progression. In order to investigate BAV hemodynamics, we performed 3D patient-specific fluid-structure interaction (FSI) simulations with fully coupled blood flow dynamics and valve motion throughout the cardiac cycle. Results showed that the hemodynamics during systole can be characterized by a systolic jet and two counter-rotating recirculation vortices. At peak systole, the jet was usually eccentric, with asymmetric recirculation vortices and helical flow motion in the ascending aorta. The flow structure at peak systole was quantified using the vorticity, flow rate reversal ratio and local normalized helicity (LNH) at four locations from the aortic root to the ascending aorta. The systolic jet was evaluated with the peak velocity, normalized flow displacement, and jet angle. It was found that peak velocity and normalized flow displacement (rather than jet angle) gave a strong correlation with the vorticity and LNH in the ascending aorta, which suggests that these two metrics could be used for clinical noninvasive evaluation of abnormal blood flow patterns in BAV patients.

Challenges and pitfalls in classification of disproportionate mitral regurgitation

The concept of disproportionate mitral regurgitation (dispropMR) has been introduced to identify patients with functional mitral regurgitation (MR) who benefit from percutaneous treatment. We aimed to examine echocardiographic characteristics behind this entity. We retrospectively included 172 consecutive patients with reduced left ventricular ejection fraction (LVEF), and more than mild MR referred to clinically indicated echocardiography. According to the proportionality ratio (effective regurgitant orifice area (EROA)/left ventricular end-diastolic volume (LVEDV)) patients were divided into dispropMR and proportionate MR (propMR) group. Potential factors which might affect proportionality definition were analyzed. 55 patients (32%) had dispropMR. Discrepant grading of MR severity was observed when using regurgitant volume (RegVol) by proximal isovelocity surface area (PISA) method or volumetric method, with significant discordance only in dispropMR (p < 0.001). Patients with dispropMR had more frequently left ventricular foreshortened images for LVEDV calculation than patients with propMR (p = 0.003), resulting in smaller LVEDV in dispropMR group. DispropMR group had more substantial dynamic variation of regurgitant flow compared to propMR. Accordingly, EROA was consistently overestimated by standard single-point PISA method compared to serial PISA method. This was more pronounced in dispropMR (bias:10.5 ± 28.3 mm2) compared to propMR group (bias:6.4 ± 12.8 mm2). DispropMR may be found in roughly one third of clinically indicated echocardiographic studies in patients with reduced LVEF and more than mild MR. EROA overestimation due to dynamic variation of regurgitant flow and LVEDV underestimation due to LV foreshortening were more frequently found in dispropMR. Our results indicate that methodological limitations of echocardiographic MR grading could not be neglected in classifying the proportionality of MR.

Effect of Nonlinear Hyperelastic Property of Arterial Tissues on the Pulse Wave Velocity Based on the Unified-Fiber-Distribution (UFD) Model

Pulse wave velocity (PWV) is a key, independent risk factor for future cardiovascular events. The Moens-Korteweg equation describes the relation between PWV and the stiffness of arterial tissue with an assumption of isotopic linear elastic property of the arterial wall. However, the arterial tissue exhibits highly nonlinear and anisotropic mechanical behaviors. There is a limited study regarding the effect of arterial nonlinear and anisotropic properties on the PWV. In this study, we investigated the impact of the arterial nonlinear hyperelastic properties on the PWV, based on our recently developed unified-fiber-distribution (UFD) model. The UFD model considers the fibers (embedded in the matrix of the tissue) as a unified distribution, which expects to be more physically consistent with the real fiber distribution than existing models that separate the fiber distribution into two/several fiber families. With the UFD model, we fitted the measured relation between the PWV and blood pressure which obtained a good accuracy. We also modeled the aging effect on the PWV based on observations that the stiffening of arterial tissue increases with aging, and the results agree well with experimental data. In addition, we did parameter studies on the dependence of the PWV on the arterial properties of fiber initial stiffness, fiber distribution, and matrix stiffness. The results indicate the PWV increases with increasing overall fiber component in the circumferential direction. The dependences of the PWV on the fiber initial stiffness, and matrix stiffness are not monotonic and change with different blood pressure. The results of this study could provide new insights into arterial property changes and disease information from the clinical measured PWV data.

Impact of assuming a circular orifice on flow error through elliptical regurgitant orifices: computational fluid dynamics and in vitro analysis of proximal flow convergence

Grounded in hydrodynamic theory, proximal isovelocity surface area (PISA) is a simplistic and practical technique widely used to quantify valvular regurgitation flow. PISA provides a relatively reasonable, though slightly underestimated flow rate for circular orifices. However, for elliptical orifices frequently seen in functional mitral regurgitation, PISA underestimates the flow rate. Based on data obtained with computational fluid dynamics (CFD) and in vitro experiments using systematically varied orifice parameters, we hypothesized that flow rate underestimation for elliptical orifices by PISA is predictable and within a clinically acceptable range. We performed 45 CFD simulations with varying orifice areas 0.1, 0.3 and 0.5 cm², orifice aspect ratios 1:1, 2:1, 3:1, 5:1, and 10:1, and peak velocities (V_max) 400, 500 and 600 cm/s. The ratio of computed effective regurgitant orifice area to true effective area (EROA_C/EROA) against the ratio of aliasing velocity to peak velocity (V_A/V_max) was analyzed for orifice shape impact. Validation was conducted with in vitro imaging in round and 3:1 elliptical orifices. Plotting EROA_C/EROA against V_A/V_max revealed marginal flow underestimation with 2:1 and 3:1 elliptical axis ratios against a circular orifice (< 10% for 8% V_A/V_max), rising to ≤ 35% for 10:1 ratio. In vitro modeling confirmed CFD findings; there was a 8.3% elliptical EROA underestimation compared to the circular orifice estimate. PISA quantification for regurgitant flow through elliptical orifices produces predictable, but generally small, underestimation deemed clinically acceptable for most regurgitant orifices.

Visualization and Quantification of the Unrepaired Complete Atrioventricular Canal Valve Using Open-Source Software

BACKGROUND

Repair of complete atrioventricular canal (CAVC) is often complicated by residual left atrioventricular valve regurgitation. The structure of the mitral and tricuspid valves in biventricular hearts has previously been shown to be associated with valve dysfunction. However, the three-dimensional (3D) structure of the entire unrepaired CAVC valve has not been quantified. Understanding the 3D structure of the CAVC may inform optimized repair.

METHODS

Novel open-source work flows were created in SlicerHeart for the modeling and quantification of CAVC valves on the basis of 3D echocardiographic images. These methods were applied to model the annulus, leaflets, and papillary muscle (PM) structure of 35 patients (29 with trisomy 21) with CAVC using transthoracic 3D echocardiography. The mean leaflet and annular shapes were calculated and visualized using shape analysis. Metrics of the complete native CAVC valve structure were compared with those of normal mitral valves using the Mann-Whitney U test. Associations between CAVC structure and atrioventricular valve regurgitation were analyzed.

RESULTS

CAVC leaflet metrics varied throughout systole. Compared with normal mitral valves, the left CAVC PMs were more acutely angled in relation to the annular plane (P < .001). In addition, the anterolateral PM was laterally and inferiorly rotated in CAVC, while the posteromedial PM was more superiorly and laterally rotated, relative to normal mitral valves (P < .001). Lower native CAVC atrioventricular valve annular height and annular height-to-valve width ratio before repair were both associated with moderate or greater left atrioventricular valve regurgitation after repair (P < .01).

CONCLUSIONS

It is feasible to model and quantify 3D CAVC structure using 3D echocardiographic images. The results demonstrate significant variation in CAVC structure across the cohort and differences in annular, leaflet, and PM structure compared with the mitral valve. These tools may be used in future studies to catalyze future research intended to identify structural associations of valve dysfunction and to optimize repair in this vulnerable and complex population.

Engineering analysis of aortic wall stress and root dilatation in the V-shape surgery for treatment of ascending aortic aneurysms

OBJECTIVES

The study objective was to evaluate the aortic wall stress and root dilatation before and after the novel V-shape surgery for the treatment of ascending aortic aneurysms and root ectasia.

METHODS

Clinical cardiac computed tomography images were obtained for 14 patients [median age, 65 years (range, 33-78); 10 (71%) males] who underwent the V-shape surgery. For 10 of the 14 patients, the computed tomography images of the whole aorta pre- and post-surgery were available, and finite element simulations were performed to obtain the stress distributions of the aortic wall at pre- and post-surgery states. For 6 of the 14 patients, the computed tomography images of the aortic root were available at 2 follow-up time points post-surgery (Post 1, within 4 months after surgery and Post 2, about 20-52 months from Post 1). We analysed the root dilatation post-surgery using change of the effective diameter of the root at the two time points and investigated the relationship between root wall stress and root dilatation.

RESULTS

The mean and peak max-principal stresses of the aortic root exhibit a significant reduction, P=0.002 between pre- and post-surgery for both root mean stress (median among the 10 patients presurgery, 285.46 kPa; post-surgery, 199.46 kPa) and root peak stress (median presurgery, 466.66 kPa; post-surgery, 342.40 kPa). The mean and peak max-principal stresses of the ascending aorta also decrease significantly from pre- to post-surgery, with P=0.004 for the mean value (median presurgery, 296.48 kPa; post-surgery, 183.87 kPa), and P=0.002 for the peak value (median presurgery, 449.73 kPa; post-surgery, 282.89 kPa), respectively. The aortic root diameter after the surgery has an average dilatation of 5.01% in total and 2.15%/year. Larger root stress results in larger root dilatation.

CONCLUSIONS

This study marks the first biomechanical analysis of the novel V-shape surgery. The study has demonstrated significant reduction in wall stress of the aortic root repaired by the surgery. The root was able to dilate mildly post-surgery. Wall stress could be a critical factor for the dilatation since larger root stress results in larger root dilatation. The dilated aortic root within 4 years after surgery is still much smaller than that of presurgery.

A Novel Approach for Semi Automated 3D Quantification of Mitral Regurgitant Volume Reflects a More Physiologic Approach to MR

BACKGROUND

Quantification of mitral regurgitation (MR) by echocardiography is an integral to assessing lesion severity, and entails integration of multiple Doppler-based parameters. These methods are primarily founded upon the principle of PISA (proximal isovelocity surface area), a 2D method known to employ several assumptions regarding MR jet characteristics. We analyzed the results of a semi-automated method of 3D-based RV estimation which accounts for jet behavior throughout the cardiac cycle, and compared it to conventional 2D PISA methods for MR.

METHODS

A total of 50 patients referred for transesophageal echocardiogram (TEE) for evaluation of primary (n= 25) and secondary MR (n=25) were included for analysis. 3D full volume color data sets were acquired, along with standard 2D methods for PISA calculation. 3D semi-automated MR flow quantification algorithm was applied offline to calculate 3D regurgitant volume (RVol), with simultaneous temporal curves generated from the 3D dataset. 3DRvol was compared to 2DRVol. 3D vena contracta area was also performed in all cases.

RESULTS

There was a modest correlation between 2DRVol and 3DRVol (r = 0.60). The semi-automated 3D approach resulted in significantly lower RV values compared to 2D PISA. Real-time and dynamic flow curve patterns were used for integral estimates of 3DRVol over the cardiac cycle, with a distinct bimodal pattern in functional MR, and brief and solitary peak in primary.

CONCLUSIONS

Using a semi-automated 3D software for quantification of mitral regurgitation allows for simultaneous calculation of 3D RVol with an automated generation of dynamic flow curves characteristic of the underlying MR mechanism. Our flow curve pattern results highlight well-known differences between MR flow dynamics in degenerative MR compared to functional MR.