The numerical analysis of non-Newtonian blood flow in human patient-specific left ventricle

Effect of RVAD Cannulation Length on Right Ventricular Thrombosis Risk: An In Silico Investigation

Left ventricular assist devices (LVADs) have been used off-label as long-term support of the right heart due to the lack of a clinically approved durable right VAD (RVAD). Whilst various techniques to reduce RVAD inflow cannula protrusion have been described, the implication of the protrusion length on right heart blood flow and subsequent risk of thrombosis remains poorly understood. This study investigates the influence of RVAD diaphragmatic cannulation length on right ventricular thrombosis risk using a patient-specific right ventricle in silico model validated with particle image velocimetry. Four cannulation lengths (5, 10, 15 and 25 mm) were evaluated in a one-way fluid-structure interaction simulation with boundary conditions generated from a lumped parameter model, simulating a biventricular supported condition. Simulation results demonstrated that the 25-mm cannulation length exhibited a lower thrombosis risk compared to 5-, 10- and 15-mm cannulation lengths due to improved flow energy distribution (25.2%, 24.4% and 17.8% increased), reduced stagnation volume (72%, 68% and 49% reduction), better washout rate (13.0%, 11.6% and 9.1% faster) and lower blood residence time (6% reduction). In the simulated scenario, our findings suggest that a longer RVAD diaphragmatic cannulation length may be beneficial in lowering thrombosis risk; however, further clinical studies are warranted.

Impact of hematocrit on pulsatile blood flow in stenosed arteries: a computational study in healthy, diabetic, and anemic models

In this study, researchers aim to enhance the realism of circulatory system simulations, focusing on factors affecting flow variations, particularly in stenotic arteries of individuals with altered hematocrit levels. Through extensive data collection and varied conditions, the goal is to attain more precise and valid results. The study conducts approximate simulations to comprehensively describe the dynamic motion of pulsatile flow. Different values of inlet velocity (UDF) are introduced, considering potential arterial distortion or occlusion due to plaque deposition, along with variations in hematocrit (Hct) levels commonly observed in patients. Three distinct types of pulsatile blood flow, corresponding to diabetes (Hct 65%), healthy (Hct 45%), and anemia (Hct 25%), are studied and compared. The research illuminates that stenosis in arteries with varying hematocrit levels significantly impacts hydrodynamic features, potentially predisposing individuals to cardiovascular diseases. Through meticulous analysis, several conclusions about hemodynamic characteristics are drawn. It is observed that both velocity and wall shear stress exhibit variation along the affected artery, influenced by stenosis and changes in hematocrit levels. Notably, the highest influence on velocity and wall shear stress is observed with Hct 65%, compared to Hct 45% and Hct 25% at the moment of stenosis. These findings hold substantial practical implications for the field of cardiovascular health, providing valuable insights into blood flow behavior in stenotic arteries with diverse hematocrit levels. Ultimately, this research contributes to more effective clinical interventions.

Prediction of Hemodynamic-Related Hemolysis in Carotid Stenosis and Aiding in Treatment Planning and Risk Stratification Using Computational Fluid Dynamics

Atherosclerosis affects human health in many ways, leading to disability or premature death due to ischemic heart disease, stroke, or limb ischemia. Poststenotic blood flow disruption may also play an essential role in artery wall impairment linked with hemolysis related to shear stress. The maximum shear stress in the atherosclerotic plaque area is the main parameter determining hemolysis risk. In our work, a 3D internal carotid artery model was built from CT scans performed on patients qualified for percutaneous angioplasty due to its symptomatic stenosis. The obtained stenosis geometries were used to conduct a series of computer simulations to identify critical parameters corresponding to the increase in shear stress in the arteries. Stenosis shape parameters responsible for the increase in shear stress were determined. The effect of changes in the carotid artery size, length, and degree of narrowing on the change in maximum shear stress was demonstrated. Then, a correlation for the quick initial diagnosis of atherosclerotic stenoses regarding the risk of hemolysis was developed. The developed relationship for rapid hemolysis risk assessment uses information from typical non-invasive tests for treated patients. Practical guidelines have been developed regarding which stenosis shape parameters pose a risk of hemolysis, which may be adapted in medical practice.

Hemolysis of red blood cells in blood vessels modeled via computational fluid dynamics

The research aims to verify the universal relationship between vessel shape and the risk of hemolysis using a rheological model of blood reflecting the physiological processes related to blood for any blood vessel. Blood is a multi-component fluid, the rheology of which depends on many factors, such as the concentration of red blood cells and local shear stress, which significantly affect the process of hemolysis. Blood rheology models used so far cannot be used for all flows and geometries. Therefore, a new rheology model has been introduced suitable for modeling hemolytic flows observed in arteries with atherosclerotic lesions in the in vivo environment. The previously presented model also has advantages in modeling local viscosity in stenosis. Geometries of the blood vessels from computed tomography scans and simplified models of the actual arteries observed during medical procedures were used in the calculations. Population Balance Based Rheology model predicts the concentration of single, deagglomerated red blood cells and the concentration and size of red blood cell agglomerates, which affect blood rheology and hemolysis. Based on the simulations carried out, a correlation was found between the shape of the vessel cavity and the risk of hemolysis. Presented results can be used in the future to create a correlation between the shape of the atherosclerotic lesions and the risk of hemolysis in the blood to make an initial risk assessment for a given patient.

