Next Article in Journal
When Four Cyclic Antipodal Pairs Are Ordered Counterclockwise in Euclidean and Hyperbolic Geometry
Previous Article in Journal
Feasibility Study on Active Structural Attenuation: Addressing Multiband Vibration in Automotive Vehicles on 2D Asymmetric Structures with a Faulty Horizontal Actuator
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mathematical Analysis and Simulation of the F-L Effect in Two-Layered Blood Flow through the Capillaries Remote from the Heart and Proximate to Human Tissue

1
Department of Physical Sciences, Mahatma Gandhi Chitrakoot Gramodaya Vishwavidyalaya Chitrakoot, Satna 485334, India
2
Department of Biological Sciences, Mahatma Gandhi Chitrakoot Gramodaya Vishwavidyalaya Chitrakoot, Satna 485334, India
3
Department of Mathematics, Dehradun Institute of Technology University, Dehradun 248009, India
4
Department of Mathematics, Government Polytechnic, Nawada 805114, India
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(6), 728; https://doi.org/10.3390/sym16060728
Submission received: 14 April 2024 / Revised: 26 May 2024 / Accepted: 5 June 2024 / Published: 11 June 2024
(This article belongs to the Section Mathematics)

Abstract

:
In this paper, we have provided a mathematical analysis of an empirical result, namely, the Fahraeus–Lindqvist effect, a phenomenon that occurs in capillary tubes with a diameter lower than 0.3 mm. According to this effect, in capillary tubes under 0.3 mm in diameter, the apparent viscosity of blood decreases as the diameter of the tube decreases, making flow possible in these vessels. A two-phase blood flow mathematical model for human capillaries has been presented here. According to Haynes’ theory, blood is separated into two layers when it flows from the capillary. It is assumed that the first layer is plasma, and the second layer is the core layer. The plasma layer flows near the wall of the capillary, and the core layer flows along the axis of the capillary. Further, the core layer is assumed to be a mixture of two phases: one is the plasma, and the other is that of RBCs. For mathematical modeling purposes, a curvilinear coordinate system has been adopted, with physical quantities used in tensorial form. Derived equations are solved to find the effective viscosity, which depends upon the radius of the capillary; that is, it reduces viscosity to make blood flow possible. A comparative study was conducted with the experimental result of this effect, and it was observed that the proposed two-phase blood flow model is much closer to the experimental data than the single-phase blood flow model, and both have the same trends. After validation of the model with the experimental result, this model was applied to human capillaries (diameter lower than 10 μ m ) to show the F-L effect, and the impact of various physiological quantities that are relevant to the flow of blood into human capillaries are also discussed here. The impact of hematocrit on various parameters has been demonstrated explicitly.

1. Introduction

