On Mitral Valve Dynamics and its Connection to Early Diastolic Flow

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.

A Hybrid Approach for Cardiac Blood Flow Vortex Ring Identification Based on Optical Flow and Lagrangian Averaged Vorticity Deviation

Objective: The measurement of cardiac blood flow vortex characteristics can help to facilitate the analysis of blood flow dynamics that regulates heart function. However, the complexity of cardiac flow along with other physical limitations makes it difficult to adequately identify the dominant vortices in a heart chamber, which play a significant role in regulating the heart function. Although the existing vortex quantification methods can achieve this goal, there are still some shortcomings: such as low precision, and ignoring the center of the vortex without the description of vortex deformation processes. To address these problems, an optical flow Lagrangian averaged vorticity deviation (Optical flow-LAVD) method is proposed.

Methodology: We examined the flow within the right atrium (RA) of the participants’ hearts, by using a single set of scans pertaining to a slice at two-chamber short-axis orientation. Toward adequate extraction of the vortex ring characteristics, a novel approach driven by the Lagrangian averaged vorticity deviation (LAVD) was implemented and applied to characterize the trajectory integral associated with vorticity deviation and the spatial mean of rings, by using phase-contrast magnetic resonance imaging (PC-MRI) datasets as a case study. To interpolate the time frames between every larger discrete frame and minimize the error caused by constructing a continuous velocity field for the integral process of LAVD, we implemented the optical flow as an interpolator and introduced the backward warping as an intermediate frame synthesis basis, which is then used to generate higher quality continuous velocity fields.

Results: Our analytical study results showed that the proposed Optical flow-LAVD method can accurately identify vortex ring and continuous velocity fields, based on optical flow information, for yielding high reconstruction outcomes. Compared with the linear interpolation and phased-based frame interpolation methods, our proposed algorithm can generate more accurate synthesized PC-MRI.

Conclusion: This study has developed a novel Optical flow-LAVD model to accurately identify cardiac vortex rings, and minimize the associated errors caused by the construction of a continuous velocity field. Our paper presents a superior vortex characteristics detection method that may potentially aid the understanding of medical experts on the dynamics of blood flow within the heart.

Deep Learning Algorithm to Detect Cardiac Sarcoidosis From Echocardiographic Movies

BACKGROUND

Because the early diagnosis of subclinical cardiac sarcoidosis (CS) remains difficult, we developed a deep learning algorithm to distinguish CS patients from healthy subjects using echocardiographic movies.Methods and Results:Among the patients who underwent echocardiography from January 2015 to December 2019, we chose 151 echocardiographic movies from 50 CS patients and 151 from 149 healthy subjects. We trained two 3D convolutional neural networks (3D-CNN) to identify CS patients using a dataset of 212 echocardiographic movies with and without a transfer learning method (Pretrained algorithm and Non-pretrained algorithm). On an independent set of 41 echocardiographic movies, the area under the receiver-operating characteristic curve (AUC) of the Pretrained algorithm was greater than that of Non-pretrained algorithm (0.842, 95% confidence interval (CI): 0.722-0.962 vs. 0.724, 95% CI: 0.566-0.882, P=0.253). The AUC from the interpretation of the same set of 41 echocardiographic movies by 5 cardiologists was not significantly different from that of the Pretrained algorithm (0.855, 95% CI: 0.735-0.975 vs. 0.842, 95% CI: 0.722-0.962, P=0.885). A sensitivity map demonstrated that the Pretrained algorithm focused on the area of the mitral valve.

CONCLUSIONS

A 3D-CNN with a transfer learning method may be a promising tool for detecting CS using an echocardiographic movie.

Numerical simulations of flow patterns in the human left ventricle model with a novel dynamic mesh morphing approach based on radial basis function

We present a new numerical simulation framework for prediction of flow patterns in the human left ventricle model. In this study, a radial basis function (RBF) mesh morphing method is developed and applied within the finite-volume computational fluid dynamics (CFD) approach. The numerical simulations are designed to closely mimic details of recent tomographic particle image velocimetry (TomoPIV) experiments. The numerically simulated dynamic motions of the left ventricle and tri-leaflet biological mitral valve are emulated through the RBF morphing method. The arbitrary Lagrangian-Eulerian (ALE) based CFD is performed with the RBF-defined deforming wall boundaries. The results obtained show a good agreement with experiments, confirming the reliability and accuracy of the developed simulation framework.

