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

Over the last couple of decades, image-based computational fluid dynamics (CFD) has revolutionized cardiovascular research by uncovering hidden features of wall strain, impact of vortices, and its use in treatment planning, as examples, that were simply not evident in the gold-standard catheterization studies done previously. In the work presented here, we have applied magnetic resonance imaging (MRI)-based CFD to study bladder voiding and to demonstrate the feasibility and potential of this approach. We used 3D dynamic MRI to image the bladder and urethra during voiding. A surface mesh processing tool was developed to process the bladder wall prior to executing a wall-motion driven CFD simulation of the bladder and urethra. The obtained flow rate and pressure were used to calculate urodynamic nomograms, which are currently used in the clinical setting to assess bladder voiding dysfunction. These nomograms concluded that our healthy volunteer has an unobstructed bladder and normal contractility. We calculated the work done to void the bladder and propose this as an additional quantitative metric to comprehensively assess bladder function. Further, we discuss the areas that would improve this relatively new methodology of image-based CFD in urodynamics.

We develop computational mechanical modeling and methods for the analysis and simulation of the motions of a human body. This type of work is crucial in many aspects of human life, ranging from comfort in riding, the motion of aged persons, sports performance and injuries, and many ergonomic issues. A prevailing approach for human motion studies is through lumped parameter models containing discrete masses for the parts of the human body with empirically determined spring, mass, damping coefficients. Such models have been effective to some extent; however, a much more faithful modeling method is to model the human body as it is, namely, as a continuum. We present this approach, and for comparison, we choose two digital CAD models of mannequins for a standing human body, one from the versatile software package LS-DYNA and another from open resources with some of our own adaptations. Our basic view in this paper is to regard human motion as a perturbation and vibration from an equilibrium position which is upright standing. A linear elastodynamic model is chosen for modal analysis, but a full nonlinear viscoelastoplastic extension is possible for full-body simulation. The motion and vibration of these two mannequin models is analyzed by modal analysis, where the normal vibration modes are determined. LS-DYNA is used as the supercomputing and simulation platform. Four sets of low-frequency modes are tabulated, discussed, visualized, and compared. Higher frequency modes are also selectively displayed. We have found that these modes of motion and vibration form intrinsic basic modes of biomechanical motion of the human body. This view is supported by our finding of the upright walking motion as a low-frequency mode in modal analysis. Such a “walking mode” shows the in-phase and out-of-phase movements between the legs and arms on the left and right sides of a human body, implying that this walking motion is spontaneous, likely not requiring any directives from the brain. Dynamic motions of CAD mannequins are also simulated by drop tests for comparisons and the validity of the models is discussed through Fourier frequency analysis. All computed modes of motion are collected in several sets of video animations for ease of visualization. Samples of LS-DYNA computer codes are also included for possible use by other researchers.

Virtual compression tests based on finite element analysis are representative noninvasive methods to evaluate bone strength. However, owing to the characteristic porous structure of bones, the material obtained from micro-computed tomography images in the finite-element model is not uniformly distributed. These characteristics cause differences in the apparent elastic moduli depending on the boundary conditions and affect the accuracy of bone-strength evaluation. Therefore, this study aimed to evaluate and compare the apparent elastic moduli under various, virtual-compression and shear-test boundary conditions. Four, nonuniform models were constructed with increasing model complexity. For representative boundary conditions, two, different, testing directions, and constrained surfaces were applied. As a result, the apparent elastic moduli of the nonuniform model varied up to 55.2% based on where the constrained surface was located in the single-end-cemented condition. Additionally, when connectivity in the test direction was lost, the accuracy of the apparent elastic moduli was low. A graphical comparison showed that the equivalent-stress distribution was more advantageous for analyzing load transferability and physical behavior than the strain-energy distribution. These results clearly show that the prediction accuracy of the apparent elastic moduli can be guaranteed if the boundary condition on the constraint and loading surfaces of the nonuniform model are applied symmetrically and the connectivity of the elements in the testing direction is well maintained. This study will aid in precision improvement of bone-strength-indicator determination for osteoporosis prevention.

