Coupled Morphological-Hemodynamic Computational Analysis of Type B Aortic Dissection: A Longitudinal Study

Morphometric and Hemodynamic Analysis of the Compressed Iliac Vein

PURPOSE

To investigate the relationship between the morphological structure and hemodynamic properties of the compressed iliac vein and explore the reason for the formation of thrombosis in the compressed iliac vein.

MATERIALS AND METHODS

A total of 11 patients with iliac vein compression syndrome (IVCS) were included in this study, and their iliac veins were reconstructed in 3 dimensions (3D). The morphological structures of the iliac veins (confluence angle, degree of stenosis) were analyzed based on the 3D model. Variations in the hemodynamic properties of the iliac vein were investigated at 4 typical moments in one cardiac cycle, and the relationship between the different morphological configurations and the pressure difference was investigated.

RESULTS

In the region of the compressed iliac vein, the blood flow velocity is accelerated and the pressure changes abruptly accompanied by the increase in pressure difference. Higher time averaged wall shear stress (TAWSS) and lower relative residence time (RRT) appeared in stenosis regions of compressed iliac vein, while TAWSS was low and RRT was large near the stenosis position. There was a strong positive correlation between the degree of stenosis and the pressure difference (r=0.894), and a positive correlation between the confluence angle of the iliac vein and the pressure difference (r=0.638).

CONCLUSION

The morphological structure of the compressed iliac vein has an obvious influence on the hemodynamic surroundings; the pressure difference becomes larger when the degree of stenosis and the confluence angle increase. The iliac vein luminal areas with low TAWSS and high RRT near the compressed location can impede blood flow and lead to accumulation of blood components, which may increase the risk of thrombosis formation and should be fully considered in the treatment of IVCS.

A Monte Carlo Sensitivity Analysis for a Dimensionally Reduced-Order Model of the Aortic Dissection

PURPOSE

Aortic dissection is associated with a high mortality rate. Although computational approaches have shed light on many aspects of the disease, a sensitivity analysis is required to determine the significance of different factors. Because of its complex geometry and high computational expense, the three-dimensional (3D) fluid-structure interaction (FSI) simulation is not a suitable approach for sensitivity analysis.

METHODS

We performed a Monte Carlo simulation (MCS) to investigate the sensitivity of hemodynamic quantities to the lumped parameters of our zero-dimensional (0D) model with numerically calculated lumped parameters. We performed local and global analyses on the effect of the model parameters on important hemodynamic quantities.

RESULTS

The MCS showed that a larger lumped resistance value for the false lumen and the tears result in a higher retrograde flow rate in the false lumen (the coefficient of variation,c v , i = 0.0183 , the sensitivityS X i σ = 0.54 , Spearman's coefficient,ρ s = 0.464 ). For the intraluminal pressure, our results show a significant role in the resistance and inertance of the true lumen (the coefficient of variation,c v , i = 0.0640 , the sensitivityS X i σ = 0.85 , and Spearman's coefficient,ρ s = 0.855 for the inertance of the true lumen).

CONCLUSION

This study highlights the necessity of comparing the results of the local and global sensitivity analyses to understand the significance of multiple lumped parameters. Because of the efficiency of the method, our approach is potentially useful to investigate and analyze medical planning.

Hemodynamic Implications of STABILISE Technique for Aortic Dissection Repair

BACKGROUND

The stent-assisted balloon-induced intimal disruption and relamination (STABILISE) technique for treatment of type B dissection has shown promising clinical results at mid-term. Computational modeling is a way of noninvasively obtaining hemodynamic effects, such as pressure and wall shear stress, leading to a better understanding of potential benefits. Particular areas of interest are (1) the effect of intimal disruption and re-lamination and (2) the effect of the bare metal stent in the visceral aortic segment.

METHODS

Single-center prospective case series. Data from 5 consecutive locally performed cases of STABILISE technique were analyzed. Included cases were type B aortic dissection with or without prior de-branching. The STABILISE procedure had to be performed without 30-day major complications. Preoperative and postoperative imaging data for each patient were transferred to the biomedical engineering team. Each case was reconstructed, meshed, and simulated with computational fluid dynamics using patient-specific data (heart rate, blood pressure, height, and weight). Hemodynamic parameters were then extracted from the simulations.

RESULTS

In all cases, computational analysis showed for postoperative patients: (1) a drop in pressure difference between lumina and (2) lower wall shear stress effects, compared to their preoperative status. These observations were most pronounced in the visceral aortic segment.

CONCLUSIONS

Computational modeling shows favourable changes in the flow dynamics of type B dissection treated using the STABILISE technique. This may suggest protective effects of this technique for long-term aortic healing and cicatrization.

Long-term survival and risk analysis of thoracic endovascular aortic repair for type B aortic dissection

To explore the safety and efficacy of thoracic endovascular aortic repair (TEVAR) in the treatment of patients with type B aortic dissection, and to evaluate the risk factors for long-term mortality. Our study retrospectively evaluated 729 patients with type B aortic dissection, who were divided into the thoracic endovascular aortic repair group and the optimal medical treatment group according to their treatment. In-hospital mortality, death within 30 days, and aortic-related mortality were lower in the thoracic endovascular aortic repair group than in the optimal medical treatment group (p < 0.05). The cumulative overall survival rates for the thoracic endovascular aortic repair group at 1 year, 5 years, and 10 years were 92.5%, 84.1%, and 73.5%, respectively. The Cox analysis found that TEVAR was beneficial in reducing mortality and that a vertical length of the dissection exceeding 150 mm was a risk factor for mortality.

A hemodynamic analysis of fenestrated physician-modified endograft repair for complicated aortic dissections involving the visceral arteries

OBJECTIVE

The aim of this study is to perform patient-specific hemodynamic simulations of the patients with complicated aortic dissection underwent Physician-modified endograft (PMEG) and evaluate the treatment outcome.

METHOD

12 patient-specific models were reconstructed from computed tomography angiography (CTA) data of 6 patients with complicated aortic dissection before and after the PMEG. Hemodynamic simulations were conducted with the same time-varying volumetric flow rate extracted from the literature and 3-element Windkessel model (3 EWM) boundary conditions were applied at the aortic outlet. Hemodynamic indicators such as time-averaged wall shear stress (TAWSS), relative residence time (RRT) and endothelial cell activation potential (ECAP) were obtained to evaluate the postoperative effect of PMEG.