Intraventricular Vortex Interaction between Transmitral Flow and Paravalvular Leak

Paravalvular leak (PVL) is a complication of transcatheter aortic valve replacement. Despite its marked clinical impact, no previous study has reported how PVL affects the intraventricular fluid dynamics. This study aims to delineate vortex interaction between PVL and transmitral flow and the influence of PVL orifice location on intraventricular fluid dynamics using Echocardiographic Particle Image Velocimetry. Three different conditions of no PVL, anterior PVL and posterior PVL were experimentally studied and clinically compared. Circulation, impulse, kinetic energy (KE) and change in KE (ΔKE) were calculated. As well, vortex formation analyses and streamline description were performed to study vortex interactions. The anterior PVL jet streamed into the LV and interfered with the transmitral flow. Posterior PVL jet formed a large clockwise vortex and collided with transmitral flow, which resulted in flow disturbance. Compared to no PVL condition, average circulation, impulse, KE and ΔKE increased in presence of PVL. In conclusion, we found that PVL jets lead to abnormal vortex formation that interfere with natural advancement of transmitral flow, and negatively affect the LV fluid dynamics parameters. PVL orifice location strongly affects the intraventricular vortex formation, and posterior PVL may have more negative effects compared to anterior PVL.

Synergy between Diastolic Mitral Valve Function and Left Ventricular Flow Aids in Valve Closure and Blood Transport during Systole

Highly resolved three-dimensional (3D) fluid structure interaction (FSI) simulation using patient-specific echocardiographic data can be a powerful tool for accurately and thoroughly elucidating the biomechanics of mitral valve (MV) function and left ventricular (LV) fluid dynamics. We developed and validated a strongly coupled FSI algorithm to fully characterize the LV flow field during diastolic MV opening under physiologic conditions. Our model revealed that distinct MV deformation and LV flow patterns developed during different diastolic stages. A vortex ring that strongly depended on MV deformation formed during early diastole. At peak E wave, the MV fully opened, with a local Reynolds number of ~5500, indicating that the flow was in the laminar-turbulent transitional regime. Our results  showed that during diastasis, the vortex structures caused the MV leaflets to converge, thus increasing mitral jet’s velocity. The vortex ring became asymmetrical, with the vortex structures on the anterior side being larger than on the posterior side. During the late diastolic stages, the flow structures advected toward the LV outflow tract, enhancing fluid transport to the aorta. This 3D-FSI study demonstrated the importance of leaflet dynamics, their effect on the vortex ring, and their influence on MV function and fluid transport within the LV during diastole.

Advanced Age Attenuates Left Ventricular Filling Efficiency Quantified Using Vortex Formation Time: A Study of Octogenarians With Normal Left Ventricular Systolic Function Undergoing Coronary Artery Surgery

OBJECTIVE

Blood flow across the mitral valve during early left ventricular (LV) filling produces a 3-dimensional rotational fluid body, known as a vortex ring, that enhances LV filling efficiency. Diastolic dysfunction is common in elderly patients, but the influence of advanced age on vortex formation is unknown. The authors tested the hypothesis that advanced age is associated with a reduction in LV filling efficiency quantified using vortex formation time (VFT) in octogenarians undergoing coronary artery bypass graft (CABG) surgery.

DESIGN

Observational study.

SETTING

Veterans Affairs medical center.

PARTICIPANTS

After institutional review board approval, octogenarians (n = 7; 82 ± 2 year [mean ± standard deviation]; ejection fraction 56% ± 7%) without valve disease or atrial arrhythmias undergoing CABG were compared with a younger cohort (n = 7; 55 ± 6 year; ejection fraction 57% ± 7%) who were undergoing coronary revascularization.

INTERVENTIONS

None.

MEASUREMENTS AND MAIN RESULTS

All patients were monitored using radial and pulmonary arterial catheters and transesophageal echocardiography. Peak early LV filling (E) and atrial systole (A) blood flow velocities and their corresponding velocity-time integrals were obtained using pulse-wave Doppler echocardiography to determine E/A, atrial filling fraction (β), and E wave deceleration time. Pulse-wave Doppler also was used to measure pulmonary venous blood flow during systole and diastole. Mitral valve diameter (D) was calculated as the average of major and minor axis lengths obtained in the midesophageal LV bicommissural and long-axis transesophageal echocardiography imaging planes, respectively. VFT was calculated as 4 × (1 - β) × SV/(πD³), where SV is the stroke volume measured using thermodilution. Systemic and pulmonary hemodynamics, LV diastolic function, and VFT were determined during steady-state conditions 30 minutes before cardiopulmonary bypass. A delayed relaxation pattern of LV filling (E/A 0.81 ± 0.16 v 1.29 ± 0.19, p = 0.00015; β 0.44 ± 0.05 v 0.35 ± 0.03, p = 0.0008; E wave deceleration time 294 ± 58 v 166 ± 28 ms, p < 0.0001; ratio of peak pulmonary venous systolic and diastolic blood flow velocity 1.42 ± 0.23 v 1.14 ± 0.20, p = 0.0255) was observed in octogenarians compared with younger patients. Mitral valve diameter was similar between groups (2.7 ± 0.2 and 2.6 ± 0.2 cm, respectively, in octogenarians v younger patients, p = 0.299). VFT was reduced in octogenarians compared with younger patients (3.0 ± 0.9 v 4.5 ± 1.2; p = 0.0171). An inverse correlation between age and VFT was shown using linear regression analysis (VFT = -0.0627 × age + 8.24; r² = 0.408; p = 0.0139).

CONCLUSION

The results indicate that LV filling efficiency quantified using VFT is reduced in octogenarians compared with younger patients undergoing coronary artery bypass grafting.

The relationship between systolic vector flow mapping parameters and left ventricular cardiac function in healthy dogs

Vector flow mapping (VFM) is a novel echocardiographic technology that shows blood flow vectors and vortexes, enabled the hydrokinetic evaluation of hemodynamics within the left ventricle. VFM provides several unique parameters: circulation, vorticity, vortex area, and energy loss. The present study aims to reveal a relationship between VFM parameters and cardiac function. Five healthy Beagle dogs were anesthetized and administered with dobutamine (0, 2, 4, 8, 12 µg/kg/min). Pressure-volume diagrams were acquired to assess cardiac function using pressure-volume conductance catheter. Systolic maximum circulation, vorticity, vortex area, and energy loss were measured using VFM. The systolic maximum circulation, systolic vorticity, systolic vortex area, and systolic energy loss were increased by dobutamine administration. There was a strongly significant correlation between the systolic maximum circulation and ejection fraction (r = 0.76), maximal positive left ventricular (LV) pressure derivatives (dP/dt max) (r = 0.80), and end-systolic LV elastance (r = 0.73). Systolic vorticity and systolic vortex area were strongly correlated with ejection fraction (r = 0.76, 0.68) and dP/dt max (r = 0.76, 0.69), and end-systolic LV elastance (r = 0.62, 0.74), respectively. Systolic energy loss was strongly correlated with dP/dt max (r = 0.78), systolic maximum circulation (r = 0.81), and systolic vorticity (r = 0.82). The present study revealed that systolic VFM parameters are associated with the LV contractility. Furthermore, systolic energy loss was susceptible to the systolic vortex parameters such as systolic vorticity and systolic maximum circulation. Systolic VFM parameters are new hydrokinetic indices reflecting LV contractility.

A calcified polymeric valve for valve-in-valve applications

The prevalence of aortic valve stenosis (AS) is increasing in the aging society. More recently, novel treatments and devices for AS, especially transcatheter aortic valve replacement (TAVR) have significantly changed the therapeutic approach to this disease. Research and development related to TAVR require testing these devices in the calcified heart valves that closely mimic a native calcific valve. However, no animal model of AS has yet been available. Alternatively, animals with normal aortic valve that are currently used for TAVR experiments do not closely replicate the aortic valve pathology required for proper testing of these devices. To solve this limitation, for the first time, we developed a novel polymeric valve whose leaflets possess calcium hydroxyapatite inclusions immersed in them. This study reports the characteristics and feasibility of these valves. Two types of the polymeric valve, i.e., moderate and severe calcified AS models were developed and tested by deploying a transcatheter valve in those and measuring the related hemodynamics. The valves were tested in a heart flow simulator, and were studied using echocardiography. Our results showed high echogenicity of the polymeric valve, that was correlated to the severity of the calcification. Aortic valve area of the polymeric valves was measured, and the severity of stenosis was defined according to the clinical guidelines. Accordingly, we showed that these novel polymeric valves closely mimic AS, and can be a desired cost-saving solution for testing the performance of the transcatheter aortic valve systems in vitro.

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.

