Two-dimensional intraventricular flow pattern visualization using the image-based computational fluid dynamics

A variable heart rate multi-compartmental coupled model of the cardiovascular and respiratory systems

Current mathematical models of the cardiovascular system that are based on systems of ordinary differential equations are limited in their ability to mimic important features of measured patient data, such as variable heart rates (HR). Such limitations present a significant obstacle in the use of such models for clinical decision-making, as it is the variations in vital signs such as HR and systolic and diastolic blood pressure that are monitored and recorded in typical critical care bedside monitoring systems. In this paper, novel extensions to well-established multi-compartmental models of the cardiovascular and respiratory systems are proposed that permit the simulation of variable HR. Furthermore, a correction to current models is also proposed to stabilize the respiratory behaviour and enable realistic simulation of vital signs over the longer time scales required for clinical management. The results of the extended model developed here show better agreement with measured bio-signals, and these extensions provide an important first step towards estimating model parameters from patient data, using methods such as neural ordinary differential equations. The approach presented is generalizable to many other similar multi-compartmental models of the cardiovascular and respiratory systems.

Computational Assessment of Valvular Dysfunction in Discrete Subaortic Stenosis: A Parametric Study

PURPOSE

Discrete subaortic stenosis (DSS) is a left-ventricular outflow tract (LVOT) obstruction caused by a membranous lesion. DSS is associated with steep aortoseptal angles (AoSAs) and is a risk factor for aortic regurgitation (AR). However, the etiology of AR secondary to DSS remains unknown. This study aimed at quantifying computationally the impact of AoSA steepening and DSS on aortic valve (AV) hemodynamics and AR.

METHODS

An LV geometry reconstructed from cine-MRI data was connected to an AV geometry to generate a unified 2D LV-AV model. Six geometrical variants were considered: unobstructed (CTRL) and DSS-obstructed LVOT (DSS), each reflecting three AoSA variations (110°, 120°, 130°). Fluid-structure interaction simulations were run to compute LVOT flow, AV leaflet dynamics, and regurgitant fraction (RF).

RESULTS

AoSA steepening and DSS generated vortex dynamics alterations and stenotic flow conditions. While the CTRL-110° model generated the highest degree of leaflet opening asymmetry, DSS preferentially altered superior leaflet kinematics, and caused leaflet-dependent alterations in systolic fluttering. LVOT steepening and DSS subjected the leaflets to increasing WSS overloads (up to 94% increase in temporal shear magnitude), while DSS also increased WSS bidirectionality on the inferior leaflet belly (+ 0.30-point in oscillatory shear index). Although AoSA steepening and DSS increased diastolic transvalvular backflow, regurgitant fractions (RF < 7%) remained below the threshold defining clinical mild AR.

CONCLUSIONS

The mechanical interactions between AV leaflets and LVOT steepening/DSS hemodynamic derangements do not cause AR. However, the leaflet WSS abnormalities predicted in those anatomies provide new support to a mechanobiological etiology of AR secondary to DSS.

Blood flow patterns estimation in the left ventricle with low-rate 2D and 3D dynamic contrast-enhanced ultrasound

BACKGROUND AND OBJECTIVE

Left ventricle (LV) dysfunction always occurs at early heart-failure stages, producing variations in the LV flow patterns. Cardiac diagnostics may therefore benefit from flow-pattern analysis. Several visualization tools have been proposed that require ultrafast ultrasound acquisitions. However, ultrafast ultrasound is not standard in clinical scanners. Meanwhile techniques that can handle low frame rates are still lacking. As a result, the clinical translation of these techniques remains limited, especially for 3D acquisitions where the volume rates are intrinsically low.

METHODS

To overcome these limitations, we propose a novel technique for the estimation of LV blood velocity and relative-pressure fields from dynamic contrast-enhanced ultrasound (DCE-US) at low frame rates. Different from other methods, our method is based on the time-delays between time-intensity curves measured at neighbor pixels in the DCE-US loops. Using Navier-Stokes equation, we regularize the obtained velocity fields and derive relative-pressure estimates. Blood flow patterns were characterized with regard to their vorticity, relative-pressure changes (dp/dt) in the LV outflow tract, and viscous energy loss, as these reflect the ejection efficiency.

RESULTS

We evaluated the proposed method on 18 patients (9 responders and 9 non-responders) who underwent cardiac resynchronization therapy (CRT). After CRT, the responder group evidenced a significant (p<0.05) increase in vorticity and peak dp/dt, and a non-significant decrease in viscous energy loss. No significant difference was found in the non-responder group. Relative feature variation before and after CRT evidenced a significant difference (p<0.05) between responders and non-responders for vorticity and peak dp/dt. Finally, the method feasibility is also shown with 3D DCE-US.

CONCLUSIONS

Using the proposed method, adequate visualization and quantification of blood flow patterns are successfully enabled based on low-rate DCE-US of the LV, facilitating the clinical adoption of the method using standard ultrasound scanners. The clinical value of the method in the context of CRT is also shown.

Colour-Doppler echocardiography flow field velocity reconstruction using a streamfunction-vorticity formulation

