Fully-coupled fluid-structure interaction simulation of the aortic and mitral valves in a realistic 3D left ventricle model

3D velocity and pressure field reconstruction in the cardiac left ventricle via physics informed neural network from echocardiography guided by 3D color Doppler

Fluid dynamics of the heart chamber can provide critical biological cues for understanding cardiac health and disease and have the potential for supporting diagnosis and prognosis. However, directly acquiring fluid dynamics information from clinical imaging remains challenging, as they are often noisy and have limited resolution, preventing accurate detailed fluid dynamics analysis. Image-based flow simulations offer high detail but are typically difficult to align with clinical velocity measurements, and as a result, may not accurately depict true fluid dynamics. Inverse-computing velocity fields from images via intra-ventricular flow mapping (VFM) has been reported, but it can become inaccurate when faced with missing or noisy measurement data, which is common with modalities such as ultrasound. Here, we propose a physics-informed neural network (PINN) framework that can accurately reconstruct detailed 3D flow fields of the cardiac left ventricle within a localized time window, using supervision from color Doppler measurements, despite their low resolution and signal-to-noise ratio. This framework couples PINN solvers at consecutive time frames with discrete temporal numerical differentiation and is thus named the "Coupled Sequential Frame PINN" or CSF-PINN. We used image-based flow simulations of fetal and adult hearts to generate synthetic color Doppler velocity data at different spatial and temporal resolution for testing the framework. Results show that CSF-PINN can accurately predict high levels of fluid dynamics details, including flow patterns, intraventricular pressure gradients, vorticity structures, and energy losses. CSF-PINN outperforms vanilla PINN in both accuracy and computational efficiency, however, its accuracy is more limited for velocity-gradient-dependent parameters, such as vorticity and wall shear stress (WSS) magnitude. CSF-PINN's accuracy is maintained even when color Doppler velocity data are spatially and temporally sparse and noisy, and when complex motions of the mitral valve are modelled. These are scenarios in which previous methodologies, including image-based flow simulations and VFM, have struggled. Additionally, we propose a scheme for advancing fluid dynamics predictions to subsequent time windows by using training from the previous time window to initialize networks for the subsequent window, further minimizing errors.

Clinical Validation of the PSCOPE Hybrid Model Prediction of Left Ventricular Assist Device Implantation Hemodynamics: Three Patient-Specific Cases

Objective

The Physiology Simulation Coupled Experiment (PSCOPE) is a hybrid modeling framework designed for mechanistic cardiovascular predictive modeling. It couples a physical fluid experiment with a lumped parameter network simulation to replicate the closed-loop feedback between simulated cardiovascular physiology and fluid dynamics in the physical experiment. This study validates PSCOPE's predictions of post-surgical physiology against clinical data in the context of HeartMate 3 left ventricular assist device implantation.

Methods

We designed a protocol to characterize the pre- and post-surgical hemodynamics of three adult HeartMate 3 patients using perioperative clinical measurements acquired from routine intensive care unit monitoring. For each patient, we tuned a lumped parameter network model to match their pre-surgical hemodynamic values, creating a patient-specific simulation of the pre-surgical physiology. The PSCOPE framework then modeled LVAD implantation by coupling these simulations to a physical HeartMate 3 device flow experiment. This hybrid model estimates physiological flow rate and pressure parameters to predict the patients' post-surgical hemodynamics.

Results

The percentage difference between PSCOPE predictions and clinical post-surgical hemodynamics ranged from 0.0% to 44.7% across different hemodynamic parameters in different patients. The predicted cardiac index, mean pulmonary arterial pressure, central venous pressure, and pulmonary arterial wedge pressure together accurately indicated the absence of post-implant right ventricular failure in all patients.

Conclusion

This validation study demonstrates the potential of PSCOPE in assisting LVAD patient management. PSCOPE hemodynamic predictions could help clinicians anticipate and manage post-implant outcomes, such as right ventricular failure, thereby improving the efficacy of surgical planning.

Transcatheter Aortic Valve Replacement in Patients With Quadricuspid Aortic Valve: A Case Series and Systematic Review

Background: Quadricuspid aortic valve (QAV) is a rare congenital cardiac anomaly associated with symptomatic aortic regurgitation (AR) or aortic stenosis (AS). Transcatheter aortic valve replacement (TAVR) for QAV remains uncertain. Methods: We retrospectively reviewed prospectively collected data from patients with QAV undergoing TAVR in our center and conducted a systematic literature review for further investigation. Results: Five patients with QAV were treated with TAVR between April 2016 and December 2023. The median age was 67 years (range: 59-86), and the median Society of Thoracic Surgeons score (STS-score) was 3.750% (range: 0.916%-11.823%). Procedural success was achieved in all cases. The median follow-up period was 3 years (from 30 days to 7 years). Four of the patients exhibited no serious complications, while one experienced delayed coronary obstruction. Our systematic review included 31 cases from 21 publications and our center. The median age of patients was 79 years (range: 57-90), including 18 males. The median STS score was 7.835%. Severe AS was present in 64.5% of the patients and severe AR in 41.9%. The most common QAV subtype was type B (48.4%). Technical success was achieved in 100% of the cases, with two cases reporting coronary obstruction and one required a permanent pacemaker implantation. During a median follow-up period of 1 year (from 30 days to 7 years), one case experienced serious complications of delayed coronary obstruction. Conclusion: The TAVR may be an alternative treatment for patients with QAV, preliminarily demonstrating feasible early and long-term results from current experience. However, extra precautions regarding coronary artery obstruction complications are necessary due to the rarity and anatomical complexity of QAV.

Image-based simulation of mitral valve dynamic closure including anisotropy

Simulation of the dynamic behavior of mitral valve closure could improve clinical treatment by predicting surgical procedures outcome. We propose here a method to achieve this goal by using the immersed boundary method. In order to go towards patient-based simulation, we tailor our method to be adapted to a valve extracted from medical image data. It includes investigating segmentation process, smoothness of geometry, case setup and the shape of the left ventricle. We also study the influence of leaflet tissue anisotropy on the quality of the valve closure by comparing with an isotropic model. As part of the anisotropy analysis, we study the influence of the principal material direction by comparing methods to obtain them without dissection. Results show that our method can be scaled to various image-based data. We evaluate the mitral valve closure quality based on measuring bulging area, contact map, and flow rate. The results show also that the anisotropic material model more precisely represents the physiological characteristics of the valve tissue. Furthermore, results indicate that the orientation of the principal material direction plays a role in the effectiveness of the valve seal.

Multi-Modal in Vitro Experiments Mimicking the Flow Through a Mitral Heart Valve Phantom

PURPOSE

Fluid-structure interaction (FSI) models are more commonly applied in medical research as computational power is increasing. However, understanding the accuracy of FSI models is crucial, especially in the context of heart valve disease in patient-specific models. Therefore, this study aimed to create a multi-modal benchmarking data set for cardiac-inspired FSI models, based on clinically important parameters, such as the pressure, velocity, and valve opening, with an in vitro phantom setup.

METHOD

An in vitro setup was developed with a 3D-printed phantom mimicking the left heart, including a deforming mitral valve. A range of pulsatile flows were created with a computer-controlled motor-and-pump setup. Catheter pressure measurements, magnetic resonance imaging (MRI), and echocardiography (Echo) imaging were used to measure pressure and velocity in the domain. Furthermore, the valve opening was quantified based on cine MRI and Echo images.

RESULT