Determining flow stasis zones in the intracranial aneurysms and the relation between these zones and aneurysms' aspect ratios after flow diversions

BACKGROUND

Flow diverter stents (FDSs) are widely used to treat aneurysms in the clinic. However, even the same flow diverter (FD) use on different patients' aneurysm sites can cause unexpected hemodynamics at the aneurysm region yielding low success rates for the overall treatment. Therefore, the present study aims to unfold why FDs do not work as they are supposed to for some patients and propose empirical correlation along with a contingency table analysis to estimate the flow stasis zones in the aneurysm sacs.

METHODS

The present work numerically evaluated the use of FRED4518 FDS on six patients' intracranial aneurysms based on patient-specific aneurysm geometries. Computational fluid dynamics (CFD) simulation results were further processed to identify the time evolution of weightless blood particles for six patients' aneurysms.

RESULTS

Stagnation zone formation, incoming and outgoing blood flow at the aneurysm neck, and statistical analysis of six patients indicated that FRED4518 showed a large flow stasis zone for an aspect ratio larger than 0.75. However, FRED4518, used for aneurysms with an aspect ratio of less than 0.65, caused small stagnant flow zones based on the number of blood particles that stayed in the aneurysm sac.

CONCLUSION

A patient-specific empirical equation is derived considering aneurysms' morphological characteristics to determine the amount of stagnated fluid flow zones and magnitude of the mean aneurysm velocity in the aneurysm sac for FRED4518 based on weightless fluid particle results for the first time in the literature. As a result, numerical simulation results and patient data-driven equation can help perceive stagnated fluid zone amount before FRED4518 placement by shedding light on neuro-interventional surgeons and radiologists.

Significance of aortoseptal angle anomalies to left ventricular hemodynamics and subaortic stenosis: A numerical study

PURPOSE

Discrete subaortic stenosis (DSS) is an obstructive cardiac disease caused by a membranous lesion in the left ventricular (LV) outflow tract (LVOT). Although its etiology is unknown, the higher prevalence of DSS in LVOT anatomies featuring a steep aortoseptal angle (AoSA) suggests a potential role for hemodynamics. Therefore, the objective of this study was to quantify the impact of AoSA steepening on the LV three-dimensional (3D) hemodynamic stress environment.

METHODS

A 3D LV model reconstructed from cardiac cine-magnetic resonance imaging was connected to four LVOT geometrical variations spanning the clinical AoSA range (115°-160°). LV hemodynamic stresses were characterized in terms of cycle-averaged pressure, temporal shear magnitude (TSM), and oscillatory shear index. The wall shear stress (WSS) topological skeleton was further analyzed by computing the scaled divergence of the WSS vector field.

RESULTS

AoSA steepening caused an increasingly perturbed subaortic flow marked by LVOT flow skewness and complex 3D secondary flow patterns. These disturbances generated WSS overloads (>45% increase in TSM vs. 160° model) on the inferior LVOT wall, and increased WSS contraction (>66% decrease in WSS divergence vs. 160° model) in regions prone to DSS membrane formation.

CONCLUSIONS

AoSA steepening generated substantial hemodynamic stress abnormalities in LVOT regions prone to DSS formation. Further studies are needed to assess the possible impact of such mechanical abnormalities on the tissue and cellular responses.

The pulsatile 3D-Hemodynamics in a doubly afflicted human descending abdominal artery with iliac branching

The study of patient-specific human arterial flow dynamics is well known to face challenges like a) apt geometric modelling, b) bifurcation zone meshing, and c) capturing the hemodynamic prone to variations with multiple disease complications. Due to aneurysms and stenosis in the same arterial network, the blood flow dynamics get affected, which needs to be explored. This study develops a new protocol for accurate geometric modelling, bifurcation zone meshing and numerically investigates the arterial network with abdominal aortic aneurysms (AAA) and right internal iliac stenosis (RIIAS). A realistic arterial model is reconstructed from the computed tomography (CT) data of a human subject. To understand the combined effect of the aneurysm and aortoiliac occlusive diseases in a patient, an arterial network with AAA, RIIAS, multiple branches tapering, and curvature has been considered. Clinically significant pulsatile blood flow simulations have been carried out to trace the alteration in the flow dynamics with multiple pathological complications under consideration. The transient blood flow dynamics are investigated via wall shear stress, wall pressure, velocity contour, streamlines, vorticity, and swirling strength. During the systolic deceleration phase, the rhythmic nested rapid secondary oscillatory WSS, adverse pressure gradients, high WSS, and high WP bands are noticed. Also, the above studies will help researchers, clinicians, and doctors understand the influence of morphological changes on hemodynamics in cardiovascular studies.

