The Effect of Vortex Formation on Left Ventricular Filling and Mitral Valve Efficiency

An Anatomically Shaped Mitral Valve for Hemodynamic Testing

In vitro modeling of the left heart relies on accurately replicating the physiological conditions of the native heart. The targeted physiological conditions include the complex fluid dynamics coming along with the opening and closing of the aortic and mitral valves. As the mitral valve possess a highly sophisticated apparatus, thence, accurately modeling it remained a missing piece in the perfect heart duplicator puzzle. In this study, we explore using a hydrogel-based mitral valve that offers a full representation of the mitral valve apparatus. The valve is tested using a custom-made mock circulatory loop to replicate the left heart. The flow analysis includes performing particle image velocimetry measurements in both left atrium and ventricle. The results showed the ability of the new mitral valve to replicate the real interventricular and atrial flow patterns during the whole cardiac cycle. Moreover, the investigated valve has a ventricular vortex formation time of 5.2, while the peak e- and a-wave ventricular velocities was 0.9 m/s and 0.4 m/s respectively.

Effect of transcatheter edge-to-edge repair device position on diastolic hemodynamic parameters: An echocardiography-based simulation study

Background

Transcatheter edge-to-edge repair (TEER) has developed from innovative technology to an established treatment strategy of mitral regurgitation (MR). The risk of iatrogenic mitral stenosis after TEER is, however, a critical factor in the conflict of interest between maximal reduction of MR and minimal impairment of left ventricular filling. We aim to investigate systematically the impact of device position on the post treatment hemodynamic outcome by involving the patient-specific segmentation of the diseased mitral valve.

Materials and methods

Transesophageal echocardiographic image data of ten patients with severe MR (age: 57 ± 8 years, 20% female) were segmented and virtually treated with TEER at three positions by using a position based dynamics approach. Pre- and post-interventional patient geometries were preprocessed for computational fluid dynamics (CFD) and simulated at peak-diastole with patient-specific blood flow boundary conditions. Simulations were performed with boundary conditions mimicking rest and stress. The simulation results were compared with clinical data acquired for a cohort of 21 symptomatic MR patients (age: 79 ± 6 years, 43% female) treated with TEER.

Results

Virtual TEER reduces the mitral valve area (MVA) from 7.5 ± 1.6 to 2.6 ± 0.6 cm2. Central device positioning resulted in a 14% smaller MVA than eccentric device positions. Furthermore, residual MVA is better predictable for central than for eccentric device positions (R2 = 0.81 vs. R2 = 0.49). The MVA reduction led to significantly higher maximal diastolic velocities (pre: 0.9 ± 0.2 m/s, post: 2.0 ± 0.5 m/s) and pressure gradients (pre: 1.5 ± 0.6 mmHg, post: 16.3 ± 9 mmHg) in spite of a mean flow rate reduction by 23% due to reduced MR after the treatment. On average, velocities were 12% and pressure gradients were 25% higher with devices in central compared to lateral or medial positions.

Conclusion

Virtual TEER treatment combined with CFD is a promising tool for predicting individual morphometric and hemodynamic outcomes. Such a tool can potentially be used to support clinical decision making, procedure planning, and risk estimation to prevent post-procedural iatrogenic mitral stenosis.

CT-Based Simulation of Left Ventricular Hemodynamics: A Pilot Study in Mitral Regurgitation and Left Ventricle Aneurysm Patients

Background

Cardiac CT (CCT) is well suited for a detailed analysis of heart structures due to its high spatial resolution, but in contrast to MRI and echocardiography, CCT does not allow an assessment of intracardiac flow. Computational fluid dynamics (CFD) can complement this shortcoming. It enables the computation of hemodynamics at a high spatio-temporal resolution based on medical images. The aim of this proposed study is to establish a CCT-based CFD methodology for the analysis of left ventricle (LV) hemodynamics and to assess the usability of the computational framework for clinical practice.

Materials and Methods

The methodology is demonstrated by means of four cases selected from a cohort of 125 multiphase CCT examinations of heart failure patients. These cases represent subcohorts of patients with and without LV aneurysm and with severe and no mitral regurgitation (MR). All selected LVs are dilated and characterized by a reduced ejection fraction (EF). End-diastolic and end-systolic image data was used to reconstruct LV geometries with 2D valves as well as the ventricular movement. The intraventricular hemodynamics were computed with a prescribed-motion CFD approach and evaluated in terms of large-scale flow patterns, energetic behavior, and intraventricular washout.

Results

In the MR patients, a disrupted E-wave jet, a fragmentary diastolic vortex formation and an increased specific energy dissipation in systole are observed. In all cases, regions with an impaired washout are visible. The results furthermore indicate that considering several cycles might provide a more detailed view of the washout process. The pre-processing times and computational expenses are in reach of clinical feasibility.

Conclusion

The proposed CCT-based CFD method allows to compute patient-specific intraventricular hemodynamics and thus complements the informative value of CCT. The method can be applied to any CCT data of common quality and represents a fair balance between model accuracy and overall expenses. With further model enhancements, the computational framework has the potential to be embedded in clinical routine workflows, to support clinical decision making and treatment planning.

Alterations in Intracardiac Flow Patterns Affect Mitral Leaflets Dynamics in a Model of Ischemic Mitral Regurgitation

PURPOSE

This study was to evaluate the effects of ischemic mitral regurgitation (IMR) on vortex formation and leaflet dynamics using an established porcine infarct model of IMR.

METHODS