RESULTS

Comparing with the preoperative models, the flow rates of most visceral arteries were increased in the postoperative models (PSMA = 0.012, PRRA = 0.013, and PLRA = 0.005). Pressure and TAWSS in visceral regions were significantly reduced (PP = 0.003 and PTAWSS = 0.017). With the false lumens (FL) covered by the stent grafts, the average TAWSS level increased in the regions of postoperative abdominal aorta (P = 0.002), and the average RRT and ECAP values decreased significantly (PRRT = 0.02 and PECAP = 0.003).

CONCLUSION

This study shows that PMEG, as a new technique for the treatment of complicated aortic dissection involving the distal tears in the visceral region, can effectively restore the abnormal blood supply of the visceral arteries, reduce the risk of aortic rupture, the formation of aortic dissection aneurysm (ADA), and thrombosis. This corresponds well with clinical retrospective studies and 1-year follow-up outcomes. The findings of this study are of great significance for the development of PMEG.

Effect of turbulence and viscosity models on wall shear stress derived biomarkers for aorta simulations

Ascending aorta simulations provide insight into patient-specific hemodynamic conditions. Numerous studies have assessed fluid biomarkers which show a potential to aid clinicians in the diagnosis process. Unfortunately, there exists a large disparity in the computational methodology used to model turbulence and viscosity. Recognizing this disparity, some authors focused on analysing the influence of either the turbulence or viscosity models on the biomarkers in order to quantify the importance of these model choices. However, no analysis has yet been done on their combined effect. In order to fully understand and quantify the effect of the computational methodology, an assessment of the combined effect of turbulence and viscosity model choice was performed. Our results show that (1) non-Newtonian viscosity has greater impact (2.9-5.0%) on wall shear stress than Large Eddy Simulation turbulence modelling (0.1-1.4%), (2) the contribution of non-Newtonian viscosity is amplified when combined with a subgrid-scale turbulence model, (3) wall shear stress is underestimated when considering Newtonian viscosity by 2.9-5.0% and (4) cycle-to-cycle variability can impact the results as much as the numerical model if insufficient cycles are performed. These results demonstrate that, when assessing the effect of computational methodologies, the resultant combined effect of the different modelling assumptions differs from the aggregated effect of the isolated modifications. Accurate aortic flow modelling requires non-Newtonian viscosity and Large Eddy Simulation turbulence modelling.

Aneurysmal growth in type-B aortic dissection: assessing the impact of patient-specific inlet conditions on key haemodynamic indices

Type-B aortic dissection is a cardiovascular disease in which a tear develops in the intimal layer of the descending aorta, allowing pressurized blood to delaminate the layers of the vessel wall. In medically managed patients, long-term aneurysmal dilatation of the false lumen (FL) is considered virtually inevitable and is associated with poorer disease outcomes. While the pathophysiological mechanisms driving FL dilatation are not yet understood, haemodynamic factors are believed to play a key role. Computational fluid dynamics (CFD) and 4D-flow MRI (4DMR) analyses have revealed correlations between flow helicity, oscillatory wall shear stress and aneurysmal dilatation of the FL. In this study, we compare CFD simulations using a patient-specific, three-dimensional, three-component inlet velocity profile (4D IVP) extracted from 4DMR data against simulations with flow rate-matched uniform and axial velocity profiles that remain widely used in the absence of 4DMR. We also evaluate the influence of measurement errors in 4DMR data by scaling the 4D IVP to the degree of imaging error detected in prior studies. We observe that oscillatory shear and helicity are highly sensitive to inlet velocity distribution and flow volume throughout the FL and conclude that the choice of IVP may greatly affect the future clinical value of simulations.

CFD Study of the Effect of the Angle Pattern on Iliac Vein Compression Syndrome

Iliac vein compression syndrome (IVCS, or May-Thurner syndrome) occurs due to the compression of the left common iliac vein between the lumbar spine and right common iliac artery. Because most patients with compression are asymptomatic, the syndrome is difficult to diagnose based on the degree of anatomical compression. In this study, we investigated how the tilt angle of the left common iliac vein affects the flow patterns in the compressed blood vessel using three-dimensional computational fluid dynamic (CFD) simulations to determine the flow fields generated after compression sites. A patient-specific iliac venous CFD model was created to verify the boundary conditions and hemodynamic parameter set in this study. Thirty-one patient-specific CFD models with various iliac venous angles were developed using computed tomography (CT) angiograms. The angles between the right or left common iliac vein and inferior vena cava at the confluence level of the common iliac vein were defined as α1 and α2. Flow fields and vortex locations after compression were calculated and compared according to the tilt angle of the veins. Our results showed that α2 affected the incidence of flow field disturbance. At α2 angles greater than 60 degrees, the incidence rate of blood flow disturbance was 90%. In addition, when α2 and α1 + α2 angles were used as indicators, significant differences in tilt angle were found between veins with laminar, transitional, and turbulent flow (p < 0.05). Using this mathematical simulation, we concluded that the tilt angle of the left common iliac vein can be used as an auxiliary indicator to determine IVCS and its severity, and as a reference for clinical decision making.

Calibration of patient-specific boundary conditions for coupled CFD models of the aorta derived from 4D Flow-MRI