Non-Newtonian blood rheology impacts left atrial stasis in patient-specific simulations

The lack of mechanically effective contraction of the left atrium (LA) during atrial fibrillation (AF) disturbs blood flow, increasing the risk of thrombosis and ischemic stroke. Thrombosis is most likely in the left atrial appendage (LAA), a small narrow sac where blood is prone to stagnate. Slow flow promotes the formation of erythrocyte aggregates in the LAA, also known as rouleaux, causing viscosity gradients that are usually disregarded in patient-specific simulations. To evaluate these non-Newtonian effects, we built atrial models derived from 4D computed tomography scans of patients and carried out computational fluid dynamics simulations using the Carreau-Yasuda constitutive relation. We examined six patients, three of whom had AF and LAA thrombosis or a history of transient ischemic attacks (TIAs). We modeled the effects of hematocrit and rouleaux formation kinetics by varying the parameterization of the Carreau-Yasuda relation and modulating non-Newtonian viscosity changes based on residence time. Comparing non-Newtonian and Newtonian simulations indicates that slow flow in the LAA increases blood viscosity, altering secondary swirling flows and intensifying blood stasis. While some of these effects are subtle when examined using instantaneous metrics like shear rate or kinetic energy, they are manifested in the blood residence time, which accumulates over multiple heartbeats. Our data also reveal that LAA blood stasis worsens when hematocrit increases, offering a potential new mechanism for the clinically reported correlation between hematocrit and stroke incidence. In summary, we submit that hematocrit-dependent non-Newtonian blood rheology should be considered when calculating patient-specific blood stasis indices by computational fluid dynamics.

Simulation of Mechanical Heart Valve Dysfunction and the Non-Newtonian Blood Model Approach

The mechanical heart valve (MHV) is commonly used for the treatment of cardiovascular diseases. Nonphysiological hemodynamic in the MHV may cause hemolysis, platelet activation, and an increased risk of thromboembolism. Thromboembolism may cause severe complications and valve dysfunction. This paper thoroughly reviewed the simulation of physical quantities (velocity distribution, vortex formation, and shear stress) in healthy and dysfunctional MHV and reviewed the non-Newtonian blood flow characteristics in MHV. In the MHV numerical study, the dysfunction will affect the simulation results, increase the pressure gradient and shear stress, and change the blood flow patterns, increasing the risks of hemolysis and platelet activation. The blood flow passes downstream and has obvious recirculation and stagnation region with the increased dysfunction severity. Due to the complex structure of the MHV, the non-Newtonian shear-thinning viscosity blood characteristics become apparent in MHV simulations. The comparative study between Newtonian and non-Newtonian always shows the difference. The shear-thinning blood viscosity model is the basics to build the blood, also the blood exhibiting viscoelastic properties. More details are needed to establish a complete and more realistic simulation.

Steady-state solutions of one-dimensional equations of non-Newtonian hemodynamics

This paper is devoted to obtaining and analysis of steady-state solutions of one-dimensional equations for the simulation of blood flow when the non-Newtonian nature of blood is taken into account. The models, based on the rheological relations, widely used for the blood, are considered. The expressions for the nonlinear frictional term are presented. For the Power Law, Simplified Cross, and Quemada models, the exact integrals of the nonlinear ordinary differential equation, obtained from the averaged momentum equation, are obtained. It is demonstrated that several solutions exist for every rheological model, but the physically relevant solutions can be selected by the appropriate value of Mach number. The effects of the velocity profile and the value of hematocrit on the steady-state solutions are analyzed. It is demonstrated that the flattening of the velocity profile, which is typical for the blood, leads to the diminishing of the length of the interval, where the solution exists. The same effect is observed when the hematocrit value is increased.

Sequential Coupling Shows Minor Effects of Fluid Dynamics on Myocardial Deformation in a Realistic Whole-Heart Model

