These findings provide significant implications for the enhancement of iliac vein stent implantation strategies and stent design. The prevalent use of stents for treating Iliac Vein Compression Syndrome (IVCS) has shown efficacy, yet the associated clinical adverse events, including stent restenosis and postoperative thrombosis, are significant concerns. Up to now, the mechanism how the stent implantation induces the restenosis and DVT is still unclear. Our study hypothesizes that these adverse outcomes arise from altered blood flow dynamics following stent implantation. Employing computational modeling and medical imaging, we simulated IVCS after various stenting procedures to assess their impact on venous blood flow characteristics, including wall shear stress (WSS), residence time (RRT), and oscillatory shear index (OSI). Our findings reveal that a stent protruding into the vena cava impedes blood circulation, with increased protrusion exacerbating this obstruction. This is particularly evident at the vein bifurcation, where low WSS and elevated OSI and RRT are observed. Moreover, a higher stent strut density further obstructs blood flow, deteriorating the hemodynamic environment. Consequently, stent protrusion into the vena cava can enhance the likelihood of adverse post-surgical events. These insights have profound implications for optimizing iliac vein stent implantation techniques and stent design.

PurposeFluid-structure interaction (FSI) models are more commonly applied in medical research as computational power is increasing. However, understanding the accuracy of FSI models is crucial, especially in the context of heart valve disease in patient-specific models. Therefore, this study aimed to create a multi-modal benchmarking data set for cardiac-inspired FSI models, based on clinically important parameters, such as the pressure, velocity, and valve opening, with an in vitro phantom setup.MethodAn in vitro setup was developed with a 3D-printed phantom mimicking the left heart, including a deforming mitral valve. A range of pulsatile flows were created with a computer-controlled motor-and-pump setup. Catheter pressure measurements, magnetic resonance imaging (MRI), and echocardiography (Echo) imaging were used to measure pressure and velocity in the domain. Furthermore, the valve opening was quantified based on cine MRI and Echo images.ResultThe experimental setup, with 0.5% cycle-to-cycle variation, was successfully built and six different flow cases were investigated. Higher velocity through the mitral valve was observed for increased cardiac output. The pressure difference across the valve also followed this trend. The flow in the phantom was qualitatively assessed by the velocity profile in the ventricle and by streamlines obtained from 4D phase-contrast MRI.ConclusionA multi-modal set of data for validation of FSI models has been created, based on parameters relevant for diagnosis of heart valve disease. All data is publicly available for future development of computational heart valve models.

To enhance the performance of machine learning (ML) models for the post-embolization recanalization of cerebral aneurysms, we evaluated the impact of hemodynamic feature derivation and selection method on six ML algorithms. We utilized computational fluid dynamics (CFD) to simulate hemodynamics in 66 cerebral aneurysms from 65 patients, including 57 stable and nine recanalized aneurysms. We derived a total of 107 features for each aneurysm, encompassing four clinical features, 12 morphological features, and 91 hemodynamic features. To investigate the influence of feature derivation and selection methods on the ML models, we employed two derivation methods, simplified and fully derived, in combination with four selection methods: all features, statistically significant analysis, stepwise multivariate logistic regression analysis (stepwise-LR), and recursive feature elimination (RFE). Model performance was assessed using the area under the receiver operating characteristic curve (AUROC) and precision-recall curve (AUPRC) on both the training and testing datasets. The AUROC values on the testing dataset exhibited a wide-ranging spectrum, spanning from 0.373 to 0.863. Fully derived features and the RFE selection method demonstrated superior performance in intra-model comparisons. The multi-layer perceptron (MLP) model, trained with RFE-selected fully derived features, achieved the best performance on the testing dataset, with an AUROC value of 0.863 (95% CI: 0.684- 1.000). Our study demonstrated the importance of feature derivation and selection in determining the performance of ML models. This enabled the development of accurate decision-making models without the need to invade the patient.

One in four deaths worldwide is due to thromboembolic disease; that is, one in four people die from blood clots first forming and then breaking off or embolizing. Once broken off, clots travel downstream, where they occlude vital blood vessels such as those of the brain, heart, or lungs, leading to strokes, heart attacks, or pulmonary embolisms, respectively. Despite clots' obvious importance, much remains to be understood about clotting and clot embolization. In our work, we take a first step toward untangling the mystery behind clot embolization and try to answer the simple question: "What makes blood clots break off?" To this end, we conducted experimentally-informed, back-of-the-envelope computations combining fracture mechanics and phase-field modeling. We also focused on deep venous clots as our model problem. Here, we show that of the three general forces that act on venous blood clots-shear stress, blood pressure, and wall stretch-induced interfacial forces-the latter may be a critical embolization force in occlusive and non-occlusive clots, while blood pressure appears to play a determinant role only for occlusive clots. Contrary to intuition and prior reports, shear stress, even when severely elevated, appears unlikely to cause embolization. This first approach to understanding the source of blood clot bulk fracture may be a critical starting point for understanding blood clot embolization. We hope to inspire future work that will build on ours and overcome the limitations of theseback-of-the-envelope computations.

PurposeThis study investigated the implications of inserting a flexible annuloplasty ring after reconstructing the entire mitral valve in a porcine model using a previously investigated tube graft design made of 2-ply small intestinal submucosa extracellular matrix (CorMatrix®).MethodsAn acute model with eight 80-kg pigs, each acting as its own control, was used. The entire mitral valve was reconstructed with a 2-ply small intestinal submucosa extracellular matrix tube graft (CorMatrix®). Subsequently, a Simulus® flexible ring was inserted. The characterization was based on mitral annular geometry and valvular dynamics with sonomicrometry and echocardiography.ResultsAfter adding the ring annuloplasty, the in-plane annular dynamics were more constant throughout the cardiac cycle compared to the reconstruction alone. However, the commissure–commissure distance was statistically significantly decreased [35.0 ± 3.4 mm vs. 27.4 ± 1.9 mm, P < 0.001, diff = − 7.6 mm, 95% CI, − 9.8 to (−5.4) mm] after ring insertion, changing the physiological annular D-shape into a circular shape which created folds at the coaptation zone resulting in a central regurgitant jet on color Doppler.ConclusionWe successfully reconstructed the entire mitral valve using 2-ply small intestinal submucosal extracellular matrix (CorMatrix®) combined with a flexible annuloplasty. The annuloplasty reduced the unphysiological systolic widening previously found with this reconstructive technique. However, the Simulus flex ring changed the physiological annular D-shape into a circular shape and hindered a correct unfolding of the leaflets. Thus, we do not recommend a flexible ring in conjunction with this reconstructive technique; further investigations are needed to discover a more suitable remodelling annuloplasty.

Intravascular endoscopy can aid in the diagnosis of coronary atherosclerosis by providing direct color images of coronary plaques. The procedure requires a blood-free optical path between the catheter and plaque, and achieving clearance safely remains an engineering challenge. In this study, we investigate the hemodynamics of saline flushing in partially occluded coronary arteries to advance the development of intravascular forward-imaging catheters that do not require balloon occlusion. In-vitro experiments and CFD simulations are used to quantify the influence of plaque size, catheter stand-off distance, saline injection flowrate, and injection orientation on the time required to achieve blood clearance. Experiments and simulation of saline injection from a dual-lumen catheter demonstrated that flushing times increase both as injection flow rate (Reynolds number) decreases and as the catheter moves distally away from the plaque. CFD simulations demonstrated that successful flushing was achieved regardless of lumen axial orientation in a 95% occluded artery. Flushing time was also found to increase as plaque size decreases for a set injection flowrate, and a lower limit for injection flowrate was found to exist for each plaques size, below which clearance was not achieved. For the three occlusion sizes investigated (90, 95, 97% by area), successful occlusion was achieved in less than 1.2s. Investigation of the pressure fields developed during injection, highlight that rapid clearance can be achieved while keeping the arterial overpressure to < 1mmHg. A dual lumen saline injection catheter was shown to produce clearance safely and effectively in models of partially occluded coronary arteries. Clearance was achieved across a range of engineering and clinical parameters without the use of a balloon occlusion, providing development guideposts for a fluid injection system in forward-imaging coronary endoscopes.

Piezoelectric energy harvesters (PEH) for cardiac pacemakers typically use animal models to assess the performance of the PEH. However, if considering multiple designs, the use of animal models and prototyping increases costs and time. To reduce the use of animal models in research for pacemaker energy harvesting applications, this study investigates the motion of a pacemaker lead wire (PLW) in vivo using fluoroscopy imaging to quantify the position and displacements as a function of time, such that the data can be used in computer simulations. The proposed technique uses fluoroscopy imaging video data of a dual chamber pacemaker implanted in a patient, and image processing allows for the motion of the PLW captured. The motion is discretized into nodes for ease of implementation in finite element software. FEA simulation is presented using a piezoelectric energy harvester design integrated in the lead wire, and the energy output is predicted by finite element computer simulation. A 2-dimensional analysis is conducted with the fluoroscopy imaging video data to characterize the PLW motion and results show close agreement with literature values. Simulations with an energy harvesting circuit using the nodal position and displacement data shows that a PEH integrated in the PLW can generate a direct current voltage of 1.12V and power output of 0.125μW, potentially extending the battery life of pacemakers by 0.75-1years. The results suggest that fluoroscopy imaging data can be effective in evaluating PEH designs rather than using animal models, saving time and costs.