The experimental setup, with 0.5% cycle-to-cycle variation, was successfully built and six different flow cases were investigated. Higher velocity through the mitral valve was observed for increased cardiac output. The pressure difference across the valve also followed this trend. The flow in the phantom was qualitatively assessed by the velocity profile in the ventricle and by streamlines obtained from 4D phase-contrast MRI.

CONCLUSION

A multi-modal set of data for validation of FSI models has been created, based on parameters relevant for diagnosis of heart valve disease. All data is publicly available for future development of computational heart valve models.

Multiscale structure and function of the aortic valve apparatus

Whereas studying the aortic valve in isolation has facilitated the development of life-saving procedures and technologies, the dynamic interplay of the aortic valve and its surrounding structures is vital to preserving their function across the wide range of conditions encountered in an active lifestyle. Our view is that these structures should be viewed as an integrated functional unit, here referred to as the aortic valve apparatus (AVA). The coupling of the aortic valve and root, left ventricular outflow tract, and blood circulation is crucial for AVA's functions: unidirectional flow out of the left ventricle, coronary perfusion, reservoir function, and support of left ventricular function. In this review, we explore the multiscale biological and physical phenomena that underlie the simultaneous fulfillment of these functions. A brief overview of the tools used to investigate the AVA, such as medical imaging modalities, experimental methods, and computational modeling, specifically fluid-structure interaction (FSI) simulations, is included. Some pathologies affecting the AVA are explored, and insights are provided on treatments and interventions that aim to maintain quality of life. The concepts explained in this article support the idea of AVA being an integrated functional unit and help identify unanswered research questions. Incorporating phenomena through the molecular, micro, meso, and whole tissue scales is crucial for understanding the sophisticated normal functions and diseases of the AVA.

Computational Fluid-Structure Interaction in Microfluidics

Micro elastofluidics is a transformative branch of microfluidics, leveraging the fluid-structure interaction (FSI) at the microscale to enhance the functionality and efficiency of various microdevices. This review paper elucidates the critical role of advanced computational FSI methods in the field of micro elastofluidics. By focusing on the interplay between fluid mechanics and structural responses, these computational methods facilitate the intricate design and optimisation of microdevices such as microvalves, micropumps, and micromixers, which rely on the precise control of fluidic and structural dynamics. In addition, these computational tools extend to the development of biomedical devices, enabling precise particle manipulation and enhancing therapeutic outcomes in cardiovascular applications. Furthermore, this paper addresses the current challenges in computational FSI and highlights the necessity for further development of tools to tackle complex, time-dependent models under microfluidic environments and varying conditions. Our review highlights the expanding potential of FSI in micro elastofluidics, offering a roadmap for future research and development in this promising area.

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

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

A novel mono-physics particle-based approach for the simulation of cardiovascular fluid-structure interaction problems

BACKGROUND AND OBJECTIVE

Fluid-structure interaction (FSI) is required in the study of several cardiovascular engineering problems were the mutual interaction between the pulsatile blood flow and the tissue structures is essential to establish the biomechanics of the system. Traditional FSI methods are partitioned approaches where two independent solvers, one for the fluid and one for the structure, are asynchronously coupled. This process results into high computational costs. In this work, a new FSI scheme which avoids the coupling of different solvers is presented in the framework of the truly incompressible smoothed particle hydrodynamics (ISPH) method.

METHODS

In the proposed FSI method, ISPH particles contribute to define both the fluid and structural domains and are solved together in a unified system. Solid particles, geometrically defined at the beginning of the simulation, are linked through spring bounds with elastic constant providing the material Young's modulus. At each iteration, internal elastic forces are calculated to restore the springs resting length. These forces are added in the predictor step of the fractional-step procedure used to solve the momentum and continuity equations for incompressible flows of all particles.

RESULTS

The method was validated with a benchmark test case consisting of a flexible beam immersed in a channel. Results showed good agreement with the system coupling approach of a well-established commercial software, ANSYS®, both in terms of fluid-dynamics and beam deformation. The approach was then applied to model a complex cardiovascular problem, consisting in the aortic valve operating function. The valve dynamics during opening and closing phases were compared qualitatively with literature results, demonstrating good consistency.

CONCLUSIONS

The method is computationally more efficient than traditional FSI strategies, and overcomes some of their main drawbacks, such as the impossibility of simulating the correct valve coaptation during the closing phase. Thanks to the incompressibility scheme, the proposed FSI method is appropriate to model biological soft tissues. The simplicity and flexibility of the approach also makes it suitable to be expanded for the modelling of thromboembolic phenomena.

Validation of fluid-structure interaction simulations of the opening phase of phantom mitral heart valves under physiologically inspired conditions

BACKGROUND AND OBJECTIVE

Atrioventricular valve disease is a common cause of heart failure, and successful surgical or interventional outcomes are crucial. Patient-specific fluid-structure interaction (FSI) modeling may provide valuable insights into valve dynamics and guidance of valve repair strategies. However, lack of validation has kept FSI modeling from clinical implementation. Therefore, this study aims to validate FSI simulations against in vitro benchmarking data, based on clinically relevant parameters for evaluating heart valve disease.

METHODS

An FSI model that mimics the left heart was developed. The domain included a deformable mitral valve of different stiffnesses run with different inlet velocities. Five different cases were simulated and compared to in vitro data based on the pressure difference across the valve, the valve opening, and the velocity in the flow domain.

RESULTS

The simulations underestimate the pressure difference across the valve by 6.8-14 % compared to catheter measurements. Evaluation of the valve opening showed an underprediction of 5.4-7.3 % when compared to cine MRI, 2D Echo, and 3D Echo data. Additionally, the simulated velocity through the valve showed a 7.9-8.4 % underprediction in relation to Doppler Echo measurements. Qualitative assessment of the velocity profile in the ventricle and the streamlines of the flow in the domain showed good agreement of the flow behavior.

CONCLUSIONS

Parameters relevant to the diagnosis of heart valve disease estimated by FSI simulations showed good agreement when compared to in vitro benchmarking data, with differences small enough not to affect the grading of heart valve disease. The FSI model is thus deemed good enough for further development toward patient-specific cases.

An In Silico Model for Predicting the Efficacy of Edge-to-Edge Repair for Mitral Regurgitation

In recent years, transcatheter edge-to-edge repair (TEER) has been widely adopted as an effective treatment for mitral regurgitation (MR). The aim of this study is to develop a personalized in silico model to predict the effect of edge-to-edge repair in advance to the procedure for each individual patient. For this purpose, we propose a combination of a valve deformation model for computing the mitral valve (MV) orifice area (MVOA) and a lumped parameter model for the hemodynamics, specifically mitral regurgitation volume (RVol). Although we cannot obtain detailed information on the three-dimensional flow field near the mitral valve, we can rapidly simulate the important medical parameters for the clinical decision support. In the present method, we construct the patient-specific pre-operative models by using the parameter optimization and then simulate the postoperative state by applying the additional clipping condition. The computed preclip MVOAs show good agreement with the clinical measurements, and the correlation coefficient takes 0.998. In addition, the MR grade in terms of RVol also has good correlation with the grade by ground truth MVOA. Finally, we try to investigate the applicability for the predicting the postclip state. The simulated valve shapes clearly show the well-known double orifice and the improvement of the MVOA, compared with the preclip state. Similarly, we confirmed the improved reverse flow and MR grade in terms of RVol. A total computational time is approximately 8 h by using general-purpose PC. These results obviously indicate that the present in silico model has good capability for the assessment of edge-to-edge repair.