Background: The human heart is a masterpiece of the highest complexity coordinating multi-physics aspects on a multi-scale range. Thus, modeling the cardiac function in silico to reproduce physiological characteristics and diseases remains challenging. Especially the complex simulation of the blood's hemodynamics and its interaction with the myocardial tissue requires a high accuracy of the underlying computational models and solvers. These demanding aspects make whole-heart fully-coupled simulations computationally highly expensive and call for simpler but still accurate models. While the mechanical deformation during the heart cycle drives the blood flow, less is known about the feedback of the blood flow onto the myocardial tissue. Methods and Results: To solve the fluid-structure interaction problem, we suggest a cycle-to-cycle coupling of the structural deformation and the fluid dynamics. In a first step, the displacement of the endocardial wall in the mechanical simulation serves as a unidirectional boundary condition for the fluid simulation. After a complete heart cycle of fluid simulation, a spatially resolved pressure factor (PF) is extracted and returned to the next iteration of the solid mechanical simulation, closing the loop of the iterative coupling procedure. All simulations were performed on an individualized whole heart geometry. The effect of the sequential coupling was assessed by global measures such as the change in deformation and-as an example of diagnostically relevant information-the particle residence time. The mechanical displacement was up to 2 mm after the first iteration. In the second iteration, the deviation was in the sub-millimeter range, implying that already one iteration of the proposed cycle-to-cycle coupling is sufficient to converge to a coupled limit cycle. Conclusion: Cycle-to-cycle coupling between cardiac mechanics and fluid dynamics can be a promising approach to account for fluid-structure interaction with low computational effort. In an individualized healthy whole-heart model, one iteration sufficed to obtain converged and physiologically plausible results.

Computational fluid dynamics of blood flow in an idealized left human heart

We construct an idealized computational model of the left human heart for the study of the blood flow dynamics in the left atrium and ventricle. We solve the Navier-Stokes equations in the ALE formulation and we prescribe the left heart wall displacement based on physiological data; moreover, we consider the presence of both the mitral and aortic valves through the resistive method. We simulate the left heart hemodynamics by means of the finite element method and we consider the variational multiscale large eddy simulation (LES) formulation to account for the transitional and nearly turbulent regimes of the blood flow in physiological conditions. The main contribution of this paper is the characterization of the blood flow in an idealized configuration of the left heart aiming at reproducing function in normal conditions. Our assessment is based on the analysis of instantaneous and phase averaged velocity fields, blood pressure, and other clinically meaningful fluid dynamics indicators. Finally, we show that our idealized computational model can be suitably used to study and critically discuss pathological scenarios like that of a regurgitant mitral valve.

Comparison of Non-Newtonian Models of One-Dimensional Hemodynamics

The paper is devoted to the comparison of different one-dimensional models of blood flow. In such models, the non-Newtonian property of blood is considered. It is demonstrated that for the large arteries, the small parameter is observed in the models, and the perturbation method can be used for the analytical solution. In the paper, the simplified nonlinear problem for the semi-infinite vessel with constant properties is solved analytically, and the solutions for different models are compared. The effects of the flattening of the velocity profile and hematocrit value on the deviation from the Newtonian model are investigated.

Multiphysics computational modelling of the cardiac ventricles

Development of cardiac multiphysics models has progressed significantly over the decades and simulations combining multiple physics interactions have become increasingly common. In this review, we summarise the progress in this field focusing on various approaches of integrating ventricular structures. electrophysiological properties, myocardial mechanics, as well as incorporating blood hemodynamics and the circulatory system. Common coupling approaches are discussed and compared, including the advantages and shortcomings of each. Currently used strategies for patient-specific implementations are highlighted and potential future improvements considered.

Effect of non-Newtonian fluid rheology on an arterial bypass graft: A numerical investigation guided by constructal design

In post-operative scenarios of arterial graft surgeries to bypass coronary artery stenosis, fluid dynamics plays a crucial role. Problems such as intimal hyperplasia have been related to fluid dynamics and wall shear stresses near the graft junction. This study focused on the question of the use of Newtonian and non-Newtonian models to represent blood in this type of problem in order to capture important flow features, as well as an analysis of the performance of geometry from the view of Constructive Theory. The objective of this study was to investigate the effects rheology on the steady-state flow and on the performance of a system consisting of an idealized version of a partially obstructed coronary artery and bypass graft. The Constructal Design Method was employed with two degrees of freedom: the ratio between bypass and artery diameters and the junction angle at the bypass inlet. The flow problem was solved numerically using the Finite Volume Method with blood modeled employing the Carreau equation for viscosity. The Computational Fluid Dynamics model associated with the Sparse Grid method generated eighteen response surfaces, each representing a severe stenosis degree of 75% for specific combinations of rheological parameters, dimensionless viscosity ratio, Carreau number and flow index at two distinct Reynolds numbers of 150 and 250. There was a considerable dependence of the pressure drop on rheological parameters. For the two Reynolds numbers studied, the Newtonian case presented the lowest value of the dimensionless pressure drop, suggesting that the choice of applying Newtonian blood may underestimate the value of pressure drop in the system by about 12.4% (Re =150) and 7.8% (Re = 250). Even so, results demonstrated that non-Newtonian rheological parameters did not influence either the shape of the response surfaces or the optimum bypass geometry, which consisted of a diameter ratio of 1 and junction angle of 30°. However, the viscosity ratio and the flow index had the greatest impact on pressure drop, recirculation zones and wall shear stress. Rheological parameters also affected the recirculation zones downstream of stenosis, where intimal hyperplasia is more prevalent. Newtonian and most non-Newtonian results had similar wall shear stresses, except for the non-Newtonian case with high viscosity ratio. In the view of Constructal Design, the geometry of best performance was independent of the rheological model. However, rheology played an important role on pressure drop and flow dynamics, allowing the prediction of recirculation zones that were not captured by a Newtonian model.

