Christopher Blum, Ulrich Steinseifer, Michael Neidlin
{"title":"Toward Uncertainty-Aware Hemolysis Modeling: A Universal Approach to Address Experimental Variance","authors":"Christopher Blum, Ulrich Steinseifer, Michael Neidlin","doi":"10.1002/cnm.70040","DOIUrl":"https://doi.org/10.1002/cnm.70040","url":null,"abstract":"<p>The purpose of this study is to address the lack of uncertainty quantification in numerical hemolysis models, which are critical for medical device evaluations. Specifically, we aim to develop a probabilistic hemolysis model, which incorporates experimental variability using the Markov Chain Monte Carlo (MCMC) method to enhance predictive accuracy and robustness. Initially, we examined the objective function landscape for fitting a Power Law hemolysis model, whose parameters are derived from inherently uncertain experimental data, by employing a grid search approach. Building on this, we applied MCMC to derive detailed stochastic distributions for the hemolysis Power Law model parameters <i>C</i>, <i>α</i>, and <i>β</i>. These distributions were then propagated through a reduced order model of the FDA benchmark pump to quantify the experimental uncertainty in hemolysis measurements with respect to the predicted pump hemolysis. Our analysis revealed a global flat minimum in the objective function landscape of the multi-parameter power law model, a phenomenon attributable to fundamental mathematical limitations in the fitting process. The probabilistic hemolysis model converged to a constant optimal <i>C</i> = 3.515 × 10<sup>−5</sup> and log normal distributions of <i>α</i> and <i>β</i> with means of 0.614 and 1.795, respectively. This probabilistic approach successfully captured both the mean and variance observed in the experimental FDA benchmark pump data. In comparison, conventional deterministic models are not able to describe experimental variation. Incorporating uncertainty quantification through MCMC enhances the robustness and predictive accuracy of hemolysis models. This method allows for better comparison of simulated hemolysis outcomes with in vitro experiments and can integrate additional datasets, potentially setting a new standard in hemolysis modeling.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"41 5","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-05-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/cnm.70040","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143944928","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}
Lei Chen, Guoqiu Liu, Lei Xian, Bo Zhang, Wen-Quan Tao
{"title":"Numerical Simulation of Blood Flow in Nutcracker Syndrome: Acquisition of Hemodynamic Parameters and Clinical Application","authors":"Lei Chen, Guoqiu Liu, Lei Xian, Bo Zhang, Wen-Quan Tao","doi":"10.1002/cnm.70031","DOIUrl":"https://doi.org/10.1002/cnm.70031","url":null,"abstract":"<div>\u0000 \u0000 <p>Adopting noninvasive techniques to provide more precise parameters related to the clinical diagnosis and treatment of nutcracker syndrome is currently the focus of research on nutcracker syndrome. However, due to individual differences in patients and limitations in monitoring techniques, there is an urgent need for a new method to obtain more accurate parameters. This work is based on imaging data from two patients with nutcracker syndrome and numerically simulates the blood flow process in the left renal vein, revealing different clinical symptoms caused by changes in flow velocity, pressure, and wall shear stress. Besides, this work also compares the dynamic changes of blood flow parameters under two different degrees of compression. The study indicates that an increase in pressure at the entrance of the left renal vein increases the risk of venous congestion. The flow separation reduces the blood flow rate by 50%–60%, causing a series of flow obstacles. The wall shear stress near the compressed area increased by 15–20 times, exacerbating the damage of blood flow to the left renal vein. The increase in the degree of compression exacerbates flow barriers and the impact of blood flow on the vascular wall. This study introduces a method of obtaining hemodynamic parameters through computational fluid dynamics and summarizes the clinical symptoms caused by abnormal changes in different blood flow parameters. This method provides a more reliable approach for the clinical diagnosis of nutcracker syndrome and the optimization design of extracorporeal stent structures since it is not limited by monitoring techniques.</p>\u0000 </div>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"41 5","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-05-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143914541","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}
Yingjie Xia, Changpeng Wang, Yan Wang, Fuyou Liang
{"title":"A Computational Model-Based Study on Trans-Stenotic Pressure Ratio of Carotid Artery Stenosis and Its Predictive Value for Cerebral Ischemia","authors":"Yingjie Xia, Changpeng Wang, Yan Wang, Fuyou Liang","doi":"10.1002/cnm.70044","DOIUrl":"https://doi.org/10.1002/cnm.70044","url":null,"abstract":"<div>\u0000 \u0000 <p>Trans-stenotic pressure ratio (herein denoted by dpPR) has been proposed as a complementary index to stenosis rate (SR) for assessing the functional severity of carotid artery stenosis (CAS); however, it remains unclear how well dpPR can indicate cerebral ischemia. In this study, a physiology-based computational model of the cerebral circulation was developed to yield a tool for generating large-scale in silico data to characterize the changes of the dpPR of the left internal CAS in response to variations in SR and various anatomical/pathophysiological factors that represent inter-patient differences. In addition, a cerebral ischemia index (CII) was defined to evaluate the predictive value of dpPR for cerebral ischemia. Results showed that dpPR was affected by many factors unrelated to the severity of stenosis, such as the anatomical structure and geometrical size of cerebral arteries, mean systemic arterial blood pressure (MAP), flow autoregulation function of cerebral microcirculation (quantified by CFAI), and coexisting contralateral CAS. In comparison with SR, dpPR exhibited a stronger correlation with CII. In particular, the relationship between dpPR and CII was found to be describable by a mathematical function if MAP and CFAI were fixed. The findings not only deepen our understanding of the physiological implications of dpPR but also provide valuable theoretical insights to guide the application of dpPR in clinical practice.</p>\u0000 </div>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"41 5","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-05-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143914252","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":"Criterion for Assessing Accumulated Neurotoxicity of Alpha-Synuclein Oligomers in Parkinson's Disease","authors":"Andrey V. Kuznetsov","doi":"10.1002/cnm.70027","DOIUrl":"https://doi.org/10.1002/cnm.70027","url":null,"abstract":"<p>The paper introduces a parameter called “accumulated neurotoxicity” of α-syn oligomers, which measures the cumulative damage these toxic species inflict on neurons over time, given the years it typically takes for such damage to manifest. A threshold value for accumulated neurotoxicity is estimated, beyond which neuron death is likely. Numerical results suggest that rapid deposition of α-syn oligomers into fibrils minimizes neurotoxicity, indicating that the formation of Lewy bodies might play a neuroprotective role. Strategies such as reducing α-syn monomer production or enhancing degradation can decrease accumulated neurotoxicity. In contrast, slower degradation (reflected by longer half-lives of monomers and free aggregates) increases neurotoxicity, supporting the idea that impaired protein degradation may contribute to Parkinson's disease progression. Accumulated neurotoxicity is highly sensitive to the half-deposition time of free α-syn aggregates into fibrils, exhibiting a sharp increase as it transitions from negligible to elevated levels, indicative of neural damage.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"41 4","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/cnm.70027","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143884208","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}
Åshild Telle, Verena Charwat, Bérénice Charrez, Henrik Finsberg, Kevin E. Healy, Samuel T. Wall
{"title":"Estimation of Active Tension in Cardiac Microtissues by Solving a PDE-Constrained Optimization Problem","authors":"Åshild Telle, Verena Charwat, Bérénice Charrez, Henrik Finsberg, Kevin E. Healy, Samuel T. Wall","doi":"10.1002/cnm.70034","DOIUrl":"https://doi.org/10.1002/cnm.70034","url":null,"abstract":"<p>Microphysiological systems (MPS) provide a highly controlled environment for the development and testing of human-induced pluripotent stem cell-based cardiac microtissues, with promising applications in disease modeling and drug development. Through optical measurements in such systems, we can quantify mechanical features such as motion and velocity during contraction. While these are useful for evaluating relative changes in muscle twitch, it remains challenging to quantify and characterize the actual active tension driving the contraction. Here, we aimed to quantify the active tension over time and space by solving an inverse problem in cardiac mechanics expressed by partial differential equations (PDEs). We formulated this as a PDE-constrained optimization problem based on a mechanical model defined for two-dimensional representations of the microtissues. Our optimization predicts active tension generated by the tissue as well as the fiber direction angle distribution. We used synthetic as well as experimental data to investigate the performance of our inversion protocol. Next, we employed the procedure to evaluate active tension changes in drug escalation studies of the inotropes omecamtiv mecarbil and Bay K8644. For both drug compounds, we observed a comparable increase in displacement, strain, and model-predicted active strain values upon higher drug doses. The estimated active tension was observed to be highest in the middle part of the tissue, and the fiber direction was mostly aligned with the longitudinal direction of the tissue. The computational framework presented here allows for spatiotemporal estimation of active tension in cardiac microtissues based on optical measurements. In the future, such methodologies might develop into valuable tools in drug development protocols.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"41 4","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/cnm.70034","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143865799","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":"Validation of a Microfluidic Device Prototype for Cancer Detection and Identification: Circulating Tumor Cells Classification Based on Cell Trajectory Analysis Leveraging Cell-Based Modeling and Machine Learning","authors":"Rifat Rejuan, Eugenio Aulisa, Wei Li, Travis Thompson, Sanjoy Kumar, Suncica Canic, Yifan Wang","doi":"10.1002/cnm.70037","DOIUrl":"https://doi.org/10.1002/cnm.70037","url":null,"abstract":"<div>\u0000 \u0000 <p>Microfluidic devices (MDs) present a novel method for detecting circulating tumor cells (CTCs), enhancing the process through targeted techniques and visual inspection. However, current approaches often yield heterogeneous CTC populations, necessitating additional processing for comprehensive analysis and phenotype identification. These procedures are often expensive, time-consuming, and need to be performed by skilled technicians. In this study, we investigate the potential of a cost-effective and efficient hyperuniform micropost MD approach for CTC classification. Our approach combines mathematical modeling of fluid–structure interactions in a simulated microfluidic channel with machine learning techniques. Specifically, we developed a cell-based modeling framework to assess CTC dynamics in erythrocyte-laden plasma flow, generating a large dataset of CTC trajectories that account for two distinct CTC phenotypes. Convolutional neural network (CNN) and recurrent neural network (RNN) were then employed to analyze the dataset and classify these phenotypes. The results demonstrate the potential effectiveness of the hyperuniform micropost MD design and analysis approach in distinguishing between different CTC phenotypes based on cell trajectory, offering a promising avenue for early cancer detection.</p>\u0000 </div>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"41 4","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143871777","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}
Laura Engelhardt, Renate Sachse, Rainer Burgkart, Wolfgang A. Wall
{"title":"Constitutive Models for Active Skeletal Muscle: Review, Comparison, and Application in a Novel Continuum Shoulder Model","authors":"Laura Engelhardt, Renate Sachse, Rainer Burgkart, Wolfgang A. Wall","doi":"10.1002/cnm.70036","DOIUrl":"https://doi.org/10.1002/cnm.70036","url":null,"abstract":"<p>The shoulder joint is one of the functionally and anatomically most sophisticated articular systems in the human body. Both complex movement patterns and the stabilization of the highly mobile joint rely on intricate three-dimensional interactions among various components. Continuum-based finite element models can capture such complexity and are thus particularly relevant in shoulder biomechanics. Considering their role as active joint stabilizers and force generators, skeletal muscles require special attention regarding their constitutive description. In this contribution, we propose a constitutive description to model active skeletal muscle within complex musculoskeletal systems, focusing on a novel continuum shoulder model. Based on a thorough review of existing material models, we select an active stress, an active strain, and a generalized active strain approach and combine the most promising and relevant features in a novel material model. We discuss the four models considering physiological, mathematical, and computational aspects, including the applied activation concepts, biophysical principles of force generation, and arising numerical challenges. To establish a basis for numerical comparison, we identify the material parameters based on experimental stress–strain data obtained under multiple active and passive loading conditions. Using the example of a fusiform muscle, we investigate force generation, deformation, and kinematics during active isometric and free contractions. Eventually, we demonstrate the applicability of the proposed material model in a novel continuum mechanical model of the human shoulder, exploring the role of rotator cuff contraction in joint stabilization.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"41 4","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-04-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/cnm.70036","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143865609","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":"New Surgical Devices for Closed Reduction of Frontal Sinus Bone Fracture","authors":"Daehan Wi, Hoyul Lee, Woo Shik Jeong, Jaesoon Choi, Youngjin Moon, Jong Woo Choi","doi":"10.1002/cnm.70042","DOIUrl":"https://doi.org/10.1002/cnm.70042","url":null,"abstract":"<p>Traditional surgical interventions for frontal sinus fractures necessitate a cut on the forehead skin, and extant closed reduction techniques aimed at enhancing accessibility continue to grapple with secure tool fixation, stable bone elevation, and screw breakage risk. To address these challenges and augment surgical efficiency, this study introduces novel surgical devices. Design parameters for models with spiral or L-shaped tips are established, considering practical medical requirements and constraints, and subsequently validated through finite element method numerical simulations using commercial software, Ansys. Four spiral-type prototypes are constructed, and three scenarios for each prototype, varying in projection distance from the device handle to the bone-device contact point, are examined via nonlinear simulation analyses. For the L-shaped type, three prototypes are developed, and static analyses are conducted for four scenarios per prototype, differing in traction force locations, based on another simulation result concerning moments of inertia calculation with a force boundary condition unlike pressure. Maximum stress results under a specific force are analyzed, and the maximum permissible force is determined under the most unfavorable force application condition. Simulation outcomes indicate that the spiral type offers greater applicability with less force to lift multiple bones, while the L-shaped type is more suitable under bone hardening conditions.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"41 4","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-04-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/cnm.70042","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143852754","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":"Comparative Analysis of Various Cavosurface Margins in Class II Restorations Using 3D Finite Element Method","authors":"Zuzanna Apel, Behzad Vafaeian, Joanna Zarzecka, Jenna Wuzinski, Derek B. Apel","doi":"10.1002/cnm.70041","DOIUrl":"https://doi.org/10.1002/cnm.70041","url":null,"abstract":"<p>The primary aim of these analyses was to evaluate the mechanical characteristics of the restored proximal surface of the lower first molar by comparing four different preparation designs: (a) slot preparation, (b) slot preparation with bevel, (c) slot preparation with bevel and rounded proximal box corners (RPBC), and (d) slot preparation with bevel, rounded proximal box corners, and gingival bevel (GB). The finite element method was utilized to assess various load scenarios applied to slot and bevelled restorations prepared using adhesive restorative materials. The numerical analysis revealed higher tensile stresses by up to 15 MPa when normal traction was applied at the interface between enamel and slot preparations than at the interface between enamel and bevelled preparations. However, the beveled restorations showed increased shear stresses in their thin beveled regions. The results imply a risk of separation for slot restorations. Conversely, incorporating a bevel (with or without RPBC and GB) significantly decreased normal stresses on the restoration edge and shifted it predominantly to compressive stresses. Thus, bevelled restorations may be less prone to debonding at their edges under occlusal loads. However, they may still be susceptible to shear debonding when locally loaded on their thin-beveled regions.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"41 4","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-04-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/cnm.70041","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143840844","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}
Vishnu S. Emani, Caglar Ozturk, Manisha Singh, Carly Long, Summer Duffy, Danielle Gottlieb Sen, Ellen T. Roche, Wesley B. Baker
{"title":"Finite Element Modeling of Abdominal Near-Infrared Spectroscopy for Infant Splanchnic Oximetry","authors":"Vishnu S. Emani, Caglar Ozturk, Manisha Singh, Carly Long, Summer Duffy, Danielle Gottlieb Sen, Ellen T. Roche, Wesley B. Baker","doi":"10.1002/cnm.70035","DOIUrl":"https://doi.org/10.1002/cnm.70035","url":null,"abstract":"<p>Abdominal near-infrared spectroscopy (NIRS) holds promise for early detection of necrotizing enterocolitis and other infant pathologies prior to irreversible injury, but the optimal NIRS sensor design is not well defined. In this study, we develop and demonstrate a computational method to evaluate NIRS sensor designs for infant splanchnic oximetry. We used a finite element (FE) approach to simulate near-infrared light transport through a 3D model of the infant abdomen constructed from computed tomography (CT) images. The simulations enable the measurement of the contrast-to-noise ratio (CNR) for splanchnic oximetry, given a specific NIRS sensor design. A key design criterion is the sensor's source–detector distance (SDD). We calculated the CNR as a function of SDD for two sensor positions near the umbilicus. Contrast-to-noise was maximal at SDDs between 4 and 5 cm, and comparable between sensor positions. Sensitivity to intestinal tissue also exceeded sensitivity to superficial adipose tissue in the 4–5 cm range. FE modeling of abdominal NIRS signals provides a means for rapid and thorough evaluation of sensor designs for infant splanchnic oximetry. By informing optimal NIRS sensor design, the computational methods presented here can improve the reliability and applicability of infant splanchnic oximetry.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":"41 4","pages":""},"PeriodicalIF":2.2,"publicationDate":"2025-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/cnm.70035","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143835926","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}