Moderate Aortic Valvular Insufficiency Invalidates Vortex Formation Time as an Index of Left Ventricular Filling Efficiency in Patients With Severe Degenerative Calcific Aortic Stenosis Undergoing Aortic Valve Replacement

OBJECTIVE

Transmitral blood flow produces a vortex ring (quantified using vortex formation time [VFT]) that enhances the efficiency of left ventricular (LV) filling. VFT is attenuated in LV hypertrophy resulting from aortic valve stenosis (AS) versus normal LV geometry. Many patients with AS also have aortic insufficiency (AI). The authors tested the hypothesis that moderate AI falsely elevates VFT by partially inhibiting mitral leaflet opening in patients with AS.

DESIGN

Observational study.

SETTING

Veterans Affairs medical center.

PARTICIPANTS

Patients with AS in the presence or absence of moderate AI (n = 8 per group) undergoing aortic valve replacement (AVR) were studied after institutional review board approval.

INTERVENTIONS

None.

MEASUREMENTS AND MAIN RESULTS

Under general anesthesia, peak early LV filling (E) and atrial systole (A) blood flow velocities and their corresponding velocity-time integrals were obtained using pulse-wave Doppler transesophageal echocardiography (TEE) to determine E/A and atrial filling fraction (beta). Mitral valve diameter (D) was calculated as the average of major and minor axis lengths obtained in the midesophageal bicommissural (transcommissural anterior-lateral-posterior medial) and LV long-axis (anterior-posterior) TEE imaging planes, respectively. VFT was calculated as 4·(1-beta)·SV/πD(3), where SV = stroke volume measured using thermodilution. Hemodynamics, diastolic function, and VFT were determined during steady-state conditions before cardiopulmonary bypass. The severity of AS (mean and peak pressure gradients, peak transvalvular jet velocity, aortic valve area) and diastolic function (E/A, beta) were similar between groups. Moderate centrally directed AI was present in 8 patients with AS (ratio of regurgitant jet width to LV outflow tract diameter of 36±6%). Pulse pressure and mean pulmonary artery pressure were elevated in patients with versus without AI, but no other differences in hemodynamics were observed. Mitral valve minor and major axis lengths, diameter, and area were reduced in the presence versus the absence of AI. VFT was increased significantly (5.7±1.7 v 3.2±0.6; p = 0.00108) in patients with AS and AI compared with AS alone.

CONCLUSION

Moderate AI falsely elevates VFT in patients with severe AS undergoing AVR by partially inhibiting mitral valve opening. VFT may be an unreliable index of LV filling efficiency with competitive diastolic flow into the LV.

In Vitro 2D PIV Measurements and Related Aperture Areas of Tricuspid Bioprosthetic Mitral Valves at the beginning of Diastole

Purpose

Besides ventricular parameters, the design and angular orientation of a prosthetic heart valve induce a specific flow field. The aim of this study was to know the inflow characteristics of a left ventricular model (LVM), investigating the behavior of tricuspid bioprosthetic mitral valves in terms of velocity profiles and related valve aperture areas at the beginning of diastole, under different conditions.

Methods

3 heart rates (HRs) were established in the LVM and each mitral bioprosthesis (27 and 31 mm diameter) was installed in 2 orientations, rotated by 180°. For each experimental setup, 2-dimensional particle image velocimetry (2D PIV) measurements and simultaneous mitral valve (MV) area detection were obtained from 50 samples.

Results

