Amadeus M Gebauer, Martin R Pfaller, Jason M Szafron, Wolfgang A Wall
{"title":"Adaptive integration of history variables in constrained mixture models for organ-scale growth and remodeling.","authors":"Amadeus M Gebauer, Martin R Pfaller, Jason M Szafron, Wolfgang A Wall","doi":"10.1002/cnm.3869","DOIUrl":"10.1002/cnm.3869","url":null,"abstract":"<p><p>In the last decades, many computational models have been developed to predict soft tissue growth and remodeling (G&R). The constrained mixture theory describes fundamental mechanobiological processes in soft tissue G&R and has been widely adopted in cardiovascular models of G&R. However, even after two decades of work, large organ-scale models are rare, mainly due to high computational costs (model evaluation and memory consumption), especially in long-range simulations. We propose two strategies to adaptively integrate history variables in constrained mixture models to enable large organ-scale simulations of G&R. Both strategies exploit that the influence of deposited tissue on the current mixture decreases over time through degradation. One strategy is independent of external loading, allowing the estimation of the computational resources ahead of the simulation. The other adapts the history snapshots based on the local mechanobiological environment so that the additional integration errors can be controlled and kept negligibly small, even in G&R scenarios with severe perturbations. We analyze the adaptively integrated constrained mixture model on a tissue patch for a parameter study and show the performance under different G&R scenarios. To confirm that adaptive strategies enable large organ-scale examples, we show simulations of different hypertension conditions with a real-world example of a biventricular heart discretized with a finite element mesh. In our example, adaptive integrations sped up simulations by a factor of three and reduced memory requirements to one-sixth. The reduction of the computational costs gets even more pronounced for simulations over longer periods. Adaptive integration of the history variables allows studying more finely resolved models and longer G&R periods while computational costs are drastically reduced and largely constant in time.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":null,"pages":null},"PeriodicalIF":2.2,"publicationDate":"2024-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142299985","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}
Sarah Iaquinta, Shahram Khazaie, Samer Albanna, Sylvain Fréour, Frédéric Jacquemin
{"title":"PREPRINT Machine Learning for the Sensitivity Analysis of a Model of the Cellular Uptake of Nanoparticles for the Treatment of Cancer.","authors":"Sarah Iaquinta, Shahram Khazaie, Samer Albanna, Sylvain Fréour, Frédéric Jacquemin","doi":"10.1002/cnm.3878","DOIUrl":"https://doi.org/10.1002/cnm.3878","url":null,"abstract":"<p><p>Experimental studies on the cellular uptake of nanoparticles (NPs), useful for the investigation of NP-based drug delivery systems, are often difficult to interpret due to the large number of parameters that can contribute to the phenomenon. It is therefore of great interest to identify insignificant parameters to reduce the number of variables used for the design of experiments. In this work, a model of the wrapping of elliptical NPs by the cell membrane is used to compare the influence of the aspect ratio of the NP, the membrane tension, the NP-membrane adhesion, and its variation during the interaction with the NP on the equilibrium state of the wrapping process. Several surrogate models, such as Kriging, Polynomial Chaos Expansion (PCE), and artificial neural networks (ANN) have been built and compared to emulate the computationally expensive model. Only the ANN-based model outperformed the other approaches by providing much better predictivity metrics and could therefore be used to compute the sensitivity indices. Our results showed that the NP's aspect ratio, the initial NP-membrane adhesion, the membrane tension, and the delay for the increase of the NP-membrane adhesion after receptor dynamics are the main contributors to the cellular internalization of the NP, while the influence of other parameters is negligible.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":null,"pages":null},"PeriodicalIF":2.2,"publicationDate":"2024-10-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142523500","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":"Prediction of Transport and Deposition of Porous Particles in the Respiratory System Using Eulerian-Lagrangian Approach.","authors":"Sajad Eshaghi, Hassan Khaleghi, Reza Maddahian","doi":"10.1002/cnm.3873","DOIUrl":"https://doi.org/10.1002/cnm.3873","url":null,"abstract":"<p><p>Deep lung delivery is crucial for respiratory disease treatment. Although nano and submicron particles exhibited a good deposition efficiency in deep regions of the lung, powder nonuniformity and particle agglomeration reduce their efficiency. Inhalation of porous particles (PPs) can overcome the mentioned challenges due to their larger size and low-density. The present study numerically investigates the deposition and penetration efficiency of orally inhaled PPs. A revised drag coefficient was implemented for PP transport. A realistic mouth-throat to the fifth generation of the lung was reconstructed from CT-scan images. A dilute suspension of uniformly distributed particles was considered at three inhalation flow rates (15, 30, and 45 L/min). Governing equations of the flow field and particle transport are solved using an Eulerian-Lagrangian approach. The results demonstrate that inhaling PPs significantly reduces the total and regional deposition of particles. There was also a critical porosity value under moderate and high inhalation flow rates for large PPs. Below this critical value, PP deposition efficiency substantially decreases. Additionally, it was also found that under low inhalation flow rates, the impact of porosity value is negligible. Almost 95% of the PPs penetrate the lower branches. These findings provide particle engineers and pharmaceutics with profound insight into developing novel inhalation techniques and drug delivery methods for deep lung delivery.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":null,"pages":null},"PeriodicalIF":2.2,"publicationDate":"2024-10-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142512250","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}
Bettine G van Willigen, M Beatrijs van der Hout-van der Jagt, Peter H M Bovendeerd, Wouter Huberts, Frans N van de Vosse
{"title":"A Multiscale Mathematical Model for the Fetal Blood Circulation of the Second Half of Pregnancy.","authors":"Bettine G van Willigen, M Beatrijs van der Hout-van der Jagt, Peter H M Bovendeerd, Wouter Huberts, Frans N van de Vosse","doi":"10.1002/cnm.3877","DOIUrl":"https://doi.org/10.1002/cnm.3877","url":null,"abstract":"<p><p>Doppler ultrasound is a commonly used method to assess hemodynamics of the fetal cardiovascular system and to monitor the well-being of the fetus. Indices based on the velocity profile are often used for diagnosis. However, precisely linking these indices to specific underlying physiology factors is challenging. Several influences, including wave reflections, fetal growth, vessel stiffness, and resistance distal to the vessel, contribute to these indices. Understanding these data is essential for making informed clinical decisions. Mathematical models can be used to investigate the relation between velocity profiles and physiological properties. This study presents a mathematical model designed to simulate velocity wave propagation throughout the fetal cardiovascular system, facilitating the assessment of factors influencing velocity-based indices. The model combines a one-fiber model of the heart with a 1D wave propagation model describing the larger vessels of the circulatory system and a lumped parameter model for the microcirculation. Fetal growth from 20 to 40 weeks of gestational age is incorporated by adjusting cardiac and circulatory parameter settings according to scaling laws. The model's results, including cardiac function, cardiac output distribution, and volume distribution, show a good agreement with literature studies for a growing healthy fetus from 20 to 40 weeks. In addition, Doppler indices are simulated in various vessels and agree with literature as well. In conclusion, this study introduces a novel closed-loop 0D-1D mathematical model that has been verified against literature studies. This model offers a valuable platform for analyzing factors influencing velocity-based indices in the fetal cardiovascular system.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":null,"pages":null},"PeriodicalIF":2.2,"publicationDate":"2024-10-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142512249","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":"Numerical Investigation of the Blast-Induced Injuries Using an Open-Source Detailed Human Model.","authors":"Alberto Morena, Lorenzo Peroni, Martina Scapin","doi":"10.1002/cnm.3879","DOIUrl":"https://doi.org/10.1002/cnm.3879","url":null,"abstract":"<p><p>Blasts are a threat both in military and civil contexts due not only to explosive devices but also to gas leakages or other accidents. Numerical models could aid to plan response strategies in the short and long term. Nevertheless, due to modeling complexities, a standardized computational framework has not been established yet. In this challenging context, the present study assesses the prediction of blast-induced traumas by using the total human model for safety (THUMS) human model, which has never been attempted before to the authors knowledge. The pedestrian model is publicly available, hence the demonstration of its suitability to predict blast injuries could benefit the establishment of a common modeling framework. Therefore, the THUMS human model was exposed to different blast scenarios both in free field and partially confined spaces and the response of vital organs was investigated. Trauma patterns to internal organs of the THUMS were consistent with available experimental data and injury thresholds. In conclusion, THUMS open-source human model demonstrated its validity to reproduce primary blast-related injuries, addressing the development of standardization of numerical simulations of human response to explosions.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":null,"pages":null},"PeriodicalIF":2.2,"publicationDate":"2024-10-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142479646","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":"Biomechanical significance of intervertebral discs on growthplate stresses in scoliotic trunks following unilateral muscle weakening: A hybrid approach of finite element and musculoskeletal modeling.","authors":"Zeinab Kamal","doi":"10.1002/cnm.3863","DOIUrl":"https://doi.org/10.1002/cnm.3863","url":null,"abstract":"<p><p>This study aimed to ascertain the relevance of intervertebral discs (IVD) in the stress distribution on growthplates (GPs) of a trunk model with adolescent idiopathic scoliosis (AIS) following a unilateral weakening of muscles. A thoracolumbar spine finite element (FE) model of a young female healthy and an AIS spine comprising GPs linked to the T12 through sacrum vertebrae. Two scenarios of including (FEI) and excluding (FEE) IVDs were considered. Then, using optimization-driven musculoskeletal models of the AIS and healthy trunks, the FE models were examined under subject-specific muscle forces and gravity loads. Results of this study demonstrate that when IVDs included in the FE model, an increase, ranging from 0.2 to 1.7 MPa, with the highest value occurring at the apex of the AIS model, in the von Mises stresses in the GPs. The ratio of 1.5 was found for the maximum von-Mises stress value on the most tilted GP in the FEI over the FEE model. Unilateral paralysis of muscles caused a reduction of 50% and 63% in the von Mises stress ratio of the concave-over-convex side of the most tilted GP in the FEI and FEE models of the AIS spine with healthy muscles, respectively. The intradiscal pressures, found for FEE and FEI models, assented to recent in-vivo investigations. Nonetheless, employing IVDs in the simulations provides an indispensable tool to anticipate the effects of neuromuscular disorders on GP stresses in an AIS spine and predict deformity progression during growth.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":null,"pages":null},"PeriodicalIF":2.2,"publicationDate":"2024-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142479645","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}
Jacques-Ezechiel N'Guessan, Muhammad Hassaan Ahmed, Matthew Leineweber, Sachin Goyal
{"title":"Piloting a Novel Computational Framework for Identifying Prosthesis-Specific Contributions to Gait Deviations.","authors":"Jacques-Ezechiel N'Guessan, Muhammad Hassaan Ahmed, Matthew Leineweber, Sachin Goyal","doi":"10.1002/cnm.3876","DOIUrl":"https://doi.org/10.1002/cnm.3876","url":null,"abstract":"<p><p>This paper introduces a novel computational framework for evaluating above-knee prostheses, addressing a major challenge in gait deviation studies: distinguishing between prosthesis-specific and patient-specific contributions to gait deviations. This innovative approach utilizes three separate computational models to quantify the changes in gait dynamics necessary to achieve a set of ideal gait kinematics across different prosthesis designs. The pilot study presented here employs a simple two-dimensional swing-phase model to conceptually demonstrate how the outcomes of this three-model framework can assess the extent to which prosthesis design impacts a user's ability to replicate the dynamics of able-bodied gait. Furthermore, this framework offers potential for optimizing passive prosthetic devices for individual patients, thereby reducing the need for real-life experiments, clinic visits, and overcoming rehabilitation challenges.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":null,"pages":null},"PeriodicalIF":2.2,"publicationDate":"2024-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142401838","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}
Rajat Mishra, Sagar Kumar Deb, Swasti Chakrabarty, Manojit Das, Monalisa Das, Sushanta Kumar Panda, Chandra Shekhar Tiwary, Amit Arora
{"title":"Human Mastication Analysis-A DEM Based Numerical Approach.","authors":"Rajat Mishra, Sagar Kumar Deb, Swasti Chakrabarty, Manojit Das, Monalisa Das, Sushanta Kumar Panda, Chandra Shekhar Tiwary, Amit Arora","doi":"10.1002/cnm.3875","DOIUrl":"https://doi.org/10.1002/cnm.3875","url":null,"abstract":"<p><p>Mastication is an essential and preliminary step of the digestion process involving fragmentation and mixing of food. Controlled muscle movement of jaws with teeth executes crushing, leading towards fragmentation of food particles. Understanding various parameters involved with the process is essential to solve any biomedical complication in the area of interest. However, exploring and analyzing such process flow through an experimental route is challenging and inefficient. Computational techniques such as discrete element numerical modeling can effectively address such problems. The current work employs the Discrete Element Method (DEM) as a numerical modeling technique to simulate the human mastication process. Tavares and Ab-T10 breakage models coupled with Gaudin Schumann and Incomplete Beta fragment distribution models have been implemented to analyze the fragmental distribution of food particles. The effect of particle shape (spherical, polyhedron, and faceted cylinder), size (aspect ratio), and orientation (vertical and horizontal) on breakage and fragment distribution is analyzed. To account for the elastic-plastic behavior and moisture content in food particles, modifications has been made in breakage models by incorporating numerical softening factor and adhesion force. The study demonstrates how numerical modeling techniques can be utilized to analyze the mastication process involving multiple process parameters.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":null,"pages":null},"PeriodicalIF":2.2,"publicationDate":"2024-10-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142401837","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":"Biomechanics of the Human Knee Joint in Maximum Voluntary Isometric Flexion: Study of Changes in Applied Moment, Agonist-Antagonist Participations, Joint Center, and Flexion Angle.","authors":"Pooya Salehi, Aboulfazl Shirazi-Adl, Farshid Ghezelbash","doi":"10.1002/cnm.3874","DOIUrl":"https://doi.org/10.1002/cnm.3874","url":null,"abstract":"<p><p>Estimation of the knee joint strength by maximum voluntary isometric contraction (MVIC) is a common practice to assess strength, coordination, safety to return to work or engage in sports after an injury, and to evaluate the efficacy of treatment modalities and rehabilitation strategies. In this study, we utilize a previously validated coupled finite element-musculoskeletal model of the lower extremity to explore the sensitivity of output measures (posterior cruciate ligament [PCL]/muscle/contact forces and passive moments) in knee MVIC flexion exercises at seated position. To do so, at three knee flexion angles (KFA), input measures (resistance moment and contribution moments of quadriceps and gastrocnemii) were varied at four levels each using the Taguchi design of experiment. Our findings reveal significant increases in PCL forces with KFA (p < 0.01), net MVIC moment (p < 0.01), and resistance moment of quadriceps (p < 0.01). In contrast, they drop at larger activity in gastrocnemii (p < 0.01). Tibiofemoral (TF) contact forces increase with the net MVIC moment (p < 0.01). The passive knee flexion moment, while highly dependent on the location at which computed, also increases with the net MVIC moment (p < 0.01). Changes in KFA, MVIC moment, and proportions thereof carried by quadriceps and/or gastrocnemii substantially affect biomechanics of the joint. Compared with level walking and stair ascent, slightly larger contact forces/stresses and much greater PCL forces are computed. This study improves our understanding of the knee joint behavior during MVIC in effective evaluation and rehabilitation interventions. Besides, it emphasizes the importance of positioning the joint center in model studies.</p>","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":null,"pages":null},"PeriodicalIF":2.2,"publicationDate":"2024-10-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142394896","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":"Effect of electromagnetic middle-ear implant simulating sites on the stapes spatial motion: A finite element analysis.","authors":"Yixiang Zhang,Houguang Liu,Lei Zhou,Jianhua Yang,Wen Liu,Shanguo Yang,Xinsheng Huang","doi":"10.1002/cnm.3871","DOIUrl":"https://doi.org/10.1002/cnm.3871","url":null,"abstract":"The electromagnetic middle-ear implant (MEI) is a new type of hearing device for addressing sensorineural and mixed hearing loss. The hearing compensation effect of the MEI varies depending on the transducer stimulation sites. This paper investigates the impact of transducer stimulation sites on MEI performance by analyzing stapes spatial motion. Firstly, we constructed a human-ear finite element model based on micro-CT scanning and inverse molding techniques. This model was validated by comparing its predictions of stapes spatial motion and cochlear response with experimental data. Then, stimulation force was applied at four common sites: umbo, incus body, incus long process and stapes to simulate the electromagnetic transducer. Results show that at low and middle frequencies, stapes-stimulating and incus-long-process-stimulating produce similar spatial motion to normal hearing; at high frequencies, incus-body-stimulating produces similar results to normal hearing. The equivalent sound pressure level generated by the stapes piston motion is less sensitive to the stimulation direction than that deduced by the stapes rocking motion.","PeriodicalId":50349,"journal":{"name":"International Journal for Numerical Methods in Biomedical Engineering","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-09-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142254002","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}