Capillaries are the smallest vessel in the human body [1]. They are tiny branches of the arterioles. Their internal diameter is about 4 to 10 μ m , less than 1/3 of a hair’s width [2,3]. They are composed of a single layer of flat cells. This layer not only forms the walls of the capillaries but also lines the other blood vessels and the heart, so all the blood in the body is contained within a single continuous envelope. Capillaries permeate every tissue of the body; indeed, no cell is more than a few thousandths of an inch from a capillary. The primary functions of the circulatory system take place in the capillary, i.e., the continuous unloading of cargo for the cells and the carrying away of the products of the cellular metabolism [3,4,5].
Capillaries are far away from the heart and narrow enough. Blood flow is already very low in arterioles due to the low heart pumping effect and high viscosity, so how is blood flow into capillaries possible? This happens due to the F-L effect [1]. It is a situation that happens in capillary-like tubes (diameter lower than 0.3 mm) [1]. The F-L effect was described independently by two Swedish scientists, Robin Fahraeus and Torsten Lindqvist, in 1931. According to this effect, when blood flow in a vessel with a diameter lower than 0.3 mm decreases in viscosity as the radius of the vessel decreases, flow is possible. When blood flows into microvessels like capillaries, resistance to blood flow increases, but the F-L effect reduces this resistance [6,7]. The F-L effect phenomenon was again repeated in 1931 by Martini et al. [8] and reinvestigated by Pries et al. [9]. This theory was improved by Haynes [10]. It describes the distribution of RBC over the vessel cross section area changes when it flows into a vessel with diameter lower than 0.3 mm. RBCs migrate to the center of the vessel and plasma close to the walls [10].
There are a large number of studies related to the F-L effect available. Most researchers have analyzed this effect and present different interpretations in their own ways. In many references, the F-L effect has been explained using Haynes’ marginal theory. Blood is assumed to be a two-layered fluid: one layer is a cell-free layer of plasma, and the other is a core layer, which is made from RBC. It is a suspended and deformable particle surrounded by a plasma layer [11,12,13,14,15,16,17,18,19,20,21]. A slower marginal layer reduces resistance and increases discharge but does not reduce the rates of energy dissipation [13], and a less-viscous, plasma-free layer of RBC flows close to the wall [19]. Sharan et al. [12] have assumed that the viscosity of the cell-free layer is different from that of plasma. This provides energy dissipation near the wall caused by RBC motion near the cell-free layer. They estimate effective viscosity in the cell-free layer and the F-L effect, depending upon the apparent viscosity, hematocrit, and tube diameter. The thickness of the marginal layer has also been determined [6]. This work was extended by Chebbi [14] with the help of a shear-induced model. The size of the core layer can be obtained with the help of this model. Hematocrit and the apparent viscosity level in the marginal layer can be determined. A fully developed hematocrit profile has also been obtained [15]. Consistency of the model with experimental results has also been given [16]. A continuum model with which to describe rheological properties, depending on the vessel’s diameter being less than 10 μ m , has been given by Medvedev [22]. For this, a generalized solution has been obtained. The effect of suspended deformability in capillaries of 8–11 μ m in diameter has been numerically simulated by Wang et al. [20]. A mesoscale simulation of blood in microvessels has been given by Bagchi [23]. About 2500 RBCs are considered as liquid capsules, allowing for deformation and migration in the microvessels with diameter 20–500 μ m . Botkin [11] assumed that a constant hematocrit level and plasma lubrication between RBCs and the capillary wall was dependent upon the wall and not the bifurcation effect. With the use of the finite element method, the analytical formula for the hydraulic resistance of the capillary depends upon the hematocrit level. Some authors have used the Casson law model to discuss flow into the capillary [18], but the Casson model is more applicable for the veins [24].
It is observed that several researchers have worked on the mathematical explanation of the F-L effect and have validated their results with experimental data. They have also discussed several physiological phenomena like the velocity profile of the core layer, the role of hematocrit in viscosity, and the thickness of the cell-free layer. They mostly assumed that blood is a two-phase fluid: one phase is the core (cell) layer, which is suspended in plasma; and the other phase is the cell-free layer, which is made from plasma. The plasma layer surrounds the suspended cell layer. Few researchers have considered the blood as single phase [24,25,26]. The constitutive equation for the single-phase blood flow model is τ = η e , where τ is stress, η is the viscosity coefficient, and e is the strain rate. The F-L effect for the single-phase blood flow model [24] is determined by
η e = η c 1 4 δ R η c η p 1 ,
where η e —effective viscosity; η c —viscosity of core layer; η p —viscosity of plasma; R—radius of capillary; and δ —thickness of plasma [24]. The constitutive equation for two-phase blood flow model is given by
τ = η e ,
where η —viscosity of blood; τ —stress tensor; and e —strain rate tensor [27]. For capillaries, the strain rate lies between 5 to 200 per second [24,28].
Blood is not an ideal fluid. In the usual case, the volume of its formed elements is 55% plasma and 45% other blood cells and particles. RBC covers 98% of the rest of the blood volume, and the remaining other 2% consists of WBC (white blood corpuscles), platelets, and other formed elements [29]. It is assumed that blood is separated into two layers when it enters the capillaries. These layers are plasma and the core layer. The plasma layer floats near the wall, and the core layer floats along the axis of the capillary [10]. The core layer is a homogeneous mixture of two phases, which are plasma and RBC. Without plasma, RBCs are not able to move; a low volume of plasma is necessary [3]. Plasma and RBCs are both liquid phases. RBCs are liquid in form, packed in a semipermeable membrane, capable of changing its shape and size, which reveals the nature of the liquid. The density of RBCs is higher than that of plasma. RBCs aggregate along the axis of the vessel and are capable of flow when mixed with plasma. Rouleax (plug) formation can be seen in the capillary. The core layer flows along the axis of the vessel, but plasma flows near the wall (see Figure 1).
A two-phase blood flow model is used to explain the F-L effect. Since capillaries are far away from the heart and near the tissue, their pumping effect can be neglected. The flow pattern has taken steady, laminar, and axial symmetry, with equal velocities of both the phases into the core layer. Physical quantities are expressed with the help of tensors, as a few quantities, like stress and strain rate, cannot be expressed in vector form.
Here, we discussed the model; extend the previous results; produce a comparative study of the results obtained from the model, including effective viscosity depending upon the radius of the capillary; and simulate results for human capillaries. In addition, the impact of various physiological parameters relevant to the flow are shown through various graphs.

2. Real Model

The whole model of blood flow in the capillary is presented in the following subsections.

2.1. Choice of Frame of Reference

Capillaries have taken approximately cylindrical tubes, although they are not perfectly cylindrical. An orthogonal curvilinear coordinate system is selected for the modeling [30]. All the physical quantities are expressed into tensorial form as biological law holds accurately in the given coordinate system.

2.2. Two-Phase Blood Flow

Let Y be the volume portion covered by the blood in unit volume, and Y can be replaced by   H a 100 , where H a is the hematocrit (the volume percentage of blood cells) [3]. Then, the volume portion covered by plasma will be 1 Y . If the mass ratio of blood cells to plasma is m r , then
m r = Y ρ c ( 1 Y ) ρ p ,
where ρ c and ρ p are the densities of the blood cells (core layer) and the plasma layer, respectively. Usually, this ratio is not constant; but in the case of the capillary, it is supposed to be constant [31].

2.3. Choice of Parameters

The following quantities are taken as parameters:
(a)
v k = v k ( x l , t ) ; l = 1 ,   2 ,   3 ; and   k = 1 ,   2 ,   3 , which is the velocity of blood at a given point x l in a space at a given time   t .
(b)
p = p ( x l , t ) ; l = 1 ,   2 ,   3 , which is the pressure of blood at a given point x l in a space at a given time   t .

2.4. Constitutive Equation

The constitutive equation of the two-layered blood flow model for capillaries is as follows:
τ l k = η m e l k p g l k , τ l k = η m g k i   v , i l p g l k ,
where η m = Y η c + ( 1 Y ) η p ;   η m , η c , and η p are the viscosity of the mixture, core layer, and plasma layer, respectively; and v , i l is the covariant derivative of the velocity tensor [32].

2.5. Boundary Conditions