On the relevance of boundary conditions and viscosity models in blood flow simulations in patient-specific aorto-coronary bypass models

Physiologically realistic results are the aim of every blood flow simulation. This is not different in aorto-coronary bypasses where the properties of the coronary circulation may significantly affect the relevance of the performed simulations. By considering three patient-specific bypass geometries, the present article focuses on two aspects of the coronary blood flow - its phasic flow pattern and its behaviour affected by blood rheology. For the phasic flow property, a multiscale modelling approach is chosen as a means to assess the ability of five different types of coronary boundary conditions (mean arterial pressure, Windkessel model and three lumped parameter models) to attain realistic coronary haemodynamics. From the analysed variants of boundary conditions, the best option in terms of physiological characteristics and its potential for use in patient-based simulations, is utilised to account for the effect of shear-dependent viscosity on the resulting haemodynamics and wall shear stress stimulation. Aside from the Newtonian model, the blood rheology is approximated by two non-Newtonian models in order to determine whether the choice of a viscosity model is important in simulations involving coronary circulation. A comprehensive analysis of obtained results demonstrated notable superiority of all lumped parameter models, especially in comparison to the constant outlet pressure, which regardless of bypass type gave overestimated and physiologically misleading results. In terms of rheology, it was noted that blood in undamaged coronary arteries behaves as a Newtonian fluid, whereas in vessels with atypical lumen geometry, such as that of anastomosis or stenosis, its shear-thinning behaviour should not be ignored.

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

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

4D flow evaluation of blood non-Newtonian behavior in left ventricle flow analysis

Blood is generally modeled as a Newtonian fluid, assuming a standard and constant viscosity; however, this assumption may not hold for the highly pulsatile and recirculating intracavitary flow in the left ventricle (LV), hampering the quantification of fluid dynamic indices of potential clinical relevance. Herein, we investigated the effect of three viscosity models on the patient-specific quantification of LV blood energetics, namely on viscous energy loss (EL), from 4D Flow magnetic resonance imaging: I) Newtonian with standard viscosity (3.7 cP), II) Newtonian with subject-specific hematocrit-dependent viscosity, III) non-Newtonian accounting for the effect of hematocrit and shear rate. Analyses were performed on 5 controls and 5 patients with cardiac light-chain amyloidosis. In Model II, viscosity ranged between 3.0 (-19%) and 4.3 cP (+16%), mildly deviating from the standard value. In the non-Newtonian model, this effect was emphasized: viscosity ranged from 3.2 to 6.0 cP, deviating maximally from the standard value in low shear rate (i.e., <100 s^-1) regions. This effect reflected on EL quantifications: in particular, as compared to Model I, Model III yielded markedly higher EL values (up to +40%) or markedly lower (down to -21%) for subjects with hematocrit higher than 39.5% and lower than 30%, respectively. Accounting for non-Newtonian blood behavior on a patient-specific basis may enhance the accuracy of intracardiac energetics assessment by 4D Flow, which may be explored as non-invasive index to discriminate between healthy and pathologic LV.

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

BACKGROUND AND OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Hemodynamic Analysis of VenaTech Convertible Vena Cava Filter Using Computational Fluid Dynamics

The VenaTech convertible filter (VTCF) has been widely used as an inferior vena cava (IVC) filter to prevent fatal pulmonary embolism in patients. However, its hemodynamics that greatly affect the filter efficacy and IVC patency are still unclear. This paper uses computational fluid dynamics with the Carreau model to simulate the non-Newtonian blood flows around the VTCF respectively deployed in the normal, reverse and three converted states in an IVC model. The results show that the prothrombotic stagnation zones are observed downstream from the normal, reverse and small open VTCFs, with the streamwise length is nearly eight times the IVC diameter. The no-slip boundary conditions of the thin-wire VTCF arms lead to the “viscous block” effect. The viscous block accelerates the blood flow by 5–15% inside the IVC and enhances the filter wall shear stress up to nearly 20 times that of the IVC only, which contributes to clot capture and thrombus lysis. The relative flow resistance is defined to evaluate the filter-induced resistance on the IVC blood flow that can be regarded as an index of IVC patency with the filter deployment. The flow resistance of the normal VTCF deployment increases dramatically by more than 60% compared with that of the IVC only and is a little higher (6%) than that of the reverse case. As the VTCF converts to a fully open configuration, the flow resistance gradually decreases to that of no filter. This work shows that even very thin VTCF arms can result in the viscous block effect and may cause significant hemodynamic impacts on clot capture, potential thrombosis and flow impedance inside the IVC. The present study also shows that CFD is a valuable and feasible in silico tool for analyzing the IVC filter hemodynamics to complement in vivo clinical and in vitro experimental studies.

