International Journal for Numerical Methods in Biomedical Engineering

Endoprosthetic hip replacement is the conventional way to treat osteoarthritis or a fracture of a dysfunctional joint. Different manufacturing methods are employed to create reliable patient-specific devices with long-term performance and biocompatibility. Recently, additive manufacturing has become a promising technique for the fabrication of medical devices, because it allows to produce complex samples with various structures of pores. Moreover, the limitations of traditional fabrication methods can be avoided...

Pulmonary hypertension is a cardiovascular disorder manifested by elevated mean arterial blood pressure (>20 mmHg) together with vessel wall stiffening and thickening due to alterations in collagen, elastin, and smooth muscle cells. Hypoxia-induced (type 3) pulmonary hypertension can be studied in animals exposed to a low oxygen environment for prolonged time periods leading to biomechanical alterations in vessel wall structure. This study introduces a novel approach to formulating a reduced order nonlinear elastic structural wall model for a large pulmonary artery...

Peristalsis is a common motion in various biological systems, especially the upper urinary tract, where it plays a pivotal role in conveying urine from the kidneys to the bladder. Using computational fluid dynamics, this study aims to investigate the effect of various peristaltic parameters on the motion of an obstacle through a two-dimensional ureter. Methodologically, Incompressible Navier-Stokes equations were utilized as the fluid domain's governing equations, and the Dynamic Mesh method (DM) was employed to simulate the peristaltic and obstacle motion...

Many transtibial amputees rate the fit between their residual limb and prosthetic socket as the most critical factor in satisfaction with using their prosthesis. This study aims to address the issue of prosthetic socket fit by reconfiguring the socket shape at the interface of the residual limb and socket. The proposed reconfigurable socket shifts pressure from sensitive areas and compensates for residual limb volume fluctuations, the most important factors in determining a good socket fit. Computed tomography scan images are employed to create the phantom limb of an amputee and to manufacture the reconfigurable socket...

To assess the health impact of inhaled aerosols, it is necessary to understand aerosol dynamics and the associated dosimetry in the human respiratory tract. Although several studies have measured or simulated the dosimetry of aerosol constituents, the respiratory tract focus areas have been limited. In particular, the aerosols generated from tobacco products are complex composites and simulating their dynamics in the respiratory tract is challenging. To assess the dosimetry of the aerosol constituents of tobacco products, we developed a revised version of the Multiple-Path Particle Dosimetry (MPPD) model, which employs (1) new geometry based on CT-scanned human respiratory tract data, (2) convective mixing in the oral cavity and deep lung, and (3) constituent partitioning between the tissue and air, and clearance...

In patients with type 2 diabetes (T2D), accurate prediction of hypoglycemic events is crucial for maintaining glycemic control and reducing their frequency. However, individuals with high blood glucose variability experience significant fluctuations over time, posing a challenge for early warning models that rely on static features. This article proposes a novel hypoglycemia early alarm framework based on dynamic feature selection. The framework incorporates domain knowledge and introduces multi-scale blood glucose features, including predicted values, essential for early warnings...

In most variance-based sensitivity analysis (SA) approaches applied to biomechanical models, statistical independence of the model input is assumed. However, often the model inputs are correlated. This might alter the interpretation of the SA results, which may severely impact the guidance provided during model development and personalization. Potential reasons for the infrequent usage of SA techniques that account for input correlation are the associated high computational costs, especially for models with many parameters, and the fact that the input correlation structure is often unknown...

We present a novel computational model for the dynamics of alveolar recruitment/derecruitment (RD), which reproduces the underlying characteristics typically observed in injured lungs. The basic idea is a pressure- and time-dependent variation of the stress-free reference volume in reduced dimensional viscoelastic elements representing the acinar tissue. We choose a variable reference volume triggered by critical opening and closing pressures in a time-dependent manner from a straightforward mechanical point of view...

Valvular heart diseases (such as stenosis and regurgitation) are recognized as a rapidly growing cause of global deaths and major contributors to disability. The most effective treatment for these pathologies is the replacement of the natural valve with a prosthetic one. Our work considers an innovative design for prosthetic aortic valves that combines the reliability and durability of artificial valves with the flexibility of tissue valves. It consists of a rigid support and three polymer leaflets which can be cut from an extruded flat sheet, and is referred to hereafter as the Wheatley aortic valve (WAV)...

Bone scaffolding is a promising approach for the treatment of critical-size bone defects. Hydroxyapatite can be used to produce highly porous scaffolds as it mimics the mineralized part of bone tissue, but its intrinsic brittleness limits its usage. Among 3D printing techniques, vat photopolymerization allows for the best printing resolution for ceramic materials. In this study, we implemented a Computed micro-Tomography based Finite Element Model of a hydroxyapatite porous scaffold fabricated by vat photopolymerization...