Two boundary conditions (see Figure 2) have been taken for modeling purposes, given as follows:
(a)
Maximum velocity of blood flow at the axes; i.e., when the radius of vessels is r = 0 , the velocity of blood flow is at a maximum,   i . e . ,   v = v 0 (say).
(b)
No slip condition (velocity on the boundary is zero); i.e., when theradius is   r = R , then   v = 0 .

3. Mathematical Formulation

The formulation of the proposed problem combines the equation of continuity and the equation of motion. Details are given in following subsections.

3.1. Equation of Continuity

When there is no source and sink in any flowing fluid region, the mass of the fluid is conserved in that region. As there is no source and sink in the human circulatory system [33], the law of conservation of mass can be applied here, which means the mass of blood that enters into the system is equal to the mass of blood that exist outside the system [28].
It is assumed that both phases in the core layer move with common velocity into the capillary. From the principle of conservation of mass, the equation of continuity for RBC phase is given as follows:
( Y ρ c ) t + ( Y ρ c v l ) , l = 0 .
The equation of continuity for plasma phase is also given as follows:
( 1 Y ) ρ p t + ( ( 1 Y ) ρ p v l ) , l = 0 ,
where v is the common velocity of both phases, ρ c and ρ p are densities of core layer and plasma layer, respectively. ( Y ρ c v l ) , l is a covariant derivative of ( Y ρ c v l ) with respect to X l . In the same way, ( ( 1 Y ) ρ p v l ) , l with respect to X l [3,34,35].
We define the uniform density ρ m as follows:
1 + r 0 ρ m = r 0 ρ c + 1 ρ p .
Above Equations (4) and (5) can be combined together as follows:
ρ m t + ( ρ m v l ) , l = 0 .
As we know, blood is an incompressible fluid; hence, ρ m is constant. Thus, the equation of continuity (7) for blood flow takes the following form:
v , l l = 0 , i . e . ,   v l X l + v l g g X l = 0 , 1 g ( g v l ) , l = 0 .

3.2. Equation of Motion

The total momentum of any system is conserved in the absence of external force. The effect of external force in our human circulatory system is zero in normal situations. So, the law of conservation of momentum can fully apply to this system. The hydrodynamical pressure p between the two phases of blood can be supposed to be uniform because both layers, core and plasma, are always in an equilibrium state in the blood. Now, applying the principle of conservation of momentum, we obtain the equation of motion for the phase of red blood cells as follows [3,34,35]:
Y ρ c v l t + Y ρ c v k v , k l = Y p , k g l k + Y η c ( g k i v , i l ) , k .
In a similar manner, the equation of motion for plasma in the core layer will be as follows:
1 Y ρ p v l t + { 1 Y ρ p v k } v , k l = ( 1 Y ) p , k g l k + ( 1 Y ) η p ( g k i v , i l ) , k .
Now, adding Equations (9) and (10) and using relation (6), the equation of motion for blood flow with the combination of both phases will be reduced as follows:
ρ m v l t + ρ m v k v , k l = p , k g l k + η m g k i v , i l , k ,
where ρ m = Y ρ c + ( 1 Y ) ρ p is the density and η m = Y η c + ( 1 Y ) η p is the viscosity of a mixture of blood (blood cells and plasma), respectively.
Equations (8) and (11) are transformed in cylindrical form.
For cylindrical coordinates,   x 1 = r , x 2 = θ ,   x 3 = z .
Components of metric tensor g l k 3 × 3 are   g 11 = 1 ,     g 22 = r 2 ,     g 33 = 1 , and the rest of the entries are zero.
The determinant of metric tensor is   g = r 2 . The conjugate metric tensor [30] is as follows:
g l k = c o f a c t o r   o f   g l k g .
Again, components of the conjugate metric tensor are   g 11 = 1 ,   g 22 = 1 r 2 ,   a n d   g 33 = 1 , and the rest of the entries are zero.
Christoffel’s symbol of the second kind, i . e . , l i   j = g l k i j , k = 1 2 g l k g j k x i + g i k x j g i j x k , is as follows:
1 2   2 = r , 2 2   1 = 2 1   2 = 1 r , and the rest are zero [30,36].
The relationship between the physical components and contravariant components [30] of the velocity of blood flow are as follows:
g 11 v 1 = v r v r = v 1 , g 22 v 2 = v θ v θ = r v 2 , and   g 33 v 3 = v z v z = v 3 , p , k g l k g l l p , k g l k .
Now, the governing equations of blood flow are taken in the following form:
Equation of continuity:
1 r   ( r v r ) r + 1 r v θ θ + v z z = 0 .
The equation of motion is given in r ,   θ , and z components, as below:
r   component is
ρ m v r t + v r v r r + v θ r v r θ v θ 2 r + v r z = p r + η m r 1 r r v r r + 1 r 2 2 v r θ 2 2 r 2 v θ θ + 2 v r z 2 ;
θ component is
ρ m v θ t + v r v θ r + v θ r v θ θ v r v θ r + v z v θ z = 1 r p θ + η m r 1 r r v θ r + 1 r 2 2 v θ θ 2 + 2 r 2 v r θ + 2 v θ z 2 ;
z component is
ρ m v z t + v r v z r + v θ r v z θ + v z v z z = p z + η m 1 r r r v z r + 1 r 2 2 v z θ 2 + 2 v   z z 2 .

4. Solution of the Problem

