Simulating impaired left ventricular-arterial coupling in aging and disease: a systematic review.

Aortic stenosis, hypertension, and left ventricular hypertrophy often coexist in the elderly, causing a detrimental mismatch in coupling between the heart and vasculature known as ventricular-vascular (VA) coupling. Impaired left VA coupling, a critical aspect of cardiovascular dysfunction in aging and disease, poses significant challenges for optimal cardiovascular performance. This systematic review aims to assess the impact of simulating and studying this coupling through computational models. By conducting a comprehensive analysis of 34 relevant articles obtained from esteemed databases such as Web of Science, Scopus, and PubMed until July 14, 2022, we explore various modeling techniques and simulation approaches employed to unravel the complex mechanisms underlying this impairment. Our review highlights the essential role of computational models in providing detailed insights beyond clinical observations, enabling a deeper understanding of the cardiovascular system. By elucidating the existing models of the heart (3D, 2D, and 0D), cardiac valves, and blood vessels (3D, 1D, and 0D), as well as discussing mechanical boundary conditions, model parameterization and validation, coupling approaches, computer resources and diverse applications, we establish a comprehensive overview of the field. The descriptions as well as the pros and cons on the choices of different dimensionality in heart, valve, and circulation are provided. Crucially, we emphasize the significance of evaluating heart-vessel interaction in pathological conditions and propose future research directions, such as the development of fully coupled personalized multidimensional models, integration of deep learning techniques, and comprehensive assessment of confounding effects on biomarkers.

Remodeling in Aortic Stenosis With Reduced and Preserved Ejection Fraction: Insight on Motion Abnormality Via 3D + Time Personalized LV Modeling in Cardiac MRI.

BACKGROUND: Increased afterload in aortic stenosis (AS) induces left ventricle (LV) remodeling to preserve a normal ejection fraction. This compensatory response can become maladaptive and manifest with motion abnormality. It is a clinical challenge to identify contractile and relaxation dysfunction during early subclinical stage to prevent irreversible deterioration. PURPOSE: To evaluate the changes of regional wall dynamics in 3D + time domain as remodeling progresses in AS. STUDY TYPE: Retrospective. POPULATION: A total of 31 AS patients with reduced and preserved ejection fraction (14 AS_rEF: 7 male, 66.5 [7.8] years old; 17 AS_pEF: 12 male, 67.0 [6.0] years old) and 15 healthy (6 male, 61.0 [7.0] years old). FIELD STRENGTH/SEQUENCE: 1.5 T Magnetic resonance imaging/steady state free precession and late-gadolinium enhancement sequences. ASSESSMENT: Individual LV models were reconstructed in 3D + time domain and motion metrics including wall thickening (TI), dyssynchrony index (DI), contraction rate (CR), and relaxation rate (RR) were automatically extracted and associated with the presence of scarring and remodeling. STATISTICAL TESTS: Shapiro-Wilk: data normality; Kruskal-Wallis: significant difference (P < 0.05); ICC and CV: variability; Mann-Whitney: effect size. RESULTS: AS_rEF group shows distinct deterioration of cardiac motions compared to AS_pEF and healthy groups (TIAS_rEF : 0.92 [0.85] mm, TIAS_pEF : 5.13 [1.99] mm, TIhealthy : 3.61 [1.09] mm, ES: 0.48-0.83; DIAS_rEF : 17.11 [7.89]%, DIAS_pEF : 6.39 [4.04]%, DIhealthy : 5.71 [1.87]%, ES: 0.32-0.85; CRAS_rEF : 8.69 [6.11] mm/second, CRAS_pEF : 16.48 [6.70] mm/second, CRhealthy : 10.82 [4.57] mm/second, ES: 0.29-0.60; RRAS_rEF : 8.45 [4.84] mm/second; RRAS_pEF : 13.49 [8.56] mm/second, RRhealthy : 9.31 [2.48] mm/second, ES: 0.14-0.43). The difference in the motion metrics between healthy and AS_pEF groups were insignificant (P-value = 0.16-0.72). AS_rEF group was dominated by eccentric hypertrophy (47.1%) with concomitant scarring. Conversely, AS_pEF group was dominated by concentric remodeling and hypertrophy (71.4%), which could demonstrate hyperkinesia with slight wall dyssynchrony than healthy. Dysfunction of LV mechanics corresponded to the presence of myocardial scarring (54.9% in AS), which reverted the compensatory mechanisms initiated and performed by LV remodeling. DATA CONCLUSION: The proposed 3D + time modeling technique may distinguish regional motion abnormalities between AS_pEF, AS_rEF, and healthy cohorts, aiding clinical diagnosis and monitoring of AS progression. Subclinical myocardial dysfunction is evident in early AS despite of normal EF. LEVEL OF EVIDENCE: 4 TECHNICAL EFFICACY: Stage 1.