Introduction: Patient-specific computational fluid dynamics (CFD) models permit analysis of complex intra-aortic hemodynamics in patients with aortic dissection (AD), where vessel morphology and disease severity are highly individualized. The simulated blood flow regime within these models is sensitive to the prescribed boundary conditions (BCs), so accurate BC selection is fundamental to achieve clinically relevant results. Methods: This study presents a novel reduced-order computational framework for the iterative flow-based calibration of 3-Element Windkessel Model (3EWM) parameters to generate patient-specific BCs. These parameters were calibrated using time-resolved flow information derived from retrospective four-dimensional flow magnetic resonance imaging (4D Flow-MRI). For a healthy and dissected case, blood flow was then investigated numerically in a fully coupled zero dimensional-three dimensional (0D-3D) numerical framework, where the vessel geometries were reconstructed from medical images. Calibration of the 3EWM parameters was automated and required ~3.5 min per branch. Results: With prescription of the calibrated BCs, the computed near-wall hemodynamics (time-averaged wall shear stress, oscillatory shear index) and perfusion distribution were consistent with clinical measurements and previous literature, yielding physiologically relevant results. BC calibration was particularly important in the AD case, where the complex flow regime was captured only after BC calibration. Discussion: This calibration methodology can therefore be applied in clinical cases where branch flow rates are known, for example, via 4D Flow-MRI or ultrasound, to generate patient-specific BCs for CFD models. It is then possible to elucidate, on a case-by-case basis, the highly individualized hemodynamics which occur due to geometric variations in aortic pathology high spatiotemporal resolution through CFD.

A numerical study of the hemodynamic behavior and gas transport in cardiovascular systems with severe cardiac or cardiopulmonary failure supported by venoarterial extracorporeal membrane oxygenation

Venoarterial extracorporeal membrane oxygenation (VA-ECMO) has been extensively demonstrated as an effective means of bridge-to-destination in the treatment of patients with severe ventricular failure or cardiopulmonary failure. However, appropriate selection of candidates and management of patients during Extracorporeal membrane oxygenation (ECMO) support remain challenging in clinical practice, due partly to insufficient understanding of the complex influences of extracorporeal membrane oxygenation support on the native cardiovascular system. In addition, questions remain as to how central and peripheral venoarterial extracorporeal membrane oxygenation modalities differ with respect to their hemodynamic impact and effectiveness of compensatory oxygen supply to end-organs. In this work, we developed a computational model to quantitatively address the hemodynamic interaction between the extracorporeal membrane oxygenation and cardiovascular systems and associated gas transport. Model-based numerical simulations were performed for cardiovascular systems with severe cardiac or cardiopulmonary failure and supported by central or peripheral venoarterial extracorporeal membrane oxygenation. Obtained results revealed that: 1) central and peripheral venoarterial extracorporeal membrane oxygenation modalities had a comparable capacity for elevating arterial blood pressure and delivering oxygenated blood to important organs/tissues, but induced differential changes of blood flow waveforms in some arteries; 2) increasing the rotation speed of extracorporeal membrane oxygenation pump (ω) could effectively improve arterial blood oxygenation, with the efficiency being especially high when ω was low and cardiopulmonary failure was severe; 3) blood oxygen indices (i.e., oxygen saturation and partial pressure) monitored at the right radial artery could be taken as surrogates for diagnosing potential hypoxemia in other arteries irrespective of the modality of extracorporeal membrane oxygenation; and 4) Left ventricular (LV) overloading could occur when ω was high, but the threshold of ω for inducing clinically significant left ventricular overloading depended strongly on the residual cardiac function. In summary, the study demonstrated the differential hemodynamic influences while comparable oxygen delivery performance of the central and peripheral venoarterial extracorporeal membrane oxygenation modalities in the management of patients with severe cardiac or cardiopulmonary failure and elucidated how the status of arterial blood oxygenation and severity of left ventricular overloading change in response to variations in ω. These model-based findings may serve as theoretical references for guiding the application of venoarterial extracorporeal membrane oxygenation or interpreting in vivo measurements in clinical practice.

Hotspots and development frontiers of postoperative complications of AD: Bibliometric analysis - a review

The research on the postoperative complications of aortic dissection (AD) has received great attention from scholars all over the world, and the number of research articles in this field has consistently increased year after year. However, no bibliometric reports have been published yet to analyze the scientific output and the current situation in this field. The Bibliometrix R-package, VOSviewer, and CiteSpace software were used to conduct a bibliometric analysis of the hotspots and development frontiers of AD. A total of 1242 articles were retrieved. The USA, China, and Japan had the highest number of publications. The five keywords with the highest frequency were "analysis," "incidence," "acute type," "graft," and "risk factor." The results also indicated that the research in related fields had shifted from surgical treatment and utilizing experience to the evidence-based exploration of risk factors and the construction of prediction models to help better manage postoperative complications of AD. This is the first bibliometric analysis of global publications on the postoperative complications of AD. The current research hotspots focus on three areas: common postoperative complications of AD, exploration of the related risk factors, and management of complications. Future research could focus on identifying risk factors through meta-analysis and using a multicenter database for AD as well as building relevant models to predict the development of complications to better facilitate the clinical management of AD patients.

Advanced risk prediction for aortic dissection patients using imaging-based computational flow analysis

Patients with either a repaired or medically managed aortic dissection have varying degrees of risk of developing late complications. High-risk patients would benefit from earlier intervention to improve their long-term survival. Currently serial imaging is used for risk stratification, which is not always reliable. On the other hand, understanding aortic haemodynamics within a dissection is essential to fully evaluate the disease and predict how it may progress. In recent decades, computational fluid dynamics (CFD) has been extensively applied to simulate complex haemodynamics within aortic diseases, and more recently, four-dimensional (4D)-flow magnetic resonance imaging (MRI) techniques have been developed for in vivo haemodynamic measurement. This paper presents a comprehensive review on the application of image-based CFD simulations and 4D-flow MRI analysis for risk prediction in aortic dissection. The key steps involved in patient-specific CFD analyses are demonstrated. Finally, we propose a workflow incorporating computational modelling for personalised assessment to aid in risk stratification and treatment decision-making.

Clinical implications of hemodynamic analysis for the three-dimension iliac vein model with different stenosis

Background

The aim of this study was to perform hemodynamic simulations of a three-dimension ideal inferior vena cava-iliac vein model with artificial stenosis to determine the degree of stenosis that requires clinical intervention.

Methods

Four three-dimension stenosis models (30%, 50%, 70%, and 90% stenosis) were constructed using commercial software (Solidworks). The inlet flow rates were acquired from previous literatures to perform the hemodynamic simulations. Changes in the old blood volume fraction, as well conventional hemodynamic parameters including pressure, differential pressure, wall shear stress, and flow patterns, over time were recorded. The pressure at the telecentric region of the stenosis increased with increasing degree of stenosis.