The blood flow is symmetric with respect to the axis of the capillary, i.e.,     v θ = 0 , and v r   ,   v z , and p are not, depending on   θ . Since only one component of velocity along the axis is effective, we have   v r = 0 ,     v θ = 0 ,   a n d   v z = v (say).
Since the flow is steady, p t = v r t = v θ t = v z t = 0 [33].
Hence, the equation of continuity (12) reduces as follows:
v z z = 0 v z = v r .
The transformed equations of motion, along with the r ,   θ , and z components, are given as follows:
The r   component of equation of motion reduces to
ρ m × 0 = p r + η m 0 , p r = 0 p = p z .
The θ component of equation of motion vanishes because blood flow is free from rotation, and applied pressure is in the normal direction.
Similarly, the z component of the equation of motion reduces to
ρ m v z 0 = p z + η m 1 r r r v r r , 0 = p z + η m 1 r r r v r r , η m 1 r d d r r d v ( r ) d r = P ,
where the pressure gradient is
p z = P .
Now, we solve Equation (18) using boundary conditions v = v 0   a t   r = 0 and v = 0 at r = R (no slip condition) and obtain the velocity of mixture of blood flow, i.e.,
v = P 4 η m R 2 r 2 .
Therefore, the velocity of the plasma layer (as if the whole capillary was filled with plasma) is
v p = P 4 η p ( R 2 r 2 ) ,   R δ r R .
Similarly, velocity of the core layer (if core the layer was filled in the whole capillary and the relative velocity of both layers was also added) is
v m = P 4 η m R 2 r 2 + P 4 η m R 2 ( R δ ) 2 η m η p 1 ,   0 r R δ ,
where R is the radius of the capillary and   δ is the thickness of the plasma layer (see in Figure 3) and is supposed to be independent of   R .

5. Result and Discussion (Biophysical Interpretation)

The flow flux of blood in the capillary [24] is defined as
Q = 0 R δ v m 2 π r   d r + R δ R v p 2 π r   d r , = 0 R δ P 4 η m R 2 r 2 + P 4 η m R 2 ( R δ ) 2 η m η p 1 2 π r   d r + R δ R P 4 η p ( R 2 r 2 ) 2 π r   d r .
Using (20) and (21) in the above expression, we obtain
Q = π P R 4 8 η p 1 1 δ R 4 1 η p η m
If the whole capillary is supposed to be filled with only one Newtonian blood, whose effective viscosity is   η e , then the flow flux of the blood [24] is
Q = π P R 4 8 η e .
Comparing Formulae (23) and (24), we obtain the effective viscosity as follows
η e = η p 1 1 δ R 4 1 η p η m 1 ,
which depends upon the radius R of the capillary. Since δ R 1 , the Formula (25) can be expanded in the higher power of δ R , taking first approximation of δ R ; then, it is reduced in the following form:
η e = η m 1 4 δ R η m η p 1   .
It is observed from Equation (26) that the effective viscosity decreases with the decrease in the radius of the capillary   R . Thus, the Fahraeus–Lindqvist effect is explained mathematically.

5.1. Numerical Simulation

A comparison and validation have conducted between the results obtained from the two-phase blood flow model and the experimental data after validation of the results simulated for the human capillary. To validate our model with the experimental result [1], the following numerical values were given to different parameters:
δ = 2 3 R ,   η c = 0.0249   P a s c a l s ,   η p = 0.00158 0.00172 Pascal s, H a = 11.6 26 [31]. For validation of the model value of the R taken from the experimental data, it is simulated on a human capillary of 4–10 μ m [3]. MATLAB2020a software is used for numerical simulation and to draw all figures.
P = p z = p r e s u r e   d r o p l e n g t h   o f   t h e   c a p i l l a r y = p f p i z f z i , p f = 2 3 S + D 2 + D 3 = p r e s s u r e   a t   v e n u l e s ,   p i = S + D 2 + D 3 = p r e s s u r e   a t   a r t e r i o l e s , S 110   m m h g S y s t o l i c   p r e s s u r e , D 70   m m h g = ( D i a s t o l i c   p r e s s u r e ) .
Table 1, Table 2, Table 3 and Table 4 for different series have been taken from the experimental data [1]. Series 1, 2, and 4 of the experiments were given by Torsten Lindqvist, and Series 3 was given by Robin Fahraeus, shown below.

5.2. Validation of Two-Phase Blood Flow Model

Figure 4 shows the graph between η e and   R . The effective viscosity η e is obtained from a single-phase model (from Equation (1)), a two-phase-model (from Equation (26)), and the experimental data of series-1. For H a = 15.6, the two-phase model is much closer to the experimental data in Series 1 than the single-phase model. So next, the impact of different physiological parameters relevant to the flow, with help of the two-phase blood flow model, was analyzed.
Figure 5 shows variation in η e and R , which are obtained from the two-phase blood flow model and experimental Series 1 to 4. In the figures of each series, η e decreases with radius   R . For H a = 15.6, 18.2, 19.8, and 27.6, the proposed model is much closer to the experimental Series 1, 2, 3, and 4, respectively (see Figure 5).
The proposed two-phase blood flow model gives a better approximation than the single-phase blood flow model. This model will be used for further studies on the human capillary. The internal diameter of the capillary lies between 4 and 10 μ m [3]. For this purpose, a capillary with radius R = 5   μ m , length Δ z = 6   μ m , and thickness of plasma layer δ = 2 3 R was used [37].

5.3. Effect of H a , η p , and P on η e