Using direct coronary ligation, five animals were subjected to a posterolateral myocardial infarction (MI) followed by an MRI at 12-weeks post MI. MR imaging consisted of 4D time-resolved left ventricular (LV) flow, full coverage 2D LV cine, and high resolution 2D cine of mitral valve dynamics. Five additional naïve animals underwent identical imaging protocols to serve as controls. Image analysis was performed to obtain mitral transvalvular flows as well as LV volumes throughout the cardiac cycle. In addition, anterior to posterior mid-leaflet tip distances were measured throughout the cardiac cycle for determination of temporal leaflet dynamics.

RESULTS

It was found IMR caused asymmetric vortex ring formation with the anterior vortex having a lower vorticity relative to its posterior counterpart. In contrast, normal ventricles create symmetric and tightly curled vortices in the basal chamber just underneath the mitral leaflets which conserve kinetic energy and aid in effective ejection. IMR animals were also evaluated for leaflet separation and were found to have a greater leaflet opening and achieved peak vorticity and peak leaflet opening later than control animals.

CONCLUSION

In conclusion, this study shows the effects that altered vortex formation, due to IMR, can have on ventricular filling and leaflet dynamics. These findings have important implications for understanding blood flow through the dilated heart and how ring annuloplasty and volume reduction interventions may influence mitral valve dynamics.

Impact of prosthetic mitral valve orientation on the ventricular flow field: Comparison using patient-specific computational fluid dynamics

Significant mitral valve regurgitation creates progressive adverse remodeling of the left ventricle (LV). Replacement of the failing valve with a prosthesis generally improves patient outcomes but leaves the patient with non-physiological intracardiac flow patterns that might contribute to their future risk of thrombus formation and embolism. It has been suggested that the angular orientation of the implanted valve might modify the postoperative distortion of the intraventricular flow field. In this study, we investigated the effect of prosthetic valve orientation on LV flow patterns by using heart geometry from a patient with LV dysfunction and a competent native mitral valve to calculate intracardiac flow fields with computational fluid dynamics (CFD). Results were validated using in vivo 4D Flow MRI. The computed flow fields were compared to calculations following virtual implantation of a mechanical heart valve oriented in four different angles to assess the effect of leaflet position. Flow patterns were visualized in long- and short-axes and quantified with flow component analysis. In comparison to a native valve, valve implantation increased the proportion of the mitral inflow remaining in the basal region and further increased the residual volume in the apical area. Only slight changes due to valve orientation were observed. Using our numerical framework, we demonstrated quantitative changes in left ventricular blood flow due to prosthetic mitral replacement. This framework may be used to improve design of prosthetic heart valves and implantation procedures to minimize the potential for apical flow stasis, and potentially assist personalized treatment planning.

Diastolic Vortex Alterations With Reducing Left Ventricular Volume: An In Vitro Study

Despite the large number of studies of intraventricular filling dynamics for potential clinical applications, little is known as to how the diastolic vortex ring properties are altered with reduction in internal volume of the cardiac left ventricle (LV). The latter is of particular importance in LV diastolic dysfunction (LVDD) and in congenital diseases such as hypertrophic cardiomyopathy (HCM), where LV hypertrophy (LVH) can reduce LV internal volume. We hypothesized that peak circulation and the rate of decay of circulation of the diastolic vortex would be altered with reducing end diastolic volume (EDV) due to increasing confinement. We tested this hypothesis on physical models of normal LV and HCM geometries, under identical prescribed inflow profiles and for multiple EDVs, using time-resolved particle image velocimetry (TR-PIV) measurements on a left heart simulator. Formation and pinch-off of the vortex ring were nearly unaffected with changes to geometry and EDV. Pinch-off occurred before the end of early filling (E-wave) in all test conditions. Peak circulation of the vortex core near the LV outflow tract (LVOT) increased with lowering EDV and was lowest for the HCM model. The rate of decay of normalized circulation in dimensionless formation time (T*) increased with decreasing EDV. When using a modified version of T* that included average LV cross-sectional area and EDV, normalized circulation of all tested EDVs collapsed closely in the normal LV model (10% maximum difference between EDVs). Collectively, our results show that LV shape and internal volume play a critical role in diastolic vortex ring dynamics.

Sex Differences in Cardiac Flow Dynamics of Healthy Volunteers

Purpose

The purpose of this study was to further understand the relationship between cardiac function and flow, on the basis of sex, by quantifying cardiac flow characteristics and relating them to cardiac muscle performance in young adults.

Materials and Methods

In this cross-sectional study, cardiac four-dimensional flow (4D flow) magnetic resonance imaging (MRI) and two-dimensional cine MRI were performed on 20 male and 19 female volunteers aged 20-35. Velocity-based metrics of flow, kinetic energy, vorticity, and efficiency indices were quantified, as well as cardiac strain metrics.

Results*

Peak systolic blood kinetic energy (male: 4.76 ± 2.66 mJ; female: 3.36 ± 1.43 mJ; p=0.047) was significantly higher in the male left ventricle (LV) than in the female LV. Peak systolic vorticity index (male: 0.008 ± 0.005 rad-m²/ml-s; female: 0.014 ± 0.007 rad-m²/ml-s; p=0.007), peak diastolic vorticity index (male: 0.007 ± 0.006 rad-m²/ml-s; female: 0.014 ± 0.010 rad-m²/ml-s; p=0.015), and cycle-average vorticity (male: 0.006 ± 0.001 rad-m²/ml-s; female: 0.011 ± 0.002 rad/s; p=0.001) were all significantly higher in the LV of women than they were in the LV of men. Radial, circumferential, and long-axis strain metrics were significantly higher in the female LV than in the male LV (p<0.05). Circumferential systolic and diastolic strain rates displayed moderate correlation to peak systolic (r=-0.38, p=0.022) and diastolic vorticity (r=0.40, p=0.015) values, respectively. *Results are reported as mean ± standard deviation.

Conclusion