Results

For the 70% stenosis model, the pressure at the telecentric region of the stenosis reached 341 Pa, and the differential pressure between the two ends of the stenosis was 363 Pa (approximately 2.7 mmHg). Moreover, in the 70% and 90% stenosis models, there was a marked change in wall shear stress in the stenosis and the proximal end region, and the flow patterns began to show the phenomenon of flow separation. Blood stasis analysis showed that the 70% stenosis model had the slowest decrease in old blood volume fraction, while the proximal end region had the largest blood residue (15%).

Conclusion

Iliac vein stenosis of approximately 70% is associated with clinically relevant hemodynamic changes, and is more closely related to DVT than other degrees of stenosis.

Fluid-structure interaction simulation of visceral perfusion and impact of different cannulation methods on aortic dissection

Hemodynamics in aortic dissection (AD) is closely associated with the risk of aortic aneurysm, rupture, and malperfusion. Altered blood flow in patients with AD can lead to severe complications such as visceral malperfusion. In this study, we aimed to investigate the effect of cannulation flow on hemodynamics in AD using a fluid-structure interaction simulation. We developed a specific-idealized AD model that included an intimal tear in the descending thoracic aorta, a re-entry tear in the left iliac artery, and nine branches. Two different cannulation methods were tested: (1) axillary cannulation (AC) only through the brachiocephalic trunk and (2) combined axillary and femoral cannulation (AFC) through the brachiocephalic trunk and the right common iliac artery. AC was found to result in the development of a pressure difference between the true lumen and false lumen, owing to the difference in the flow rate through each lumen. This pressure difference collapsed the true lumen, disturbing blood flow to the celiac and superior mesenteric arteries. However, in AFC, the pressure levels between the two lumens were similar, and no collapse occurred. Moreover, the visceral flow was higher than that in AC. Lastly, the stiffness of the intimal flap affected the true lumen's collapse.

Non-invasive estimation of the parameters of a three-element windkessel model of aortic arch arteries in patients undergoing thoracic endovascular aortic repair

Background: Image-based computational hemodynamic modeling and simulations are important for personalized diagnosis and treatment of cardiovascular diseases. However, the required patient-specific boundary conditions are often not available and need to be estimated. Methods: We propose a pipeline for estimating the parameters of the popular three-element Windkessel (WK3) models (a proximal resistor in series with a parallel combination of a distal resistor and a capacitor) of the aortic arch arteries in patients receiving thoracic endovascular aortic repair of aneurysms. Pre-operative and post-operative 1-week duplex ultrasound scans were performed to obtain blood flow rates, and intra-operative pressure measurements were also performed invasively using a pressure transducer pre- and post-stent graft deployment in arch arteries. The patient-specific WK3 model parameters were derived from the flow rate and pressure waveforms using an optimization algorithm reducing the error between simulated and measured pressure data. The resistors were normalized by total resistance, and the capacitor was normalized by total resistance and heart rate. The normalized WK3 parameters can be combined with readily available vessel diameter, brachial blood pressure, and heart rate data to estimate WK3 parameters of other patients non-invasively. Results: Ten patients were studied. The medians (interquartile range) of the normalized proximal resistor, distal resistor, and capacitor parameters are 0.10 (0.07-0.15), 0.90 (0.84-0.93), and 0.46 (0.33-0.58), respectively, for common carotid artery; 0.03 (0.02-0.04), 0.97 (0.96-0.98), and 1.91 (1.63-2.26) for subclavian artery; 0.18 (0.08-0.41), 0.82 (0.59-0.92), and 0.47 (0.32-0.85) for vertebral artery. The estimated pressure showed fairly high tolerance to patient-specific inlet flow rate waveforms using the WK3 parameters estimated from the medians of the normalized parameters. Conclusion: When patient-specific outflow boundary conditions are not available, our proposed pipeline can be used to estimate the WK3 parameters of arch arteries.

An artificial intelligence-based platform for automatically estimating time-averaged wall shear stress in the ascending aorta

Aims

Aortopathies are a series of disorders requiring multiple indicators to assess risk. Time-averaged wall shear stress (TAWSS) is currently considered as the primary indicator of aortopathies progression, which can only be calculated by Computational Fluid Dynamics (CFD). However, CFD's complexity and high computational cost, greatly limit its application. The study aimed to construct a deep learning platform which could accurately estimate TAWSS in ascending aorta.

Methods and results

A total of 154 patients who had thoracic computed tomography angiography were included and randomly divided into two parts: training set (90%, n = 139) and testing set (10%, n = 15). TAWSS were calculated via CFD. The artificial intelligence (AI)-based model was trained and assessed using the dice coefficient (DC), normalized mean absolute error (NMAE), and root mean square error (RMSE). Our AI platform brought into correspondence with the manual segmentation (DC = 0.86) and the CFD findings (NMAE, 7.8773% ± 4.7144%; RMSE, 0.0098 ± 0.0097), while saving 12000-fold computational cost.

Conclusion

The high-efficiency and robust AI platform can automatically estimate value and distribution of TAWSS in ascending aorta, which may be suitable for clinical applications and provide potential ideas for CFD-based problem solving.

Influence of MRI-based boundary conditions on type B aortic dissection simulations in false lumen with or without abdominal aorta involvement

Most computational hemodynamic studies of aortic dissections rely on idealized or general boundary conditions. However, numerical simulations that ignore the characteristics of the abdominal branch arteries may not be conducive to accurately observing the hemodynamic changes below the branch arteries. In the present study, two men (M-I and M-II) with type B aortic dissection (TBAD) underwent arterial-phase computed tomography angiography and four-dimensional flow magnetic resonance imaging (MRI) before and after thoracic endovascular aortic repair (TEVAR). The finite element method was used to simulate the computational fluid dynamic parameters of TBAD [false lumen (FL) with or without visceral artery involvement] under MRI-specific and three idealized boundary conditions in one cardiac cycle. Compared to the results of zero pressure and outflow boundary conditions, the simulations with MRI boundary conditions were closer to the initial MRI data. The pressure difference between true lumen and FL after TEVAR under the other three boundary conditions was lower than that of the MRI-specific results. The results of the outflow boundary conditions could not characterize the effect of the increased wall pressure near the left renal artery caused by the impact of Tear-1, which raised concerns about the distal organ and limb perfused by FL. After TEVAR, the flow velocity and wall pressure in the FL and the distribution areas of high time average wall shear stress and oscillating shear index were reduced. The difference between the calculation results for different boundary conditions was lower in M-II, wherein FL did not involve the abdominal aorta branches than in M-I. The boundary conditions of the abdominal branch arteries from MRI data might be valuable in elucidating the hemodynamic changes of the descending aorta in TBAD patients before and after treatment, especially those with FL involving the branch arteries.