We introduce a new method (Doppler Velocity Reconstruction or DoVeR), for reconstructing two-component velocity fields from colour Doppler scans. DoVeR employs the streamfunction-vorticity equation, which satisfies mass conservation while accurately approximating the flow rate of rotation. We validated DoVeR using artificial colour Doppler images generated from computational fluid dynamics models of left ventricle (LV) flow. We compare DoVeR against the conventional intraventricular vector flow mapping (iVFM1D) and reformulated iVFM (iVFM2D). LV model error analysis showed that DoVeR is more robust to noise and probe placement, with noise RMS errors (nRMSE) between 3.81% and 6.67%, while the iVFM methods delivered 4.16-24.17% for iVFM1D and 4.06-400.21% for iVFM2D. We test the DoVeR and iVFM methods using in vivo mouse LV ultrasound scans. DoVeR yielded more haemodynamically accurate reconstructions, suggesting that it can provide a more reliable approach for robust quantification of cardiac flow.

Impact of Aortoseptal Angle Abnormalities and Discrete Subaortic Stenosis on Left-Ventricular Outflow Tract Hemodynamics: Preliminary Computational Assessment

Discrete subaortic stenosis (DSS) is an obstruction of the left ventricular outflow tract (LVOT) due to the formation of a fibromuscular membrane upstream of the aortic valve. DSS is a major risk factor for aortic regurgitation (AR), which often persists after surgical resection of the membrane. While the etiology of DSS and secondary AR is largely unknown, the frequent association between DSS and aortoseptal angle (AoSA) abnormalities has supported the emergence of a mechanobiological pathway by which hemodynamic stress alterations on the septal wall could trigger a biological cascade leading to fibrosis and membrane formation. The resulting LVOT flow disturbances could activate the valve endothelium and contribute to AR. In an effort to assess this hypothetical mechano-etiology, this study aimed at isolating computationally the effects of AoSA abnormalities on septal wall shear stress (WSS), and the impact of DSS on LVOT hemodynamics. Two-dimensional computational fluid dynamics models featuring a normal AoSA (N-LV), a steep AoSA (S-LV), and a steep AoSA with a DSS lesion (DSS-LV) were designed to compute the flow in patient-specific left ventricles (LVs). Boundary conditions consisted of transient velocity profiles at the mitral inlet and LVOT outlet, and patient-specific LV wall motion. The deformation of the DSS lesion was computed using a two-way fluid-structure interaction modeling strategy. Turbulence was accounted for via implementation of the k-ω turbulence model. While the N-LV and S-LV models generated similar LVOT flow characteristics, the DSS-LV model resulted in an asymmetric LVOT jet-like structure, subaortic stenotic conditions (up to 2.4-fold increase in peak velocity, 45% reduction in effective jet diameter vs. N-LV/S-LV), increased vorticity (2.8-fold increase) and turbulence (5- and 3-order-of-magnitude increase in turbulent kinetic energy and Reynolds shear stress, respectively). The steep AoSA subjected the septal wall to a 23% and 69% overload in temporal shear magnitude and gradient, respectively, without any substantial change in oscillatory shear index. This study reveals the existence of WSS overloads on septal wall regions prone to DSS lesion formation in steep LVOTs, and the development of highly turbulent, stenotic and asymmetric flow in DSS LVOTs, which support a possible mechano etiology for DSS and secondary AR.

Investigation of left heart flow using a numerical correlation to model heart wall motion

A three-dimensional computational fluid dynamics (CFD) method has been developed to model the flow in the left heart including atrium and ventricle. Since time resolution of the medical scans does not fit the requirements of the CFD calculations, the main challenge in a numerical simulation of heart chambers is wall motion modeling. This study employs a novel three-dimensional approximation scheme to correlate the wall boundary and grid movement in systole and diastole. It uses a geometry extracted from medical images in the literature and deformed based on the reported flow rates. The opening and closing of the mitral (MV) and the aortic valve (AV) considered as simultaneous events. Unstructured tetragonal grids were used for the meshing of the domain. The calculation was performed by a Navier-Stokes solver using the arbitrary Lagrange-Euler (ALE) formulation. Results show that the proposed correlation for the wall motion could predict the main features of heart flows.

Application of Patient-Specific Computational Fluid Dynamics in Coronary and Intra-Cardiac Flow Simulations: Challenges and Opportunities

The emergence of new cardiac diagnostics and therapeutics of the heart has given rise to the challenging field of virtual design and testing of technologies in a patient-specific environment. Given the recent advances in medical imaging, computational power and mathematical algorithms, patient-specific cardiac models can be produced from cardiac images faster, and more efficiently than ever before. The emergence of patient-specific computational fluid dynamics (CFD) has paved the way for the new field of computer-aided diagnostics. This article provides a review of CFD methods, challenges and opportunities in coronary and intra-cardiac flow simulations. It includes a review of market products and clinical trials. Key components of patient-specific CFD are covered briefly which include image segmentation, geometry reconstruction, mesh generation, fluid-structure interaction, and solver techniques.