CFD investigations of a shape-memory polymer foam-based endovascular embolization device for the treatment of intracranial aneurysms

CFD investigations of a shape-memory polymer foam-based endovascular embolization device for the treatment of intracranial aneurysms | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article CFD investigations of a shape-memory polymer foam-based endovascular embolization device for the treatment of intracranial aneurysms Tanner Cabaniss, Ryan Bodlak, Yingtao Liu, Geoffrey Colby, Hyowon Lee, and 4 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-5014601/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 25 Nov, 2024 Read the published version in Biomechanics and Modeling in Mechanobiology → Version 1 posted 9 You are reading this latest preprint version Abstract The hemodynamic and convective heat transfer effects of a patient-specific endovascular therapeutic agent based on shape memory polymer foam (SMPf) are evaluated using computational fluid dynamics studies for six patient-specific aneurysm geometries. The SMPf device is modeled as a continuous porous medium with full expansion for the flow studies and with various degrees of expansion for the heat transfer studies. The flow simulation parameters were qualitatively validated based on the existing literature. Further, a mesh independence study was conducted to verify an optimal cell size and reduce the computational costs. For convective heat transfer, a worst-case scenario is evaluated where the minimum volumetric flow rate is applied alongside the zero-flux boundary conditions. In the flow simulations, we found a reduction of the average intra-aneurysmal flow of > 85% and a reduction of the maximum intra-aneurysmal flow of > 45% for all presented geometries. These findings were compared with the literature on numerical simulations of hemodynamic and heat transfer of SMPf devices. The results obtained from this study can serve as a guide for optimizing the design and development of patient-specific SMPf devices aimed at personalized endovascular embolization of intracranial aneurysms. Computational fluid dynamics intracranial aneurysms cerebral aneurysms shape memory polymers patient-specific therapeutics Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Figure 9 1. Introduction Intracranial aneurysms (ICAs), which affect approximately 4% of the American population (Keedy et al. 2006), are saccular protrusions that arise form focal dilation of the arterial wall (Toth et al. 2018). It has been found that 85% of ICAs develop within the circle of Willis , the complex central cranial vasculature (Keedy et al. 2006). Although there are different classifications of ICAs, namely saccular, fusiform, dissecting, and mycotic, this study focuses exclusively on saccular aneurysms as they account for approximately 90% of ICAs (Keedy 2006 ). The primary risk of ICAs is the possibility of rupture, which can further lead to subarachnoid hemorrhage: a bleeding in the intracranial cavity that leads to a variety of life-threatening conditions such as severe brain damage, long-term morbidity and/or mortality (Schievink et al. 1997). While the causative physiological mechanisms for ICA rupture are not yet fully understood, current literature presents correlations with multiple factors, including age, hypertension, smoking history, alcohol abuse, hypercholesterolemia, and carotid stenosis (Wagner et al. 2005). The success rate of ICA treatment is mainly measured using the Raymond-Roy Occlusion Classification scale. This classification evaluates the degree of occlusion of the aneurysm sac into three distinct classes: Class I – complete obliteration; Class II – residual neck; and Class III – residual aneurysm (Roy et al. 2001 ). Complete obliteration, also called complete occlusion, is defined as the cessation of intra-aneurysmal blood circulation, with the majority of blood flowing through the parent artery as if the aneurysm were not present. Guglielmi detachable coils (GDCs), developed in the 1990s (Guglielmi et al. 1992 ), have long been considered the gold standard for endovascular treatment of ICAs. According to a retrospective study by Gonda et al. ( 2014 ) over 2,509 unruptured aneurysms, while GDCs have a 30-day postoperative mortality rate lower than that of highly invasive microsurgical clipping (1.1% vs. 2.3%), in the long term they have unsatisfactorily lower rates of complete occlusion (i.e., Raymond-Roy Occlusion Classification Class I). Another retrospective study by Murayama et al. ( 2003 ) reported that traditional coiling had an astonishing complete closure rate of 55%, compared to complete occlusion rates of > 90% with microsurgical clipping (Natarajan et al. 2008 ; Schwartz et al. 2018 ; Choi et al. 2020 ). To address this critical shortcoming in aneurysm occlusion, new endovascular devices have been developed over the past 2–3 decades. These devices consist of direct improvements to GDCs through the introduction of hydrogel and shape memory polymer (SMP)-based coatings as well as novel endovascular devices such as stents and flow diverters (Alderazi et al. 2014 ; Taschner et al. 2018 ; Herting et al. 2019 ). While stent-based devices have demonstrated improved complete occlusion rates of up to 83.3% (Hanel et al. 2023 ), flow diversion via stents is only useful in the treatment of sidewall aneurysms, so bifurcation aneurysms continue to be treated primarily via coil-based techniques. Although these results provide significant room for improvement, new endovascular technology may be the answer to achieving comparable rates of complete occlusion to microsurgical clipping procedures without the associated procedural risks of a craniotomy. This emerging technology is based on the use of shape memory polymer foams (SMPf) – a subgroup of SMP materials that possess a high degree of porosity while maintaining the shape recovery properties of typical shape memory polymers (Pineda-Castillo et al. 2022 ). These properties allow the SMPf to be compressed into a microcatheter, implanted into an aneurysm cavity, and subsequently expand and recover to its original programmed geometry via shape recovery activation, thereby occluding the aneurysm sac. The SMPf have been fabricated with both closed and open cell architectures, the latter enabling higher compression ratios required for endovascular techniques (Wang et al. 2019 ). Moreover, recent work has used discrete lattice architectures to enable further degrees of structural optimization (Pineda-Castillo et al. 2023 ). To trigger the shape recovery of SMPf through Joule heating, conductive polymer coatings and nanomaterials such as carbon nanotubes were introduced to increase the electrical and thermal conductivity of the fabricated SMPf (Pineda-Castillo et al. 2020 ; Pineda-Castillo et al. 2021 ). Horn et al. ( 2017 ) conducted in vivo studies on porcine side-wall aneurysms that demonstrated the promise and efficacy of this therapeutic technique. However, several inquiries remain regarding the optimal geometric considerations of these foams. These salient challenges revolve around crucial aspects of foam design such as mechanical stability, optimal configurations to reduce intra-aneurysmal residual circulation, and promotion of platelet aggregation to improve coagulation processes. One of the earliest applications of computational fluid dynamics (CFD) in aneurysm research came from Gonzalez et al. ( 1992 ), who used pulsatile flow and a non-Newtonian fluid model to simulate blood flow through an aneurysm in the middle cerebral artery. By 2004, only eleven additional articles had been published on the topic, but research efforts have increased rapidly since then, as a study by Saqr et al. ( 2020 ) found 795 publications on ICA hemodynamics using CFD as the primary in silico modeling tool. While there is a wealth of information about hemodynamics in ICAs and the effect of other devices (e.g., stents, coils, etc.), only limited efforts have been made to investigate the impact of SMPf-based devices for personalized treatment of ICAs. The aim of this study is therefore to use CFD to examine the hemodynamic influence of SMPf implantation in “patient-specific” aneurysm geometries. Specifically, our aim is to quantify the reduction in intra-aneurysmal blood flow caused by the SMPf device and to evaluate the convective heat transfer to the surrounding arterial tissues during initiation of shape recovery as well as successive device expansion. The results of this work could also lead us to a deeper understanding of the relevant intricacies of SMPf devices for patient-specific endovascular ICA therapeutics. The remainder of this paper is organized as follows: in Section 2 , we describe the development of a CFD-based approach to analyze the hemodynamics of untreated ICAs and aneurysms implanted with an SMPf-based device, as well as the heat transfer mechanisms associated with activation of shape recovery. In Section 3 , we present the results of the models and validate their efficacy using an example from the established literature. In Section 4 , we will further assess the advantages and limitations of the patient-specific SMPf device when compared to other established endovascular techniques, followed by some concluding remarks. 2. Methods 2.1 Geometry Preparation Eleven anonymized, patient-specific 3D rotational angiography datasets were acquired from the University of Oklahoma Health Sciences Center in 2018 and 2019 (IRB #7932). Of these datasets, nine had sufficient resolution and were selected for image segmentation. These were successively narrowed down to five unique geometries from the collected data; three bifurcation and two sidewall aneurysms as well as another geometry from the literature (Valen-Sendstad et al. 2018 ) for our numerical modeling/simulation. The four segmented geometries excluded from this study were structurally similar to the geometries presented. Image segmentation was performed in Amira (Thermo Fisher Scientific, Waltham, MA) with \(\:0.108\times\:0.108\times\:0.5\) mm rectangular prism voxels. The segmented patient-specific ICA lumen surface was then imported into Meshmixer (Autodesk, San Rafael, CA) as an STL file and slightly smoothed to remove localized peaks and kinks due to noise in the image data. The vessels were then cut normal to their cross sections to form inlets and outlets, while the outlets were maintained at approximately twice the vessel diameter once their respective directions were established, and the inlets were cut at least seven vessel diameters from the aneurysm sac, as suggested by the literature convention (Moyle et al. 2006 ; Hua et al. 2015 ). These cuts were made to reduce the computational costs while maintaining validity of the flow characteristics throughout the studied system (Hodis et al. 2015 ). The full extent of the modeled vasculature for each geometry and the respective more detailed views of the ICAs are shown in Fig. 1 , while the main geometric characteristics of the modeled aneurysms are summarized in Table 1 . Here, the bifurcation aneurysm B1 was taken from the 2015 International Aneurysm CFD Challenge (Valen-Sendstad et al. 2018 ). Table 1 Characteristics of the six patient-specific aneurysms studied (B1-B4 and S1, S2). Location l max (mm) AR SR V (mm 3 ) A ostium (mm 2 ) B1 MCA bifurc. 7.7 1.0 1.8 76 28.2 B2 MCA bifurc. 7.8 1.5 2.1 50 12.4 B3 MCA bifurc. 5.6 1.6 1.9 41 21.0 B4 Basilar apex 11.1 1.3 2.5 313 36.1 S1 Sup. hypophyseal 9.0 2.3 2.0 167 24.4 S2 Ant. choroidal 2.9 0.8 0.4 10 6.6 MCA birfurc. = middle cerebral artery bifurcation; Sup. hypophyseal = superior hypophyseal artery; Ant. choroidal = anterior choroidal artery; l max = facet diameter; AR = aspect ratio (sac depth: neck width); SR = size ratio (sac depth: parent artery diameter); V = volume of aneurysm sac; A ostium = interface area between aneurysm sac and main vessel. 2.2 Simulation Parameters Since one of the goals of this study is to support the development of individually-optimized SMPf devices (Kunkel et al. 2018 ;Wang et al. 2019 ), steady-state simulations were performed in OpenFOAM 8 (Weller et al. 1998 ), assuming that the greatest degree of blood flow into the SMPf device occurs during systole. Moreover, because steady-state simulations with systolic flow tend to overestimate the velocity field (Geers et al. 2014 ), this consideration provided an upper bound on the worst-case flow that the SMPf device and the treated aneurysm would experience. Accordingly, steady-state solutions were achieved using the SIMPLEC (Semi-Implicit Method for Pressure Linked Equations-Consistent) algorithm and intra-aneurysmal blood flow residuals were monitored to ensure convergence. The meshes were constructed predominantly from hexahedron cells, while prismatic and tetrahedral cells were incorporated along the boundaries with a growth factor of 1.2 in the 12–20 inflation layers, resulting in mesh sizes from \(\:4.4\times\:{10}^{6}\) to \(\:7.9\times\:{10}^{6}\) cells. The results of the mesh selection process and convergence study are presented in Sections 2.6 and 3.1 . In these simulations, the blood flow was modeled to be an incompressible Newtonian fluid with a dynamic viscosity \(\:\mu\:=3.65\times\:{10}^{-3}\) Pa·s, a density \(\:\rho\:=1050\) kg/m 3 and a kinematic viscosity \(\:\upsilon\:=3.5\times\:{10}^{-6}\) m 2 /s (Hasgall et al. 2018 ). No material was assigned to the vessel walls, but they were treated as rigid according to convention, with a no-slip boundary condition applied along the walls (Berg et al. 2019 ). Next, the inlet boundary conditions were defined using the peak systolic volumetric flow rate (VFR) values estimated from population-averaged flow waveforms for the ICA and basilar artery reported in the literature for a normal, young cohort: 9.26 mL/s and 5.63 mL/s, respectively (Horn et al. 2017 ). Generalized waveforms for older adults were avoided due to the age range of the present cohort was unknown, although both datasets would likely have produced similar flow patterns and wall shear stresses (WSS) (Xu et al. 2018 ). According to Chnafa et al. ( 2018 ), using the splitting method outlined above, VFRs were prescribed on all but one outlet, typically the smallest. The unassigned outlet was left in the zero pressure condition while a zero-gradient pressure boundary was applied to all other outlets. Finally, there is increasing evidence that transitional or turbulent flow occurs in physiological blood flow (Saqr et al. 2020 ) and that turbulent flow has been observed in the aneurysm space after systole (Valen-Sendstad et al. 2011 ). Accordingly, the \(\:k-{\omega\:}\) shear stress transport turbulence model (Menter et al. 1994) was adopted to realize the turbulence characteristics in the simulated intracranial circulation environment. 2.3 Simulations of an SMPf-Based Device Implementation Treatment efficacy was evaluated by virtually implanting SMPf devices into each patient-specific aneurysm geometry. First, the geometry of each SMPf device was generated by manually excising the aneurysm sac from the surrounding vasculature (Fig. 2 ). The devices were then idealized and represented the maximum possible foam volume for the specified geometries. Manufacturing considerations or defects/discrepancies were not taken into account. Because in this study we focused on providing general design suggestions rather than seeking to evaluate a patient-specific case in a clinical setting, we chose a continuum porous media approach to model the SMP foam. The porous medium specifications considered in this study followed the previous work by Muschenborn et al. ( 2013 ) with a porosity \(\:\varPhi\:=98\%\) , permeability \(\:k=3.25\times\:{10}^{-9}\sim1.24\times\:{10}^{-8}\) m 2 , and form factor \(\:C=2720-7260\) m −1 . These parameter ranges were verified by evaluating the pressure gradient of the flow through the porous media model at different flow velocities and comparing them with the original study (Muschenborn et al. 2013 ). For the results presented in this work, the lowest permeability model was selected because it demonstrated the greatest performance in reducing intra-aneurysmal flow. 2.4 Simulations of Heat Transfer Process During Initial Shape Recovery that Triggers Heating To quantify the thermal effects of device deployment and the heating associated with the initiation of shape recovery, an idealized, rigid-walled solid cylinder with a diameter of 1 mm was used. The use of a compressed solid model was justified because the pores of the closed cell SMPf were almost completely closed upon compression (Fig. 3 ). The length of the cylinder models varied between ICA cases, but in all cases ranged from the aneurysm neck to a 1 mm offset from the vessel walls perpendicular to the neck, centered in the intra-aneurysmal space. The potential thermally induced tissue damage can be described by the cumulative equivalent minutes at 43°C (CEM 43 ), defined by $$\:CE{M}_{43}\left(t\right)={\int\:}_{{t}_{0}}^{{t}_{f}}{R}^{\left(43-T\left(t\right)\right)}\text{d}t$$ 1 , where \(\:{t}_{0}\) is the initial time, \(\:{t}_{f}\) is the end time of exposure in minutes, \(\:R\) is a constant of 0.5 for \(\:T\left(t\right)\ge\:43\:℃\) , 0.25 for \(\:39\:℃\le\:T\left(t\right)<43\:℃\) , and 0.0 for \(\:T\left(t\right)<39\:℃\) (Sapareto et al. 1984). According to van Rhoon et al. ( 2013 ), irreversible damage to brain tissues occurs with as little as CEM 43 =10. Assuming a constant exposure temperature, \(\:40\:℃\) and \(\:45\:℃\) stimuli would result in tissue damage after > 600 min and 2.5 min, respectively. Therefore, the SMPf device in this study was considered to have a glass transition temperature T g of \(\:39\:℃\) (Pineda-Castillo et al. 2021 ). Pineda-Castillo et al. ( 2021 ) have shown that the SMPf device can heat from well below body temperature to its transition temperature in less than a minute while completely occluding an in vitro aneurysm sac in the same amount of time. Thus, a fixed temperature of \(\:40\:℃\) was used in the heating portion of this numerical study, just above the T g of the SMPf material. A zero-heat flux (perfectly insulated) boundary condition was imposed on the vessel to model the worst-case temperature rise (Ortega et al. 2007 ). The flow through the inlet was set to body temperature ( \(\:37\:℃\) ), with the blood being assigned a specific heat capacity of \(\:{c}_{p}=3600\) J/kg· \(\:℃\) and a thermal conductivity of \(\:{k}_{t}=0.52\) W/m· \(\:℃\) (Hasgall et al. 2018 ). The inlet flow for this study was modeled at the diastolic VFR of 2.98 mL/s and 1.77 mL/s for the internal carotid artery and basilar artery, respectively (Gwilliam et al. 2009 ) to minimize convection from the region and to simulate the highest possible temperature rise on the walls. The simulations were considered steady-state and included the Reynolds-averaged Navier-Stokes equations, although the slow flow rates significantly reduced the instabilities observed in the previous systolic flow simulations. 2.5 Simulations of Heat Transfer During Device Expansion Next, transient heat transfer simulations were performed on a representative patient geometry, Case B2, which demonstrated the highest level of heat transfer during the initial heating stage of the deployment. To limit the computational expense of the additional solution steps, the domain upstream of the ICA was shortened to approximately three diameters of the parent artery (Fig. 4 ) and remeshed with a cell size of 0.20 mm to generate approximately 41,000 cells. To provide a more realistic flow condition than the previous diastolic study, outflow VFRs were calculated from the mean ICA velocity. All other boundary conditions were retained. Although the simulation was transient, the flow was still modeled as steady-state rather than pulsatile, using a time step of 1 ms and a PIMPLE algorithm, an iterative solver which is a combination of PISO (Pressure Implicit with Splitting of Operator) and SIMPLE (Semi-Implicit Method for Pressure-Linked Equations). For computational efficiency reasons, the total time steps were limited to 1,000. The geometry was constructed as a partially expanded foam with an offset of 0.5 mm from the vessel wall. As in the flow simulations, the model was defined as a porous medium with a permeability of \(\:5.7\times\:{10}^{-9}\) m 2 and form factor of 7260 m − 1 . The entire zone of the porous medium was assigned a temperature of \(\:40\:℃\) and \(\:45\:℃\) for the two respective shape recovery-induced heating processes to provide an optimal and worst-case heating scenario. The surrounding fluid volume was modeled similarly to that described in Section 2.4 . 2.6 Simulation Verification A mesh independence study was performed on the aneurysm B1 geometry (Fig. 1 ), with simulations covering a range of 0.05–0.16 mm cell sizes ( \(\:6.87\times\:{10}^{5}-1.25\times\:{10}^{7}\) cells). The number of cells was approximately doubled for each successive cell size, without the finest mesh containing 1.5 times as many cells as the second finest mesh. The maximum and average flow velocities for 1,000 points evenly distributed throughout the aneurysm domain were calculated in ParaView (Ahrens et al. 2005 ). To reduce the influence of insignificant errors/discrepancies in the simulation, the slowest 1% of the sampled points, corresponding to a velocity of 0.03 m/s, were excluded from the analyses. To further verify the simulation parameters and results, the mean internal carotid artery flow rate of 5.11 mL/s (Gwilliam et al. 2009 ) was scaled and applied to a 0.08 mm mesh. This enabled a direct comparison of WSS values between this study and the time-averaged WSS results of the six most experienced teams at the 2015 International Aneurysm CFD Challenge (Valen-Sendstad et al. 2018 ). However, despite their experience, there was no clear consensus regarding the WSS magnitude numerically determined between the six teams: i.e., four teams achieved comparable values and the other two teams reported drastically different values. While a direct quantitative and numerical comparison may not have a high degree of validity in this case, qualitative contrasts can be drawn based on the WSS distributions presented by the teams as they are holistically more similar. 2.7 Simulation Metrics and Data Analysis For the flow velocity studies, the WSS distribution was beneficial to qualitatively verify the validity of the flow within the model. Intra-aneurysmal velocity is also an important metric to evaluate since reducing flow is a necessity for rebuilding and recovering arterial tissue. Intra-aneurysmal velocity was evaluated for maximum and mean quantities for the interior of the aneurysm space. For the heat transfer studies, the temperature for the intra-aneurysmal space and surrounding vessels was assessed, and visualized using isosurfaces and contours. For all studies presented, ParaView (Ahrens et al. 2005 ) was utilized to extract the above metrics from the simulated data and generate visualizations. The Python library Matplotlib (Hunter 2007 ) was used to generate the CEM 43 plots. 3. Results 3.1 Study on Mesh Size The mesh size was determined to be \(\:s=0.08\) mm as this was the maximum cell size with less than 5% mean absolute error (MAE) compared to the control ( \(\:s=0.05\) mm), resulting in a cell density of approximately 1,700-1,800 cells/mm 3 , which is consistent with the cell densities presented in the literature (Marzo et al. 2009 ; Hodis et al. 2012 ; Aranda et al. 2018). The mesh convergence study results are summarized in Table 2 . Further qualitative comparisons of the WSS contours with those of the 2015 International Aneurysm CFD Challenge (Valen-Sendstad et al. 2018 ) showed good agreement in the WSS distributions along the aneurysm surface (Fig. 5 ). Table 2 Detailed results of the mesh independence study. (Velocity values are calculated over the entire aneurysm domain.) s (mm) Cells ( \(\:\times\:{10}^{3}\) ) IA cells ( \(\:\times\:{10}^{3}\) ) u max (m/s) u avg (m/s) MAE (%) 0.05 12450 938 1.106 0.494 – 0.06 7780 612 1.103 0.494 1.5 0.08 4160 329 1.097 0.494 4.4 0.10 2230 194 1.086 0.494 9.0 0.13 1150 106 1.061 0.490 15.9 0.16 687 67 1.054 0.490 18.3 s = mesh size; IA = intra-aneurysmal; u max = maximum velocity; u avg = average velocity; MAE = mean absolute error. 3.2 Simulations of Untreated and Treated ICAs As discussed in Section 2.3 , porous media models were used to define the permeability of the SMPf device. For the results presented here, the lowest permeability of the foam was adopted, defined by a permeability of \(\:5.7\times\:{10}^{-9}\) m 2 and a form factor of 7260 m − 1 . In each case, treatment of the aneurysm showed significant decrease in blood flow into the aneurysm cavity. Overall, in each case, the maximum flow velocity was reduced by > 40% and the average flow velocity was reduced by up to 95%, as shown in Table 3 . The metric used to quantify this reduction was the velocity reduction ratio VRR, i.e. $$\:VRR=\frac{{u}_{\text{p}\text{r}\text{e}}-{u}_{\text{p}\text{o}\text{s}\text{t}}}{{u}_{\text{p}\text{r}\text{e}}}\:,$$ 2 where u is the predicted fluid velocity and the subscripts, pre and post, denote the respective velocities before and after device implantation, respectively. Table 3 Simulation results of intracranial flow velocity before and after SMPf implantation. All velocity values are given in m/s. u max,pre u max,post VRR max u avg,pre u avg,post VRR avg B1 1.107 0.456 58.8% 0.494 0.054 89.0% B2 0.746 0.387 48.2% 0.172 0.017 90.0% B3 0.775 0.333 57.1% 0.328 0.033 90.1% B4 1.180 0.488 58.7% 0.333 0.022 93.3% S1 0.997 0.428 57.1% 0.219 0.011 95.0% S2 0.825 0.329 60.1% 0.403 0.048 88.0% u = predicted fluid velocity; pre = before device implantation; post = after device implantation; VRR = velocity reduction ratio; Furthermore, this residual flow reduction can be visualized using the velocity contours of the central aneurysmal plane and inlet vessel for both the untreated and treated cases (Fig. 6 ). More profound differences were found in the bifurcation cases (B1, B2, B3, and B4), where high-velocity inflow jets were observed prior to treatment. In addition to reducing flow within the intra-aneurysmal space, device implantation appeared to have stabilized the downstream flow, as indicated by the smoother velocity contours in these areas. While the flow reduction results are promising for high efficacy of a patient-specific SMPf-based device, the deployment of these devices introduces several biocompatibility concerns, with thermal damage to surrounding arterial tissue being the greatest concern and requiring further investigation. 3.3 Heat Transfer Simulations for Initial Heating of the SMPf Device – Triggering Shape Recovery As described previously, arterial tissue damage occurs near CEM 43 of 10 minutes (van Rhoon et al. 2013 ). Therefore, the thermal effects of the device deployment process on the aneurysm and artery wall tissue must be carefully investigated to identify possible thermal damage with respect to the appropriate activation temperature for shape recovery. The first case considered was the initial thermal activation heat transfer, as at this point the convection patterns from the SMPf device to the surrounding fluid became visible. Temperature isosurfaces were plotted to observe heat convection and the distribution of the initial activation heat to the surrounding fluid (Fig. 7 ). The selected isosurface was chosen at \(\:37.5\:℃\) just \(\:0.5\:℃\) above body temperature, to maximize interpretable data of the underlying heat transfer process and impact. In general, the \(\:37.5\:\text{℃}\) isosurfaces (Fig. 7 ) did not contact any walls within the aneurysm space or downstream vasculature, which was expected given the minimal temperature gradient generated during device activation at \(\:40\:℃\) . However, there are two exceptions. In case B1, a thin stream of heated fluid impinges on the wall at a downstream bifurcation, although heat transfer thereby is negligible. In addition, given the low fluid velocities compared to the other cases (Fig. 6 ), it was not surprising that Case B2 was somewhat more concerning as it was found that convection within the dome was significantly hindered. This allows fluid at a higher temperature level to come into contact with the vessel wall at the apex of the treated aneurysm. This resulted in an open isosurface at \(\:37.5\:℃\) , indicating local wall temperatures above body temperature. Nevertheless, with the \(\:39\:℃\) isosurface superimposed, the temperatures on the vessel wall were still lower than those of the device ( \(\:40\:℃\) ). Because Case B2 showed the greatest concern for exposing vessel walls to potentially degenerative temperatures, this case was utilized to further explore the impact of heat transfer during the subsequent expansion phase of the device. 3.4 Transient Heat Transfer During SMPf Device Expansion and Assessment of Tissue Damage The maximum wall temperature was observed to increase to \(\:40\:℃\) within 1 s at both stimulus temperatures, which is due to the relatively slow flow within the aneurysm dome and the small distance between the SMP foam and the tissue wall (Fig. 8 ). It was found that the blood flow pattern is such that only certain regions of the aneurysm sac reach elevated temperatures, while almost all of the heated blood from those areas drains into only one of the bifurcation vessels, significantly increasing its wall temperature. The time histories for the heating and cooling process of the SMPf device shape recovery can be combined to calculate the CEM 43 dosage for Case B2. Assuming a total heating time of 1 minute (Pineda-Castillo et al. 2021 ), the maximum vessel wall temperature curves were generated throughout the entire shape recovery process (Fig. 9 (a)). Since the blood flow temperature for the \(\:40\:℃\) stimulus was found to never exceed \(\:40\:℃\) , the CEM 43 measurement is almost negligible for the entire heating process (i.e., 0.02 min). In contrast, the \(\:45\:℃\) stimulus resulted in a worst-case CEM 43 of 4.9 min. Even in this (over-heating) worst-case scenario, there would be no thermal damage to the surrounding arteries/aneurysm wall (van Rhoon et al. 2013 ). The general relationship between time and temperature for CEM 43 is shown in Fig. 9 (b). 4. Discussion In this study, we examined the viability of a patient-specific SMPf device for the treatment of ICAs. We explored the hemodynamical effects of device implantation on intra-aneurysmal flow as well as the implications of a thermal heating mechanism within the aneurysm cavity, taking into account both the initial heating stages and transient heating during device expansion. Below we will discuss the potential advantages and disadvantages regarding the efficacy of the proposed SMPf patient-specific embolic device and the validity of the numerical studies presented here. 4.1 Flow Simulations In this work, we presented baseline and treatment flow simulations for six different patient-specific geometries, with the SMPf device modeled as a porous medium. In the untreated cases, we found an average maximum velocity in the aneurysm space of \(\:0.94\pm\:0.18\) m/s and an overall average velocity of \(\:0.32\pm\:0.12\) m/s across all patient geometries. After treatment, we found an average maximum velocity of \(\:0.40\pm\:0.07\) m/s and an overall average velocity of \(\:0.031\pm\:0.017\) m/s resulting in an approximate reduction of 60% and 90% in the maximum and average flow velocities, respectively. While the low permeability foam used in our numerical simulations showed promising results, in real practical applications higher permeability may be required to meet the required compression constraints as the device compresses radially and is inserted into an endovascular micro-catheter with an inner diameter of < 1 mm for use in clinical setting. In our internal pilot studies, not presented here, we found that higher permeability foams do not perform as well as the idealized low permeability foam presented in this manuscript. A study by Ouared et al. ( 2016 ) found that with vascular stents, even a 30% reduction in intra-aneurysmal flow can promote thrombosis in the aneurysm sac, which can ultimately lead to effective aneurysm occlusion. Therefore, even with higher permeability foam, there is still evidence for the viability of the SMPf device presented. In addition, Ortega et al. ( 2013 ) evaluated an SMPf device with higher permeability ( \(\:k\approx\:2\times\:{10}^{-8}\) m 2 ) and found a flow reduction in the intra-aneurysmal region furthest from the aneurysm neck up to 95%. Their results are consistent with our results presented in this work. While endovascular coil embolization remains one of the most prominent methods for the endovascular (minimally invasive) treatment of ICAs, the SMPf device offers potential solutions to some of the more common limitations of GDC coiling therapeutic. Most importantly, the patient-specific SMPf aims to maximally reduce the residual volume of the aneurysm, which is a cause of recurrence in patients treated with coils and requires repeated surgery (Griessenauer et al. 2016 ). Another shortcoming of coil embolization that the proposed SMPf device aims to improve is post-treatment compaction. Because the soft platinum wires used in coiling are delicate and flexibly deformable, clinical evidence has shown that blood inflow into the aneurysm sac can gradually compact the packed coils, further leading to residual volume near the aneurysm neck and causing recurrence (Yu et al. 2019 ). On the other hand, the SMPf device would have a higher degree of structural integrity due to its highly interconnected cellular structure. 4.2 Device Heating and Deployment Heat transfer from SMPf device deployment and initiation of shape recovery is highly dependent on the specific aneurysm geometry and flow speed (Fig. 7 ). For example, transient analysis of Case B2, which showed the highest level of heat transfer from the device to the aneurysmal and arterial walls, indicated that the wall temperature would approach that of the expanding SMPf within seconds. While this was most clearly observed in Case B2, it is likely that this condition applies to all study cases during further stages of device expansion as fluid is forced from the intra-aneurysmal space and convection capacity is reduced. However, the potential thermal damage to tissue from this heating process depends heavily on the temperature required for the device to expand (e.g., a thermal activation threshold very close to body temperature would still not cause damage to surrounding tissue regardless of the duration of the activation/expansion process). Furthermore, the heterogeneity of the arterial tissue composing the aneurysm sac, both in terms of microstructural composition and tissue thickness (Laurence et al., 2021 ), may result in a non-uniform response to the heating stimulus. This can potentially lead to local thermal damage below the previously mentioned stimulus thresholds. Other possible solutions to the heating problem include reducing the time required for full expansion, which is also related to the material properties of the SMPf, as the discrepancy between the temperature provided to the device and the T g of the device regulates the rate at which expansion occurs (Kunkel et al. 2018 ; Shahi et al. 2022). This positive correlation between shape recovery rate and stimulus temperature is related to the increased kinetic energy within the polymer system at higher temperatures, which allows the polymer chains to move more freely and reestablish their initial geometric state more quickly (Davidson et al. 2015). Furthermore, although the boundary conditions presented in this study provided the stimulus temperature ( \(\:40\:℃\) and \(\:45\:℃\) ) on the exterior of the SMPf model, in practice with a central heating element, the low thermal conductivity of a SMPf (Soloveva et al. 2022 ) would likely have a strong thermal gradient that reduces heat transfer to the blood and surrounding tissues during all stages of deployment. 4.3 Study Limitations and Future Work Although the results presented in this study are promising, there are notable limitations that should also be discussed. First, no experimental validation was performed to verify the accuracy of the presented results. Second, the porous models presented in the study were all isotropic in nature, ignoring the potential effects of anisotropic foam, which are more common in practice due to the foaming process (Ortega et al. 2013 ). Furthermore, the steady-state assumption neglected potential benefits of the foam-based device, e.g., previous research studies have shown that SMPf can filter out high-frequency velocity fluctuations in the aneurysm space (Ortega et al. 2013 ). In short, the SMPf device promises to reduce intra-aneurysmal flow velocities and offers potential solutions to the limitations encountered with traditional coil-based methods. Now we will discuss some possible drawbacks arising from the thermal activation mechanism for SMPf deployment. The most urgent extension of this work is the construction of an in vitro flow loop with an integrated particle image velocimetry system for experimental validation of the studies reported herein. Moreover, extensions of computational simulations with pulsating, time-dependent conditions can provide insights into flow behavior caused by inertial effects not foreseen in this work. An analysis of the deployment of the device would be helpful in determining the viability of the underlying concepts of a SMPf device (e.g., examining the stresses that the expanding foam exerts on the vessel wall in the event of device misalignment, similar to the work of Hwang et al. 2012 ). Further analysis of the thermal properties of the specific SMPf material and inclusion in the aforementioned heat transfer simulations would provide a more realistic description of the plausibly harmful heat generation of the device during expansion in the intra-aneurysmal space. Another possible extension is the inclusion of a fluid-structure interaction (FSI) approach to model blood flow-induced deformations in the aneurysm and foam. With such an approach, the surfaces are not rigid and can respond to fluid motion. Accordingly, FSI modeling is more computationally expensive and requires a higher degree of boundary conditions and model data since both the fluid and solid domains are modeled. If properly performed, FSI modeling would offer a more accurate physiological representation of the blood flow patterns in the ICA and the effects of the proposed SMPf device on it. 4.4 Conclusion In this work, we presented results on the hemodynamic effects and possible heat transfer impairments of a SMPf device for the treatment of ICAs. We have demonstrated a significant flow reduction in the intra-aneurysmal space, effects of initial shape recovery on various aneurysm geometries, and, in the worst case, transient heat transfer to the surrounding tissue. Although these results were promising and prospective solutions have been developed to address concerns about device heat transfer during shape recovery induction, future work is still required to address the full complexity of treating ICAs and develop a highly effective treatment device. However, this work adds to the growing body of research assessing the effectiveness of SMPf devices for the endovascular treatment of ICAs, thereby providing additional evidence for further investigation and the design and development of personalized SMPf-based endovascular devices. Declarations conflict of interest The authors of this paper have no financial or personal relationships with other people or organizations that could inappropriately influence (bias) our work. Use of human subject data Deidentified angiographic imaging data from patients were collected and used in this study. The study protocol for human subjects was reviewed and approved by the Institution Review Board of the University of Oklahoma Health Sciences Center (IRB #7932). Author Contribution T.L.C.: Methodology, Figure Preparation, Manuscript Writing/EditingR.B.: Conceptualization, Methodology, Figure Preparation, Manuscript WritingY.L.: Conceptualization, Resources, Student/Researcher Mentoring, Dissemination, Manuscript Review/EditingG.P.C.: Conceptualization, Manuscript Review/EditingH.L.: Conceptualization, Manuscript Review/EditingB.N.B.: Conceptualization, Methodology, Project Administration, Manuscript Review/EditingR.G.: Conceptualization, Resources, Manuscript Review/EditingG.A.H.: Conceptualization, Methodology, Manuscript Review/EditingC.-H.L.: Conceptualization, Methodology, Project Administration, Resources, Student/Researcher Mentoring, Dissemination, Manuscript Review/Editing Acknowledgements This project was supported by grants from the Oklahoma Center for the Advancement of Science and Technology (OCAST, grant no. HR18-002), and the Oklahoma Shared Translational and Clinical Research (OSCTR, NIGMS U54GM104938). In addition, this publication reports research supported by the National Heart, Lung, And Blood Institute of the National Institutes of Health under Award Number R01HL159475, and by the University of Parma through the action Bando di Ateneo 2021 per la ricerca co-funded by MUR-Italian Ministry of Universities and Research – D.M. 737/2021-PNR-PNRR – NextGenerationEU. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Data Availability Data is provided within the manuscript or supplementary information files. References Ahrens J, Geveci B, Law C, Hansen C, Johnson C (2005) 36-paraview: An end-user tool for large-data visualization. The Visualization Handbook, 717, 50038-50031. https://doi.org/10.1016/B978-012387582-2/50038-1 Alderazi YJ, Shastri D, Kass-Hout T, Prestigiacomo CJ, Gandhi CD (2014) Flow diverters for intracranial aneurysms. Stroke Res. Treat, 2014, 415653. https://doi.org/10.1155/2014/415653 Aranda A, Valencia A (2018) Study on the relationship between wall shear stress and aspect ratio of cerebral aneurysms with different pressure differences using CFD simulations. J Mech Med Biol, 18(05), 1850055. https://doi.org/10.1142/S0219519418500550 Berg P, Saalfeld S, Voß S, Beuing O, Janiga G (2019) A review on the reliability of hemodynamic modeling in intracranial aneurysms: Why computational fluid dynamics alone cannot solve the equation. Neurosurg Focus, 47(1), E15. https://doi.org/10.3171/2019.4.FOCUS19181 Chnafa C, Brina O, Pereira V, Steinman D (2018) Better than nothing: A rational approach for minimizing the impact of outflow strategy on cerebrovascular simulations. AJNR Am J Neuroradiol, 39(2), 337-343. https://doi.org/10.3174/ajnr.A5484 Choi JH, Lee KS, Kim B-S, Shin YS (2020) Treatment outcomes of large and giant intracranial aneurysms according to various treatment modalities. Acta Neurochir, 162, 2745-2752. https://doi.org/10.1007/s00701-020-04540-1 Davidson JD, Goulbourne NC (2015) Microscopic mechanisms of the shape memory effect in crosslinked polymers. Smart Mater Struct, 24(5), 055014. https://doi.org/10.1088/0964-1726/24/5/055014 Geers A, Larrabide I, Morales H, Frangi A (2014) Approximating hemodynamics of cerebral aneurysms with steady flow simulations. J Biomech, 47(1), 178-185. https://doi.org/10.1016/j.jbiomech.2013.09.033 Gonda DD, Khalessi AA, McCutcheon BA, Marcus LP, Noorbakhsh A, Chen CC, Chang DC, Carter BS (2014) Long-term follow-up of unruptured intracranial aneurysms repaired in California. J Neurosurg, 120(6), 1349-1357. https://doi.org/10.3171/2014.3.JNS131159 Gonzalez CF, Cho YI, Ortega HV, Moret J (1992) Intracranial aneurysms: Flow analysis of their origin and progression. AJNR Am J Neuroradiol, 13(1), 181-188. Griessenauer CJ, Adeeb N, Foreman PM, Gupta R, Patel AS, Moore J, Abud TG, Thomas AJ, Ogilvy CS, Baccin CE (2016) Impact of coil packing density and coiling technique on occlusion rates for aneurysms treated with stent-assisted coil embolization. World Neurosurg, 94, 157-166. https://doi.org/10.1016/j.wneu.2016.06.127 Guglielmi G, Viñuela F, Briganti F, Duckwiler G (1992) Carotid-cavernous fistula caused by a ruptured intracavernous aneurysm: Endovascular treatment by electrothrombosis with detachable coils. Neurosurgery, 31(3), 591-596. Gwilliam MN, Hoggard N, Capener D, Singh P, Marzo A, Verma PK, Wilkinson ID (2009) MR derived volumetric flow rate waveforms at locations within the common carotid, internal carotid, and basilar arteries. J Cereb Blood Flow Metab, 29(12), 1975-1982. https://doi.org/10.1038/jcbfm.2009.176 Hanel RA, Cortez GM, Lopes DM, Nelson PK, Siddiqui AH, Jabbour P, Mendes Pereira V, Istvan IS, Zaidat OO, Bettegowda C, Colby GP, Mokin M, Schirmer CM, Hellinger FR, Given C, Krings T, Taussky P, Toth G, Fraser JF, Chen M, Priest R, Kan P, Fiorella D, Frei D, Aagaard-Kienitz B, Diaz O, Malek AM, Cawley CM, Puri AS, Kallmes DF (2023) Prospective study on embolization of intracranial aneurysms with the pipeline device (PREMIER study): 3-year results with the application of a flow diverter specific occlusion classification. J Neurointerv Surg, 15(3), 248-254. https://doi.org/10.1136/neurintsurg-2021-018501 Hasgall PA, Di Gennaro F, Baumgartner C, Neufeld E, Lloyd B, Gosselin M, Payne D, Klingenböck A, Kuster N (2018) IT’IS Database for thermal and electromagnetic parameters of biological tissues. Version 4.0. https://www.scienceopen.com/document?vid=a95fbaa4-efd8-429a-a59e-5e208fea2e45 Herting SM, Ding Y, Boyle AJ, Dai D, Nash LD, Asnafi S, Jakaitis DR, Johnson CR, Graul LM, Yeh C (2019) In vivo comparison of shape memory polymer foam‐coated and bare metal coils for aneurysm occlusion in the rabbit elastase model. J Biomed Mater Res B Appl Biomater, 107(8), 2466-2475. https://doi.org/10.1002/jbm.b.34337 Hodis S, Kargar S, Kallmes DF, Dragomir-Daescu D (2015) Artery length sensitivity in patient-specific cerebral aneurysm simulations. AJNR Am J Neuroradiol, 36(4), 737-743. https://doi.org/10.3174/ajnr.A4179 Hodis S, Uthamaraj S, Smith AL, Dennis KD, Kallmes DF, Dragomir-Daescu D (2012) Grid convergence errors in hemodynamic solution of patient-specific cerebral aneurysms. J Biomech, 45(16), 2907-2913. https://doi.org/10.1016/j.jbiomech.2012.07.030 Horn J, Hwang W, Jessen SL, Keller BK, Miller MW, Tuzun E, Hartman J, Clubb Jr FJ, Maitland DJ (2017) Comparison of shape memory polymer foam versus bare metal coil treatments in an in vivo porcine sidewall aneurysm model. J Biomed Mater Res B Appl Biomater, 105(7), 1892-1905. https://doi.org/10.1002/jbm.b.33725 Hua Y, Oh JH, Kim YB (2015) Influence of parent artery segmentation and boundary conditions on hemodynamic characteristics of intracranial aneurysms. Yonsei Med J, 56(5), 1328-1337. https://doi.org/10.3349/ymj.2015.56.5.1328 Hunter JD (2007) Matplotlib: A 2D graphics environment. Comput Sci Eng, 9(3), 90-95. https://doi.org/10.1109/MCSE.2007.55 Hwang W, Volk BL, Akberali F, Singhal P, Criscione JC, Maitland DJ (2012) Estimation of aneurysm wall stresses created by treatment with a shape memory polymer foam device. Biomech Model Mechanobiol, 11(5), 715-729. https://doi.org/10.1007/s10237-011-0345-8 Keedy A (2006) An overview of intracranial aneurysms. McGill J Med, 9(2), 141-146. Kunkel R, Laurence D, Wang J, Robinson D, Scherrer J, Wu Y, Bohnstedt BN, Chien A, Liu Y, Lee C-H (2018) Synthesis and characterization of bio-compatible shape memory polymers with potential applications to endovascular embolization of intracranial aneurysms. J Mech Behav Biomed Mater, 88, 422-430. https://doi.org/10.1016/j.jmbbm.2018.08.037 Laurence DW, Homburg H, Yan F, Tang Q, Fung K-M, Bohnstedt BN, Holzapfel GA, Lee C-H (2021) A pilot study on biaxial mechanical, collagen microstructural, and morphological characteristics of a resected human intracranial aneurysm tissue. Sci Rep, 11(1), 3525. https://doi.org/10.1038/s41598-021-82991-x Marzo A, Singh P, Reymond P, Stergiopulos N, Patel U, Hose R (2009) Influence of inlet boundary conditions on the local haemodynamics of intracranial aneurysms. Comput Meth Biomech Biomed Eng, 12(4), 431-444. https://doi.org/10.1080/10255840802654335 Menter FR (1994) Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J, 32(8), 1598-1605. https://doi.org/10.2514/3.12149 Moyle KR, Antiga L, Steinman DA (2006) Inlet conditions for image-based CFD models of the carotid bifurcation: Is it reasonable to assume fully developed flow? J Biomech Eng, 128(3), 371-379. https://doi.org/10.1115/1.2187035 Murayama Y, Nien YL, Duckwiler G, Gobin YP, Jahan R, Frazee J, Martin N, Viñuela F (2003) Guglielmi detachable coil embolization of cerebral aneurysms: 11 years’ experience. J Neurosurg, 98(5), 959-966. https://doi.org/10.3171/jns.2003.98.5.0959 Muschenborn AD, Ortega JM, Szafron JM, Szafron DJ, Maitland DJ (2013) Porous media properties of reticulated shape memory polymer foams and mock embolic coils for aneurysm treatment. Biomed Eng Online, 12, 1-13. https://doi.org/10.1186/1475-925X-12-103 Natarajan S, Sekhar LN, Ghodke B, Britz GW, Bhagawati D, Temkin N (2008) Outcomes of ruptured intracranial aneurysms treated by microsurgical clipping and endovascular coiling in a high-volume center. AJNR Am J Neuroradiol, 29(4), 753-759. https://doi.org/10.3174/ajnr.A0895 Ortega JM, Maitland DJ, Wilson T, Tsai W, Savaş Ö, Saloner D (2007) Vascular dynamics of a shape memory polymer foam aneurysm treatment technique. Ann Biomed Eng, 35, 1870-1884. https://doi.org/10.1007/s10439-007-9358-y Ortega JM, Hartman J, Rodriguez JN, Maitland DJ (2013) Virtual treatment of basilar aneurysms using shape memory polymer foam. Ann Biomed Eng, 41(4), 725-743. https://doi.org/10.1007/s10439-012-0719-9 Ouared R, Larrabide I, Brina O, Bouillot P, Erceg G, Yilmaz H, Lovblad KO, Mendes Pereira V (2016) Computational fluid dynamics analysis of flow reduction induced by flow-diverting stents in intracranial aneurysms: A patient-unspecific hemodynamics change perspective. J Neurointerv Surg, 8(12), 1288-1293. https://doi.org/10.1136/neurintsurg-2015-012154 Pineda-Castillo SA, Luo J, Bohnstedt BN, Lee C-H, Liu Y (2020) Shape Memory Polymer Foam with Tunable Properties for Treatment of Intracranial Aneurysm. ASME 2020 International Mechanical Engineering Congress and Exposition, IMECE2020-24291, V005T05A010. https://doi.org/10.1115/IMECE2020-24291 Pineda-Castillo SA, Luo J, Lee H, Bohnstedt BN, Liu Y, Lee C-H (2021) Effects of carbon nanotube infiltration on a shape memory polymer‐based device for brain aneurysm therapeutics: Design and characterization of a Joule‐heating triggering mechanism. Adv Eng Mater, 23(6), 2100322. https://doi.org/10.1002/adem.202100322 Pineda-Castillo SA, Stiles AM, Bohnstedt BN, Lee H, Liu Y, Lee C-H (2022) Shape memory polymer-based endovascular devices: Design criteria and future perspective. Polymers, 14(13), 2526. https://doi.org/10.3390/polym14132526 Pineda-Castillo SA, Cabaniss TL, Aboukeila H, Grady BP, Lee H, Bohnstedt BN, Liu Y, Lee C-H (2023) Toward the development of a shape memory polymer for individualized endovascular therapy of intracranial aneurysms using a 3D-printing/leaching method. Adv Eng Mater, 25(23), 2300683. https://doi.org/10.1002/adem.202300683 Roy D, Milot G, Raymond J (2001) Endovascular treatment of unruptured aneurysms. Stroke, 32(9), 1998-2004. https://doi.org/10.1161/hs0901.095600 Sapareto SA, Dewey WC (1984) Thermal dose determination in cancer therapy. Int J Radiat Oncol Biol Phys, 10(6), 787-800. https://doi.org/10.1016/0360-3016(84)90379-1 Saqr KM, Rashad S, Tupin S, Niizuma K, Hassan T, Tominaga T, Ohta M (2020) What does computational fluid dynamics tell us about intracranial aneurysms? A meta-analysis and critical review. J Cereb Blood Flow Metab, 40(5), 1021-1039. https://doi.org/10.1177/0271678X19854640 Saqr KM, Tupin S, Rashad S, Endo T, Niizuma K, Tominaga T, Ohta M (2020) Physiologic blood flow is turbulent. Sci Rep, 10(1), 15492. https://doi.org/10.1038/s41598-020-72309-8 Schievink WI (1997) Intracranial aneurysms. New Engl J Med, 336(1), 28-40. https://doi.org/10.1056/NEJM199701023360106 Schwartz C, Aster H-C, Al-Schameri R, Müller-Thies-Broussalis E, Griessenauer CJ, Killer-Oberpfalzer M (2018) Microsurgical clipping and endovascular treatment of middle cerebral artery aneurysms in an interdisciplinary treatment concept: Comparison of long-term results. Interv Neuroradiol, 24(6), 608-614. https://doi.org/10.1177/1591019918792231 Shahi K, Ramachandran V (2022) Theoretical and experimental investigation of shape memory polymers programmed below glass transition temperature. Polymers, 14(13), 2753. https://doi.org/10.3390/polym14132753 Soloveva OV, Solovev SA, Vankov YV, Shakurova RZ (2022) Experimental studies of the effective thermal conductivity of polyurethane foams with different morphologies. Processes, 10(11), 2257. https://doi.org/10.3390/pr10112257 Taschner CA, Chapot R, Costalat V, Machi P, Courthéoux P, Barreau X, Berge J, Pierot L, Kadziolka K, Jean B (2018) Second-generation hydrogel coils for the endovascular treatment of intracranial aneurysms: A randomized controlled trial. Stroke, 49(3), 667-674. https://doi.org/10.1161/STROKEAHA.117.018707 Toth G, Cerejo R (2018) Intracranial aneurysms: Review of current science and management. Vasc Med, 23(3), 276-288. https://doi.org/10.1177/1358863X18754693 Valen-Sendstad K, Mardal K-A, Mortensen M, Reif BAP, Langtangen HP (2011) Direct numerical simulation of transitional flow in a patient-specific intracranial aneurysm. J Biomech, 44(16), 2826-2832. https://doi.org/10.1016/j.jbiomech.2011.08.015 Valen-Sendstad K, Bergersen AW, Shimogonya Y, Goubergrits L, Bruening J, Pallares J, Cito S, Piskin S, Pekkan K, Geers AJ, Larrabide I, Rapaka S, Mihalef V, Fu W, Qiao A, Jain K, Roller S, Mardal K-A, Kamakoti R, et al. (2018) Real-world variability in the prediction of intracranial aneurysm wall shear stress: The 2015 International Aneurysm CFD Challenge. Cardiovasc Eng Technol, 9(4), 544-564. https://doi.org/10.1007/s13239-018-00374-2 van Rhoon GC, Samaras T, Yarmolenko PS, Dewhirst MW, Neufeld E, Kuster N (2013) CEM43°C thermal dose thresholds: A potential guide for magnetic resonance radiofrequency exposure levels? Eur J Radiol, 23(8), 2215-2227. https://doi.org/10.1007/s00330-013-2825-y Wagner M, Stenger K (2005) Unruptured intracranial aneurysms: using evidence and outcomes to guide patient teaching. Crit Care Nurs Quart, 28(4), 341-354. https://doi.org/10.1097/00002727-200510000-00007 Wang J, Luo J, Kunkel R, Saha M, Bohnstedt BN, Lee C-H, Liu Y (2019). Development of shape memory polymer nanocomposite foam for treatment of intracranial aneurysms. Mater Lett, 250, 38-41. https://doi.org/10.1016/j.matlet.2019.04.112 Weller HG, Tabor G, Jasak H, Fureby C (1998) A tensorial approach to computational continuum mechanics using object-oriented techniques. Comput Phys, 12(6), 620-631. Xu L, Liang F, Zhao B, Wan J, Liu H (2018) Influence of aging-induced flow waveform variation on hemodynamics in aneurysms present at the internal carotid artery: A computational model-based study. Comput Biol Med, 101, 51-60. https://doi.org/10.1016/j.compbiomed.2018.08.004 Yu L-B, Yang X-J, Zhang Q, Zhang S-S, Zhang Y, Wang R, Zhang D (2019). Management of recurrent intracranial aneurysms after coil embolization: A novel classification scheme based on angiography. J Neurosurg, 131(5), 1455-1461. https://doi.org/10.3171/2018.6.JNS181046 Additional Declarations No competing interests reported. Supplementary Files DescriptionofEnclosedData.docx Cite Share Download PDF Status: Published Journal Publication published 25 Nov, 2024 Read the published version in Biomechanics and Modeling in Mechanobiology → Version 1 posted Editorial decision: Revision requested 22 Sep, 2024 Reviews received at journal 19 Sep, 2024 Reviews received at journal 17 Sep, 2024 Reviewers agreed at journal 09 Sep, 2024 Reviewers agreed at journal 08 Sep, 2024 Reviewers invited by journal 07 Sep, 2024 Editor assigned by journal 06 Sep, 2024 Submission checks completed at journal 01 Sep, 2024 First submitted to journal 01 Sep, 2024 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-5014601","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":357518631,"identity":"cdf2bc40-a438-44ff-a4c8-cf6af7ffefa6","order_by":0,"name":"Tanner Cabaniss","email":"","orcid":"","institution":"The University of Oklahoma Norman","correspondingAuthor":false,"prefix":"","firstName":"Tanner","middleName":"","lastName":"Cabaniss","suffix":""},{"id":357518632,"identity":"668c4af9-b317-4512-a394-1df3e9978de9","order_by":1,"name":"Ryan Bodlak","email":"","orcid":"","institution":"The University of Oklahoma Norman","correspondingAuthor":false,"prefix":"","firstName":"Ryan","middleName":"","lastName":"Bodlak","suffix":""},{"id":357518633,"identity":"58c2597b-2a1d-40ef-bbb2-04712feaffc7","order_by":2,"name":"Yingtao Liu","email":"","orcid":"","institution":"The University of Oklahoma Norman","correspondingAuthor":false,"prefix":"","firstName":"Yingtao","middleName":"","lastName":"Liu","suffix":""},{"id":357518634,"identity":"a27c6bb6-1c57-44aa-959f-e0e2851f19ee","order_by":3,"name":"Geoffrey Colby","email":"","orcid":"","institution":"UCLA Health","correspondingAuthor":false,"prefix":"","firstName":"Geoffrey","middleName":"","lastName":"Colby","suffix":""},{"id":357518635,"identity":"082d716b-630b-4a3e-a45d-63f25b607e39","order_by":4,"name":"Hyowon Lee","email":"","orcid":"","institution":"Purdue University West Lafayette","correspondingAuthor":false,"prefix":"","firstName":"Hyowon","middleName":"","lastName":"Lee","suffix":""},{"id":357518636,"identity":"f2f3340d-d08f-4bf7-9619-1123958f0fea","order_by":5,"name":"Bradley Bohnstedt","email":"","orcid":"","institution":"Indiana University School of Medicine Indianapolis","correspondingAuthor":false,"prefix":"","firstName":"Bradley","middleName":"","lastName":"Bohnstedt","suffix":""},{"id":357518637,"identity":"f8437bfe-64a5-4500-baa9-fdc87f365021","order_by":6,"name":"Rinaldo Garziera","email":"","orcid":"","institution":"University of Parma","correspondingAuthor":false,"prefix":"","firstName":"Rinaldo","middleName":"","lastName":"Garziera","suffix":""},{"id":357518638,"identity":"660bcbb6-2be9-4386-9c95-8f8c28e08d44","order_by":7,"name":"Gerhard Holzapfel","email":"","orcid":"","institution":"Graz University of Technology","correspondingAuthor":false,"prefix":"","firstName":"Gerhard","middleName":"","lastName":"Holzapfel","suffix":""},{"id":357518639,"identity":"f2aa8978-656a-42dd-8bf1-05e43ec3edec","order_by":8,"name":"Chung-Hao Lee","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA20lEQVRIiWNgGAWjYDACCSBmbAAS7ED8ASJmQKQWZiA9g2QtzDzEaOGf3WP28OsOOzkGZuaDj23b7PIZ2Ju3SeC15M4Zc2PZM8nGDMxsyca5bcmWDTzHyvBqMZDIMZOWbGNObGDmMZPObTtgwAAUIUZLPVAL//ffliAt8m8Ia5H82HYYZAsbMyPYFh78WiRupJVJM7YdN2ZjZjOW7DmXbMDGk1ZsgU8L/4zkbZI/26rl+NmbH374UWZnwM9+eOMNfFpAABwdbDAeGx6VcMD4gxhVo2AUjIJRMHIBAJQdOswT6cbuAAAAAElFTkSuQmCC","orcid":"","institution":"University of California Riverside","correspondingAuthor":true,"prefix":"","firstName":"Chung-Hao","middleName":"","lastName":"Lee","suffix":""}],"badges":[],"createdAt":"2024-09-02 00:36:21","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-5014601/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-5014601/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1007/s10237-024-01910-x","type":"published","date":"2024-11-25T15:57:31+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":66947310,"identity":"d792f8a7-f29f-4327-b4f6-824f67cdd5f2","added_by":"auto","created_at":"2024-10-18 09:53:06","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":1130995,"visible":true,"origin":"","legend":"\u003cp\u003ePatient-specific vasculature and aneurysms: (a) segmented vasculature for each geometry with inflow (red arrows) and outflow (blue arrows); (b) magnified view of the patient-specific aneurysms (sizes are to scale relative to each other; B = bifurcation aneurysm; S = side-wall aneurysm).\u003c/p\u003e","description":"","filename":"Figure1.png","url":"https://assets-eu.researchsquare.com/files/rs-5014601/v1/a8587ca1a098d7e9707014e0.png"},{"id":66944537,"identity":"6d985965-da83-4ae6-a8b4-426b9164695b","added_by":"auto","created_at":"2024-10-18 09:37:07","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":663516,"visible":true,"origin":"","legend":"\u003cp\u003eSMPf device geometry for each of the six patient-specific ICA cases examined. Devices are shown with slight transparency to represent the full neck region (views are shown at size for clarity but are not a consistent scale).\u003c/p\u003e","description":"","filename":"Figure214.png","url":"https://assets-eu.researchsquare.com/files/rs-5014601/v1/c5893c4b0f67e331aa96aa32.png"},{"id":66944531,"identity":"50248ee4-a199-4dfa-bbbf-4b804b0d5b59","added_by":"auto","created_at":"2024-10-18 09:37:06","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":3116416,"visible":true,"origin":"","legend":"\u003cp\u003eSchematic of the compression of a representative SMPf device in its original state at 50% and 90% radial compressive strain.\u003c/p\u003e","description":"","filename":"Figure3.png","url":"https://assets-eu.researchsquare.com/files/rs-5014601/v1/18175c2bffdba2c275ccc969.png"},{"id":66946387,"identity":"ad81a609-4832-47d4-a823-4da5c2ad50d4","added_by":"auto","created_at":"2024-10-18 09:45:06","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":1525356,"visible":true,"origin":"","legend":"\u003cp\u003eGeometry shortening for Case B2 in the original, flow simulation and heat transfer simulation states.\u003c/p\u003e","description":"","filename":"Figure4.png","url":"https://assets-eu.researchsquare.com/files/rs-5014601/v1/d9e4061b667bf3bca7925c74.png"},{"id":66944535,"identity":"4ce9daca-820c-4a4f-ba5a-c25d4df687f6","added_by":"auto","created_at":"2024-10-18 09:37:06","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":3230874,"visible":true,"origin":"","legend":"\u003cp\u003eComparison of WSS values to representative results from the 2015 International Aneurysm CFD Challenge (Valen-Sendstad \u003cem\u003eet al.\u003c/em\u003e 2018): (a) WSS distributions for the six highly experienced teams (units are in Pa, and layout rearranged); (b) WSS distribution calculated for mean flow rate (\u003cem\u003etop\u003c/em\u003e) and systolic flow rate (\u003cem\u003ebottom\u003c/em\u003e).\u003c/p\u003e","description":"","filename":"Figure5.png","url":"https://assets-eu.researchsquare.com/files/rs-5014601/v1/1acfca6f0faaf0c07d23ce05.png"},{"id":66944530,"identity":"657d5fb1-d110-4b02-8d9e-d929175735dd","added_by":"auto","created_at":"2024-10-18 09:37:06","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":1660038,"visible":true,"origin":"","legend":"\u003cp\u003eComparison of intra-aneurysmal velocities before and after treatment in the central plane of each of the 6 ICAs studied.\u003c/p\u003e","description":"","filename":"Figure6.png","url":"https://assets-eu.researchsquare.com/files/rs-5014601/v1/6ae698132da7a73d1c9611a3.png"},{"id":66944533,"identity":"9b0de37f-e9af-4a15-af6d-c1f69870a448","added_by":"auto","created_at":"2024-10-18 09:37:06","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":1576015,"visible":true,"origin":"","legend":"\u003cp\u003eTemperature isosurface of 37.5 °C for steady-state heating in patient-specific aneurysm geometries. Note the inclusion of a second isosurface (\u003cem\u003eT\u003c/em\u003e= 39 °C) for Case B2 and the indicated downstream isosurface contact for Case B1.\u003c/p\u003e","description":"","filename":"Figure7.png","url":"https://assets-eu.researchsquare.com/files/rs-5014601/v1/241b9c9c4aed425986bbbc37.png"},{"id":66946389,"identity":"d2994174-962e-4266-bcf2-60b43d2b26ad","added_by":"auto","created_at":"2024-10-18 09:45:06","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":1046147,"visible":true,"origin":"","legend":"\u003cp\u003eTemperature at the vessel walls over 1 s of the expansion and heating process, starting with a 40% expanded foam: (a) 40 °C stimulus; (b) 45 °C stimulus.\u003c/p\u003e","description":"","filename":"Figure8.png","url":"https://assets-eu.researchsquare.com/files/rs-5014601/v1/9b9f24c22e5c898c747481e0.png"},{"id":66947311,"identity":"87412332-5f21-49e0-8f2d-920119814e90","added_by":"auto","created_at":"2024-10-18 09:53:06","extension":"png","order_by":9,"title":"Figure 9","display":"","copyAsset":false,"role":"figure","size":646500,"visible":true,"origin":"","legend":"\u003cp\u003eAnalysis of potential thermal damage caused by the shape recovery heating stimulus: (a) time history of the maximum wall temperature of the simulated stimulus for the triggering of the SMPf shape recovery and the resulting cooling; (b) representative temperature \u003cem\u003eT\u003c/em\u003e and time \u003cem\u003et\u003c/em\u003e combinations related to CEM\u003csub\u003e43\u003c/sub\u003e = 10 min (blue = no thermal damage; red = thermal damage).\u003c/p\u003e","description":"","filename":"Figure9.png","url":"https://assets-eu.researchsquare.com/files/rs-5014601/v1/222845de93644a75d5954753.png"},{"id":70388701,"identity":"5d91ca6c-29df-4a5d-87af-c4ac8f6de006","added_by":"auto","created_at":"2024-12-02 17:26:59","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":14397760,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-5014601/v1/db7b6e47-c152-4fb2-a709-da31e0f6a654.pdf"},{"id":66944526,"identity":"ff6e7dbb-2409-4961-93eb-cc9e3b8b1fa2","added_by":"auto","created_at":"2024-10-18 09:37:06","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":59881,"visible":true,"origin":"","legend":"","description":"","filename":"DescriptionofEnclosedData.docx","url":"https://assets-eu.researchsquare.com/files/rs-5014601/v1/ada4bf44aeb60be6100f18cd.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"CFD investigations of a shape-memory polymer foam-based endovascular embolization device for the treatment of intracranial aneurysms","fulltext":[{"header":"1. Introduction","content":"\u003cp\u003eIntracranial aneurysms (ICAs), which affect approximately 4% of the American population (Keedy \u003cem\u003eet al.\u003c/em\u003e 2006), are saccular protrusions that arise form focal dilation of the arterial wall (Toth \u003cem\u003eet al.\u003c/em\u003e 2018). It has been found that 85% of ICAs develop within \u003cem\u003ethe circle of Willis\u003c/em\u003e, the complex central cranial vasculature (Keedy \u003cem\u003eet al.\u003c/em\u003e 2006). Although there are different classifications of ICAs, namely saccular, fusiform, dissecting, and mycotic, this study focuses exclusively on saccular aneurysms as they account for approximately 90% of ICAs (Keedy \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2006\u003c/span\u003e). The primary risk of ICAs is the possibility of rupture, which can further lead to subarachnoid hemorrhage: a bleeding in the intracranial cavity that leads to a variety of life-threatening conditions such as severe brain damage, long-term morbidity and/or mortality (Schievink \u003cem\u003eet al.\u003c/em\u003e 1997). While the causative physiological mechanisms for ICA rupture are not yet fully understood, current literature presents correlations with multiple factors, including age, hypertension, smoking history, alcohol abuse, hypercholesterolemia, and carotid stenosis (Wagner \u003cem\u003eet al.\u003c/em\u003e 2005).\u003c/p\u003e \u003cp\u003eThe success rate of ICA treatment is mainly measured using the Raymond-Roy Occlusion Classification scale. This classification evaluates the degree of occlusion of the aneurysm sac into three distinct classes: Class I \u0026ndash; complete obliteration; Class II \u0026ndash; residual neck; and Class III \u0026ndash; residual aneurysm (Roy et al. \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2001\u003c/span\u003e). Complete obliteration, also called complete occlusion, is defined as the cessation of intra-aneurysmal blood circulation, with the majority of blood flowing through the parent artery as if the aneurysm were not present.\u003c/p\u003e \u003cp\u003eGuglielmi detachable coils (GDCs), developed in the 1990s (Guglielmi et al. \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e1992\u003c/span\u003e), have long been considered the gold standard for \u003cem\u003eendovascular\u003c/em\u003e treatment of ICAs. According to a retrospective study by Gonda et al. (\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2014\u003c/span\u003e) over 2,509 unruptured aneurysms, while GDCs have a 30-day postoperative mortality rate lower than that of highly invasive microsurgical clipping (1.1% vs. 2.3%), in the long term they have unsatisfactorily lower rates of complete occlusion (i.e., Raymond-Roy Occlusion Classification Class I). Another retrospective study by Murayama et al. (\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2003\u003c/span\u003e) reported that traditional coiling had an astonishing complete closure rate of 55%, compared to complete occlusion rates of \u0026gt;\u0026thinsp;90% with microsurgical clipping (Natarajan et al. \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2008\u003c/span\u003e; Schwartz et al. \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Choi et al. \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eTo address this critical shortcoming in aneurysm occlusion, new endovascular devices have been developed over the past 2\u0026ndash;3 decades. These devices consist of direct improvements to GDCs through the introduction of hydrogel and shape memory polymer (SMP)-based coatings as well as novel endovascular devices such as stents and flow diverters (Alderazi et al. \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Taschner et al. \u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Herting et al. \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). While stent-based devices have demonstrated improved complete occlusion rates of up to 83.3% (Hanel et al. \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), flow diversion via stents is only useful in the treatment of sidewall aneurysms, so bifurcation aneurysms continue to be treated primarily via coil-based techniques. Although these results provide significant room for improvement, new endovascular technology may be the answer to achieving comparable rates of complete occlusion to microsurgical clipping procedures without the associated procedural risks of a craniotomy.\u003c/p\u003e \u003cp\u003eThis emerging technology is based on the use of shape memory polymer foams (SMPf) \u0026ndash; a subgroup of SMP materials that possess a high degree of porosity while maintaining the shape recovery properties of typical shape memory polymers (Pineda-Castillo et al. \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). These properties allow the SMPf to be compressed into a microcatheter, implanted into an aneurysm cavity, and subsequently expand and recover to its original programmed geometry via shape recovery activation, thereby occluding the aneurysm sac. The SMPf have been fabricated with both closed and open cell architectures, the latter enabling higher compression ratios required for endovascular techniques (Wang et al. \u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). Moreover, recent work has used discrete lattice architectures to enable further degrees of structural optimization (Pineda-Castillo et al. \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). To trigger the shape recovery of SMPf through Joule heating, conductive polymer coatings and nanomaterials such as carbon nanotubes were introduced to increase the electrical and thermal conductivity of the fabricated SMPf (Pineda-Castillo et al. \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Pineda-Castillo et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Horn et al. (\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2017\u003c/span\u003e) conducted \u003cem\u003ein vivo\u003c/em\u003e studies on porcine side-wall aneurysms that demonstrated the promise and efficacy of this therapeutic technique. However, several inquiries remain regarding the optimal geometric considerations of these foams. These salient challenges revolve around crucial aspects of foam design such as mechanical stability, optimal configurations to reduce intra-aneurysmal residual circulation, and promotion of platelet aggregation to improve coagulation processes.\u003c/p\u003e \u003cp\u003eOne of the earliest applications of computational fluid dynamics (CFD) in aneurysm research came from Gonzalez et al. (\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e1992\u003c/span\u003e), who used pulsatile flow and a non-Newtonian fluid model to simulate blood flow through an aneurysm in the middle cerebral artery. By 2004, only eleven additional articles had been published on the topic, but research efforts have increased rapidly since then, as a study by Saqr et al. (\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) found 795 publications on ICA hemodynamics using CFD as the primary \u003cem\u003ein silico\u003c/em\u003e modeling tool. While there is a wealth of information about hemodynamics in ICAs and the effect of other devices (e.g., stents, coils, etc.), only limited efforts have been made to investigate the impact of SMPf-based devices for \u003cem\u003epersonalized\u003c/em\u003e treatment of ICAs.\u003c/p\u003e \u003cp\u003eThe aim of this study is therefore to use CFD to examine the hemodynamic influence of SMPf implantation in \u0026ldquo;patient-specific\u0026rdquo; aneurysm geometries. Specifically, our aim is to quantify the reduction in intra-aneurysmal blood flow caused by the SMPf device and to evaluate the convective heat transfer to the surrounding arterial tissues during initiation of shape recovery as well as successive device expansion. The results of this work could also lead us to a deeper understanding of the relevant intricacies of SMPf devices for patient-specific endovascular ICA therapeutics.\u003c/p\u003e \u003cp\u003eThe remainder of this paper is organized as follows: in Section \u003cspan refid=\"Sec2\" class=\"InternalRef\"\u003e2\u003c/span\u003e, we describe the development of a CFD-based approach to analyze the hemodynamics of untreated ICAs and aneurysms implanted with an SMPf-based device, as well as the heat transfer mechanisms associated with activation of shape recovery. In Section \u003cspan refid=\"Sec10\" class=\"InternalRef\"\u003e3\u003c/span\u003e, we present the results of the models and validate their efficacy using an example from the established literature. In Section \u003cspan refid=\"Sec15\" class=\"InternalRef\"\u003e4\u003c/span\u003e, we will further assess the advantages and limitations of the patient-specific SMPf device when compared to other established endovascular techniques, followed by some concluding remarks.\u003c/p\u003e"},{"header":"2. Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003e2.1 Geometry Preparation\u003c/h2\u003e \u003cp\u003eEleven anonymized, patient-specific 3D rotational angiography datasets were acquired from the University of Oklahoma Health Sciences Center in 2018 and 2019 (IRB #7932). Of these datasets, nine had sufficient resolution and were selected for image segmentation. These were successively narrowed down to five unique geometries from the collected data; three bifurcation and two sidewall aneurysms as well as another geometry from the literature (Valen-Sendstad et al. \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2018\u003c/span\u003e) for our numerical modeling/simulation. The four segmented geometries excluded from this study were structurally similar to the geometries presented. Image segmentation was performed in Amira (Thermo Fisher Scientific, Waltham, MA) with \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:0.108\\times\\:0.108\\times\\:0.5\\)\u003c/span\u003e\u003c/span\u003e mm rectangular prism voxels. The segmented patient-specific ICA lumen surface was then imported into Meshmixer (Autodesk, San Rafael, CA) as an STL file and slightly smoothed to remove localized peaks and kinks due to noise in the image data. The vessels were then cut normal to their cross sections to form inlets and outlets, while the outlets were maintained at approximately twice the vessel diameter once their respective directions were established, and the inlets were cut at least seven vessel diameters from the aneurysm sac, as suggested by the literature convention (Moyle et al. \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2006\u003c/span\u003e; Hua et al. \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). These cuts were made to reduce the computational costs while maintaining validity of the flow characteristics throughout the studied system (Hodis et al. \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). The full extent of the modeled vasculature for each geometry and the respective more detailed views of the ICAs are shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e, while the main geometric characteristics of the modeled aneurysms are summarized in Table\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e. Here, the bifurcation aneurysm B1 was taken from the 2015 International Aneurysm CFD Challenge (Valen-Sendstad et al. \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2018\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eCharacteristics of the six patient-specific aneurysms studied (B1-B4 and S1, S2).\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"7\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eLocation\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003el\u003c/em\u003e\u003csub\u003emax\u003c/sub\u003e (mm)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eAR\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eSR\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cem\u003eV\u003c/em\u003e (mm\u003csup\u003e3\u003c/sup\u003e)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cem\u003eA\u003c/em\u003e\u003csub\u003eostium\u003c/sub\u003e (mm\u003csup\u003e2\u003c/sup\u003e)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eB1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMCA bifurc.\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e7.7\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e1.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e76\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e28.2\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eB2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMCA bifurc.\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e7.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e2.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e50\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e12.4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eB3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eMCA bifurc.\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e5.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.6\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e1.9\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e41\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e21.0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eB4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eBasilar apex\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e11.1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e2.5\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e313\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e36.1\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eS1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSup. hypophyseal\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e9.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e2.3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e2.0\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e167\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e24.4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eS2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eAnt. choroidal\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e2.9\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e0.8\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e6.6\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"BlockQuote\"\u003e \u003cp\u003eMCA birfurc. = middle cerebral artery bifurcation; Sup. hypophyseal\u0026thinsp;=\u0026thinsp;superior hypophyseal artery; Ant. choroidal\u0026thinsp;=\u0026thinsp;anterior choroidal artery; \u003cem\u003el\u003c/em\u003e\u003csub\u003emax\u003c/sub\u003e = facet diameter; AR\u0026thinsp;=\u0026thinsp;aspect ratio (sac depth: neck width); SR\u0026thinsp;=\u0026thinsp;size ratio (sac depth: parent artery diameter); \u003cem\u003eV\u003c/em\u003e\u0026thinsp;=\u0026thinsp;volume of aneurysm sac; \u003cem\u003eA\u003c/em\u003e\u003csub\u003eostium\u003c/sub\u003e = interface area between aneurysm sac and main vessel.\u003c/p\u003e \u003c/div\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003e2.2 Simulation Parameters\u003c/h2\u003e \u003cp\u003eSince one of the goals of this study is to support the development of individually-optimized SMPf devices (Kunkel et al. \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2018\u003c/span\u003e;Wang et al. \u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e2019\u003c/span\u003e), steady-state simulations were performed in OpenFOAM 8 (Weller et al. \u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e1998\u003c/span\u003e), assuming that the greatest degree of blood flow into the SMPf device occurs during systole. Moreover, because steady-state simulations with systolic flow tend to overestimate the velocity field (Geers et al. \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2014\u003c/span\u003e), this consideration provided an upper bound on the worst-case flow that the SMPf device and the treated aneurysm would experience. Accordingly, steady-state solutions were achieved using the SIMPLEC (Semi-Implicit Method for Pressure Linked Equations-Consistent) algorithm and intra-aneurysmal blood flow residuals were monitored to ensure convergence. The meshes were constructed predominantly from hexahedron cells, while prismatic and tetrahedral cells were incorporated along the boundaries with a growth factor of 1.2 in the 12\u0026ndash;20 inflation layers, resulting in mesh sizes from \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:4.4\\times\\:{10}^{6}\\)\u003c/span\u003e\u003c/span\u003e to \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:7.9\\times\\:{10}^{6}\\)\u003c/span\u003e\u003c/span\u003e cells. The results of the mesh selection process and convergence study are presented in Sections \u003cspan refid=\"Sec8\" class=\"InternalRef\"\u003e2.6\u003c/span\u003e and \u003cspan refid=\"Sec11\" class=\"InternalRef\"\u003e3.1\u003c/span\u003e.\u003c/p\u003e \u003cp\u003eIn these simulations, the blood flow was modeled to be an incompressible Newtonian fluid with a dynamic viscosity \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\mu\\:=3.65\\times\\:{10}^{-3}\\)\u003c/span\u003e\u003c/span\u003e Pa\u0026middot;s, a density \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\rho\\:=1050\\)\u003c/span\u003e\u003c/span\u003e kg/m\u003csup\u003e3\u003c/sup\u003e and a kinematic viscosity \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\upsilon\\:=3.5\\times\\:{10}^{-6}\\)\u003c/span\u003e\u003c/span\u003e m\u003csup\u003e2\u003c/sup\u003e/s (Hasgall et al. \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). No material was assigned to the vessel walls, but they were treated as rigid according to convention, with a no-slip boundary condition applied along the walls (Berg et al. \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2019\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eNext, the inlet boundary conditions were defined using the peak systolic volumetric flow rate (VFR) values estimated from population-averaged flow waveforms for the ICA and basilar artery reported in the literature for a normal, young cohort: 9.26 mL/s and 5.63 mL/s, respectively (Horn et al. \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Generalized waveforms for older adults were avoided due to the age range of the present cohort was unknown, although both datasets would likely have produced similar flow patterns and wall shear stresses (WSS) (Xu et al. \u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). According to Chnafa et al. (\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2018\u003c/span\u003e), using the splitting method outlined above, VFRs were prescribed on all but one outlet, typically the smallest. The unassigned outlet was left in the zero pressure condition while a zero-gradient pressure boundary was applied to all other outlets.\u003c/p\u003e \u003cp\u003eFinally, there is increasing evidence that transitional or turbulent flow occurs in physiological blood flow (Saqr et al. \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2020\u003c/span\u003e) and that turbulent flow has been observed in the aneurysm space after systole (Valen-Sendstad et al. \u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Accordingly, the \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:k-{\\omega\\:}\\)\u003c/span\u003e\u003c/span\u003e shear stress transport turbulence model (Menter \u003cem\u003eet al.\u003c/em\u003e 1994) was adopted to realize the turbulence characteristics in the simulated intracranial circulation environment.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003e2.3 Simulations of an SMPf-Based Device Implementation\u003c/h2\u003e \u003cp\u003eTreatment efficacy was evaluated by virtually implanting SMPf devices into each patient-specific aneurysm geometry. First, the geometry of each SMPf device was generated by manually excising the aneurysm sac from the surrounding vasculature (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e). The devices were then idealized and represented the maximum possible foam volume for the specified geometries. Manufacturing considerations or defects/discrepancies were not taken into account.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eBecause in this study we focused on providing general design suggestions rather than seeking to evaluate a patient-specific case in a clinical setting, we chose a continuum porous media approach to model the SMP foam. The porous medium specifications considered in this study followed the previous work by Muschenborn et al. (\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2013\u003c/span\u003e) with a porosity \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\varPhi\\:=98\\%\\)\u003c/span\u003e\u003c/span\u003e, permeability \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:k=3.25\\times\\:{10}^{-9}\\sim1.24\\times\\:{10}^{-8}\\)\u003c/span\u003e\u003c/span\u003e m\u003csup\u003e2\u003c/sup\u003e, and form factor \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:C=2720-7260\\)\u003c/span\u003e\u003c/span\u003e m\u003csup\u003e\u0026minus;1\u003c/sup\u003e. These parameter ranges were verified by evaluating the pressure gradient of the flow through the porous media model at different flow velocities and comparing them with the original study (Muschenborn et al. \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). For the results presented in this work, the lowest permeability model was selected because it demonstrated the greatest performance in reducing intra-aneurysmal flow.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003e2.4 Simulations of Heat Transfer Process During Initial Shape Recovery that Triggers Heating\u003c/h2\u003e \u003cp\u003eTo quantify the thermal effects of device deployment and the heating associated with the initiation of shape recovery, an idealized, rigid-walled solid cylinder with a diameter of 1 mm was used. The use of a compressed solid model was justified because the pores of the closed cell SMPf were almost completely closed upon compression (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e). The length of the cylinder models varied between ICA cases, but in all cases ranged from the aneurysm neck to a 1 mm offset from the vessel walls perpendicular to the neck, centered in the intra-aneurysmal space.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe potential thermally induced tissue damage can be described by the cumulative equivalent minutes at 43\u0026deg;C (CEM\u003csub\u003e43\u003c/sub\u003e), defined by\u003cdiv id=\"Equ1\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ1\" name=\"EquationSource\"\u003e\n$$\\:CE{M}_{43}\\left(t\\right)={\\int\\:}_{{t}_{0}}^{{t}_{f}}{R}^{\\left(43-T\\left(t\\right)\\right)}\\text{d}t$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e1\u003c/div\u003e\u003c/div\u003e,\u003c/p\u003e \u003cp\u003ewhere \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{t}_{0}\\)\u003c/span\u003e\u003c/span\u003e is the initial time, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{t}_{f}\\)\u003c/span\u003e\u003c/span\u003e is the end time of exposure in minutes, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:R\\)\u003c/span\u003e\u003c/span\u003e is a constant of 0.5 for \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:T\\left(t\\right)\\ge\\:43\\:℃\\)\u003c/span\u003e\u003c/span\u003e, 0.25 for \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:39\\:℃\\le\\:T\\left(t\\right)\u0026lt;43\\:℃\\)\u003c/span\u003e\u003c/span\u003e, and 0.0 for \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:T\\left(t\\right)\u0026lt;39\\:℃\\)\u003c/span\u003e\u003c/span\u003e (Sapareto \u003cem\u003eet al.\u003c/em\u003e 1984). According to van Rhoon et al. (\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2013\u003c/span\u003e), irreversible damage to brain tissues occurs with as little as CEM\u003csub\u003e43\u003c/sub\u003e=10. Assuming a constant exposure temperature, \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:40\\:℃\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:45\\:℃\\)\u003c/span\u003e\u003c/span\u003e stimuli would result in tissue damage after \u0026gt;\u0026thinsp;600 min and 2.5 min, respectively. Therefore, the SMPf device in this study was considered to have a glass transition temperature \u003cem\u003eT\u003c/em\u003e\u003csub\u003eg\u003c/sub\u003e of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:39\\:℃\\)\u003c/span\u003e\u003c/span\u003e (Pineda-Castillo et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Pineda-Castillo et al. (\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2021\u003c/span\u003e) have shown that the SMPf device can heat from well below body temperature to its transition temperature in less than a minute while completely occluding an \u003cem\u003ein vitro\u003c/em\u003e aneurysm sac in the same amount of time. Thus, a fixed temperature of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:40\\:℃\\)\u003c/span\u003e\u003c/span\u003e was used in the heating portion of this numerical study, just above the \u003cem\u003eT\u003c/em\u003e\u003csub\u003eg\u003c/sub\u003e of the SMPf material. A zero-heat flux (perfectly insulated) boundary condition was imposed on the vessel to model the worst-case temperature rise (Ortega et al. \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2007\u003c/span\u003e). The flow through the inlet was set to body temperature (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:37\\:℃\\)\u003c/span\u003e\u003c/span\u003e), with the blood being assigned a specific heat capacity of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{c}_{p}=3600\\)\u003c/span\u003e\u003c/span\u003e J/kg\u0026middot;\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:℃\\)\u003c/span\u003e\u003c/span\u003e and a thermal conductivity of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{k}_{t}=0.52\\)\u003c/span\u003e\u003c/span\u003e W/m\u0026middot;\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:℃\\)\u003c/span\u003e\u003c/span\u003e (Hasgall et al. \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). The inlet flow for this study was modeled at the diastolic VFR of 2.98 mL/s and 1.77 mL/s for the internal carotid artery and basilar artery, respectively (Gwilliam et al. \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2009\u003c/span\u003e) to minimize convection from the region and to simulate the highest possible temperature rise on the walls. The simulations were considered steady-state and included the Reynolds-averaged Navier-Stokes equations, although the slow flow rates significantly reduced the instabilities observed in the previous systolic flow simulations.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003e2.5 Simulations of Heat Transfer During Device Expansion\u003c/h2\u003e \u003cp\u003eNext, transient heat transfer simulations were performed on a representative patient geometry, Case B2, which demonstrated the highest level of heat transfer during the initial heating stage of the deployment. To limit the computational expense of the additional solution steps, the domain upstream of the ICA was shortened to approximately three diameters of the parent artery (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e) and remeshed with a cell size of 0.20 mm to generate approximately 41,000 cells.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eTo provide a more realistic flow condition than the previous diastolic study, outflow VFRs were calculated from the mean ICA velocity. All other boundary conditions were retained. Although the simulation was transient, the flow was still modeled as steady-state rather than pulsatile, using a time step of 1 ms and a PIMPLE algorithm, an iterative solver which is a combination of PISO (Pressure Implicit with Splitting of Operator) and SIMPLE (Semi-Implicit Method for Pressure-Linked Equations). For computational efficiency reasons, the total time steps were limited to 1,000. The geometry was constructed as a partially expanded foam with an offset of 0.5 mm from the vessel wall. As in the flow simulations, the model was defined as a porous medium with a permeability of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:5.7\\times\\:{10}^{-9}\\)\u003c/span\u003e\u003c/span\u003e m\u003csup\u003e2\u003c/sup\u003e and form factor of 7260 m\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e. The entire zone of the porous medium was assigned a temperature of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:40\\:℃\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:45\\:℃\\)\u003c/span\u003e\u003c/span\u003e for the two respective shape recovery-induced heating processes to provide an optimal and worst-case heating scenario. The surrounding fluid volume was modeled similarly to that described in Section \u003cspan refid=\"Sec6\" class=\"InternalRef\"\u003e2.4\u003c/span\u003e.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003e2.6 Simulation Verification\u003c/h2\u003e \u003cp\u003eA mesh independence study was performed on the aneurysm B1 geometry (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e), with simulations covering a range of 0.05\u0026ndash;0.16 mm cell sizes (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:6.87\\times\\:{10}^{5}-1.25\\times\\:{10}^{7}\\)\u003c/span\u003e\u003c/span\u003e cells). The number of cells was approximately doubled for each successive cell size, without the finest mesh containing 1.5 times as many cells as the second finest mesh. The maximum and average flow velocities for 1,000 points evenly distributed throughout the aneurysm domain were calculated in ParaView (Ahrens et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2005\u003c/span\u003e). To reduce the influence of insignificant errors/discrepancies in the simulation, the slowest 1% of the sampled points, corresponding to a velocity of 0.03 m/s, were excluded from the analyses.\u003c/p\u003e \u003cp\u003eTo further verify the simulation parameters and results, the mean internal carotid artery flow rate of 5.11 mL/s (Gwilliam et al. \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2009\u003c/span\u003e) was scaled and applied to a 0.08 mm mesh. This enabled a direct comparison of WSS values between this study and the time-averaged WSS results of the six most experienced teams at the 2015 International Aneurysm CFD Challenge (Valen-Sendstad et al. \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). However, despite their experience, there was no clear consensus regarding the WSS magnitude numerically determined between the six teams: i.e., four teams achieved comparable values and the other two teams reported drastically different values. While a direct quantitative and numerical comparison may not have a high degree of validity in this case, qualitative contrasts can be drawn based on the WSS distributions presented by the teams as they are holistically more similar.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003e2.7 Simulation Metrics and Data Analysis\u003c/h2\u003e \u003cp\u003eFor the flow velocity studies, the WSS distribution was beneficial to qualitatively verify the validity of the flow within the model. Intra-aneurysmal velocity is also an important metric to evaluate since reducing flow is a necessity for rebuilding and recovering arterial tissue. Intra-aneurysmal velocity was evaluated for maximum and mean quantities for the interior of the aneurysm space. For the heat transfer studies, the temperature for the intra-aneurysmal space and surrounding vessels was assessed, and visualized using isosurfaces and contours. For all studies presented, ParaView (Ahrens et al. \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2005\u003c/span\u003e) was utilized to extract the above metrics from the simulated data and generate visualizations. The Python library Matplotlib (Hunter \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2007\u003c/span\u003e) was used to generate the CEM\u003csub\u003e43\u003c/sub\u003e plots.\u003c/p\u003e \u003c/div\u003e"},{"header":"3. Results","content":"\u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003e3.1 Study on Mesh Size\u003c/h2\u003e \u003cp\u003eThe mesh size was determined to be \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:s=0.08\\)\u003c/span\u003e\u003c/span\u003e mm as this was the maximum cell size with less than 5% mean absolute error (MAE) compared to the control (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:s=0.05\\)\u003c/span\u003e\u003c/span\u003e mm), resulting in a cell density of approximately 1,700-1,800 cells/mm\u003csup\u003e3\u003c/sup\u003e, which is consistent with the cell densities presented in the literature (Marzo et al. \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2009\u003c/span\u003e; Hodis et al. \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2012\u003c/span\u003e; Aranda \u003cem\u003eet al.\u003c/em\u003e 2018). The mesh convergence study results are summarized in Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e. Further qualitative comparisons of the WSS contours with those of the 2015 International Aneurysm CFD Challenge (Valen-Sendstad et al. \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2018\u003c/span\u003e) showed good agreement in the WSS distributions along the aneurysm surface (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eDetailed results of the mesh independence study. (Velocity values are calculated over the entire aneurysm domain.)\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"6\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u003cem\u003es\u003c/em\u003e (mm)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCells (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\times\\:{10}^{3}\\)\u003c/span\u003e\u003c/span\u003e)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eIA cells (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\times\\:{10}^{3}\\)\u003c/span\u003e\u003c/span\u003e)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003e\u003cem\u003eu\u003c/em\u003e\u003csub\u003emax\u003c/sub\u003e (m/s)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cem\u003eu\u003c/em\u003e\u003csub\u003eavg\u003c/sub\u003e (m/s)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eMAE (%)\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e0.05\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e12450\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e938\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.106\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.494\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u0026ndash;\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e0.06\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e7780\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e612\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.103\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.494\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e1.5\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e0.08\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e4160\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e329\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.097\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.494\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e4.4\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e0.10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e2230\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e194\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.086\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.494\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e9.0\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e0.13\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e1150\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e106\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.061\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.490\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e15.9\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e0.16\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e687\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e67\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e1.054\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.490\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e18.3\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003e \u003cdiv class=\"BlockQuote\"\u003e \u003cp\u003e \u003cem\u003es\u003c/em\u003e\u0026thinsp;=\u0026thinsp;mesh size; IA\u0026thinsp;=\u0026thinsp;intra-aneurysmal; \u003cem\u003eu\u003c/em\u003e\u003csub\u003emax\u003c/sub\u003e = maximum velocity; \u003cem\u003eu\u003c/em\u003e\u003csub\u003eavg\u003c/sub\u003e = average velocity; MAE\u0026thinsp;=\u0026thinsp;mean absolute error.\u003c/p\u003e \u003c/div\u003e \u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003e3.2 Simulations of Untreated and Treated ICAs\u003c/h2\u003e \u003cp\u003eAs discussed in Section \u003cspan refid=\"Sec5\" class=\"InternalRef\"\u003e2.3\u003c/span\u003e, porous media models were used to define the permeability of the SMPf device. For the results presented here, the lowest permeability of the foam was adopted, defined by a permeability of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:5.7\\times\\:{10}^{-9}\\)\u003c/span\u003e\u003c/span\u003e m\u003csup\u003e2\u003c/sup\u003e and a form factor of 7260 m\u003csup\u003e\u0026minus;\u0026thinsp;1\u003c/sup\u003e. In each case, treatment of the aneurysm showed significant decrease in blood flow into the aneurysm cavity. Overall, in each case, the maximum flow velocity was reduced by \u0026gt;\u0026thinsp;40% and the average flow velocity was reduced by up to 95%, as shown in Table\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e. The metric used to quantify this reduction was the velocity reduction ratio VRR, i.e.\u003cdiv id=\"Equ2\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ2\" name=\"EquationSource\"\u003e\n$$\\:VRR=\\frac{{u}_{\\text{p}\\text{r}\\text{e}}-{u}_{\\text{p}\\text{o}\\text{s}\\text{t}}}{{u}_{\\text{p}\\text{r}\\text{e}}}\\:,$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e2\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere \u003cem\u003eu\u003c/em\u003e is the predicted fluid velocity and the subscripts, pre and post, denote the respective velocities before and after device implantation, respectively.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab3\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 3\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eSimulation results of intracranial flow velocity before and after SMPf implantation. All velocity values are given in m/s.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"7\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e\u0026nbsp;\u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cem\u003eu\u003c/em\u003e\u003csub\u003emax,pre\u003c/sub\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cem\u003eu\u003c/em\u003e\u003csub\u003emax,post\u003c/sub\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003e\u003cem\u003eVRR\u003c/em\u003e\u003csub\u003emax\u003c/sub\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cem\u003eu\u003c/em\u003e\u003csub\u003eavg,pre\u003c/sub\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003e\u003cem\u003eu\u003c/em\u003e\u003csub\u003eavg,post\u003c/sub\u003e\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003e\u003cem\u003eVRR\u003c/em\u003e\u003csub\u003eavg\u003c/sub\u003e\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eB1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e1.107\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e0.456\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e58.8%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.494\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e0.054\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e89.0%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eB2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.746\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e0.387\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e48.2%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.172\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e0.017\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e90.0%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eB3\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.775\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e0.333\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e57.1%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.328\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e0.033\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e90.1%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eB4\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e1.180\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e0.488\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e58.7%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.333\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e0.022\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e93.3%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eS1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.997\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e0.428\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e57.1%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.219\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e0.011\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e95.0%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eS2\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e0.825\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e0.329\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e60.1%\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.403\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e0.048\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c7\"\u003e \u003cp\u003e88.0%\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003ctfoot\u003e \u003ctr\u003e\u003ctd colspan=\"7\"\u003e\u003cem\u003eu\u003c/em\u003e\u0026thinsp;=\u0026thinsp;predicted fluid velocity; pre\u0026thinsp;=\u0026thinsp;before device implantation;\u003c/td\u003e\u003c/tr\u003e \u003ctr\u003e\u003ctd colspan=\"7\"\u003epost\u0026thinsp;=\u0026thinsp;after device implantation; \u003cem\u003eVRR\u003c/em\u003e\u0026thinsp;=\u0026thinsp;velocity reduction ratio;\u003c/td\u003e\u003c/tr\u003e \u003c/tfoot\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eFurthermore, this residual flow reduction can be visualized using the velocity contours of the central aneurysmal plane and inlet vessel for both the untreated and treated cases (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e). More profound differences were found in the bifurcation cases (B1, B2, B3, and B4), where high-velocity inflow jets were observed prior to treatment. In addition to reducing flow within the intra-aneurysmal space, device implantation appeared to have stabilized the downstream flow, as indicated by the smoother velocity contours in these areas. While the flow reduction results are promising for high efficacy of a patient-specific SMPf-based device, the deployment of these devices introduces several biocompatibility concerns, with thermal damage to surrounding arterial tissue being the greatest concern and requiring further investigation.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003e3.3 Heat Transfer Simulations for Initial Heating of the SMPf Device \u0026ndash; Triggering Shape Recovery\u003c/h2\u003e \u003cp\u003eAs described previously, arterial tissue damage occurs near CEM\u003csub\u003e43\u003c/sub\u003e of 10 minutes (van Rhoon et al. \u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). Therefore, the thermal effects of the device deployment process on the aneurysm and artery wall tissue must be carefully investigated to identify possible thermal damage with respect to the appropriate activation temperature for shape recovery. The first case considered was the initial thermal activation heat transfer, as at this point the convection patterns from the SMPf device to the surrounding fluid became visible. Temperature isosurfaces were plotted to observe heat convection and the distribution of the initial activation heat to the surrounding fluid (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e). The selected isosurface was chosen at \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:37.5\\:℃\\)\u003c/span\u003e\u003c/span\u003e just \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:0.5\\:℃\\)\u003c/span\u003e\u003c/span\u003e above body temperature, to maximize interpretable data of the underlying heat transfer process and impact.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eIn general, the \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:37.5\\:\\text{℃}\\)\u003c/span\u003e\u003c/span\u003e isosurfaces (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e) did not contact any walls within the aneurysm space or downstream vasculature, which was expected given the minimal temperature gradient generated during device activation at \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:40\\:℃\\)\u003c/span\u003e\u003c/span\u003e. However, there are two exceptions. In case B1, a thin stream of heated fluid impinges on the wall at a downstream bifurcation, although heat transfer thereby is negligible. In addition, given the low fluid velocities compared to the other cases (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e), it was not surprising that Case B2 was somewhat more concerning as it was found that convection within the dome was significantly hindered. This allows fluid at a higher temperature level to come into contact with the vessel wall at the apex of the treated aneurysm. This resulted in an open isosurface at \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:37.5\\:℃\\)\u003c/span\u003e\u003c/span\u003e, indicating local wall temperatures above body temperature. Nevertheless, with the \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:39\\:℃\\)\u003c/span\u003e\u003c/span\u003e isosurface superimposed, the temperatures on the vessel wall were still lower than those of the device (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:40\\:℃\\)\u003c/span\u003e\u003c/span\u003e). Because Case B2 showed the greatest concern for exposing vessel walls to potentially degenerative temperatures, this case was utilized to further explore the impact of heat transfer during the subsequent expansion phase of the device.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section2\"\u003e \u003ch2\u003e3.4 Transient Heat Transfer During SMPf Device Expansion and Assessment of Tissue Damage\u003c/h2\u003e \u003cp\u003eThe maximum wall temperature was observed to increase to \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:40\\:℃\\)\u003c/span\u003e\u003c/span\u003e within 1 s at both stimulus temperatures, which is due to the relatively slow flow within the aneurysm dome and the small distance between the SMP foam and the tissue wall (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e). It was found that the blood flow pattern is such that only certain regions of the aneurysm sac reach elevated temperatures, while almost all of the heated blood from those areas drains into only one of the bifurcation vessels, significantly increasing its wall temperature.\u003c/p\u003e \u003cp\u003eThe time histories for the heating and cooling process of the SMPf device shape recovery can be combined to calculate the CEM\u003csub\u003e43\u003c/sub\u003e dosage for Case B2. Assuming a total heating time of 1 minute (Pineda-Castillo et al. \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e2021\u003c/span\u003e), the maximum vessel wall temperature curves were generated throughout the entire shape recovery process (Fig.\u0026nbsp;\u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e9\u003c/span\u003e(a)). Since the blood flow temperature for the \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:40\\:℃\\)\u003c/span\u003e\u003c/span\u003e stimulus was found to never exceed \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:40\\:℃\\)\u003c/span\u003e\u003c/span\u003e, the CEM\u003csub\u003e43\u003c/sub\u003e measurement is almost negligible for the entire heating process (i.e., 0.02 min). In contrast, the \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:45\\:℃\\)\u003c/span\u003e\u003c/span\u003e stimulus resulted in a worst-case CEM\u003csub\u003e43\u003c/sub\u003e of 4.9 min. Even in this (over-heating) worst-case scenario, there would be no thermal damage to the surrounding arteries/aneurysm wall (van Rhoon et al. \u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). The general relationship between time and temperature for CEM\u003csub\u003e43\u003c/sub\u003e is shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig9\" class=\"InternalRef\"\u003e9\u003c/span\u003e(b).\u003c/p\u003e \u003c/div\u003e"},{"header":"4. Discussion","content":"\u003cp\u003eIn this study, we examined the viability of a patient-specific SMPf device for the treatment of ICAs. We explored the hemodynamical effects of device implantation on intra-aneurysmal flow as well as the implications of a thermal heating mechanism within the aneurysm cavity, taking into account both the initial heating stages and transient heating during device expansion. Below we will discuss the potential advantages and disadvantages regarding the efficacy of the proposed SMPf patient-specific embolic device and the validity of the numerical studies presented here.\u003c/p\u003e \u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003e4.1 Flow Simulations\u003c/h2\u003e \u003cp\u003eIn this work, we presented baseline and treatment flow simulations for six different patient-specific geometries, with the SMPf device modeled as a porous medium. In the untreated cases, we found an average maximum velocity in the aneurysm space of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:0.94\\pm\\:0.18\\)\u003c/span\u003e\u003c/span\u003e m/s and an overall average velocity of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:0.32\\pm\\:0.12\\)\u003c/span\u003e\u003c/span\u003e m/s across all patient geometries. After treatment, we found an average maximum velocity of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:0.40\\pm\\:0.07\\)\u003c/span\u003e\u003c/span\u003e m/s and an overall average velocity of \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:0.031\\pm\\:0.017\\)\u003c/span\u003e\u003c/span\u003e m/s resulting in an approximate reduction of 60% and 90% in the maximum and average flow velocities, respectively. While the low permeability foam used in our numerical simulations showed promising results, in real practical applications higher permeability may be required to meet the required compression constraints as the device compresses radially and is inserted into an endovascular micro-catheter with an inner diameter of \u0026lt;\u0026thinsp;1 mm for use in clinical setting.\u003c/p\u003e \u003cp\u003eIn our internal pilot studies, not presented here, we found that higher permeability foams do not perform as well as the idealized low permeability foam presented in this manuscript. A study by Ouared et al. (\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2016\u003c/span\u003e) found that with vascular stents, even a 30% reduction in intra-aneurysmal flow can promote thrombosis in the aneurysm sac, which can ultimately lead to effective aneurysm occlusion. Therefore, even with higher permeability foam, there is still evidence for the viability of the SMPf device presented. In addition, Ortega et al. (\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2013\u003c/span\u003e) evaluated an SMPf device with higher permeability (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:k\\approx\\:2\\times\\:{10}^{-8}\\)\u003c/span\u003e\u003c/span\u003e m\u003csup\u003e2\u003c/sup\u003e) and found a flow reduction in the intra-aneurysmal region furthest from the aneurysm neck up to 95%. Their results are consistent with our results presented in this work.\u003c/p\u003e \u003cp\u003eWhile endovascular coil embolization remains one of the most prominent methods for the endovascular (minimally invasive) treatment of ICAs, the SMPf device offers potential solutions to some of the more common limitations of GDC coiling therapeutic. Most importantly, the patient-specific SMPf aims to maximally reduce the residual volume of the aneurysm, which is a cause of recurrence in patients treated with coils and requires repeated surgery (Griessenauer et al. \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). Another shortcoming of coil embolization that the proposed SMPf device aims to improve is post-treatment compaction. Because the soft platinum wires used in coiling are delicate and flexibly deformable, clinical evidence has shown that blood inflow into the aneurysm sac can gradually compact the packed coils, further leading to residual volume near the aneurysm neck and causing recurrence (Yu et al. \u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). On the other hand, the SMPf device would have a higher degree of structural integrity due to its highly interconnected cellular structure.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec17\" class=\"Section2\"\u003e \u003ch2\u003e4.2 Device Heating and Deployment\u003c/h2\u003e \u003cp\u003eHeat transfer from SMPf device deployment and initiation of shape recovery is highly dependent on the specific aneurysm geometry and flow speed (Fig.\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e). For example, transient analysis of Case B2, which showed the highest level of heat transfer from the device to the aneurysmal and arterial walls, indicated that the wall temperature would approach that of the expanding SMPf within seconds. While this was most clearly observed in Case B2, it is likely that this condition applies to all study cases during further stages of device expansion as fluid is forced from the intra-aneurysmal space and convection capacity is reduced. However, the potential thermal damage to tissue from this heating process depends heavily on the temperature required for the device to expand (e.g., a thermal activation threshold very close to body temperature would still not cause damage to surrounding tissue regardless of the duration of the activation/expansion process). Furthermore, the heterogeneity of the arterial tissue composing the aneurysm sac, both in terms of microstructural composition and tissue thickness (Laurence et al., \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2021\u003c/span\u003e), may result in a non-uniform response to the heating stimulus. This can potentially lead to local thermal damage below the previously mentioned stimulus thresholds.\u003c/p\u003e \u003cp\u003eOther possible solutions to the heating problem include reducing the time required for full expansion, which is also related to the material properties of the SMPf, as the discrepancy between the temperature provided to the device and the \u003cem\u003eT\u003c/em\u003e\u003csub\u003eg\u003c/sub\u003e of the device regulates the rate at which expansion occurs (Kunkel et al. \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Shahi \u003cem\u003eet al.\u003c/em\u003e 2022). This positive correlation between shape recovery rate and stimulus temperature is related to the increased kinetic energy within the polymer system at higher temperatures, which allows the polymer chains to move more freely and reestablish their initial geometric state more quickly (Davidson \u003cem\u003eet al.\u003c/em\u003e 2015). Furthermore, although the boundary conditions presented in this study provided the stimulus temperature (\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:40\\:℃\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:45\\:℃\\)\u003c/span\u003e\u003c/span\u003e) on the exterior of the SMPf model, in practice with a central heating element, the low thermal conductivity of a SMPf (Soloveva et al. \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) would likely have a strong thermal gradient that reduces heat transfer to the blood and surrounding tissues during all stages of deployment.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec18\" class=\"Section2\"\u003e \u003ch2\u003e4.3 Study Limitations and Future Work\u003c/h2\u003e \u003cp\u003eAlthough the results presented in this study are promising, there are notable limitations that should also be discussed. First, no experimental validation was performed to verify the accuracy of the presented results. Second, the porous models presented in the study were all isotropic in nature, ignoring the potential effects of anisotropic foam, which are more common in practice due to the foaming process (Ortega et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2013\u003c/span\u003e). Furthermore, the steady-state assumption neglected potential benefits of the foam-based device, e.g., previous research studies have shown that SMPf can filter out high-frequency velocity fluctuations in the aneurysm space (Ortega et al. \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2013\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eIn short, the SMPf device promises to reduce intra-aneurysmal flow velocities and offers potential solutions to the limitations encountered with traditional coil-based methods. Now we will discuss some possible drawbacks arising from the thermal activation mechanism for SMPf deployment. The most urgent extension of this work is the construction of an \u003cem\u003ein vitro\u003c/em\u003e flow loop with an integrated particle image velocimetry system for experimental validation of the studies reported herein. Moreover, extensions of computational simulations with pulsating, time-dependent conditions can provide insights into flow behavior caused by inertial effects not foreseen in this work. An analysis of the deployment of the device would be helpful in determining the viability of the underlying concepts of a SMPf device (e.g., examining the stresses that the expanding foam exerts on the vessel wall in the event of device misalignment, similar to the work of Hwang et al. \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2012\u003c/span\u003e). Further analysis of the thermal properties of the specific SMPf material and inclusion in the aforementioned heat transfer simulations would provide a more realistic description of the plausibly harmful heat generation of the device during expansion in the intra-aneurysmal space.\u003c/p\u003e \u003cp\u003eAnother possible extension is the inclusion of a fluid-structure interaction (FSI) approach to model blood flow-induced deformations in the aneurysm and foam. With such an approach, the surfaces are not rigid and can respond to fluid motion. Accordingly, FSI modeling is more computationally expensive and requires a higher degree of boundary conditions and model data since both the fluid and solid domains are modeled. If properly performed, FSI modeling would offer a more accurate physiological representation of the blood flow patterns in the ICA and the effects of the proposed SMPf device on it.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec19\" class=\"Section2\"\u003e \u003ch2\u003e4.4 Conclusion\u003c/h2\u003e \u003cp\u003eIn this work, we presented results on the hemodynamic effects and possible heat transfer impairments of a SMPf device for the treatment of ICAs. We have demonstrated a significant flow reduction in the intra-aneurysmal space, effects of initial shape recovery on various aneurysm geometries, and, in the worst case, transient heat transfer to the surrounding tissue. Although these results were promising and prospective solutions have been developed to address concerns about device heat transfer during shape recovery induction, future work is still required to address the full complexity of treating ICAs and develop a highly effective treatment device. However, this work adds to the growing body of research assessing the effectiveness of SMPf devices for the endovascular treatment of ICAs, thereby providing additional evidence for further investigation and the design and development of personalized SMPf-based endovascular devices.\u003c/p\u003e \u003c/div\u003e"},{"header":"Declarations","content":"\u003ch2\u003econflict of interest\u003c/h2\u003e \u003cp\u003eThe authors of this paper have no financial or personal relationships with other people or organizations that could inappropriately influence (bias) our work.\u003c/p\u003e \u003ch2\u003eUse of human subject data\u003c/h2\u003e \u003cp\u003eDeidentified angiographic imaging data from patients were collected and used in this study. The study protocol for human subjects was reviewed and approved by the Institution Review Board of the University of Oklahoma Health Sciences Center (IRB #7932).\u003c/p\u003e \u003c/p\u003e\u003ch2\u003eAuthor Contribution\u003c/h2\u003e\u003cp\u003eT.L.C.: Methodology, Figure Preparation, Manuscript Writing/EditingR.B.: Conceptualization, Methodology, Figure Preparation, Manuscript WritingY.L.: Conceptualization, Resources, Student/Researcher Mentoring, Dissemination, Manuscript Review/EditingG.P.C.: Conceptualization, Manuscript Review/EditingH.L.: Conceptualization, Manuscript Review/EditingB.N.B.: Conceptualization, Methodology, Project Administration, Manuscript Review/EditingR.G.: Conceptualization, Resources, Manuscript Review/EditingG.A.H.: Conceptualization, Methodology, Manuscript Review/EditingC.-H.L.: Conceptualization, Methodology, Project Administration, Resources, Student/Researcher Mentoring, Dissemination, Manuscript Review/Editing\u003c/p\u003e\u003ch2\u003eAcknowledgements\u003c/h2\u003e \u003cp\u003eThis project was supported by grants from the Oklahoma Center for the Advancement of Science and Technology (OCAST, grant no. HR18-002), and the Oklahoma Shared Translational and Clinical Research (OSCTR, NIGMS U54GM104938). In addition, this publication reports research supported by the National Heart, Lung, And Blood Institute of the National Institutes of Health under Award Number R01HL159475, and by the University of Parma through the action Bando di Ateneo 2021 per la ricerca co-funded by MUR-Italian Ministry of Universities and Research \u0026ndash; D.M. 737/2021-PNR-PNRR \u0026ndash; NextGenerationEU. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.\u003c/p\u003e\u003ch2\u003eData Availability\u003c/h2\u003e\u003cp\u003eData is provided within the manuscript or supplementary information files.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n \u003cli\u003eAhrens J, Geveci B, Law C, Hansen C, Johnson C (2005) 36-paraview: An end-user tool for large-data visualization. The Visualization Handbook, 717, 50038-50031. https://doi.org/10.1016/B978-012387582-2/50038-1\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eAlderazi YJ, Shastri D, Kass-Hout T, Prestigiacomo CJ, Gandhi CD (2014) Flow diverters for intracranial aneurysms. Stroke Res. Treat, 2014, 415653. https://doi.org/10.1155/2014/415653\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eAranda A, Valencia A (2018) Study on the relationship between wall shear stress and aspect ratio of cerebral aneurysms with different pressure differences using CFD simulations. J Mech Med Biol, 18(05), 1850055. https://doi.org/10.1142/S0219519418500550\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eBerg P, Saalfeld S, Vo\u0026szlig; S, Beuing O, Janiga G (2019) A review on the reliability of hemodynamic modeling in intracranial aneurysms: Why computational fluid dynamics alone cannot solve the equation. Neurosurg Focus, 47(1), E15. https://doi.org/10.3171/2019.4.FOCUS19181\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eChnafa C, Brina O, Pereira V, Steinman D (2018) Better than nothing: A rational approach for minimizing the impact of outflow strategy on cerebrovascular simulations. AJNR Am J Neuroradiol, 39(2), 337-343. https://doi.org/10.3174/ajnr.A5484\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eChoi JH, Lee KS, Kim B-S, Shin YS (2020) Treatment outcomes of large and giant intracranial aneurysms according to various treatment modalities. Acta Neurochir, 162, 2745-2752. https://doi.org/10.1007/s00701-020-04540-1\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eDavidson JD, Goulbourne NC (2015) Microscopic mechanisms of the shape memory effect in crosslinked polymers. Smart Mater Struct, 24(5), 055014. https://doi.org/10.1088/0964-1726/24/5/055014\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eGeers A, Larrabide I, Morales H, Frangi A (2014) Approximating hemodynamics of cerebral aneurysms with steady flow simulations. J Biomech, 47(1), 178-185. https://doi.org/10.1016/j.jbiomech.2013.09.033\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eGonda DD, Khalessi AA, McCutcheon BA, Marcus LP, Noorbakhsh A, Chen CC, Chang DC, Carter BS (2014) Long-term follow-up of unruptured intracranial aneurysms repaired in California. J Neurosurg, 120(6), 1349-1357. https://doi.org/10.3171/2014.3.JNS131159 \u0026nbsp;\u003c/li\u003e\n \u003cli\u003eGonzalez CF, Cho YI, Ortega HV, Moret J (1992) Intracranial aneurysms: Flow analysis of their origin and progression. AJNR Am J Neuroradiol, 13(1), 181-188.\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eGriessenauer CJ, Adeeb N, Foreman PM, Gupta R, Patel AS, Moore J, Abud TG, Thomas AJ, Ogilvy CS, Baccin CE (2016) Impact of coil packing density and coiling technique on occlusion rates for aneurysms treated with stent-assisted coil embolization. World Neurosurg, 94, 157-166. https://doi.org/10.1016/j.wneu.2016.06.127 \u0026nbsp;\u003c/li\u003e\n \u003cli\u003eGuglielmi G, Vi\u0026ntilde;uela F, Briganti F, Duckwiler G (1992) Carotid-cavernous fistula caused by a ruptured intracavernous aneurysm: Endovascular treatment by electrothrombosis with detachable coils. Neurosurgery, 31(3), 591-596.\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eGwilliam MN, Hoggard N, Capener D, Singh P, Marzo A, Verma PK, Wilkinson ID (2009) MR derived volumetric flow rate waveforms at locations within the common carotid, internal carotid, and basilar arteries. J Cereb Blood Flow Metab, 29(12), 1975-1982. https://doi.org/10.1038/jcbfm.2009.176\u003c/li\u003e\n \u003cli\u003eHanel RA, Cortez GM, Lopes DM, Nelson PK, Siddiqui AH, Jabbour P, Mendes Pereira V, Istvan IS, Zaidat OO, Bettegowda C, Colby GP, Mokin M, \u0026nbsp;Schirmer CM, Hellinger FR, Given C, Krings T, Taussky P, Toth G, Fraser JF, Chen M, Priest R, Kan P, Fiorella D, Frei D, Aagaard-Kienitz B, Diaz O, Malek AM, Cawley CM, Puri AS, Kallmes DF (2023) Prospective study on embolization of intracranial aneurysms with the pipeline device (PREMIER study): 3-year results with the application of a flow diverter specific occlusion classification. J Neurointerv Surg, 15(3), 248-254. https://doi.org/10.1136/neurintsurg-2021-018501\u003c/li\u003e\n \u003cli\u003eHasgall PA, Di Gennaro F, Baumgartner C, Neufeld E, Lloyd B, Gosselin M, Payne D, Klingenb\u0026ouml;ck A, Kuster N (2018) IT\u0026rsquo;IS Database for thermal and electromagnetic parameters of biological tissues. Version 4.0. https://www.scienceopen.com/document?vid=a95fbaa4-efd8-429a-a59e-5e208fea2e45\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eHerting SM, Ding Y, Boyle AJ, Dai D, Nash LD, Asnafi S, Jakaitis DR, Johnson CR, Graul LM, Yeh C (2019) In vivo comparison of shape memory polymer foam‐coated and bare metal coils for aneurysm occlusion in the rabbit elastase model. J Biomed Mater Res B Appl Biomater, 107(8), 2466-2475. https://doi.org/10.1002/jbm.b.34337\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eHodis S, Kargar S, Kallmes DF, Dragomir-Daescu D (2015) Artery length sensitivity in patient-specific cerebral aneurysm simulations. AJNR Am J Neuroradiol, 36(4), 737-743. https://doi.org/10.3174/ajnr.A4179 \u0026nbsp;\u003c/li\u003e\n \u003cli\u003eHodis S, Uthamaraj S, Smith AL, Dennis KD, Kallmes DF, Dragomir-Daescu D (2012) Grid convergence errors in hemodynamic solution of patient-specific cerebral aneurysms. J Biomech, 45(16), 2907-2913. https://doi.org/10.1016/j.jbiomech.2012.07.030\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eHorn J, Hwang W, Jessen SL, Keller BK, Miller MW, Tuzun E, Hartman J, Clubb Jr FJ, Maitland DJ (2017) Comparison of shape memory polymer foam versus bare metal coil treatments in an in vivo porcine sidewall aneurysm model. J Biomed Mater Res B Appl Biomater, 105(7), 1892-1905. https://doi.org/10.1002/jbm.b.33725\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eHua Y, Oh JH, Kim YB (2015) Influence of parent artery segmentation and boundary conditions on hemodynamic characteristics of intracranial aneurysms. Yonsei Med J, 56(5), 1328-1337. https://doi.org/10.3349/ymj.2015.56.5.1328 \u0026nbsp;\u003c/li\u003e\n \u003cli\u003eHunter JD (2007) Matplotlib: A 2D graphics environment. Comput Sci Eng, 9(3), 90-95. https://doi.org/10.1109/MCSE.2007.55 \u0026nbsp;\u003c/li\u003e\n \u003cli\u003eHwang W, Volk BL, Akberali F, Singhal P, Criscione JC, Maitland DJ (2012) Estimation of aneurysm wall stresses created by treatment with a shape memory polymer foam device. Biomech Model Mechanobiol, 11(5), 715-729. https://doi.org/10.1007/s10237-011-0345-8 \u0026nbsp;\u003c/li\u003e\n \u003cli\u003eKeedy A (2006) An overview of intracranial aneurysms. McGill J Med, 9(2), 141-146.\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eKunkel R, Laurence D, Wang J, Robinson D, Scherrer J, Wu Y, Bohnstedt BN, Chien A, Liu Y, Lee C-H (2018) Synthesis and characterization of bio-compatible shape memory polymers with potential applications to endovascular embolization of intracranial aneurysms. J Mech Behav Biomed Mater, 88, 422-430. https://doi.org/10.1016/j.jmbbm.2018.08.037\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eLaurence DW, Homburg H, Yan F, Tang Q, Fung K-M, Bohnstedt BN, Holzapfel GA, Lee C-H (2021) A pilot study on biaxial mechanical, collagen microstructural, and morphological characteristics of a resected human intracranial aneurysm tissue. Sci Rep, 11(1), 3525. https://doi.org/10.1038/s41598-021-82991-x\u003c/li\u003e\n \u003cli\u003eMarzo A, Singh P, Reymond P, Stergiopulos N, Patel U, Hose R (2009) Influence of inlet boundary conditions on the local haemodynamics of intracranial aneurysms. Comput Meth Biomech Biomed Eng, 12(4), 431-444. https://doi.org/10.1080/10255840802654335\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eMenter FR (1994) Two-equation eddy-viscosity turbulence models for engineering applications. AIAA J, 32(8), 1598-1605. https://doi.org/10.2514/3.12149\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eMoyle KR, Antiga L, Steinman DA (2006) Inlet conditions for image-based CFD models of the carotid bifurcation: Is it reasonable to assume fully developed flow? J Biomech Eng, 128(3), 371-379. https://doi.org/10.1115/1.2187035 \u0026nbsp;\u003c/li\u003e\n \u003cli\u003eMurayama Y, Nien YL, Duckwiler G, Gobin YP, Jahan R, Frazee J, Martin N, Vi\u0026ntilde;uela F (2003) Guglielmi detachable coil embolization of cerebral aneurysms: 11 years\u0026rsquo; experience. J Neurosurg, 98(5), 959-966. https://doi.org/10.3171/jns.2003.98.5.0959\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eMuschenborn AD, Ortega JM, Szafron JM, Szafron DJ, Maitland DJ (2013) Porous media properties of reticulated shape memory polymer foams and mock embolic coils for aneurysm treatment. Biomed Eng Online, 12, 1-13. https://doi.org/10.1186/1475-925X-12-103\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eNatarajan S, Sekhar LN, Ghodke B, Britz GW, Bhagawati D, Temkin N (2008) Outcomes of ruptured intracranial aneurysms treated by microsurgical clipping and endovascular coiling in a high-volume center. AJNR Am J Neuroradiol, 29(4), 753-759. https://doi.org/10.3174/ajnr.A0895\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eOrtega JM, Maitland DJ, Wilson T, Tsai W, Savaş \u0026Ouml;, Saloner D (2007) Vascular dynamics of a shape memory polymer foam aneurysm treatment technique. Ann Biomed Eng, 35, 1870-1884. https://doi.org/10.1007/s10439-007-9358-y\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eOrtega JM, Hartman J, Rodriguez JN, Maitland DJ (2013) Virtual treatment of basilar aneurysms using shape memory polymer foam. Ann Biomed Eng, 41(4), 725-743. https://doi.org/10.1007/s10439-012-0719-9 \u0026nbsp;\u003c/li\u003e\n \u003cli\u003eOuared R, Larrabide I, Brina O, Bouillot P, Erceg G, Yilmaz H, Lovblad KO, Mendes Pereira V (2016) Computational fluid dynamics analysis of flow reduction induced by flow-diverting stents in intracranial aneurysms: A patient-unspecific hemodynamics change perspective. J Neurointerv Surg, 8(12), 1288-1293. https://doi.org/10.1136/neurintsurg-2015-012154\u0026nbsp;\u003c/li\u003e\n \u003cli\u003ePineda-Castillo SA, Luo J, Bohnstedt BN, Lee C-H, Liu Y (2020) Shape Memory Polymer Foam with Tunable Properties for Treatment of Intracranial Aneurysm. ASME 2020 International Mechanical Engineering Congress and Exposition, IMECE2020-24291, V005T05A010. https://doi.org/10.1115/IMECE2020-24291\u0026nbsp;\u003c/li\u003e\n \u003cli\u003ePineda-Castillo SA, Luo J, Lee H, Bohnstedt BN, Liu Y, Lee C-H (2021) Effects of carbon nanotube infiltration on a shape memory polymer‐based device for brain aneurysm therapeutics: Design and characterization of a Joule‐heating triggering mechanism. Adv Eng Mater, 23(6), 2100322. https://doi.org/10.1002/adem.202100322\u0026nbsp;\u003c/li\u003e\n \u003cli\u003ePineda-Castillo SA, Stiles AM, Bohnstedt BN, Lee H, Liu Y, Lee C-H (2022) Shape memory polymer-based endovascular devices: Design criteria and future perspective. Polymers, 14(13), 2526. https://doi.org/10.3390/polym14132526\u0026nbsp;\u003c/li\u003e\n \u003cli\u003ePineda-Castillo SA, Cabaniss TL, Aboukeila H, Grady BP, Lee H, Bohnstedt BN, Liu Y, Lee C-H (2023) Toward the development of a shape memory polymer for individualized endovascular therapy of intracranial aneurysms using a 3D-printing/leaching method. Adv Eng Mater, 25(23), 2300683. https://doi.org/10.1002/adem.202300683 \u0026nbsp;\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eRoy D, Milot G, Raymond J (2001) Endovascular treatment of unruptured aneurysms. Stroke, 32(9), 1998-2004. https://doi.org/10.1161/hs0901.095600\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eSapareto SA, Dewey WC (1984) Thermal dose determination in cancer therapy. Int J Radiat Oncol Biol Phys, 10(6), 787-800. https://doi.org/10.1016/0360-3016(84)90379-1\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eSaqr KM, Rashad S, Tupin S, Niizuma K, Hassan T, Tominaga T, Ohta M (2020) What does computational fluid dynamics tell us about intracranial aneurysms? A meta-analysis and critical review. J Cereb Blood Flow Metab, 40(5), 1021-1039. https://doi.org/10.1177/0271678X19854640\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eSaqr KM, Tupin S, Rashad S, Endo T, Niizuma K, Tominaga T, Ohta M (2020) Physiologic blood flow is turbulent. Sci Rep, 10(1), 15492. https://doi.org/10.1038/s41598-020-72309-8\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eSchievink WI (1997) Intracranial aneurysms. New Engl J Med, 336(1), 28-40. https://doi.org/10.1056/NEJM199701023360106\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eSchwartz C, Aster H-C, Al-Schameri R, M\u0026uuml;ller-Thies-Broussalis E, Griessenauer CJ, Killer-Oberpfalzer M (2018) Microsurgical clipping and endovascular treatment of middle cerebral artery aneurysms in an interdisciplinary treatment concept: Comparison of long-term results. Interv Neuroradiol, 24(6), 608-614. https://doi.org/10.1177/1591019918792231\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eShahi K, Ramachandran V (2022) Theoretical and experimental investigation of shape memory polymers programmed below glass transition temperature. Polymers, 14(13), 2753. https://doi.org/10.3390/polym14132753\u003c/li\u003e\n \u003cli\u003eSoloveva OV, Solovev SA, Vankov YV, Shakurova RZ (2022) Experimental studies of the effective thermal conductivity of polyurethane foams with different morphologies. Processes, 10(11), 2257. https://doi.org/10.3390/pr10112257\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eTaschner CA, Chapot R, Costalat V, Machi P, Courth\u0026eacute;oux P, Barreau X, Berge J, Pierot L, Kadziolka K, Jean B (2018) Second-generation hydrogel coils for the endovascular treatment of intracranial aneurysms: A randomized controlled trial. Stroke, 49(3), 667-674. https://doi.org/10.1161/STROKEAHA.117.018707\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eToth G, Cerejo R (2018) Intracranial aneurysms: Review of current science and management. Vasc Med, 23(3), 276-288. https://doi.org/10.1177/1358863X18754693\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eValen-Sendstad K, Mardal K-A, Mortensen M, Reif BAP, Langtangen HP (2011) Direct numerical simulation of transitional flow in a patient-specific intracranial aneurysm. J Biomech, 44(16), 2826-2832. https://doi.org/10.1016/j.jbiomech.2011.08.015\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eValen-Sendstad K, Bergersen AW, Shimogonya Y, Goubergrits L, Bruening J, Pallares J, Cito S, Piskin S, Pekkan K, Geers AJ, Larrabide I, Rapaka S, Mihalef V, Fu W, Qiao A, Jain K, Roller S, Mardal K-A, Kamakoti R, et al. (2018) Real-world variability in the prediction of intracranial aneurysm wall shear stress: The 2015 International Aneurysm CFD Challenge. Cardiovasc Eng Technol, 9(4), 544-564. https://doi.org/10.1007/s13239-018-00374-2 \u0026nbsp;\u0026nbsp;\u003c/li\u003e\n \u003cli\u003evan Rhoon GC, Samaras T, Yarmolenko PS, Dewhirst MW, Neufeld E, Kuster N (2013) CEM43\u0026deg;C thermal dose thresholds: A potential guide for magnetic resonance radiofrequency exposure levels? Eur J Radiol, 23(8), 2215-2227. https://doi.org/10.1007/s00330-013-2825-y \u0026nbsp;\u003c/li\u003e\n \u003cli\u003eWagner M, Stenger K (2005) Unruptured intracranial aneurysms: using evidence and outcomes to guide patient teaching. Crit Care Nurs Quart, 28(4), 341-354. https://doi.org/10.1097/00002727-200510000-00007\u003c/li\u003e\n \u003cli\u003eWang J, Luo J, Kunkel R, Saha M, Bohnstedt BN, Lee C-H, Liu Y (2019). Development of shape memory polymer nanocomposite foam for treatment of intracranial aneurysms. Mater Lett, 250, 38-41. https://doi.org/10.1016/j.matlet.2019.04.112\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eWeller HG, Tabor G, Jasak H, Fureby C (1998) A tensorial approach to computational continuum mechanics using object-oriented techniques. Comput Phys, 12(6), 620-631.\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eXu L, Liang F, Zhao B, Wan J, Liu H (2018) Influence of aging-induced flow waveform variation on hemodynamics in aneurysms present at the internal carotid artery: A computational model-based study. Comput Biol Med, 101, 51-60. https://doi.org/10.1016/j.compbiomed.2018.08.004\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eYu L-B, Yang X-J, Zhang Q, Zhang S-S, Zhang Y, Wang R, Zhang D (2019). Management of recurrent intracranial aneurysms after coil embolization: A novel classification scheme based on angiography. J Neurosurg, 131(5), 1455-1461. https://doi.org/10.3171/2018.6.JNS181046\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"biomechanics-and-modeling-in-mechanobiology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"bmmb","sideBox":"Learn more about [Biomechanics and Modeling in Mechanobiology](http://link.springer.com/journal/10237)","snPcode":"10237","submissionUrl":"https://submission.nature.com/new-submission/10237/3","title":"Biomechanics and Modeling in Mechanobiology","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"Computational fluid dynamics, intracranial aneurysms, cerebral aneurysms, shape memory polymers, patient-specific therapeutics","lastPublishedDoi":"10.21203/rs.3.rs-5014601/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-5014601/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eThe hemodynamic and convective heat transfer effects of a patient-specific endovascular therapeutic agent based on shape memory polymer foam (SMPf) are evaluated using computational fluid dynamics studies for six patient-specific aneurysm geometries. The SMPf device is modeled as a continuous porous medium with full expansion for the flow studies and with various degrees of expansion for the heat transfer studies. The flow simulation parameters were qualitatively validated based on the existing literature. Further, a mesh independence study was conducted to verify an optimal cell size and reduce the computational costs. For convective heat transfer, a worst-case scenario is evaluated where the minimum volumetric flow rate is applied alongside the zero-flux boundary conditions. In the flow simulations, we found a reduction of the average intra-aneurysmal flow of \u0026gt;\u0026thinsp;85% and a reduction of the maximum intra-aneurysmal flow of \u0026gt;\u0026thinsp;45% for all presented geometries. These findings were compared with the literature on numerical simulations of hemodynamic and heat transfer of SMPf devices. The results obtained from this study can serve as a guide for optimizing the design and development of patient-specific SMPf devices aimed at personalized endovascular embolization of intracranial aneurysms.\u003c/p\u003e","manuscriptTitle":"CFD investigations of a shape-memory polymer foam-based endovascular embolization device for the treatment of intracranial aneurysms","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-10-18 09:37:01","doi":"10.21203/rs.3.rs-5014601/v1","editorialEvents":[{"type":"communityComments","content":0},{"type":"decision","content":"Revision requested","date":"2024-09-22T17:39:03+00:00","index":"","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-09-19T19:13:08+00:00","index":"hide","fulltext":""},{"type":"editorInvitedReview","content":"","date":"2024-09-17T10:50:38+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"198235223747969342541331818816033672247","date":"2024-09-09T12:25:06+00:00","index":"hide","fulltext":""},{"type":"reviewerAgreed","content":"281120690169005976481450325491888944370","date":"2024-09-08T10:17:46+00:00","index":"hide","fulltext":""},{"type":"reviewersInvited","content":"","date":"2024-09-07T12:02:46+00:00","index":"","fulltext":""},{"type":"editorAssigned","content":"","date":"2024-09-06T11:55:43+00:00","index":"","fulltext":""},{"type":"checksComplete","content":"","date":"2024-09-02T01:27:43+00:00","index":"","fulltext":""},{"type":"submitted","content":"Biomechanics and Modeling in Mechanobiology","date":"2024-09-02T00:33:57+00:00","index":"","fulltext":""}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"biomechanics-and-modeling-in-mechanobiology","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"bmmb","sideBox":"Learn more about [Biomechanics and Modeling in Mechanobiology](http://link.springer.com/journal/10237)","snPcode":"10237","submissionUrl":"https://submission.nature.com/new-submission/10237/3","title":"Biomechanics and Modeling in Mechanobiology","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"em","reportingPortfolio":"Springer Hybrid","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"f0fc16d7-2818-4d97-8b0c-5e479e57019d","owner":[],"postedDate":"October 18th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[],"tags":[],"updatedAt":"2024-12-02T17:21:53+00:00","versionOfRecord":{"articleIdentity":"rs-5014601","link":"https://doi.org/10.1007/s10237-024-01910-x","journal":{"identity":"biomechanics-and-modeling-in-mechanobiology","isVorOnly":false,"title":"Biomechanics and Modeling in Mechanobiology"},"publishedOn":"2024-11-25 15:57:31","publishedOnDateReadable":"November 25th, 2024"},"versionCreatedAt":"2024-10-18 09:37:01","video":"","vorDoi":"10.1007/s10237-024-01910-x","vorDoiUrl":"https://doi.org/10.1007/s10237-024-01910-x","workflowStages":[]},"version":"v1","identity":"rs-5014601","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-5014601","identity":"rs-5014601","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}