Patients presenting with coarctation of the aorta (CoA) may also suffer from co-existing transverse arch hypoplasia (TAH). Depending on the risks associated with the surgery and the severity of TAH, clinicians may decide to repair only CoA, and monitor the TAH to see if it improves as the patient grows. While acutely successful, eventually hemodynamics may become suboptimal if TAH is left untreated. The objective of this work aims to develop a patient-specific surgical planning framework for predicting and assessing postoperative outcomes of simple CoA repair and comprehensive repair of CoA and TAH. The surgical planning framework consisted of virtual clamp placement, stenosis resection, and design and optimization of patient-specific aortic grafts that involved geometrical modeling of the graft and computational fluid dynamics (CFD) simulation for evaluating various surgical plans. Time-dependent CFD simulations were performed using Windkessel boundary conditions at the outlets that were obtained from patient-specific non-invasive pressure and flow data to predict hemodynamics before and after the virtual repairs. We applied the proposed framework to investigate optimal repairs for six patients (n = 6) diagnosed with both CoA and TAH. Design optimization was performed by creating a combination of a tubular graft and a waterslide patch to reconstruct the aortic arch. The surfaces of the designed graft were parameterized to optimize the shape. Peak systolic pressure drop (PSPD) and time-averaged wall shear stress (TAWSS) were used as performance metrics to evaluate surgical outcomes of various graft designs and implantation. The average PSPD improvements were 28% and 44% after the isolated CoA repair and comprehensive repair, respectively. Maximum values of TAWSS were decreased by 60% after CoA repair and further improved by 22% after the comprehensive repair. The oscillatory shear index was calculated and the values were confirmed to be in the normal range after the repairs. The results showed that the comprehensive repair outperforms the simple CoA repair and may be more advantageous in the long term in some patients. We demonstrated that the surgical planning and patient-specific flow simulations could potentially affect the selection and outcomes of aorta repairs.

Bio-impedance analysis (BIA) has been widely investigated for hemodynamic monitoring. However, previous works rarely modelled two synchronously pulsatile arteries (representing the radial and ulnar arteries) in the wrist/forearm model. This work aims to clarify and quantify the influences of two pulsatile arteries on BIA. First, two blood-filled arteries were structured in a 3D wrist segment using the finite element method (FEM). Afterwards, an easy-to-produce two-arteries artificial wrist was fabricated with two components: gelatine-based surrounding tissue phantom and saline blood phantom. A syringe driver was utilised to constrict the arteries, and the impedance signals were measured using a Multi-frequency Impedance Analyser (MFIA). Both simulation and experimental results demonstrated the non-negligible influences of the ulnar artery on the overall BIA, inducing unwanted resistance changes to the acquired signals from the radial artery. The phantom experiments revealed the summation of the individual resistance changes caused by a single pulsatile artery was approximately equal to the measured resistance change caused by two synchronously pulsatile arteries, confirming the measured impedance signal at the wrist contains the pulsatile information from both arteries. This work is the first simulation and phantom investigation into two synchronously pulsatile arteries under BIA in the distal forearm, providing a better insight and understanding in the morphology of measured impedance signals. Future research can accordingly select either a small spacing 4-spot electrode configuration for a single artery sensing or a band electrode configuration for overall pulsatile arteries sensing. A more accurate estimation of blood volume change and pulse wave analysis (PWA) could help to develop cuffless blood pressure measurement (BPM).

The impact of the distribution in space of the inlet velocity in the numerical simulations of the hemodynamics in the thoracic aorta is systematically investigated. A real healthy aorta geometry, for which in-vivo measurements are available, is considered. The distribution is modeled through a truncated cone shape, which is a suitable approximation of the real one downstream of a trileaflet aortic valve during the systolic part of the cardiac cycle. The ratio between the upper and the lower base of the truncated cone and the position of the center of the upper base are selected as uncertain parameters. A stochastic approach is chosen, based on the generalized Polynomial Chaos expansion, to obtain accurate response surfaces of the quantities of interest in the parameter space. The selected parameters influence the velocity distribution in the ascending aorta. Consequently, effects on the wall shear stress are observed, confirming the need to use patient-specific inlet conditions if interested in the hemodynamics of this region. The surface base ratio is globally the most important parameter. Conversely, the impact on the velocity and wall shear stress in the aortic arch and descending aorta is almost negligible.