Preprocedural Planning of Left Atrial Appendage Occlusion: A Review of the Use of Additive Manufacturing

Stroke is a significant public health problem, with non-valvular atrial fibrillation (NVAF) being one of its main causes. This cardiovascular arrhythmia predisposes to the production of intracardiac thrombi, mostly formed in the left atrial appendage (LAA). When there are contraindications to treatment with oral anticoagulants, another therapeutic option to reduce the possibility of thrombus formation in the LAA is the implantation of an occlusion device by cardiac catheterization. The effectiveness of LAA occlusion is dependent on accurate preprocedural device sizing and proper device positioning at the LAA ostium, to ensure sufficient device anchoring and avoid peri-device leaks. Additive manufacturing, commonly known as three-dimensional printing (3DP), of LAA models is beginning to emerge in the scientific literature to address these challenges through procedural simulation. This review aims at clarifying the impact of 3DP on preprocedural planning of LAA occlusion, specifically in the training of cardiac surgeons and in the assessment of the perfect adjustment between the LAA and the biomedical implant.

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

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

Comparison of Immersed Boundary Simulations of Heart Valve Hemodynamics Against In Vitro 4D Flow MRI Data

The immersed boundary (IB) method is a mathematical framework for fluid-structure interaction problems (FSI) that was originally developed to simulate flows around heart valves. Direct comparison of FSI simulations around heart valves against experimental data is challenging, however, due to the difficulty of performing robust and effective simulations, the complications of modeling a specific physical experiment, and the need to acquire experimental data that is directly comparable to simulation data. Such comparators are a necessary precursor for further formal validation studies of FSI simulations involving heart valves. In this work, we performed physical experiments of flow through a pulmonary valve in an in vitro pulse duplicator, and measured the corresponding velocity field using 4D flow MRI (4-dimensional flow magnetic resonance imaging). We constructed a computer model of this pulmonary artery setup, including modeling valve geometry and material properties via a technique called design-based elasticity, and simulated flow through it with the IB method. The simulated flow fields showed excellent qualitative agreement with experiments, excellent agreement on integral metrics, and reasonable relative error in the entire flow domain and on slices of interest. These results illustrate how to construct a computational model of a physical experiment for use as a comparator.

Research on fatigue optimization simulation of polymeric heart valve based on the iterative sub-regional thickened method

Prosthetic polymeric heart valves (PHVs) have the potential to overcome the inherent material and design limitations of traditional valves in the treatment of valvular heart disease; however, their durability remains limited. Optimal design of the valve structure is necessary to improve their durability. This study aimed to enhance the fatigue resistance of PHVs by improving the stress distribution. Iterative subregional thickening of the leaflets was used, and the mechanical stress distribution and hemodynamics of these polymeric tri-leaflet valves were characterized using a fluid-structure interaction approach. Subregional thickening led to a reduction in stress concentration on the leaflet, with the effective orifice area still meeting ISO 5840-3 and the regurgitant volume achieving a similar value to those in previous studies. The maximum stress in the final iteration was reduced by 28% compared with that of the prototype. The proposed method shows potential for analyzing the stress distribution and hemodynamic performance of subregional thickened valves and can further improve the durability of PHVs.

Intra-Cardiac Flow from Geometry Prescribed Computational Fluid Dynamics: Comparison with Ultrasound Vector Flow Imaging

PURPOSE

This paper investigates the accuracy of blood flow velocities simulated from a geometry prescribed computational fluid dynamics (CFD) pipeline by applying it to a dynamic heart phantom. The CFD flow patterns are compared to a direct flow measurement by ultrasound vector flow imaging (VFI). The hypothesis is that the simulated velocity magnitudes are within one standard deviation of the measured velocities.

METHODS

The CFD pipeline uses computed tomography angiography (CTA) images with 20 volumes per cardiac cycle as geometry input. Fluid domain movement is prescribed from volumetric image registration using the CTA image data. Inlet and outlet conditions are defined by the experimental setup. VFI is systematically measured in parallel planes, and compared to the corresponding planes in the simulated time dependent three dimensional fluid velocity field.

RESULTS

The measured VFI and simulated CFD have similar flow patterns when compared qualitatively. A quantitative comparison of the velocity magnitude is also performed at specific regions of interest. These are evaluated at 11 non-overlapping time bins and compared by linear regression giving R2 = 0.809, SD = 0.060 m/s, intercept = - 0.039 m/s, and slope = 1.09. Excluding an outlier at the inlet, the correspondence between CFD and VFI improves to: R2 = 0.823, SD = 0.048 m/s, intercept = -0.030 m/s, and slope = 1.01.

CONCLUSION

The direct comparison of flow patterns shows that the proposed CFD pipeline provide realistic flow patterns in a well-controlled experimental setup. The demanded accuracy is obtained close to the inlet and outlet, but not in locations far from these.

An Image-Based Computational Fluid Dynamics Study of Mitral Regurgitation in Presence of Prolapse

PURPOSE

In this work we performed an imaged-based computational study of the systolic fluid dynamics in presence of mitral valve regurgitation (MVR). In particular, we compared healthy and different regurgitant scenarios with the aim of quantifying different hemodynamic quantities.

METHODS

We performed computational fluid dynamic (CFD) simulations in the left ventricle, left atrium and aortic root, with a resistive immersed method, a turbulence model, and with imposed systolic wall motion reconstructed from Cine-MRI images, which allowed us to segment also the mitral valve. For the regurgitant scenarios we considered an increase of the heart rate and a dilation of the left ventricle.

RESULTS

Our results highlighted that MVR gave rise to regurgitant jets through the mitral orifice impinging against the atrial walls and scratching against the mitral valve leading to high values of wall shear stresses (WSSs) with respect to the healthy case.

CONCLUSION

CFD with prescribed wall motion and immersed mitral valve revealed to be an effective tool to quantitatively describe hemodynamics in case of MVR and to compare different regurgitant scenarios. Our findings highlighted in particular the presence of transition to turbulence in the atrium and allowed us to quantify some important cardiac indices such as cardiac output and WSS.

A CFD study on the interplay of torsion and vortex guidance by the mitral valve on the left ventricular wash-out making use of overset meshes (Chimera technique)

Cardiovascular disease often occurs with silent and gradual alterations of cardiac blood flow that can lead to the onset of chronic pathological conditions. Image-based patient-specific Computational Fluid Dynamics (CFD) models allow for an extensive quantification of the flow field beyond the direct capabilities of medical imaging techniques that could support the clinicians in the early diagnosis, follow-up, and treatment planning of patients. Nonetheless, the large and impulsive kinematics of the left ventricle (LV) and the mitral valve (MV) pose relevant modeling challenges. Arbitrary Lagrangian-Eulerian (ALE) based computational fluid dynamics (CFD) methods struggle with the complex 3D mesh handling of rapidly moving valve leaflets within the left ventricle (LV). We, therefore, developed a Chimera-based (overset meshing) method to build a patient-specific 3D CFD model of the beating LV which includes a patient-inspired kinematic model of the mitral valve (LVMV). Simulations were performed with and without torsion. In addition, to evaluate how the intracardiac LV flow is impacted by the MV leaflet kinematics, a third version of the model without the MV was generated (LV with torsion). For all model versions, six cardiac cycles were simulated. All simulations demonstrated cycle-to-cycle variations that persisted after six cycles but were albeit marginal in terms of the magnitude of standard deviation of velocity and vorticity which may be related to the dissipative nature of the numerical scheme used. The MV was found to have a crucial role in the development of the intraventricular flow by enhancing the direct flow, the apical washout, and the propagation of the inlet jet towards the apical region. Consequently, the MV is an essential feature in the patient-specific CFD modeling of the LV. The impact of torsion was marginal on velocity, vorticity, wall shear stress, and energy loss, whereas it resulted to be significant in the evaluation of particle residence times. Therefore, including torsion could be considered in patient-specific CFD models of the LV, particularly when aiming to study stasis and residence time. We conclude that, despite some technical limitations encountered, the Chimera technique is a promising alternative for ALE methods for 3D CFD models of the heart that include the motion of valve leaflets.