Left ventricular vorticity metrics were observed to be higher in women than in men and displayed moderate correlation to cardiac strain metrics. The methods and results of this study may be used to further understand the sex-based cardiac efficiency relationship between cardiac function and flow.

Thrust and Hydrodynamic Efficiency of the Bundled Flagella

The motility mechanism of prokaryotic organisms has inspired many untethered microswimmers that could potentially perform minimally invasive medical procedures in stagnant fluid regions inside the human body. Some of these microswimmers are inspired by bacteria with single or multiple helical flagella to propel efficiently and fast. For multiple flagella configurations, the direct measurement of thrust and hydrodynamic propulsion efficiency has been challenging due to the ambiguous mechanical coupling between the flow field and mechanical power input. To address this challenge and to compare alternative micropropulsion designs, a methodology based on volumetric velocity field acquisition is developed to acquire the key propulsive performance parameters from scaled-up swimmer prototypes. A digital particle image velocimetry (PIV) analysis protocol was implemented and experiments were conducted with the aid of computational fluid dynamics (CFD). First, this methodology was validated using a rotating single-flagellum similitude model. In addition to the standard PIV error assessment, validation studies included 2D vs. 3D PIV, axial vs. lateral PIV and simultaneously acquired direct thrust force measurement comparisons. Compatible with typical micropropulsion flow regimes, experiments were conducted both for very low and higher Reynolds (Re) number regimes (up to a Re number = 0.01) than that are reported in the literature. Finally, multiple flagella bundling configurations at 0°, 90° and 180° helical phase-shift angles were studied using scaled-up multiple concentric flagella thrust elements. Thrust generation was found to be maximal for the in-phase (0°) bundling configuration but with ~50% lower hydrodynamic efficiency than the single flagellum. The proposed measurement protocol and static thrust test-bench can be used for bio-inspired microscale propulsion methods, where direct thrust and efficiency measurement are required.

Analysis of cavopulmonary and cardiac flow characteristics in fontan Patients: Comparison with healthy volunteers

BACKGROUND

Characterizing the flow of the Fontan circuit, and correlating flow characteristics with the development of complications, is an important clinical challenge. Past work has analyzed the flow characteristics of Fontan circulation on a component-by-component basis. 4D flow MRI with radial projections allows for large volumetric coverage, and therefore can be used to analyze the flow through many codependent cardiovascular components in a single imaging session.

PURPOSE

To describe flow characteristics across the entire Fontan circuit and to compare these with the flow characteristics in healthy volunteers.

STUDY TYPE

Prospective.

SUBJECTS

Eleven single ventricle patients with a Fontan connection and 15 healthy controls.

SEQUENCE

Phase contrast with vastly undersampled isotropic projection reconstruction (PC-VIPR) at a field strength of 3 T.

ASSESSMENT

Cavopulmonary and ventricular flow distributions, blood flow kinetic energy, vorticities, efficiency indices, and other flow parameters were analyzed using Ensight and MatLab.

STATISTICAL TESTS

The results were compared across Fontan subjects, between respiratory phases, and between Fontan subjects and healthy volunteers using a Student's t-test for unequal sample sizes and linear regression.

RESULTS

Cava-specific pulmonary flow distributions of Fontan patients varied significantly between respiratory phases (P < 0.05). Ventricular kinetic energy (KE) was significantly higher in Fontan patients than it was in healthy controls, leading to a lower cardiac efficiency metric in the Fontan group. A significant diastolic KE time-shift was also observed in the Fontan patient group. Peak diastolic KE was significantly higher in the single ventricle of patients with right ventricle morphology than it was in left ventricle morphology patients.

DATA CONCLUSION

Radial 4D flow MRI can be used for comprehensive analysis of single ventricle Fontan flow characteristics.

LEVEL OF EVIDENCE

2 Technical Efficacy: Stage 2 J. Magn. Reson. Imaging 2019.

Computational medical imaging and hemodynamics framework for functional analysis and assessment of cardiovascular structures

Cardiac dysfunction constitutes common cardiovascular health issues in the society, and has been an investigation topic of strong focus by researchers in the medical imaging community. Diagnostic modalities based on echocardiography, magnetic resonance imaging, chest radiography and computed tomography are common techniques that provide cardiovascular structural information to diagnose heart defects. However, functional information of cardiovascular flow, which can in fact be used to support the diagnosis of many cardiovascular diseases with a myriad of hemodynamics performance indicators, remains unexplored to its full potential. Some of these indicators constitute important cardiac functional parameters affecting the cardiovascular abnormalities. With the advancement of computer technology that facilitates high speed computational fluid dynamics, the realization of a support diagnostic platform of hemodynamics quantification and analysis can be achieved. This article reviews the state-of-the-art medical imaging and high fidelity multi-physics computational analyses that together enable reconstruction of cardiovascular structures and hemodynamic flow patterns within them, such as of the left ventricle (LV) and carotid bifurcations. The combined medical imaging and hemodynamic analysis enables us to study the mechanisms of cardiovascular disease-causing dysfunctions, such as how (1) cardiomyopathy causes left ventricular remodeling and loss of contractility leading to heart failure, and (2) modeling of LV construction and simulation of intra-LV hemodynamics can enable us to determine the optimum procedure of surgical ventriculation to restore its contractility and health This combined medical imaging and hemodynamics framework can potentially extend medical knowledge of cardiovascular defects and associated hemodynamic behavior and their surgical restoration, by means of an integrated medical image diagnostics and hemodynamic performance analysis framework.

Aortic Regurgitation Generates a Kinematic Obstruction Which Hinders Left Ventricular Filling