An epicardial approach is often used in radiofrequency (RF) catheter ablation to ablate ventricular tachycardia when an endocardial approach fails. Our objective was to analyze the effect of the position of the dispersive patch (DP) on lesion size using computer modeling during epicardial approach. We compared the posterior position (patient's back), commonly used in clinical practice, to the anterior position (patient's chest). The model considered ventricular wall thicknesses between 4 and 8 mm, and electrode insertion depths between .3 and .7 mm. RF pulses were simulated with 20 W of power for 30 s duration. Statistically significant differences (P < .001) were found between both DP positions in terms of baseline impedance, RF current (at 15 s) and thermal lesion size. The anterior position involved lower impedance (130.8 ± 4.7 vs. 146.2 ± 4.9 Ω) and a higher current (401.5 ± 5.6 vs. 377.5 ± 5.1 mA). The anterior position created lesion sizes larger than the posterior position: 8.9 ± 0.4 vs. 8.4 ± 0.4 mm in maximum width, 8.6 ± 0.4 vs. 8.1 ± 0.4 mm in surface width, and 4.5 ± 0.4 vs. 4.3 ± 0.4 mm in depth. Our results suggest that: (1) the redirection of the RF currents due to repositioning the PD has little impact on lesion size and only affects baseline impedance, and (2) the differences in lesion size are only 0.5 mm wider and 0.2 mm deeper for the anterior position, which does not seem to have a clinical impact in the context of VT ablation.

Infrared thermography is gaining relevance in breast cancer assessment. For this purpose, breast segmentation in thermograms is an important task for performing automatic image analysis and detecting possible temperature changes that indicate the presence of malignancy. However, it is not a simple task since the breast limit borders, especially the top borders, often have low contrast, making it difficult to isolate the breast area. Several algorithms have been proposed for breast segmentation, but these highly depend on the contrast at the lower breast borders and on filtering algorithms to remove false edges. This work focuses on taking advantage of the distinctive inframammary shape to simplify the definition of the lower breast border, regardless of the contrast level, which indeed also provides a strong anatomical reference to support the definition of the poorly marked upper boundary of the breasts, which has been one of the major challenges in the literature. In order to demonstrate viability of the proposed technique for an automatic breast segmentation, we applied it to a database with 180 thermograms and compared their results with those reported by others in the literature. We found that our approach achieved a high performance, in terms of Intersection over Union of 0.934, even higher than that reported by artificial intelligence algorithms. The performance is invariant to breast sizes and thermal contrast of the images.

Intracranial aneurysms (IAs) pose severe health risks influenced by hemodynamics. This study focuses on the intricate characterization of hemodynamic conditions within the IA walls and their influence on bleb development, aiming to enhance understanding of aneurysm stability and the risk of rupture. The methods emphasized utilizing a comprehensive dataset of 359 IAs and 213 IA blebs from 268 patients to reconstruct patient-specific vascular models, analyzing blood flow using finite element methods to solve the unsteady Navier–Stokes equations, the segmentation of aneurysm wall subregions and the hemodynamic metrics wall shear stress (WSS), its metrics, and the critical points in WSS fields were computed and analyzed across different aneurysm subregions defined by saccular, streamwise, and topographical divisions. The results revealed significant variations in these metrics, correlating distinct hemodynamic environments with wall features on the aneurysm walls, such as bleb formation. Critical findings indicated that regions with low WSS and high OSI, particularly in the body and central regions of aneurysms, are prone to conditions that promote bleb formation. Conversely, areas exposed to high WSS and positive divergence, like the aneurysm neck, inflow, and outflow regions, exhibited a different but substantial risk profile for bleb development, influenced by flow impingements and convergences. These insights highlight the complexity of aneurysm behavior, suggesting that both high and low-shear environments can contribute to aneurysm pathology through distinct mechanisms.

The aortic valve (AV) is crucial for cardiovascular (CV) hemodynamic, impacting cardiac output (CO) and left ventricular volumetric flow rate (LVQ). Its nonlinear behavior challenges standard LVQ prediction methods as well as CO one. This study presents a novel approach for modeling the AV in the CV system, offering an improved method for estimating crucial parameters like LVQ across various AV conditions, including aortic stenosis (AS). The model, based on AV channel length during the entire cardiac phase, introduces a time-varying AV resistance (TV-AVR) parameterized by the pressure ratio across the AV and LVQ, enabling the simulation of both healthy and AS-related conditions. To validate this model, in vitro measurements are compared using a hybrid mock circulatory loop device. An unconventional use of a convolutional neural network (CNN) corrects the model's estimates, eliminating the need for labeled datasets. This approach, incorporating real-time learning and transforming 1-D CV signals into 2-D tensors, significantly improves the accuracy of LVQ measurements, achieving an error rate of less than 3.41 ± 4.84% for CO in healthy conditions and 2.83 ± 1.35% in AS cases—a 33.13% enhancement over linear diode models. These results underscore the potential of this approach for enhancing the diagnosis, prediction, and treatment of AV diseases. The key contributions of the proposed method encompass nonlinear TV-AVR estimation, investigation of transient CV responses, prediction of instantaneous CO, development of a flexible framework for noninvasive measurements integration, and the introduction of an adjustable resistance model using an extended Kalman filter (EKF) and CNN combination, all without requiring labeled data.