Interplay between Artificial Intelligence and Biomechanics Modeling in the Cardiovascular Disease Prediction

Cardiovascular disease (CVD) is the most common cause of morbidity and mortality worldwide, and early accurate diagnosis is the key point for improving and optimizing the prognosis of CVD. Recent progress in artificial intelligence (AI), especially machine learning (ML) technology, makes it possible to predict CVD. In this review, we first briefly introduced the overview development of artificial intelligence. Then we summarized some ML applications in cardiovascular diseases, including ML-based models to directly predict CVD based on risk factors or medical imaging findings and the ML-based hemodynamics with vascular geometries, equations, and methods for indirect assessment of CVD. We also discussed case studies where ML could be used as the surrogate for computational fluid dynamics in data-driven models and physics-driven models. ML models could be a surrogate for computational fluid dynamics, accelerate the process of disease prediction, and reduce manual intervention. Lastly, we briefly summarized the research difficulties and prospected the future development of AI technology in cardiovascular diseases.

Left Atrial Appendage Occlusion: Current Stroke Prevention Strategies and a Shift Toward Data-Driven, Patient-Specific Approaches

The left atrial appendage (LAA) is a complex structure with unknown physiologic function protruding from the main body of the left atrium. In patients with atrial fibrillation, the left atrium does not contract effectively. Insufficient atrial and LAA contractility predisposes the LAA morphology to hemostasis and thrombus formation, leading to an increased risk of cardioembolic events. Oral anticoagulation therapies are the mainstay of stroke prevention options for patients; however, not all patients are candidates for long-term oral anticoagulation. Percutaneous occlusion devices are an attractive alternative to long-term anticoagulation therapy, although they are not without limitations, such as peri-implant leakage and device-related thrombosis. Although efforts have been made to reduce these risks, significant interpatient heterogeneity inevitably yields some degree of device-anatomy mismatch that is difficult to resolve using current devices and can ultimately lead to insufficient occlusion and poor patient outcomes. In this state-of-the-art review, we evaluated the anatomy of the LAA as well as the current pathophysiologic understanding and stroke prevention strategies used in the management of the risk of stroke associated with atrial fibrillation. We highlighted recent advances in computed tomography imaging, preprocedural planning, computational modeling, and novel additive manufacturing techniques, which represent the tools needed for a paradigm shift toward patient-centric LAA occlusion. Together, we envisage that these techniques will facilitate a pipeline from the imaging of patient anatomy to patient-specific computational and bench-top models that enable customized, data-driven approaches for LAA occlusion that are engineered specifically to meet each patient's unique needs.

Fluid structure computational model of simulating mitral valve motion in a contracting left ventricle

BACKGROUND

Fluid structure interaction simulations h hold promise in studying normal and abnormal cardiac function, including the effect of fluid dynamics on mitral valve (MV) leaflet motion. The goal of this study was to develop a 3D fluid structure interaction computational model to simulate bileaflet MV when interacting with blood motion in left ventricle (LV).

METHODS

The model consists of ideal geometric-shaped MV leaflets and the LV, with MV dimensions based on human anatomical measurements. An experimentally-based hyperelastic isotropic material was used to model the mechanical behaviour of the MV leaflets, with chordae tendineae and papillary muscle tips also incorporated. LV myocardial tissue was prescribed using a transverse isotropic hyperelastic formulation. Incompressible Navier-Stokes fluid formulations were used to govern the blood motion, and the Arbitrary Lagrangian Eulerian (ALE) method was employed to determine the mesh deformation of the fluid and solid domains due to trans-valvular pressure on MV boundaries and the resulting leaflet movement.

RESULTS

The LV-MV generic model was able to reproduce physiological MV leaflet opening and closing profiles resulting from the time-varying atrial and ventricular pressures, as well as simulating normal and prolapsed MV states. Additionally, the model was able to simulate blood flow patterns after insertion of a prosthetic MV with and without left ventricular outflow tract flow obstruction. In the MV-LV normal model, the regurgitant blood flow fraction was 10.1 %, with no abnormality in cardiac function according to the mitral regurgitation severity grades reported by the American Society of Echocardiography.

CONCLUSION

Our simulation approach provides insights into intraventricular fluid dynamics in a contracting LV with normal and prolapsed MV function, as well as aiding in the understanding of possible complications after transcatheter MV implantation prior to clinical trials.

Clinical Impact of Computational Heart Valve Models

This paper provides a review of engineering applications and computational methods used to analyze the dynamics of heart valve closures in healthy and diseased states. Computational methods are a cost-effective tool that can be used to evaluate the flow parameters of heart valves. Valve repair and replacement have long-term stability and biocompatibility issues, highlighting the need for a more robust method for resolving valvular disease. For example, while fluid-structure interaction analyses are still scarcely utilized to study aortic valves, computational fluid dynamics is used to assess the effect of different aortic valve morphologies on velocity profiles, flow patterns, helicity, wall shear stress, and oscillatory shear index in the thoracic aorta. It has been analyzed that computational flow dynamic analyses can be integrated with other methods to create a superior, more compatible method of understanding risk and compatibility.

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.

Computational Methods for Fluid-Structure Interaction Simulation of Heart Valves in Patient-Specific Left Heart Anatomies

Given the complexity of human left heart anatomy and valvular structures, the fluid–structure interaction (FSI) simulation of native and prosthetic valves poses a significant challenge for numerical methods. In this review, recent numerical advancements for both fluid and structural solvers for heart valves in patient-specific left hearts are systematically considered, emphasizing the numerical treatments of blood flow and valve surfaces, which are the most critical aspects for accurate simulations. Numerical methods for hemodynamics are considered under both the continuum and discrete (particle) approaches. The numerical treatments for the structural dynamics of aortic/mitral valves and FSI coupling methods between the solid Ωs and fluid domain Ωf are also reviewed. Future work toward more advanced patient-specific simulations is also discussed, including the fusion of high-fidelity simulation within vivo measurements and physics-based digital twining based on data analytics and machine learning techniques.

A design-based model of the aortic valve for fluid-structure interaction

This paper presents a new method for modeling the mechanics of the aortic valve and simulates its interaction with blood. As much as possible, the model construction is based on first principles, but such that the model is consistent with experimental observations. We require that tension in the leaflets must support a pressure, then derive a system of partial differential equations governing its mechanical equilibrium. The solution to these differential equations is referred to as the predicted loaded configuration; it includes the loaded leaflet geometry, fiber orientations and tensions needed to support the prescribed load. From this configuration, we derive a reference configuration and constitutive law. In fluid-structure interaction simulations with the immersed boundary method, the model seals reliably under physiological pressures and opens freely over multiple cardiac cycles. Further, model closure is robust to extreme hypo- and hypertensive pressures. Then, exploiting the unique features of this model construction, we conduct experiments on reference configurations, constitutive laws and gross morphology. These experiments suggest the following conclusions: (1) The loaded geometry, tensions and tangent moduli primarily determine model function. (2) Alterations to the reference configuration have little effect if the predicted loaded configuration is identical. (3) The leaflets must have sufficiently nonlinear material response to function over a variety of pressures. (4) Valve performance is highly sensitive to free edge length and leaflet height. These conclusions suggest appropriate gross morphology and material properties for the design of prosthetic aortic valves. In future studies, our aortic valve modeling framework can be used with patient-specific models of vascular or cardiac flow.

Intraventricular Flow Simulations in Singular Right Ventricles Reveal Deteriorated Washout and Low Vortex Formation

Purpose

Patients with a functionally univentricular heart represent one of the most common severe cardiac lesions with a prevalence of 3 per 10,000 live births. Hemodynamics of the singular ventricle is a major research topic in cardiology and there exists a relationship between fluid dynamical features and cardiac behavior in health and disease. The aim of the present work was to compare intraventricular flow in single right ventricle (SRV) patients and subjects with healthy left hearts (LV) through patient-specific CFD simulations.

Methods

Three-dimensional real-time echocardiographic images were obtained for five SRV patients and two healthy subjects and CFD simulations with a moving mesh methodology were performed. Intraventricular vortex formation and vortex formation time (VFT) as well as the turbulent kinetic energy (TKE) and ventricular washout were evaluated.

Results

The results show significantly lower values for the VFT and the TKE in SRV patients compared with healthy LV subjects. Furthermore, vortex formation does not progress to the apex in SRV patients. These findings were confirmed by a significantly lower washout in SRV patients.

Conclusions

The study pinpoints the intriguing role of intraventricular flows to characterize performance of SRVs that goes beyond standard clinical metrics such as ejection fraction.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13239-021-00598-9.

Computational Analysis of Wall Shear Stress Patterns on Calcified and Bicuspid Aortic Valves: Focus on Radial and Coaptation Patterns

Calcification and bicuspid valve formation are important aortic valve disorders that disturb the hemodynamics and the valve function. The detailed analysis of aortic valve hemodynamics would lead to a better understanding of the disease’s etiology. We computationally modeled the aortic valve using simplified three-dimensional geometry and inlet velocity conditions obtained via echocardiography. We examined various calcification severities and bicuspid valve formation. Fluid-structure interaction (FSI) analyses were adapted using ANSYS Workbench to incorporate both flow dynamics and leaflet deformation accurately. Simulation results were validated by comparing leaflet movements in B-mode echo recordings. Results indicate that the biomechanical environment is significantly changed for calcified and bicuspid valves. High flow jet velocities are observed in the calcified valves which results in high transvalvular pressure difference (TPG). Wall shear stresses (WSS) increased with the calcification on both fibrosa (aorta side) and ventricularis (left ventricle side) surfaces of the leaflet. The WSS distribution is regular on the ventricularis, as the WSS values proportionally increase from the base to the tip of the leaflet. However, WSS patterns are spatially complex on the fibrosa side. Low WSS levels and spatially complex WSS patterns on the fibrosa side are considered as promoting factors for further calcification and valvular diseases.

A Hybrid Microstructural-Continuum Multiscale Approach for Modeling Hyperelastic Fibrous Soft Tissue

The heterogeneous, nonlinear, anisotropic material behavior of biological tissues makes precise definition of an accurate constitutive model difficult. One possible solution to this issue would be to define microstructural elements and perform fully coupled multiscale simulation. However, for complex geometries and loading scenarios, the computational costs of such simulations can be prohibitive. Ideally then, we should seek a method that contains microstructural detail, but leverages the speed of classical continuum-based finite-element (FE) modeling. In this work, we demonstrate the use of the Holzapfel-Gasser-Ogden (HGO) model [1, 2] to fit the behavior of microstructural network models. We show that Delaunay microstructural networks can be fit to the HGO strain energy function by calculating fiber network strain energy and average fiber stretch ratio. We then use the HGO constitutive model in a FE framework to improve the speed of our hybrid model, and demonstrate that this method, combined with a material property update scheme, can match a full multiscale simulation. This method gives us flexibility in defining complex FE simulations that would be impossible, or at least prohibitively time consuming, in multiscale simulation, while still accounting for microstructural heterogeneity.

Transcatheter Heart Valve Implantation in Bicuspid Patients with Self-Expanding Device

Bicuspid aortic valve (BAV) patients are conventionally not treated by transcathether aortic valve implantation (TAVI) because of anatomic constraint with unfavorable outcome. Patient-specific numerical simulation of TAVI in BAV may predict important clinical insights to assess the conformability of the transcathether heart valves (THV) implanted on the aortic root of members of this challenging patient population. We aimed to develop a computational approach and virtually simulate TAVI in a group of n.6 stenotic BAV patients using the self-expanding Evolut Pro THV. Specifically, the structural mechanics were evaluated by a finite-element model to estimate the deformed THV configuration in the oval bicuspid anatomy. Then, a fluid–solid interaction analysis based on the smoothed-particle hydrodynamics (SPH) technique was adopted to quantify the blood-flow patterns as well as the regions at high risk of paravalvular leakage (PVL). Simulations demonstrated a slight asymmetric and elliptical expansion of the THV stent frame in the BAV anatomy. The contact pressure between the luminal aortic root surface and the THV stent frame was determined to quantify the device anchoring force at the level of the aortic annulus and mid-ascending aorta. At late diastole, PVL was found in the gap between the aortic wall and THV stent frame. Though the modeling framework was not validated by clinical data, this study could be considered a further step towards the use of numerical simulations for the assessment of TAVI in BAV, aiming at understanding patients not suitable for device implantation on an anatomic basis.

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

BACKGROUND

Mitral regurgitation (MR) quantification by the proximal isovelocity surface area (PISA) method remains challenging. Using computer models, the authors evaluated the accuracy of different PISA methods and quantified their errors.

METHODS

Five functional MR computer models of different geometric and tethering abnormalities were created, validated, and treated as phantom models, from which the reference values were directly obtained. Virtual two-dimensional (2D) PISA and three-dimensional (3D) PISA (both peak and integrated values) were performed on these phantom models. By comparing virtual PISA results with reference values, the accuracy of different PISA methods was evaluated, and their sources of errors were quantified.

RESULTS

Compared with reference values of regurgitant flow rate, excellent correlations were found for true PISA (r = 0.99, bias = 32.3 ± 35.3 mL/sec), 3D PISA (r = 0.97, bias = -24.4 ± 55.5 mL/sec), followed by multiplane 2D hemicylindrical PISA (r = 0.88, bias = -24.1 ± 85.4 mL/sec) and hemiellipsoidal PISA (r = 0.91, bias = -55.7 ± 96.6 mL/sec). Weaker correlations were found for single-plane 2D hemispherical PISA (parasternal long-axis: r = 0.71, bias = -77.6 ± 124.5 mL/sec; apical two-chamber: r = 0.69, bias = -52.0 ± 122.0 mL/sec; apical four-chamber: r = 0.82, bias = -65.5 ± 107.3 mL/sec). For regurgitant volume quantification, integrated PISA was more accurate than peak PISA. The bias of 3D PISA improved from -12.7 ± 7.8 mL (peak PISA) to -2.1 ± 5.3 mL (integrated PISA).

CONCLUSIONS

