International Journal for Numerical Methods in Biomedical Engineering最新文献

{"title":"Influence of Fiber Dispersion Representation on the Accuracy of the Mechanical Response of Healthy and Aneurysmal Aortic Wall Tissue","authors":"Omid Ghorbani, Victorien Prot, Gerhard A. Holzapfel, Bjørn Skallerud","doi":"10.1002/cnm.70186","DOIUrl":"10.1002/cnm.70186","url":null,"abstract":"<div>\u0000 \u0000 <p>This study investigates the influence of fiber dispersion representation on the accuracy of the mechanical response of anisotropic hyperelastic material models for healthy and aneurysmal human aortic wall tissue under planar biaxial loading. Five fitting strategies were compared: (i) detailed integration of fiber orientation and dispersion data across the thickness; (ii) a single, thickness-averaged fiber distribution; (iii) two representative layers (media and adventitia); (iv) multiple layers discretized in the radial direction; and (v) a single layer with a symmetrical in-plane distribution of collagen fibers. The goodness of fit was calculated using the coefficient of determination averaged over all experimental tests. The most accurate fits were achieved with strategy (i), with a goodness of fit <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msubsup>\u0000 <mi>R</mi>\u0000 <mi>avg</mi>\u0000 <mn>2</mn>\u0000 </msubsup>\u0000 <mo>≈</mo>\u0000 <mn>0.96</mn>\u0000 <mo>±</mo>\u0000 <mn>0.02</mn>\u0000 </mrow>\u0000 <annotation>$$ {R}_{mathrm{avg}}^2approx 0.96pm 0.02 $$</annotation>\u0000 </semantics></math> for the combined cohort of aneurysm patients and the healthy cohort; comparable accuracy was achieved by the two-layer model (iii), <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msubsup>\u0000 <mi>R</mi>\u0000 <mi>avg</mi>\u0000 <mn>2</mn>\u0000 </msubsup>\u0000 <mo>≈</mo>\u0000 <mn>0.91</mn>\u0000 <mo>±</mo>\u0000 <mn>0.06</mn>\u0000 </mrow>\u0000 <annotation>$$ {R}_{mathrm{avg}}^2approx 0.91pm 0.06 $$</annotation>\u0000 </semantics></math>, but with significantly less complexity in model implementation. As a secondary objective of this study, regional parameters that correlate with the two-layer model (iii) and are intended for finite element analyses are presented. The method using a symmetrical single layer in model (v) yielded the lowest accuracy, <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msubsup>\u0000 <mi>R</mi>\u0000 <mi>avg</mi>\u0000 <mn>2</mn>\u0000 </msubsup>\u0000 <mo>≈</mo>\u0000 <mn>0.73</mn>\u0000 <mo>±</mo>\u0000 <mn>0.14</mn>\u0000 </mrow>\u0000 <annotation>$$ {R}_{mathrm{avg}}^2approx 0.73pm 0.14 $$</annotation>\u0000 </semantics></math>, and highlighted the need to capture non-symmetric fiber families. The modeling of elastin as a separate anisotropic material was also investigated by comparing two different constitutive modeling approaches. While the healthy cohort showed the best fit by includin","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"42 6","pages":""},"PeriodicalIF":2.4,"publicationDate":"2026-06-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148165337","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Investigation of Contrast Retention in PCOM Aneurysms as a Possible Marker for Instability and Rupture","authors":"Laurel M. M. Marsh,&nbsp;Fernando Mut,&nbsp;Alireza Chitsaz,&nbsp;Rainald Lohner,&nbsp;Anne M. Robertson,&nbsp;Naoki Kaneko,&nbsp;Juan R. Cebral","doi":"10.1002/cnm.70184","DOIUrl":"10.1002/cnm.70184","url":null,"abstract":"<p>Identification of cerebral aneurysms at risk of destabilization and rupture remains challenging. Although flow stagnation has been recognized as a risk factor, the underlying mechanisms and its predictive power are poorly understood. The purpose was to investigate possible associations between flow stagnation determined from computational fluid dynamics (CFD)-based virtual angiograms and aneurysm instability and rupture. A total of 312 aneurysms at the posterior communicating (PCOM) artery were analyzed with image-based CFD. Virtual angiograms were constructed by simulating contrast injections, and contrast retention metrics were calculated from time density curves of the aneurysm and their parent artery. These metrics were compared between 23 stable and 18 unstable aneurysms from a longitudinal dataset and between 129 unruptured and 142 ruptured aneurysms from a cross-sectional dataset. Aneurysmal contrast retention was significantly larger (<i>p</i> &lt; 0.05) in ruptured compared to unruptured aneurysms as well as when considering only small aneurysms (size &lt; 7 mm). Similar associations were obtained when comparing stable and growing aneurysms. Multivariate analysis showed that aneurysm instability can be forecasted from flow stagnation with good accuracy (AUC = 0.88), while aneurysm rupture with only moderate accuracy (AUC = 0.72) or marginal accuracy for small aneurysms (AUC = 0.66). Mean residence time was identified as a good predictor of aneurysm instability and rupture. Flow stagnation, determined by contrast retention in angiographies (virtual or real), is associated with PCOM aneurysm instability (growth or symptoms development) and can be used to identify aneurysms prone to further growth and potentially rupture. However, aneurysms can also rupture without exhibiting signs of flow stagnation.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"42 6","pages":""},"PeriodicalIF":2.4,"publicationDate":"2026-06-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/cnm.70184","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148165259","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A Novel Decomposition-Based Growth Model for Simulating Stress-Modulated Spinal Growth","authors":"Serhat Onur Çakmak,&nbsp;Ercan Gürses","doi":"10.1002/cnm.70187","DOIUrl":"10.1002/cnm.70187","url":null,"abstract":"<p>Understanding the mechanobiological mechanisms of spinal growth is essential for modeling progressive deformities, particularly adolescent idiopathic scoliosis (AIS). In this study, we present a novel finite element methodology for simulating mechanobiological growth governed by the Hueter-Volkmann law. The approach is developed in the large deformation framework and uses the Hueter-Volkmann law for the first time in a decomposition-based framework in which the growth is incorporated into the constitutive model via multiplicative decomposition of the deformation gradient. For comparison purposes, we also implement the sequential method widely used in the literature. Although both approaches are based on the Hueter-Volkmann law, which relates the normal compressive stress on the growth plate to local growth inhibition, they differ significantly in how this principle is integrated into the finite element framework. Unlike the sequential method, which requires successive growth computation and mechanical solution steps, along with repeated updates of the growth direction, the proposed decomposition-based approach solves stress-induced deformations and growth simultaneously. It is formulated in the Lagrangian configuration, eliminating the need to update the growth direction throughout the simulation. The proposed formulation not only offers superior numerical stability but also exhibits lower computational complexity compared to the sequential method, particularly under large deformations. The method is evaluated using a finite element model of a simplified Functional Spinal Unit (FSU) under symmetric and asymmetric compressive loading over a two-year growth period. While both approaches produce comparable wedge angle progression, the decomposition-based formulation, which more realistically resembles the actual growth process, results in smoother deformations, less element distortion, and substantially reduced computational cost. Additionally, a robust wedge angle calculation technique is introduced using least plane fitting to full endplate nodes of the vertebral bodies, improving geometric accuracy and reducing mesh sensitivity. Together, these features establish a reliable and efficient framework for long-term spinal growth simulations and offer a promising foundation for future patient-specific modeling and clinical applications.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"42 6","pages":""},"PeriodicalIF":2.4,"publicationDate":"2026-06-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/cnm.70187","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148165308","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Automated Centerline Extraction From Meshed Vascular Models","authors":"Gala Sanchez Van Moer,&nbsp;Shawn C. Shadden","doi":"10.1002/cnm.70183","DOIUrl":"10.1002/cnm.70183","url":null,"abstract":"<p>Centerlines of vascular structures are essential for the analysis of vascular anatomy as well as the development of computational models for simulating biomechanics. In this paper, we present a centerline extraction framework for vascular geometries implemented on piece-wise linear surface meshes, from which tetrahedral volumetric meshes are generated. Our approach involves solving the Eikonal equation using the finite element method to describe wave propagation for automatic endpoint identification and centerline tracing. We evaluated the framework on 19 vascular meshes with varying anatomies and construction methods, comparing our results to those obtained using the widely used Vascular Modeling Toolkit (VMTK) as a benchmark. The proposed method produces well-centered centerlines within a reasonable time frame without the need for tedious outlet labeling required by VMTK. Furthermore, it considerably outperforms VMTK when extracting centerlines from less smooth surface meshes, such as those generated using machine learning-based segmentation techniques.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"42 6","pages":""},"PeriodicalIF":2.4,"publicationDate":"2026-06-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/cnm.70183","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148145376","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Leaflets Morphology Optimization of Valved Pulmonary Arterial Conduit Based on Experimental Design","authors":"Xiaofan Zheng,&nbsp;Yong Feng,&nbsp;Yunhan Cai,&nbsp;Wentao Yan,&nbsp;Qiqi Shi,&nbsp;Huifeng Zhang,&nbsp;Shengzhang Wang","doi":"10.1002/cnm.70185","DOIUrl":"https://doi.org/10.1002/cnm.70185","url":null,"abstract":"<div>\u0000 \u0000 <p>This study investigated the influence of leaflet geometry on the hemodynamic performance of pulmonary valved conduits and optimized leaflet trimming parameters using a design of experiments (DOE) approach combined with in vitro hydrodynamic testing. In vitro hydrodynamic experiments were performed to validate the feasibility of the valved conduit. A computational model was then established based on the experimental setup. DOE was applied to systematically assess and optimize leaflet tailoring parameters. Fluid–structure interaction simulations showed that leaflet tailoring significantly affects hemodynamic outcomes. Reducing leaflet height improved conduit performance. At a cardiac output of 3.5 L/min, the optimized valve exhibited a 3.79% decrease in peak velocity, a 45.18% reduction in maximum equivalent strain, and a 66.57% decrease in coaptation area. These findings suggest that scaling down the free edge length relative to graft diameter and moderately reducing leaflet height may help reduce diastolic regurgitation caused by redundant leaflet overlap while lowering leaflet strain, thereby enhancing valve durability.</p>\u0000 </div>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"42 6","pages":""},"PeriodicalIF":2.4,"publicationDate":"2026-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148099670","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Evaluation of Hemodynamic Environment of Intracranial Aneurysms After Flow Disruption Based on Angiographic Signatures","authors":"Fernando Mut,&nbsp;Laurel Marsh,&nbsp;David Kallmes,&nbsp;Ramanathan Kadirvel,&nbsp;Juan R. Cebral","doi":"10.1002/cnm.70179","DOIUrl":"https://doi.org/10.1002/cnm.70179","url":null,"abstract":"<p>Incomplete aneurysm occlusion after flow disruption treatment strongly depends on the persistence of strong flow into the aneurysm, which in turn depends on the aneurysm and parent artery geometry as well as the device placement. However, at the time of treatment it remains challenging to predict future outcomes (complete or incomplete occlusion). This paper describes a strategy to evaluate the local hemodynamic environment created immediately after flow-modifying device deployment to treat cerebral aneurysms and connect it to the subsequent long-term outcomes. The approach is based on the concept of angiographic signatures that can be extracted from DSA sequences or computational models and represented as images that can be used to compare different aneurysmal environments. The methodology is demonstrated on two small series of four experimental rabbit aneurysms each, treated with either intrasaccular flow disruptors or endoluminal flow diverters. Noticeable differences in the angiographic signatures of completely and incompletely occluded aneurysms at follow-up were observed. Specifically, larger values of mean transit times and longer time-to-peaks were observed deeper into the aneurysms that remained incompletely occluded in both series. The findings suggest that angiographic signatures can be used to assess the underlying hemodynamic environment immediately after device implantation and prognosticate the likelihood of future complete or incomplete occlusion.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"42 6","pages":""},"PeriodicalIF":2.4,"publicationDate":"2026-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/cnm.70179","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148099669","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Numerical Inverse Design of Patient-Specific Dental Implants: Accelerating FEA-Based Optimization via Evolutionary Neural Surrogates","authors":"María Prados-Privado","doi":"10.1002/cnm.70182","DOIUrl":"10.1002/cnm.70182","url":null,"abstract":"<div>\u0000 \u0000 <p>Mechanical complications in dental implantology often arise from a mismatch between standardized geometries and patient-specific anatomical constraints. While high-fidelity finite element analysis (FEA) is the gold standard for identifying these risks, its high computational cost creates a “bottleneck” that prevents its use in routine surgical planning. To develop a numerical inverse design framework that overcomes this limitation by accelerating the generation of optimal implant geometries through a neural-surrogate-assisted optimization engine. A high-fidelity dataset of 3000 high-fidelity 3D FEA simulations was used to train a multilayer perceptron (MLP) regressor (<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msup>\u0000 <mi>R</mi>\u0000 <mn>2</mn>\u0000 </msup>\u0000 <mo>=</mo>\u0000 <mn>0.9673</mn>\u0000 </mrow>\u0000 <annotation>$$ {R}^2=0.9673 $$</annotation>\u0000 </semantics></math>). This model served as a real-time inference engine replacing expensive iterative simulations within surrogate-assisted evolutionary optimization framework. The framework's accuracy and clinical validity were tested on a synthetic cohort of 50 virtual patients representing diverse bone qualities and loading conditions. The accelerated inverse design system achieved a statistically significant and substantial reduction in peak von Mises stress at the bone–implant interface compared with standard clinical protocols (<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>p</mi>\u0000 <mo>&lt;</mo>\u0000 <mn>0.001</mn>\u0000 </mrow>\u0000 <annotation>$$ p&lt;0.001 $$</annotation>\u0000 </semantics></math>, Cohen's <i>d</i> = 3.22). The optimization successfully identified patient-specific compensatory strategies, such as prescribing wider diameters for low-density bone to maximize load distribution. By integrating deep learning surrogates to accelerate complex numerical optimization, this framework enables real-time, patient-specific implant prescription. This methodology effectively bypasses the computational burden of traditional FEA, offering a scalable numerical solution for personalized surgical planning.</p>\u0000 </div>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"42 6","pages":""},"PeriodicalIF":2.4,"publicationDate":"2026-05-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148058158","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Influence of Input Data Composition and Measurement Errors on Computational Model-Mediated Assessment of Cardiovascular Properties: An In Silico Study","authors":"Zhicheng Zhang,&nbsp;Zhaojun Li,&nbsp;Di Sun,&nbsp;Fuyou Liang","doi":"10.1002/cnm.70181","DOIUrl":"10.1002/cnm.70181","url":null,"abstract":"<div>\u0000 \u0000 <p>The integration of computational cardiovascular models with clinical data through parameter estimation algorithms represents a promising approach for quantitative cardiovascular assessment. However, the sensitivity of assessment to variations in the composition of and measurement errors in input clinical data remains insufficiently understood. In this study, we utilized a cardiovascular model to generate in silico datasets representing diverse virtual subjects. From these datasets, hemodynamic data mimicking noninvasive clinical measurements were extracted as inputs for the parameter estimation algorithm (herein the Levenberg–Marquardt algorithm), while the corresponding values of model parameters representing key cardiovascular properties (e.g., arterial stiffness, cardiac function, and vascular resistance) served as the ground truth for evaluating the accuracy of parameter estimation. The obtained results showed that input datasets comprising brachial arterial blood pressures and flow velocities in the ascending aorta and four peripheral arteries (supplying major organs/tissues) enabled accurate parameter estimation. Cardiovascular models reassigned with the estimated parameter values precisely replicated the input data and predicted other important hemodynamic variables (e.g., aortic blood pressure). While these capabilities were largely maintained upon the omission of flow velocity data in a single peripheral artery, they were significantly compromised when flow velocity data in two peripheral arteries were removed. Furthermore, the introduction of measurement errors into the input data substantially increased errors in both parameter estimation and hemodynamic prediction. In summary, this study underscores the importance of rigorous clinical data selection and high-precision in vivo measurements for enhancing the reliability of computational model-mediated cardiovascular assessment.</p>\u0000 </div>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"42 6","pages":""},"PeriodicalIF":2.4,"publicationDate":"2026-05-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148058121","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"","authors":"","doi":"","DOIUrl":"","url":null,"abstract":"","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"42 5","pages":""},"PeriodicalIF":2.4,"publicationDate":"2026-05-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148089711","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"","authors":"","doi":"","DOIUrl":"","url":null,"abstract":"","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"42 5","pages":""},"PeriodicalIF":2.4,"publicationDate":"2026-05-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148083735","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}