Blood flow analysis with computational fluid dynamics and 4D-flow MRI for vascular diseases

Both computational fluid dynamics (CFD) and time-resolved, three-dimensional, phase-contrast, magnetic resonance imaging (4D-flow MRI) enable visualization of time-varying blood flow structures and quantification of blood flow in vascular diseases. However, they are totally different. CFD is a method to calculate blood flow by solving the governing equations of fluid mechanics, so the obtained flow field is somewhat virtual. On the other hand, 4D-flow MRI measures blood flow in vivo, thus the flow is real. Recently, with the development and enhancement of computers, medical imaging techniques, and related software, blood flow analysis has become more accessible to clinicians and its usefulness in vascular diseases has been demonstrated. In this review, we have outlined the methods and characteristics of CFD and 4D-flow MRI, respectively. We have discussed the differences in the characteristics between both methods; reviewed the milestones achieved by blood flow analysis in various vascular diseases; and discussed the usefulness, challenges, and limitations of blood flow analysis. We have discussed the difficulties and limitations of current blood flow analysis. We have also discussed our views on future directions.

The Necessity to Seal the Re-Entry Tears of Aortic Dissection After TEVAR: A Hemodynamic Indicator

Thoracic endovascular aortic repair (TEVAR) is a common treatment for Stanford type B aortic dissection (TBAD). However, re-entry tears might be found distal to the stented region which transports blood between the true and false lumens. Sealing the re-entry tears, especially for the thoracic tears, could further reduce blood perfusion to the false lumen; however, it might also bring risks by re-intervention or surgery. Wise determination of the necessity to seal the re-entry tears is needed. In this study, patient-specific models of TBAD were reconstructed, and the modified models were established by virtually excluding the thoracic re-entries. Computational hemodynamics was investigated, and the variation of the functional index and first balance position (FBP) of the luminal pressure difference, due to the sealing of the re-entries, was reported. The results showed that the direction of the net flow through the unstented thoracic re-entries varied among cases. Excluding the re-entries with the net flow toward the false lumen may induce the FBP moving distally and the relative particle residence time increasing in the false lumen. This study preliminarily demonstrated that the hemodynamic status of the re-entry tears might serve as an indicator to the necessity of sealing. By quantifying the through-tear flow exchange and shift of FBP, one can predict the hemodynamic benefit by sealing the thoracic re-entries and thus wisely determine the necessity of further interventional management.

Numerical modeling of residual type B aortic dissection: longitudinal analysis of favorable and unfavorable evolution

Residual type B aortic dissection was numerically investigated to highlight the contribution of biomechanical parameters to the pathology's evolution. Patient-specific geometries from cases involving both favorable and unfavorable evolution were modeled to assess their hemodynamic features. This original approach was supported by a longitudinal study confirming the association between morphological changes, hemodynamic features, adverse clinical outcomes, and CT-angioscan observations on the same patient. Comparing one patient with unfavorable evolution with one with favorable one, we identify potential biomechanical indicators predictive of unfavorable evolution: (i) a patent false lumen with a flow rate above 50% of inlet flow rate; (ii) high wall shear stress above 18 Pa at entry tears, and above 10 Pa at some regions of the false lumen wall; (iii) low time-averaged wall shear stress in distal false lumen below 0.5 Pa; (iv) vortical structure dynamics. Although these comparisons could only be conducted on 2 patients and need to be confirmed by a larger number of cases, our findings point to these hemodynamic markers as possible candidates for early evaluation of the pathology's evolution towards an unfavorable scenario. Graphical Abstract Correlation between hemodynamics index and thrombus initiation for unfavorable case. ET₂ and ET₃ are entry tear numbers 2 and 3 respectively. WSS is wall shear stress. TAWSS is time average shear stress.

Evaluation of Computational Methodologies for Accurate Prediction of Wall Shear Stress and Turbulence Parameters in a Patient-Specific Aorta

Background: Recent studies suggest that blood flow in main arteries is intrinsically disturbed, even under healthy conditions. Despite this, many computational fluid dynamics (CFD) analyses of aortic haemodynamics make the assumption of laminar flow, and best practices surrounding appropriate modelling choices are lacking. This study aims to address this gap by evaluating different modelling and post-processing approaches in simulations of a patient-specific aorta. Methods: Magnetic resonance imaging (MRI) and 4D flow MRI from a patient with aortic valve stenosis were used to reconstruct the aortic geometry and derive patient-specific inlet and outlet boundary conditions. Three different computational approaches were considered based on assumed laminar or assumed disturbed flow states including low-resolution laminar (LR-Laminar), high-resolution laminar (HR-Laminar) and large-eddy simulation (LES). Each simulation was ran for 30 cardiac cycles and post-processing was conducted on either the final cardiac cycle, or using a phase-averaged approach which utilised all 30 simulated cycles. Model capabilities were evaluated in terms of mean and turbulence-based parameters. Results: All simulation types, regardless of post-processing approach could correctly predict velocity values and flow patterns throughout the aorta. Lower resolution simulations could not accurately predict gradient-derived parameters including wall shear stress and viscous energy loss (largest differences up to 44.6% and 130.3%, respectively), although phase-averaging these parameters improved predictions. The HR-Laminar simulation produced more comparable results to LES with largest differences in wall shear stress and viscous energy loss parameters up to 5.1% and 11.6%, respectively. Laminar-based parameters were better estimated than turbulence-based parameters. Conclusion: Our findings suggest that well-resolved laminar simulations can accurately predict many laminar-based parameters in disturbed flows, but there is no clear benefit to running a HR-Laminar simulation over an LES simulation based on their comparable computational cost. Additionally, post-processing "typical" laminar simulation results with a phase-averaged approach is a simple and cost-effective way to improve accuracy of lower-resolution simulation results.

