Predicting false lumen thrombosis in patient-specific models of aortic dissection

Mechanisms of aortic dissection: From pathological changes to experimental and in silico models

Aortic dissection continues to be responsible for significant morbidity and mortality, although recent advances in medical data assimilation and in experimental and in silico models have improved our understanding of the initiation and progression of the accumulation of blood within the aortic wall. Hence, there remains a pressing necessity for innovative and enhanced models to more accurately characterize the associated pathological changes. Early on, experimental models were employed to uncover mechanisms in aortic dissection, such as hemodynamic changes and alterations in wall microstructure, and to assess the efficacy of medical implants. While experimental models were once the only option available, more recently they are also being used to validate in silico models. Based on an improved understanding of the deteriorated microstructure of the aortic wall, numerous multiscale material models have been proposed in recent decades to study the state of stress in dissected aortas, including the changes associated with damage and failure. Furthermore, when integrated with accessible patient-derived medical data, in silico models prove to be an invaluable tool for identifying correlations between hemodynamics, wall stresses, or thrombus formation in the deteriorated aortic wall. They are also advantageous for model-guided design of medical implants with the aim of evaluating the deployment and migration of implants in patients. Nonetheless, the utility of in silico models depends largely on patient-derived medical data, such as chosen boundary conditions or tissue properties. In this review article, our objective is to provide a thorough summary of medical data elucidating the pathological alterations associated with this disease. Concurrently, we aim to assess experimental models, as well as multiscale material and patient data-informed in silico models, that investigate various aspects of aortic dissection. In conclusion, we present a discourse on future perspectives, encompassing aspects of disease modeling, numerical challenges, and clinical applications, with a particular focus on aortic dissection. The aspiration is to inspire future studies, deepen our comprehension of the disease, and ultimately shape clinical care and treatment decisions.

Exploring hemodynamic mechanisms and re-intervention strategies for partial false lumen thrombosis in Stanford type B aortic dissection after thoracic aortic endovascular repair

OBJECTIVES

False lumen (FL) thrombosis status for Stanford type B aortic dissection (TBAD) after thoracic endovascular aortic repair (TEVAR) is critical for evaluating aortic remodeling and long-term prognosis. This study aimed to monitor the morphology evolution of partial FL thrombosis (PFLT) and its hemodynamic conditions through an innovative approach, providing a re-intervention strategy from both morphologic and hemodynamic perspectives.

METHODS

Three-dimensional geometries are extracted from a five-year follow-up of CTA images for TBAD after TEVAR. The morphology and hemodynamics of PFLT are comprehensively analyzed based on patient-specific reconstructions and computational fluid dynamics (CFD). The impact of various strategies treating risk factors of PFLT, including proximal entry closure, left renal artery stenting, or accessory renal artery embolism on hemodynamics is assessed.

RESULTS

The introduced morphologic approaches appropriately reflected the evolution of PFLT. Gradual dilation of FL (surface area from 82.63cm2 to 98.84cm2, volume from 45.12 mL to 63.40 mL, increase in distal tear (from 3.72 cm to 4.32 cm), and fluctuation of thrombosis-blood lumen boundary are observed. For further surgical preparation in the absence of unanimously recognized re-intervention indicators, velocity and wall shear stress distributions reveal different simulated re-interventions have distinctly suppressive effects on hemodynamic conditions within the PFLT, providing valuable insights for further surgical preparation.

CONCLUSIONS

The present study demonstrates a re-intervention strategy for PFLT in TBAD patients after TEVAR utilizing morphologic and hemodynamic analyses. Acknowledging the deterioration of PFLT may result in adverse long-term outcomes, this strategy might offer an alternative approach for clinical monitoring and management of related patients.

Effects of inlet boundary conditions on blood flow and thrombosis modelling in patients with chronic thromboembolic pulmonary hypertension

To investigate the impact of patient-specific boundary conditions (BC) on blood flow and thrombosis modelling for patients with chronic thromboembolic pulmonary hypertension (CTEPH), three types of BCs were utilized to construct CTEPH models based on computed tomography pulmonary angiography images. First BC type is the patient-specific velocity profiles at the main pulmonary artery using phase contrast MRI (PC-MRI). The other two simplified types are the pulsatile BC and steady BC, which are obtained by spatially and temporally averaging the PC-MRI BC. Hemodynamic features including helical density, time-averaged wall shear stress (TAWSS) and oscillatory shear index (OSI), and thrombosis were compared for the three types BCs. The results indicated that, compared to the MRI BC, steady BC overestimated helical density and TAWSS in the pulmonary arteries by approximately 63.1% and 60%, respectively. The impact of simplified pulsatile BC on TAWSS and OSI in most regions of the pulmonary arteries was negligible with differences within 5%. Regarding thrombosis, the area predicted under pulsatile BC was approximately 80% smaller than that under PC-MRI BC. In conclusion, compared to PC-MRI BC, steady inlet BC tend to overestimate hemodynamic parameters, while pulsatile inlet BC yield similar wall shear stress based on parameters in most regions of the pulmonary artery. Patient-specific PC-MRI inlet BC should be used for accurate predictions of helical flow pattern and thrombus formation.

Longitudinal investigation of aortic dissection in mice with computational fluid dynamics

Predicting late adverse events in aortic dissections is challenging. One commonly observed risk factor is partial thrombosis of the false lumen. In this study we investigated false lumen thrombus progression over 7 days in four mice with angiotensin II-induced aortic dissection. We performed computational fluid dynamic simulations with subject-specific boundary conditions from velocity and pressure measurements. We investigated endothelial cell activation potential, mean velocity, thrombus formation potential, and other hemodynamic factors. Our findings support the hypothesis that flow stagnation is the predominant hemodynamic factor driving a large thrombus ratio in false lumina, particularly those with a single fenestration.

Numerical investigation on the impact of different coronary aneurysms morphologies on thrombus formation and hemodynamics: a comparative study

Coronary artery aneurysms (CAAs) are morphologically classified as saccular and fusiform. There is still a great deal of clinical controversy as to which types of CAA are more likely to cause thrombosis. Therefore, the main objective of this study was to evaluate the trend of thrombus growth in CAAs with different morphologies and to assess the risk of possible long-term complications based on hemodynamic parameters. Utilizing computed tomography angiography (CTA) data from eight healthy coronary arteries, two distinct morphologies of coronary artery aneurysms (CAAs) were reconstructed. Distribution of four wall shear stress (WSS)-based indicators and three helicity indicators was analyzed in this study. Meanwhile, a thrombus growth model was introduced to analyze the thrombus formation in CAAs with different morphologies. The research results showed the distribution of most WSS indicators between saccular and fusiform CAAs was not statistically significant. However, due to the presence of a more pronounced helical flow pattern, irregular helical flow structure and longer time of flow stagnation in saccular CAAs during the cardiac cycle, the mean and maximum relative residence time (RRT) were significantly higher in saccular CAAs than in fusiform CAAs (P < 0.05). This may increase the risk of saccular coronary arteries leading to aneurysmal dilatation or even rupture. Although the two CAAs had similar rates of thrombosis, fusiform CAAs may more early cause obstruction of the main coronary flow channel where the aneurysm is located due to thrombosis growth. Thus, the risk of thrombosis in fusiform coronary aneurysms may warrant greater clinical concern.

Predictive Methods for Thrombus Formation in the Treatment of Aortic Dissection and Cerebral Aneurysms: A Comprehensive Review

Thrombus formation plays a crucial role in the clinical treatment of certain diseases. In conditions such as aortic dissection and cerebral aneurysm, complete thrombus occlusion in the affected region is desired to reduce blood flow into the false lumen or aneurysm sac, leading to a decrease in the tension exerted on the vascular wall and making it less likely to rupture. However, desired thrombosis sometimes fails to occur. Predicting thrombus formation can provide valuable information in such cases. This article offers a comprehensive review of conventional methods for predicting thrombus formation. In reviews conducted from the year 2000 to the present, the number of published related papers every five years has increased more than tenfold. We also found that the predictive methods can be classified into two categories: those based on the hemodynamic evaluation parameters and those based on hemodynamic and mathematical models that simulate the transport and reaction of blood components. Through our discussions, we identified several challenges that need to be resolved, including predictions based on patient-specific condition, model validation, multi-scale problems, the mechanisms of thrombus formation, and ensuring cost effectiveness. This review aims to guide researchers interested in exploring thrombus formation prediction within clinical treatments.

Haemodynamic effects of non-Newtonian fluid blood on the abdominal aorta before and after double tear rupture

OBJECTIVES

Intimal tears caused by aortic dissection can weaken the arterial wall and lead to aortic aneurysms. However, the effect of different tear states on the blood flow behaviour remains complex. This study uses a novel approach that combines numerical haemodynamic simulation with in vitro experiments to elucidate the effect of arterial dissection rupture on the complex blood flow state within the abdominal aneurysm and the endogenous causes of end-organ malperfusion.

MATERIALS AND METHODS

Based on the CT imaging data and clinical physiological parameters, the overall arterial models including aortic dissection and aneurysm with single tear and double tear were established, and the turbulence behaviours and haemodynamic characteristics of arterial dissection and aneurysm under different blood pressures were simulated by using non-Newtonian flow fluids with the pulsatile blood flow rate of the clinical patients as a cycle, and the results of the numerical simulation were verified by in vitro simulation experiments.

RESULTS

Hemodynamic simulations revealed that the aneurysm and single-tear false lumen generated a maximum pressure of 320.591 mmHg, 267 % over the 120 mmHg criterion. The pressure differential generates reflux, leading to a WSS of 2247.9 Pa at the TL inlet and blood flow velocities of up to 6.41 m/s inducing extend of the inlet. DTD Medium FL instantaneous WP above 120 mmHg Standard 151 % Additionally, there was 82.5 % higher flow in the right iliac aorta than in the left iliac aorta, which triggered malperfusion. Thrombus was accumulated distal to the tear and turbulence. These results are consistent with the findings of the in vitro experiments.

CONCLUSIONS

This study reveals the haemodynamic mechanisms by which aortic dissection induces aortic aneurysms to produce different risk states. This will contribute to in vitro simulation studies as a new fulcrum in the process of moving from numerical simulation to clinical trials.

Assessing the impact of tear direction in coronary artery dissection on thrombosis development: A hemodynamic computational study

OBJECTIVE

Iatrogenic coronary artery dissection is a complication of coronary intimal injury and dissection due to improper catheter manipulation. The impact of tear direction on the prognosis of coronary artery dissection (CAD) remains unclear. This study examines the hemodynamic effects of different tear directions (transverse and longitudinal) of CAD and evaluates the risk of thrombosis, rupture and further dilatation of CAD.

METHODS

Two types of CAD models (Type I: transverse tear, Type II: longitudinal tear) were reconstructed from the aorto-coronary CTA dataset of 8 healthy cases. Four WSS-based indicators were analyzed, including time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), and cross flow index (CFI). A thrombus growth model was also introduced to predict the trend of thrombus growth in CAD with two different tear directions.

