Biomechanics and Modeling in Mechanobiology最新文献

{"title":"A Hamilton principle-based model for diffusion-driven biofilm growth","authors":"Felix Klempt,&nbsp;Meisam Soleimani,&nbsp;Peter Wriggers,&nbsp;Philipp Junker","doi":"10.1007/s10237-024-01883-x","DOIUrl":"10.1007/s10237-024-01883-x","url":null,"abstract":"<div><p>Dense communities of bacteria, also known as biofilms, are ubiquitous in all of our everyday life. They are not only always surrounding us, but are also active inside our bodies, for example in the oral cavity. While some biofilms are beneficial or even necessary for human life, others can be harmful. Therefore, it is highly important to gain an in-depth understanding of biofilms which can be achieved by <i>in vitro</i> or <i>in vivo</i> experiments. Since these experiments are often time-consuming or expensive, <i>in silico</i> models have proven themselves to be a viable tool in assisting the description and analysis of these complicated processes. Current biofilm growth simulations are using mainly two approaches for describing the underlying models. The volumetric approach splits the deformation tensor into a growth and an elastic part. In this approach, the mass never changes, unless some additional constraints are enforced. The density-based approach, on the other hand, uses an evolution equation to update the growing tissue by adding mass. Here, the density stays constant, and no pressure is exerted. The <i>in silico</i> model presented in this work combines the two approaches. Thus, it is possible to capture stresses inside of the biofilm while adding mass. Since this approach is directly derived from Hamilton’s principle, it fulfills the first and second law of thermodynamics automatically, which other models need to be checked for separately. In this work, we show the derivation of the model as well as some selected numerical experiments. The numerical experiments show a good phenomenological agreement with what is to be expected from a growing biofilm. The numerical behavior is stable, and we are thus capable of solving complicated boundary value problems. In addition, the model is very reactive to different input parameters, thereby different behavior of various biofilms can be captured without modifying the model.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"2091 - 2113"},"PeriodicalIF":3.0,"publicationDate":"2024-09-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10237-024-01883-x.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142338777","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Multiscale computational analysis of the steady fluid flow through a lymph node","authors":"Alberto Girelli,&nbsp;Giulia Giantesio,&nbsp;Alessandro Musesti,&nbsp;Raimondo Penta","doi":"10.1007/s10237-024-01879-7","DOIUrl":"10.1007/s10237-024-01879-7","url":null,"abstract":"<div><p>Lymph Nodes (LNs) are crucial to the immune and lymphatic systems, filtering harmful substances and regulating lymph transport. LNs consist of a lymphoid compartment (LC) that forms a porous bulk region, and a subcapsular sinus (SCS), which is a free-fluid region. Mathematical and mechanical challenges arise in understanding lymph flow dynamics. The highly vascularized lymph node connects the lymphatic and blood systems, emphasizing its essential role in maintaining the fluid balance in the body. In this work, we describe a mathematical model in a steady setting to describe the lymph transport in a lymph node. We couple the fluid flow in the SCS governed by an incompressible Stokes equation with the fluid flow in LC, described by a model obtained by means of asymptotic homogenisation technique, taking into account the multiscale nature of the node and the fluid exchange with the blood vessels inside it. We solve this model using numerical simulations and we analyze the lymph transport inside the node to elucidate its regulatory mechanisms and significance. Our results highlight the crucial role of the microstructure of the lymph node in regularising its fluid balance. These results can pave the way to a better understanding of the mechanisms underlying the lymph node’s multiscale functionalities which can be significantly affected by specific physiological and pathological conditions, such as those characterising malignant tissues.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"2005 - 2023"},"PeriodicalIF":3.0,"publicationDate":"2024-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10237-024-01879-7.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142338778","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"3D growth and remodeling theory supports the hypothesis of staphyloma formation from local scleral weakening under normal intraocular pressure","authors":"Fabian A. Braeu,&nbsp;Stéphane Avril,&nbsp;Michaël J. A. Girard","doi":"10.1007/s10237-024-01885-9","DOIUrl":"10.1007/s10237-024-01885-9","url":null,"abstract":"<div><p>The purpose of this study was to assess whether growth and remodeling (G&amp;R) theory could explain staphyloma formation from a local scleral weakening—as could occur from age-related elastin degradation, myopia progression, or other factors. A finite element model of a healthy eye was reconstructed, including the lamina cribrosa, the peripapillary sclera, and the peripheral sclera. The homogenized constrained mixture model was employed to simulate the adaptation of the sclera to alterations in its biomechanical environment over a duration of 13.7 years. G&amp;R processes were triggered by reducing the shear stiffness of the ground matrix in the peripapillary sclera and lamina cribrosa by 85%. Three distinct G&amp;R scenarios were investigated: (1) low mass turnover rate in combination with transmural volumetric growth; (2) high mass turnover rate in combination with transmural volumetric growth; and (3) high mass turnover rate in combination with mass density growth. In <b>scenario 1</b>, we observed a significant outpouching of the posterior pole, closely resembling the shape of a Type-III staphyloma. Additionally, we found a notable change in scleral curvature and a thinning of the peripapillary sclera by 84%. In contrast, <b>scenario 2</b> and <b>3</b> exhibited less drastic deformations, with stable posterior staphylomas after approximately 7 years. Our proposed framework suggests that local scleral weakening is sufficient to trigger staphyloma formation under a normal level of intraocular pressure. Our model also reproduced characteristics of Type-III staphylomas. With patient-specific scleral geometries (as could be obtained with wide-field optical coherence tomography), our framework could be clinically translated to help us identify those at risks of developing posterior staphylomas.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"2137 - 2154"},"PeriodicalIF":3.0,"publicationDate":"2024-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142338776","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Prediction of shear stress imposed on alveolar epithelium of healthy and diseased lungs","authors":"Alexandros Livanos,&nbsp;Konstantinos Bouchoris,&nbsp;Kyriaki-Evangelia Aslani,&nbsp;Konstantinos Gourgoulianis,&nbsp;Vasilis Bontozoglou","doi":"10.1007/s10237-024-01889-5","DOIUrl":"10.1007/s10237-024-01889-5","url":null,"abstract":"<div><p>Lung alveoli are modeled as spherical caps, lined internally by a thin surfactant-laden liquid film, and the periodic wall shear stress exerted along the epithelium during small-amplitude radial oscillations of their wall is computed. A novel set of boundary conditions, applied at the rim, reveals the dominant role of Marangoni stresses. These stresses develop along the air/liquid interface due to spatial gradients of interfacial surfactant concentration and are transported to the wall by the action of viscosity. The effect of a variety of geometric and functional characteristics, including rim interstitial thickness, alveolar opening angle and liquid film thickness and viscosity, is interrogated, and the results are discussed in relation to the onset and evolution of acute and chronic lung diseases, such as asthmatic attacks, pulmonary emphysema and pulmonary fibrosis.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"2213 - 2227"},"PeriodicalIF":3.0,"publicationDate":"2024-09-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142278428","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical analysis of blood flow in a branched modular stent-graft for aneurysms covering all zones of the aortic arch","authors":"Mário Luis Ferreira da Silva,&nbsp;Matheus Carvalho Barbosa Costa,&nbsp;Saulo de Freitas Gonçalves,&nbsp;Rudolf Huebner,&nbsp;Túlio Pinho Navarro","doi":"10.1007/s10237-024-01887-7","DOIUrl":"10.1007/s10237-024-01887-7","url":null,"abstract":"<div><p>Due to the anatomical complexity of the aortic arch for the development of stent-grafts for total repair, this region remains without a validated and routinely used endovascular option. In this work, a modular stent-graft for aneurysms that covers all aortic arch zones, proposed by us and previously structurally evaluated, was evaluated from the point of view of haemodynamics using fluid-structural numerical simulations. Blood was assumed to be non-Newtonian shear-thinning using the Carreau model, and the arterial wall was assumed to be anisotropic hyperelastic using the Holzapfel model. Nitinol and expanded polytetrafluoroethylene (PTFE-e) were used as materials for the stents and the graft, respectively. Nitinol was modelled as a superelastic material with shape memory by the Auricchio model, and PTFE-e was modelled as an isotropic linear elastic material. To validate the numerical model, a silicone model representative of the aneurysmal aorta was subjected to tests on an experimental bench representative of the circulatory system. The numerical results showed that the stent-graft restored flow behaviour, making it less oscillatory, but increasing the strain rate, turbulence kinetic energy, and viscosity compared to the pathological case. Taking the mean of the entire cycle, the increase in turbulence kinetic energy was 198.82% in the brachiocephalic trunk, 144.63% in the left common carotid artery and 284.03% in the left subclavian artery after stent-graft implantation. Based on wall shear stress parameters, it was possible to identify that the internal branches of the stent-graft and the stent-graft fixation sites in the artery were the most favourable regions for the deposition and accumulation of thrombus. In these regions, the oscillating shear index reached the maximum value of 0.5 and the time-averaged wall shear stress was close to zero, which led the relative residence time to reach values above 15 Pa⁻¹. The stent-graft was able to preserve flow in the supra-aortic branches.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"2177 - 2191"},"PeriodicalIF":3.0,"publicationDate":"2024-09-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142278427","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Simulating the impact of tumor mechanical forces on glymphatic networks in the brain parenchyma","authors":"Saeed Siri,&nbsp;Alice Burchett,&nbsp;Meenal Datta","doi":"10.1007/s10237-024-01890-y","DOIUrl":"10.1007/s10237-024-01890-y","url":null,"abstract":"&lt;div&gt;&lt;p&gt;The brain glymphatic system is currently being explored in the context of many neurological disorders and diseases, including traumatic brain injury, Alzheimer’s disease, and ischemic stroke. However, little is known about the impact of brain tumors on glymphatic function. Mechanical forces generated during tumor development and growth may be responsible for compromised glymphatic transport pathways, reducing waste clearance and cerebrospinal fluid (CSF) transport in the brain parenchyma. One such force is solid stress, i.e., growth-induced forces from cell hyperproliferation and excess matrix deposition. Because there are no prior studies assessing the impact of tumor-derived solid stress on glymphatic system structure and performance in the brain parenchyma, this study serves to fill an important gap in the field. We adapted a previously developed Electrical Analog Model using MATLAB Simulink for glymphatic transport coupled with Finite Element Analysis for tumor mechanical stresses and strains in COMSOL. This allowed simulation of the impact of tumor mechanical force generation on fluid transport within brain parenchymal glymphatic units—which include perivascular spaces, astrocytic networks, interstitial spaces, and capillary basement membranes. We conducted a parametric analysis to compare the contributions of tumor size, tumor proximity, and ratio of glymphatic subunits to the stress and strain experienced by the glymphatic unit and corresponding reduction in flow rate of CSF. Mechanical stresses intensify with proximity to the tumor and increasing tumor size, highlighting the vulnerability of nearby glymphatic units to tumor-derived forces. Our stress and strain profiles reveal compressive deformation of these surrounding glymphatics and demonstrate that varying the relative contributions of astrocytes vs. interstitial spaces impact the resulting glymphatic structure significantly under tumor mechanical forces. Increased tumor size and proximity caused increased stress and strain across all glymphatic subunits, as does decreased astrocyte composition. Indeed, our model reveals an inverse correlation between extent of astrocyte contribution to the composition of the glymphatic unit and the resulting mechanical stress. This increased mechanical strain across the glymphatic unit decreases the venous efflux rate of CSF, dependent on the degree of strain and the specific glymphatic subunit of interest. For example, a 20% mechanical strain on capillary basement membranes does not significantly decrease venous efflux (2% decrease in flow rates), while the same magnitude of strain on astrocyte networks and interstitial spaces decreases efflux flow rates by 7% and 22%, respectively. Our simulations reveal that solid stress from growing brain tumors directly reduces glymphatic fluid transport, independently from biochemical effects from cancer cells. Understanding these pathophysiological implications is crucial for developing targeted inter","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"2229 - 2241"},"PeriodicalIF":3.0,"publicationDate":"2024-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10237-024-01890-y.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142249689","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Computational analysis of heart valve growth and remodeling after the Ross procedure","authors":"Elmer Middendorp,&nbsp;Fabian Braeu,&nbsp;Frank P. T. Baaijens,&nbsp;Jay D. Humphrey,&nbsp;Christian J. Cyron,&nbsp;Sandra Loerakker","doi":"10.1007/s10237-024-01874-y","DOIUrl":"10.1007/s10237-024-01874-y","url":null,"abstract":"<div><p>During the Ross procedure, an aortic heart valve is replaced by a patient’s own pulmonary valve. The pulmonary autograft subsequently undergoes substantial growth and remodeling (G&amp;R) due to its exposure to increased hemodynamic loads. In this study, we developed a homogenized constrained mixture model to understand the observed adaptation of the autograft leaflets in response to the changed hemodynamic environment. This model was based on the hypothesis that tissue G&amp;R aims to preserve mechanical homeostasis for each tissue constituent. To model the Ross procedure, we simulated the exposure of a pulmonary valve to aortic pressure conditions and the subsequent G&amp;R of the valve. Specifically, we investigated the effects of assuming either stress- or stretch-based mechanical homeostasis, the use of blood pressure control, and the effect of root dilation. With this model, we could explain different observations from published clinical studies, such as the increase in thickness, change in collagen organization, and change in tissue composition. In addition, we found that G&amp;R based on stress-based homeostasis could better capture the observed changes in tissue composition than G&amp;R based on stretch-based homeostasis, and that root dilation or blood pressure control can result in more leaflet elongation. Finally, our model demonstrated that successful adaptation can only occur when the mechanically induced tissue deposition is sufficiently larger than tissue degradation, such that leaflet thickening overrules leaflet dilation. In conclusion, our findings demonstrated that G&amp;R based on mechanical homeostasis can capture the observed heart valve adaptation after the Ross procedure. Finally, this study presents a novel homogenized mixture model that can be used to investigate other cases of heart valve G&amp;R as well.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"1889 - 1907"},"PeriodicalIF":3.0,"publicationDate":"2024-09-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10237-024-01874-y.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142225457","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A morpho-viscoelasticity theory for growth in proliferating aggregates","authors":"Prakhar Bandil,&nbsp;Franck J. Vernerey","doi":"10.1007/s10237-024-01886-8","DOIUrl":"10.1007/s10237-024-01886-8","url":null,"abstract":"<div><p>Despite significant research efforts in the continuum modeling of biological growth, certain aspects have been overlooked. For instance, numerous investigations have examined the influence of morphogenetic cell behaviors, like division and intercalation, on the mechanical response of passive (non-growing) tissues. Yet, their impact on active growth dynamics remains inadequately explored. A key reason for this inadequacy stems from challenges in the continuum treatment of cell-level processes. While some coarse-grained models have been proposed to address these shortcomings, a focus on cell division and cell expansion has been missing, rendering them unusable when it comes to modeling growth. Moreover, existing studies are limited to two-dimensional tissues and are yet to be formally extended to three-dimensional multicellular systems. To address these limitations, we here present a generalized multiscale model for three-dimensional aggregates that accounts for complex morphogenetic movements that include division, expansion, and intercalation. The proposed continuum theory thus allows for a comprehensive exploration into the growth and dissipation mechanics of proliferating aggregates, such as spheroids and organoids. \u0000</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"2155 - 2176"},"PeriodicalIF":3.0,"publicationDate":"2024-09-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142103128","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Personalized computational electro-mechanics simulations to optimize cardiac resynchronization therapy","authors":"Emilia Capuano,&nbsp;Francesco Regazzoni,&nbsp;Massimiliano Maines,&nbsp;Silvia Fornara,&nbsp;Vanessa Locatelli,&nbsp;Domenico Catanzariti,&nbsp;Simone Stella,&nbsp;Fabio Nobile,&nbsp;Maurizio Del Greco,&nbsp;Christian Vergara","doi":"10.1007/s10237-024-01878-8","DOIUrl":"10.1007/s10237-024-01878-8","url":null,"abstract":"<div><p>In this study, we present a computational framework designed to evaluate virtual scenarios of cardiac resynchronization therapy (CRT) and compare their effectiveness based on relevant clinical biomarkers. Our approach involves electro-mechanical numerical simulations personalized, for patients with left bundle branch block, by means of a calibration obtained using data from Electro-Anatomical Mapping System (EAMS) measures acquired by cardiologists during the CRT procedure, as well as ventricular pressures and volumes, both obtained pre-implantation. We validate the calibration by using EAMS data coming from right pacing conditions. Three patients with fibrosis and three without are considered to explore various conditions. Our virtual scenarios consist of personalized numerical experiments, incorporating different positions of the left electrode along reconstructed epicardial veins; different locations of the right electrode; different ventriculo-ventricular delays. The aim is to offer a comprehensive tool capable of optimizing CRT efficiency for individual patients. We provide preliminary answers on optimal electrode placement and delay, by computing some relevant biomarkers such as <span>(dP/dt_{max})</span>, ejection fraction, stroke work. From our numerical experiments, we found that the latest activated segment during sinus rhythm is an effective choice for the non-fibrotic cases for the location of the left electrode. Also, our results showed that the activation of the right electrode before the left one seems to improve the CRT performance for the non-fibrotic cases. Last, we found that the CRT performance seems to improve by positioning the right electrode halfway between the base and the apex. This work is on the line of computational works for the study of CRT and introduces new features in the field, such as the presence of the epicardial veins and the movement of the right electrode. All these studies from the different research groups can in future synergistically flow together in the development of a tool which clinicians could use during the procedure to have quantitative information about the patient’s propagation in different scenarios.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"1977 - 2004"},"PeriodicalIF":3.0,"publicationDate":"2024-08-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10237-024-01878-8.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142078723","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A comprehensive numerical approach to coil placement in cerebral aneurysms: mathematical modeling and in silico occlusion classification","authors":"Fabian Holzberger,&nbsp;Markus Muhr,&nbsp;Barbara Wohlmuth","doi":"10.1007/s10237-024-01882-y","DOIUrl":"10.1007/s10237-024-01882-y","url":null,"abstract":"<div><p>Endovascular coil embolization is one of the primary treatment techniques for cerebral aneurysms. Although it is a well-established and minimally invasive method, it bears the risk of suboptimal coil placement which can lead to incomplete occlusion of the aneurysm possibly causing recurrence. One of the key features of coils is that they have an imprinted natural shape supporting the fixation within the aneurysm. For the spatial discretization, our mathematical coil model is based on the discrete elastic rod model which results in a dimension-reduced 1D system of differential equations. We include bending and twisting responses to account for the coils natural curvature and allow for the placement of several coils having different material parameters. Collisions between coil segments and the aneurysm wall are handled by an efficient contact algorithm that relies on an octree based collision detection. In time, we use a standard symplectic semi-implicit Euler time stepping method. Our model can be easily incorporated into blood flow simulations of embolized aneurysms. In order to differentiate optimal from suboptimal placements, we employ a suitable in silico Raymond–Roy-type occlusion classification and measure the local packing density in the aneurysm at its neck, wall region and core. We investigate the impact of uncertainties in the coil parameters and embolization procedure. To this end, we vary the position and the angle of insertion of the micro-catheter, and approximate the local packing density distributions by evaluating sample statistics.</p></div>","PeriodicalId":489,"journal":{"name":"Biomechanics and Modeling in Mechanobiology","volume":"23 6","pages":"2063 - 2089"},"PeriodicalIF":3.0,"publicationDate":"2024-08-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10237-024-01882-y.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142003295","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}