A two-fluid blood stasis model for false lumen thrombosis after type B dissection repair

The formation of thrombosis is a major concern in the false lumen (FL) for post-TEVAR (thoracic endovascular aortic repair) patients. Blood stasis is one of the key factors which lead to the formation of thrombosis in the arterial systems. This study proposed a computational model for blood stasis, using a two-fluid principle to track the locations of blood residual over time. The current study applied this novel model to evaluate blood stasis and thrombosis potential in four patient-specific post-TEVAR FLs of type B aortic dissection, with their follow-up in-vivo observations two years after TEVAR. The locations and topologies of residual blood in the FL predicted by the model agreed well with the in-vivo observations of thrombus. In addition, the results corresponded better with clinical observations in terms of interpatient comparison of degree of thrombosis, compared with conventional hemodynamic parameters. The blood stasis model serves as a valuable addition to conventional metrics to better predict thrombosis potential. Collectively, these metrics can provide an efficient non-invasive method for evaluating blood stasis and thrombosis potential in arterial system, and useful guidance for clinicians' operative planning and postoperative evaluation.

Computational Fluid Dynamic Technique for Assessment of How Changing Character of Blood Flow and Different Value of Hct Influence Blood Hemodynamic in Dissected Aorta

Using computer tomography angiography (CTA) and computational structural analysis, we present a non-invasive method of mass flow rate/velocity and wall stress analysis in type B aortic dissection. Three-dimensional (3D) computer models of the aorta were calculated using pre-operative (baseline) and post-operative CT data from 12 male patients (aged from 51 to 64 years) who were treated for acute type B dissection. A computational fluid dynamics (CFD) technique was used to quantify the displacement forces acting on the aortic wall in the areas of endografts placement. The mass flow rate and wall stress were measured and quantified using the CFD technique. The CFD model indicated the places with a lower value of blood velocity and shear rate, which corelated with higher blood viscosity and a probability of thrombus appearance. Moreover, with the increase in Hct, blood viscosity also increased, while the intensity of blood flow provoked changing viscosity values in these areas. Furthermore, the velocity gradient near the tear surface caused high wall WSS; this could lead to a decreased resistance in the aorta’s wall with further implications to a patient.

Analysis of Turbulence Effects in a Patient-Specific Aorta with Aortic Valve Stenosis

Blood flow in the aorta is often assumed laminar, however aortic valve pathologies may induce transition to turbulence and our understanding of turbulence effects is incomplete. The aim of the study was to provide a detailed analysis of turbulence effects in aortic valve stenosis (AVS).

METHODS

Large-eddy simulation (LES) of flow through a patient-specific aorta with AVS was conducted. Magnetic resonance imaging (MRI) was performed and used for geometric reconstruction and patient-specific boundary conditions. Computed velocity field was compared with 4D flow MRI to check qualitative and quantitative consistency. The effect of turbulence was evaluated in terms of fluctuating kinetic energy, turbulence-related wall shear stress (WSS) and energy loss.

RESULTS

Our analysis suggested that turbulence was induced by a combination of a high velocity jet impinging on the arterial wall and a dilated ascending aorta which provided sufficient space for turbulence to develop. Turbulent WSS contributed to 40% of the total WSS in the ascending aorta and 38% in the entire aorta. Viscous and turbulent irreversible energy losses accounted for 3.9 and 2.7% of the total stroke work, respectively.

CONCLUSIONS

This study demonstrates the importance of turbulence in assessing aortic haemodynamics in a patient with AVS. Neglecting the turbulent contribution to WSS could potentially result in a significant underestimation of the total WSS. Further work is warranted to extend the analysis to more AVS cases and patients with other aortic valve diseases.

Computational numerical analysis of different cannulation methods during cardiopulmonary bypass of type A aortic dissection model based on computational fluid dynamics

Background

The aim of the present study was to use a numerical simulation based on computational fluid dynamics (CFD) to analyze the difference of different cannulation methods on hemodynamics characteristic in a type A aortic dissection (TAAD) model.

Methods

A finite-element analysis based on the CFD model of a TAAD patient was used, and axillary artery cannulation (AAC), innominate artery cannulation (IAC), and femoral artery cannulation (FAC) were analyzed under different situations, including a cardiac output (CO) of 2.5 L/min and cardiopulmonary bypass (CPB) of 2.5 L/min (partial CPB before cross-clamping aorta, defined as condition A), and a CO of 0 L/min and CPB of 5 L/min (aortic cross-clamping phase, defined as condition B). The insertion of an 8-mm cannula into the different models was simulated. Hemodynamic characteristics, including wall shear stress, wall stress, blood flow, and velocity were analyzed.

Results

In condition A, the total flow of branches of the aortic arch was 2,009.5 mL/min (AAC), 1,855.47 mL/min (IAC), and 1,648.03 mL/min (FAC). All cannulation methods improved left renal blood perfusion. However, in relation to blood flow in the right renal artery, FAC showed the highest blood flow (105 mL/min). The results in condition B were similar to those of condition A. The velocity, shear stress, and stress of entry tear via AAC and IAC decreased in condition B compared with condition A. The velocity, shear stress, stress of tear via AAC was lower than that of IAC.

Conclusions

Different cannulation modes have an effect on the hemodynamic characteristic of the tear, but this effect is related to different states of CPB. AAC was found to superior to IAC, especially in reducing velocity, stress, and shear stress of site of tear. However, IAC and AAC are more conductive to blood supply than FAC in branch vessels of the aortic arch without being affected by the CPB state.

Influence of Distal Re-entry Tears on False Lumen Thrombosis After Thoracic Endovascular Aortic Repair in Type B Aortic Dissection Patients: A Computational Fluid Dynamics Simulation

PURPOSE

