In small capillary hemodynamics, the cell-free layer is a near-wall layer of plasma absent of red blood cells since they are subject to migration to the capillary center in Poiseuille flow[1] .Cell-free marginal layer model is a mathematical model which tries to explain Fåhræus–Lindqvist effect mathematically.

Mathematical modeling

Governing equations

Consider steady flow of blood through a capillary of radius \(R\). The capillary cross section can be divided into a core region and cell-free plasma region near the wall. The governing equations for both regions can be given by the following equations[2] :

\[ \frac{ -\Delta P}{ L }=\frac{1}{r}\frac{d}{dr}(\mu_c r \frac{du_c}{dr});\] \( 0 \le r\ \le R-\delta\,\) \[ \frac{ -\Delta P}{ L }=\frac{1}{r}\frac{d}{dr}(\mu_p r \frac{du_p}{dr});\] \( R-\delta\le r\ \le R\ \,\)

where:

\[\Delta P \] is the pressure drop across the capillary \[L\] is the length of capillary \[ u_c \] is velocity in core region \[ u_p \] is velocity of plasma in cell-free region \[ \mu_{c} \] is viscosity in core region \[ \mu_{p} \] is viscosity of plasma in cell-free region \[\delta\] is the cell-free plasma layer thickness

Boundary conditions

The boundary conditions to obtain the solution for the two differential equations presented above are that the velocity gradient is zero in the tube center, no slip occurs at the tube wall and the velocity and the shear stress are continuous at the interface between the two zones. These boundary conditions can be expressed mathematically as:

  • \(\left. \frac{du_c}{dr}\right|_{r= 0}=0\)
  • \(\left.u_p\right|_{r= R}=0\)
  • \(\left.u_p\right|_{r= R-\delta}=\left.u_c\right|_{r= R-\delta}\)
  • \(\left.\tau_p\right|_{r= R-\delta}=\left.\tau_c\right|_{r= R-\delta}\)

Velocity profiles

Integrating governing equations with respect to r and applying the above discussed boundary conditions will result in:

\[ u_c=\frac{ \Delta P R^2}{ 4\mu_c L }[1-(\frac{ R-\delta}{R})^2-\frac{\mu_p}{\mu_c}(\frac{r}{R})^2+\frac{\mu_p}{\mu_c}(\frac{ R-\delta}{R})^2]\]

\[ u_p=\frac{ \Delta P R^2}{ 4\mu_p L }[1-(\frac{r}{R})^2]\]

Volumetric flow rate

Total volumetric flow rate is the algebraic sum of the flow rates in core and plasma region. The expression for the total volumetric flow rate can be written as:

\[ Q=\frac{ \pi \Delta P R^4}{ 8\mu_p L }[1-(1-\frac{\delta}{R})^4(1-\frac{\mu_p}{\mu_c})]\]

Comparison with the viscosity which applies in the Poiseuille flow yields effective viscosity, \( \mu_{e} \) as:

\[ \mu_{e}=\frac{\mu_p}{[1-(1-\frac{\delta}{R})^4(1-\frac{\mu_p}{\mu_c})]} \]

It can be realized when the radius of the blood vessel is very larger than the thickness of the cell-free plasma layer, the effective viscosity is equal to bulk blood viscosity \( \mu_c \) at high shear rates (Newtonian fluid).

See also

References

  1. Script error
  2. Script error