RESULTS

For most of the WSS-based indicators, including TAWSS, RRT, and CFI, no statistically significant differences were observed across the CAD models with varying tear directions, except for OSI, where a significant difference was noted (p < 0.05). Meanwhile, in terms of thrombus growth, the thrombus growing at the tear of the Type I (transverse tear) CAD model extended into the true lumen earlier than that of the Type II (longitudinal tear) model.

CONCLUSIONS

Numerical simulations suggest that: (1) The CAD with transverse tear have a high risk of further tearing of the dissection at the distal end of the tear. (2) The CAD with longitudinal tear create a hemodynamic environment characterized by low TAWSS and high OSI in the false lumen, which may additionally increase the risk of vessel wall injury. (3) The CAD with transverse tear may have a higher risk of thrombosis and coronary obstruction and myocardial ischemia in the early phase of the dissection.

A new computational fluid dynamics based noninvasive assessment of portacaval pressure gradient

Accurate assessment of portacaval pressure gradient (PCG) in patients with portal hypertension (PH) is of great significance both for diagnosis and treatment. This study aims to develop a noninvasive method for assessing PCG in PH patients and evaluate its accuracy and effectiveness. This study recruited 37 PH patients treated with transjugular intrahepatic portosystemic shunt (TIPS). computed tomography angiography was used to create three dimension (3D) models of each patient before and after TIPS. Doppler ultrasound examinations were conducted to obtain the patient's portal vein flow (or splenic vein and superior mesenteric vein). Using computational fluid dynamics (CFD) simulation, the patient's pre-TIPS and post-TIPS PCG was determined by the 3D models and ultrasound measurements. The accuracy of these noninvasive results was then compared to clinical invasive measurements. The results showed a strong linear correlation between the PCG simulated by CFD and the clinical invasive measurements both before and after TIPS (R2 = 0.998, P < 0.001 and R2 = 0.959, P < 0.001). The evaluation accuracy of this noninvasive method reached 94 %, and the influence of ultrasound result errors on the numerical accuracy was found to be marginal if the error was less than 20 %. Furthermore, the information about the hemodynamic environment in the portal system was obtained by this numerical method. Spiral flow patterns were observed in the portal vein of some patients. In a conclusion, this study proposes a noninvasive numerical method for assessing PCG in PH patients before and after TIPS. This method can assist doctors in accurately diagnosing patients and selecting appropriate treatment plans. Additionally, it can be used to further investigate potential biomechanical causes of complications related to TIPS in the future.

A numerical study of the effect of thrombus breakdown on predicted thrombus formation and growth

Thrombosis is a complex biological process which involves many biochemical reactions and is influenced by blood flow. Various computational models have been developed to simulate natural thrombosis in diseases such as aortic dissection (AD), and device-induced thrombosis in blood-contacting biomedical devices. While most hemodynamics-based models consider the role of low shear stress in the initiation and growth of thrombus, they often ignore the effect of thrombus breakdown induced by elevated shear stress. In this study, a new shear stress-induced thrombus breakdown function is proposed and implemented in our previously published thrombosis model. The performance of the refined model is assessed by quantitative comparison with experimental data on thrombus formation in a backward-facing step geometry, and qualitative comparison with in vivo data obtained from an AD patient. Our results show that incorporating thrombus breakdown improves accuracy in predicted thrombus volume and captures the same pattern of thrombus evolution as measured experimentally and in vivo. In the backward-facing step geometry, thrombus breakdown impedes growth over the step and downstream, allowing a stable thrombus to be reached more quickly. Moreover, the predicted thrombus volume, height and length are in better agreement with the experimental measurements compared to the original model which does not consider thrombus breakdown. In the patient-specific AD, the refined model outperforms the original model in predicting the extent and location of thrombosis. In conclusion, the effect of thrombus breakdown is not negligible and should be included in computational models of thrombosis.

Impact on hemodynamics in carotid arteries with carotid webs at different locations: A Numerical Study Integrating Thrombus Growth Model

OBJECTIVE

Carotid webs (CWs), lesions in the carotid arteries, are gaining research interest due to the unclear link to ischemic stroke. Similarity to atherosclerosis in lesion location adds the complexity. The main purpose of study is to investigate the hemodynamic effects of CWs at different locations in carotid arteries.

METHODS

Three types of models with CWs were reconstructed from the CTA dataset of 8 healthy carotid arteries (Models A: CWs at the common carotid artery; B: at the origin of internal carotid artery; C: at the carotid sinus). Wall shear stress (WSS)-based parameters, including time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), and endothelial cell activation potential (ECAP) were analyzed. A thrombus growth model was also incorporated to assess long-term thrombus formation across different carotid webs locations.

RESULTS

Models A exhibited helical flow, whereas models B and C showed disturbed flow in the carotid sinus. Recirculation in Models A and B was mainly downstream of CWs, while Models C had both upstream and downstream recirculation. In addition, models A had higher overall TAWSS levels, with the smallest region of TAWSS < 0.4 pa (7.78 ± 8.35%). In contrast, Models C had larger areas with TAWSS < 0.4 pa, RRT > 100, and ECAP > 1.5, accounting for 14.18 ± 5.28%, 1.51 ± 1.17%, and 10.36 ± 4.10%, respectively. Noting that thrombus volume was highest in Models C (7.20 ± 3.95%).