An Algorithm for the Noninvasive and Personalized Measurement of Microvascular Blood Viscosity Using Physiological Parameters

Blood viscosity is one of the important parameters to characterize hemorheological properties of the human body. Its real-time and dynamic measurement has important physiological significance for studying the development and prevention of chronic diseases. This study researched noninvasive and personalized measurement of microvascular blood viscosity. In the microcirculation capillary network blood flow model, combined with pulse wave parameters, multiple regression analysis was used to fit the simulated radius of personalized physiological blood vessels to calculate the microvascular blood viscosity. The fitted value related to the simulated radius of the physiological blood vessel had a high correlation with the corresponding theoretically derived value (correlation coefficient: 0.904, P ≤ 0.001). The calculated value of the microvascular blood viscosity had a certain correlation with the clinical whole blood viscosity at a low shear rate (correlation coefficient: 0.443, P < 0.05). This algorithm could provide effective means for noninvasive and long-term individual monitoring and family health care.

Understanding the effect of effective metal surface area of flow diverter stent's on the patient-specific intracranial aneurysm numerical model using Lagrangian coherent structures

The effective metal surface area (EMSA) of flow diversions plays an essential role in the occlusion mechanism inside the aneurysm since the value of EMSA determines the amount of blood flow into the aneurysm sac. In the present study, three different models of a flow diverter stent, namely FRED 4017, FRED 4038, and FRED 4539, were virtually placed at the aneurysm neck of a 52-years-old female patient to identify the effect of EMSA on stagnation region formation inside the aneurysm sac. Lagrangian coherent structures (LCSs), hyperbolic time, and particle tracking analysis were employed to the velocity vectors obtained from computational fluid dynamics (CFD). It is noticed that use of FRED 4017 stent with 0.42 EMSA value caused nearly 40% of the weightless blood flow particles (more than FRED 4038 and FRED 4539) to stay inside the aneurysm while only 0.35% of the blood flow was remaining inside the aneurysm sac when no stent was placed into the aneurysm site. Furthermore, hyperbolic time computations illustrated the formation of stagnation fluid flow zones that can be associated with the residence time of the blood flow particles. Lastly, the results of hyperbolic time analysis are in good agreement with digital subtraction angiography (DSA) images taken in the clinic a few minutes after a FRED 4017 implantation.

Non-Newtonian Effects on Patient-Specific Modeling of Fontan Hemodynamics

The Fontan procedure is a common palliative surgery for congenital single ventricle patients. In silico and in vitro patient-specific modeling approaches are widely utilized to investigate potential improvements of Fontan hemodynamics that are related to long-term complications. However, there is a lack of consensus regarding the use of non-Newtonian rheology, warranting a systematic investigation. This study conducted in silico patient-specific modeling for twelve Fontan patients, using a Newtonian and a non-Newtonian model for each patient. Differences were quantified by examining clinically relevant metrics: indexed power loss (iPL), indexed viscous dissipation rate (iVDR), hepatic flow distribution (HFD), and regions of low wall shear stress (A_WSS). Four sets of "non-Newtonian importance factors" were calculated to explore their effectiveness in identifying the non-Newtonian effect. No statistical differences were observed in iPL, iVDR, and HFD between the two models at the population-level, but large inter-patient variations exist. Significant differences were detected regarding A_WSS, and its correlations with non-Newtonian importance factors were discussed. Additionally, simulations using the non-Newtonian model were computationally faster than those using the Newtonian model. These findings distinguish good importance factors for identifying non-Newtonian rheology and encourage the use of a non-Newtonian model to assess Fontan hemodynamics.

Effects of different non-Newtonian models on unsteady blood flow hemodynamics in patient-specific arterial models with in-vivo validation

The aim of this study is to demonstrate the implications of using different blood rheological models in the simulation of blood flow dynamics in atherosclerotic coronary arteries. Computational fluid dynamics simulation was performed using three-dimensional (3D) patient-specific models of diseased left anterior descending (LAD) coronary arteries with varying degrees of stenosis severity. The three-dimensional arterial models were reconstructed from 3D quantitative coronary angiography, and input flow conditions were prescribed with blood flow conditions measured in-vivo. Different blood viscosity models were used for the simulations, and they include Newtonian and also non-Newtonian models such as Bingham, Carreau, Carreau-Yasuda, Casson, modified Casson, Cross, modified Cross, simplified Cross, Herschel Bulkley, Kuang-Luo (K-L), PowellErying, modified PowellErying, Power-law, Quemada and Walburn-Schneck models. Results from this study show that the time-averaged velocity at the centre of the arteries produced in the CFD simulations that uses the Carreau, modified Casson or Quemada blood viscosity models corresponded exceptionally well with the clinical measurements regardless of stenosis severities and hence, highlights the usefulness of these models to determine the potential determinants of blood vessel wall integrity such as dynamic blood viscosity, blood velocity and wall shear stress.

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

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

