Using Discrete Multiphysics Modelling to Assess the Effect of Calcification on Hemodynamic and Mechanical Deformation of Aortic Valve

Development of a digital twin of a tablet that mimics a real solid dosage form: Differences in the dissolution profile in conventional mini-USP II and a biorelevant colon model

The performance of colon-targeted solid dosage forms is commonly assessed using standardised pharmacopeial dissolution apparatuses like the USP II or the miniaturised replica, the mini-USP II. However, these fail to replicate the hydrodynamics and shear stresses in the colonic environment, which is crucial for the tablet's drug release process. In this work, computer simulations are used to create a digital twin of a dissolution apparatus and to develop a method to create a digital twin of a tablet that behaves realistically. These models are used to investigate the drug release profiles and shear rates acting on a tablet at different paddle speeds in the mini-USP II and biorelevant colon models to understand how the mini-USP II can be operated to achieve more realistic (i.e., in vivo) hydrodynamic conditions. The behaviour of the tablet and the motility patterns used in the simulations are derived from experimental and in vivo data, respectively, to obtain profound insights into the tablet's disintegration/drug release processes. We recommend an "on-off" operating mode in the mini-USP II to generate shear rate peaks, which would better reflect the in vivo conditions of the human colon instead of constant paddle speed.

Fluid-structure interaction simulation of calcified aortic valve stenosis

Calcified aortic valve stenosis (CAVS) is caused by calcium buildup and tissue thickening that impede the blood flow from left ventricle (LV) to aorta. In recent years, CAVS has become one of the most common cardiovascular diseases. Therefore, it is necessary to study the mechanics of aortic valve (AV) caused by calcification. In this paper, based on a previous idealized AV model, the hybrid immersed boundary/finite element method (IB/FE) is used to study AV dynamics and hemodynamic performance under normal and calcified conditions. The computational CAVS model is realized by dividing the AV leaflets into a calcified region and a healthy region, and each is described by a specific constitutive equation. Our results show that calcification can significantly affect AV dynamics. For example, the elasticity and mobility of the leaflets decrease due to calcification, leading to a smaller opening area with a high forward jet flow across the valve. The calcified valve also experiences an increase in local stress and strain. The increased loading due to AV stenosis further leads to a significant increase in left ventricular energy loss and transvalvular pressure gradients. The model predicted hemodynamic parameters are in general consistent with the risk classification of AV stenosis in the clinic. Therefore, mathematical models of AV with calcification have the potential to deepen our understanding of AV stenosis-induced ventricular dysfunction and facilitate the development of computational engineering-assisted medical diagnosis in AV related diseases.

Simulating the Hydrodynamic Conditions of the Human Ascending Colon: A Digital Twin of the Dynamic Colon Model

The performance of solid oral dosage forms targeting the colon is typically evaluated using standardised pharmacopeial dissolution apparatuses. However, these fail to replicate colonic hydrodynamics. This study develops a digital twin of the Dynamic Colon Model; a physiologically representative in vitro model of the human proximal colon. Magnetic resonance imaging of the Dynamic Colon Model verified that the digital twin robustly replicated flow patterns under different physiological conditions (media viscosity, volume, and peristaltic wave speed). During local contractile activity, antegrade flows of 0.06-0.78 cm s^-1 and backflows of -2.16--0.21 cm s^-1 were measured. Mean wall shear rates were strongly time and viscosity dependent although peaks were measured between 3.05-10.12 s^-1 and 5.11-20.34 s^-1 in the Dynamic Colon Model and its digital twin respectively, comparable to previous estimates of the USPII with paddle speeds of 25 and 50 rpm. It is recommended that viscosity and shear rates are considered when designing future dissolution test methodologies for colon-targeted formulations. In the USPII, paddle speeds >50 rpm may not recreate physiologically relevant shear rates. These findings demonstrate how the combination of biorelevant in vitro and in silico models can provide new insights for dissolution testing beyond established pharmacopeial methods.

A Simplified Framework for Modelling Viscoelastic Fluids in Discrete Multiphysics

A simplified modelling technique for modelling viscoelastic fluids is proposed from the perspective of Discrete Multiphysics. This approach, based on the concept of linear additive composition of energy potentials, aims to integrate Smooth Particle Hydrodynamics (SPH) with an equivalent elastic potential tailored for fluid flow simulations. The model was implemented using a particle-based software, explored and thoroughly validated with results from numerical experiments on three different flow conditions. The model was able to successfully capture a large extent of viscoelastic responses to external forcing, ranging from pure viscous flows to creep-dominated Bingham type of behaviour. It is concluded that, thanks to the modularity and tunable characteristics of the parameters involved, the proposed modelling approach can be a powerful simulation tool for modelling or mimicking the behaviour of viscoelastic substances.