CONCLUSIONS

Numerical simulations indicate that: 1) CWs have less hemodynamic impact when located in the CCA, but may increase flow resistance leading to distal branch ischemia; 2) CWs contribute to thrombus formation, primarily downstream in the common carotid artery and internal carotid artery origin, and both upstream and downstream in the sinus; 3) CWs at the origin of the ICA are more likely to result in disturbed blood flow patterns and thrombus aggregation than the other two locations, which may increase the risk of ischemic stroke in distal cerebral arteries.

Decoding thrombosis through code: a review of computational models

From the molecular level up to a blood vessel, thrombosis and hemostasis involves many interconnected biochemical and biophysical processes over a wide range of length and time scales. Computational modeling has gained eminence in offering insights into these processes beyond what can be obtained from in vitro or in vivo experiments, or clinical measurements. The multiscale and multiphysics nature of thrombosis has inspired a wide range of modeling approaches that aim to address how a thrombus forms and dismantles. Here, we review recent advances in computational modeling with a focus on platelet-based thrombosis. We attempt to summarize the diverse range of modeling efforts straddling the wide-spectrum of physical phenomena, length scales, and time scales; highlighting key advancements and insights from existing studies. Potential information gleaned from models is discussed, ranging from identification of thrombus-prone regions in patient-specific vasculature to modeling thrombus deformation and embolization in response to fluid forces. Furthermore, we highlight several limitations of current models, future directions in the field, and opportunities for clinical translation, to illustrate the state-of-the-art. There are a plethora of opportunity areas for which models can be expanded, ranging from topics of thromboinflammation to platelet production and clearance. Through successes demonstrated in existing studies described here, as well as continued advancements in computational methodologies and computer processing speeds and memory, in silico investigations in thrombosis are poised to bring about significant knowledge growth in the years to come.

Irregular anatomical features can alter hemodynamics in Takayasu arteritis

Objective

Takayasu arteritis (TA) is a difficult disease to deal with because there are neither reliable clinical signs, laboratory biomarkers, nor a single noninvasive imaging technique that can be used for early diagnosis and disease activity monitoring. Knowledge of aortic hemodynamics in TA is lacking. This study aimed to fill this gap by assessing hemodynamics in patients with TA using image-based computational fluid dynamics (CFD) simulations.

Methods

Eleven patients with TA were included in the present study. Patient-specific geometries were reconstructed from either clinical aortic computed tomography angiography or magnetic resonance angiography studies and coupled with physiological boundary conditions for CFD simulations. Key anatomical and hemodynamic parameters were compared with a control group consisting of 18 age- and sex-matched adults without TA who had healthy aortas.

Results

Compared with controls, patients with TA had significantly higher aortic velocities (0.9 m/s [0.7, 1.1 m/s] vs 0.6 m/s [0.5, 0.7 m/s]; P = .002), maximum time-averaged wall shear stress (14.2 Pa [9.8, 20.9 Pa] vs 8.0 Pa [6.2, 10.3 Pa]; P = .004), and maximum pressure drops between the ascending and descending aorta (36.9 mm Hg [29.0, 49.3 mm Hg] vs 28.5 mm Hg [25.8, 31.5 mm Hg]; P = .004). These significant hemodynamic alterations in patients with TA might result from abnormal anatomical features including smaller arch diameter (20.0 mm [13.8, 23.3 mm] vs 25.2 mm [23.3, 26.8 mm]; P = .003), supra-aortic branch diameters (21.9 mm [18.5, 24.6 mm] vs 25.7 mm [24.3, 28.3 mm]; P = .003) and descending aorta diameter (14.7 mm [12.2, 16.8 mm] vs 22.5 mm [19.8, 24.0 mm]; P < .001).

Conclusions

CFD analysis reveals hemodynamic changes in the aorta of patients with TA. The applicability of CFD technique coupled with standard imaging assessments in predicting disease progression of such patients will be explored in future studies. Future large cohort study with outcome correlation is also warranted.

Clinical Relevance

Based on patient-specific computational fluid dynamics simulations, the present retrospective study revealed significant difference in aortic hemodynamics between the patients with and without Takayasu arteritis (TA). To the best of our knowledge, this study is the first to evaluate hemodynamic conditions within TA, demonstrating the potential of computational flow modeling in capturing abnormal hemodynamic forces, such as high wall shear stress, resulted from irregular morphological changes. In the future, assessing the hemodynamic parameters within patients with TA during the prestenotic period, together with longitudinal computational fluid dynamics studies may allow better monitoring and management of TA.

A computational study of the effects of size, location, and number of tears on haemodynamics in surgically repaired type A aortic dissection

Objective

This study aimed to comprehensively examine the roles of size, location, and number of tears in the progression of surgically repaired type A aortic dissection (TAAD) by assessing haemodynamic changes through patient-specific computational fluid dynamic (CFD) simulations.

Methods

Two patient-specific TAAD geometries with replaced ascending aorta were reconstructed based upon computed 15 tomography (CT) scans, after which 10 hypothetical models (5 per patient) with different tear configurations were artificially created. CFD simulations were performed on all the models under physiologically realistic boundary conditions.

Results

Our simulation results showed that increasing either the size or number of the re-entry tears reduced the luminal pressure difference (LPD) and maximum time-averaged wall shear stress (TAWSS), as well as areas exposed to abnormally high or low TAWSS values. Models with a large re-entry tear outperformed the others by reducing the maximum LPD by 1.88 mmHg and 7.39 mmHg, for patients 1 and 2, respectively. Moreover, proximally located re-entry tears in the descending aorta were more effective at reducing LPD than distal re-entry tears.

Discussion

These computational results indicate that the presence of a relatively large re-entry tear in the proximal descending aorta might help stabilize post-surgery aortic growth. This finding has important implications for the management and risk stratification of surgically repaired TAAD patients. Nevertheless, further validation in a large patient cohort is needed.

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.

A Computational Study of Dynamic Obstruction in Type B Aortic Dissection

A serious complication in aortic dissection is dynamic obstruction of the true lumen (TL). Dynamic obstruction results in malperfusion, a blockage of blood flow to a vital organ. Clinical data reveal that increases in central blood pressure promote dynamic obstruction. However, the mechanisms by which high pressures result in TL collapse are underexplored and poorly understood. Here, we developed a computational model to investigate biomechanical and hemodynamical factors involved in Dynamic obstruction. We hypothesize that relatively small pressure gradient between TL and false lumen (FL) are sufficient to displace the flap and induce obstruction. An idealized fluid-structure interaction model of type B aortic dissection was created. Simulations were performed under mean cardiac output while inducing dynamic changes in blood pressure by altering FL outflow resistance. As FL resistance increased, central aortic pressure increased from 95.7 to 115.3 mmHg. Concurrent with blood pressure increase, flap motion was observed, resulting in TL collapse, consistent with clinical findings. The maximum pressure gradient between TL and FL over the course of the dynamic obstruction was 4.5 mmHg, consistent with our hypothesis. Furthermore, the final stage of dynamic obstruction was very sudden in nature, occurring over a short time (<1 s) in our simulation, consistent with the clinical understanding of this dramatic event. Simulations also revealed sudden drops in flow and pressure in the TL in response to the flap motion, consistent with first stages of malperfusion. To our knowledge, this study represents the first computational analysis of potential mechanisms driving dynamic obstruction in aortic dissection.

Monitoring of false lumen thrombosis in type B aortic dissection by impedance cardiography - A multiphysics simulation study