The mechanics of bile flow in the biliary system plays an important role in studying bile stasis and gallstone formation. Bile duct stricture is an abnormal phenomenon that refers to the bile duct getting smaller or narrower. The main objective of this study is to study the influence of stricture on bile flow dynamics using numerical methods. We employed a numerical Computational Fluid Dynamics model of the bile flow within a strictured hepatic duct. We studied and compared the influence of stricture severity, stricture length, eccentricity, and bile flow property on the bile flow dynamics...

The work presented here provides the guidelines and results for designing and implementing a highly sensitive modified Van der Pol - Duffing oscillator with a trigonometric damping function (VTD). This VTD can exhibit periodic and quasi-chaotic behavior necessary for application in weak signal detection. Here, we present two proposals: (1) A method based on a quasi-chaotic intermittent array (ANLIOA), whose all VTD parameters are calculated and fine-tuned toward a critical state between chaotic and periodic state through a Lyapunov exponent procedure, and (2) A method based on a single oscillator in an adaptive stopping oscillation system (ANLSOS), where VTD is established within an oscillatory regime...

We present our continuous efforts from a modeling and numerical viewpoint to develop a powerful and flexible mathematical and computational framework called Ocular Mathematical Virtual Simulator (OMVS). The OMVS aims to solve problems arising in biomechanics and hemodynamics within the human eye. We discuss our contribution towards improving the reliability and reproducibility of computational studies by performing a thorough validation of the numerical predictions against experimental data. The OMVS proved capable of simulating complex multiphysics and multiscale scenarios motivated by the study of glaucoma...

The study was carried out to investigate the effect of the artery with different pulse values and stenosis rates on the pressure drop, the peristaltic pump outlet pressure, fractional flow reserve (FFR) and most importantly the amount of power consumed by the peristaltic pump. For this purpose, images taken from the clinical environment were produced as models (10 mm inlet diameter) with 0% and 70% percent areal stenosis rates (PSR) on a three-dimensional (3D) printer. In the experimental system, pure water was used as the fluid at 54, 84, 114, 132, and 168 bpm pulse values...

In silico trials are a promising way to increase the efficiency of the development, and the time to market of cardiovascular implantable devices. The development of transcatheter aortic valve implantation (TAVI) devices, could benefit from in silico trials to overcome frequently occurring complications such as paravalvular leakage and conduction problems. To be able to perform in silico TAVI trials virtual cohorts of TAVI patients are required. In a virtual cohort, individual patients are represented by computer models that usually require patient-specific aortic valve geometries...

Breast cancer is a significant public health issue affecting women worldwide. While advancements in treatment options have led to improved survival rates, the impact of breast cancer and its treatments on bone health cannot be overlooked. Bone remodeling is a complex process regulated by the delicate balance between bone formation and resorption. Any disruption to this balance can lead to decreased bone density, increased fracture risk, and compromised physical function. To investigate the effects of breast cancer and its treatments on bone remodeling, a finite element model was developed in this study...

We present a reduced order model for efficient nonlinear homogenization of bones, accounting for strength difference effects and containing some well-known plasticity models (like von Mises or Drucker-Prager) as special cases. The reduced order homogenization is done by using a cluster-based model order reduction technique, called cluster-based nonuniform transformation field analysis. For an offline phase, a space-time decomposition is performed on the mesoscopic plastic strain fields, while a clustering analysis is employed for a spatial decomposition of the mesoscale RVE model...

Reducing the computational time required by high-fidelity, full-order models (FOMs) for the solution of problems in cardiac mechanics is crucial to allow the translation of patient-specific simulations into clinical practice. Indeed, while FOMs, such as those based on the finite element method, provide valuable information on the cardiac mechanical function, accurate numerical results can be obtained at the price of very fine spatio-temporal discretizations. As a matter of fact, simulating even just a few heartbeats can require up to hours of wall time on high-performance computing architectures...

In this work, we develop numerical methods for the solution of blood flow and coagulation on dynamic adaptive moving meshes. We consider the blood flow as a flow of incompressible Newtonian fluid governed by the Navier-Stokes equations. The blood coagulation is introduced through the additional Darcy term, with a permeability coefficient dependent on reactions. To this end, we introduce moving mesh collocated finite-volume methods for the Navier-Stokes equations, advection-diffusion equations, and a method for the stiff cascade of reactions...

When modeling complex fluid networks using one-dimensional (1D) approaches, boundary conditions can be imposed using zero-dimensional (0D) models. An application case is the modeling of the entire human circulation using closed-loop models. These models can be considered as a tool to investigate short-term transient and stationary hemodynamic responses to postural changes. The first shortcoming of existing 1D modeling methods in simulating these sudden maneuvers is their inability to deal with rapid variations in flow conditions, as they are limited to the subsonic case...