The effect of chordae tendineae on systolic flow

When using Computational Fluid Dynamics to simulate ventricular blood flow in the heart, it has been common practice to neglect the effect of the sub-valvular apparatus and the trabeculae on the flow conditions. In this study, we analyze the effect of neglecting the chordae tendineae on the fluid flow and pressure drop. To test the assumption we use a previously developed dynamic 3D model of the left ventricle, aorta and valves that is based on 3D echocardiographic recordings. To this model we add the chordae tendineae as a sub-grid model. The previously developed 3D model for the left ventricle during systole is based on real-time three-dimensional echocardiography (RT3DE) recordings of a 30 years old female volunteer. The segmented ventricular wall does not include details of the aorta and the mitral valve, so these were reconstructed. The subgrid model for the flow across the chordae tendineae is based on the Actuator Line Method, which means that they are represented by drag coefficients. The analysis shows that the effect of the chordae tendineae on the pressure drop and work efficiency of the normal heart during systole is minor, and it seems that for simulating ventricular fluid flow and pressure drop during systole, one can follow the current practice and ignore the chordae. However, there can be local effects such as small vortices behind the chordae. Whether such effects are important for a particular application must be evaluated for the given case.

Left Ventricular Trabeculations Decrease the Wall Shear Stress and Increase the Intra-Ventricular Pressure Drop in CFD Simulations

The aim of the present study is to characterize the hemodynamics of left ventricular (LV) geometries to examine the impact of trabeculae and papillary muscles (PMs) on blood flow using high performance computing (HPC). Five pairs of detailed and smoothed LV endocardium models were reconstructed from high-resolution magnetic resonance images (MRI) of ex-vivo human hearts. The detailed model of one LV pair is characterized only by the PMs and few big trabeculae, to represent state of art level of endocardial detail. The other four detailed models obtained include instead endocardial structures measuring ≥1 mm² in cross-sectional area. The geometrical characterizations were done using computational fluid dynamics (CFD) simulations with rigid walls and both constant and transient flow inputs on the detailed and smoothed models for comparison. These simulations do not represent a clinical or physiological scenario, but a characterization of the interaction of endocardial structures with blood flow. Steady flow simulations were employed to quantify the pressure drop between the inlet and the outlet of the LVs and the wall shear stress (WSS). Coherent structures were analyzed using the Q-criterion for both constant and transient flow inputs. Our results show that trabeculae and PMs increase the intra-ventricular pressure drop, reduce the WSS and disrupt the dominant single vortex, usually present in the smoothed-endocardium models, generating secondary small vortices. Given that obtaining high resolution anatomical detail is challenging in-vivo, we propose that the effect of trabeculations can be incorporated into smoothed ventricular geometries by adding a porous layer along the LV endocardial wall. Results show that a porous layer of a thickness of 1.2·10⁻² m with a porosity of 20 kg/m² on the smoothed-endocardium ventricle models approximates the pressure drops, vorticities and WSS observed in the detailed models.

Hemodynamic forces in the left and right ventricles of the human heart using 4D flow magnetic resonance imaging: Phantom validation, reproducibility, sensitivity to respiratory gating and free analysis software

Purpose

To investigate the accuracy, reproducibility and sensitivity to respiratory gating, field strength and ventricle segmentation of hemodynamic force quantification in the left and right ventricles of the heart (LV and RV) using 4D-flow magnetic resonance imaging (MRI), and to provide free hemodynamic force analysis software.

Materials and methods

A pulsatile flow phantom was imaged using 4D flow MRI and laser-based particle image velocimetry (PIV). Cardiac 4D flow MRI was performed in healthy volunteers at 1.5T (n = 23). Reproducibility was investigated using MR scanners from two different vendors on the same day (n = 8). Subsets of volunteers were also imaged without respiratory gating (n = 17), at 3T on the same day (n = 6), and 1–12 days later on the same scanner (n = 9, median 6 days). Agreement was measured using the intraclass correlation coefficient (ICC).

Results

Phantom validation showed good accuracy for both scanners (Scanner 1: bias -14±9%, y = 0.82x+0.08, R² = 0.96, Scanner 2: bias -12±8%, y = 0.99x-0.08, R² = 1.00). Force reproducibility was strong in the LV (0.09±0.07 vs 0.09±0.07 N, bias 0.00±0.04 N, ICC = 0.87) and RV (0.09±0.06 vs 0.09±0.05 N, bias 0.00±0.03, ICC = 0.83). Strong to very strong agreement was found for scans with and without respiratory gating (LV/RV: ICC = 0.94/0.95), scans on different days (ICC = 0.92/0.87), and 1.5T and 3T scans (ICC = 0.93/0.94).