Aortic dissection is caused by a tear on the aortic wall that allows blood to flow through the wall layers. Usually, this tear involves the intimal and partly the medial layer of the aortic wall. As a result, a new false lumen develops besides the original aorta, denoted then as the true lumen. The local hemodynamic conditions such as flow disturbances, recirculations and low wall shear stress may cause thrombus formation and growth in the false lumen. Since the false lumen status is a significant predictor for late-dissection-related deaths, it is of great importance in the medical management of patients with aortic dissection. The hemodynamic changes in the aorta also alter the electrical conductivity of blood. Since the blood is much more conductive than other tissues in the body, such changes can be identified with non-invasive methods such as impedance cardiography. Therefore, in this study, the capability of impedance cardiography in monitoring thrombosis in the false lumen is studied by multiphysics simulations to assist clinicians in the medical management of patients under treatment. To tackle this problem, a 3D computational fluid dynamics simulation has been set up to model thrombosis in the false lumen and its impact on the blood flow-induced conductivity changes. The electrical conductivity changes of blood have been assigned as material properties of the blood-filled aorta in a 3D finite element electric simulation model to investigate the impact of conductivity changes on the measured impedance from the body's surface. The results show remarkable changes in the electrical conductivity distribution in the measurement region due to thrombosis in the false lumen, which significantly impacts the morphology of the impedance cardiogram. Thus, frequent monitoring of impedance cardiography signals may allow tracking the thrombus formation and growth in the false lumen.

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.

Imaging and biophysical modelling of thrombogenic mechanisms in atrial fibrillation and stroke

Atrial fibrillation (AF) underlies almost one third of all ischaemic strokes, with the left atrial appendage (LAA) identified as the primary thromboembolic source. Current stroke risk stratification approaches, such as the CHA₂DS₂-VASc score, rely mostly on clinical comorbidities, rather than thrombogenic mechanisms such as blood stasis, hypercoagulability and endothelial dysfunction-known as Virchow's triad. While detection of AF-related thrombi is possible using established cardiac imaging techniques, such as transoesophageal echocardiography, there is a growing need to reliably assess AF-patient thrombogenicity prior to thrombus formation. Over the past decade, cardiac imaging and image-based biophysical modelling have emerged as powerful tools for reproducing the mechanisms of thrombogenesis. Clinical imaging modalities such as cardiac computed tomography, magnetic resonance and echocardiographic techniques can measure blood flow velocities and identify LA fibrosis (an indicator of endothelial dysfunction), but imaging remains limited in its ability to assess blood coagulation dynamics. In-silico cardiac modelling tools-such as computational fluid dynamics for blood flow, reaction-diffusion-convection equations to mimic the coagulation cascade, and surrogate flow metrics associated with endothelial damage-have grown in prevalence and advanced mechanistic understanding of thrombogenesis. However, neither technique alone can fully elucidate thrombogenicity in AF. In future, combining cardiac imaging with in-silico modelling and integrating machine learning approaches for rapid results directly from imaging data will require development under a rigorous framework of verification and clinical validation, but may pave the way towards enhanced personalised stroke risk stratification in the growing population of AF patients. This Review will focus on the significant progress in these fields.

Morphological parameters affecting false lumen thrombosis following type B aortic dissection: a systematic study based on simulations of idealized models

Type B aortic dissection (TBAD) carries a high risk of complications, particularly with a partially thrombosed or patent false lumen (FL). Therefore, uncovering the risk factors leading to FL thrombosis is crucial to identify high-risk patients. Although studies have shown that morphological parameters of the dissected aorta are related to FL thrombosis, often conflicting results have been reported. We show that recent models of thrombus evolution in combination with sensitivity analysis methods can provide valuable insights into how combinations of morphological parameters affect the prospect of FL thrombosis. Based on clinical data, an idealized geometry of a TBAD is generated and parameterized. After implementing the thrombus model in computational fluid dynamics simulations, a global sensitivity analysis for selected morphological parameters is performed. We then introduce dimensionless morphological parameters to scale the results to individual patients. The sensitivity analysis demonstrates that the most sensitive parameters influencing FL thrombosis are the FL diameter and the size and location of intimal tears. A higher risk of partial thrombosis is observed when the FL diameter is larger than the true lumen diameter. Reducing the ratio of the distal to proximal tear size increases the risk of FL patency. In summary, these parameters play a dominant role in classifying morphologies into patent, partially thrombosed, and fully thrombosed FL. In this study, we point out the predictive role of morphological parameters for FL thrombosis in TBAD and show that the results are in good agreement with available clinical studies.

Effects of cardiac function alterations on the risk of postoperative thrombotic complications in patients receiving endovascular aortic repair

Introduction: Chronic heart disease (CHD) is a common comorbidity of patients receiving endovascular aneurysm repair (EVAR) for abdominal aortic aneurysms (AAA). The explicit relationship between ventricular systolic function and EVAR complication of thrombotic events is unknown. Methods: In this study, we proposed a three-dimensional numerical model coupled with the lumped-elements heart model, which is capable of simulating thrombus formation in diverse systolic functions. The relation of cardiac functions and the predicted risk of thrombus formation in the aorta and/or endograft of 4 patients who underwent EVAR was investigated. Relative risks for thrombus formation were identified using machine-learning algorithms. Results: The computational results demonstrate that thrombus tended to form on the interior side of the aorta arch and iliac branches, and cardiac function can affect blood flow field and affect thrombus formation, which is consistent with the four patients' post-operative imaging follow-up. We also found that RRT, OSI, TAWSS in thrombosis area are lower than whole average. In addition, we found that the thrombus formation has negative correlations with the maximum ventricular contractile force (r = -.281 ± .101) and positive correlations with the minimum ventricular contractile force (r = .238 ± .074), whereas the effect of heart rate (r = -.015 ± .121) on thrombus formation is not significant. Conclusion: In conclusion, changes in ventricular systolic function may alter the risk of thrombotic events after EVAR repair, which could provide insight into the selection of adjuvant therapy strategies for AAA patients with CHD.

Effects of non-Newtonian viscosity on arterial and venous flow and transport

It is well known that blood exhibits non-Newtonian viscosity, but it is generally modeled as a Newtonian fluid. However, in situations of low shear rate, the validity of the Newtonian assumption is questionable. In this study, we investigated differences between Newtonian and non-Newtonian hemodynamic metrics such as velocity, vorticity, and wall shear stress. In addition, we investigated cardiovascular transport using two different approaches, Eulerian mass transport and Lagrangian particle tracking. Non-Newtonian solutions revealed important differences in both hemodynamic and transport metrics relative to the Newtonian model. Most notably for the hemodynamic metrics, in-plane velocity and vorticity were consistently larger in the Newtonian approximation for both arterial and venous flows. Conversely, wall shear stresses were larger for the non-Newtonian case for both the arterial and venous models. Our results also indicate that for the Lagrangian metrics, the history of accumulated shear was consistently larger for both arterial and venous flows in the Newtonian approximation. Lastly, our results also suggest that the Newtonian model produces larger near wall and luminal mass transport values compared to the non-Newtonian model, likely due to the increased vorticity and recirculation. These findings demonstrate the importance of accounting for non-Newtonian behavior in cardiovascular flows exhibiting significant regions of low shear rate and recirculation.

Computational simulation of stent thrombosis induced by various degrees of stent malapposition

Percutaneous coronary intervention with stent implantation is one of the most commonly used approaches to treat coronary artery stenosis. Stent malapposition (SM) can increase the incidence of stent thrombosis, but the quantitative association between SM distance and stent thrombosis is poorly clarified. The objective of this study is to determine the biomechanical reaction mechanisms underlying stent thrombosis induced by SM and to quantify the effect of different SM severity grades on thrombosis. The thrombus simulation was performed in a continuous model based on the diffusion-convection response of blood substance transport. Simulated models included well-apposed stents and malapposed stents with various severities where the detachment distances ranged from 0 to 400 μm. The abnormal shear stress induced by SM was considered a critical contributor affecting stent thrombosis, which was dependent on changing SM distances in the simulation. The results illustrate that the proportion of thrombus volume was 1.88% at a SM distance of 75 μm (mild), 3.46% at 150 μm, and 3.93% at 400 μm (severe), but that a slight drop (3.18%) appeared at the detachment distance of 225 μm (intermediate). The results indicate that when the SM distance was less than 150 μm, the thrombus rose notably as the gap distance increased, whereas the progression of thrombogenicity weakened when it exceeded 150 μm. Therefore, more attention should be paid when SM is present at a gap distance of 150 μm. Moreover, when the SM length of stents are the same, thrombus tends to accumulate downstream towards the distal end of the stent as the SM distance increases.