Distal re-entry tears play a significant role in false lumen (FL) thrombosis, which will strongly affect the postoperative long-term survival of patients with type B aortic dissection (TBAD) after thoracic endovascular aortic repair (TEVAR). This study aimed to investigate the influence of a peculiar morphological parameter of the residual re-entry tears in TBAD patients after TEVAR on long-term FL thrombosis using the computational fluid dynamics.

METHODS

Ideal population-based three-dimensional models of post-operative TBAD were established. Numerical simulation was performed to investigate the hemodynamic differences caused by different tear features, including the tear count, the maximum distance between tears, and the tear area.

RESULTS

Although the low relative residence time (RRT) area did not change significantly when the tear distance was fixed, the area of oscillatory shear index (OSI) > 0.45 and endothelial cell activation potential (ECAP) > 1.5 decreased significantly with the tear count and area increased and a dramatic increase in blood flow into the FL was also observed. When tear count and total area were fixed, for each 10-mm increase in the maximum distance between tears, the area of low RRT in the FL increased significantly, while the average pressure difference increased by 10.85%.

CONCLUSION

The different morphology of the re-entry tears had different effects on the thrombosis-related hemodynamic parameters in FL following TEVAR. and the number of re-entry tears was most crucial to the potential thrombosis in the post-TEVAR FL of TBAD patients.

Flow analysis of aortic dissection: comparison of inflow boundary conditions for computational models based on 4D PCMRI and Doppler ultrasound

Computational hemodynamics quantifying the flow environment is an important tool in understanding aortic dissection. In this study, various inflow boundaries were applied on a patient-specific model and compared to the individualized velocimetry. The results indicated that the computations generally overestimated the flow volume and underestimated the wall shear stress. By quantifying the accuracy of the simulation results, two inflow settings were suggested. One was individualized, the PCMRI-extracted 4D flow information, and the other was averaged by healthy data, the ultrasound-extracted averaged flow waveform with parabolic velocity profile. This study might contribute to improving the precise computation of aortic dissection hemodynamics.

Global Sensitivity Analysis for Patient-Specific Aortic Simulations: The Role of Geometry, Boundary Condition and Large Eddy Simulation Modeling Parameters

Numerical simulations for computational hemodynamics in clinical settings require a combination of many ingredients, mathematical models, solvers and patient-specific data. The sensitivity of the solutions to these factors may be critical, particularly when we have a partial or noisy knowledge of data. Uncertainty quantification is crucial to assess the reliability of the results. We present here an extensive sensitivity analysis in aortic flow simulations, to quantify the dependence of clinically relevant quantities to the patient-specific geometry and the inflow boundary conditions. Geometry and inflow conditions are generally believed to have a major impact on numerical simulations. We resort to a global sensitivity analysis, (i.e., not restricted to a linearization around a working point), based on polynomial chaos expansion (PCE) and the associated Sobol' indices. We regard the geometry and the inflow conditions as the realization of a parametric stochastic process. To construct a physically consistent stochastic process for the geometry, we use a set of longitudinal-in-time images of a patient with an abdominal aortic aneurysm (AAA) to parametrize geometrical variations. Aortic flow is highly disturbed during systole. This leads to high computational costs, even amplified in a sensitivity analysis -when many simulations are needed. To mitigate this, we consider here a large Eddy simulation (LES) model. Our model depends in particular on a user-defined parameter called filter radius. We borrowed the tools of the global sensitivity analysis to assess the sensitivity of the solution to this parameter too. The targeted quantities of interest (QoI) include: the total kinetic energy (TKE), the time-average wall shear stress (TAWSS), and the oscillatory shear index (OSI). The results show that these indexes are mostly sensitive to the geometry. Also, we find that the sensitivity may be different during different instants of the heartbeat and in different regions of the domain of interest. This analysis helps to assess the reliability of in silico tools for clinical applications.

A two-phase flow approach for modeling blood stasis and estimating the thrombosis potential of a ventricular assist device

Thrombosis and its related events have become a major concern during the development and optimization of ventricular assist devices (VADs, also called blood pumps), and limit their clinical use and economic benefits. Attempts have been made to model the thrombosis formation, considering hemodynamic and biochemical processes. However, the complexities and computational expenses are prohibitive. Blood stasis is one of the key factors which may lead to the formation of thrombosis and excessive thromboembolic risks for patients. This study proposed a novel approach for modeling blood stasis, based on a two-phase flow principle. The locations of blood residual can be tracked over time, so that regions of blood stasis can be identified. The blood stasis in an axial blood pump is simulated under various working conditions, the results agree well with the experimental results. In contrast, conventional hemodynamic metrics such as velocity, time-averaged wall shear stress (TAWSS), and relative residence time (RRT), were contradictory in judging risk of blood stasis and thrombosis, and inconsistent with experimental results. We also found that the pump operating at the designed rotational speed is less prone to blood stasis. The model provides an efficient and fast alternative for evaluating blood stasis and thrombosis potential in blood pumps, and will be a valuable addition to the tools to support the design and improvement of VADs.

A Combined In Vivo, In Vitro, In Silico Approach for Patient-Specific Haemodynamic Studies of Aortic Dissection

The optimal treatment of Type-B aortic dissection (AD) is still a subject of debate, with up to 50% of the cases developing late-term complications requiring invasive intervention. A better understanding of the patient-specific haemodynamic features of AD can provide useful insights on disease progression and support clinical management. In this work, a novel in vitro and in silico framework to perform personalised studies of AD, informed by non-invasive clinical data, is presented. A Type-B AD was investigated in silico using computational fluid dynamics (CFD) and in vitro by means of a state-of-the-art mock circulatory loop and particle image velocimetry (PIV). Both models not only reproduced the anatomical features of the patient, but also imposed physiologically-accurate and personalised boundary conditions. Experimental flow rate and pressure waveforms, as well as detailed velocity fields acquired via PIV, are extensively compared against numerical predictions at different locations in the aorta, showing excellent agreement. This work demonstrates how experimental and numerical tools can be developed in synergy to accurately reproduce patient-specific AD blood flow. The combined platform presented herein constitutes a powerful tool for advanced haemodynamic studies for a range of vascular conditions, allowing not only the validation of CFD models, but also clinical decision support, surgical planning as well as medical device innovation.