An incompetent aortic valve (AV) results in aortic regurgitation (AR), where retrograde flow of blood into the left ventricle (LV) is observed. In this work, we parametrically characterized the detailed changes in intra-ventricular flow during diastole as a result of AR in a physiological in vitro left-heart simulator (LHS). The loss of energy within the LV as the level of AR increased was also assessed. The validated LHS consisted of an optically-clear, flexible wall LV and a modular AV holder. Two-component, planar, digital particle image velocimetry was used to visualize and quantify intra-ventricular flow. A large coherent vortical structure which engulfed the whole LV was observed under control conditions. In the cases with AR, the regurgitant jet was observed to generate a "kinematic obstruction" between the mitral valve and the LV apex, preventing the trans-mitral jet from generating a coherent vortical structure. The regurgitant jet was also observed to impinge on the inferolateral wall of the LV. Energy dissipation rate (EDR) for no, trace, mild, and moderate AR were found to be 1.15, 2.26, 3.56, and 5.99 W/m³, respectively. This study has, for the first time, performed an in vitro characterization of intra-ventricular flow in the presence of AR. Mechanistically, the formation of a "kinematic obstruction" appears to be the cause of the increased EDR (a metric quantifiable in vivo) during AR. EDR increases non-linearly with AR fraction and could potentially be used as a metric to grade severity of AR and develop clinical interventional timing strategies for patients.

Role of Mitral Annulus Diastolic Geometry on Intraventricular Filling Dynamics

The mitral valve (MV) is a bileaflet valve positioned between the left atrium and ventricle of the heart. The annulus of the MV has been observed to undergo geometric changes during the cardiac cycle, transforming from a saddle D-shape during systole to a flat (and less eccentric) D-shape during diastole. Prosthetic MV devices, including heart valves and annuloplasty rings, are designed based on these two configurations, with the circular design of some prosthetic heart valves (PHVs) being an approximation of the less eccentric, flat D-shape. Characterizing the effects of these geometrical variations on the filling efficiency of the left ventricle (LV) is required to understand why the flat D-shaped annulus is observed in the native MV during diastole in addition to optimizing the design of prosthetic devices. We hypothesize that the D-shaped annulus reduces energy loss during ventricular filling. An experimental left heart simulator (LHS) consisting of a flexible-walled LV physical model was used to characterize the filling efficiency of the two mitral annular geometries. The strength of the dominant vortical structure formed and the energy dissipation rate (EDR) of the measured fields, during the diastolic period of the cardiac cycle, were used as metrics to quantify the filling efficiency. Our results indicated that the O-shaped annulus generates a stronger (25% relative to the D-shaped annulus) vortical structure than that of the D-shaped annulus. It was also found that the O-shaped annulus resulted in higher EDR values throughout the diastolic period of the cardiac cycle. The results support the hypothesis that a D-shaped mitral annulus reduces dissipative energy losses in ventricular filling during diastole and in turn suggests that a symmetric stent design does not provide lower filling efficiency than an equivalent asymmetric design.

A D-Shaped Bileaflet Bioprosthesis which Replicates Physiological Left Ventricular Flow Patterns

Prior studies have shown that in a healthy heart, there exist a large asymmetric vortex structure that aids in establishing a steady flow field in the left ventricle. However, the implantation of existing artificial heart valves at the mitral position is found to have a negative effect on this physiological flow pattern. In light of this, a novel D-shaped bileaflet porcine bioprosthesis (GD valve) has been designed based on the native geometry mitral valve, with the hypothesis that biomimicry in valve design can restore physiological left ventricle flow patterns after valve implantation. An in-vitro experiment using two dimensional particle velocimetry imaging was carried out to determine the hemodynamic performance of the new bileaflet design and then compared to that of the well-established St. Jude Epic valve which functioned as a control in the experiment. Although both valves were found to have similar Reynolds shear stress and Turbulent Kinetic Energy levels, the novel D-shape valve was found to have lower turbulence intensity and greater mean kinetic energy conservation.

Atrial systole enhances intraventricular filling flow propagation during increasing heart rate

Diastolic fluid dynamics in the left ventricle (LV) has been examined in multiple clinical studies for understanding cardiac function in healthy humans and developing diagnostic measures in disease conditions. The question of how intraventricular filling vortex flow pattern is affected by increasing heart rate (HR) is still unanswered. Previous studies on healthy subjects have shown a correlation between increasing HR and diminished E/A ratio of transmitral peak velocities during early filling (E-wave) to atrial systole (A-wave). We hypothesize that with increasing HR under constant E/A ratio, E-wave contribution to intraventricular vortex propagation is diminished. A physiologic in vitro flow phantom consisting of a LV physical model was used for this study. HR was varied across 70, 100 and 120 beats per minute (bpm) with E/A of 1.1-1.2. Intraventricular flow patterns were characterized using 2D particle image velocimetry measured across three parallel longitudinal (apical-basal) planes in the LV. A pair of counter-rotating vortices was observed during E-wave across all HRs. With increasing HR, diminished vortex propagation occurred during E-wave and atrial systole was found to amplify secondary vorticity production. The diastolic time point where peak vortex circulation occurred was delayed with increasing HR, with peak circulation for 120bpm occurring as late as 90% into diastole near the end of A-wave. The role of atrial systole is elevated for higher HR due to the limited time available for filling. Our baseline findings and analysis approach can be applied to studies of clinical conditions where impaired exercise tolerance is observed.

Cardiovascular magnetic resonance compatible physical model of the left ventricle for multi-modality characterization of wall motion and hemodynamics

BACKGROUND

The development of clinically applicable fluid-structure interaction (FSI) models of the left heart is inherently challenging when using in vivo cardiovascular magnetic resonance (CMR) data for validation, due to the lack of a well-controlled system where detailed measurements of the ventricular wall motion and flow field are available a priori. The purpose of this study was to (a) develop a clinically relevant, CMR-compatible left heart physical model; and (b) compare the left ventricular (LV) volume reconstructions and hemodynamic data obtained using CMR to laboratory-based experimental modalities.

METHODS

The LV was constructed from optically clear flexible silicone rubber. The geometry was based off a healthy patient's LV geometry during peak systole. The LV phantom was attached to a left heart simulator consisting of an aorta, atrium, and systemic resistance and compliance elements. Experiments were conducted for heart rate of 70 bpm. Wall motion measurements were obtained using high speed stereo-photogrammetry (SP) and cine-CMR, while flow field measurements were obtained using digital particle image velocimetry (DPIV) and phase-contrast magnetic resonance (PC-CMR).

RESULTS

The model reproduced physiologically accurate hemodynamics (aortic pressure = 120/80 mmHg; cardiac output = 3.5 L/min). DPIV and PC-CMR results of the center plane flow within the ventricle matched, both qualitatively and quantitatively, with flow from the atrium into the LV having a velocity of about 1.15 m/s for both modalities. The normalized LV volume through the cardiac cycle computed from CMR data matched closely to that from SP. The mean difference between CMR and SP was 5.5 ± 3.7%.

CONCLUSIONS

The model presented here can thus be used for the purposes of: (a) acquiring CMR data for validation of FSI simulations, (b) determining accuracy of cine-CMR reconstruction methods, and

The Influence of Mitral Annuloplasty on Left Ventricular Flow Dynamics

BACKGROUND

Mitral valve (MV) repair using annuloplasty rings is the preferred method of treatment for MV regurgitation, but the impact of annuloplasty ring placement on left ventricular intraventricular flow has not been studied.

METHODS

Annuloplasty rings of varying sizes were placed in 5 healthy sheep (intercommissural ring sizes were 24, 26, 28, 30, and 32 mm), and three-dimensional phase contrast magnetic resonance imaging (4D flow MRI) was performed before and 1 week after ring placement.

RESULTS

Normal diastolic flow consisted of diastolic intraventricular vortices that naturally unwound during systole. Postsurgical intraventricular flow was highly disturbed in all sheep, and the disturbance was greatest for undersized rings. Ring size was highly correlated with the diastolic inflow angle (Pearson's r = -0.62, p < 0.1, 95% confidence interval: -0.92 to 0.14). There was a mean angle increase of mean diastolic inflow angle increase of 12.3 degrees (< 30 mm, p < 0.01, 95% confidence interval: 4.8 to 19.6) for rings less than 30 mm. There was an inverse relationship between peak velocity and annuloplasty ring area (Pearson's r = -0.80, p < 0.05, 95% confidence interval: -0.96 to -0.2). Transmitral pressure gradients increased significantly from baseline 0.73 ± 0.18 mm Hg to after annuloplasty 2.31 ± 1.04 mm Hg (p < 0.05).

CONCLUSIONS

Mitral valve annuloplasty ring placement disturbs normal left ventricular intraventricular flow patterns, and the degree of disturbance is closely associated with annuloplasty ring size.

Numerical Modeling of Intraventricular Flow during Diastole after Implantation of BMHV

This work presents a numerical simulation of intraventricular flow after the implantation of a bileaflet mechanical heart valve at the mitral position. The left ventricle was simplified conceptually as a truncated prolate spheroid and its motion was prescribed based on that of a healthy subject. The rigid leaflet rotation was driven by the transmitral flow and hence the leaflet dynamics were solved using fluid-structure interaction approach. The simulation results showed that the bileaflet mechanical heart valve at the mitral position behaved similarly to that at the aortic position. Sudden area expansion near the aortic root initiated a clockwise anterior vortex, and the continuous injection of flow through the orifice resulted in further growth of the anterior vortex during diastole, which dominated the intraventricular flow. This flow feature is beneficial to preserving the flow momentum and redirecting the blood flow towards the aortic valve. To the best of our knowledge, this is the first attempt to numerically model intraventricular flow with the mechanical heart valve incorporated at the mitral position using a fluid-structure interaction approach. This study facilitates future patient-specific studies.

Towards a fast and efficient approach for modelling the patient-specific ventricular haemodynamics

Computer modelling of the heart has emerged over the past decade as a powerful technique to explore the cardiovascular pathophysiology and inform clinical diagnosis. The current state-of-the-art in biophysical modelling requires a wealth of, potentially invasive, clinical data for the parametrisation and validation of the models, a process that is still too long and complex to be compatible with the clinical decision-making time. Therefore, there remains a need for models that can be quickly customised to reconstruct physical processes difficult to measure directly in patients. In this paper, we propose a less resource-intensive approach to modelling, whereby computational fluid-dynamics (CFD) models are constrained exclusively by boundary motion derived from imaging data through a validated wall tracking algorithm. These models are generated and parametrised based solely on ultrasound data, whose acquisition is fast, inexpensive and routine in all patients. To maximise the time and computational efficiency, a semi-automated pipeline is embedded in an image processing workflow to personalise the models. Applying this approach to two patient cases, we demonstrate this tool can be directly used in the clinic to interpret and complement the available clinical data by providing a quantitative indication of clinical markers that cannot be easily derived from imaging, such as pressure gradients and the flow energy.

Computational Study of Hemodynamic Effects of Abnormal E/A Ratio on Left Ventricular Filling