Assessing Complex Left Ventricular Adaptations in Aortic Stenosis Using Personalized 3D + time Cardiac MRI Modeling.

Left ventricular adaptations can be a complex process under the influence of aortic stenosis (AS) and comorbidities. This study proposed and assessed the feasibility of using a motion-corrected personalized 3D + time LV modeling technique to evaluate the adaptive and maladaptive LV response to aid treatment decision-making. A total of 22 AS patients were analyzed and compared against 10 healthy subjects. The 3D + time analysis showed a highly distinct and personalized pattern of remodeling in individual AS patients which is associated with comorbidities and fibrosis. Patients with AS alone showed better wall thickening and synchrony than those comorbid with hypertension. Ischemic heart disease in AS caused impaired wall thickening and synchrony and systolic function. Apart from showing significant correlations to echocardiography and clinical MRI measurements (r: 0.70-0.95; p < 0.01), the proposed technique helped in detecting subclinical and subtle LV dysfunction, providing a better approach to evaluate AS patients for specific treatment, surgical planning, and follow-up recovery.

Comparison of diametric and volumetric changes in Stanford type B aortic dissection patients in assessing aortic remodeling post-stent graft treatment.

BACKGROUND: The study aims to analyze the correlation between the maximal diameter (both axial and orthogonal) and volume changes in the true (TL) and false lumens (FL) after stent-grafting for Stanford type B aortic dissection. METHOD: Computed tomography angiography was performed on 13 type B aortic dissection patients before and after procedure, and at 6 and 12 months follow-up. The lumens were divided into three regions: the stented area (Region 1), distal to the stent graft to the celiac artery (Region 2), and between the celiac artery and the iliac bifurcation (Region 3). Changes in aortic morphology were quantified by the increase or decrease of diametric and volumetric percentages from baseline measurements. RESULTS: At Region 1, the TL diameter and volume increased (pre-treatment: volume =51.4±41.9 mL, maximal axial diameter =22.4±6.8 mm, maximal orthogonal diameter =21.6±7.2 mm; follow-up: volume =130.7±69.2 mL, maximal axial diameter =40.1±8.1 mm, maximal orthogonal diameter =31.9+2.6 mm, P<0.05 for all comparisons), while FL decreased (pre-treatment: volume =129.6±150.5 mL; maximal axial diameter =43.0±15.8 mm; maximal orthogonal diameter =28.3±12.6 mm; follow-up: volume =66.6±95.0 mL, maximal axial diameter =24.5±19.9 mm, maximal orthogonal diameter =16.9±13.7, P<0.05 for all comparisons). Due to the uniformity in size throughout the vessel, high concordance was observed between diametric and volumetric measurements in the stented region with 93% and 92% between maximal axial diameter and volume for the true/false lumens, and 90% and 92% between maximal orthogonal diameter and volume for the true/false lumens. Large discrepancies were observed between the different measurement methods at regions distal to the stent graft, with up to 46% differences between maximal orthogonal diameter and volume. CONCLUSIONS: Volume measurement was shown to be a much more sensitive indicator in identifying lumen expansion/shrinkage at the distal stented region.

The study of myocardial ischemia-reperfusion treatment through computational modelling.

Reperfusion of the blood flow to ischemic myocardium is the standard treatment for patients suffering myocardial infarction. However, the reperfusion itself can also induce myocardial injury, in which the actual mechanism and its risk factors remain unclear. This work aims to study the mechanism of ischemia-reperfusion treatment using a three-dimensional (3D) oxygen diffusion model. An electrical model is then coupled to an oxygen model to identify the possible region of myocardial damage. Our findings show that the value of oxygen exceeds its optimum (>1.0) at the ischemic area during early reperfusion period. This complication was exacerbated in a longer ischemic period. While a longer reperfusion time causes a continuous excessive oxygen supply to the ischemic area throughout the reperfusion time. This work also suggests the use of less than 0.8 of initial oxygen concentration in the reperfusion treatment to prevent undesired upsurge at the early reperfusion period and further myocardial injury. We also found the region at risk for myocardial injury is confined in the ischemic vicinity revealed by its electrical conductivity impairment. Although there is a risk that reperfusion leads to myocardial injury for excessive oxygen accumulation, the reperfusion treatment is helpful in reducing the infarct size.