For a capillary F-L effect, as shown in Figure 6a–c, the radii were taken from 0 < r R . η e decreases rapidly in the capillary with a radius less than 3   μ m . In Figure 6a,b, η e decreases with decreasing values of r . Figure 6a shows variation in η e with radius 0 < r R for different values of H a . From Figure 6a, it is observed that η e decreases with the higher concentration of H a . The concentration of RBC in the capillaries increases but it decreases η e along with r ; hence, flow is possible. Figure 6b shows the impact of   η p on   η e ; it decreases with higher values of viscosity of the plasma layer. Figure 6c shows the impact of pressure drop P on   η e . It is assumed that the capillaries are far away from the heart and proximate to the tissue; hence, its pumping effect has been negligible. This effect can also be seen in Figure 6c. Increasing the value of P does not much affect η e .

5.4. Effect of H a , η p , and P   on   v , v m ,   a n d   v p

The velocity profile v for the mixture of blood is shown in Figure 7a–c. In each figure ,   v increases with the decreasing radius of the capillary, tending to zero at r = R . Figure 7a shows the effect of   H a on v . It can be observed from Figure 7a that a higher volume percentage of RBC decreases   v . The impact of pressure drop P on whole blood is shown in Figure 7b. Although the heart pumping effect in not very effective into the capillaries, its larger value increases blood velocity. Figure 7c shows the effect of plasma viscosity on   v . η p slightly affects velocity v . Higher values of η p decrease v , but they do not greatly influence it.
Blood is separated into two layers when it flows from the capillary. Graphs of combined velocity profiles have been plotted for the core layer v m (for radius 0 r R δ ) and the plasma layer v p   R δ r R (see Figure 8a–c). Figure 8a shows the impact of the volume percentage of RBC on both layers. As the plasma layer is free from RBC, which flows near the wall of the capillary, the impact of H a cannot be seen in the velocity profile of the plasma layer. Its effect can be seen only on the core layer. v m decreases as H a increases. The pressure drop also increases the velocity of both layers (see Figure 8b). The impact of the viscosity of the plasma layer on both layers is shown in Figure 8c. The velocity of both layers increases as plasma viscosity decreases.

5.5. Effect of H a , η p , and P   on Q

Variation in Q with respect to radius r is shown in Figure 9a–c. In each figure, Q decreases rapidly with the decreasing radius r , and after the radius 2   μ m changes in Q cannot be seen explicitly; it becomes stable for radii less than 2   μ m . The impact of H a is shown in Figure 9a. The graphs between Q and r for different values of H a overlap, which that means Q is not much affected by   H a . Figure 9b shows the effect of the pressure drop on   Q . Increasing values of P also increase Q , but for radii less than 2   μ m , the curves coincide, which means that for small- radius vessels, changes in Q become static. Figure 9c shows the impact of plasma viscosity on   Q . The effect of η p on Q can be seen explicitly for the vessel with a radius larger than 2   μ m . Q increases as   η p decreases, but for a radius less than 2   μ m , Q does not change as much.
Figure 8. Variation in v m with radius ( 0 r R δ ) and v p with radius ( R δ r R ) , respectively, for different values of (a) hematocrit H a , (b) pressure drop   P , and (c) plasma viscosity η p .
Figure 8. Variation in v m with radius ( 0 r R δ ) and v p with radius ( R δ r R ) , respectively, for different values of (a) hematocrit H a , (b) pressure drop   P , and (c) plasma viscosity η p .
Symmetry 16 00728 g008

6. Conclusions

A two-phase microcirculation model that explains the F-L effect has been presented here. For the validation of the model, the obtained results are compared with the experimental data and single-phase blood flow model, respectively. Further, we simulated these results on the human capillary. The effects of various physiological parameters have been shown through different graphs. We present the conclusions of this study as follows:
  • A mathematical expression has been obtained between η e and R which successfully explains the F-L effect.
  • We have obtained effective, rather than apparent, viscosity.
  • With the help of this model, extension of Haynes theory is possible because for the relative momentum of RBC, the plasma is mixed in the core layer, and both phases are homogeneously distributed in the core layer.
  • Experimental data and the results obtained from the presented model have the same trends, but this model provides an extension of the F-L effect.
  • The presented model successfully explains all the physical phenomena which occur in the capillary.
  • The flow flux and velocity of the blood decreases in the capillary due to the higher concentration of H a , but flow is possible because effective viscosity decreases as the radius of the capillary decreases.
  • The pressure gradient can increase the effective viscosity and flow flux, but it is not very effective in terms of blood flowing into capillaries as they are remote from the heart and proximate to the human tissue.
There are several two-phase blood flow models of the F-L effect, but uses of tensors and the RBCs of the core layer considered are liquid in form and packed in a semipermeable membrane, capable of changing its shape and size, which demonstrates the nature of the liquid, not of a suspended particle, moving the model closer to the reality. This model not only successfully explains the F-L effect but is also able to determine the impact of several physiological parameters on blood flow in the capillary.

Author Contributions

P.M.: conceptualization, data curation, formal analysis, investigation, methodology, project administration, resources, software, supervision, validation, visualization, writing—original draft, writing—review and editing; V.U.: conceptualization, formal analysis, investigation, methodology supervision, validation, visualization; S.K.C.: data curation; formal analysis; V.C.: software supervision, visualization, formal analysis, investigation; D.K.: data curation, formal analysis, investigation, software supervision, visualization, writing—review and editing, project administration, resources. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

This research paper does not contain new data; the results of the experimental data are presented within the article.

Acknowledgments

