Liang Zhong
Jun-Mei Zhang
Boyang Su
Ru San Tan
John C. Allen
Ghassan S. Kassab- 1National Heart Centre Singapore, National Heart Research Institute of Singapore, Singapore, Singapore
- 2Duke-NUS Medical School, Singapore, Singapore
- 3California Medical Innovations Institute, San Diego, CA, United States
The emergence of new cardiac diagnostics and therapeutics of the heart has given rise to the challenging field of virtual design and testing of technologies in a patient-specific environment. Given the recent advances in medical imaging, computational power and mathematical algorithms, patient-specific cardiac models can be produced from cardiac images faster, and more efficiently than ever before. The emergence of patient-specific computational fluid dynamics (CFD) has paved the way for the new field of computer-aided diagnostics. This article provides a review of CFD methods, challenges and opportunities in coronary and intra-cardiac flow simulations. It includes a review of market products and clinical trials. Key components of patient-specific CFD are covered briefly which include image segmentation, geometry reconstruction, mesh generation, fluid-structure interaction, and solver techniques.
Introduction—Coronary Artery Disease, Cardiac Dysfunction, and Diagnosis
In coronary artery disease (CAD) atherosclerotic build-up can narrow the arterial lumen, resulting in myocardial ischemia. Prevalence of CAD is 6% in the general population and up to 20% in those aged over 65 years. About 13% of deaths are due to CAD. By 2030, it is projected that 15% of male deaths will be attributable to CAD (World Health Organization, 2011).
CAD can be diagnosed by means of either an anatomic parameter, such as diameter stenosis or a functional parameter linked to coronary territory myocardial ischemia. Stenosis does not invariably impair distal coronary flow, and this is particularly true with regard to the intermediate coronary artery lesions (i.e., diameter stenosis between 30 and 70%). Non-invasive tests of myocardial ischemia (e.g., nuclear myocardial perfusion imaging, stress echocardiography) identify areas of the most severely reduced relative coronary flow reserve. They are fairly accurate for myocardial ischemia detection on a per-patient basis, but these perform less well in quantifying severity of individual coronary territory ischemia. The latter is relevant in multi-vessel percutaneous coronary intervention (PCI), where coronary physiological information, overlaid on detailed maps of patient-specific coronary artery anatomy, dictates management decisions. Fractional flow reserve (FFR) measured during invasive coronary angiography (ICA) under adenosine-induced hyperemia has emerged as the gold standard for assessment of coronary flow physiology and coronary territory ischemia (Johnson et al., 2012).
Diagnosis of heart contractile dysfunction requires demonstration of either diastolic or systolic function abnormalities. The gold standard for determining diastolic dysfunction is an increase in invasively measured ventricular end-diastolic pressure–>15 mmHg in the case of the left ventricle (LV) (Nishimura and Tajik, 1997). Systolic dysfunction is assessed by the change in maximal ventricular pressure (P) during isovolumic contraction, dP/dtmax (Yamada et al., 1998). Multiple ventricular pressure-volume loops assayed using a conductance catheter under varying loading conditions can yield end-diastolic (Eed) and end-systolic elastances (Ees) that characterize ventricular diastolic and systolic dysfunction, respectively (Burkhoff et al., 2005). Emerging noninvasive echocardiographic and cardiac magnetic resonance (MRI) imaging techniques enable corroborative assessment of regional and global cardiac chamber dysfunction involving strain and strain rate (Zhong et al., 2012), curvedness (a descriptor of three-dimensional ventricular shape) and curvedness rate (Zhong et al., 2009), ventricular contractility dσ*/dtmax, where σ* is pressure-normalized wall stress (Zhong et al., 2014), and atrio-ventricular velocities (Leng et al., 2016).
Patient-specific computational fluid dynamics (CFD) modeling is a recent development. Non-invasive FFR (FFRCT) is derived from CFD modeling of images acquired using computed tomography coronary angiography (CTCA). With invasive FFR as the gold standard, FFRCT ≤ 0.80 is superior to both CTCA and ICA determined diameter stenosis for ascertaining ischemia on a coronary artery territory basis (Min et al., 2012). FFRCT analysis is solely available via a centralized commercial web-based service of the HeartFlow® company. Time-consuming computational demands and high costs−6 h and $2000 USD to process a case (Kimura et al., 2015)—hamper widespread clinical adoption. The requisite offsite handling of sensitive confidential patient information and associated medical conditions is a highly delicate issue involving IT-security, potential for data abuse, etc.
Unlike coronary blood flow simulation, CFD studies on intra-cardiac flows are primarily confined to research purposes owing to the complexity of modeling intra-cardiac flows. In truth, coronary and intra-cardiac flows are closely connected. Coronary artery dysfunction leads to myocardium ischemia, and intra-cardiac flows provide blood for circulation throughout the body, including the coronary circulation. Prolonged and untreated myocardial ischemia could increase the risk for death or myocardial infarction (Iskander and Iskandrian, 1998). Future integration of both coronary and intra-cardiac flow simulations is desirable to enable a comprehensive assessment of cardiac circulatory pathophysiology. This paper aims to pave the way for integrated simulations and focuses on a progress review of CFD applications in modeling coronary and intra-cardiac flow. Other applications of CFD in cardiovascular disease can be found in Morris et al. (2015).
Coronary Flow Simulation
Challenges and Opportunities in Patient-Specific Simulation Techniques for Studying Blood Flow in Coronary Arteries
In general, the tasks involved in performing CFD simulation for a patient-specific coronary artery tree are as follows: (1) Image acquisition and segmentation to reconstruct a 3D patient-specific coronary model; (2) CFD preprocessing to discretize the domain with meshes and define the boundary conditions; and (3) Solving the fluid governing equations using a fluid solver and post-processing to visualize the flow field (Figure 1). If fluid structure interaction (FSI) is considered, an additional solid solver is used, and coupling between fluid and solid solvers is implemented. Table 1 summarizes the challenges involved.
Figure 1. Schematic drawings for the procedures of patient-specific simulations for blood flow in coronary arteries, including (1) acquisition of CTA or ICA images (2) segmentation of acquired images (3) reconstruction of 3D model (4) mesh generation to discretize the 3D model (5) solving the mass and momentum conservation equations to simulate the blood flow in coronary arteries, if FSI is taken into consideration, solid solver will also be activated, and (6) presentation of simulations.
Table 1. Current challenges and opportunities in numerical simulation of coronary arteries.
Coronary anatomy can be imaged using ICA, intravascular ultrasound (IVUS), optical coherence tomography (OCT), CTCA, and MRI (Zhang et al., 2014). Invasive IVUS and OCT yields high-resolution cross-sectional views of the coronary arteries, and can be used in conjunction with biplane ICA to reconstruct the 3D vessel model. Since non-invasive CTCA possesses higher spatial resolution than MRI and echocardiography, it is widely used for 3D patient-specific coronary model reconstruction. However present CTCA has a spatial resolution of about 0.3 mm, which limits its use to coronary arteries of 1 mm or greater in diameter. Although CAD is not generally characterized nor is FFR measured in such small vessels, the latter is essential for characterizing coronary microcirculation.
Sophisticated segmentation approaches such as level-set segmentation (Bekkers and Taylor, 2008) have been applied to reconstruct 2D and 3D patient-specific coronary models, either by fusion of biplane ICA with IVUS images (Papafaklis et al., 2007) or directly from CTCA images (Torii et al., 2009). Commercial (e.g., 3D Doctor, Mimics, SliceOmatic, Amira) and open-source general image processing tools [e.g., VTK, ITK, ITK-SNAP, VTK, Analyze, and ImageJ (or Fiji)], make reconstruction of patient-specific models from medical images possible. Furthermore, their plug-in capability allows easy customization of segmentation tools. Artifacts such as calcification, motion and mis-registration are not easily overcome by segmentation techniques, and remain challenging.
In terms of simulation tools, commercial software such as ANSYS (including ICEM, FLUENT, CFX), STAR-CCM, and open-source tools (e.g., OpenFOAM) is applicable to general CFD simulations, including simulating the blood flow in coronary arteries. SimVascular (Schmidt et al., 2008) is a special tool designed for simulating the blood flow in vessels. These tools allow users specifications on mesh generation, boundary conditions settings and etc. As regards meshes, mesh generation schemes can be classified as structured or unstructured meshes. Structured grid generators, including “block-structured” techniques (used in ICEM CFD, TrueGrid, and IA-FEMesh) generally require complex iterative smoothing procedures to align elements with boundaries or physical domain. For the complex 3D coronary artery models reconstructed from medical images, unstructured meshes are commonly needed, which are built based on node coordinates and the connections between nodes to form elements. Commercial packages (e.g., ANSYS, TGrid) and open-source (e.g., TetGen, gmesh) allow automatic discretization of complex geometry with tetrahedral meshes. However quality control of tetrahedral meshes can be challenging and varies according to the mesh generation method employed (Wittek et al., 2016). The advancing front method, such as the Delaunay triangulation method, can provide better control of the mesh quality, but at the expense of prolonged computational time. In addition, 4-noded tetrahedral elements are involved with artificial stiffening, which presents challenges in modeling soft tissue, such as the coronary artery wall (Wittek et al., 2010). Higher-order and mixed-formulation tetrahedral elements can assist in overcoming these challenges. Nevertheless their computational cost is about four times higher than the 4-noded tetrahedral elements (Bourdin et al., 2007). To overcome the difficulty of generating good quality meshes for complex geometry with limited time and the convergence difficulties in modeling structures with large deformations, meshless methods have been recognized as one possible solution (Doblare et al., 2005). Meshless methods discretize the computational domain into a cloud of nodal points (Belinha, 2016). This discretization flexibility may allow direct model generation from CTCA or MRI images. However, meshless methods also have substantial shortcomings: (i) limited strict mathematical proof, (ii) incomplete theory, and (iii) lower computation efficiency compared to traditional computational methods using meshes (Zhang et al., 2012).
To achieve an acceptable simulation of blood flow in the coronary arteries, proper boundary conditions are paramount. Although flow and pressure waveforms can be obtained from the literature, in-vitro and in-vivo measurements are necessary for accurate simulations. For many CFD applications, it is virtually impossible to know flow and/or pressure waveforms a priori due to the difficulty of obtaining simultaneous measurements in coronary arteries. To solve this problem, multi-scale simulations have been developed that couple 3D simulation with reduced-order (1 or 0-dimensional) models at the boundaries. These models characterize pressure and flow rate in upstream and downstream vasculatures (Kim et al., 2010) as resistance, compliance, and impedance. How to determine the values of these patient-specific parameters remains a dilemma. Although morphologic information (e.g., scaling law) is widely used, numerical optimization may be necessary to tune these patient-specific parameters (Spilker et al., 2007).
Another obstacle in multi-scale simulation is how to solve the flow-governing equations using reduced-order models as boundary conditions. Both robust implicit approaches and explicit methods have been used. Explicit methods do not require changing the numerical algorithms to solve the governing equations (Sankaran et al., 2012), although they are generally less robust compared to the implicit fully coupled methods (Kim et al., 2010).
Recent progress in FSI techniques has allowed treating the coronary vessels as compliant. In the traditional Arbitrary Lagrangian-Eulerian (ALE) method (Malvè et al., 2012), boundaries and interfaces of both fluid and structural computational domains are precisely tracked during the iterations. When taking heart movement into consideration, re-meshing computational domains is often necessary to maintain mesh quality, which substantially increases computational cost. Over the years, stability of the ALE method has been improved (Hecht and Pironneau, 2017).
Alternative FSI techniques used to simulate flow in the presence of a moving boundary include immersed boundary and coupled momentum methods. These use fixed fluid meshes with boundaries defined by a set of moving Lagrangian points (Peskin, 2003) or linear membrane (Figueroa et al., 2006). Although prescribed heart motion has been used for simulating the blood flow in left (Prosi et al., 2004) and right (Torii et al., 2010) coronary arteries, more efficient, and robust FSI techniques are necessary to model large deformations of coronary arteries when taking heart motion during a cardiac cycle into account.
Recent development of fluid-solid-growth modeling incorporates vascular wall geometry, structure and properties governed by stress-mediated growth and remodeling (G&R) into FSI simulation (Figueroa et al., 2009). In this way, biofluid mechanics, biosolid mechanics, and biotransport phenomena, such as arterial growth, remodeling, and adaptation (Valentín et al., 2013) can be better understood. Recently, machine learning has been adopted into the CFD simulation to reduce the computational time incorporating features trained by a database generated from a set of offline CFD simulations (Sankaran et al., 2015; Itu et al., 2016). Sankaran and Marsden (2011) developed an adaptive collocation algorithm to quantify the effect of input uncertainties on cardiovascular simulation. Using a data-driven framework, Sankaran et al. (2016) studied the impact of anatomic and physiologic uncertainty (e.g., various boundary conditions and blood viscosity) on blood flow simulation. These data-driven modeling approaches wherein CFD data is used to enrich and refine the models may become popular in the future.
In addition, the assumption of blood as a Newtonian fluid is only valid for shear rate higher than 100 s−1, which may not be true when a flow recirculation region is formed near a coronary stenosis. The influence of non-Newtonian properties of blood on the velocity distribution and shear-thinning has been studied via various single and multi-phase non-Newtonian hemodynamic models (Jung et al., 2006). An effect of non-Newtonian properties on overall pressure drop across the arterial stenosis was exhibited at a flow with the Reynolds number of 100 or less (Cho and Kensey, 1991).
Challenges and Opportunities in Applying Patient-Specific CFD Technologies for Diagnosing the Severity of Coronary Artery Disease
Recent patient-specific CFD simulations provide detailed hemodynamic parameter (HP) information, such as pressure (P), wall pressure gradient (WPG) (Liu et al., 2012), wall shear stress (WSS) (Papafaklis et al., 2007; Stone et al., 2012), oscillatory shear index (OSI), relative residence time (RRT), and stress phase angle (SPA) (Knight et al., 2010), enabling characterization of HP distributions on coronary vessel walls. Definitions of these HPs (Table 2) and illustrations of their distributions in a left coronary artery tree (Figure 2) are presented. In the simulations, pressure and resistance boundary conditions were assigned to the inlet and outlets, respectively. A non-slip condition was imposed on the wall (Zhang et al., 2015).
Table 2. Hemodynamic parameters (HPs) predicted by CFD (Computational Fluid Dynamics) to link with CAD (Coronary Artery Disease).
Figure 2. Distributions of (A) P (Pressure), (B) WPG (wall pressure gradient), (C) WSS (wall shear stress), (D) OSI (oscillatory shear index), (E) RRT (relative residence time), and (F) SPA (stress phase angle) on the virtually healthy and diseased left coronary artery trees respectively. Labels of (A–D) indicate the stenosis locations. In the virtually healthy artery model, low WSS, and high RRT was exhibited in three of the four locations, where the stenoses were formed, and high WSS with low RRT was exhibited in the fourth. These findings suggest that coronary plaque is more likely to form in locations with low- WSS- and- high- RRT or high- WSS- and- low- RRT. From Zhang et al. (2015).
Among the HPs, pressure was closely related to FFR (Johnson et al., 2012), where FFR is the ratio of pressure distal to stenosis and aortic pressure at hyperemia, obtained during ICA. Alternative approaches for calculating FFR non-invasively have been attempted by several groups.
The Heart Flow company is a pioneer to combine CTCA images with CFD for calculating non-invasive FFRCT in CAD diagnosis. Kim et al (Kim et al., 2010) pioneered non-invasive FFRCT by reconstructing 3D patient-specific coronary artery model from CTCA images and coupling lumped parameter models to an implicit solver of fluid-governing equations (Figure 3A). Multicenter clinical trials DISCOVER-FLOW, DeFACTO, and NXT (Koo et al., 2011; Min et al., 2012; Nørgaard et al., 2014) demonstrated superior diagnostic accuracy for FFRCT vs. CTCA alone. PLATFORM trial (Douglas et al., 2015) (Lu et al., 2017) demonstrated the feasibility and safety of FFRCT as a diagnostic strategy in triage of patients with suspected CAD compared to standard of care. REAL-FFRCT (Kawaji et al., 2017) demonstrated good diagnostic performance of FFRCT even in patients with severely calcified vessels. Recently HeartFlow's FFRCT software has been approved by FDA for measuring coronary blockages non-invasively. However it used cloud technology for uploading CTCA images and downloading FFRCT report. This could involve data security issues that would hinder the on-site computation.
Figure 3. Diagrams and characteristics for calculating non-invasive FFR through combining CTCA with CFD by the companies of (A) Heart Flow: FFRCT (Gaur et al., 2013); (B) Siemens (1st generation: cFFR) (Coenen et al., 2015); (C) Siemens (2nd generation: cFFRML) (Itu et al., 2016); and (D) Toshiba: CT-FFR (Ko et al., 2017).
To facilitate on-site FFR computation supported by Siemens Company, Coenen et al. modeled coronary vessel as 1D segment and employed a reduced-order model for simulating the coronary circulation. In this way, the computational time was reduced to 5–10 min per patient (Coenen et al., 2015). Calculated pressure information, viz. computational FFR (cFFR), was mapped onto a 3D model reconstructed from CTCA images (Figure 3B). The correlation between cFFR and FFR, however, was poor (r = 0.59). To further reduce the computational time, Itu and colleagues from Siemens applied machine-learning approach for computing cFFRML with features extracted from training database (Itu et al., 2016). The database consisted of synthetically generated coronary artery models and corresponding FFR values computed from the CFD algorithm (Figure 3C).
Another attempt to reduce the computation time to be < 30 min was conducted by Ko and colleagues from Toshiba Medical Systems Corp. core laboratory (Ko et al., 2017). Differing from the above studies, 4 CTCA images were acquired and reconstructed at phases of 70, 80, 90, and 99% of the R-R interval. The arterial luminal deformation was taken into consideration and a reduced-order fluid model was used to simulate a 1-dimensional pressure and flow distribution in coronary tree (Figure 3D). In this approach, the interaction between fluid and structure is taken into account to some extent, although it is not a fully-coupled FSI simulation.
With the exception of CTCA images, attempts have been made to derive FFRQCA using CFD simulation in patient-specific coronary artery models reconstructed from ICA In a study involving 77 coronary vessels (Tu et al., 2014), FFRQCA correlated well with the gold standard FFR (r = 0.81, p < 0.001). Invasive FFRQCA obviates the need for pressure wire/catheter and adenosine. QFR was further derived from 3 flow models using fixed empiric flow velocity; modeled hyperemic flow velocity derived from measured angiography without administration of adenosine, and measured hyperemic flow velocity, respectively (Tu et al., 2016). Diagnostic accuracy of QFR was tested in the FAVOR II China (Xu et al., 2017) and WIFI II (Westra et al., 2018) studies. Based on this method, Medis QAngio 3D XA software was developed to calculate QFR.
Combining coronary angiogram images with CFD simulation was also studied by Morris et al. (2013) to estimate virtual fractional flow reserve (vFFR) with generic boundary conditions. VIRTUheart software was therefore developed to facilitate the calculation of vFFR. CathWorks is another tool available for FFR simulation-based service through the combination of coronary angiograms and CFD simulation.
Infusion OCT with a coronary angiogram was used by Poon et al. (2015, 2017) in an attempt to reconstruct the vessel and calculate virtually derived FFR.
Other HPs (Table 2) have been considered in relation to CAD based on biomechanical forces involved in regulation of blood vessel structure (Langille and O'Donnell, 1986). Among them, WSS has been the most studied. WSS is the frictional force of blood exerted tangentially to the luminal surface of the blood vessel per unit area. WSS is typically within the range of 15–20 dynes/cm2 (Kassab and Navia, 2006) for normal arteries; abnormal WSS outside this range promotes atherogenesis. There are two competing theories. In in-vitro experiments on canine thoracic aorta endothelial cells (ECs), ECs became abnormal for WSS >379 dynes/cm2 (Fry, 1968). This implied that high WSS might injure and denudate the ECs resulting in atherogenesis. Conversely, Caro et al. found intimal thickening and atherosclerosis with WSS < 6 dynes/cm2 (Caro et al., 1969). They conjectured that low WSS was associated with prolonged particle retention time, and increased intimal accumulation of lipids, leading to atherogenesis. Indeed, Rutsch et al. proposed several signaling pathways relating disturbed WSS with EC mechanico-chemical transductions (Rutsch et al., 2011).
CFD has been applied in the study of HP distributions, particularly WSS, in both idealized and patient-specific coronary artery models (Papafaklis et al., 2007). The PREDICTION study (Stone et al., 2012) enrolled 506 PCI patients and studied the natural history of plaque development in a subset of 374 (74%) over a 6- to 10-month period post-PCI. WSS was calculated for 3D coronary artery models reconstructed by combining intracoronary IVUS and biplane ICA images. Large plaque burden and low local WSS were found to be independent and additive predictors of plaque progression and luminal narrowing. Other clinical trials are currently underway investigating potential links between HPs and coronary atherogenesis (Antoniadis et al., 2015). These presage exciting possibilities for identifying indicators that may be useful for CAD diagnosis and patient management.
Intra-Cardiac Flow Simulation
Challenges and Opportunities in Patient-Specific Simulations for Studying Intra-Cardiac Flow
Figure 4 illustrates the general process of CFD simulation of blood flow in the LV based on cardiac MRI images (Su et al., 2015). Images comprise the long-axis and a stack of 12–14 short-axis images covering the LV from apex to base. Typical slice thickness is 8 mm; and typical frame rate 20–40 per cardiac cycle. Short-axis images are segmented either manually or automatically wherein blood pool is distinguished from heart muscle, and papillary muscles are included in blood pool. Long-axis images are used to track the mitral annulus at the intersection of the LV and the left atrium, which cannot be segmented easily on short-axis images due to through-plane displacement. A similar method is applied to construct the aortic annulus at the LV outflow tract. A patient-specific model based on segmentations is thereby generated (steps are highlighted by red rectangle in 1.1 of Figure 4).
Figure 4. Flow chart of CFD simulation of patient-specific intra-cardiac flow. Black box: each step of numerical simulation (Su et al., 2016). Red box: left ventricle reconstruction from MRI images. Green box: CFD mesh generation. Violet box: Mapping between reconstructed model and CFD surface mesh. Dotted blue box: CFD mesh resulted from mapping. Blue box: A series of CFD meshes at each frame. Elements (blue) and grids (green) are for reconstructed geometry and CFD mesh, respectively.
Next, the complex 3D models are discretized into tetrahedral, hexahedral or polyhedral grids using either in-house or commercial mesh generators such as ANSYS ICEM CFD, STAR-CCM+. A tetrahedral mesh is frequently adopted, which requires re-meshing when spring-based smoothing fails to cope with large deformations. Polyhedral meshes confer superior convergence speed (Spiegel et al., 2011) and can be easily implemented for FSI simulation. To factor in wall motion during numerical simulation, surface mesh numbers and their connectivity must match at each time step. By exploiting consistent topology within the patient-specific model, surface meshes at other time frames are generated based on the corresponding LV geometry, and the correlation between the 3D model and the surface mesh at the first time step. Cubic-spline interpolation is commonly applied, as the frame rate of MRI is inadequate for numerical simulation.
In simulations that use the immersed boundary method (Peskin, 2003), a Cartesian mesh is typically applied to the blood domain, where the LV surface consists of triangular facets. In other words, the discretization step can be skipped by using the patient-specific model directly.
Boundary conditions at inlet and outlet are initialized with appropriate blood properties prior to solving the Navier-Stokes equations. Cardiac flow profiles measured from velocity-encoded imaging modalities such as phase-contrast MRI scans can be used to specify inlet and outlet boundary conditions (Wong et al., 2017). Aortic/pulmonary pressure can only be measured invasively, which might not always be possible.
In addition, the heart is a complex multi-scale system involving the interaction of cardiac electrophysiology with muscle tissue, rapid valve opening and closure, and large wall deformation during the cardiac cycle (Quarteroni et al., 2017). To consider the coupling between electrophysiology and heart mechanics, fiber orientation/ architecture information is necessary for modeling electrical conduction and associated force generation (Crozier et al., 2016). As in-vivo acquisition of the fiber architecture is difficult (Toussaint et al., 2013), “Rule-based” methods was used to assign fiber orientation to ventricular cardiac models, which did not include the fiber structure within endocardial and intramural structures (Bayer et al., 2012). “Atlas-based” approaches used a mesh warping process to map the meshes associated with diffusion tensor MRI fiber data on an idealized template mesh. In this way, the fiber architecture information can be automatically incorporated into new patient-specific model (Vadakkumpadan et al., 2012).
The appropriate choice of constitutive models and material parameters to represent valves in computational studies is another topic related to the heart valve mechanics. In the study of Rausch et al. (2013), mitral valves were modeled as follows: (i) the Neo-Hooken isotropic nonlinear hyperplastic model (neglecting the anisotropic microstructure of mitral leaflet tissue; (ii) the coupled May-Newman model to characterize the heterogeneous response of the entire mitral valve complex; (iii) the decoupled Holzapfel model to represent the anisotropic properties of arterial tissue. The last two models resulted in smaller local displacement errors relative to the first model.
Cardiac tissue properties are the other important parameters for modeling the cardiac multi-scale interaction. However in-vitro measurement of tissue properties using a cardiac tissue mechanics testing system (Golob et al., 2014) is not applicable to the circumference in most cases. Estimating the material stiffness from pressure-volume loop analysis might be a practical way (Pironet et al., 2013). Modeling flow-mediated thrombus generation becomes feasible by coupling the hemodynamic equations to the biochemical convection–diffusion–reaction (CDR) equations (Mittal et al., 2016). Challenges of patient-specific intra-cardiac flow simulation are summarized in Table 3.
Table 3. Current challenges and opportunities in numerical simulation of intra-cardiac flow.
Figure 5 summarizes various post-processing analyses of intra-cardiac flow. The flow mapping facilitates visualization of flow patterns throughout the cardiac cycle. Characteristic of diastolic flow is the clockwise anterior vortex, which is believed to preserve momentum and direct the flow toward the outflow tract (Chnafa et al., 2014). The flow mapping depicts one 2D slice of the complex 3D LV flow at a time. Therefore, 3D vortex structure is superior, and has been widely applied in the literature.
Figure 5. Post-processing of intra-cardiac flow showing (A) flow mapping (Seo et al., 2014); (B) vortex structure (Seo et al., 2014); (C) kinetic energy (Seo and Mittal, 2013); and (D) flow component (Svalbring et al., 2016).
In a normal subject, a vortex ring is observed during the rapid filling phase, which results from flow separation at the mitral valve tip (Seo et al., 2014). The vortex ring interacts with endocardium and dissipates (i.e., breaking into smaller eddies during diastasis). Another vortex ring is again generated during atrial contraction, albeit weaker in terms of penetration distance. Both 2D flow mapping and 3D vortex structure provide quantitative analyses of LV flow.
A number of parameters have been proposed to quantify kinetic energy and energy loss during cardiac motion. Pumping efficiency during systole is defined as the ratio of total flux of mechanical energy through the aorta and power expended on blood flow by heart muscle (Seo et al., 2014). In contrast, myocardial efficiency based on the pressure-volume loop is defined as the ratio of systole work and energy consumption of heart muscle. Kinetic energy dissipation is closely related to flow field, and is the energy dissipated into heat. Relatively higher energy dissipation is expected in regions with complex flow fields due to vortex formation. Since energy dissipation is sensitive to flow field, it has been used to diagnose ventricular dysfunction (Mangual et al., 2013).
To quantify blood transit across the LV, numerous trajectories of massless particles are obtained over a few cardiac cycles (Pedrizzetti and Domenichini, 2014). According to flow trajectories, end-diastolic volume is categorized into volumes of four subgroups: (i) direct flow, blood that passes through ventricle during one heartbeat; (ii) retained flow, blood that enters during diastole but does not exit at end of systole in one heartbeat; (iii) delayed ejection flow, blood that enters ventricle in earlier cardiac cycles but exits in current cycle; and (iv) residual volume, blood that resides in ventricle over a number of cardiac cycles. In-vivo studies demonstrate direct and residual flow (≥2 cycles) constitute about 37 and 30%, respectively, of end-diastolic volume, respectively (Eriksson et al., 2011). Although it is believed that larger residual volume may promote thrombogenesis, clinical significance of the relative distribution of the four volumes of subgroups has not been elucidated. The considerable differences among numerical studies (Mangual et al., 2013; Seo and Mittal, 2013) could be due to uncertainties, differing assumptions and the methodologies adopted.
Challenges and Opportunities in Use of Patient-Specific CFD Technologies for Diagnosing Heart Dysfunction
A promising application of numerical models is found in surgical and interventional planning for predicting procedural outcomes. Corsini and coworkers compared two surgical options for a patient with single ventricle malformation (Corsini et al., 2014). A secondary application is to advance the understanding of the effects of myocardial disease and surgery on ventricular flow. Su and colleagues deduced that HCM retarded the formation of vortex ring during diastole because the narrowed LV chamber delayed flow separation and escalated energy dissipation (Su et al., 2014a). Although various effects of myocardial disease on ventricular flow have been investigated, sample sizes in these numerical studies were too small to provide meaningful insights.
Since the 1990s, numerical studies have been conducted to study heart function (Hunter et al., 2003; Sugiura et al., 2012) including intraventricular flows. The early studies were focused on ideal models due to limitations of imaging techniques, particularly for segmenting images from noisy data. More contemporary published patient-specific numerical studies are summarized in Table 4, in which the majority of studies focus on the LV rather than the RV (where the geometry is more complex than LV). As shown in Table 4, very few studies employed reconstructed 3D models from echocardiographic images (Mangual et al., 2013). Among imaging modalities, MRI is the most widely used owing to superior soft tissue contrast, and absence of ionizing radiation. In our experience, the accuracy of current commercial (e.g., 3D Doctor, Mimics, SliceOmatic, Amira) and open (e.g., ITK, ITK-SNAP, VTK, Analyze) segmentation software packages depends on the experience of the user who processes the MRI images. This is a hurdle to wider adoption of numerical studies. Most studies focus on normal subjects, with relatively few myocardial disease cases, for example, right ventricle (RV) dysfunction (Mangual et al., 2012), dilated cardiomyopathy (Mangual et al., 2013), hypertrophic cardiomyopathy (HCM) (Su et al., 2014a), single ventricle (Corsini et al., 2014), and myocardial infarction (Khalafvand et al., 2012) (Figure 6).
Table 4. Published patient-specific CFD simulations of heart ventricles.
Figure 6. CFD simulation of intra-cardiac flow showing (A) vortex formation in dilated cardiomyopathy (DCM) (Mangual et al., 2013); (B) the distribution of blood velocity before and after surgical ventricular restoration (SVR) (Khalafvand et al., 2014); (C) vortex formation in right ventricle (Mangual et al., 2012); (D) vortex formation in hypertrophic cardiomyopathy (HCM); (E) distribution of blood velocity in a single ventricle (SV) (SVC, superior vena cava; RPA, right pulmonary artery; LPA, left pulmonary artery); (F) the distribution of blood velocity in the LV with myocardial infarction (MI) (Khalafvand et al., 2012).
Although valve leaflet dynamics influence intraventricular flow development relative to vortex formation mechanisms and penetration depth, they are not taken into consideration in most studies as current spatiotemporal resolution of MRI and CT images are inadequate for assessing the rapid opening and closing of thin heart valves. Although the mitral valve was incorporated by Seo et al. (2014), the model was not a full patient-specific model. Only the initial geometry was based on patient-specific CT data, while the remaining geometries were simply dilated according to an ideal model. In addition, the mitral valve motion was pre-defined rather than based on patient-specific data. Mihalef et al. (2011) studied a patient with severe mitral stenosis and regurgitation. The expected strong forward jet toward the apex during diastole and the reversal toward the atrium during systole balanced out, preempting a mismatch between the leaflet dynamics and LV volumetric changes in this selected case.
One feasible solution to model valvular motion is the FSI (Khalafvand et al., 2011; Domenichini and Pedrizzetti, 2014; Doost et al., 2016), which to date has been applied only in 2D studies (Table 4). Basically, there are two methods for modeling myocardium deformation during a cardiac cycle: dynamic mesh and immersed boundary method. Volumetric mesh deforms to cope with the motion of the boundary (e.g., the ventricle) in the dynamic mesh method; while in the immersed boundary method, modeling is accomplished using stationary grids and adding force near the boundary in the Navier-Stokes equations to take ventricular wall effects into account. Although the immersed boundary method avoids the issues of potential meshing failure during mesh generation and deformation, additional functions must be added to obtain the solution, which results in extra computational cost.
Conclusions
Regulatory authorities such as US Food and Drug Administration have recognized the value of computer modeling and simulation in the regulatory approval process (Malinauskas et al., 2017). Authorities encourage use of the simulation to complement bench, animal and human testing. There are an increasing number of FDA applications that include simulations. The Heart Flow's FFRCT software has been approved by the US FDA for measuring coronary blockages non-invasively. Because most healthcare practitioners and organizations are unfamiliar to the technical, computational and simulation methodologies of patient-specific CFD, and the methodology is not yet fully developed, there is an understandable hesitancy to embrace the approach. This is in addition to the presently unrealized ability of researchers and practitioners advocates to effectively communicate the potential benefits of patient-specific CFD. One possible reason for the lack of interest by clinicians could be the lack of a validation protocol in general. A general validation protocol would stipulate procedures and methods for measuring a specified clinical quantity using a standard technique and comparing it to a CFD computation. However, the virtual absence of specific clinical quantities with recognized links to the vast majority of pathologies—with the exception of FFR for PCI—upon which decisions could be predicted, creates a particularly challenging obstacle in the validation of patient-specific models. In addition, development of patient-specific models is a time-consuming task that requires patient-specific geometries, material properties, and realistic boundary conditions. These represent formidable challenges, but at the same time significant opportunities to interject precision medicine into clinical practice for improved clinical outcomes.
Author Contributions
LZ, J-MZ and BS: conception or design of the work; LZ, J-MZ, BS, RT, JA, and GK: draft the work or revise it critically for important intellectual content. All authors have seen and approved the final version of manuscript.
Funding
This work has been partially supported by National Medical Research Council through project (grant number NMRC/BnB/0017/2015) and (grant number NMRC/OFIRG/0018/2016). This research was also supported by Biomedical Research Council through research grant (14/1/32/24/002).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
Alkadhi, H., Stolzmann, P., Scheffel, H., Desbiolles, L., Baumüller, S., Plass, A., et al. (2008). Radiation dose of cardiac dual-source CT: the effect of tailoring the protocol to patient-specific parameters. Eur. J. Radiol. 68, 385–391. doi: 10.1016/j.ejrad.2008.08.015
Antoniadis, A., Mortier, P., Kassab, G., Dubini, G., Foin, N., Murasato, Y., et al. (2015). Biomechanical modeling to improve coronary artery bifurcation stenting: expert review document on techniques and clinical implementation. JACC Cardiovasc. Interv. 8, 1281–1296. doi: 10.1016/j.jcin.2015.06.015
Aueron, F. M., and Gruentzig, A. R. (1984). Percutaneous transluminal coronary angioplasty: indication and current status. Prim Cardiol. 10, 97–107.
Bavo, A. M., Pouch, A. M., Degroote, J., Vierendeels, J., Gorman, J. H., Gorman, R. C., et al. (2016). Patient-specific CFD models for intraventricular flow analysis from 3D ultrasound imaing: comparison of three clinical cases. J. Biomech. 50, 144–150. doi: 10.1016/j.jbiomech.2016.11.039
Bayer, J. D., Blake, R. C., Plank, G., and Trayanova, N. (2012). A novel rule-based algorithm for assigning myocardial fiber orientation to computational heart models. Ann. Biomed. Eng. 40, 2243–2254. doi: 10.1007/s10439-012-0593-5
Bekkers, E. J., and Taylor, C. A. (2008). Multiscale vascular surface model generation from medical imaging data using hierarchical features. IEEE Trans. Med. Imaging 27, 331–341. doi: 10.1109/TMI.2007.905081
Belinha, J. (2016). Meshless methods: the future of computational biomechanical simulation. J. Biom. Biostat. 7:325. doi: 10.4172/2155-6180.1000325
Bourdin, X. X., Torsseille, P., Petit, P., and Beillas, P. (2007). “Comparison of tetrahedral and hexahedral meshes for organ finite element modeling: an application kidney impact,” in 20th International Technical Conference on the Enhanced Safety of Vehicles (Lyon: NHTSA).
Burkhoff, D., Mirsky, I., and Suga, H. (2005). Assessment of systolic and diastolic ventricular properties via pressure-volume analysis: a guide for clinical, translational, and basic researchers. Am. J. Physiol. Heart Circ. Physiol. 289, H501–H512. doi: 10.1152/ajpheart.00138.2005
Caro, C. G., Fitz-Geraldm, J. M., and Schroter, R. C. (1969). Arterial wall shear and distribution of early atheroma in man. Nature 223, 1159–1160. doi: 10.1038/2231159a0
Chnafa, C., Mendez, S., and Nicoud, F. (2014). Image-based large-eddy simulation in a realistic left heart. Comput. Fluids 94, 173–187. doi: 10.1016/j.compfluid.2014.01.030
Chnafa, C., Mendez, S., and Nicoud, F. (2016). Image-based simulations show important flow fluctuations in a normal left ventricle: what could be the implications? Ann. Biomed. Eng. 44, 3346–3358. doi: 10.1007/s10439-016-1614-6
Cho, Y. I., and Kensey, K. R. (1991). Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: steady flows. Biorheology 28, 241–262. doi: 10.3233/BIR-1991-283-415
Coenen, A., Lubbers, M., Kurata, A., Kono, A., Dedic, A., Chelu, R. G., et al. (2015). Fractional flow reserve computed from noninvasive CT angiography data: diagnostic performance of an on-site clinician-operated computational fluid dynamics algorithm. Radiology 274, 674–683. doi: 10.1148/radiol.14140992
Corsini, C., Baker, C., Kung, E., Schievano, S., Arbia, G., Baretta, A., et al. (2014). An integrated approach to patient-specific predictive modeling for single ventricle heart palliation. Comput. Methods Biomech. Biomed. Engin. 17, 1572–1589. doi: 10.1080/10255842.2012.758254
Crozier, A., Augustin, C. M., Neic, A., Prassl, A. J., Holler, M., Fastl, T. E., et al. (2016). Image-based personalization of cardiac anatomy for coupled electromechanical modeling. Ann. Biomed. Eng. 44, 58–70. doi: 10.1007/s10439-015-1474-5
Dahl, S., and Vierendeels, J. (2012). FSI simulation of asymmetric mitral valve dynamics during diastolic filling. Comput. Methods Biomech. Biomed. Engin. 15, 121–130. doi: 10.1080/10255842.2010.517200
De Vecchi, C., Caudron, J., Dubourg, B., Pirot, N., Lefebvre, V., Bauer, F., et al. (2014). Effect of the ellipsoid shape of the left ventricular outflow tract on the echocardiographic assessment of aortic valve area in aortic stenosis. J. Cardiovasc. Comput. Tomogr. 8, 52–57. doi: 10.1016/j.jcct.2013.12.006
Doblare, M., Cueto, E., Calvo, B., Martinez, M. A., Garcia, J. M., and Cegonino, J. (2005). On the employ of meshless methods in biomechanics. Comput. Methods Appl. Mech. 194, 801–821. doi: 10.1016/j.cma.2004.06.031
Doenst, T., Spiegel, K., Reik, M., Markl, M., Hennig, J., Nitzsche, S., et al. (2009). Fluid-dynamic modeling of the human left ventricle: methodology and application to surgical ventricular reconstruction. Ann. Thorac. Surg. 87, 1187–1195. doi: 10.1016/j.athoracsur.2009.01.036
Domenichini, F., and Pedrizzetti, G. (2014). Asymptotic model of fluid–tissue interaction for mitral valve dynamics. Cardiovasc. Engin. Technol. 6, 95–104. doi: 10.1007/s13239-014-0201-y
Doost, S. N., Ghista, D., Su, B., Zhong, L., and Morsi, Y. S. (2016). Heart blood flow simulation: a perspective review. BioMed. Engin. Online 15:101. doi: 10.1186/s12938-016-0224-8
Doost, S. N., Zhong, L., Su, B., and Morsi, Y. S. (2017). Two-dimensional intraventricular flow pattern visualization using the image-based computational fluid dynamics. Comput. Methods Biomechanics Biomed. Engin. 20, 492–507. doi: 10.1080/10255842.2016.1250891
Douglas, P. S., Pontone, G., Hlatky, M. A., Patel, M. R., Norgaard, B. L., Byrne, R. A., et al. (2015). Clinical outcomes of fractional flow reserve by computed tomographic angiography-guided diagnostic strategies vs. usual care in patients with suspected coronary artery disease: the prospective longitudinal trial of FFRct: outcome and resource impacts stud. Eur. Heart J. 36, 3359–3367. doi: 10.1093/eurheartj/ehv444
Eriksson, J., Dyverfeldt, P., Engvall, J., Bolger, A. F., Ebbers, T., and Carlhäll, C. J. (2011). Quantification of presystolic blood flow organization and energetics in the human left ventricle. Am. J. Physiol. Heart Circ. Physiol. 300, H2135–H2141. doi: 10.1152/ajpheart.00993.2010
Eshtehardi, P., McDaniel, M. C., Suo, J., Dhawan, S. S., Timmins, L. H., Binongo, J. N. G., et al. (2012). Association of coronary wall shear stress with atherosclerotic plaque burden, composition, and distribution in patients with coronary artery disease. J. Am. Heart Assoc. 1:e002543. doi: 10.1161/JAHA.112.002543
Figueroa, C. A., Baek, S., Taylor, C. A., and Humphrey, J. D. (2009). A computational framework for fluid-solid-growth modeling in cardiovascular simulations. Comput. Methods Appl. Mech. Eng. 198, 3583–3602. doi: 10.1016/j.cma.2008.09.013
Figueroa, C. A., Vignon-clementel, I. E., Jansen, K. E., Hughes, T. J. R., and Taylor, C. A. (2006). A coupled momentum method for modeling blood flow in three-dimensional deformable arteries. Comput. Methods Appl. Mech. Engin. 195, 5685–5706. doi: 10.1016/j.cma.2005.11.011
Fry, D. L. (1968). Acute vascular endothelial changes associated with increased blood velocity gradients. Circ. Res. 22, 165–197. doi: 10.1161/01.RES.22.2.165
Galassi, F., Alkhalil, M., Lee, R., Martindale, P., Kharbanda, R. K., Channon, K. M., et al. (2018). 3D reconstruction of coronary arteries from 2D angiographic projections using nonuniform rational basis splines (NURBS) for accurate modelling of coronary stenoses. PLoS ONE 13:e0190650. doi: 10.1371/journal.pone.0190650
Gaur, S., Achenbach, S., Leipsic, J., Mauri, L., Bezerra, H. G., Jensen, J. M., et al. (2013). Rationale and design of the HeartFlowNXT (HeartFlow analysis of coronary blood flow using CT angiography: NeXt sTeps) study. J. Cardiovasc. Comput. Tomogr. 7, 279–288. doi: 10.1016/j.jcct.2013.09.003
Gijsen, F. J., Wentzel, J. J., Thury, A., Mastik, F., Schaar, J. A., Schuurbiers, J. C. H., et al. (2008). Strain distribution over plaques in human coronary arteries relates to shear stress. Am. J. Physiol. Heart Circ. Physiol. 295, H1608–H1614. doi: 10.1152/ajpheart.01081.2007
Glagov, S., Zarins, C. K., Giddens, D. P., and Ku, D. N. (1989). “Mechanical factors in the pathogenesis, localization and evolution of atherosclerotic plaques,” in Diseases of the Arterial Wall, eds B. Camilleri and B. Fiessinger (Berlin: Springer-Verlag), 217–239.
Golob, M., Moss, R. L., and Chesler, N. C. (2014). Cardiac tissue structure, properties, and performance: a materials science perspective. Ann. Biomed. Eng. 42, 2003–2013. doi: 10.1007/s10439-014-1071-z
Hecht, F., and Pironneau, O. (2017). An energy stable monolithic Eulerian fluid-structure finite element method. Int. J. Numer. Methods Fluids 85, 430–446. doi: 10.1002/fld.4388
Himburg, H. A., Grzybowski, D. M., Hazel, A. L., LaMack, J. A., and Friedman, M. H. (2004). Spatial comparison between wall shear stress measures and porcine arterial endothelial permeability. Am. J. Physiol. Heart. Circ. Physiol. 286, 1916–1922. doi: 10.1152/ajpheart.00897.2003
Hunter, P. J., Pullan, A. J., and Smaill, B. H. (2003). Modeling total heart function. Annu. Rev. Biomed. Eng. 5, 147–177. doi: 10.1146/annurev.bioeng.5.040202.121537
Imanparast, A., Fatouraee, N., and Sharif, F. (2016). The impact of valve simplifications on left ventricular hemodynamics in a three dimensional simulation based on in vivo MRI data. J. Biomech. 49, 1482–1489. doi: 10.1016/j.jbiomech.2016.03.021
Iskander, S., and Iskandrian, A. E. (1998). Risk assessment using single-photon emission computed tomographic technetium-99m sestamibi imaging. J. Am. Coll. Cardiol. 32, 57–62. doi: 10.1016/S0735-1097(98)00177-6
Itu, L., Rapaka, S., Passerini, T., Georgescu, B., Schwemmer, C., Schoebinger, M., et al. (2016). A machine-learning approach for computation of fractional flow reserve from coronary computed tomography. J. Appl. Physiol. 121, 42–52. doi: 10.1152/japplphysiol.00752.2015
Johnson, N. P., Kirkeeide, R. L. L., and Gould, K. L. L. (2012). Is discordance of coronary flow reserve and fractional flow reserve due to methodology or clinically relevant coronary pathophysiology? JACC Cardiovasc. Imaging 5, 193–202. doi: 10.1016/j.jcmg.2011.09.020
Jung, J., Lyczkowski, R. W., Panchal, C. B., and Hassanein, A. (2006). Multiphase hemodynamic simulation of pulsatile flow in a coronary artery. J. Biomech. 39, 2064–2073. doi: 10.1016/j.jbiomech.2005.06.023
Kantor, B., Kuzo, R. S., and Gerber, T. C. (2007). Coronary computed tomographic angiography: current and future uses. Heart Metab. 34, 5–9.
Kassab, G. S., and Navia, J. A. (2006). Biomechanical considerations in the design of graft: the homeostasis hypothesis. Annu. Rev. Biomed. Eng. 8, 499–535. doi: 10.1146/annurev.bioeng.8.010506.105023
Kawaji, T., Shiomi, H., Morishita, H., Morimoto, T., Taylor, C. A., Kanao, S., et al. (2017). Feasibility and diagnostic performance of fractional flow reserve measurement derived from coronary computed tomography angiography in real clinical practice. Int. J. Cardiovasc. Imaging 33, 271–281. doi: 10.1007/s10554-016-0995-9
Khalafvand, S. S., Ng, E. Y. K., and Zhong, L. (2011). CFD simulation of flow through heart: a perspective review. Comput. Methods Biomech. Biomed. Engin. 14, 113–132. doi: 10.1080/10255842.2010.493515
Khalafvand, S. S., Ng, E. Y. K., Zhong, L., and Hung, T. K. (2012). Fluid-dynamics modelling of the human left ventricle with dynamic mesh for normal and myocardial infarction: preliminary study. Comput. Biol. Med. 42, 863–870. doi: 10.1016/j.compbiomed.2012.06.010
Khalafvand, S. S., Zhong, L., and Ng, E. Y. K. (2014). Three-dimensional CFD/MRI modeling reveals that ventricular surgical restoration improves ventricular function by modifying intraventricular blood flow. Int. J. Numer. Method. Biomed. Eng. 30, 1044–1056. doi: 10.1002/cnm.2643
Kim, H. J., Vignon-Clementel, I. E., Coogan, J. S., Figueroa, C. A., Jansen, K. E., and Taylor, C. A. (2010). Patient-specific modeling of blood flow and pressure in human coronary arteries. Ann. Biomed. Eng. 38, 3195–3209. doi: 10.1007/s10439-010-0083-6
Kimura, T., Shiomi, H., Kuribayashi, S., Isshiki, T., Kanazawa, S., Ito, H., et al. (2015). Cost analysis of non-invasive fractional flow reserve derived from coronary computed tomographic angiography in Japan. Cardiovasc. Interv. Ther. 30, 38–44. doi: 10.1007/s12928-014-0285-1
Kleinstreuer, C., Hyun, S., Buchanan, J. R., Longest, P. W., Archie, J. P., and Truskey, G. A. (2001). Hemodynamic parameters and early intimal thickening in branching blood vessels. Crit. Rev. Biomed. Eng. 29, 1–64. doi: 10.1615/CritRevBiomedEng.v29.i1.10
Knight, J., Olgac, U., Saur, S. C., Poulikakos, D., Marshall, W. Jr., Cattin, P. C., et al. (2010). Choosing the optimal wall shear parameter for the prediction of plaque location-A patient-specific computational study in human right coronary arteries. Atherosclerosis 211, 445–450. doi: 10.1016/j.atherosclerosis.2010.03.001
Ko, B. S., Cameron, J. D., Munnur, R. K., Wong, D. T. L., Fujisawa, Y., Sakaguchi, T., et al. (2017). Noninvasive CT-derived FFR based on structural and fluid analysis: a comparison with invasive FFR for detection of functionally significant stenosis. JACC Cardiovasc. Imaging 10, 663–673. doi: 10.1016/j.jcmg.2016.07.005
Koo, B. K., Erglis, A., Doh, J.-H., Daniels, D. V., Jegere, S., Kim, H. S., et al. (2011). Diagnosis of ischemia-causing coronary stenoses by noninvasive fractional flow reserve computed from coronary computed tomographic angiograms. Results from the prospective multicenter DISCOVER-FLOW (Diagnosis of Ischemia-Causing Stenoses Obtained via Noninvasive Fractional Flow Reserve) study.J. Am. Coll. Cardiol. 58, 1989–1997. doi: 10.1016/j.jacc.2011.06.066
Krittian, S., Janoske, U., Oertel, H., and Böhlke, T. (2010). Partitioned fluid-solid coupling for cardiovascular blood flow: left-ventricular fluid mechanics. Ann. Biomed. Eng. 38, 1426–1441. doi: 10.1007/s10439-009-9895-7
Ku, D. N., Giddens, D. P., Zarins, C. K., and Glagov, S. (1985). Pulsatile flow and atherosclerosis in the human carotid bifurcation. Positive correlation between plaque location and low oscillating shear stress. Arteriosclerosis 5, 293–302. doi: 10.1161/01.ATV.5.3.293
Langille, B. L., and O'Donnell, F. (1986). Reductions in arterial diameter produced by chronic decreases in blood flow are endothelium-dependent. Science 23, 405–407.
Le, T., and Sotiropoulos, F. (2012). On the three-dimensional vortical structure of early distolic flow in a patient-specific left ventricle. Eur. J. Mech. B Fluids 35, 20–24. doi: 10.1016/j.euromechflu.2012.01.013
Leng, S., Jiang, M., Zhao, X. D., Allen, J. C., Kassab, G. S., Ouyang, R. Z., et al. (2016). Three-dimensional tricuspid annular motion analysis from cardiac magnetic resonance feature-tracking. Ann. Biomed. Eng. 44, 3522–3538. doi: 10.1007/s10439-016-1695-2
Lewis, M. A., Pascoal, A., Keevil, S. F., and Lewis, C. A. (2016). Selecting a CT scanner for cardiac imaging: the heart of the matter. Br. J. Radiol. 89, 569–590. doi: 10.1259/bjr.20160376
Liu, B., Zheng, J., Bach, R., and Tang, D. (2012). Correlations of coronary plaque wall thickness with wall pressure and wall pressure gradient: a representative case study. Biomed. Eng. Online 11:43. doi: 10.1186/1475-925X-11-43
Long, Q., Merrifield, R., Xu, X. Y., Kilner, P., Firmin, D. N., and G-Z, Y. (2008). Subject-specific computational simulation of left ventricular flow based on magnetic resonance imaging. Proc. Inst. Mech. Eng. H. 222, 475–485. doi: 10.1243/09544119JEIM310
Lu, M. T., Ferencik, M., Roberts, R. S., Lee, K. L., Ivanov, A., Adami, E., et al. (2017). Noninvasive FFR derived from coronary CT angiography: management and outcomes in the PROMISE trial. JACC Cardiovasc. Imaging 10, 1350–1358. doi: 10.1016/j.jcmg.2016.11.024
Malinauskas, R. A., Hariharan, P., Day, S. W., Herbertson, L. H., Buesen, M., Steinseifer, U., et al. (2017). FDA benchmark medical device flow models for CFD validation. ASAIO J. 63, 150–160. doi: 10.1097/MAT.0000000000000499
Malvè, M., García, A., Ohayon, J., and Martínez, M. A. (2012). Unsteady blood flow and mass transfer of a human left coronary artery bifurcation: FSI vs. CFD. Int. Commun. Heat Mass Transf. 39, 745–751. doi: 10.1016/j.icheatmasstransfer.2012.04.009
Mangual, J. O., Domenichini, F., and Pedrizzetti, G. (2012). Three dimensional numerical assessment of the right ventricular flow using 4D echocardiography boundary data. Eur. J. Mech. B Fluids 35, 25–30. doi: 10.1016/j.euromechflu.2012.01.022
Mangual, J. O., Kraigher-Krainer, E., De Luca, A., Toncelli, L., Shah, A., Solomon, S., et al. (2013). Comparative numerical study on left ventricular fluid dynamics after dilated cardiomyopathy. J. Biomech. 46, 1611–1617. doi: 10.1016/j.jbiomech.2013.04.012
Mihalef, V., Ionasec, R. I., Sharma, P., Georgescu, B., Voigt, I., Suehling, M., et al. (2011). Patient-specific modelling of whole heart anatomy, dynamics and haemodynamics from four-dimensional cardiac CT images. Interface Focus. 1, 286–296. doi: 10.1098/rsfs.2010.0036
Min, J. K., Leipsic, J., Pencina, M. J., Berman, D. S., Koo, B. K., van Mieghem, C., et al. (2012). Diagnostic accuracy of fractional flow reserve from anatomic CT angiography. JAMA. 308, 1237–1245. doi: 10.1001/2012.jama.11274
Mittal, R., Seo, J. H., Vedula, V., Choi, Y. J., Liu, H., Huang, H. H., et al. (2016). Computational modeling of cardiac hemodynamics: current status and future outlook. J. Comput. Phys. 305, 1065–1082. doi: 10.1016/j.jcp.2015.11.022
Moosavi, M. H., Fatouraee, N., Katoozian, H., Pashaei, A., Camara, O., and Frangi, A. E. (2014). Numerical simulation of blood flow in the left ventricle and aortic sinus using magnetic resonance imaging and computational fluid dyanmics. Comput. Methods Biomech. Biomed. Engin. 17, 740–749. doi: 10.1080/10255842.2012.715638
Morris, P. D., Narracott, A., von Tengg-Kobligk, H., Soto, D. A. S., Hsiao, S., Lungu, A., et al. (2015). Computational fluid dynamics modelling in cardiovascular medicine. Heart 102, 18–28. doi: 10.1136/heartjnl-2015-308044
Morris, P. D., Ryan, D., Morton, A. C., Lycett, R., Lawford, P. V., Hose, D. R., et al. (2013). Virtual fractional flow reserve from coronary angiography: modeling the significance of coronary lesions. Results from the VIRTU-1 (VIRTUal fractional flow reserve from coronary angiography) study. JACC Cardiovasc. Interv. 6, 149–157. doi: 10.1016/j.jcin.2012.08.024
Nguyen, V. T., Loon, C. J., Nguyen, H. H., Liang, Z., and Leo, H. L. (2013). A semi-automated method for patient-specific computational flow modelling of left ventricles. Comput. Methods Biomech. Biomed. Engin. 18, 401–413. doi: 10.1080/10255842.2013.803534
Nishimura, R. A., and Tajik, A. J. (1997). Evaluation of diastolic filling of left ventricle in health and disease: doppler echocardiography is the clinician's Rosetta Stone. J. Am. Coll. Cardiol. 30, 8–18.
Nørgaard, B., Leipsic, J., Gaur, S., Seneviratne, S., Ko, B. S., Ito, H., et al. (2014). Diagnostic performance of non-invasive fractional flow reserve derived from coronary CT angiography in suspected coronary artery disease: the NXT trial. J. Am. Coll. Cardiol. 63, 1145–1155. doi: 10.1016/j.jacc.2013.11.043
Papafaklis, M. I., Bourantas, C. V., Theodorakis, P. E., Katsouras, C. S., Fotiadis, D. I., and Michalis, L. K. (2007). Association of endothelial shear stress with plaque thickness in a real three-dimensional left main coronary artery bifurcation model. Int. J. Cardiol. 115, 276–278. doi: 10.1016/j.ijcard.2006.04.030
Pedrizzetti, G., and Domenichini, F. (2014). Left ventricular fluid mechanics: the long way from theoretical models to clinical applications. Ann. Biomed. Eng. 43, 26–40. doi: 10.1007/s10439-014-1101-x
Peskin, C. S. (2003). The immersed boundary method. Acta. Numer. 11, 479–517. doi: 10.1017/S0962492902000077
Pironet, A., Desaive, T., Kosta, S., Lucas, A., Paeme, S., Collet, A., et al. (2013). A multi-scale cardiovascular system model can account for the load-dependence of the end-systolic pressure-volume relationship. Biomed. Eng. Online 12:8. doi: 10.1186/1475-925X-12-8
Poon, E. K. W., Hayat, U., Thondapu, V., Ooi, A. S. H., Haq, M. A. U., Moore, S., et al. (2015). Advances in three-dimensional coronary imaging and computational fluid dynamics: Is virtual fractional flow reserve more than just a pretty picture? Coron. Artery Dis. 26, e43–e54. doi: 10.1097/MCA.0000000000000219
Poon, E. K. W., Thondapu, B., Revalor, E., Ooi, A., and Barlis, P. (2017). Coronary optical coherence tomography-derived virtual fractional flow reserve (FFR): anatomy and physiology all-in-one. Eur. Heart J. 38, 3604–3605. doi: 10.1093/eurheartj/ehx594
Prosi, M., Perktold, K., Ding, Z., and Friedman, M. H. (2004). Influence of curvature dynamics on pulsatile coronary artery flow in a realistic bifurcation model. J. Biomech. 37, 1767–1775. doi: 10.1016/j.jbiomech.2004.01.021
Qiu, Y., and Tarbell, J. M. (1996). Computational simulation of flow in the end-to-end anastomosis of a rigid graft and a compliant artery. ASAIO J. 42, M702–M709. doi: 10.1097/00002480-199609000-00078
Quarteroni, A., Lassila, T., Rossi, S., and Ruiz-Baier, R. (2017). Integrated heart—coupling multiscale and multiphysics models for the simulation of the cardiac function. Comput. Methods Appl. Mech. Eng. 314, 345–407. doi: 10.1016/j.cma.2016.05.031
Rausch, M. K., Famaey, N., Shultz, T. O. B., Bothe, W., Miller, D. C., and Kuhl, E. (2013). Mechanics of the mitral valve: a critical review, an in vivo parameter identification, and the effect of prestrain. Biomech. Model. Mechanobiol. 12, 1053–1071. doi: 10.1007/s10237-012-0462-z
Rutsch, F., Nitschke, Y., and Terkeltaub, R. (2011). Genetics in arterial calcification: pieces of a puzzle and cogs in a wheel. Circ. Res. 109, 578–592. doi: 10.1161/CIRCRESAHA.111.247965
Saber, N. R., Gosman, A. D., Wood, N. B., Kilner, P. J., Charrier, C. L., and Firmin, D. N. (2001). Computational flow modeling of the left ventricle based on in vivo MRI data: initial experience. Ann. Biomed. Eng. 29, 275–283. doi: 10.1114/1.1359452
Saber, N. R., Wood, N. B., Gosman, A. D., Merrifield, R. D., Yang, G. Z., Charrier, C. L., et al. (2003). Progress towards patient-specific computational flow modeling of the left heart via combination of magnetic resonance imaging with computational fluid dynamics. Ann. Biomed. Eng. 31, 42–52. doi: 10.1114/1.1533073
Saeed, M., Van, T. A., Krug, R., Hetts, S. W., and Wilson, M. W. (2015). Cardiac MR imaging: current status and future direction. Cardiovasc. Diagn. Ther. 5, 290–310. doi: 10.3978/j.issn.2223-3652.2015.06.07
Sankaran, S., and Marsden, A. L. (2011). A stochastic collocation method for uncertainty quantification and propagation in cardiovascular simulations. J. Biomech. Eng. 133:031001. doi: 10.1115/1.4003259
Sankaran, S., Grady, L., and Taylor, C. A. (2015). Impact of geometric uncertainty on hemodynamic simulations using machine learning. Comput. Methods Appl. Mech. Eng. 297, 167–190. doi: 10.1016/j.cma.2015.08.014
Sankaran, S., Kim, H. J., Choi, G., and Taylor, C. A. (2016). Uncertainty quantification in coronary blood flow simulations: impact of geometry, boundary conditions and blood viscosity. J. Biomech. 49, 2540–2547. doi: 10.1016/j.jbiomech.2016.01.002
Sankaran, S., Moghadam, M. E., Kahn, A. M., Tseng, E. E., Guccione, J. M., and Marsden, A. L. (2012). Patient-specific multiscale modeling of blood flow for coronary artery bypass graft surgery. Ann. Biomed. Eng. 40, 2228–2242. doi: 10.1007/s10439-012-0579-3
Schenkel, T., Malve, M., Reik, M., Markl, M., Jung, B., and Oertel, H. (2009). MRI-based CFD analysis of flow in a human left ventricle: methodology and application to a healthy heart. Ann. Biomed. Eng. 37, 503–515. doi: 10.1007/s10439-008-9627-4
Schmidt, J. P., Delp, S. L., Sherman, M. A., Taylor, C. A., Pande, V. S., and Altman, R. B. (2008). The simbios national center: systems biology in motion. Proc. IEEE. 96, 1266–1280. doi: 10.1109/JPROC.2008.925454
Seo, J. H., and Mittal, R. (2013). Effect of diastolic flow patterns on the function of the left ventricle. Phys. Fluids 25:110801. doi: 10.1063/1.4819067
Seo, J. H., Vedula, V., Abraham, T., Lardo, A. C., Dawoud, F., Luo, H., et al. (2014). Effect of the mitral valve on diastolic flow patterns. Phys. Fluids 26:121901. doi: 10.1063/1.4904094
Spiegel, M., Redel, T., Zhang, Y. J., Struffert, T., Hornegger, J., Grossman, R. G., et al. (2011). Tetrahedral vs. polyhedral mesh size evaluation on flow velocity and wall shear stress for cerebral hemodynamic simulation. Comput. Methods Biomech. Biomed. Engin. 14, 9–22. doi: 10.1080/10255842.2010.518565
Spilker, R. L., Feinstein, J. A., Parker, D. W., Reddy, V. M., and Taylor, C. A. (2007). Morphometry-based impedance boundary conditions for patient-specific modeling of blood flow in pulmonary arteries. Ann. Biomed. Eng. 35, 546–559. doi: 10.1007/s10439-006-9240-3
Stone, P. H., Coskun, A. U., Kinlay, S., Clark, M. E., Sonka, M., Wahle, A., et al. (2003). Effect of endothelial shear stress on the progression of coronary artery disease, vascular remodeling, and in-stent restenosis in humans: in vivo 6-month follow-up study. Circulation 108, 438–444. doi: 10.1161/01.CIR.0000080882.35274.AD
Stone, P. H., Saito, S., Takahashi, S., Makita, Y., Nakamura, S., Kawasaki, T., et al. (2012). Prediction of progression of coronary artery disease and clinical outcomes using vascular profiling of endothelial shear stress and arterial plaque characteristics: the PREDICTION study. Circulation 126, 172–181. doi: 10.1161/CIRCULATIONAHA.112.096438
Su, B., Kabinejadian, F., Phang, H. Q., Kumar, G. P., Cui, F., Kim, S., et al. (2015). Numerical modeling of intraventricular flow during diastole after implantation of BMHV. PLoS ONE 10:e0126315. doi: 10.1371/journal.pone.0126315
Su, B., Tan, R. S., Tan, J. L., Guo, K. W. Q., Zhang, J. M., Leng, S., et al. (2016). Cardiac MRI based numerical modeling of left ventricular fluid dynamics with mitral valve incorporated. J. Biomech. 49, 1199–1205. doi: 10.1016/j.jbiomech.2016.03.008
Su, B., Zhang, J. M., Tang, H. C., Wan, M. Lim, C. C. W., Su, Y., et al. (2014a). “Patient-specific blood flows and vortex formations in patients with hypertrophic cardiomyopathy using computational fluid dynamics,” in Biomedical Engineering and Sciences (IECBES), Conference on IEEE, (Sarawak) 276–280.
Su, B., Zhong, L., Wang, X. K., Zhang, J. M., Tan, R. S., Allen, J. C., et al. (2014b). Numerical simulation of patient-specific left ventricular model with both mitral and aortic valves by FSI approach. Comput. Methods Programs Biomed. 113, 474–482. doi: 10.1016/j.cmpb.2013.11.009
Sugiura, S., Washio, T., Hatano, A., Okada, J., Watanabe, H., and Hisada, T. (2012). Multi-scale simulations of cardiac electrophysiology and mechanics using the University of Tokyo heart simulator. Prog. Biophys. Mol. Biol. 110, 380–389. doi: 10.1016/j.pbiomolbio.2012.07.001
Svalbring, E., Fredriksson, A., Eriksson, J., Dyverfeldt, P., Ebbers, T., Bolger, A. F., et al. (2016). Altered diastolic flow patterns and kinetic energy in subtle left ventricular remodeling and dysfunction detected by 4D flow MRI. PLoS ONE 11:e0161391. doi: 10.1371/journal.pone.0161391.
Tay, W. B., Tseng, Y. H., Lin, L. Y., and Tseng, W. Y. (2011). Towards patient-specific cardiovascular modeling system using the immersed boundary technique. Biomed. Eng. Online 10:52. doi: 10.1186/1475-925X-10-52
Torii, R., Keegan, J., Wood, N. B., Dowsey, A. W., Hughes, A. D., Yang, G. Z., et al. (2010). MR image-based geometric and hemodynamic investigation of the right coronary artery with dynamic vessel motion. Ann. Biomed. Eng. 38, 2606–2620. doi: 10.1007/s10439-010-0008-4
Torii, R., Wood, N. B., Hadjiloizou, N., Dowsey, A. W., Wright, A. R., Hughes, A. D., et al. (2009). Stress phase angle depicts differences in coronary artery hemodynamics due to changes in flow and geometry after percutaneous coronary intervention. Am. J. Physiol. Heart Circ. Physiol. 296, H765–H776. doi: 10.1152/ajpheart.01166.2007
Toussaint, N., Stoeck, C. T., Schaeffter, T., Kozerke, S., Sermesant, M., and Batchelor, P. G. (2013). In vivo human cardiac fibre architecture estimation using shape-based diffusion tensor processing. Med. Image Anal. 17, 1243–1255. doi: 10.1016/j.media.2013.02.008
Tu, S., Barbato, E., Köszegi, Z., Yang, J., Sun, Z., Holm, N. R., et al. (2014). Fractional flow reserve calculation from 3-dimensional quantitative coronary angiography and TIMI frame count: a fast computer model to quantify the functional significance of moderately obstructed coronary arteries. JACC Cardiovasc. Interv. 7, 768–777. doi: 10.1016/j.jcin.2014.03.004
Tu, S., Westra, J., Yang, J., von Birgelen, C., Ferrara, A., Pellicano, M., et al. (2016). Diagnostic accuracy of fast computational approaches to derive fractional flow reserve from diagnostic coronary angiography: the international multicenter FAVOR pilot study. JACC Cardiovasc. Interv. 9, 2024–2035. doi: 10.1016/j.jcin.2016.07.013
Vadakkumpadan, F., Arevalo, H., Ceritoglu, C., Miller, M., and Trayanova, N. (2012).Image-based estimation of ventricular fiber orientations for personalized modeling of cardiac electrophysiology. IEEE Trans. Med. Imaging 31, 1051–1060. doi: 10.1109/TMI.2012.2184799
Valentín, A., Humphrey, J. D., and Holzapfel, G. A. (2013). A finite element-based constrained mixture implementation for arterial growth, remodeling, and adaptation: theory and numerical verification. Int. J. Numer. Method. Biomed. Eng. 29, 822–849. doi: 10.1002/cnm.2555
Vedula, V., George, R., Younes, L., and Mittal, R. (2015). Hemodynamics in the left atrium and its effect on ventricular flow patterns. J. Biomech. Eng. 137:111003. doi: 10.1115/1.4031487
Westra, J., Tu, S., Winther, S., Nissen, L., Vestergaard, M.-B., Andersen, B. K., et al. (2018). Evaluation of coronary artery stenosis by quantitative flow ratio during invasive coronary angiography. Circ. Cardiovasc. Imaging 11:e007107. doi: 10.1161/CIRCIMAGING.117.007107
Wittek, A., Grosland, N. M., Joldes, G. R., Magnotta, V., and Miller, K. (2016). From finite element meshes to clouds of points: a review of methods for generation of computational biomechanics models for patient-specific applications. Ann. Biomed. Eng. 44, 3–15. doi: 10.1007/s10439-015-1469-2
Wittek, A., Joldes, G., Couton, M., Warfield, S. K., and Miller, K. (2010). Patient-specific non-linear finite element modelling for predicting soft organ deformation in real-time; application to non-rigid neuroimage registration. Prog. Biophys. Mol. Biol. 103, 292–303. doi: 10.1016/j.pbiomolbio.2010.09.001
Wong, K. K., Wang, D., Ko, J. K. L., Mazumdar, J., Le, T. T., and Ghista, D. (2017). Computational medical imaging and hemodynamics framework for functional analysis and assessment of cardiovascular structures. Biomed. Eng. Online 16:35. doi: 10.1186/s12938-017-0326-y
World Health Organization (2011). Burden: Mortality, Morbidity and Risk Factors. Global Status Report on Noncommunicable Diseases 2010, 9–31.
Xu, B., Tu, S., Qiao, S., Qu, X., Chen, Y., Yang, J., et al. (2017). Diagnostic accuracy of angiography-Based quantitative flow ratio measurements for online assessment of coronary stenosis. J. Am. Coll. Cardiol. 70, 3077–3087. doi: 10.1016/j.jacc.2017.10.035
Yamada, H., Oki, T., Tabata, T., Iuchi, A., and Ito, S. (1998). Assessment of left ventricular systolic wall motion velocity with pulsed tissue Doppler imaging: comparison with peak dP/dt of the left ventricular pressure curve. J. Am. Soc. Echocardiogr. 11, 442–449. doi: 10.1016/S0894-7317(98)70024-0
Zhang, J. M., Luo, T., Tan, S. Y., Lomarda, A. M., Wong, A. S., Keng, F. Y., et al. (2015). Hemodynamic analysis of patient-specific coronary artery tree. Int. J. Numer. Methods Biomed. Eng. 31:e02708. doi: 10.1002/cnm.2708
Zhang, J. M., Shuang, D., Baskaran, L., Wu, W., Teo, S. K., Huang, W., et al. (2018). Advanced analyses of computed tomography coronary angiography can help discriminate ischemic lesions. Int. J. Cardiol. doi: 10.1016/j.ijcard.2018.04.020. [Epub ahead print].
Zhang, J. M., Zhong, L., Luo, T., Lomarda, A. M., Huo, Y., Yap, J., et al. (2016). Simplified models of non-invasive fractional flow reserve based on CT images. PLoS ONE 11:e0153070. doi: 10.1371/journal.pone.0153070
Zhang, J. M., Zhong, L., Su, B., Wan, M., Yap, J. S., Tham, J. P., et al. (2014). Perspective on CFD studies of coronary artery disease lesions and hemodynamics: a review. Int. J. Numer. Methods Biomed. Eng. 30, 659–680. doi: 10.1002/cnm.2625
Zhang, X., Zhang, P., and Zhang, L. (2012). A simple technique to improve computational efficiency of meshless methods. Proc. Engin. 31, 1102–1107. doi: 10.1016/j.proeng.2012.01.1149
Zhong, L., Gobeawan, L., Su, Y., Tan, J. L., Ghista, D., Chua, T., et al. (2012). Right ventricular regional wall curvedness and area strain in patients with repaired tetralogy of fallot. Am. J. Physiol. Heart Circ. Physiol. 302, H1306–H1316. doi: 10.1152/ajpheart.00679.2011
Zhong, L., Huang, F. Q., Tan, L. K., Allen, J. C., Ding, Z. P., Kassab, G., et al. (2014). Age and gender-specific changes in left ventricular systolic function in human volunteers. Int. J. Cardiol. 172, e102–e105. doi: 10.1016/j.ijcard.2013.12.128
Keywords: blood flow, computational fluid dynamics (CFD), patient-specific, cardiovascular, coronary, intra-cardiac flow simulation
Citation: Zhong L, Zhang J-M, Su B, Tan RS, Allen JC and Kassab GS (2018) Application of Patient-Specific Computational Fluid Dynamics in Coronary and Intra-Cardiac Flow Simulations: Challenges and Opportunities. Front. Physiol. 9:742. doi: 10.3389/fphys.2018.00742
Received: 04 December 2017; Accepted: 28 May 2018;
Published: 26 June 2018.
Edited by:
Viatcheslav Gurev, IBM, United StatesReviewed by:
Stephen Michael Moore, IBM Research, AustraliaSethuraman Sankaran, HeartFlow, Inc., United States
Copyright © 2018 Zhong, Zhang, Su, Tan, Allen and Kassab. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Liang Zhong, zhong.liang@nhcs.com.sg
Ru San Tan, tan.ru.san@singhealth.com.sg
†These authors have contributed equally to this work and are joint first author.