For functional MR quantification, 2D hemispherical PISA had significant underestimation, multiplane 2D hemiellipsoidal and hemicylindrical PISA showed improved accuracy, and 3D PISA was the most accurate. The PISA method is subject to both systematic underestimation due to the Doppler angle effect and systematic overestimation when regurgitant flow is not perpendicular to PISA contour. Integrated PISA is able to capture dynamic MR and is therefore more accurate than peak PISA. The sum of regurgitant flow rates is the most feasible way to perform integrated PISA.

Fluid Flow Characteristics of Healthy and Calcified Aortic Valves Using Three-Dimensional Lagrangian Coherent Structures Analysis

Aortic valve calcification is an important cardiovascular disorder that deteriorates the accurate functioning of the valve leaflets. The increasing stiffness due to the calcification prevents the complete closure of the valve and therefore leads to significant hemodynamic alterations. Computational fluid dynamics (CFD) modeling enables the investigation of the entire flow domain by processing medical images from aortic valve patients. In this study, we computationally modeled and simulated a 3D aortic valve using patient-specific dimensions of the aortic root and aortic sinus. Leaflet stiffness is deteriorated in aortic valve disease due to calcification. In order to investigate the influence of leaflet calcification on flow dynamics, three different leaflet-stiffness values were considered for healthy, mildly calcified, and severely calcified leaflets. Time-dependent CFD results were used for applying the Lagrangian coherent structures (LCS) technique by performing finite-time Lyapunov exponent (FTLE) computations along with Lagrangian particle residence time (PRT) analysis to identify unique vortex structures at the front and backside of the leaflets. Obtained results indicated that the peak flow velocity at the valve orifice increased with the calcification rate. For the healthy aortic valve, a low-pressure field was observed at the leaflet tips. This low-pressure field gradually expanded through the entire aortic sinus as the calcification level increased. FTLE field plots of the healthy and calcified valves showed a variety of differences in terms of flow structures. When the number of fluid particles in the healthy valve model was taken as reference, 1.59 and 1.74 times more particles accumulated in the mildly and severely calcified valves, respectively, indicating that the calcified valves were not sufficiently opened to allow normal mass flow rates.

Personalized intervention cardiology with transcatheter aortic valve replacement made possible with a non-invasive monitoring and diagnostic framework

One of the most common acute and chronic cardiovascular disease conditions is aortic stenosis, a disease in which the aortic valve is damaged and can no longer function properly. Moreover, aortic stenosis commonly exists in combination with other conditions causing so many patients suffer from the most general and fundamentally challenging condition: complex valvular, ventricular and vascular disease (C3VD). Transcatheter aortic valve replacement (TAVR) is a new less invasive intervention and is a growing alternative for patients with aortic stenosis. Although blood flow quantification is critical for accurate and early diagnosis of C3VD in both pre and post-TAVR, proper diagnostic methods are still lacking because the fluid-dynamics methods that can be used as engines of new diagnostic tools are not well developed yet. Despite remarkable advances in medical imaging, imaging on its own is not enough to quantify the blood flow effectively. Moreover, understanding of C3VD in both pre and post-TAVR and its progression has been hindered by the absence of a proper non-invasive tool for the assessment of the cardiovascular function. To enable the development of new non-invasive diagnostic methods, we developed an innovative image-based patient-specific computational fluid dynamics framework for patients with C3VD who undergo TAVR to quantify metrics of: (1) global circulatory function; (2) global cardiac function as well as (3) local cardiac fluid dynamics. This framework is based on an innovative non-invasive Doppler-based patient-specific lumped-parameter algorithm and a 3-D strongly-coupled fluid-solid interaction. We validated the framework against clinical cardiac catheterization and Doppler echocardiographic measurements and demonstrated its diagnostic utility by providing novel analyses and interpretations of clinical data in eleven C3VD patients in pre and post-TAVR status. Our findings position this framework as a promising new non-invasive diagnostic tool that can provide blood flow metrics while posing no risk to the patient. The diagnostic information, that the framework can provide, is vitally needed to improve clinical outcomes, to assess patient risk and to plan treatment.

A Simplified In Silico Model of Left Ventricular Outflow in Patients After Transcatheter Mitral Valve Replacement with Anterior Leaflet Laceration

In silico modeling has been proposed as a tool to simulate left ventricular (LV) outflow tract (LVOT) obstruction in patients undergoing transcatheter mitral valve replacement (TMVR). This study validated a simplified approach to simulate LV outflow hemodynamics in the setting of TMVR with anterior leaflet laceration, a clinical technique used to mitigate the risk of LVOT obstruction. Personalized, 3-dimensional computational fluid dynamics models were developed from computed tomography images of six patients who underwent TMVR with anterior leaflet laceration. LV outflow hemodynamics were simulated using the patient-specific anatomy and the peak systolic flow rate as boundary conditions. The peak outflow velocity, a clinically relevant hemodynamic metric, was extracted from each simulation (v_sim-peak) and compared with the clinical measurement from Doppler echocardiography (v_clin-peak) for validation. In silico models were successfully developed and implemented for all patients. The pre-processing time was 2 h per model and the simulation could be completed within 3 h. In three patients, the lacerated anterior leaflet exposed open cells of the transcatheter valve to flow. Good agreement was obtained between v_sim-peak and v_clin-peak (r = 0.97, p < 0.01) with average discrepancies of 5 ± 2% and 14 ± 1% for patients with exposed and unexposed cells of the transcatheter valve, respectively. The proposed in silico modeling paradigm therefore simulated LV outflow hemodynamics in a time-efficient manner and demonstrated good agreement with clinical measurements. Future studies should investigate the ability of this paradigm to support clinical applications.

Fluid-Structure Interaction Analyses of Biological Systems Using Smoothed-Particle Hydrodynamics

Due to the inherent complexity of biological applications that more often than not include fluids and structures interacting together, the development of computational fluid-structure interaction models is necessary to achieve a quantitative understanding of their structure and function in both health and disease. The functions of biological structures usually include their interactions with the surrounding fluids. Hence, we contend that the use of fluid-structure interaction models in computational studies of biological systems is practical, if not necessary. The ultimate goal is to develop computational models to predict human biological processes. These models are meant to guide us through the multitude of possible diseases affecting our organs and lead to more effective methods for disease diagnosis, risk stratification, and therapy. This review paper summarizes computational models that use smoothed-particle hydrodynamics to simulate the fluid-structure interactions in complex biological systems.

Computational investigation of left ventricular hemodynamics following bioprosthetic aortic and mitral valve replacement

The left ventricle of the heart is a fundamental structure in the human cardiac system that pumps oxygenated blood into the systemic circulation. Several valvular conditions can cause the aortic and mitral valves associated with the left ventricle to become severely diseased and require replacement. However, the clinical outcomes of such operations, specifically the postoperative ventricular hemodynamics of replacing both valves, are not well understood. This work uses computational fluid-structure interaction (FSI) to develop an improved understanding of this effect by modeling a left ventricle with the aortic and mitral valves replaced with bioprostheses. We use a hybrid Arbitrary Lagrangian-Eulerian/immersogeometric framework to accommodate the analysis of cardiac hemodynamics and heart valve structural mechanics in a moving fluid domain. The motion of the endocardium is obtained from a cardiac biomechanics simulation and provided as an input to the proposed numerical framework. The results from the simulations in this work indicate that the replacement of the native mitral valve with a tri-radially symmetric bioprosthesis dramatically changes the ventricular hemodynamics. Most significantly, the vortical motion in the left ventricle is found to reverse direction after mitral valve replacement. This study demonstrates that the proposed computational FSI framework is capable of simulating complex multiphysics problems and can provide an in-depth understanding of the cardiac mechanics.

