ترغب بنشر مسار تعليمي؟ اضغط هنا

Spatio-temporal dynamics of dilute red blood cell suspensions in a microchannel flow at low Reynolds number

114   0   0.0 ( 0 )
 نشر من قبل Timm Kr\\\"uger
 تاريخ النشر 2019
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

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.



قيم البحث

اقرأ أيضاً

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 (R BCs) 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.
A new approach to turbulence simulation, based on a combination of large-eddy simulation (LES) for the whole flow and an array of non-space-filling quasi-direct numerical simulations (QDNS), which sample the response of near-wall turbulence to large- scale forcing, is proposed and evaluated. The technique overcomes some of the cost limitations of turbulence simulation, since the main flow is treated with a coarse-grid LES, with the equivalent of wall functions supplied by the near-wall sampled QDNS. Two cases are tested, at friction Reynolds number Re$_tau$=4200 and 20,000. The total grid node count for the first case is less than half a million and less than two million for the second case, with the calculations only requiring a desktop computer. A good agreement with published DNS is found at Re$_tau$=4200, both in terms of the mean velocity profile and the streamwise velocity fluctuation statistics, which correctly show a substantial increase in near-wall turbulence levels due to a modulation of near-wall streaks by large-scale structures. The trend continues at Re$_tau$=20,000, in agreement with experiment, which represents one of the major achievements of the new approach. A number of detailed aspects of the model, including numerical resolution, LES-QDNS coupling strategy and sub-grid model are explored. A low level of grid sensitivity is demonstrated for both the QDNS and LES aspects. Since the method does not assume a law of the wall, it can in principle be applied to flows that are out of equilibrium.
Soft hydraulics, which addresses the interaction between an internal flow and a compliant conduit, is a central problem in microfluidics. We analyze Newtonian fluid flow in a rectangular duct with a soft top wall at steady state. The resulting fluid- -structure interaction (FSI) is formulated for both vanishing and finite flow inertia. At the leading-order in the small aspect ratio, the lubrication approximation implies that the pressure only varies in the streamwise direction. Meanwhile, the compliant walls slenderness makes the fluid--solid interface behave like a Winkler foundation, with the displacement fully determined by the local pressure. Coupling flow and deformation and averaging across the cross-section leads to a one-dimensional reduced model. In the case of vanishing flow inertia, an effective deformed channel height is defined rigorously to eliminate the spanwise dependence of the deformation. It is shown that a previously-used averaged height concept is an acceptable approximation. From the one-dimensional model, a friction factor and the corresponding Poiseuille number are derived. Unlike the rigid duct case, the Poiseuille number for a compliant duct is not constant but varies in the streamwise direction. Compliance can increase the Poiseuille number by a factor of up to four. The model for finite flow inertia is obtained by assuming a parabolic vertical variation of the streamwise velocity. To satisfy the displacement constraints along the edges of the channel, weak tension is introduced in the streamwise direction to regularize the Winkler-foundation-like model. Matched asymptotic solutions of the regularized model are derived.
We present guidelines to estimate the effect of electrostatic repulsion in sedimenting dilute particle suspensions. Our results are based on combined Langevin dynamics and lattice Boltzmann simulations for a range of particle radii, Debye lengths and particle concentrations. They show a simple relationship between the slope $K$ of the sedimentation velocity over the concentration versus the range $chi$ of the electrostatic repulsion normalized by the average particle-particle distance. When $chi to 0$, the particles are too far away from each other to interact electrostatically and $K=6.55$ as predicted by the theory of Batchelor. As $chi$ increases, $K$ likewise increases up to a maximum around $chi=0.5$ and then decreases again to a concentration-dependent constant over the range $chi=0.5-1$, while the particles transition from a disordered gas-like distribution to a liquid-like state with a narrow distribution of the interparticle spacing.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا