Multiscale Modeling of Blood Flow: Coupling Finite Elements with Smoothed Dissipative Particle Dynamics

Computational mesoscale framework for biological clustering and fractal aggregation

Hierarchical clustering due to diffusion and reaction is a widespread occurrence in natural phenomena, displaying fractal behavior with non-integer size scaling. The study of this phenomenon has garnered interest in both biological systems such as morphogenesis and blood clotting, and synthetic systems such as colloids and polymers. The modeling of biological clustering can be difficult, as it can occur on a variety of scales and involve multiple mechanisms, necessitating the use of various methods to capture its behavior. Here, we propose a novel framework, the generalized-mesoscale-clustering (GMC), for the study of complex hierarchical clustering phenomena in biological systems. The GMC framework incorporates the effects of hydrodynamic interactions, bonding, and surface tension, and allows for the analysis of both static and dynamic states of cluster development. The framework is applied to a range of biological clustering mechanisms, with a focus on blood-related clustering from fibrin network formation to platelet aggregation. Our study highlights the importance of a comprehensive characterization of the structural properties of the cluster, including fractal dimension, pore-scale diffusion, initiation time, and consolidation time, in fully understanding the behavior of biological clustering systems. The GMC framework also provides the potential to investigate the temporal evolution and mechanical properties of the clusters by tracking bond density and including hydrodynamic interactions.

An arbitrary Lagrangian Eulerian smoothed particle hydrodynamics (ALE-SPH) method with a boundary volume fraction formulation for fluid-structure interaction

We present a new weakly-compressible smoothed particle hydrodynamics (SPH) method capable of modeling non-slip fixed and moving wall boundary conditions. The formulation combines a boundary volume fraction (BVF) wall approach with the transport-velocity SPH method. The resulting method, named SPH-BVF, offers detection of arbitrarily shaped solid walls on-the-fly, with small computational overhead due to its local formulation. This simple framework is capable of solving problems that are difficult or infeasible for standard SPH, namely flows subject to large shear stresses or at moderate Reynolds numbers, and mass transfer in deformable boundaries. In addition, the method extends the transport-velocity formulation to reaction-diffusion transport of mass in Newtonian fluids and linear elastic solids, which is common in biological structures. Taken together, the SPH-BVF method provides a good balance of simplicity and versatility, while avoiding some of the standard obstacles associated with SPH: particle penetration at the boundaries, tension instabilities and anisotropic particle alignments, that hamper SPH from being applied to complex problems such as fluid-structure interaction in a biological system.

Development of a novel hybrid method combining finite difference method and dissipative particle dynamics to simulate thrombus formation on orifice flow

In our previous works, the transport of activated platelets (APs) on orifice flow has been simulated by finite difference method (FDM). And the distribution of AP concentration on the flow was obtained. However, the effect of platelet aggregation on the distribution of AP concentration can't be investigated by FDM because FDM can't simulate platelet aggregation. On the other hand, platelet aggregation has been simulated by dissipative particle dynamics (DPD). In this paper, a hybrid method combining FDM and DPD is proposed to investigate the effect of platelet aggregation on the distribution of AP concentration. And the hybrid method is used to simulate thrombus formation on orifice flow. As for the effect of platelet aggregation, it is found that the distribution of AP concentration in the hybrid method is different from the distribution in FDM at the places of platelet aggregation. It is considered that the difference is induced by platelet aggregation. As for the distribution of thrombus, higher AP concentration and more aggregated APs are found around the reattachment point and in the recirculation area. It is considered that thrombus is mainly distributed at these places in the simulation. And according to our previous experimental results, thrombus is mainly distributed around the reattachment point and in the recirculation area. It is concluded that the effect of platelet aggregation on the distribution of AP concentration can be investigated by the hybrid method, and the computational results agree with our previous experimental results.

Recent advances in computational modeling of fibrin clot formation: A review