MRI in CFD for chronic type B aortic dissection: Ready for prime time?

OBJECTIVES

Better tools are needed for risk assessment of Type B aortic dissection (TBAD) to determine optimal treatment for patients with uncomplicated disease. Magnetic resonance imaging (MRI) has the potential to inform computational fluid dynamics (CFD) simulations for TBAD by providing individualised quantification of haemodynamic parameters, for assessment of complication risks. This systematic review aims to present an overview of MRI applications for CFD studies of TBAD.

METHODS

Following PRISMA guidelines, a search in Medline, Embase, and the Scopus Library identified 136 potentially relevant articles. Studies were included if they used MRI to inform CFD simulation in TBAD.

RESULTS

There were 20 articles meeting the inclusion criteria. 19 studies used phase contrast MRI (PC-MRI) to provide data for CFD flow boundary conditions. In 12 studies, CFD haemodynamic parameter results were validated against PC-MRI. In eight studies, geometric models were developed from MR angiography. In three studies, aortic wall or intimal flap motion data were derived from PC/cine MRI.

CONCLUSIONS

MRI provides complementary patient-specific information in CFD haemodynamic studies for TBAD that can be used for personalised care. MRI provides structural, dynamic and flow data to inform CFD for pre-treatment planning, potentially advancing its integration into clinical decision-making. The use of MRI to inform CFD in TBAD surgical planning is promising, however further validation and larger cohort studies are required.

Model reduction of coagulation cascade based on genetic algorithm

Fibrin is an important product of the coagulation cascade, and plays an eminent role in platelet stabilization. Since coagulation cascade models typically involve the reaction kinetics of dozens of proteins, which will incur burdensome computational costs when coupled to blood flow in complex geometries, researchers often ignore this process when constructing thrombosis models. However, previous studies have shown that fundamental aspects of coagulation can be reproduced with simpler models, which motivated us to obtain a reduced-order model of fibrin generation through a systematic approach. Therefore, we introduced a semi-automatic framework to perform model-reduction of cascade reactions in this study, which consisted of two processes. Specifically, the retained protein species and cascade reactions were determined based on published studies and simulation results from the full cascade model, while the optimal reaction rates for the new cascade network were determined using a genetic algorithm. The framework has been applied to a 19-species coagulation model that triggers fibrin generation in internal fields via reactive boundaries, and a 10-species reduced-order model was obtained to reproduce the kinetics of fibrinogenesis in the full cascade model at different boundary tissue factor concentrations. This reduced-order model of fibrinogenesis would be valuable for thrombosis modeling that considers both the coagulation cascade and platelet activity. Furthermore, the framework proposed herein can also be applied to the reductions of other cascade reaction models.

Shear-driven modelling of thrombus formation in type B aortic dissection

Background: Type B aortic dissection (TBAD) is a dangerous pathological condition with a high mortality rate. TBAD is initiated by an intimal tear that allows blood to flow between the aortic wall layers, causing them to separate. As a result, alongside the original aorta (true lumen), a false lumen (FL) develops. TBAD compromises the whole cardiovascular system, in the worst case resulting in complete aortic rupture. Clinical studies have shown that dilation and rupture of the FL are related to the failure of the FL to thrombose. Complete FL thrombosis has been found to improve the clinical outcomes of patients with chronic TBAD and is the desired outcome of any treatment. Partial FL thrombosis has been associated with late dissection-related deaths and the requirement for re-intervention, thus the level of FL thrombosis is dominant in classifying the risk of TBAD patients. Therefore, it is important to investigate and understand under which conditions complete thrombosis of the FL occurs.

Method: Local FL hemodynamics play an essential role in thrombus formation and growth. In this study, we developed a simplified phenomenological model to predict FL thrombosis in TBAD under physiological flow conditions. Based on an existing shear-driven thrombosis model, a comprehensive model reduction study was performed to improve computational efficiency. The reduced model has been implemented in Ansys CFX and applied to a TBAD case following thoracic endovascular aortic repair (TEVAR) to test the model. Predicted thrombus formation based on post-TEVAR geometry at 1-month was compared to actual thrombus formation observed on a 3-year follow-up CT scan.

Results: The predicted FL status is in excellent agreement with the 3-year follow-up scan, both in terms of thrombus location and total volume, thus validating the new model. The computational cost of the new model is significantly lower than the previous thrombus model, with an approximate 65% reduction in computational time. Such improvement means the new model is a significant step towards clinical applicability.

Conclusion: The thrombosis model developed in this study is accurate and efficient at predicting FL thrombosis based on patient-specific data, and may assist clinicians in choosing individualized treatments in the future.

Patient specific approach to analysis of shear-induced platelet activation in haemodialysis arteriovenous fistula

Shear-induced platelet activation (SIPAct) is an important mechanism of thrombosis initiation under high blood flow. This mechanism relies on the interaction of platelets with the von Willebrand factor (VWF) capable of unfolding under high shear stress. High shear stress occurs in the arteriovenous fistula (AVF) commonly used for haemodialysis. A novel patient-specific approach for the modelling of SIPAct in the AVF was proposed. This enabled us to estimate the SIPAct level via computational fluid dynamics. The suggested approach was applied for the SIPAct analysis in AVF geometries reconstructed from medical images. The approach facilitates the determination of the SIPAct level dependence on both biomechanical (AVF flow rate) and biochemical factors (VWF multimer size). It was found that the dependence of the SIPAct level on the AVF flow rate can be approximated by a power law. The critical flow rate was a decreasing function of the VWF multimer size. Moreover, the critical AVF flow rate highly depended on patient-specific factors, e.g., the vessel geometry. This indicates that the approach may be adopted to elucidate patient-specific thrombosis risk factors in haemodialysis patients.

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.

Predicting the risk of postsplenectomy thrombosis in patients with portal hypertension using computational hemodynamics models: A proof-of-concept study

BACKGROUND

The high incidence of thrombosis in the portal venous system following splenectomy (a frequently adopted surgery for treating portal hypertension in patients with splenomegaly and hypersplenism) is a critical clinical issue. The aim of this study was to address whether quantification of postsplenectomy hemodynamics has potential value for assessing the risk of postsplenectomy thrombosis.

METHODS

Computational models were constructed for three portal hypertensive patients treated with splenectomy based on their preoperative clinical data to quantify hemodynamics in the portal venous system before and after splenectomy, respectively. Each patient was followed up for three or five months after surgery and examined with CT to screen potential thrombosis.

FINDINGS

The area ratio of wall regions exposed to low wall shear stress was small before splenectomy in all patients, which increased markedly after splenectomy and exhibited enlarged inter-patient differences. The largest area ratio of low wall shear stress and most severe flow stagnation after splenectomy were predicted for the patient suffering from postsplenectomy thrombosis, with the wall regions exposed to low wall shear stress corresponding well with the CT-detected distribution of thrombus. Further analyses revealed that postoperative hemodynamic characteristics were considerably influenced by the anatomorphological features of the portal venous system.

INTERPRETATION

Postoperative hemodynamic conditions in the portal venous system are highly patient-specific and have a potential link to postsplenectomy thrombosis, which indicates that patient-specific hemodynamic studies may serve as a complement to routine clinical assessments for refining risk stratification and postoperative patient management.

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.

Fluid-Structure Interaction Simulations of Repaired Type A Aortic Dissection: a Comprehensive Comparison With Rigid Wall Models

This study aimed to evaluate the effect of aortic wall compliance on intraluminal hemodynamics within surgically repaired type A aortic dissection (TAAD). Fully coupled two-way fluid-structure interaction (FSI) simulations were performed on two patient-specific post-surgery TAAD models reconstructed from computed tomography angiography images. Our FSI model incorporated prestress and different material properties for the aorta and graft. Computational results, including velocity, wall shear stress (WSS) and pressure difference between the true and false lumen, were compared between the FSI and rigid wall simulations. It was found that the FSI model predicted lower blood velocities and WSS along the dissected aorta. In particular, the area exposed to low time-averaged WSS ( ≤ 0.2   P a ) was increased from 21 cm2 (rigid) to 38 cm2 (FSI) in patient 1 and from 35 cm2 (rigid) to 144 cm2 (FSI) in patient 2. FSI models also produced more disturbed flow where much larger regions presented with higher turbulence intensity as compared to the rigid wall models. The effect of wall compliance on pressure difference between the true and false lumen was insignificant, with the maximum difference between FSI and rigid models being less than 0.25 mmHg for the two patient-specific models. Comparisons of simulation results for models with different Young's moduli revealed that a more compliant wall resulted in further reduction in velocity and WSS magnitudes because of increased displacements. This study demonstrated the importance of FSI simulation for accurate prediction of low WSS regions in surgically repaired TAAD, but a rigid wall computational fluid dynamics simulation would be sufficient for prediction of luminal pressure difference.