Conclusion

Software for quantification of hemodynamic forces in 4D-flow MRI was developed, and results show high accuracy and strong to very strong reproducibility for both the LV and RV, supporting its use for research and clinical investigations. The software including source code is released freely for research.

Patient-Specific Left Ventricular Flow Simulations From Transthoracic Echocardiography: Robustness Evaluation and Validation Against Ultrasound Doppler and Magnetic Resonance Imaging

The combination of medical imaging with computational fluid dynamics (CFD) has enabled the study of 3-D blood flow on a patient-specific level. However, with models based on gated high-resolution data, the study of transient flows, and any model implementation into routine cardiac care, is challenging. This paper presents a novel pathway for patient-specific CFD modelling of the left ventricle (LV), using 4-D transthoracic echocardiography (TTE) as input modality. To evaluate the clinical usability, two sub-studies were performed. First, a robustness evaluation was performed, where repeated models with alternating input variables were generated for six subjects and changes in simulated output quantified. Second, a validation study was carried out, where the pathway accuracy was evaluated against pulsed-wave Doppler (100 subjects), and 2-D through-plane phase-contrast magnetic resonance imaging measurements over seven intraventricular planes (6 subjects). The robustness evaluation indicated a model deviation of <12%, with highest regional and temporal deviations at apical segments and at peak systole, respectively. The validation study showed an error of <11% (velocities <10 cm/s) for all subjects, with no significant regional or temporal differences observed. With the patient-specific pathway shown to provide robust output with high accuracy, and with the pathway dependent only on 4-D TTE, the method has a high potential to be used within future clinical studies on 3-D intraventricular flow patterns. To this, future model developments in the form of e.g., anatomically accurate LV valves may further enhance the clinical value of the simulations.

Modeling Left Ventricular Blood Flow Using Smoothed Particle Hydrodynamics

This study aims to investigate the capability of smoothed particle hydrodynamics (SPH), a fully Lagrangian mesh-free method, to simulate the bulk blood flow dynamics in two realistic left ventricular (LV) models. Three dimensional geometries and motion of the LV, proximal left atrium and aortic root are extracted from cardiac magnetic resonance imaging and multi-slice computed tomography imaging data. SPH simulation results are analyzed and compared with those obtained using a traditional finite volume-based numerical method, and to in vivo phase contrast magnetic resonance imaging and echocardiography data, in terms of the large-scale blood flow phenomena usually clinically measured. A quantitative comparison of the velocity fields and global flow parameters between the in silico models and the in vivo data shows a reasonable agreement, given the inherent uncertainties and limitations in the modeling and imaging techniques. The results indicate the capability of SPH as a promising tool for predicting clinically relevant large-scale LV flow information.

Heart function and hemodynamic analysis for zebrafish embryos

The Zebrafish has emerged to become a powerful vertebrate animal model for cardiovascular research in recent years. Its advantages include easy genetic manipulation, transparency, small size, low cost, and the ability to survive without active circulation at early stages of development. Sequencing the whole genome and identifying ortholog genes with human genome made it possible to induce clinically relevant cardiovascular defects via genetic approaches. Heart function and disturbed hemodynamics need to be assessed in a reliable manner for these disease models in order to reveal the mechanobiology of induced defects. This effort requires precise determination of blood flow patterns as well as hemodynamic stress (i.e., wall shear stress and pressure) levels within the developing heart. While traditional approach involves time-lapse brightfield microscopy to track cell and tissue movements, in more recent studies fast light-sheet fluorescent microscopes are utilized for that purpose. Integration of more complicated techniques like particle image velocimetry and computational fluid dynamics modeling for hemodynamic analysis holds a great promise to the advancement of the Zebrafish studies. Here, we discuss the latest developments in heart function and hemodynamic analysis for Zebrafish embryos and conclude with our future perspective on dynamic analysis of the Zebrafish cardiovascular system. Developmental Dynamics 246:868-880, 2017. © 2017 Wiley Periodicals, Inc.

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

The image-based computational fluid dynamics (IB-CFD) technique, as the combination of medical images and the CFD method, is utilized in this research to analyze the left ventricle (LV) hemodynamics. The research primarily aims to propose a semi-automated technique utilizing some freely available and commercial software packages in order to simulate the LV hemodynamics using the IB-CFD technique. In this research, moreover, two different physiological time-resolved 2D models of a patient-specific LV with two different types of aortic and mitral valves, including the orifice-type valves and integrated with rigid leaflets, are adopted to visualize the process of developing intraventricular vortex formation and propagation. The blood flow pattern over the whole cardiac cycle of two models is also compared to investigate the effect of utilizing different valve types in the process of the intraventricular vortex formation. Numerical findings indicate that the model with integrated valves can predict more complex intraventricular flow that can match better the physiological flow pattern in comparison to the orifice-type model.

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.