The second author is thankful to CSIR-UGC, New Delhi, India, for the financial support under the JRF (UGC Ref. No.:1107/(CSIR-UGC NET DEC. 2017)) scheme.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

τ l k   Stress tensor [N/m2]
e l k Strain rate tensor [S−1]
g l k Metric tensor
g l k Conjugate metric tensor
l i   j Christoffel’s symbol of second kind
pPressure [N/m2]
η Viscosity constant [Pascal s]
Y Volume portion covered by the blood cells in unit volume
H a   Hematocrit [volume percentage of RBC]
δ Thickness of plasma layer
ρ c Density of core layer [kg/m3]
ρ p Density of plasma layer [kg/m3]
ρ m Density of mixture of blood [kg/m3]
m r Mass ratio of blood cells to plasma
η m Viscosity of mixture [Pascal s]
η c Viscosity of core layer [Pascal s]
η p Viscosity plasma layer [Pascal s]
η e Effective viscosity [Pascal s]
v Velocity of mixture of blood [m/s]
v p Velocity of plasma [m/s]
v m Velocity of core layer [m/s]
v r Radial component of velocity [m/s]
v θ Angular component of velocity [m/s]
v z Axial component of velocity [m/s]
D Diameter of capillary [m]
P Pressure gradient [Pascal/m]
p Pressure drop [Pascal]
RRadius of capillary [m]
LLength of capillary [m]
Abbreviations
F-LFahraeus–Lindqvist
T.L.Torsten Lindqvist
R.F.Robin Fahreaeus
R B C Red blood corpuscle
WBCWhite blood corpuscle