Three-dimensional numerical simulations are employed to investigate the hemodynamic effects of abnormal E/A ratios on left ventricular filling. The simulations are performed in a simplified geometric model of the left ventricle (LV) in conjunction with a specified endocardial motion. The model has been carefully designed to match the important geometric and flow parameters under the physiological conditions. A wide range of E/A ratios from 0 to infinity is employed with the aim to cover all the possible stages of left ventricle diastolic dysfunction (DD). The effects of abnormal E/A ratios on vortex dynamics, flow propagation velocity, energy consumption as well as flow transport and mixing are extensively discussed. Our results are able to confirm some common findings reported by the previous studies, and also uncover some interesting new features.

Numerical simulation of patient-specific left ventricular model with both mitral and aortic valves by FSI approach

Intraventricular flow is important in understanding left ventricular function; however, relevant numerical simulations are limited, especially when heart valve function is taken into account. In this study, intraventricular flow in a patient-specific left ventricle has been modelled in two-dimension (2D) with both mitral and aortic valves integrated. The arbitrary Lagrangian-Eulerian (ALE) approach was employed to handle the large mesh deformation induced by the beating ventricular wall and moving leaflets. Ventricular wall deformation was predefined based on MRI data, while leaflet dynamics were predicted numerically by fluid-structure interaction (FSI). Comparisons of simulation results with in vitro and in vivo measurements reported in the literature demonstrated that numerical method in combination with MRI was able to predict qualitatively the patient-specific intraventricular flow. To the best of our knowledge, we are the first to simulate patient-specific ventricular flow taking into account both mitral and aortic valves.

Intraventricular flow patterns and stasis in the LVAD-assisted heart

Left ventricular assist device (LVAD) support disrupts the natural blood flow path through the heart, introducing flow patterns associated with thrombosis, especially in the presence of medical devices. The aim of this study was to quantitatively evaluate the flow patterns in the left ventricle (LV) of the LVAD-assisted heart, with a focus on alterations in vortex development and stasis. Particle image velocimetry of a LVAD-supported LV model was performed in a mock circulatory loop. In the Pre-LVAD flow condition, a vortex ring initiating from the LV base migrated toward the apex during diastole and remained in the LV by the end of ejection. During LVAD support, vortex formation was relatively unchanged although vortex circulation and kinetic energy increased with LVAD speed, particularly in systole. However, as pulsatility decreased and aortic valve opening ceased, a region of fluid stasis formed near the left ventricular outflow tract. These findings suggest that LVAD support does not substantially alter vortex dynamics unless cardiac function is minimal. The altered blood flow introduced by the LVAD results in stasis adjacent to the LV outflow tract, which increases the risk of thrombus formation in the heart.

Simulation of left ventricle flow dynamics with dilated cardiomyopathy during the filling phase

Dilated cardiomyopathy (DCM) is a common cardiac disease which leads to the deterioration in cardiac performance. A computational fluid dynamics (CFD) approach can be used to enhance our understanding of the disease, by providing us with a detailed map of the intraventricular flow and pressure distributions. In the present work, effect of ventricular size on the intraventricular flow dynamics and intraventricular pressure gradients (IVPGs) was studied using two different implementation methods, i.e. the geometry-prescribed and the fluid structure interaction (FSI) methods. Results showed that vortex strength and IVPGs are significantly reduced in a dilated heart, leading to an increased risk of thrombus formation and impaired ventricular filling. We suggest FSI method as the ultimate method in studying ventricular dysfunction as it provides additional cardiac disease prognostic factors and more realistic model implementation.

Fluid-structure interaction of an aortic heart valve prosthesis driven by an animated anatomic left ventricle

We develop a novel large-scale kinematic model for animating the left ventricle (LV) wall and use this model to drive the fluid-structure interaction (FSI) between the ensuing blood flow and a mechanical heart valve prosthesis implanted in the aortic position of an anatomic LV/aorta configuration. The kinematic model is of lumped type and employs a cell-based, FitzHugh-Nagumo framework to simulate the motion of the LV wall in response to an excitation wavefront propagating along the heart wall. The emerging large-scale LV wall motion exhibits complex contractile mechanisms that include contraction (twist) and expansion (untwist). The kinematic model is shown to yield global LV motion parameters that are well within the physiologic range throughout the cardiac cycle. The FSI between the leaflets of the mechanical heart valve and the blood flow driven by the dynamic LV wall motion and mitral inflow is simulated using the curvilinear immersed boundary (CURVIB) method [1, 2] implemented in conjunction with a domain decomposition approach. The computed results show that the simulated flow patterns are in good qualitative agreement with in vivo observations. The simulations also reveal complex kinematics of the valve leaflets, thus, underscoring the need for patient-specific simulations of heart valve prosthesis and other cardiac devices.

FSI simulation of intra-ventricular flow in patient-specific ventricular model with both mitral and aortic valves

Investigating the intra-ventricular flow is the most important to understand the left ventricular function. In this study, we proposed a fluid-structure interaction (FSI) approach to simulate the blood flow in patient-specific model by combining both mitral and aortic valves. To accommodate the large mesh deformation, moving arbitrary Lagrangian-Eulerian (ALE) meshes were used for moving ventricular wall and rotating leaflets of valves. The left ventricular wall was predescribed according to the points acquired from magnetic resonance image (MRI). Mitral and aortic valves were integrated into the model by assuming each leaflet as a rigid body. Fluid-structure interaction (FSI) approach was adopted to capture the rapid motion of leaflets. The simulation results were qualitatively similar to the measurements reported in literatures. To the best of our knowledge, this is the first to simulate the patient-specific ventricular flow with the presence of both mitral and aortic valves.

MEDICAL IMAGING AND COMPUTER-AIDED FLOW ANALYSIS OF A HEART WITH ATRIAL SEPTAL DEFECT