The results from the velocity profiles immediately downstream of mitral bioprostheses showed the influence of valve orientation for moderate HRs, although for a similar magnitude of mean velocity vectors. The geometries of MV open areas for each HR were similar regardless of valve orientation, except for the 27-mm valve at 90 beats per minute (bpm), and for the 31-mm valve at 60 bpm. Moreover, for each HR, similar percentages of valve open area were obtained regardless of MV nominal diameters.

Conclusions

In conclusion, the experimental setup for the 2D PIV measurements synchronized with the MV area detection was a useful tool for knowing the inflow characteristics of the LVM.

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

A Reconstruction Method of Blood Flow Velocity in Left Ventricle Using Color Flow Ultrasound

Vortex flow imaging is a relatively new medical imaging method for the dynamic visualization of intracardiac blood flow, a potentially useful index of cardiac dysfunction. A reconstruction method is proposed here to quantify the distribution of blood flow velocity fields inside the left ventricle from color flow images compiled from ultrasound measurements. In this paper, a 2D incompressible Navier-Stokes equation with a mass source term is proposed to utilize the measurable color flow ultrasound data in a plane along with the moving boundary condition. The proposed model reflects out-of-plane blood flows on the imaging plane through the mass source term. The boundary conditions to solve the system of equations are derived from the dimensions of the ventricle extracted from 2D echocardiography data. The performance of the proposed method is evaluated numerically using synthetic flow data acquired from simulating left ventricle flows. The numerical simulations show the feasibility and potential usefulness of the proposed method of reconstructing the intracardiac flow fields. Of particular note is the finding that the mass source term in the proposed model improves the reconstruction performance.

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.

2-D left ventricular flow estimation by combining speckle tracking with Navier-Stokes-based regularization: an in silico, in vitro and in vivo study

Despite the availability of multiple ultrasound approaches to left ventricular (LV) flow characterization in two dimensions, this technique remains in its childhood and further developments seem warranted. This article describes a new methodology for tracking the 2-D LV flow field based on ultrasound data. Hereto, a standard speckle tracking algorithm was modified by using a dynamic kernel embedding Navier-Stokes-based regularization in an iterative manner. The performance of the proposed approach was first quantified in synthetic ultrasound data based on a computational fluid dynamics model of LV flow. Next, an experimental flow phantom setup mimicking the normal human heart was used for experimental validation by employing simultaneous optical particle image velocimetry as a standard reference technique. Finally, the applicability of the approach was tested in a clinical setting. On the basis of the simulated data, pointwise evaluation of the estimated velocity vectors correlated well (mean r = 0.84) with the computational fluid dynamics measurement. During the filling period of the left ventricle, the properties of the main vortex obtained from the proposed method were also measured, and their correlations with the reference measurement were also calculated (radius, r = 0.96; circulation, r = 0.85; weighted center, r = 0.81). In vitro results at 60 bpm during one cardiac cycle confirmed that the algorithm properly measures typical characteristics of the vortex (radius, r = 0.60; circulation, r = 0.81; weighted center, r = 0.92). Preliminary qualitative results on clinical data revealed physiologic flow fields.

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.

Current clinical application of intracardiac flow analysis using echocardiography

In evaluating the cardiac function, it is important to have a comprehensive assessment of structural factors, such as the myocardial or valvular function and intracardiac flow dynamics that pass the heart. Vortex flow that form during left ventricular filling have specific geometry and anatomical location that are critical determinants of directed blood flow during ejection. The formation of abnormal vortices relates to the abnormal cardiac function. Therefore, vortex flow may offer a novel index of cardiac dysfunction. Intracardiac flow visualization using ultrasound technique has definite advantages with a higher temporal resolution and availability in real time clinical setting. Vector flow mapping based on color-Doppler and contrast echocardiography using particle image velocimetry is currently being used for visualizing the intracardiac flow. The purpose of this review is to provide readers with an update on the current method for analyzing intracardiac flow using echocardiography and its clinical applications.

The vortex formation time to diastolic function relation: assessment of pseudonormalized versus normal filling