References

  1. Fahraeus, R.; Lindqvist, J.T. The viscosity of the blood in narrow capillary tubes. Am. J. Physiol. 1931, 96, 562–568. [Google Scholar] [CrossRef]
  2. Wilson, R. Anatomy and Physiology in Health and Illness, 11th ed.; Churchill Livingstone: London, UK, 2014. [Google Scholar]
  3. Guyton, A.C.; Hall, J.E. Textbook of Medical Physiological, 11th ed.; Elsevier Inc.: New Delhi, India, 2006. [Google Scholar]
  4. Dalkara, T. Cerebral Microcirculation: An Introduction; Springer: Berlin/Heidelberg, Germany, 2015; pp. 655–680. [Google Scholar]
  5. Cipolla, M.J. The Cerebral Circulation; Morgan and Claypool Life Sciences: San Rafael, CA, USA, 2009; ISBN 9781615040131. [Google Scholar]
  6. Toksvang, L.N.; Berg, R.M.G. Using a classic paper by Robin Fåhræus and Torsten Lindqvist to teach basic hemorheology. Adv. Physiol. Educ. 2013, 37, 129–133. [Google Scholar] [CrossRef] [PubMed]
  7. Secomb, T.W.; Pries, A.R. Blood viscosity in microvessels: Experiment and theory. C. R. Phys. 2013, 14, 470–478. [Google Scholar] [CrossRef] [PubMed]
  8. Martini, P.; Pierach, A.; Scheryer, E. Die stromung des blutes in engen gefäßen. Eine abweichung vom poiseuille’schen gesetz. Dtsch. Arch. Klin. Med. 1931, 169, 212–222. [Google Scholar]
  9. Pries, A.R.; Neuhas, D.; Gaehtgens, P. Blood viscosity in tube flow: Dependence on diameter and hematocrit. Am. J. Physiol. Heart Circ. Physiol. 1992, 263, H1770–H1778. [Google Scholar] [CrossRef] [PubMed]
  10. Haynes, R.F. Physical basis of the dependence of blood viscosity on tube radius. Am. J. Physiol. 1960, 198, 1193–1200. [Google Scholar] [CrossRef] [PubMed]
  11. Botkin, N.D.; Kovtanyuk, A.E.; Turova, V.L.; Sidorenko, I.N.; Lampe, R. Accounting for tube haematocrit in modeling of blood flow in cerebral capillary networks. Comput. Math. Methods Med. 2019, 2019, 4235937. [Google Scholar] [CrossRef] [PubMed]
  12. Sharan, M.; Popel, A.S. A two-phase model for flow of blood in narrow tubes with increased effective viscosity near the wall. Biorheology 2001, 38, 415–428. [Google Scholar] [PubMed]
  13. Ascolese, M.; Farina, A.; Fasano, A. The Fåhræus-Lindqvist effect in small blood vessels: How does it help the heart. J. Biol. Phys. 2019, 45, 379–394. [Google Scholar] [CrossRef]
  14. Chebbi, R. Dynamics of blood flow: Modeling of the Fåhræus–Lindqvist effect. J. Biol. Phys. 2015, 41, 313–326. [Google Scholar] [CrossRef]
  15. Chebbi, R. Dynamics of blood flow: Modeling of Fahraeus and Fahraeus-Lindqvist effects using a shear-induced red blood cell migration model. J. Biol. Phys. 2018, 44, 591–603. [Google Scholar] [CrossRef] [PubMed]
  16. Chebbi, R. A two-zone shear-induced red blood cell migration model for blood flow in microvessels. Front. Phys. 2019, 7, 206. [Google Scholar] [CrossRef]
  17. Farina, A.; Rosso, F.; Fasano, A. A Continuum mechanics model for the Fahraeus-Lindqvist effect. J. Biol. Phys. 2021, 47, 253–270. [Google Scholar] [CrossRef] [PubMed]
  18. Farina, A.; Fasano, A.; Rosso, F. Mathematical models for some aspects of blood microcirculation. Symmetry 2021, 13, 1020. [Google Scholar] [CrossRef]
  19. Possenti, L.; Di Gregorio, S.; Gerosa, F.M.; Raimondi, G.; Casagrande, G.; Costantino, M.L.; Zunino, P. A computational model for microcirculation including Fahraeus-Lindqvist effect, plasma skimming and fluid exchange with the tissue interstitium. Int. J. Numer. Methods Biomed. Eng. 2019, 35, e3165. [Google Scholar] [CrossRef] [PubMed]
  20. Wang, T.; Xing, Z. Characterization of blood flow in capillaries by numerical simulation. J. Mod. Phys. 2010, 1, 349. [Google Scholar] [CrossRef]
  21. Medvedev, A.E.; Fomin, V.M. Two-phase blood-flow model in large and small vessels. Dokl. Phys. 2011, 56, 610–613. [Google Scholar] [CrossRef]
  22. Medvedev, A.E. Two-phase blood-flow model. Russ. J. Biomech. 2013, 17, 18–32. [Google Scholar]
  23. Bagchi, P. Mesoscale simulation of blood flow in small vessels. Biophys. J. 2007, 92, 1858–1877. [Google Scholar] [CrossRef]
  24. Kapoor, J.N. Mathematical Models in Biology and Medicine; E.W.P.: New Delhi, India, 1992; ISBN 8185336822. [Google Scholar]
  25. Debnath, L. On a micro-continuum model pulsatile blood flow. Acta Mech. 1976, 24, 165–177. [Google Scholar] [CrossRef]
  26. Devi, G.; Devanathan, R. Peristaltic motion of a micropolar fluid. Proc. Indian Acad. Sci. 1975, 81, 149–163. [Google Scholar] [CrossRef]
  27. Upadhyay, V.; Chaturvedi, S.K.; Upadhyay, A. A mathematical model on effect of stenosis in two phase blood flow in arteries remote from the heart. Int. Acad. Phys. Sci. 2012, 16, 247–257. [Google Scholar]
  28. Fung, Y.C. Biomathematics Mechanical Properties of Living Tissues; Springer: New York, NY, USA, 1981; pp. 118–182. [Google Scholar]
  29. Sherman, I.W.; Sherman, V.G. Biology—A Human Approach; Oxford University Press: New York, NY, USA; Oxford, UK, 1989; pp. 276–277. [Google Scholar]
  30. Murray, R.S. Vector Analysis and an Introduction to Tensor Analysis; Schaum’s Outline Series; McGraw-Hill: New York, NY, USA, 2021. [Google Scholar]
  31. Maurya, P.; Upadhyay, V.; Chaturvedi, S.K.; Kumar, D. Mathematical study and simulation on stenosed carotid arteries with the help of two-phase blood flow model. Can. J. Chem. Eng. 2023, 101, 5468–5481. [Google Scholar] [CrossRef]
  32. De, U.C.; Abosos, S.A.; Joydeep, G.S. Tensor Calculus, 2nd ed.; Alpha Science International, Ltd.: Oxford, UK, 2008. [Google Scholar]
  33. Landau, L.D.; Lipchitz, E.M. Fluid Mechanics; Pergamon Press: Oxford, UK, 1959. [Google Scholar]
  34. Wiat, L.; Fine, J. Applied Bio Fluid Mechanics; McGraw-Hill Companies: New York, NY, USA, 2007. [Google Scholar]
  35. Mazumdar, J.N. Bio Fluid Mechanics; World Scientific: Singapore, 2004. [Google Scholar]
  36. Spain, B. Tensor Calculus: A Concise Course; Courier Corporation: North Chelmsford, MA, USA, 2003. [Google Scholar]
  37. Batra, S.; Rakusan, K. Capillary length, tortuosity, and spacing in rat myocardium during cardiac cycle. Am. J. Physiol.-Heart Circ. Physiol. 1992, 263, H1369–H1376. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of blood flow in the capillary.