Computer-aided magnetic resonance (MR) fluid motion tracking and cardiac vorticity quantification of the right atrial flow is implemented in this study to suggest a new method for the diagnosis of an atrial septal defect (ASD). MR signals of blood moving in a cardiac chamber can be represented as an image and vary in intensity at every consecutive cardiac phase. A method was devised to perform flow analysis using MR imaging without modification of scan mode or protocol that allows velocity encoding. A single vortex or multiple vortices may appear in the cardiac chamber. However, velocity fields in any flow scenario are normally unable to reveal information for a concise analysis; therefore, in addition to velocity maps, vorticity flow maps on which the velocity field is superimposed are presented. Through a case study, the difference in vortex strengths pre- and post-atrial septal occlusion can be examined, and the results can be verified using computational fluid dynamics. Based on this framework, the degree of vortical flow was assessed for the right atrium of a subject with atrial septal defect. A relationship can be established between right atrial vorticity and the ASD. As such, there is clear utility of the developed system in its potential as a prognostic and investigative tool for the quantitative assessment of cardiac abnormalities parallel to examining magnetic resonance images.

Assessment of left ventricular hemodynamics and function of patients with uremia by vortex formation using vector flow mapping

A novel echocardiographic method, vector flow mapping (VFM), acquires velocity vector from color Doppler velocity data. The purpose of this study was to evaluate whether VFM could provide useful information on intracardiac flow and helpful to evaluate left ventricular (LV) function. Thirty-eight patients with uremia undergoing hemodialysis and 30 healthy volunteers were enrolled. The maximum vector velocity, maximum diameter and duration of the intracardiac vortex were measured using VFM software during systole and diastole. The maximum vector velocity of the vortex and the peak velocities at the basal septum and lateral mitral annulus measured by tissue Doppler imaging (TDI) were correlated. The maximum diameter and duration of vortex formation were significantly higher in uremic patients compared with the control group during the ejection phase (40.6 ± 7.9 cm/sec vs. 28.1 ± 3.9 cm/sec; 297.1 ± 22.1 msec vs. 145.4 ± 19.3 msec, all P < 0.001). The maximal diameters of the vortex were higher in uremic patients compared with the control group during diastole (25.6 ± 3.4 mm vs. 16.4 ± 2.1 mm; 34.3 ± 3.1 mm vs. 26.8 ± 3.9 mm; 37.5 ± 2.4 mm vs. 20.9 ± 2.1 mm; all P < 0.001). The maximum vector velocities were lower in mid-diastole and late diastole (23.6 ± 2.3 cm/sec vs. 45.2 ± 3.7 cm/sec; 31.9 ± 2.9 cm/sec vs. 54.7 ± 3.2 cm/sec, all P < 0.001). There was a correlation between the maximum vector velocity of the vortex in mid-diastole and E'/A' at the septum and lateral mitral annulus (r = 0.70, r = 0.76, P < 0.001). Vortex can be utilized to provide intracardiac dynamic information using VFM and it may be a good supplement for evaluating LV function.

Medical image diagnostics based on computer-aided flow analysis using magnetic resonance images

Most of the cardiac abnormalities have an implication on hemodynamics and affect cardiovascular health. Diagnostic imaging modalities such as computed tomography and magnetic resonance imaging provide excellent anatomical information on myocardial structures, but fail to show the cardiac flow and detect heart defects in vivo condition. The computerized technique for fluid motion estimation by pixel intensity tracking based on magnetic resonance signals represents a promising technique for functional assessment of cardiovascular disease, as it can provide functional information of the heart in addition to analysis of its anatomy. Cardiovascular flow characteristics can be measured in both normal controls and patients with cardiac abnormalities such as atrial septal defect, thus, enabling identification of the underlying causes of these flow phenomena. This review paper focuses on an overview of a flow analysis scheme based on computer-aided evaluation of magnetic resonance intensity images, in comparison with other commonly used medical imaging modalities. Details of the proposed technique are provided with validations being conducted at selected abnormal cardiovascular patients. It is expected that this new technique can potentially extend applications for characterizing cardiovascular defects and their hemodynamic behavior.

Towards patient-specific cardiovascular modeling system using the immersed boundary technique

BACKGROUND

Previous research shows that the flow dynamics in the left ventricle (LV) reveal important information about cardiac health. This information can be used in early diagnosis of patients with potential heart problems. The current study introduces a patient-specific cardiovascular-modelling system (CMS) which simulates the flow dynamics in the LV to facilitate physicians in early diagnosis of patients before heart failure.

METHODS

The proposed system will identify possible disease conditions and facilitates early diagnosis through hybrid computational fluid dynamics (CFD) simulation and time-resolved magnetic resonance imaging (4-D MRI). The simulation is based on the 3-D heart model, which can simultaneously compute fluid and elastic boundary motions using the immersed boundary method. At this preliminary stage, the 4-D MRI is used to provide an appropriate comparison. This allows flexible investigation of the flow features in the ventricles and their responses.

RESULTS

The results simulate various flow rates and kinetic energy in the diastole and systole phases, demonstrating the feasibility of capturing some of the important characteristics of the heart during different phases. However, some discrepancies exist in the pulmonary vein and aorta flow rate between the numerical and experimental data. Further studies are essential to investigate and solve the remaining problems before using the data in clinical diagnostics.

CONCLUSIONS

The results show that by using a simple reservoir pressure boundary condition (RPBC), we are able to capture some essential variations found in the clinical data. Our approach establishes a first-step framework of a practical patient-specific CMS, which comprises a 3-D CFD model (without involving actual hemodynamic data yet) to simulate the heart and the 4-D PC-MRI system. At this stage, the 4-D PC-MRI system is used for verification purpose rather than input. This brings us closer to our goal of developing a practical patient-specific CMS, which will be pursued next. We anticipate that in the future, this hybrid system can potentially identify possible disease conditions in LV through comprehensive analysis and facilitates physicians in early diagnosis of probable cardiac problems.