In early diastole, the suction pump feature of the left ventricle opens the mitral valve and aspirates atrial blood. The ventricle fills via a blunt profiled cylindrical jet of blood that forms an asymmetric toroidal vortex ring inside the ventricle whose growth has been quantified by the standard (dimensionless) expression for vortex formation time, VFT_standard = {transmitral velocity time integral}/{mitral orifice diameter}. It can differentiate between hearts having distinguishable early transmitral (Doppler E-wave) filling patterns. An alternative validated expression, VFT_kinematic reexpresses VFT_standard by incorporating left heart, near “constant-volume pump” physiology thereby revealing VFT_kinematic's explicit dependence on maximum rate of longitudinal chamber expansion (E′). In this work, we show that VFT_kinematic can differentiate between hearts having indistinguishable E-wave patterns, such as pseudonormal (PN; 0.75 < E/A < 1.5 and E/E′ > 8) versus normal. Thirteen age-matched normal and 12 PN data sets (738 total cardiac cycles), all having normal LVEF, were selected from our Cardiovascular Biophysics Laboratory database. Doppler E-, lateral annular E′-waves, and M-mode data (mitral leaflet separation, chamber dimension) was used to compute VFT_standard and VFT_kinematic. VFT_standard did not differentiate between groups (normal [3.58 ± 1.06] vs. PN [4.18 ± 0.79], P = 0.13). In comparison, VFT_kinematic for normal (3.15 ± 1.28) versus PN (4.75 ± 1.35) yielded P = 0.006. Hence, the applicability of VFT_kinematic for diastolic function quantitation has been broadened to include analysis of PN filling patterns in age-matched groups.

Intracardiac flow visualization: current status and future directions

Non-invasive cardiovascular imaging initially focused on heart structures, allowing the visualization of their motion and inferring its functional status from it. Colour-Doppler and cardiac magnetic resonance (CMR) have allowed a visual approach to intracardiac flow behaviour, as well as measuring its velocity at single selected spots. Recently, the application of new technologies to medical use and, particularly, to cardiology has allowed, through different algorithms in CMR and applications of ultrasound-related techniques, the description and analysis of flow behaviour in all points and directions of the selected region, creating the opportunity to incorporate new data reflecting cardiac performance to cardiovascular imaging. The following review provides an overview of the currently available imaging techniques that enable flow visualization, as well as its present and future applications based on the available literature and on-going works.

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.

Chronic pressure-overload hypertrophy attenuates vortex formation time in patients with severe aortic stenosis and preserved left ventricular systolic function undergoing aortic valve replacement

OBJECTIVE

Transmitral blood flow produces a vortex ring that enhances the hydraulic efficiency of early left ventricular (LV) filling. The effect of pressure-overload hypertrophy on the duration of LV vortex ring formation (vortex formation time [VFT]) is unknown. The current investigation tested the hypothesis that chronic LV pressure-overload hypertrophy produced by severe aortic stenosis (AS) reduces VFT in patients with preserved LV systolic function undergoing aortic valve replacement.

DESIGN

Observational study.

SETTING

Veterans Affairs Medical Center.

PARTICIPANTS

After the Institutional Review Board's approval, 8 patients (7 men and 1 woman; age, 62±5 y; and ejection fraction, 59%±5%) with AS (peak pressure gradient, 81±22 mmHg; aortic valve area, 0.78±0.25 cm(2)) scheduled for aortic valve replacement were compared with 8 patients (all men; age, 63±3 y; and ejection fraction, 60%±7%) without AS undergoing coronary artery bypass graft surgery.

INTERVENTIONS

None.

MEASUREMENTS AND MAIN RESULTS

Under general anesthesia, peak early LV filling (E) and atrial systole (A) blood flow velocities and their corresponding velocity-time integrals were obtained using pulse-wave Doppler echocardiography to determine E/A and atrial filling fraction (β). Mitral valve diameter (D) was calculated as the average of minor and major axis lengths obtained in the midesophageal bicommissural and long-axis transesophageal echocardiography imaging planes, respectively. Posterior wall thickness (PWT) was measured at end-diastole using M-mode echocardiography. VFT was calculated as 4×(1-β)×SV/πD(3), where SV = stroke volume measured using thermodilution. Systemic and pulmonary hemodynamics, LV diastolic function, PWT, and VFT were determined during steady-state conditions 30 minutes before cardiopulmonary bypass. Early LV filling was attenuated in patients with AS (eg, E/A, 0.77±0.11 compared with 1.23±0.13; β, 0.43±0.09 compared with 0.35±0.02; p<0.05 for each). LV hypertrophy was observed (PWT, 1.4±0.1 cm compared with 1.1±0.2 cm; p<0.05) and VFT was lower (3.0±0.9 v 4.3±0.5; p<0.05) in patients with versus without AS. Linear regression analysis showed a significant correlation between VFT and PWT (VFT = -2.57 ×PWT + 6.81; r(2) = 0.345; p = 0.017).