Figure 1. Schematic diagram of blood flow in the capillary.
Symmetry 16 00728 g001
Figure 2. Boundary conditions on the vessel.
Figure 2. Boundary conditions on the vessel.
Symmetry 16 00728 g002
Figure 3. Schematic diagram of blood flow in the capillary.
Figure 3. Schematic diagram of blood flow in the capillary.
Symmetry 16 00728 g003
Figure 4. Comparison of single-phase and two-phase blood flow models with experimental data [1] when   P = 0.2591   p a s c a l , η p = 0.00158   p a s c a l s , η c = 0.0216   p a s c a l s ,   a n d   H a = 15.6 .
Figure 4. Comparison of single-phase and two-phase blood flow models with experimental data [1] when   P = 0.2591   p a s c a l , η p = 0.00158   p a s c a l s , η c = 0.0216   p a s c a l s ,   a n d   H a = 15.6 .
Symmetry 16 00728 g004
Figure 5. Graph of effective viscosity vs. radius of capillary with experimental result [1] for different values of hematocrit H a when   P = 0.2591   p a s c a ,   η p = 0.00158   p a s c a l , a n d   η c = 0.0216   p a s c a l .
Figure 5. Graph of effective viscosity vs. radius of capillary with experimental result [1] for different values of hematocrit H a when   P = 0.2591   p a s c a ,   η p = 0.00158   p a s c a l , a n d   η c = 0.0216   p a s c a l .
Symmetry 16 00728 g005
Figure 6. Variation in effective viscosity η e   with radius 0 < r R of the capillary for different values of (a) hematocrit H a , (b) plasma viscosity η p , and (c) pressure drop P .
Figure 6. Variation in effective viscosity η e   with radius 0 < r R of the capillary for different values of (a) hematocrit H a , (b) plasma viscosity η p , and (c) pressure drop P .
Symmetry 16 00728 g006
Figure 7. Variation in the velocity of blood mixture with radius ( 0 r R ) for different values of (a) hematocrit H a , (b) pressure drop P , and (c) plasma viscosity η p .
Figure 7. Variation in the velocity of blood mixture with radius ( 0 r R ) for different values of (a) hematocrit H a , (b) pressure drop P , and (c) plasma viscosity η p .
Symmetry 16 00728 g007
Figure 9. Variation in flow flux Q with radius ( 0 r R δ ): (a) hematocrit H a ; (b) pressure drop P ; (c) plasma viscosity η p .
Figure 9. Variation in flow flux Q with radius ( 0 r R δ ): (a) hematocrit H a ; (b) pressure drop P ; (c) plasma viscosity η p .
Symmetry 16 00728 g009
Table 1. Experiment Series 1, Torsten Lindqvist (T.L.).
Table 1. Experiment Series 1, Torsten Lindqvist (T.L.).
Series 1
Relative   Viscosity   of   Plasma   η p = 1.63 × 10 3 P a s c a l S
S. No.Length of Tube (m)Diameter (m)Relative Viscosity (Pascal-S)
10.1265 5.05 × 10 4 4.60 × 10 3
20.1210 3.11 × 10 4 4.54 × 10 3
30.0456 1.75 × 10 4 4.22 × 10 3
40.0500 1.20 × 10 4 3.98 × 10 3
50.0326 8.10 × 10 5 3.65 × 10 3
60.0118 6.50 × 10 5 3.38 × 10 3
70.0118 4.70 × 10 5 2.95 × 10 3
80.0126 4.00 × 10 5 2.79 × 10 3
Table 2. Experiment Series 2, Torsten Lindqvist (T.L.).
Table 2. Experiment Series 2, Torsten Lindqvist (T.L.).
Series 2
Relative   Viscosity   of   Plasma   η p = 1.65 × 10 3 P a s c a l S
S. No.Length of Tube (m)Diameter (m)Relative Viscosity (Pascal-S)
10.1255 2.77 × 10 4 4.96 × 10 3
20.0456 1.75 × 10 4 4.54 × 10 3
30.0500 1.20 × 10 4 4.22 × 10 3
40.0118 4.70 × 10 5 3.98 × 10 3
Table 3. Experiment Series 3, Robin Fahraeus (R.F.).
Table 3. Experiment Series 3, Robin Fahraeus (R.F.).
Series 3
Relative   Viscosity   of   Plasma   η p = 1.60 × 10 3 P a s c a l S
S. No.Length of Tube (m)Diameter (m)Relative Viscosity (Pascal-S)
1---- 4.60 × 10 4 5.21 × 10 3
20.1000 2.77 × 10 4 5.10 × 10 3
30.0850 1.05 × 10 4 4.28 × 10 3
4---- 5.20 × 10 5 3.52 × 10 3
Table 4. Experiment Series 4, Torsten Lindqvist (T.L.).
Table 4. Experiment Series 4, Torsten Lindqvist (T.L.).
Series 4
Relative   Viscosity   of   Plasma   η p = 1.72 × 10 3 P a s c a l S
S. No.Length of Tube (m)Diameter (m)Relative Viscosity (Pascal-S)
10.1225 2.77 × 10 4 6.80 × 10 3
20.0456 1.75 × 10 4 6.28 × 10 3
30.0500 1.20 × 10 4 5.59 × 10 3
40.0279 8.10 × 10 5 5.08 × 10 3
50.0070 6.50 × 10 5 4.71 × 10 3
60.0118 4.70 × 10 5 4.28 × 10 3
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Upadhyay, V.; Maurya, P.; Chaturvedi, S.K.; Chaurasiya, V.; Kumar, D. A Mathematical Analysis and Simulation of the F-L Effect in Two-Layered Blood Flow through the Capillaries Remote from the Heart and Proximate to Human Tissue. Symmetry 2024, 16, 728. https://doi.org/10.3390/sym16060728

AMA Style

Upadhyay V, Maurya P, Chaturvedi SK, Chaurasiya V, Kumar D. A Mathematical Analysis and Simulation of the F-L Effect in Two-Layered Blood Flow through the Capillaries Remote from the Heart and Proximate to Human Tissue. Symmetry. 2024; 16(6):728. https://doi.org/10.3390/sym16060728

Chicago/Turabian Style

Upadhyay, Virendra, Pooja Maurya, Surya Kant Chaturvedi, Vikas Chaurasiya, and Dinesh Kumar. 2024. "A Mathematical Analysis and Simulation of the F-L Effect in Two-Layered Blood Flow through the Capillaries Remote from the Heart and Proximate to Human Tissue" Symmetry 16, no. 6: 728. https://doi.org/10.3390/sym16060728

APA Style

Upadhyay, V., Maurya, P., Chaturvedi, S. K., Chaurasiya, V., & Kumar, D. (2024). A Mathematical Analysis and Simulation of the F-L Effect in Two-Layered Blood Flow through the Capillaries Remote from the Heart and Proximate to Human Tissue. Symmetry, 16(6), 728. https://doi.org/10.3390/sym16060728

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop