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

Computational cardiac physiology for new modelers: Origins, foundations, and future

Mathematical models of the cardiovascular system have come a long way since they were first introduced in the early 19th century. Driven by a rapid development of experimental techniques, numerical methods, and computer hardware, detailed models that describe physical scales from the molecular level up to organs and organ systems have been derived and used for physiological research. Mathematical and computational models can be seen as condensed and quantitative formulations of extensive physiological knowledge and are used for formulating and testing hypotheses, interpreting and directing experimental research, and have contributed substantially to our understanding of cardiovascular physiology. However, in spite of the strengths of mathematics to precisely describe complex relationships and the obvious need for the mathematical and computational models to be informed by experimental data, there still exist considerable barriers between experimental and computational physiological research. In this review, we present a historical overview of the development of mathematical and computational models in cardiovascular physiology, including the current state of the art. We further argue why a tighter integration is needed between experimental and computational scientists in physiology, and point out important obstacles and challenges that must be overcome in order to fully realize the synergy of experimental and computational physiological research.

A Study of the Fluid–Structure Interaction of the Plaque Circumferential Distribution in the Left Coronary Artery

Atherosclerotic plaques within the coronary arteries can prevent blood from flowing to downstream tissues, causing coronary heart disease and a myocardial infarction over time. The degree of stenosis is an important reference point during percutaneous coronary intervention (PCI). However, clinically, patients with the same degree of stenosis exhibit different degrees of disease severity. To investigate the connection between this phenomenon and the plaque circumferential distribution, in this paper, four models with different plaque circumferential locations were made based on the CT data. The blood in the coronary arteries was simulated using the fluid–structure interaction method in ANSYS Workbench software. The results showed that the risk of plaque rupture was less affected by the circumferential distribution of plaque, and the distribution of blood in each branch was affected by the circumferential distribution of plaque. Low TAWSS areas were found posterior to the plaque, and the TAWSS < 0.4 Pa area was ranked from highest to lowest in each model species: plaque on the side away from the left circumflex branch, plaque on the side away from the heart; plaque on the side close to the heart; and plaque on the side close to the left circumflex branch. The same trend was also found in the OSI. It was concluded that the circumferential distribution of plaques affects their further development. This finding will be useful for clinical treatment.

A parallel non-nested two-level domain decomposition method for simulating blood flows in cerebral artery of stroke patient

Numerical simulation of blood flows in patient-specific arteries can be useful for the understanding of vascular diseases, as well as for surgery planning. In this paper, we simulate blood flows in the full cerebral artery of stroke patients. To accurately resolve the flow in this rather complex geometry with stenosis is challenging and it is also important to obtain the results in a short amount of computing time so that the simulation can be used in pre- and/or post-surgery planning. For this purpose, we introduce a highly scalable, parallel non-nested two-level domain decomposition method for the three-dimensional unsteady incompressible Navier-Stokes equations with an impedance outlet boundary condition. The problem is discretized with a stabilized finite element method on unstructured meshes in space and a fully implicit method in time, and the large nonlinear systems are solved by a preconditioned parallel Newton-Krylov method with a two-level Schwarz method. The key component of the method is a non-nested coarse problem solved using a subset of processor cores and its solution is interpolated to the fine space using radial basis functions. To validate and verify the proposed algorithm and its highly parallel implementation, we consider a case with available clinical data and show that the computed result matches with the measured data. Further numerical experiments indicate that the proposed method works well for realistic geometry and parameters of a full size cerebral artery of an adult stroke patient on a supercomputers with thousands of processor cores.

Exogenous spermine attenuates myocardial fibrosis in diabetic cardiomyopathy by inhibiting endoplasmic reticulum stress and the canonical Wnt signaling pathway