Phenotyping of hypertensive heart disease and hypertrophic cardiomyopathy using personalized 3D modelling and cardiac cine MRI.

Differential diagnosis of hypertensive heart disease (HHD) and hypertrophic cardiomyopathy (HCM) is clinically challenging but important for treatment management. This study aims to phenotype HHD and HCM in 3D + time domain by using a multiparametric motion-corrected personalized modeling algorithm and cardiac magnetic resonance (CMR). 44 CMR data, including 12 healthy, 16 HHD and 16 HCM cases, were examined. Multiple CMR phenotype data consisting of geometric and dynamic variables were extracted globally and regionally from the models over a full cardiac cycle for comparison against healthy models and clinical reports. Statistical classifications were used to identify the distinctive characteristics and disease subtypes with overlapping functional data, providing insights into the challenges for differential diagnosis of both types of disease. While HCM is characterized by localized extreme hypertrophy of the LV, wall thickening/contraction/strain was found to be normal and in sync, though it was occasionally exaggerated at normotrophic/less severely hypertrophic regions during systole to preserve the overall ejection fraction (EF) and systolic functionality. Additionally, we observed that hypertrophy in HHD could also be localized, although at less extreme conditions (i.e. more concentric). While fibrosis occurs mostly in those HCM cases with aortic obstruction, only minority of HHD patients were found affected by fibrosis. We demonstrate that subgroups of HHD (i.e. preserved and reduced EF: HHDpEF & HHDrEF) have different 3D + time CMR characteristics. While HHDpEF has cardiac functions in normal range, dilation and heart failure are indicated in HHDrEF as reflected by low LV wall thickening/contraction/strain and synchrony, as well as much reduced EF.

Effect of intimal flap motion on flow in acute type B aortic dissection by using fluid-structure interaction.

A monolithic, fully coupled fluid-structure interaction (FSI) computational framework was developed to account for dissection flap motion in acute type B aortic dissection (TBAD). Analysis of results included wall deformation, pressure, flow, wall shear stress (WSS), von Mises stress and comparison of hemodynamics between rigid wall and FSI models. Our FSI model mimicked realistic wall deformation that resulted in maximum compression of the distal true lumen (TL) by 21.4%. The substantial movement of intimal flap mostly affected flow conditions in the false lumen (FL). Flap motion facilitated more flow entering the FL at peak systole, with the TL to FL flow split changing from 88:12 in the rigid model to 83:17 in the FSI model. There was more disturbed flow in the FL during systole (5.8% FSI vs 5.2% rigid) and diastole (13.5% FSI vs 9.8% rigid), via a λ2 -criterion. The flap-induced disturbed flow near the tears in the FSI model caused an increase of local WSS by up to 70.0% during diastole. This resulted in a significant reduction in the size of low time-averaged WSS (TAWSS) regions in the FL (113.11 cm2 FSI vs 177.44 cm2 rigid). Moreover, the FSI model predicted lower systolic pressure, higher diastolic pressure, and hence lower pulse pressure. Our results provided new insights into the possible impact of flap motion on flow in aortic dissections, which are particularly important when evaluating hemodynamics of acute TBAD. NOVELTY STATEMENT: Our monolithic fully coupled FSI computational framework is able to reproduce experimentally measured range of flap deformation in aortic dissection, thereby providing novel insights into the influence of physiological flap motion on the flow and pressure distributions. The drastic flap movement increases the flow resistance in the true lumen and causes more disturbed flow in the false lumen, as visualized through the λ2 criterion. The flap-induced luminal pressure is dampened, thereby affecting pressure measures, which may serve as potential prognostic indicators for late complications in acute uncomplicated TBAD patients.

The role of end-diastolic myocardial fibre stretch on infarct extension.