CONCLUSION

The results indicated that pressure-overload hypertrophy produced by AS reduced VFT in patients with normal LV systolic function undergoing aortic valve replacement.

The evolution of intraventricular vortex during ejection studied by using vector flow mapping

AIMS

The purpose of this study was to assess the evolution of intraventricular vortex during left ventricular (LV) ejection.

METHODS

Vector flow mapping was performed in 51 patients with coronary artery disease and LV ejection fraction (EF) >50%, 70 patients with EF <50% (13 with coronary artery disease and 57 with dilated cardiomyopathy), and 62 healthy volunteers.

RESULTS

In normals and patients with EF >50%, the intraventricular vortex dissipated quickly during early ejection. In patients with EF <50%, the vortex stayed mainly at apex and persisted for a significantly longer time. The evolution of vortex during ejection was significantly correlated with QRS width, EF, fractional shortening, LV outflow velocity time integral, wall motion score index (WMSI), LV dimensions, left atrial diameter, and diastolic mitral annular velocities. LV end-diastolic short diameter and WMSI were the independent determinants of the duration of vortex (R(2) = 0.482, P < 0.001). End-systolic short diameter and apical WMSI were the independent determinants of duration of vortex corrected for ejection time (R(2) = 0.565, P < 0.001). End-systolic short diameter was the independent determinant of percentage change in vortex area during early ejection (R(2) = 0.355, P < 0.001). End-systolic short diameter and ejection time were the independent determinants of percentage change in vortex flow volume (R(2) = 0.415, P < 0.001).

CONCLUSIONS

In patients with LV systolic dysfunction, the vortex persists during ejection and stays mainly at apex. The vortex evolution during ejection is closely associated with LV dimensions and functions.

Vortex ring formation in the left ventricle of the heart: analysis by 4D flow MRI and Lagrangian coherent structures

Recent studies suggest that vortex ring formation during left ventricular (LV) rapid filling is an optimized mechanism for blood transport, and that the volume of the vortex ring is an important measure. However, due to lack of quantitative methods, the volume of the vortex ring has not previously been studied. Lagrangian Coherent Structures (LCS) is a new flow analysis method, which enables in vivo quantification of vortex ring volume. Therefore, we aimed to investigate if vortex ring volume in the human LV can be reliably quantified using LCS and magnetic resonance velocity mapping (4D PC-MR). Flow velocities were measured using 4D PC-MR in 9 healthy volunteers and 4 patients with dilated ischemic cardiomyopathy. LV LCS were computed from flow velocities and manually delineated in all subjects. Vortex volume in the healthy volunteers was 51 ± 6% of the LV volume, and 21 ± 5% in the patients. Interobserver variability was -1 ± 13% and interstudy variability was -2 ± 12%. Compared to idealized flow experiments, the vortex rings showed additional complexity and asymmetry, related to endocardial trabeculation and papillary muscles. In conclusion, LCS and 4D PC-MR enables measurement of vortex ring volume during rapid filling of the LV.

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.

Mitral valve annuloplasty: a quantitative clinical and mechanical comparison of different annuloplasty devices