Myocardial fibrosis is one of the main pathological manifestations of diabetic cardiomyopathy (DCM). Spermine (SPM), a product of polyamine metabolism, plays an important role in many cardiac diseases including hypertrophy, ischemia, and infarction, but its role in diabetic myocardial fibrosis has not been clarified. This study aimed to investigate the role of polyamine metabolism, specifically SPM, in diabetic myocardial fibrosis and to explore the related mechanisms. We used intraperitoneal injection of streptozotocin (STZ, 60 mg/kg) in Wistar rats and high glucose (HG, 40 mM) stimulated cardiac fibroblasts (CFs) to established a type 1 diabetes (T1D) model in vivo and in vitro, which were pretreated with exogenous SPM (5 mg/kg per day and 5 μM). The results showed that hyperglycemia induced the expression of the key polyamine synthesis enzyme ornithine decarboxylase (ODC) decreased and the key catabolic enzyme spermidine/spermine N¹ -acetyltransferase (SSAT) increased compared with those in the control group. The body weight, blood insulin level, and cardiac ejection function were decreased, while blood glucose, heart weight, the ratio of heart weight to body weight, myocardial interstitial collagen deposition, and endoplasmic reticulum stress (ERS)-related protein (glucose-regulated protein-78, glucose-regulated protein-94, activating transcription factor-4, and C/EBP homology protein) expression in the T1D group were all significantly increased. HG also caused an increased expression of Wnt3, β-catenin (in cytoplasm and nucleus), while Axin2 and phosphorylated β-catenin decreased. Exogenous SPM improved the above changes caused by polyamine metabolic disorders. In conclusion, polyamine metabolism disorder occurs in the myocardial tissue of diabetic rats, causing myocardial fibrosis and ERS. Exogenous SPM plays a myocardial protective role via inhibiting of ERS and the canonical Wnt/β-catenin signaling pathway.

A review study on blood in human coronary artery: Numerical approach

Computational fluid dynamics (CFD) study of blood flow in human coronary artery is one of the emerging fields of Biomed- ical engineering. In present review paper, Finite Volume Method with governing equations and boundary conditions are briefly discussed for different coronary models. Many researchers have come up with astonishing results related to the various factors (blood viscosity, rate of blood flow, shear stress on the arterial wall, Reynolds number, etc.) affecting the hemodynamic of blood in the right/left coronary artery. The aim of this paper is to present an overview of all those work done by the researchers to justify their work related to factors which hampers proper functioning of heart and lead to Coronary Artery Disease (CAD). Governing equations like Navier-stokes equations, continuity equations etc. are widely used and are solved using CFD solver to get a clearer view of coronary artery blockage. Different boundary conditions and blood properties published in the last ten years are summarized in the tabulated form. This table will help new researchers to work on this area.

Human fetal hearts with tetralogy of Fallot have altered fluid dynamics and forces

Studies have suggested the effect of blood flow forces in pathogenesis and progression of some congenital heart malformations. It is therefore of interest to study the fluid mechanic environment of the malformed prenatal heart, such as the tetralogy of Fallot (TOF), especially when little is known about fetal TOF. In this study, we performed patient-specific ultrasound-based flow simulations of three TOF and seven normal human fetal hearts. TOF right ventricles (RVs) had smaller end-diastolic volumes (EDVs) but similar stroke volumes (SVs), whereas TOF left ventricles (LVs) had similar EDVs but slightly increased SVs compared with normal ventricles. Simulations showed that TOF ventricles had elevated systolic intraventricular pressure gradient (IVPG) and required additional energy for ejection but IVPG elevations were considered to be mild relative to arterial pressure. TOF RVs and LVs had similar pressures because of equalization via ventricular septal defect (VSD). Furthermore, relative to normal, TOF RVs had increased diastolic wall shear stresses (WSS) but TOF LVs were not. This was caused by high tricuspid inflow that exceeded RV SV, leading to right-to-left shunting and chaotic flow with enhanced vorticity interaction with the wall to elevate WSS. Two of the three TOF RVs but none of the LVs had increased thickness. As pressure elevations were mild, we hypothesized that pressure and WSS elevation could play a role in the RV thickening, among other causative factors. Finally, the endocardium surrounding the VSD consistently experienced high WSS because of RV-to-LV flow shunt and high flow rate through the over-riding aorta. NEW & NOTEWORTHY Blood flow forces are thought to cause congenital heart malformations and influence disease progression. We performed novel investigations of intracardiac fluid mechanics of tetralogy of Fallot (TOF) human fetal hearts and found essential differences from normal hearts. The TOF right ventricle (RV) and left ventricle had similar and elevated pressure but only the TOF RV had elevated wall shear stress because of elevated tricuspid inflow, and this may contribute to the observed RV thickening. TOF hearts also expended more energy for ejection.

High-speed, high-frequency ultrasound, in utero vector-flow imaging of mouse embryos

Real-time imaging of the embryonic murine cardiovascular system is challenging due to the small size of the mouse embryo and rapid heart rate. High-frequency, linear-array ultrasound systems designed for small-animal imaging provide high-frame-rate and Doppler modes but are limited in regards to the field of view that can be imaged at fine-temporal and -spatial resolution. Here, a plane-wave imaging method was used to obtain high-speed image data from in utero mouse embryos and multi-angle, vector-flow algorithms were applied to the data to provide information on blood flow patterns in major organs. An 18-MHz linear array was used to acquire plane-wave data at absolute frame rates ≥10 kHz using a set of fixed transmission angles. After beamforming, vector-flow processing and image compounding, effective frame rates were on the order of 2 kHz. Data were acquired from the embryonic liver, heart and umbilical cord. Vector-flow results clearly revealed the complex nature of blood-flow patterns in the embryo with fine-temporal and -spatial resolution.

Comprehensive computational assessment of blood flow characteristics of left ventricle based on in-vivo MRI in presence of artificial myocardial infarction

BACKGROUND

Understanding the effects of cardiac diseases on the heart's functionality which is the purpose of many biomedical researches, directly affects the diagnostic and therapeutic methods. Myocardial infarction (MI) is a common complication of cardiac ischemia, however, the impact of MI on the left ventricle (LV) flow patterns has not been widely considered by computational fluid dynamics studies thus far.

METHODS

In this study, we present an insightful numerical method that creates an artificial MI on an image-based fluid-structure interactional model of normal LV to investigate its influence on the flow in comparison with the normal case. Seventeen different models were developed to evaluate the effects of location, percentage, myocardial material properties and dilation size of MI on the LV's performance, area strain, wall displacement, pressure-volume loop, wall shear stress and velocity field.

RESULTS

The results show that MI considerably changes blood flow features which are fully dependent on MI parameters. For the case of constant MI location, the effect of a decrease of infarcted myocardium stiffness, increase of dilation size and increase of MI percentage are mostly similar. Although the location differences of MI under other constant conditions have similar impact on the ejection fraction, they also lead to dissimilar variations in the LV flow pattern and other indicators.

CONCLUSIONS

The presented model showed a capable computational method for investigating various mechanical MI conditions with respect to cardiac flow pattern. The perspective of this model development seems to be an applicable tool for MI clinical diagnosis and prediction of complications related to MI.

Pulsatile hemodynamics in patient-specific thoracic aortic dissection models constructed from computed tomography angiography

PURPOSE

Thoracic aortic dissection (TAD) is considered one of the most catastrophic and non-traumatic cardiovascular diseases associated with high morbidity and mortality rates in clinical treatment. The purpose of this paper is to investigate the pulsatile hemodynamics changes throughout a cardiac cycle in a Stanford Type B TAD model with the aid of computational fluid dynamics (CFD) method.

METHODS

A patient-specific dissected aorta geometry was reconstructed from the three-dimensional (3D) computed tomography angiography (CTA) scanning. The realistic time-dependent pulsatile boundary conditions were prescribed for our 3D patient-specific TAD model. Blood was considered to be an incompressible, Newtonian fluid. The aortic wall was assumed to be rigid, and a no-slip boundary condition was applied at the wall. CFD simulations were processed using the finite volume (FV) method to investigate the pulsatile hemodynamics in terms of blood flow velocity, aortic wall pressure, wall shear stress and flow vorticity. In the experiments, blood velocity, pressure, wall shear stress and vorticity distributions were analyzed qualitatively and quantitatively.

RESULTS

The experimental results demonstrated a high wall shear stress and strong vertical flow at dissection initiation. The results also indicated that wall shear progressed along the false lumen, which is a possible cause of blood flow between aortic wall layers.

Heart blood flow simulation: a perspective review

Cardiovascular disease (CVD), the leading cause of death today, incorporates a wide range of cardiovascular system malfunctions that affect heart functionality. It is believed that the hemodynamic loads exerted on the cardiovascular system, the left ventricle (LV) in particular, are the leading cause of CVD initiation and propagation. Moreover, it is believed that the diagnosis and prognosis of CVD at an early stage could reduce its high mortality and morbidity rate. Therefore, a set of robust clinical cardiovascular assessment tools has been introduced to compute the cardiovascular hemodynamics in order to provide useful insights to physicians to recognize indicators leading to CVD and also to aid the diagnosis of CVD. Recently, a combination of computational fluid dynamics (CFD) and different medical imaging tools, image-based CFD (IB-CFD), has been widely employed for cardiovascular functional assessment by providing reliable hemodynamic parameters. Even though the capability of CFD to provide reliable flow dynamics in general fluid mechanics problems has been widely demonstrated for many years, up to now, the clinical implications of the IB-CFD patient-specific LVs have not been applicable due to its limitations and complications. In this paper, we review investigations conducted to numerically simulate patient-specific human LV over the past 15 years using IB-CFD methods. Firstly, we divide different studies according to the different LV types (physiological and different pathological conditions) that have been chosen to reconstruct the geometry, and then discuss their contributions, methodologies, limitations, and findings. In this regard, we have studied CFD simulations of intraventricular flows and related cardiology insights, for (i) Physiological patient-specific LV models, (ii) Pathological heart patient-specific models, including myocardial infarction, dilated cardiomyopathy, hypertrophic cardiomyopathy and hypoplastic left heart syndrome. Finally, we discuss the current stage of the IB-CFD LV simulations in order to mimic realistic hemodynamics of patient-specific LVs. We can conclude that heart flow simulation is on the right track for developing into a useful clinical tool for heart function assessment, by (i) incorporating most of heart structures' (such as heart valves) operations, and (ii) providing useful diagnostic indices based hemodynamic parameters, for routine adoption in clinical usage.

