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

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.

Computational models of the cardiovascular system are increasingly used for the diagnosis, treatment, and prevention of cardiovascular disease. Before being used for translational applications, the predictive abilities of these models need to be thoroughly demonstrated through verification, validation, and uncertainty quantification. When results depend on multiple uncertain inputs, sensitivity analysis is typically the first step required to separate relevant from unimportant inputs, and is key to determine an initial reduction on the problem dimensionality that will significantly affect the cost of all downstream analysis tasks. For computationally expensive models with numerous uncertain inputs, sample-based sensitivity analysis may become impractical due to the substantial number of model evaluations it typically necessitates. To overcome this limitation, we consider recently proposed Multifidelity Monte Carlo estimators for Sobol’ sensitivity indices, and demonstrate their applicability to an idealized model of the common carotid artery. Variance reduction is achieved combining a small number of three-dimensional fluid–structure interaction simulations with affordable one- and zero-dimensional reduced-order models. These multifidelity Monte Carlo estimators are compared with traditional Monte Carlo and polynomial chaos expansion estimates. Specifically, we show consistent sensitivity ranks for both bi- (1D/0D) and tri-fidelity (3D/1D/0D) estimators, and superior variance reduction compared to traditional single-fidelity Monte Carlo estimators for the same computational budget. As the computational burden of Monte Carlo estimators for Sobol’ indices is significantly affected by the problem dimensionality, polynomial chaos expansion is found to have lower computational cost for idealized models with smooth stochastic response.

Microwave ablation has become a viable alternative for cancer treatment for patients who cannot undergo surgery. During this procedure, a single-slot coaxial antenna is employed to effectively deliver microwave energy to the targeted tissue. The success of the treatment was measured by the amount of ablation zone created during the ablation procedure. The significantly large blood vessel placed near the antenna causes heat dissipation by convection around the blood vessel. The heat sink effect could result in insufficient ablation, raising the risk of local tumor recurrence. In this study, we investigated the heat loss due to large blood vessels and the relationship between blood velocity and temperature distribution. The hepatic artery, with a diameter of 4 mm and a height of 50 mm and two branches, is considered in the computational domain. The temperature profile, localized tissue contraction, and ablation zones were simulated for initial blood velocities 0.05, 0.1, and 0.16 m/s using the 3D Pennes bio-heat equation, temperature–time dependent model, and cell death model, respectively. Temperature-dependent blood velocity is modeled using the Navier–Stokes equation, and the fluid–solid interaction boundary is treated as a convective boundary. For discretization, we utilized $��H�\left(\mathrm{curl},\Omega \right)�$ elements for the wave propagation model, $��H^1�\left(\Omega \right)�$ elements for the Pennes bio-heat model, and $��{\left(�H^1�\left(\Omega \right)�\right)}^3�\times �L_0^2�\left(\Omega \right)�$ elements for the Navier–Stokes equation, where $��\Omega �$ represents the computational domain. The simulated results show that blood vessels and blood velocity have a significant impact on temperature distribution, tissue contraction, and the volume of the ablation zone.

We present a 3D discrete-continuum model to simulate blood pressure in large microvascular tissues in the absence of known capillary network architecture. Our hybrid approach combines a 1D Poiseuille flow description for large, discrete arteriolar and venular networks coupled to a continuum-based Darcy model, point sources of flux, for transport in the capillary bed. We evaluate our hybrid approach using a vascular network imaged from the mouse brain medulla/pons using multi-fluorescence high-resolution episcopic microscopy (MF-HREM). We use the fully-resolved vascular network to predict the hydraulic conductivity of the capillary network and generate a fully-discrete pressure solution to benchmark against. Our results demonstrate that the discrete-continuum methodology is a computationally feasible and effective tool for predicting blood pressure in real-world microvascular tissues when capillary microvessels are poorly defined.

Strokes are one of the leading causes of death in the United States. Stroke treatment involves removal or dissolution of the obstruction (usually a clot) in the blocked artery by catheter insertion. A computer simulation to systematically plan such patient-specific treatments needs a network of about 10⁵ blood vessels including collaterals. The existing computational fluid dynamic (CFD) solvers are not employed for stroke treatment planning as they are incapable of providing solutions for such big arterial trees in a reasonable amount of time. This work presents a novel one-dimensional mathematical formulation for blood flow modeling in an elastic blood vessel with a centrally placed rigid catheter. The governing equations are first-order hyperbolic partial differential equations, and the hypergeometric function needs to be computed to obtain the characteristic system of these hyperbolic equations. We employed the Discontinuous Galerkin method to solve the hyperbolic system and validated the implementation by comparing it against a well-established 3D CFD solver using idealized vessels and a realistic truncated arterial network. The results showed clinically insignificant differences in steady flow cases, with overall variations between 1D and 3D models remaining below 10%. Additionally, the solver accurately captured wave reflection phenomena at domain discontinuities in unsteady cases. A primary advantage of this model over 3D solvers is its ease in obtaining a discretized geometry of complex vasculatures with multiple arterial branches. Thus, the 1D computational model offers good accuracy and applicability in simulating complex vasculatures, demonstrating promising potential for investigating patient-specific endovascular interventions in strokes.

Thyroid nodules are commonly diagnosed with ultrasonography, which includes internal characteristics, varying looks, and hazy boundaries, making it challenging for a clinician to differentiate between malignant and benign forms based only on visual identification. The advancement of AI, particularly DL, provides significant breakthroughs in the domain of medical image identification. Yet, there are certain obstacles to achieving accuracy as well as efficacy in thyroid nodule detection. The thyroid nodules in this study are detected and classified using an inventive hybrid deep learning-assisted multi-classification method. The median blur method is applied in this work to eliminate the salt and pepper noise from the image. Then MPIU-Net-based segmentation is utilized to segment the image. The LGBPNP-based features are retrieved from the segmented image to obtain a single histogram sequence of the LGBP pattern in addition to other features like extraction of multi-texton and LTP-based features. After the feature extraction, the data augmentation process is applied and then the features are fed to the hybrid classification-based nodule classification model that comprises Deep Maxout and CNN, this hybrid classification trains the features and predicts the thyroid nodule. Additionally, the TIRADS score classification is used for the projected malignant thyroid nodule coupled with statistical features collected from the segmented. The DBNAAF with transfer learning model is employed to classify the grading of malignant thyroid nodules, where the weights of the model are learned with transfer learning. The MCC of the Hybrid Model is 0.9445, whereas the DCNN is 0.6858, YOLOV3-DMRF is 0.7229, CNN is 0.7780, DBN is 0.7601, Bi-GRU is 0.7038, Deep Maxout is 0.7528, and RNN is 0.8522, respectively.