Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userOn the effects of membrane viscosity on transient red blood cell dynamics
109
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Computational Fluid Dynamics (CFD) is currently used to design and improve the hydraulic properties of biomedical devices, wherein the large scale blood circulation needs to be simulated by accounting for the mechanical response of red blood cells (RBCs) at mesoscales. In many practical instances, biomedical devices work on time-scales comparable to the intrinsic relaxation time of RBCs: thus, a systematic understanding of the time-dependent response of erythrocyte membranes is crucial for the effective design of such devices. So far, this information has been deduced from experimental data, which do not necessarily adapt to the the broad variety of the fluid dynamic conditions that can be encountered in practice. This work explores the novel possibility of studying the time-dependent response of an erythrocyte membrane to external mechanical loads via ab-initio, mesoscale numerical simulations, with a primary focus on the detailed characterisation of the RBC relaxation time $t_c$ following the arrest of the external mechanical load. The adopted mesoscale model exploits a hybrid Immersed Boundary-Lattice Boltzmann Method (IB-LBM), coupled with the Standard Linear Solid model (SLS) to account for the RBC membrane viscosity. We underscore the key importance of the 2D membrane viscosity $mu_{m}$ to correctly reproduce the relaxation time of the RBC membrane. A detailed assessment of the dependencies on the typology and strength of the applied mechanical loads is also provided. Overall, our findings open interesting future perspectives for the study of the non-linear response of RBCs immersed in time-dependent strain fields.
rate research
Read More
Microfluidic technologies are commonly used for the manipulation of red blood cell (RBC) suspensions and analyses of flow-mediated biomechanics. To enhance the performance of microfluidic devices, understanding the dynamics of the suspensions processed within is crucial. We report novel aspects of the spatio-temporal dynamics of RBC suspensions flowing through a typical microchannel at low Reynolds number. Through experiments with dilute RBC suspensions, we find an off-centre two-peak (OCTP) profile of cells contrary to the centralised distribution commonly reported for low-inertia flows. This is reminiscent of the well-known tubular pinch effect which arises from inertial effects. However, given the conditions of negligible inertia in our experiments, an alternative explanation is needed for this OCTP profile. Our massively-parallel simulations of RBC flow in real-size microfluidic dimensions using the immersed-boundary-lattice-Boltzmann method (IB-LBM) confirm the experimental findings and elucidate the underlying mechanism for the counterintuitive RBC pattern. By analysing the RBC migration and cell-free layer (CFL) development within a high-aspect-ratio channel, we show that such a distribution is co-determined by the spatial decay of hydrodynamic lift and the global deficiency of cell dispersion in dilute suspensions. We find a CFL development length greater than 46 and 28 hydraulic diameters in the experiment and simulation, respectively, exceeding typical lengths of microfluidic designs. Our work highlights the key role of transient cell distribution in dilute suspensions, which may negatively affect the reliability of experimental results if not taken into account.
Flux of rigid or soft particles (such as drops, vesicles, red blood cells, etc.) in a channel is a complex function of particle concentration, which depends on the details of induced dissipation and suspension structure due to hydrodynamic interactions with walls or between neighboring particles. Through two-dimensional and three-dimensional simulations and a simple model that reveals the contribution of the main characteristics of the flowing suspension, we discuss the existence of an optimal volume fraction for cell transport and its dependence on the cell mechanical properties. The example of blood is explored in detail, by adopting the commonly used modeling of red blood cells dynamics. We highlight the complexity of optimization at the level of a network, due to the antagonist evolution of local volume fraction and optimal volume fraction with the channels diameter. In the case of the blood network, the most recent results on the size evolution of vessels along the circulatory network of healthy organs suggest that the red blood cell volume fraction (hematocrit) of healthy subjects is close to optimality, as far as transport only is concerned. However, the hematocrit value of patients suffering from diverse red blood cel pathologies may strongly deviate from optimality.
Thrombosis is a common complication following the surgical implantation of blood contacting devices, and is strongly influenced by the phenomenon of near-wall enrichment of platelets. This paper describes a multi-constituent continuum approach to study this phenomenon. A mixture-theory model is used to describe the motion of the plasma and the red blood cells (RBCs) and the interactions between the two components. A transport model is developed to study the influence of the RBC field on the platelets. The model is used to study blood flow in a rectangular micro-channel, a sudden expansion micro-channel, and a channel containing micro crevices (representing a practical problem encountered in most blood-wetted devices). The simulations show that in the rectangular channel the concentration of platelets near the walls is about five times higher than the concentration near the centerline of the channel. It is also noticed that in the channel with crevices, extremely a large number of platelets accumulate in the deep part of the crevices and this may serve as the nidus for thrombosis occurring in medical devices. Keywords: Platelets; Blood flow; RBCs; Micro-channel; Mixture theory; Two-fluids; Thrombosis
A hybrid computational method coupling the lattice-Boltzmann (LB) method and a Langevin-dynamics (LD) method is developed to simulate nanoscale particle and polymer (NPP) suspensions in the presence of both thermal fluctuation and long-range many-body hydrodynamic interactions (HI). Brownian motion of the NPP is explicitly captured by a stochastic forcing term in the LD method. The LD method is two-way coupled to the non-fluctuating LB fluid through a discrete LB forcing source distribution to capture the long-range HI. To ensure intrinsically linear scalability with respect to the number of particles, an Eulerian-host algorithm for short-distance particle neighbor search and interaction is developed and embedded to LB-LD framework. The validity and accuracy of the LB-LD approach are demonstrated through several sample problems. The simulation results show good agreements with theory and experiment. The LB-LD approach can be favorably incorporated into complex multiscale computational frameworks for efficiently simulating multiscale, multicomponent particulate suspension systems such as complex blood suspensions.
The effects of viscoelasticity on the dynamics and break-up of fluid threads in microfluidic T-junctions are investigated using numerical simulations of dilute polymer solutions at changing the Capillary number ($mbox {Ca}$), i.e. at changing the balance between the viscous forces and the surface tension at the interface, up to $mbox{Ca} approx 3 times 10^{-2}$. A Navier-Stokes (NS) description of the solvent based on the lattice Boltzmann models (LBM) is here coupled to constitutive equations for finite extensible non-linear elastic dumbbells with the closure proposed by Peterlin (FENE-P model). We present the results of three-dimensional simulations in a range of $mbox{Ca}$ which is broad enough to characterize all the three characteristic mechanisms of breakup in the confined T-junction, i.e. ${it squeezing}$, ${it dripping}$ and ${it jetting}$ regimes. The various model parameters of the FENE-P constitutive equations, including the polymer relaxation time $tau_P$ and the finite extensibility parameter $L^2$, are changed to provide quantitative details on how the dynamics and break-up properties are affected by viscoelasticity. We will analyze cases with ${it Droplet ~Viscoelasticity}$ (DV), where viscoelastic properties are confined in the dispersed (d) phase, as well as cases with ${it Matrix ~Viscoelasticity}$ (MV), where viscoelastic properties are confined in the continuous (c) phase. Moderate flow-rate ratios $Q approx {cal O}(1)$ of the two phases are considered in the present study. Overall, we find that the effects are more pronounced in the case with MV, as the flow driving the break-up process upstream of the emerging thread can be sensibly perturbed by the polymer stresses.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