A New Mathematical Numerical Model to Evaluate the Risk of Thrombosis in Three Clinical Ventricular Assist Devices

(1) Background: Thrombosis is the main complication in patients supported with ventricular assist devices (VAD). Models that accurately predict the risk of thrombus formation in VADs are still lacking. When VADs are clinically assisted, their complex geometric configuration and high rotating speed inevitably generate complex flow fields and high shear stress. These non-physiological factors can damage blood cells and proteins, release coagulant factors and trigger thrombosis. In this study, a more accurate model for thrombus assessment was constructed by integrating parameters such as shear stress, residence time and coagulant factors, so as to accurately assess the probability of thrombosis in three clinical VADs. (2) Methods: A mathematical model was constructed to assess platelet activation and thrombosis within VADs. By solving the transport equation, the influence of various factors such as shear stress, residence time and coagulation factors on platelet activation was considered. The diffusion equation was applied to determine the role of activated platelets and substance deposition on thrombus formation. The momentum equation was introduced to describe the obstruction to blood flow when thrombus is formed, and finally a more comprehensive and accurate model for thrombus assessment in patients with VAD was obtained. Numerical simulations of three clinically VADs (CH-VAD, HVAD and HMII) were performed using this model. The simulation results were compared with experimental data on platelet activation caused by the three VADs. The simulated thrombogenic potential in different regions of MHII was compared with the frequency of thrombosis occurring in the regions in clinic. The regions of high thrombotic risk for HVAD and HMII observed in experiments were compared with the regions predicted by simulation. (3) Results: It was found that the percentage of activated platelets within the VAD obtained by solving the thrombosis model developed in this study was in high agreement with the experimental data (r² = 0.984), the likelihood of thrombosis in the regions of the simulation showed excellent correlation with the clinical statistics (r² = 0.994), and the regions of high thrombotic risk predicted by the simulation were consistent with the experimental results. Further study revealed that the three clinical VADs (CH-VAD, HVAD and HMII) were prone to thrombus formation in the inner side of the secondary flow passage, the clearance between cone and impeller, and the corner region of the inlet pipe, respectively. The risk of platelet activation and thrombus formation for the three VADs was low to high for CH-VAD, HVAD, and HM II, respectively. (4) Conclusions: In this study, a more comprehensive and accurate thrombosis model was constructed by combining parameters such as shear stress, residence time, and coagulation factors. Simulation results of thrombotic risk received with this model showed excellent correlation with experimental and clinical data. It is important for determining the degree of platelet activation in VAD and identifying regions prone to thrombus formation, as well as guiding the optimal design of VAD and clinical treatment.

Computational Prediction of Thrombosis in Food and Drug Administration's Benchmark Nozzle

Thrombosis seriously threatens human cardiovascular health and the safe operation of medical devices. The Food and Drug Administration’s (FDA) benchmark nozzle model was designed to include the typical structure of medical devices. However, the thrombosis in the FDA nozzle has yet not been investigated. The objective of this study is to predict the thrombus formation process in the idealized medical device by coupling computational fluid dynamics and a macroscopic hemodynamic-based thrombus model. We developed the hemodynamic-based thrombus model by considering the effect of platelet consumption. The thrombus model was quantitatively validated by referring to the latest thrombosis experiment, which was performed in a backward-facing step with human blood flow. The same setup was applied in the FDA nozzle to simulate the thrombus formation process. The thrombus shaped like a ring was firstly observed in the FDA benchmark nozzle. Subsequently, the accuracy of the shear-stress transport turbulence model was confirmed in different turbulent flow conditions. Five scenarios with different Reynolds numbers were carried out. We found that turbulence could change the shape of centrosymmetric thrombus to axisymmetric and high Reynolds number blood flow would delay or even prevent thrombosis. Overall, the present study reports the thrombosis process in the FDA benchmark nozzle using the numerical simulation method, and the primary findings may shed light on the effect of turbulence on thrombosis.

Functional Evaluation of Embedded Modular Single-Branched Stent Graft: Application to Type B Aortic Dissection With Aberrant Right Subclavian Artery

Background

Endovascular repair of type B aortic dissection (TBAD) with aberrant right subclavian artery (ARSA) is challenging due to anatomical complexity. The embedded modular single-branched stent graft (EMSBSG) could solve this problem. However, the hemodynamic efficacy of this innovative technique has not been fully assessed. This study aimed to propose morphometric and functional indicators to quantify the outcomes of EMSBSG in treating TBAD with ARSA.

Material and Methods

A patient who had TBAD with ARSA underwent EMSBSG implantation was admitted. Computational fluid dynamics (CFD) and three-dimensional structural analyses were conducted based on CTA datasets before the operation (Pre-1) and at 4 and 25 days after EMSBSG implantation (Post-1 and Post-2). Quantitative and qualitative functional analyses were conducted via pressure-, velocity- and wall shear stress (WSS) -based parameters, such as the luminal pressure difference (LPD), total energy loss, and flow distribution ratio. By precisely registering the aortas at the three time points, parameter variations in the EMSBSG region were also computed to investigate the prognostic improvement after EMSBSG implantation.

Results

The first balance point of LPD distally shifted to the abdominal aorta in Post-1 by a distance of 20.172 cm, and shifted out of the dissected region in Post-2, indicating positive pressure recovery post EMSBSG. The flow distribution ratios of all aortic arch branches increased after EMSBSG implantation. A positive normal deformation index in the EMSBSG region confirmed true lumen expansion; dominant AR_N (area ratio of negative value) of pressure and WSS-based parameters indicated an improved prognosis from Post-1 to Post-2.

Conclusions

The short-term results of EMSBSG in treating TBAD with ARSA proved to be promising, especially in EMSBSG region. Comprehensive evaluation could provide new insight into the therapy of TBAD with ARSA. Thus, it might guide the further management of complex aortic arch lesions.

Uncertainty quantification of a thrombosis model considering the clotting assay PFA-100®

Mathematical models of thrombosis are currently used to study clinical scenarios of pathological thrombus formation. As these models become more complex to predict thrombus formation dynamics high computational cost must be alleviated and inherent uncertainties must be assessed. Evaluating model uncertainties allows to increase the confidence in model predictions and identify avenues of improvement for both thrombosis modeling and anti-platelet therapies. In this work, an uncertainty quantification analysis of a multi-constituent thrombosis model is performed considering a common assay for platelet function (PFA-100®). The analysis is facilitated thanks to time-evolving polynomial chaos expansions used as a parametric surrogate for the full thrombosis model considering two quantities of interest; namely, thrombus volume and occlusion percentage. The surrogate is thoroughly validated and provides a straightforward access to a global sensitivity analysis via computation of Sobol' coefficients. Six out of 15 parameters linked to thrombus consitution, vWF activity, and platelet adhesion dynamics were found to be most influential in the simulation variability considering only individual effects; while parameter interactions are highlighted when considering the total Sobol' indices. The influential parameters are related to thrombus constitution, vWF activity, and platelet to platelet adhesion dynamics. The surrogate model allowed to predict realistic PFA-100® closure times of 300,000 virtual cases that followed the trends observed in clinical data. The current methodology could be used including common anti-platelet therapies to identify scenarios that preserve the hematological balance.

An integrated fluid-structure interaction and thrombosis model for type B aortic dissection

False lumen thrombosis (FLT) in type B aortic dissection has been associated with the progression of dissection and treatment outcome. Existing computational models mostly assume rigid wall behavior which ignores the effect of flap motion on flow and thrombus formation within the FL. In this study, we have combined a fully coupled fluid-structure interaction (FSI) approach with a shear-driven thrombosis model described by a series of convection-diffusion reaction equations. The integrated FSI-thrombosis model has been applied to an idealized dissection geometry to investigate the interaction between vessel wall motion and growing thrombus. Our simulation results show that wall compliance and flap motion can influence the progression of FLT. The main difference between the rigid and FSI models is the continuous development of vortices near the tears caused by drastic flap motion up to 4.45 mm. Flap-induced high shear stress and shear rates around tears help to transport activated platelets further to the neighboring region, thus speeding up thrombus formation during the accelerated phase in the FSI models. Reducing flap mobility by increasing the Young's modulus of the flap slows down the thrombus growth. Compared to the rigid model, the predicted thrombus volume is 25% larger using the FSI-thrombosis model with a relatively mobile flap. Furthermore, our FSI-thrombosis model can capture the gradual effect of thrombus growth on the flow field, leading to flow obstruction in the FL, increased blood viscosity and reduced flap motion. This model is a step closer toward simulating realistic thrombus growth in aortic dissection, by taking into account the effect of intimal flap and vessel wall motion.

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.

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.

Proximal true lumen collapse in a chronic type B aortic dissection patient: A case report

BACKGROUND

In the context of aortic dissection, increasing pressure within the newly formed false lumen can result in the progressive compression of the true aortic channel. However, true lumen collapse in chronic type B aortic dissection (cTBAD) patients is rare, with few clinical or experimental studies to date having explored the causes of such collapse.

CASE SUMMARY

In the present report, we describe a rare case of true-lumen collapse in an 83-year-old patient diagnosed with cTBAD, and we discuss potential therapeutic interventions for such cases. Following thoracic endovascular aortic repair (TEVAR), computed tomography angiography revealed satisfactory stent-graft positioning, no endoleakage, true lumen enlargement, thrombus formation in the false lumen, and slight enlargement of the true lumen distal to the stent-graft. Computational hemodynamic analyses indicated that the wall shear stress and pressure within the false lumen were significantly reduced following TEVAR.

CONCLUSION

TEVAR treatment of cTBAD patients suffering from proximal true lumen collapse can facilitate some degree of effective remodeling.

Delivery of Nitric Oxide in the Cardiovascular System: Implications for Clinical Diagnosis and Therapy

Nitric oxide (NO) is a key molecule in cardiovascular homeostasis and its abnormal delivery is highly associated with the occurrence and development of cardiovascular disease (CVD). The assessment and manipulation of NO delivery is crucial to the diagnosis and therapy of CVD, such as endothelial dysfunction, atherosclerotic progression, pulmonary hypertension, and cardiovascular manifestations of coronavirus (COVID-19). However, due to the low concentration and fast reaction characteristics of NO in the cardiovascular system, clinical applications centered on NO delivery are challenging. In this tutorial review, we first summarized the methods to estimate the in vivo NO delivery process, based on computational modeling and flow-mediated dilation, to assess endothelial function and vulnerability of atherosclerotic plaque. Then, emerging bioimaging technologies that have the potential to experimentally measure arterial NO concentration were discussed, including Raman spectroscopy and electrochemical sensors. In addition to diagnostic methods, therapies aimed at controlling NO delivery to regulate CVD were reviewed, including the NO release platform to treat endothelial dysfunction and atherosclerosis and inhaled NO therapy to treat pulmonary hypertension and COVID-19. Two potential methods to improve the effectiveness of existing NO therapy were also discussed, including the combination of NO release platform and computational modeling, and stem cell therapy, which currently remains at the laboratory stage but has clinical potential for the treatment of CVD.

Computational modeling of blood component transport related to coronary artery thrombosis in Kawasaki disease

Coronary artery thrombosis is the major risk associated with Kawasaki disease (KD). Long-term management of KD patients with persistent aneurysms requires a thrombotic risk assessment and clinical decisions regarding the administration of anticoagulation therapy. Computational fluid dynamics has demonstrated that abnormal KD coronary artery hemodynamics can be associated with thrombosis. However, the underlying mechanisms of clot formation are not yet fully understood. Here we present a new model incorporating data from patient-specific simulated velocity fields to track platelet activation and accumulation. We use a system of Reaction-Advection-Diffusion equations solved with a stabilized finite element method to describe the evolution of non-activated platelets and activated platelet concentrations [AP], local concentrations of adenosine diphosphate (ADP) and poly-phosphate (PolyP). The activation of platelets is modeled as a function of shear-rate exposure and local concentration of agonists. We compared the distribution of activated platelets in a healthy coronary case and six cases with coronary artery aneurysms caused by KD, including three with confirmed thrombosis. Results show spatial correlation between regions of higher concentration of activated platelets and the reported location of the clot, suggesting predictive capabilities of this model towards identifying regions at high risk for thrombosis. Also, the concentration levels of ADP and PolyP in cases with confirmed thrombosis are higher than the reported critical values associated with platelet aggregation (ADP) and activation of the intrinsic coagulation pathway (PolyP). These findings suggest the potential initiation of a coagulation pathway even in the absence of an extrinsic factor. Finally, computational simulations show that in regions of flow stagnation, biochemical activation, as a result of local agonist concentration, is dominant. Identifying the leading factors to a pro-coagulant environment in each case—mechanical or biochemical—could help define improved strategies for thrombosis prevention tailored for each patient.

Computational studies aiming to model thrombosis often rely on an arterial wall injury. Collagen and other extracellular matrix components are exposed to the bloodstream, which facilitates platelet adhesion to the wall and subsequent clot formation. However, these models are not adequate to explain thrombosis in other settings where even in the absence of a focal lesion, clots may still form under certain flow conditions. Coronary artery aneurysm thrombosis following KD is an example of the need to understand the mechanisms of thrombus initiation in the absence of an extrinsic factor. This study provides a new framework to investigate thrombus initiation in KD from a patient-specific perspective, which integrates fluid mechanics and biochemistry and which could help quantify the pro-coagulant environment induced by the aneurysm and become a predictive tool. The work presented here has broad relevance to other clinical situations where flow stagnation and transport are driving factors in thrombus formation.

Outflow Through Aortic Side Branches Drives False Lumen Patency in Type B Aortic Dissection

Objective: Thoracic endovascular aortic repair (TEVAR) for type B aortic dissection (TBAD) aims to induce false lumen (FL) thrombosis by sealing intimal tears between the true (TL) and the FL, and blocking the inflow into the FL. Incomplete thrombosis of the FL is correlated with poor clinical outcome. We hypothesize that the number of major and minor branches arising from the FL affects FL patency and may negatively influence TEVAR induced FL thrombosis.

Methods: Computed tomography (CT)-scans from 89 patients diagnosed with TBAD [best medical treatment (BMT) n = 52, TEVAR n = 37] from two high-volume vascular surgery centers were analyzed retrospectively. Analysis included evaluation of the FL patency status, the number, location and size of intimal tears, and the presence of minor and major side branches originating from the FL. Multiple regression analysis was conducted to evaluate obtained parameters as predictors for FL thrombosis status.

Results: In univariate analysis, the strongest correlation for FL patency was found for the number of major (R = 0.79) and minor (R = 0.86) side branches originating from the FL. When applying a multiple linear regression model, the number of major (normalized beta 0.37; P < 0.001) and minor (normalized beta 0.41; P < 0.01) side branches arising from the FL were valid predictors for the axial length of the patent and non-patent FL, and additionally determined the length of the patent FL at 12-month follow-up in patients that underwent TEVAR.

Conclusions: Our data suggest that the number of minor side branches that originate from the FL in TBAD is an important determinant of FL patency, to a greater degree than previously assumed.

Risk factors for incomplete thrombosis in false lumen in sub-acute type B aortic dissection post-TEVAR

