Software that combines deep learning, 3D reconstruction and CFD to analyze the state of carotid arteries from ultrasound imaging.

BACKGROUND: Ultrasound is one of the non-invasive techniques that are used in clinical diagnostics of carotid artery disease. OBJECTIVE: This paper presents software methodology that can be used in combination with this imaging technique to provide additional information about the state of patient-specific artery. METHODS: Overall three modules are combined within the proposed methodology. A clinical dataset is used within the deep learning module to extract the contours of the carotid artery. This data is then used within the second module to perform the three-dimensional reconstruction of the geometry of the carotid bifurcation and ultimately this geometry is used within the third module, where the hemodynamic analysis is performed. The obtained distributions of hemodynamic quantities enable a more detailed analysis of the blood flow and state of the arterial wall and could be useful to predict further progress of present abnormalities in the carotid bifurcation. RESULTS: The performance of the deep learning module was demonstrated through the high values of relevant common classification metric parameters. Also, the accuracy of the proposed methodology was shown through the validation of results for the reconstructed parameters against the clinically measured values. CONCLUSION: The presented methodology could be used in combination with standard clinical ultrasound examination to quickly provide additional quantitative and qualitative information about the state of the patient's carotid bifurcation and thus ensure a treatment that is more adapted to the specific patient.

Computational Study of Abdominal Aortic Aneurysm Walls Accounting for Patient-Specific Non-Uniform Intraluminal Thrombus Thickness and Distinct Material Models: A Pre- and Post-Rupture Case.

An intraluminal thrombus (ILT) is present in the majority of abdominal aortic aneurysms, playing a crucial role in their growth and rupture. Although most computational studies do not include the ILT, in the present study, this is taken into account, laying out the whole simulation procedure, namely, from computed tomography scans to medical image segmentation, geometry reconstruction, mesh generation, biomaterial modeling, finite element analysis, and post-processing, all carried out in open software. By processing the tomography scans of a patient's aneurysm before and after rupture, digital twins are reconstructed assuming a uniform aortic wall thickness. The ILT and the aortic wall are assigned different biomaterial models; namely, the first is modeled as an isotropic linear elastic material, and the second is modeled as the Mooney-Rivlin hyperelastic material as well as the transversely isotropic hyperelastic Holzapfel-Gasser-Ogden nonlinear material. The implementation of the latter requires the designation of local Cartesian coordinate systems in the aortic wall, suitably oriented in space, for the proper orientation of the collagen fibers. The composite aneurysm geometries (ILT and aortic wall structures) are loaded with normal and hypertensive static intraluminal pressure. Based on the calculated stress and strain distributions, ILT seems to be protecting the aneurysm from a structural point of view, as the highest stresses appear in the thrombus-free areas of the aneurysmal wall.

Effect of Ascending Aortic Curvature on Flow in the Sinus and Neo-sinus Following TAVR: A Patient-Specific Study.

Patient-specific aortic geometry and its influence on the flow in the vicinity of Transcatheter Aortic Valve (TAV) has been highlighted in numerous studies using both in silico and in vitro experiments. However, there has not yet been a detailed Particle Image Velocimetry (PIV) experiment conducted to quantify the relationship between the geometry, flow downstream of TAV, and the flow in the sinus and the neo-sinus. We tested six different patient-specific aorta models with a 26-mm SAPIEN 3 valve (Edwards Lifesciences, Irvine, CA, USA) in a left heart simulator with coronary flow. Velocities in all three cusps and circulation downstream of TAV were computed to evaluate the influence of the ascending aorta curvature on the flow field. The in vitro analysis showed that the patient-specific aortic curvature had positive correlation to the circulation in the ascending aorta (p = 0.036) and circulation had negative correlation to the particle washout time in the cusps (p = 0.011). These results showed that distinct vortical flow patterns in the ascending aorta as the main jet impinges on the aortic wall causes a recirculation region that facilitates the flow back into the sinus and the neo-sinus, thus reducing the risk of flow stagnation and washout time.

The effects of plaque morphological characteristics on the post-stenotic flow in left main coronary artery bifurcation.

Local post-stenotic hemodynamics has critical influence in the atherosclerotic plaque progression occurring in susceptible arterial sites, in particular the left main coronary artery (LMCA) bifurcation. Understanding the effects of plaque morphological characteristics: stenosis severity (SS), eccentricity index (EI) and lesion length (LL) on the post-stenotic flow behavior can significantly improve treatment planning. In order to investigate these effects, we have employed computational fluid dynamics (CFD) simulations in twenty computer-generated and five patient-specific LMCA models and the hemodynamic parameters: velocity, pressure (P), wall pressure gradient (WPG), wall shear stress (WSS), time averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT) and helicity intensity (h2) were analyzed. Our results revealed that the effect of stenosis eccentricity varied significantly for different values of stenosis severity and lesion length. Regions with low WSS, low TAWSS and high RRT were more prominent in models having higher stenosis severity. For smaller lesion length, at low and moderate stenosis severity, surface area with low TAWSS and high RRT decreased with increasing eccentricity index, whereas for high stenosis severity models, low TAWSS region and average RRT values increased with eccentricity. However, for models with longer lesion length, regions with high OSI and RRT overall increased gradually with eccentricity. The helicity intensity (h2) of all models remained very low except at the most eccentric model with longer lesion length. The presence of very high helical flow in this model suggests the possibility of atheroprotective flow. It can be concluded that all plaque morphological characteristics covered under this investigation play an important role in plaque progression.

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

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

Towards optimal flow diverter porosity for the treatment of intracranial aneurysm.

PURPOSE: Low-porosity endovascular stents, known as flow diverters (FDs), have been proposed as an effective and minimally invasive treatment for sidewall intracranial aneurysms (IAs). Although it has been reported that the efficacy of a FD is substantially influenced by its porosity, clinical doctors would clearly prefer to do their interventions optimally based on refined quantitative data. This study focuses on the association between the porosity configurations and the FD efficacy, in order to provide practical data to help the clinical doctors optimize the interventions. METHOD: Numerical simulations in fluid dynamics were performed using four patient-specific IA geometries, pulsatile velocity profiles and braided fully resolved FDs. The variation of velocity and wall shear stress within the IAs, were investigated in this study. Lattice Boltzmann method (LBM) was used to solve the main challenge centered on the diversity of spatial scales since the typical diameter of struts of FDs is only 25µm while the artery normally can be larger by a hundred times. RESULTS: Numerical simulations revealed that the blood flow within IA sac was substantially reduced when the porosity is less than 86%. In particular, the flow condition within each IA sac is favorite to initialize thrombus formation when porosity is less than 70%. CONCLUSION: Our study suggests the existence of a porosity threshold below which the efficacy of a FD will be sufficient for the patients to initialize the thrombus formation. Therefore, by estimating the porosity of FD on patient-specific information, it may be potentially to predict whether or the blood flow condition will successfully become prothrombotic after the FD intervention.

Geometry reconstruction method for patient-specific finite element models for the assessment of tibia fracture risk in osteogenesis imperfecta.

Lower limb deformation in children with osteogenesis imperfecta (OI) impairs ambulation and may lead to fracture. Corrective surgery is based on empirical assessment criteria. The objective was to develop a reconstruction method of the tibia for OI patients that could be used as input of a comprehensive finite element model to assess fracture risks. Data were obtained from three children with OI and tibia deformities. Four pQCT scans were registered to biplanar radiographs, and a template mesh was deformed to fit the bone outline. Cortical bone thickness was computed. Sensitivity of the model to missing slices of pQCT was assessed by calculating maximal von Mises stress for a vertical hopping load case. Sensitivity of the model to ±5 % of cortical thickness measurements was assessed by calculating loads at fracture. Difference between the mesh contour and bone outline on the radiographs was below 1 mm. Removal of one pQCT slice increased maximal von Mises stress by up to 10 %. Simulated ±5 % variation of cortical bone thickness leads to variations of up to 4.1 % on predicted fracture loads. Using clinically available tibia imaging from children with OI, the developed reconstruction method allowed the building of patient-specific finite element models.

Functional assessment of the stenotic carotid artery by CFD-based pressure gradient evaluation.

The functional assessment of a hemodynamic significant stenosis base on blood pressure variation has been applied for evaluation of the myocardial ischemic event. This functional assessment shows great potential for improving the accuracy of the classification of the severity of carotid stenosis. To explore the value of grading the stenosis using a pressure gradient (PG)-we had reconstructed patient-specific carotid geometries based on MRI images-computational fluid dynamics were performed to analyze the PG in their stenotic arteries. Doppler ultrasound image data and the corresponding MRI image data of 19 patients with carotid stenosis were collected. Based on these, 31 stenotic carotid arterial geometries were reconstructed. A combinatorial boundary condition method was implemented for steady-state computer fluid dynamics simulations. Anatomic parameters, including tortuosity (T), the angle of bifurcation, and the cross-sectional area of the remaining lumen, were collected to investigate the effect on the pressure distribution. The PG is highly correlated with the severe stenosis (r = 0.902), whereas generally, the T and the angle of the bifurcation negatively correlate to the pressure drop of the internal carotid artery stenosis. The calculation required <10 min/case, which made it prepared for the fast diagnosis of the severe stenosis. According to the results, we had proposed a potential threshold value for distinguishing severe stenosis from mild-moderate stenosis (PG = 0.88). In conclusion, the PG could serve as the additional factor for improving the accuracy of grading the severity of the stenosis.