Geometric Analysis of Type B Aortic Dissections Shows Aortic Remodeling After Intervention Using Multilayer Stents

Recently, multilayer stents for type B aortic dissections (TBAD) have been proposed to decrease false lumen flow, increase and streamline true lumen flow, and retain branch vessel patency. We aimed to provide a protocol with standardized techniques to investigate aortic remodeling of TBAD by multilayer flow modulators (MFM) in static geometric and hemodynamic analyses. Combining existing literature and new insights, a standardized protocol was designed. Using pre- and postoperative CT scans, geometric models were constructed, lumen dimensions were calculated, computational fluid dynamics (CFD) models were composed, and velocity and pressures were calculated. Sixteen TBAD cases treated with MFM were included for analysis. For each case, aortic remodeling was analyzed using post-processing medical imaging software. After 3D models were created, geometrical anatomical measurements were performed, and meshes for finite element analysis were generated. MFM cases were compared pre- and postoperatively; true lumen volumes increased (p < 0.001), false lumen volumes decreased (p = 0.001), true lumen diameter at the plane of maximum compression (PMC) increased (p < 0.001), and false lumen index decreased (p = 0.008). True lumen flow was streamlined, and the overall fluid velocity and pressures decreased (p < 0.001 and p = 0.006, respectively). This protocol provided a standardized method to evaluate the effects of MFM treatments in TBAD on geometric analyses, PMC, and CFD outcomes.

The application of computational modeling for risk prediction in type B aortic dissection

OBJECTIVE

New tools are urgently needed to help with surgical decision-making in type B aortic dissection (TBAD) that is uncomplicated at the time of initial presentation. This narrative review aims to answer the clinical question, Can computational modeling be used to predict risk in acute and chronic Stanford TBAD?

METHODS

The review (PROSPERO 2018 CRD42018104472) focused on risk prediction in TBAD. A comprehensive search of the Ovid MEDLINE database, using terms related to computational modeling and aortic dissection, was conducted to find studies of any form published between 1998 and 2018. Cohort studies, case series, and case reports of adults (older than 18 years) with computed tomography or magnetic resonance imaging diagnosis of TBAD were included. Computational modeling was applied in all selected studies.

RESULTS

There were 37 studies about computational modeling of TBAD identified from the search, and the findings were synthesized into a narrative review. Computational modeling can produce numerically calculated values of stresses, pressures, and flow velocities that are difficult to measure in vivo. Hemodynamic parameters-high or low wall shear stress, high pressure gradient between lumens during the cardiac cycle, and high false lumen flow rate-have been linked to the pathogenesis of branch malperfusion and aneurysm formation by numerous studies. Considering the major outcomes of end-organ failure, aortic rupture, and stabilization and remodeling, hypotheses have been generated about inter-relationships of measurable parameters in computational models with observable anatomic and pathologic changes, resulting in specific clinical outcomes.

CONCLUSIONS

There is consistency in study findings about computational modeling in TBAD, although a limited number of patients have been analyzed using various techniques. The mechanistic patterns of association found in this narrative review should be investigated in larger cohort prospective studies to further refine our understanding. It highlights the importance of patient-specific computational hemodynamic parameters in clinical decision-making algorithms. The current challenge is to develop and to test a risk assessment method that can be used by clinicians for TBAD.

Highly integrated workflows for exploring cardiovascular conditions: Exemplars of precision medicine in Alzheimer's disease and aortic dissection

For precision medicine to be implemented through the lens of in silico technology, it is imperative that biophysical research workflows offer insight into treatments that are specific to a particular illness and to a particular subject. The boundaries of precision medicine can be extended using multiscale, biophysics-centred workflows that consider the fundamental underpinnings of the constituents of cells and tissues and their dynamic environments. Utilising numerical techniques that can capture the broad spectrum of biological flows within complex, deformable and permeable organs and tissues is of paramount importance when considering the core prerequisites of any state-of-the-art precision medicine pipeline. In this work, a succinct breakdown of two precision medicine pipelines developed within two Virtual Physiological Human (VPH) projects are given. The first workflow is targeted on the trajectory of Alzheimer's Disease, and caters for novel hypothesis testing through a multicompartmental poroelastic model which is integrated with a high throughput imaging workflow and subject-specific blood flow variability model. The second workflow gives rise to the patient specific exploration of Aortic Dissections via a multi-scale and compliant model, harnessing imaging, computational fluid-dynamics (CFD) and dynamic boundary conditions. Results relating to the first workflow include some core outputs of the multiporoelastic modelling framework, and the representation of peri-arterial swelling and peri-venous drainage solution fields. The latter solution fields were statistically analysed for a cohort of thirty-five subjects (stratified with respect to disease status, gender and activity level). The second workflow allowed for a better understanding of complex aortic dissection cases utilising both a rigid-wall model informed by minimal and clinically common datasets as well as a moving-wall model informed by rich datasets.

Patient-specific CFD modelling in the thoracic aorta with PC-MRI-based boundary conditions: A least-square three-element Windkessel approach

The increasing use of computational fluid dynamics for simulating blood flow in clinics demands the identification of appropriate patient-specific boundary conditions for the customization of the mathematical models. These conditions should ideally be retrieved from measurements. However, finite resolution of devices as well as other practical/ethical reasons prevent the construction of complete data sets necessary to make the mathematical problems well posed. Available data need to be completed by modelling assumptions, whose impact on the final solution has to be carefully addressed. Focusing on aortic vascular districts and related pathologies, we present here a method for efficiently and robustly prescribing phase contrast MRI-based patient-specific data as boundary conditions at the domain of interest. In particular, for the outlets, the basic idea is to obtain pressure conditions from an appropriate elaboration of available flow rates on the basis of a 3D/0D dimensionally heterogeneous modelling. The key point is that the parameters are obtained by a constrained optimization procedure. The rationale is that pressure conditions have a reduced impact on the numerical solution compared with velocity conditions, yielding a simulation framework less exposed to noise and inconsistency of the data, as well as to the arbitrariness of the underlying modelling assumptions. Numerical results confirm the reliability of the approach in comparison with other patient-specific approaches adopted in the literature.