The field of thrombosis and hemostasis is crucial for understanding and developing new therapies for pathologies such as deep vein thrombosis, diabetes related strokes, pulmonary embolisms, and hemorrhaging related diseases. In the last two decades, an exponential growth in studies related to fibrin clot formation using computational tools has been observed. Despite this growth, the complete mechanism behind thrombus formation and hemostasis has been long and rife with obstacles; however, significant progress has been made in the present century. The computational models and methods used in this context are diversified into different spatiotemporal scales, yet there is no single model which can predict both physiological and mechanical properties of fibrin clots. In this review, we list the major strategies employed by researchers in modeling fibrin clot formation using recent and existing computational techniques. This review organizes the computational strategies into continuum level, system level, discrete particle (DPD), and multi-scale methods. We also discuss strengths and weaknesses of various methods and future directions in which computational modeling of fibrin clots can advance.

Multiscale simulation of ideal mixtures using smoothed dissipative particle dynamics

Smoothed dissipative particle dynamics (SDPD) [P. Español and M. Revenga, Phys. Rev. E 67, 026705 (2003)] is a thermodynamically consistent particle-based continuum hydrodynamics solver that features scale-dependent thermal fluctuations. We obtain a new formulation of this stochastic method for ideal two-component mixtures through a discretization of the advection-diffusion equation with thermal noise in the concentration field. The resulting multicomponent approach is consistent with the interpretation of the SDPD particles as moving volumes of fluid and reproduces the correct fluctuations and diffusion dynamics. Subsequently, we provide a general multiscale multicomponent SDPD framework for simulations of molecularly miscible systems spanning length scales from nanometers to the non-fluctuating continuum limit. This approach reproduces appropriate equilibrium properties and is validated with simulation of simple one-dimensional diffusion across multiple length scales.

Size consistency in smoothed dissipative particle dynamics

Smoothed dissipative particle dynamics (SDPD) is a mesoscopic method that allows one to select the level of resolution at which a fluid is simulated. In this work, we study the consistency of the resulting thermodynamic properties as a function of the size of the mesoparticles, both at equilibrium and out of equilibrium. We also propose a reformulation of the SDPD equations in terms of energy variables. This increases the similarities with dissipative particle dynamics with energy conservation and opens the way for a coupling between the two methods. Finally, we present a numerical scheme for SDPD that ensures the conservation of the invariants of the dynamics. Numerical simulations illustrate this approach.

Hybrid molecular-continuum simulations using smoothed dissipative particle dynamics

We present a new multiscale simulation methodology for coupling a region with atomistic detail simulated via molecular dynamics (MD) to a numerical solution of the fluctuating Navier-Stokes equations obtained from smoothed dissipative particle dynamics (SDPD). In this approach, chemical potential gradients emerge due to differences in resolution within the total system and are reduced by introducing a pairwise thermodynamic force inside the buffer region between the two domains where particles change from MD to SDPD types. When combined with a multi-resolution SDPD approach, such as the one proposed by Kulkarni et al. [J. Chem. Phys. 138, 234105 (2013)], this method makes it possible to systematically couple atomistic models to arbitrarily coarse continuum domains modeled as SDPD fluids with varying resolution. We test this technique by showing that it correctly reproduces thermodynamic properties across the entire simulation domain for a simple Lennard-Jones fluid. Furthermore, we demonstrate that this approach is also suitable for non-equilibrium problems by applying it to simulations of the start up of shear flow. The robustness of the method is illustrated with two different flow scenarios in which shear forces act in directions parallel and perpendicular to the interface separating the continuum and atomistic domains. In both cases, we obtain the correct transient velocity profile. We also perform a triple-scale shear flow simulation where we include two SDPD regions with different resolutions in addition to a MD domain, illustrating the feasibility of a three-scale coupling.

A clinically oriented introduction and review on finite element models of the human cochlea

Due to the inaccessibility of the inner ear, direct in vivo information on cochlear mechanics is difficult to obtain. Mathematical modelling is a promising way to provide insight into the physiology and pathology of the cochlea. Finite element method (FEM) is one of the most popular discrete mathematical modelling techniques, mainly used in engineering that has been increasingly used to model the cochlea and its elements. The aim of this overview is to provide a brief introduction to the use of FEM in modelling and predicting the behavior of the cochlea in normal and pathological conditions. It will focus on methodological issues, modelling assumptions, simulation of clinical scenarios, and pathologies.