There is scarce information about the risk factors for incomplete false lumen thrombosis (FLT) among type B aortic dissection (AD) patients, particularly in the sub-acute phase following thoracic endovascular aortic repair (TEVAR). We enrolled consecutive sub-acute type B AD patients at Xinqiao Hospital (Chongqing, China) from May 2010 to December 2019. Patients with severe heart failure, cancer, and myocardial infarction were excluded. The postoperative computed tomography angiography (CTA) data were extracted from the most recent follow-up aortic CTA. Multivariate logistic regressions were applied to identify the association between FLT and clinical or imaging factors. Fifty-five subjects were enrolled in our study. Twelve participants showed complete FLT, and 2 of these died during the follow-up, while 8 patients died in incomplete FLT group. In the multivariate analysis, maximum abdominal aorta diameter (OR 1.20, 95% CI 1.016-1.417 p = 0.032) and the number of branches arising from the false lumen (FL) (OR 15.062, 95% 1.681-134.982 p = 0.015) were significantly associated with incomplete FLT. The C-statistics was 0.873 (95% CI 0.773-0.972) for the model. The FL diameter (p < 0.001) was significantly shorter following TEVAR, while the true lumen diameter (p < 0.001) and maximum abdominal aorta diameter (p = 0.011) were larger after the aortic repair. There were 21.8% of sub-acute type B AD patients presented complete FLT post-TEVAR. Maximum abdominal aorta diameter and the number of branches arising from the FL were independent risk factors for incomplete FLT. The true lumen diameter, maximum abdominal aorta diameter, and FL diameter changed notably following TEVAR.

The effect of beta-blockers on hemodynamic parameters in patient-specific blood flow simulations of type-B aortic dissection: a virtual study

Aortic dissection (AD) is one of the fatal and complex conditions. Since there is a lack of a specific treatment guideline for type-B AD, a better understanding of patient-specific hemodynamics and therapy outcomes can potentially control the progression of the disease and aid in the clinical decision-making process. In this work, a patient-specific geometry of type-B AD is reconstructed from computed tomography images, and a numerical simulation using personalised computational fluid dynamics (CFD) with three-element Windkessel model boundary condition at each outlet is implemented. According to the physiological response of beta-blockers to the reduction of left ventricular contractions, three case studies with different heart rates are created. Several hemodynamic features, including time-averaged wall shear stress (TAWSS), highly oscillatory, low magnitude shear (HOLMES), and flow pattern are investigated and compared between each case. Results show that decreasing TAWSS, which is caused by the reduction of the velocity gradient, prevents vessel wall at entry tear from rupture. Additionally, with the increase in HOLMES value at distal false lumen, calcification and plaque formation in the moderate and regular-heart rate cases are successfully controlled. This work demonstrates how CFD methods with non-invasive hemodynamic metrics can be developed to predict the hemodynamic changes before medication or other invasive operations. These consequences can be a powerful framework for clinicians and surgical communities to improve their diagnostic and pre-procedural planning.

An integrated fluid-chemical model toward modeling the thrombus formation in an idealized model of aortic dissection

Type B aortic dissection is a major aortic catastrophe that can be acutely complicated by rapid expansion, rupture, and malperfusion syndromes. The separation of the intima from aortic walls will form a second blood-filled lumen defined as "false lumen (FL)", where the thrombus is more likely to form due to the local stasis hemodynamic conditions. Complete thrombosis of FL is associated with a beneficial outcome while patency and partial thrombosis will lead to later complications. However, the thrombosis mechanism is still unclear and little is known about the impact of chemical species transported by blood flow on this process. The proteins involved in the coagulation cascade (CC) may play an important role in the process of thrombosis, especially in the activation and stabilization of platelets. Based on this hypothesis, a reduced-order fluid-chemical model was established to simulate CC in an aortic dissection phantom with two tears. A high level of fibrin is continuously observed at the top of the FL and some time-varying areas between two tears, indicating a high likelihood of thrombus formation there. This finding is consistent with the clinical observation. The time evolution of coagulation factors is greatly affected by local hemodynamics, especially in the high disturbance zone where the evolution has characteristics of periodic changes consistent with the flow field. The ability of the proposed model to reproduce the CC response provides a potential application to integrate with a model that can simulate platelet activities, forming a biochemical-based model which would help unveil the mechanisms of thrombosis in FL and the clinical decision of appropriate treatment.

Sequential numerical simulation of vascular remodeling and thrombosis in unconventional hybrid repair of ruptured middle aortic syndrome

Unconventional surgical procedures may be utilized in treating complicated middle aortic syndrome (MAS), the outcome and prognosis of which remain largely undetermined due to limited numbers and significant heterogeneity of this population. Using computational fluid dynamics (CFD) analysis, this study aimed to assess the dynamic changes of postoperative aortic flow in seeking to unveil the relationship between hemodynamics and vascular remodeling and thrombotic events. One patient with middle aortic syndrome complicated with aortic rupture was treated with hybrid repair of extra-anatomic bypass and fenestrated endovascular aortic repair. The patient was followed-up for 8 months by computational tomography angiography and Doppler ultrasound. Thoracoabdominal aortic blood flow and locations with ongoing thrombosis at 1, 3, and 6 months postoperatively were simulated and analyzed. Remodeling processes, including low wall shear-mediated constrictive remodeling of non-stented aorta, neointimal hyperplasia at suture lines, and minimal thrombosis at various locations, were evident. Meanwhile, abdominal blood flow was tri-phasic at 1 month after surgery, and was reversed and stabilized at 6 months. The distribution of newly formed thrombus vary at different follow-up stages, which were in line with the numerical simulation of thrombosis from different postoperative time points. CFD-based sequential monitoring is of promising value in capturing dynamic changes of vascular outcome.

The effect of coarctation degrees on wall shear stress indices

Coarctation of the aorta (CoA) is a congenital tightening of the proximal descending aorta. Flow quantification can be immensely valuable for an early and accurate diagnosis. However, there is a lack of appropriate diagnostic approaches for a variety of cardiovascular diseases, such as CoA. An accurate understanding of the disease depends on measurements of the global haemodynamics (criteria for heart function) and also the local haemodynamics (detailed data on the dynamics of blood flow). Playing a significant role in clinical processes, wall shear stress (WSS) cannot be measured clinically; thus, computation tools are needed to give an insight into this crucial haemodynamic parameter. In the present study, in order to enable the progress of non-invasive approaches that quantify global and local haemodynamics for different CoA severities, innovative computational blueprint simulations that include fluid-solid interaction models are developed. Since there is no clear approach for managing the CoA regarding its severity, this study proposes the use of WSS indices and pressure gradient to better establish a framework for treatment procedures in CoA patients with different severities. This provides a platform for improving CoA therapy on a patient-specific level, in which physicians can perform treatment methods based on WSS indices on top of using a mere experience. Results show how severe CoA affects the aorta in comparison to the milder cases, which can give the medical community valuable information before and after any intervention.

The characteristics of distal tears affect false lumen thrombosis rate after thoracic endovascular aortic repair for acute type B dissection

OBJECTIVES

A low false lumen thrombosis rate (FLTR) is common in patients with type B aortic dissection after they have undergone thoracic endovascular aortic repair, which indicates a poor long-term prognosis. This study aimed to establish a quantitative linear regression model to predict false lumen (FL) thrombosis accurately using morphological parameters.

METHODS

In this retrospective study, we included 59 patients diagnosed with acute type B aortic dissection between 2014 and 2017. Morphological parameters were measured. Univariable and multivariable linear regression analyses were performed, and a linear regression model relating FLTR with the number of re-entry tears was proposed. Ten patients were further chosen to validate the linear relationship, and idealized aortic dissection models were adopted for haemodynamic analysis.

RESULTS

Only the total area and number of re-entry tears were negatively correlated with FL thrombosis (P < 0.001). Moreover, based on the univariable regression, the number of re-entry tears played a more crucial role in FLTR (R2 = 0.509 vs R2 = 0.298), and the linear relationship model was created as follows: thrombosis rate (%) = -11.25 × distal tear number + 105.24. This model was perfectly matched in 10 patients (concordance correlation coefficient = 0.880, P = 0.947). Moreover, when the total area of re-entry tears was constant, the net blood flow increased rapidly with an increase in the tear count.

CONCLUSIONS

The number of re-entry tears could be a crucial related factor of FL thrombosis; the larger the number of re-entry tears is, the lower the possibility of thrombosis is.

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.