This study aims to investigate how inert gas affects the partial pressure of alveolar and venous blood using a fast and accurate operator splitting method (OSM). Unlike previous complex methods, such as the finite element method (FEM), OSM effectively separates governing equations into smaller sub-problems, facilitating a better understanding of inert gas transport and exchange between blood capillaries and surrounding tissue. The governing equations were discretized with a fully implicit finite difference method (FDM), which enables the use of larger time steps. The model employed partial differential equations, considering convection-diffusion in blood and only diffusion in tissue. The study explores the impact of initial arterial pressure, breathing frequency, blood flow velocity, solubility, and diffusivity on the partial pressure of inert gas in blood and tissue. Additionally, the effects of anesthetic inert gas and oxygen on venous blood partial pressure were analyzed. Simulation results demonstrate that the high solubility and diffusivity of anesthetic inert gas lead to its prolonged presence in blood and tissue, resulting in lower partial pressure in venous blood. These findings enhance our understanding of inert gas interaction with alveolar/venous blood, with potential implications for medical diagnostics and therapies.

A high failure rate is associated with fracture plates in proximal humerus fractures. The causes of failure remain unclear due to the complexity of the problem including the number and position of the screws, their length and orientation in the space. Finite element (FE) analysis has been used for the analysis of plating of proximal humeral fractures, but due to computational costs is unable to fully explore all potential screw combinations. Surrogate modelling is a viable solution, having the potential to significantly reduce the computational cost whilst requiring a moderate number of training sets. This study aimed to develop adaptive neural network (ANN)-based surrogate models to predict the strain in the humeral bone as a result of changing the length of the screws. The ANN models were trained using data from FE simulations of a single humerus, and after defining the best training sample size, multiple and single-output models were developed. The best performing ANN model was used to predict all the possible screw length configurations. The ANN predictions were compared with the FE results of unseen data, showing a good correlation (R² = 0.99) and low levels of error (RMSE = 0.51%–1.83% strain). The ANN predictions of all possible screw length configurations showed that the screw that provided the medial support was the most influential on the predicted strain. Overall, the ANN-based surrogate model accurately captured bone strains and has the potential to be used for more complex problems with a larger number of variables.

The mechanisms behind intracranial aneurysm formation and rupture are not fully understood, with factors such as location, patient demographics, and hemodynamics playing a role. Additionally, the significance of anatomical features like blebs in ruptures is debated. This highlights the necessity for comprehensive research that combines patient-specific risk factors with a detailed analysis of local hemodynamic characteristics at bleb and rupture sites. Our study analyzed 359 intracranial aneurysms from 268 patients, reconstructing patient-specific models for hemodynamic simulations based on 3D rotational angiographic images and intraoperative videos. We identified aneurysm subregions and delineated rupture sites, characterizing blebs and their regional overlap, employing statistical comparisons across demographics, and other risk factors. This work identifies patterns in aneurysm rupture sites, predominantly at the dome, with variations across patient demographics. Hypertensive and anterior communicating artery (ACom) aneurysms showed specific rupture patterns and bleb associations, indicating two pathways: high-flow in ACom with thin blebs at impingement sites and low-flow, oscillatory conditions in middle cerebral artery (MCA) aneurysms fostering thick blebs. Bleb characteristics varied with gender, age, and smoking, linking rupture risks to hemodynamic factors and patient profiles. These insights enhance understanding of the hemodynamic mechanisms leading to rupture events. This analysis elucidates the role of localized hemodynamics in intracranial aneurysm rupture, challenging the emphasis on location by revealing how flow variations influence stability and risk. We identify two pathways to wall failure—high-flow and low-flow conditions—highlighting the complexity of aneurysm behavior. Additionally, this research advances our knowledge of how inherent patient-specific characteristics impact these processes, which need further investigation.