No Arabic abstract
The microvascular networks in the body of vertebrates consist of the smallest vessels such as arterioles, capillaries, and venules. The flow of RBCs through these networks ensures the gas exchange in as well as the transport of nutrients to the tissues. Any alterations in this blood flow may have severe implications on the health state. Since the vessels in these networks obey dimensions similar to the diameter of RBCs, dynamic effects on the cellular scale play a key role. The steady progression in the numerical modeling of RBCs, even in complex networks, has led to novel findings in the field of hemodynamics, especially concerning the impact and the dynamics of lingering events, when a cell meets a branch of the network. However, these results are yet to be matched by a detailed analysis of the lingering experiments in vivo. To quantify this lingering effect in in vivo experiments, this study analyzes branching vessels in the microvasculature of Syrian golden hamsters via intravital microscopy and the use of an implanted dorsal skinfold chamber. It also presents a detailed analysis of these lingering effects of cells at the apex of bifurcating vessels, affecting the temporal distribution of cell-free areas of blood flow in the branches, even causing a partial blockage in severe cases.
Using a multiscale blood flow solver, the complete diffusion tensor of nanoparticle (NP) in sheared cellular blood flow is calculated over a wide range of shear rate and haematocrit. In the short-time regime, NPs exhibit anomalous dispersive behaviors under high shear and high haematocrit due to the transient elongation and alignment of the red blood cells (RBCs). In the long-time regime, the NP diffusion tensor features high anisotropy. Particularly, there exists a critical shear rate ($sim$100 $s^{-1}$) around which the shear-rate dependence of the diffusivity tensor changes from linear to nonlinear scale. Above the critical shear rate, the cross-stream diffusivity terms vary sublinearly with shear rate, while the longitudinal term varies superlinearly. The dependence on haematocrit is linear in general except at high shear rates, where a sublinear scale is found for the vorticity term and a quadratic scale for the longitudinal term. Through analysis of the suspension microstructure and numerical experiments, the nonlinear hemorheological dependence of the NP diffusion tensor is attributed to the streamwise elongation and cross-stream contraction of RBCs under high shear, quantified by a Capillary number. The RBC size is shown to be the characteristic length scale affecting the RBC-enhanced shear-induced diffusion (RESID), while the NP size at submicron exhibits negligible influence on the RESID. Based on the observed scaling behaviors, empirical correlations are proposed to bridge the NP diffusion tensor to specific shear rate and haematocrit. The characterized NP diffusion tensor provides a constitutive relation that can lead to more effective continuum models to tackle large-scale NP biotransport applications.
Various biological processes such as transport of oxygen and nutrients, thrombus formation, vascular angiogenesis and remodeling are related to cellular/subcellular level biological processes, where mesoscopic simulations resolving detailed cell dynamics provide a key to understanding and identifying the cellular basis of disease. To break this bottleneck and achieve a biologically meaningful timescale, we propose a multiscale parareal algorithm in which a continuum-based solver supervises a mesoscopic simulation in the time-domain. Using an iterative prediction-correction strategy, the parallel-in-time mesoscopic simulation supervised by its continuum-based counterpart can converge fast. The effectiveness of the proposed method is first verified in a time-dependent flow with a sinusoidal flowrate through a Y-shaped bifurcation channel. Physical quantities of interest including velocity, wall shear stress and flowrate are computed to compare against those of reference solutions, showing a less than 1% relative error on flowrate in the Newtonian flow and a less than 3% relative error in the non-Newtonian blood flow. The proposed method is then applied to a large-scale mesoscopic simulation of microvessel blood flow in a zebrafish hindbrain for temporal acceleration. The time-dependent blood flow from heartbeats in this realistic vascular network of zebrafish hindbrain is simulated using dissipative particle dynamics as the mesoscopic model, which is supervised by a one-dimensional blood flow model (continuum-based model) in multiple temporal sub-domains. The computational analysis shows that the resulting microvessel blood flow converges to the reference solution after only two iterations. The proposed method is suitable for long-time mesoscopic simulations with complex fluids and geometries.
Blood flowing through microvascular bifurcations has been an active research topic for many decades, while the partitioning pattern of nanoscale solutes in the blood remains relatively unexplored. Here, we demonstrate a multiscale computational framework for direct numerical simulation of the nanoparticle (NP) partitioning through physiologically-relevant vascular bifurcations in the presence of red blood cells (RBCs). The computational framework is established by embedding a newly-developed particulate suspension inflow/outflow boundary condition into a multiscale blood flow solver. The computational framework is verified by recovering a tubular blood flow without a bifurcation and validated against the experimental measurement of an intravital bifurcation flow. The classic Zweifach-Fung (ZF) effect is shown to be well captured by the method. Moreover, we observe that NPs exhibit a ZF-like heterogeneous partition in response to the heterogeneous partition of the RBC phase. The NP partitioning prioritizes the high-flow-rate daughter branch except for extreme (large or small) suspension flow partition ratios under which the complete phase separation tends to occur. By analyzing the flow field and the particle trajectories, we show that the ZF-like heterogeneity in NP partition can be explained by the RBC-entrainment effect caused by the deviation of the flow separatrix preceded by the tank-treading of RBCs near the bifurcation junction. The recovery of homogeneity in the NP partition under extreme flow partition ratios is due to the plasma skimming of NPs in the cell-free layer. These findings, based on the multiscale computational framework, provide biophysical insights to the heterogeneous distribution of NPs in microvascular beds that are observed pathophysiologically.
Bubbles introduced to the arterial circulation during invasive medical procedures can have devastating consequences for brain function but their effects are currently difficult to quantify. Here we present a Monte-Carlo simulation investigating the impact of gas bubbles on cerebral blood flow. For the first time, this model includes realistic adhesion forces, bubble deformation, fluid dynamical considerations, and bubble dissolution. This allows investigation of the effects of buoyancy, solubility, and blood pressure on embolus clearance. Our results illustrate that blockages depend on several factors, including the number and size distribution of incident emboli, dissolution time and blood pressure. We found it essential to model the deformation of bubbles to avoid overestimation of arterial obstruction. Incorporation of buoyancy effects within our model slightly reduced the overall level of obstruction but did not decrease embolus clearance times. We found that higher blood pressures generate lower levels of obstruction and improve embolus clearance. Finally, we demonstrate the effects of gas solubility and discuss potential clinical applications of the model.
We present experiments on RBCs that flow through microcapillaries under physiological conditions. We show that the RBC clusters form as a subtle imbrication between hydrodynamics interaction and adhesion forces because of plasma proteins. Clusters form along the capillaries and macromolecule-induced adhesion contribute to their stability. However, at high yet physiological flow velocities, shear stresses overcome part of the adhesion forces, and cluster stabilization due to hydrodynamics becomes stronger. For the case of pure hydrodynamic interaction, cell-to-cell distances have a pronounced bimodal distribution. Our 2D-numerical simulations on vesicles captures the transition between adhesive and non-adhesive clusters at different flow velocities.