Computational fluid dynamics of the right atrium: Assessment of modelling criteria for the evaluation of dialysis catheters

Central venous catheters are widely used in haemodialysis therapy, having to respect design requirements for appropriate performance. These are placed within the right atrium (RA); however, there is no prior computational study assessing different catheter designs while mimicking their native environment. Here, a computational fluid dynamics model of the RA, based on realistic geometry and transient physiological boundary conditions, was developed and validated. Symmetric, split and step catheter designs were virtually placed in the RA and their performance was evaluated by: assessing their interaction with the RA haemodynamic environment through prediction of flow vorticity and wall shear stress (WSS) magnitudes (1); and quantifying recirculation and tip shear stress (2). Haemodynamic predictions from our RA model showed good agreement with the literature. Catheter placement in the RA increased average vorticity, which could indicate alterations of normal blood flow, and altered WSS magnitudes and distribution, which could indicate changes in tissue mechanical properties. All designs had recirculation and elevated shear stress values, which can induce platelet activation and subsequently thrombosis. The symmetric design, however, had the lowest associated values (best performance), while step design catheters working in reverse mode were associated with worsened performance. Different tip placements also impacted on catheter performance. Our findings suggest that using a realistically anatomical RA model to study catheter performance and interaction with the haemodynamic environment is crucial, and that care needs to be given to correct tip placement within the RA for improved recirculation percentages and diminished shear stress values.

Precision medicine in human heart modeling : Perspectives, challenges, and opportunities

Precision medicine is a new frontier in healthcare that uses scientific methods to customize medical treatment to the individual genes, anatomy, physiology, and lifestyle of each person. In cardiovascular health, precision medicine has emerged as a promising paradigm to enable cost-effective solutions that improve quality of life and reduce mortality rates. However, the exact role in precision medicine for human heart modeling has not yet been fully explored. Here, we discuss the challenges and opportunities for personalized human heart simulations, from diagnosis to device design, treatment planning, and prognosis. With a view toward personalization, we map out the history of anatomic, physical, and constitutive human heart models throughout the past three decades. We illustrate recent human heart modeling in electrophysiology, cardiac mechanics, and fluid dynamics and highlight clinically relevant applications of these models for drug development, pacing lead failure, heart failure, ventricular assist devices, edge-to-edge repair, and annuloplasty. With a view toward translational medicine, we provide a clinical perspective on virtual imaging trials and a regulatory perspective on medical device innovation. We show that precision medicine in human heart modeling does not necessarily require a fully personalized, high-resolution whole heart model with an entire personalized medical history. Instead, we advocate for creating personalized models out of population-based libraries with geometric, biological, physical, and clinical information by morphing between clinical data and medical histories from cohorts of patients using machine learning. We anticipate that this perspective will shape the path toward introducing human heart simulations into precision medicine with the ultimate goals to facilitate clinical decision making, guide treatment planning, and accelerate device design.

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.

Replication of left ventricular haemodynamics with a simple planar mitral valve model

Tools for the numerical prediction of haemodynamics in multi-disciplinary integrated heart simulations have to be based on computational models that can be solved with low computational effort and still provide physiological flow characteristics. In this context the mitral valve model is important since it strongly influences the flow kinematics, especially during the diastolic phase. In contrast to a 3D valve, a vastly simplified valve model in form of a simple diode is known to be unable to reproduce the characteristic vortex formation and unable to promote a proper ventricular washout. In the present study, an adaptation of the widely used simplest modelling approach for the mitral valve is employed and compared to a physiologically inspired 3D valve within the same ventricular geometry. The adapted approach shows enhanced vortex formation and an improved ventricular washout in comparison to the diode type model. It further shows a high potential in reproducing the main flow characteristics and related particle residence times generated by a 3D valve.

Computer simulations of transapical mitral valve repair with neochordae implantation: Clinical implications

OBJECTIVES

Transapical beating heart neochordae implantation is an innovative mitral valve repair technique that has demonstrated promising clinical results in patients with primary mitral regurgitation. However, as clinical experience continues to increase, neochordae implantation criteria have not been fully standardized. The aim of this study was to investigate the biomechanical effects of selecting an antero-lateral apical access site compared with a postero-lateral site, and suboptimal neochordae length compared with optimal suture length, on restoring physiologic left heart dynamics.

METHODS

Transapical neochordae implantation using 3 and 4 sutures was computer simulated under 3 posterior mitral leaflet prolapse conditions: isolated P2, multiscallop P2/P3 and multiscallop P2/P1. Physiologic, pre- and postrepair left heart dynamics were evaluated using a fluid-structure interaction modeling framework.

RESULTS

Despite the absence of residual mitral regurgitation in all postrepair models with optimal neochordae length, selecting an antero-lateral apical entry site for the treatment of P2/P3 prolapse generated a significant increase (>80%) in neochordae tension and P2 peak stress, with respect to a postero-lateral entry site. During isolated P2 prolapse repair, although neochordae overtension by 5% led to minimal hemodynamic changes in the regurgitant volume compared with using an optimal suture length, a significant increase in systolic and diastolic neochordae tension (>300%) and posterior leaflet average stress (70%-460%) was quantified. On the other hand, neochordae undertension by 5% led to worsening of regurgitation severity.

CONCLUSIONS

This parametric computer study represents a further step toward an improved understanding of the biomechanical outcomes of transapical neochordae technologies.

3D Printing, Computational Modeling, and Artificial Intelligence for Structural Heart Disease

Structural heart disease (SHD) is a new field within cardiovascular medicine. Traditional imaging modalities fall short in supporting the needs of SHD interventions, as they have been constructed around the concept of disease diagnosis. SHD interventions disrupt traditional concepts of imaging in requiring imaging to plan, simulate, and predict intraprocedural outcomes. In transcatheter SHD interventions, the absence of a gold-standard open cavity surgical field deprives physicians of the opportunity for tactile feedback and visual confirmation of cardiac anatomy. Hence, dependency on imaging in periprocedural guidance has led to evolution of a new generation of procedural skillsets, concept of a visual field, and technologies in the periprocedural planning period to accelerate preclinical device development, physician, and patient education. Adaptation of 3-dimensional (3D) printing in clinical care and procedural planning has demonstrated a reduction in early-operator learning curve for transcatheter interventions. Integration of computation modeling to 3D printing has accelerated research and development understanding of fluid mechanics within device testing. Application of 3D printing, computational modeling, and ultimately incorporation of artificial intelligence is changing the landscape of physician training and delivery of patient-centric care. Transcatheter structural heart interventions are requiring in-depth periprocedural understanding of cardiac pathophysiology and device interactions not afforded by traditional imaging metrics.

A hybrid echocardiography-CFD framework for ventricular flow simulations

Image-based CFD is a powerful tool to study cardiovascular flows while 2D echocardiography (echo) is the most widely used noninvasive imaging modality for the diagnosis of heart disease. Here, echo is combined with CFD, that is, an echo-CFD framework, to study ventricular flows. To achieve this, the previous 3D reconstruction from multiple 2D echo at standard cross sections is extended by: (a) reconstructing aortic and mitral valves from 2D echo and closing the left-ventricle (LV) geometry by approximating a superior wall; (b) incorporating the physiological assumption of the fixed apex as a reference (fixed) point in the 3D reconstruction; and (c) incorporating several smoothing algorithms to remove the nonphysical oscillations (ringing) near the basal section. The method is applied to echo from a baseline LV and one after inducing acute myocardial ischemia (AMI). The 3D reconstruction is validated by comparing it against a reference reconstruction from many echo sections while flow simulations are validated against the Doppler ultrasound velocity measurements. The sensitivity study shows that the choice of the smoothing algorithm does not change the flow pattern inside the LV. However, the presence of the mitral valve can significantly change the flow pattern during the diastole phase. In addition, the abnormal shape of a LV with AMI can drastically change the flow during diastole. Furthermore, the hemodynamic energy loss, as an indicator of the LV pumping performance, for different test cases is calculated, which shows a larger energy loss for a LV with AMI compared to the baseline one.

A Comprehensive Engineering Analysis of Left Heart Dynamics After MitraClip in a Functional Mitral Regurgitation Patient

Percutaneous edge-to-edge mitral valve (MV) repair using MitraClip has been recently established as a treatment option for patients with heart failure and functional mitral regurgitation (MR), which significantly expands the number of patients that can be treated with this device. This study aimed to quantify the morphologic, hemodynamic and structural changes, and evaluate the biomechanical interaction between the MitraClip and the left heart (LH) complex of a heart failure patient with functional MR using a fluid-structure interaction (FSI) modeling framework. MitraClip implantation using lateral, central and double clip positions, as well as combined annuloplasty procedures were simulated in a patient-specific LH model that integrates detailed anatomic structures, incorporates age- and gender-matched non-linear elastic material properties, and accounts for mitral chordae tethering. Our results showed that antero-posterior distance, mitral annulus spherecity index, anatomic regurgitant orifice area, and anatomic opening orifice area decreased by up to 28, 39, 52, and 71%, respectively, when compared to the pre-clip model. MitraClip implantation immediately decreased the MR severity and improved the hemodynamic profile, but imposed a non-physiologic configuration and loading on the mitral apparatus, with anterior and posterior leaflet stress significantly increasing up to 210 and 145% during diastole, respectively. For this patient case, while implanting a combined central clip and ring resulted in the highest reduction in the regurgitant volume (46%), this configuration also led to mitral stenosis. Patient-specific computer simulations as used here can be a powerful tool to examine the complex device-host biomechanical interaction, and may be useful to guide device positioning for potential favorable clinical outcomes.

New Insights into Valve Hemodynamics

Heart valve diseases are common disorders with five million annual diagnoses being made in the United States alone. All heart valve disorders alter cardiac hemodynamic performance; therefore, treatments aim to restore normal flow. This paper reviews the state-of-the-art clinical and engineering advancements in heart valve treatments with a focus on hemodynamics. We review engineering studies and clinical literature on the experience with devices for aortic valve treatment, as well as the latest advancements in mitral valve treatments and the pulmonic and tricuspid valves on the right side of the heart. Upcoming innovations will potentially revolutionize treatment of heart valve disorders. These advancements, and more gradual enhancements in the procedural techniques and imaging modalities, could improve the quality of life of patients suffering from valvular disease who currently cannot be treated.

The Impact of Self-Expandable Transcatheter Aortic Valve Replacement on Concomitant Functional Mitral Regurgitation: A Comprehensive Engineering Analysis

Background

Mitral regurgitation (MR) is present in a large proportion of patients who undergo transcatheter aortic valve replacement (TAVR). However, existing clinical data on the impact of TAVR on early post-procedural MR severity are contradictory. Using a comprehensive computational engineering methodology, this study aimed to evaluate quantitatively the structural and hemodynamic impact of TAVR on aortic-mitral continuity and MR severity in a rigorously developed and validated patient-specific left heart (LH) computer model with aortic stenosis and concomitant functional MR.

Methods

TAVR procedure was virtually simulated using a self-expandable valve (SEV) at three implantation heights. Pre- and post-TAVR LH dynamics as well as intra-operative biomechanics were analyzed.

Results

No significant differences in early MR improvement (<10%) were noted at the three implantation depths when compared to the pre-TAVR state. The high deployment model resulted in the highest stress in the native aortic leaflets, lowest stent-tissue contact force, highest aortic-mitral angle, and highest MR reduction for this patient case. When comparing SEV vs. balloon-expandable valve (BEV) performance at an optimal implantation height, the SEV gave a higher regurgitant volume ⋅ than the pre-TAVR model (40.49 vs 37.59 ml), while the BEV model gave the lowest regurgitant volume (33.84 vs 37.59 ml).

Conclusions

Contact force, aortic-mitral angle, and valve annuli compression were identified as possible mechanistic parameters that may suggest avenues for acute MR improvement. Albeit a single patient parametric study, it is our hope that such detailed engineering analysis could shed some light into the underlying biomechanical mechanisms of TAVR impact on MR.

Transapical mitral valve repair with neochordae implantation: FSI analysis of neochordae number and complexity of leaflet prolapse

Transapical mitral valve repair with neochordae implantation is a relatively new minimally invasive technique to treat primary mitral regurgitation. Quantifying the complex biomechanical interaction and interdependence between the left heart structures and the neochordae during this procedure is technically challenging. The aim of this parametric computational study is to investigate the immediate effects of neochordae number and complexity of leaflet prolapse on restoring physiologic left heart dynamics after optimal transapical neochordae repair procedures. Neochordae implantation using three and four sutures was modeled under three clinically relevant prolapse conditions: isolated P2, multi-scallop P2/P3, and multi-scallop P2/P1. A fluid-structure interaction (FSI) modeling framework was used to evaluate the left heart dynamics under baseline, prerepair, and postrepair states. Despite immediate restoration of leaflet coaptation and no residual mitral regurgitation in all postrepair models, the average and peak stresses in the repaired scallop(s) increased >40% and >100%, respectively, compared with the baseline state. Additionally, anterior mitral leaflet marginal chordae tension increased >30%, while posterior mitral leaflet chordae tension decreased at least 30%. No marked differences in hemodynamic performance, in native and neochordae forces, and in leaflet stress were found when implanting three or four sutures. We report, to our knowledge, the first set of time-dependent in silico FSI human neochordae tension measurements during transapical neochordae repair. This work represents a further step towards an improved understanding of the biomechanical outcomes of minimally invasive mitral valve repair procedures.

A biomechanical model of the pathological aortic valve: simulation of aortic stenosis

Aortic stenosis (AS) disease is a narrowing of the aortic valve (AV) opening which reduces blood flow from the heart causing several health complications. Although a lot of work has been done in AV simulations, most of the efforts have been conducted regarding healthy valves. In this article, a new three-dimensional patient-specific biomechanical model of the valve, based on a parametric formulation of the stenosis that permits the simulation of different degrees of pathology, is presented. The formulation is based on a double approach: the first one is done from the geometric point of view, reducing the effective ejection area of the AV by joining leaflets using a zipper effect to sew them; the second one, in terms of functionality, is based on the modification of AV tissue properties due to the effect of calcifications. Both healthy and stenotic valves were created using patient-specific data and results of the numerical simulation of the valve function are provided. Analysis of the results shows a variation in the first principal stress, geometric orifice area, and blood velocity which were validated against clinical data. Thus, the possibility to create a pipeline which allows the integration of patient-specific data from echocardiographic images and iFR studies to perform finite elements analysis is proved.