Myocardial infarct extension, a process involving the enlargement of infarct and border zone, leads to progressive degeneration of left ventricular (LV) function and eventually gives rise to heart failure. Despite carrying a high risk, the causation of infarct extension is still a subject of much speculation. In this study, patient-specific LV models were developed to investigate the correlation between infarct extension and impaired regional mechanics. Subsequently, sensitivity analysis was performed to examine the causal factors responsible for the impaired regional mechanics observed in regions surrounding the infarct and border zone. From our simulations, fibre strain, fibre stress and fibre stress-strain loop (FSSL) were the key biomechanical variables affected in these regions. Among these variables, only FSSL was correlated with infarct extension, as reflected in its work density dissipation (WDD) index value, with high WDD indices recorded at regions with infarct extension. Impaired FSSL is caused by inadequate contraction force generation during the isovolumic contraction and ejection phases. Our further analysis revealed that the inadequacy in contraction force generation is not necessarily due to impaired myocardial intrinsic contractility, but at least in part, due to inadequate muscle fibre stretch at end-diastole, which depresses the ability of myocardium to generate adequate contraction force in the subsequent systole (according to the Frank-Starling law). Moreover, an excessively stiff infarct may cause its neighbouring myocardium to be understretched at end-diastole, subsequently depressing the systolic contractile force of the neighbouring myocardium, which was found to be correlated with infarct extension.

Impact of myocardial infarction on intraventricular vortex and flow energetics assessed using computational simulations.