Left Ventricular Vortex Under Mitral Valve Edge-to-Edge Repair

Mitral valve (MV) edge-to-edge repair (ETER) changes MV geometry by approximation of MV leaflets, and impacts left ventricle (LV) filling fluid mechanics. The purpose of this study was to investigate LV vortex with MV ETER during diastole. A computational MV-LV model was developed with MV ETER at the central free edges of the anterior and posterior leaflets. It was supposed that LV would elongate apically during diastole. The elongation deformation was controlled by the intraventricular flow rate. MV leaflets were modeled as a semi-prolate sphere with two symmetrical circular orifices and fixed at the maximum valve opening. MV chordae were neglected. FLUENT was used to simulate blood flow through the MV and in the LV. MV ETER generated two jets deflected laterally toward the LV wall in rapid LV filling. The jets impinged the LV wall obliquely and moved apically along the LV wall. Jet energy was primarily lost near the impingement. The jet from each MV orifice was surrounded by a vortex ring. The two vortex rings dissipated at the end of diastole. The total energy loss increased inversely with the MV orifice area. The atrio-ventricular pressure gradient was adverse near the end of diastole and possibly in diastasis. Reduction of the total orifice area led to more increment in the transmitral pressure drop than in the transmitral velocity. In conclusion, during diastole, two deflected jets from the MV under ETER impinged the LV wall. Major energy loss occurred around the jet impingement. Two vortex rings dissipated at the end of diastole with little storage of inflow energy for blood ejection in the following process of systole. MV ETER increased energy loss and lowered LV filling efficiency. The maintaining of a larger orifice area after ETER might not significantly increase energy loss in the LV during diastole and the transmitral pressure drop. The adverse pressure gradient from the atrium to the LV might be the mechanism of MV closure in the late diastole.

Two-dimensional intraventricular flow mapping by digital processing conventional color-Doppler echocardiography images

Doppler echocardiography remains the most extended clinical modality for the evaluation of left ventricular (LV) function. Current Doppler ultrasound methods, however, are limited to the representation of a single flow velocity component. We thus developed a novel technique to construct 2D time-resolved (2D+t) LV velocity fields from conventional transthoracic clinical acquisitions. Combining color-Doppler velocities with LV wall positions, the cross-beam blood velocities were calculated using the continuity equation under a planar flow assumption. To validate the algorithm, 2D Doppler flow mapping and laser particle image velocimetry (PIV) measurements were carried out in an atrio-ventricular duplicator. Phase-contrast magnetic resonance (MR) acquisitions were used to measure in vivo the error due to the 2D flow assumption and to potential scan-plane misalignment. Finally, the applicability of the Doppler technique was tested in the clinical setting. In vitro experiments demonstrated that the new method yields an accurate quantitative description of the main vortex that forms during the cardiac cycle (mean error for vortex radius, position and circulation). MR image analysis evidenced that the error due to the planar flow assumption is close to 15% and does not preclude the characterization of major vortex properties neither in the normal nor in the dilated LV. These results are yet to be confirmed by a head-to-head clinical validation study. Clinical Doppler studies showed that the method is readily applicable and that a single large anterograde vortex develops in the healthy ventricle while supplementary retrograde swirling structures may appear in the diseased heart. The proposed echocardiographic method based on the continuity equation is fast, clinically-compliant and does not require complex training. This technique will potentially enable investigators to study of additional quantitative aspects of intraventricular flow dynamics in the clinical setting by high-throughput processing conventional color-Doppler images.

Cardiac flow component analysis

In a chamber of the heart, large-scale vortices are shown to exist as the result of the dynamic blood flow and unique morphological changes of the chamber wall. As the cardiovascular flow varies over a cardiac cycle, there is a need for a robust quantification method to analyze its vorticity and circulation. We attempt to measure vortex characteristics by means of two-dimensional vorticity maps and vortex circulation. First, we develop vortex component analysis by segmenting the vortices using an data clustering algorithm before histograms of their vorticity distribution are generated. The next stage is to generate the statistics of the vorticity maps for each phase of the cardiac cycle to allow analysis of the flow. This is followed by evaluating the circulation of each segmented vortex. The proposed approach is dedicated to examining vortices within the human heart chamber. The vorticity field can indicate the strength and number of large-scale vortices in the chamber. We provide the results of the flow analysis after vorticity map segmentation and the statistical properties that characterize the vorticity components. The success of the cardiac measurement and analysis is illustrated by a case study of the right atrium. Our investigation shows that it is possible to utilize a data clustering algorithm to segment vortices after vorticity mapping, and that the vorticity and circulation analysis of a chamber vorticity can provide new insights into the blood flow within the cardiovascular structure.

Recent advances in the application of computational mechanics to the diagnosis and treatment of cardiovascular disease

During the last 30 years, research into the pathogenesis and progression of cardiovascular disease has had to employ a multidisciplinary approach involving a wide range of subject areas, from molecular and cell biology to computational mechanics and experimental solid and fluid mechanics. In general, research was driven by the need to provide answers to questions of critical importance for disease management. Ongoing improvements in the spatial resolution of medical imaging equipment coupled to an exponential growth in the capacity, flexibility and speed of computational techniques have provided a valuable opportunity for numerical simulations and complex experimental techniques to make a contribution to improving the diagnosis and clinical management of many forms of cardiovascular disease. This paper contains a review of recent progress in the numerical simulation of cardiovascular mechanics, focusing on three particular areas: patient-specific modeling and the optimization of surgery in pediatric cardiology, evaluating the risk of rupture in aortic aneurysms, and noninvasive characterization of intraventricular flow in the management of heart failure.