Fluid-Structure Interaction in Coronary Stents: A Discrete Multiphysics Approach

Stenting is a common method for treating atherosclerosis. A metal or polymer stent is deployed to open the stenosed artery or vein. After the stent is deployed, the blood flow dynamics influence the mechanics by compressing and expanding the structure. If the stent does not respond properly to the resulting stress, vascular wall injury or re-stenosis can occur. In this work, a Discrete Multiphysics modelling approach is used to study the mechanical deformation of the coronary stent and its relationship with the blood flow dynamics. The major parameters responsible for deforming the stent are sorted in terms of dimensionless numbers and a relationship between the elastic forces in the stent and pressure forces in the fluid is established. The blood flow and the stiffness of the stent material contribute significantly to the stent deformation and affect its rate of deformation. The stress distribution in the stent is not uniform with the higher stresses occurring at the nodes of the structure. From the relationship (correlation) between the elastic force and the pressure force, depending on the type of material used for the stent, the model can be used to predict whether the stent is at risk of fracture or not after deployment.

Modelling Particle Agglomeration on through Elastic Valves under Flow

This work proposes a model of particle agglomeration in elastic valves replicating the geometry and the fluid dynamics of a venous valve. The fluid dynamics is simulated with Smooth Particle Hydrodynamics, the elastic leaflets of the valve with the Lattice Spring Model, while agglomeration is modelled with a 4-2 Lennard-Jones potential. All the models are combined together within a single Discrete Multiphysics framework. The results show that particle agglomeration occurs near the leaflets, supporting the hypothesis, proposed in previous experimental work, that clot formation in deep venous thrombosis is driven by the fluid dynamics in the valve.

Numerical Simulations of Red-Blood Cells in Fluid Flow: A Discrete Multiphysics Study

In this paper, we present a methodological study of modelling red blood cells (RBCs) in shear-induced flows based on the discrete multiphysics (DMP) approach. The DMP is an alternative approach from traditional multiphysics based on meshless particle-based methods. The proposed technique has been successful in modelling multiphysics and multi-phase problems with large interfacial deformations such as those in biological systems. In this study, we present the proposed method and introduce an accurate geometrical representation of the RBC. The results were validated against available data in the literature. We further illustrate that the proposed method is capable of modelling the rupture of the RBC membrane with minimum computational difficulty.

How to Modify LAMMPS: From the Prospective of a Particle Method Researcher

LAMMPS is a powerful simulator originally developed for molecular dynamics that, today, also accounts for other particle-based algorithms such as DEM, SPH, or Peridynamics. The versatility of this software is further enhanced by the fact that it is open-source and modifiable by users. This property suits particularly well Discrete Multiphysics and hybrid models that combine multiple particle methods in the same simulation. Modifying LAMMPS can be challenging for researchers with little coding experience. The available material explaining how to modify LAMMPS is either too basic or too advanced for the average researcher. In this work, we provide several examples, with increasing level of complexity, suitable for researchers and practitioners in physics and engineering, who are familiar with coding without been experts. For each feature, step by step instructions for implementing them in LAMMPS are shown to allow researchers to easily follow the procedure and compile a new version of the code. The aim is to fill a gap in the literature with particular reference to the scientific community that uses particle methods for (discrete) multiphysics.

Modelling and Simulation of the Drug Release from a Solid Dosage Form in the Human Ascending Colon: The Influence of Different Motility Patterns and Fluid Viscosities

For colonic drug delivery, the ascending part of the colon is the most favourable site as it offers the most suitable environmental conditions for drug dissolution. Commonly, the performance of a drug formulation is assessed using standardised dissolution apparatus, which does not replicate the hydrodynamics and shear stress evoked by wall motion in the colon. In this work, computer simulations are used to analyse and understand the influence of different biorelevant motility patterns on the disintegration/drug release of a solid dosage form (tablet) under different fluid conditions (viscosities) to mimic the ascending colonic environment. Furthermore, the ability of the motility pattern to distribute the drug in the ascending colon luminal environment is analysed to provide data for a spatiotemporal concentration profile. The motility patterns used are derived from in vivo data representing different motility patterns in the human ascending colon. The applied motility patterns show considerable differences in the drug release rate from the tablet, as well as in the ability to distribute the drug along the colon. The drug dissolution/disintegration process from a solid dosage form is primarily influenced by the hydrodynamic and shear stress it experiences, i.e., a combination of motility pattern and fluid viscosity. Reduced fluid motion leads to a more pronounced influence of diffusion in the tablet dissolution process. The motility pattern that provoked frequent single shear stress peaks seemed to be more effective in achieving a higher drug release rate. The ability to simulate drug release profiles under biorelevant colonic environmental conditions provides valuable feedback to better understand the drug formulation and how this can be optimised to ensure that the drug is present in the desired concentration within the ascending colon.