Mitral valve annuloplasty is a common surgical technique used in the repair of a leaking valve by implanting an annuloplasty device. To enhance repair durability, these devices are designed to increase leaflet coaptation, while preserving the native annular shape and motion; however, the precise impact of device implantation on annular deformation, strain, and curvature is unknown. In this article, we quantify how three frequently used devices significantly impair native annular dynamics. In controlled in vivo experiments, we surgically implanted 11 flexible-incomplete, 11 semi-rigid-complete, and 12 rigid-complete devices around the mitral annuli of 34 sheep, each tagged with 16 equally spaced tantalum markers. We recorded four-dimensional marker coordinates using biplane videofluoroscopy, first with device and then without, which were used to create mathematical models using piecewise cubic splines. Clinical metrics (characteristic anatomical distances) revealed significant global reduction in annular dynamics upon device implantation. Mechanical metrics (strain and curvature fields) explained this reduction via a local loss of anterior dilation and posterior contraction. Overall, all three devices unfavorably caused reduction in annular dynamics. The flexible-incomplete device, however, preserved native annular dynamics to a larger extent than the complete devices. Heterogeneous strain and curvature profiles suggest the need for heterogeneous support, which may spawn more rational design of annuloplasty devices using design concepts of functionally graded materials.

The left ventricular intracavitary vortex during the isovolumic contraction period as detected by vector flow mapping

AIMS

The purpose of this study was to characterize left ventricular (LV) intracavitary flow during the isovolumic contraction (IVC) period in humans using vector flow mapping.

METHODS

Color flow Doppler imaging was performed from the apical long-axis view in 61 patients with heart failure and 58 healthy volunteers. Doppler flow data obtained during IVC were analyzed offline with vector flow mapping.

RESULTS

A large vortex was formed from the LV inflow toward the outflow during IVC. In normal subjects, the area of the vortex was sustained, but the flow volume decreased significantly during IVC (P < 0.001). A significant apex-to-base flow velocity gradient was shown along the outflow axis on aortic valve opening. However, both the area and flow volume of the vortex decreased more severely during IVC in the patients (P < 0.001). The apex-to-base flow velocity gradient along the outflow axis disappeared and a reversed velocity gradient was observed at the basal-mid level on aortic valve opening. In multivariate models, a decreased LV ejection fraction was the only independent predictor of the percentage decrease in area of the vortex during the IVC (P < 0.001), and a larger QRS width (P = 0.028) and LV end-systolic long diameter (P = 0.002) were independent predictors of the percentage decrease in flow volume of the vortex.

CONCLUSIONS

The vortex across the LV inflow-outflow region during IVC facilitates the ejection of blood during early systole, and an unsustained vortex may be associated with impaired cardiac function.

Vortex formation time-to-left ventricular early rapid filling relation: model-based prediction with echocardiographic validation

During early rapid filling, blood aspirated by the left ventricle (LV) generates an asymmetric toroidal vortex whose development has been quantified using vortex formation time (VFT), a dimensionless index defined by the length-to-diameter ratio of the aspirated (equivalent cylindrical) fluid column. Since LV wall motion generates the atrioventricular pressure gradient resulting in the early transmitral flow (Doppler E-wave) and associated vortex formation, we hypothesized that the causal relation between VFT and diastolic function (DF), parametrized by stiffness, relaxation, and load, can be elucidated via kinematic modeling. Gharib et al. (Gharib M, Rambod E, Kheradvar A, Sahn DJ, Dabiri JO. Proc Natl Acad Sci USA 103: 6305–6308, 2006) approximated E-wave shape as a triangle and calculated VFTGharib as triangle (E-wave) area (cm) divided by peak (Doppler M-mode derived) mitral orifice diameter (cm). We used a validated kinematic model of filling for the E-wave as a function of time, parametrized by stiffness, viscoelasticity, and load. To calculate VFTkinematic, we computed the curvilinear E-wave area (using the kinematic model) and divided it by peak effective orifice diameter. The derived VFT-to-LV early rapid filling relation predicts VFT to be a function of peak E-wave-to-peak mitral annular tissue velocity (Doppler E′-wave) ratio as (E/E′)3/2. Validation utilized 262 cardiac cycles of simultaneous echocardiographic high-fidelity hemodynamic data from 12 subjects. VFTGharib and VFTkinematic were calculated for each subject and were well-correlated ( R2 = 0.66). In accordance with prediction, VFTkinematic to (E/E′)3/2 relationship was validated ( R 2 = 0.63). We conclude that VFTkinematic is a DF index computable in terms of global kinematic filling parameters of stiffness, viscoelasticity, and load. Validation of the fluid mechanics-to-chamber kinematics relation unites previously unassociated DF assessment methods and elucidates the mechanistic basis of the strong correlation between VFT and (E/E′)3/2.