Multiphysics and multiscale modelling, data-model fusion and integration of organ physiology in the clinic: ventricular cardiac mechanics

With heart and cardiovascular diseases continually challenging healthcare systems worldwide, translating basic research on cardiac (patho)physiology into clinical care is essential. Exacerbating this already extensive challenge is the complexity of the heart, relying on its hierarchical structure and function to maintain cardiovascular flow. Computational modelling has been proposed and actively pursued as a tool for accelerating research and translation. Allowing exploration of the relationships between physics, multiscale mechanisms and function, computational modelling provides a platform for improving our understanding of the heart. Further integration of experimental and clinical data through data assimilation and parameter estimation techniques is bringing computational models closer to use in routine clinical practice. This article reviews developments in computational cardiac modelling and how their integration with medical imaging data is providing new pathways for translational cardiac modelling.

Computational fluid dynamics modelling in cardiovascular medicine

This paper reviews the methods, benefits and challenges associated with the adoption and translation of computational fluid dynamics (CFD) modelling within cardiovascular medicine. CFD, a specialist area of mathematics and a branch of fluid mechanics, is used routinely in a diverse range of safety-critical engineering systems, which increasingly is being applied to the cardiovascular system. By facilitating rapid, economical, low-risk prototyping, CFD modelling has already revolutionised research and development of devices such as stents, valve prostheses, and ventricular assist devices. Combined with cardiovascular imaging, CFD simulation enables detailed characterisation of complex physiological pressure and flow fields and the computation of metrics which cannot be directly measured, for example, wall shear stress. CFD models are now being translated into clinical tools for physicians to use across the spectrum of coronary, valvular, congenital, myocardial and peripheral vascular diseases. CFD modelling is apposite for minimally-invasive patient assessment. Patient-specific (incorporating data unique to the individual) and multi-scale (combining models of different length- and time-scales) modelling enables individualised risk prediction and virtual treatment planning. This represents a significant departure from traditional dependence upon registry-based, population-averaged data. Model integration is progressively moving towards ‘digital patient’ or ‘virtual physiological human’ representations. When combined with population-scale numerical models, these models have the potential to reduce the cost, time and risk associated with clinical trials. The adoption of CFD modelling signals a new era in cardiovascular medicine. While potentially highly beneficial, a number of academic and commercial groups are addressing the associated methodological, regulatory, education- and service-related challenges.

4D flow cardiovascular magnetic resonance consensus statement

Pulsatile blood flow through the cavities of the heart and great vessels is time-varying and multidirectional. Access to all regions, phases and directions of cardiovascular flows has formerly been limited. Four-dimensional (4D) flow cardiovascular magnetic resonance (CMR) has enabled more comprehensive access to such flows, with typical spatial resolution of 1.5×1.5×1.5 - 3×3×3 mm(3), typical temporal resolution of 30-40 ms, and acquisition times in the order of 5 to 25 min. This consensus paper is the work of physicists, physicians and biomedical engineers, active in the development and implementation of 4D Flow CMR, who have repeatedly met to share experience and ideas. The paper aims to assist understanding of acquisition and analysis methods, and their potential clinical applications with a focus on the heart and greater vessels. We describe that 4D Flow CMR can be clinically advantageous because placement of a single acquisition volume is straightforward and enables flow through any plane across it to be calculated retrospectively and with good accuracy. We also specify research and development goals that have yet to be satisfactorily achieved. Derived flow parameters, generally needing further development or validation for clinical use, include measurements of wall shear stress, pressure difference, turbulent kinetic energy, and intracardiac flow components. The dependence of measurement accuracy on acquisition parameters is considered, as are the uses of different visualization strategies for appropriate representation of time-varying multidirectional flow fields. Finally, we offer suggestions for more consistent, user-friendly implementation of 4D Flow CMR acquisition and data handling with a view to multicenter studies and more widespread adoption of the approach in routine clinical investigations.