Flow energetics have been proposed as early indicators of progressive left ventricular (LV) functional impairment in patients with myocardial infarction (MI), but its correlation with individual MI parameters has not been fully explored. Using electro-fluid-structure interaction LV models, this study investigated the correlation between four MI parameters: infarct size, infarct multiplicity, regional enhancement of contractility at the viable myocardium area (RECVM), and LV mechanical dyssynchrony (LVMD) with intraventricular vortex and flow energetics. In LV with small infarcts, our results showed that infarct appearance amplified the energy dissipation index (DI), where substantial viscous energy loss was observed in areas with high flow velocity and near the infarct-vortex interface. The LV with small multiple infarcts and RECVM showed remarkable DI increment during systole and diastole. In correlation analysis, the systolic kinetic energy fluctuation index (E') was positively related to ejection fraction (EF) (R2  = 0.982) but negatively correlated with diastolic E' (R2  = 0.970). Diastolic E' was inversely correlated with vortex kinetic energy (R2  = 0.960) and vortex depth (R2  = 0.876). We showed an excessive systolic DI could differentiate infarcted LV with normal EF from healthy LV. Strong flow acceleration, LVMD, and vortex-infarct interactions were predominant factors that induced excessive DI in infarcted LVs. Instead of causing undesired flow turbulence, high systolic E' suggested the existence of energetic flow acceleration, while high diastolic E' implied an inefficient diastolic filling. Thus, systolic E' is not a suitable early indicator for progressive LV dysfunction in MI patients, while diastolic E' may be a useful index to indicate diastolic impairment in these patients.

Assessing regional left ventricular thickening dysfunction and dyssynchrony via personalized modeling and 3D wall thickness measurements for acute myocardial infarction.

BACKGROUND: Existing clinical diagnostic and assessment methods could be improved to facilitate early detection and treatment of cardiac dysfunction associated with acute myocardial infarction (AMI) to reduce morbidity and mortality. PURPOSE: To develop 3D personalized left ventricular (LV) models and thickening assessment framework for assessing regional wall thickening dysfunction and dyssynchrony in AMI patients. STUDY TYPE: Retrospective study, diagnostic accuracy. SUBJECTS: Forty-four subjects consisting of 15 healthy subjects and 29 AMI patients. FIELD STRENGTH/SEQUENCE: 1.5T/steady-state free precession cine MRI scans; LGE MRI scans. ASSESSMENT: Quantitative thickening measurements across all cardiac phases were correlated and validated against clinical evaluation of infarct transmurality by an experienced cardiac radiologist based on the American Heart Association (AHA) 17-segment model. STATISTICAL TEST: Nonparametric 2-k related sample-based Kruskal-Wallis test; Mann-Whitney U-test; Pearson's correlation coefficient. RESULTS: Healthy LV wall segments undergo significant wall thickening (P < 0.05) during ejection and have on average a thicker wall (8.73 ± 1.01 mm) compared with infarcted wall segments (2.86 ± 1.11 mm). Myocardium with thick infarct (ie, >50% transmurality) underwent remarkable wall thinning during contraction (thickening index [TI] = 1.46 ± 0.26 mm) as opposed to healthy myocardium (TI = 4.01 ± 1.04 mm). For AMI patients, LV that showed signs of thinning were found to be associated with a significantly higher percentage of dyssynchrony as compared with healthy subjects (dyssynchrony index [DI] = 15.0 ± 5.0% vs. 7.5 ± 2.0%, P < 0.01). Also, a strong correlation was found between our TI and left ventricular ejection fraction (LVEF) (r = 0.892, P < 0.01), and moderate correlation between DI and LVEF (r = 0.494, P < 0.01). DATA CONCLUSION: The extracted regional wall thickening and DIs are shown to be strongly correlated with infarct severity, therefore suggestive of possible practical clinical utility. LEVEL OF EVIDENCE: 2 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;49:1006-1019.

Treatment of chronic back pain using indirect vibroacoustic therapy: A pilot study.

BACKGROUND: Low frequency sound wave stimulation therapy has become increasingly popular in the rehabilitation fields, due to its ease, less fatiguing and time efficient application. OBJECTIVE: This 12-week pilot study examines the efficacy of applying low frequency sound wave stimulation (between 16-160 Hz) through both hands and feet on relieving pain and improving functional ability in patients with chronic back pain. METHODS: Twenty-three participants with chronic shoulder (eleven participants) or low back pain (twelve participants) underwent a 12-week vibration therapy program of three sessions per week. A low frequency sound wave device comprising four piezoelectric vibration-type tactile tranducers enclosed in separate 5-cm diameter circular plates, which generate sinusoidal vibratory stimuli at a frequency of 16-160 Hz, was used in this study. Primary outcome measure was pain sensation measured using the Visual Analogue Scale (P-VAS). The secondary outcome measures were pain-related disability measured using the pain disability index (PDI) and quality of life measured using the SF-12. RESULTS: At week 12, significant reductions in pain sensation and pain-related disability were observed, with mean reductions of 3.5 points in P-VAS and 13.5 points in the PDI scores. Sixty-five percent of the participants had a reduction of at least 3 points on the P-VAS score, while 52% participants showed a decrease of at least 10 points in the PDI score. Significant improvement was observed in the SF-12 physical composite score but not the mental composite score. CONCLUSIONS: The preliminary findings showed that passive application of low frequency sound wave stimulation therapy through both hands and feet was effective in alleviating pain and improving functional ability in patients with chronic back pain.

Quantitative analysis of intraventricular flow-energetics and vortex in ischaemic hearts.

OBJECTIVE: This study investigated the intraventricular flow dynamics in ischaemic heart disease patients. PATIENTS AND METHODS: Fourteen patients with normal ejection fraction and 16 patients with reduced ejection fraction were compared with 20 healthy individuals. Phase-contrast MRI was used to assess intraventricular flow variables and speckle-tracking echocardiography to assess myocardial strain and left ventricular (LV) dyssynchrony. Infarct size was acquired using delayed-enhancement MRI. RESULTS: The results obtained showed no significant differences in intraventricular flow variables between the healthy group and the patients with normal ejection fraction group, whereas considerable reductions in kinetic energy (KE) fluctuation index, E' (P<0.001) and vortex KE (P=0.003) were found in the patients with reduced ejection fraction group. In multivariate analysis, only vortex KE and infarct size were significantly related to LV ejection fraction (P<0.001); furthermore, vortex KE was correlated negatively with energy dissipation, energy dissipation index (r=-0.44, P=0.021). CONCLUSION: This study highlights that flow energetic indices have limited applicability as early predictors of LV progressive dysfunction, whereas vortex KE could be an alternative to LV performance.

Left ventricular flow propagation velocity measurement: Is it cast in stone?

This study aims to investigate the measurement of left ventricular flow propagation velocity, V p, using phase contrast magnetic resonance imaging and to assess the discrepancies resulting from inflow jet direction and individual left ventricular size. Three V p measuring techniques, namely non-adaptive (NA), adaptive positions (AP) and adaptive vectors (AV) method, were suggested and compared. We performed the comparison on nine healthy volunteers and nine post-infarct patients at four measurement positions, respectively, at one-third, one-half, two-thirds and the conventional 4 cm distances from the mitral valve leaflet into the left ventricle. We found that the V p measurement was affected by both the inflow jet direction and measurement positions. Both NA and AP methods overestimated V p, especially in dilated left ventricles, while the AV method showed the strongest correlation with the isovolumic relaxation myocardial strain rate (r = 0.53, p < 0.05). Using the AV method, notable difference in mean V p was also observed between healthy volunteers and post-infarct patients at positions of: one-half (81 ± 31 vs. 58 ± 25 cm/s), two-thirds (89 ± 32 vs. 45 ± 15 cm/s) and 4 cm (98 ± 23 vs. 47 ± 13 cm/s) distances. The use of AV method and measurement position at one-half distance was found to be the most suitable method for assessing diastolic dysfunction given varying left ventricular sizes and inflow jet directions.

Effect of spatial inlet velocity profiles on the vortex formation pattern in a dilated left ventricle.

Despite the advancement of cardiac imaging technologies, these have traditionally been limited to global geometrical measurements. Computational fluid dynamics (CFD) has emerged as a reliable tool that provides flow field information and other variables essential for the assessment of the cardiac function. Extensive studies have shown that vortex formation and propagation during the filling phase acts as a promising indicator for the diagnosis of the cardiac health condition. Proper setting of the boundary conditions is crucial in a CFD study as they are important determinants, that affect the simulation results. In this article, the effect of different transmitral velocity profiles (parabolic and uniform profile) on the vortex formation patterns during diastole was studied in a ventricle with dilated cardiomyopathy (DCM). The resulting vortex evolution pattern using the uniform inlet velocity profile agreed with that reported in the literature, which revealed an increase in thrombus risk in a ventricle with DCM. However the application of a parabolic velocity profile at the inlet yields a deviated vortical flow pattern and overestimates the propagation velocity of the vortex ring towards the apex of the ventricle. This study highlighted that uniform inlet velocity profile should be applied in the study of the filling dynamics in a left ventricle because it produces results closer to that observed experimentally.

Sensitivity analysis of left ventricle with dilated cardiomyopathy in fluid structure simulation.

Dilated cardiomyopathy (DCM) is the most common myocardial disease. It not only leads to systolic dysfunction but also diastolic deficiency. We sought to investigate the effect of idiopathic and ischemic DCM on the intraventricular fluid dynamics and myocardial wall mechanics using a 2D axisymmetrical fluid structure interaction model. In addition, we also studied the individual effect of parameters related to DCM, i.e. peak E-wave velocity, end systolic volume, wall compliance and sphericity index on several important fluid dynamics and myocardial wall mechanics variables during ventricular filling. Intraventricular fluid dynamics and myocardial wall deformation are significantly impaired under DCM conditions, being demonstrated by low vortex intensity, low flow propagation velocity, low intraventricular pressure difference (IVPD) and strain rates, and high-end diastolic pressure and wall stress. Our sensitivity analysis results showed that flow propagation velocity substantially decreases with an increase in wall stiffness, and is relatively independent of preload at low-peak E-wave velocity. Early IVPD is mainly affected by the rate of change of the early filling velocity and end systolic volume which changes the ventriculo:annular ratio. Regional strain rate, on the other hand, is significantly correlated with regional stiffness, and therefore forms a useful indicator for myocardial regional ischemia. The sensitivity analysis results enhance our understanding of the mechanisms leading to clinically observable changes in patients with DCM.

Review on CFD simulation in heart with dilated cardiomyopathy and myocardial infarction.

The heart is a sophisticated functional organ that plays a crucial role in the blood circulatory system. Hemodynamics within the heart chamber can be indicative of exert cardiac health. Due to the limitations of current cardiac imaging modalities, computational fluid dynamics (CFD) have been widely used for the purposes of cardiac function assessment and heart disease diagnosis, as they provide detailed insights into the cardiac flow field. An understanding of ventricular hemodynamics and pathological severities can be gained through studies that employ the CFD method. In this research the hemodynamics of two common myocardial diseases, dilated cardiomyopathy (DCM) and myocardial infarction (MI) were investigated, during both the filling phase and the whole cardiac cycle, through a prescribed geometry and fluid structure interaction (FSI) approach. The results of the research indicated that early stage disease identification and the improvement of cardiac assisting devices and therapeutic procedures can be facilitated through the use of the CFD method.