Micromorphic theory as a model for blood in the microcirculation: correction and analysis

This paper analyzes the applicability of Eringen’s Generalized Continuum Theories as a model for human blood in the microcirculation. The applied theory considers a fluid with a fully deformable substructure, namely a micromorphic fluid. This analysis is motivated by the fact that blood itself can be considered a suspension of deformable particles, i.e., red blood cells (RBCs), suspended in a Newtonian fluid, i.e., blood plasma. As a consequence, non-Newtonian phenomena such as shear-thinning are observed in blood. To test the micromorphic fluid as a model for blood, the solution for the velocity and the motion of substructure is determined for a cylindrical pipe flow and compared to experimental results of blood flow through narrow glass capillaries representing idealized blood vessels. A similar analysis was also conducted by Kang and Eringen in 1976, but it contains some misprints and minor errors regarding the mathematical expressions and subsequent discussion which are corrected in this paper. For certain material parameters, the micromorphic fluid models capture high-shear blood flow in narrow glass capillaries very well. This concerns both the velocity profiles and the shear-thinning behavior. Furthermore, a parameter study reveals that the flexibility of substructure governs the micromorphic shear-thinning. In this regard, parallels can be drawn to the shear-thinning of human blood, which is also induced by the deformability of RBCs. This makes the micromorphic fluid a complex but accurate model for human blood, at least for the considered experiments.

blood in narrow blood vessels, developed by Eringen in 1964. For a better understanding of the considered approach and the related experiments, a brief introduction to the two fields relevant for this paper, namely hemodynamics and Eringen's micro-continuum theories, is given below. Then, the governing equations of a micromorphic fluid are stated. They are subsequently used to determine the motion and the micromotion of a micromorphic fluid in a cylindrical tube. After validating the derived results, the motion of the micromorphic fluid is compared to human blood flowing through narrow glass capillaries in Sect. 4. This problem was already solved and discussed by Kang and Eringen in 1976, leading to promising conclusions [22]. However, their solution contains some misprints and the derivation of their results was not documented, so the results are not reproducible. The purpose of this paper is to evaluate, whether the conclusions made in [22] are still valid.

Hemodynamics
Human blood can be construed as a tissue, an emulsion or a two-phase fluid, depending on the field of study. However, in the rheological context of this paper, it is considered as a fluid in which particles are suspended. Most importantly, it can, in general, not be considered as a Newtonian fluid. This is because it has many properties that are not covered by classic fluid mechanics. The aim of this section is to point out the differences to a Newtonian fluid by listing and-as far as possible-explaining several phenomena leading to these differences. To be able to do so, it is useful to understand, what exactly blood is, and thus, what components it is composed of.
Human blood and that of vertebrates, in general, consists of an aqueous solution, namely plasma, and different particles, i.e., RBCs, white blood cells (WBCs) and platelets. Of these, the RBCs have the highest volume fraction, or Hematocrit (Ht), of about 40% to 45% [33]. The blood plasma makes up the highest volume fraction of blood and can be considered a Newtonian fluid. It consists mainly of water, but also contains "organic molecules, proteins, and salts," [3, p. 435]. The plasma viscosity at 37 • C usually takes values of 1.1mPa s to 1.3mPa s [23], whereas the viscosity of whole blood cannot be described by a single value, because it varies with changing flow situations, as discussed below. As a consequence of the high hematocrit, this behavior is primarily attributed to the RBCs [26]. Under special circumstances, WBCs and platelets can also have a severe impact on the motion of blood, for example through platelet adhesion and aggregation due to injured blood vessel and high shear rates leading to thrombosis, as explained in [28]. But these special circumstances are not considered in this work.
As mentioned above, the hematocrit has a major influence on blood flow properties: For example, an increased hematocrit yields an increased viscosity for all vessel sizes and shear rates [20]. The larger viscosity due to more cells suspended in plasma is a logical consequence of rigid particles disturbing the streamlines and thereby increasing the flow resistance [5]. At a sufficiently large shear rate, a migration of RBCs toward the center of the vessel is observed, resulting in a high local hematocrit in the inner vessel and a RBC-free plasmatic zone near the vascular wall. Consequently, the local viscosity near the vessel wall corresponds to that of plasma, whose value is close to that of water and many times smaller than that of blood in large vessels. As a result, the blood flows without any particle-wall-interaction-induced flow resistance through the vessel, which promotes a decreased viscosity for whole blood. While the formation of a plasmatic layer occurs for every vascular radius, experiments indicate that the absolute thickness of this layer does not change along the vessel radius but is fixed at about 1.8µm [31,32]. Hence, the reduced viscosity for larger radii is negligible. This effect is generally referred to as shear-thinning or, in the context of blood Fåhraeus-Lindqvist effect, see [13]. Shear-thinning is characterized by the apparent viscosity, which is derived from the Hagen-Poiseuille equation [25,33]: (1.1) Here, p/L, R and Q denote the pressure drop along the distance L, the vessel radius and the volume flow rate, respectively. For sufficiently large shear rates, the volume flow rate Q is approximately proportional to the pressure drop p/L multiplied by a factor M(R), which only depends on the vessel radius R. Hence, the asymptotic apparent viscosity is defined as a function of R [19,21]: . (1. 2) The error of this approximation was found to be negligible for pseudo-shear rates v > 50s −1 which is a measure for the actually occurring shear rates [29]: Here, v m denotes the volumetrically averaged velocity. It has been shown that for R < 300 µm the apparent viscosity decreases significantly, monotonously with a decreasing vessel radius [13]. Consequently, in smaller vessels blood becomes continuously less viscous. Since most of RBCs are packed near the vessel center, they flow with similar velocity as the surrounding particles, forming a plug-shaped velocity profile, cf., [15]. It can be proven that the cell-free layer formation and the Fåhraeus-Lindqvist effect are a result of RBC deformability inducing axial migration, whereas the exact mechanisms leading to fluidity in the core region are not clear yet [16]. While recent studies indicate that a simple deformation of RBCs along the streamlines might be the cause of shear-thinning in the vessel core region [24], another possible explanation is an extraordinary deformation called tank-treading that provokes a drop-like behavior [37]. The RBC deformation in the core region further reduces the viscosity of blood.
A direct consequence of the cell-free layer is the so-called Fåhraeus effect, which was discovered in 1929: The majority of RBCs in the cross section of a narrow blood vessel are located near the vessel core. Obviously, the average of the velocity in this RBC-rich core is significantly higher than the average plasma velocity along the vessel cross section, especially because a larger fraction of plasma is located near the vessel walls as compared to a homogeneous solution. Thus, in average in a narrow vessel RBCs travel faster than the suspending plasma. Consequently, the average volume fraction of RBCs in a cross section, namely the true Hematocrit, must be smaller than the Hematocrit of the whole blood [12]. This variable relation between true tube Hematocrit and feed Hematocrit is referred to as the Fåhraeus effect and was shown to decrease with decreasing vessel radius and with decreasing feed hematocrit [2].

Micro-continuum approach
In this study, the micromorphic theory is used as a continuum model for blood. This theory was introduced by Eringen and Suhubi for solid materials and adjusted by Eringen to model fluent media [6,8]. The approach aims at describing media with a deformable substructure by considering the actual rotation and deformation of substructure on a continuum level. Note that there are other continuum mechanical approaches applied to non-Newtonian fluid problems, where the non-Newtonian behavior is modeled by identifying mathematical instead of phenomenological relations governing the material behavior. For instance, the gradient theory from [34] considers a stress tensor containing higher gradients of the strain-rate, whereas the stress tensor in the micromorphic approach contains the micromotion instead, as shown below.
For a micromorphic material, each material point is considered the center of mass of a deformable particle of infinitesimal size. The position of a material point is described by the vector X or x in the reference or current configuration, respectively. Analogously, on the micro-level, the infinitesimal vectors Ξ and ξ describe the position of a material point within a particle with respect to its center of mass in the reference and current configuration, respectively. The inverse motion and inverse micromotion are required to exist and to be unique, because two material points cannot be located at the same spatial point simultaneously. In order to relate differential line elements in the two configurations, the deformation tensor F and the micromotion tensor χ are introduced [8]: As indicated by the lower and upper case index letters, F and χ depend on both the current and reference configuration and are two-point tensors. Similar to the macro-scale deformation tensor, the micromotion tensor describes the particle deformation on the micro-scale. Thereby, it approximately links the microscopic position vectors in the two configurations via [8]: Furthermore, the determinants of both F and χ are assumed to be positive, in order to 'retain the right-hand screw orientations of the frames-of-references,' [10, p. 6]. This is sufficient for the inverse tensors F −1 and χ −1 to exist. The Eulerian velocity gradient L and gyration tensor G are defined as [8]: where v is the velocity field, ∇ denotes the del or nabla operator, and a superscript dot indicates the Lagrangian or material derivative:φ The micromotion induces a particle deformation on the micro-scale which is accompanied by a change of the mass distribution, characterized by the microinertia tensor i [8]: Here, and dv denote the density and an infinitesimal volume element in the current configuration, respectively, while the superposed prime indicates the corresponding quantity on the micro-scale. Definition (1.8) slightly deviates from the definition of the inertia tensor commonly used in classical mechanics [17]. However, the two formulations can be transformed into each other, as described in [36]. By the use of the inertia tensor i and gyration tensor G, the spin inertia per unit mass σ [6] is defined: The spin inertia itself is not a time-rate, though it contains the material derivative of the gyration tensor. It characterizes the change in time of a particle's spin angular momentum (cf. [36]). The field quantities introduced in this section are used to formulate the governing balance and constitutive equations of a micromorphic fluid below.

Micromorphic fluid theory
This section comprises the governing equations to describe the motion of a micromorphic fluid in the current configuration. Apart from fields of velocity, temperature, and pressure known from classical fluid mechanics, micromorphic theory additionally depends on the gyration tensor G which consists of nine independent components as introduced above. For a better interpretation, this tensor can be decomposed into its symmetric part sym(G) representing the microshearing and stretching, and its skew-symmetric part skw(G) representing the microrotation. The evolution of these fields is also explained by the balance equations stated below.

Balance equations
The balance equations are presented in their local form since they are required to hold in every arbitrary infinitesimal control volume dv within the domain Ω. On the boundary ∂Ω, certain jump conditions resulting from boundary integrals in the balances' global form must hold, if no Dirichlet boundary conditions are specified. These will not be stated in this article, but further information on them can be found in [6]. In addition, all balance equations stated in this subsection can also be found in [6]. First, the usual local form of the balance of mass is given by The momentum equation reads Here, f and t are the external body force and stress tensor, respectively. In contrary to classical continuum mechanics, the stress tensor is usually not symmetric for microfluids. The micro-motion is governed by the so-called micro-spin balance, and depends on the symmetric micro-stress average tensor s, the first moments stress tensor of third-order 3 λ, and the external moments tensor l. Here, the micro-stress average tensor corresponds to the averaged stress tensor of a point particle. It can be obtained from the momentum equation within a particle, see [8]. The balance of angular momentum of a point particle is fulfilled, if the corresponding micro-stress is symmetric, because a point particle itself does not contain a microstructure. Since the spin inertia σ contains the material derivative of G, Eq. (2.3) corresponds to a balance equation describing the micromotion's change in time. The microinertia is assumed to be conserved, which is formalized by the following identity: Furthermore, the balance of internal energy is given by where ε, q and h denote the specific internal energy density, the heat vector field, and a heat source, respectively. Moreover, the double dot operator " ·· " is defined through Finally, the Clausius-Duhem inequality in differential representation assures that the second law of thermodynamics is not violated [10]. Here, η and θ denote the specific entropy and the absolute temperature, respectively, while ψ = − θη represents the Helmholtz free energy. Since in this article only non-heat-conducting materials are considered, the energy balance will not be solved.

Constitutive equations
In order to obtain a set of Partial Differential Equations (PDEs) governing the microfluid fields, constitutive equations need to be employed. In this article, the linear constitutive equations from [6] are used if not stated otherwise. First, the non-symmetric stress tensor links the stress to the velocity gradient, the gyration tensor, and the thermodynamic pressure π, and contains the viscosity coefficients λ v , λ 0 , μ v , μ 0 , and μ 1 as material parameters. Moreover, the velocity gradient was decomposed into its symmetric part, i.e., the strain-rate tensor d = sym(v ⊗ ∇), and its skew-symmetric part, namely the vorticity tensor w = skw(v ⊗ ∇). The average micro-stress tensor depends on the velocity gradient and on the gyration tensor as well as on the material parameters η v , η 0 and ζ 2 . By equating these parameters with the corresponding ones in the stress tensor, i.e., η v = λ v and η 0 = λ 0 , Fig. 1 Underlying domain the governing equations reduce to the Navier-Stokes equations in the absence of micromotion, i.e., if G = w [6]. Finally, the first stress moments tensor introduces 15 new gyration viscosity coefficients γ i and, furthermore, depends on the gyration tensor, the metric tensor g, and the transposition tensor 4 T:=e i ⊗ e j ⊗ e j ⊗ e i . The first stress moments tensor represents the "surface tractions acting on the [particle] surfaces," [8, p. 194], and thereby affects the stress on the macro-scale via Eq. (2.3). It should be mentioned that in more recent works, a different representation of the constitutive equations is chosen. For example in [11], the first moments stress tensor only depends on 11 viscosity coefficients but is also independent of the metric tensor. Usually, the viscosity coefficients are functions of the temperature, but due to the neglect of non-heat-conductive media, they can be treated as material constants.

Blood as a micromorphic in idealized narrow blood vessels
In this section, blood flow in a cylindrical pipe, which can be regarded as an idealized blood vessel is modeled as a micromorphic fluid. In order to assess how suitable the micromorphic model is, the hemodynamic properties of a micromorphic fluid are analyzed. This is achieved by deriving the solution to the laminar flow problem. Precisely, the solution of the Boundary Value Problem (BVP) is found for the laminar incompressible stationary case. Such a solution is not only valuable to determine restrictions on material parameters for blood but also serves as a reference solution for the validation of numerical solutions. A solution of the problem was already presented in [22], but it contains minor errors which also affect the final results. Therefore, the correct solution is presented and also validated in the following.

Assumptions and boundary conditions
The domain corresponds to a circular pipe with radius R, and the problem will therefore be solved in cylindrical coordinates. Thus, the domain can be expressed as A sketch of this domain is shown in Fig. 1. Furthermore, the assumption of rotational symmetry allows to solve the problem in the r -z-plane: Accordingly, the former physical boundary at the pipe wall Γ 1 :={[r, z] ∈ R 2 | r = R} and the cylinder axis Γ 2 :={[r, z] ∈ R 2 | r = 0} are defined, and make up the boundary of the computational domain: Another consequence of rotational symmetry is that the microrotation tensor contains only two components spanned by e r and e z , and therefore, the micro-rotation reads Analogously, the micro-shearing only occurs due to motion in the r -z-plane. This means where the migration velocities G rr and G zz were neglected, because they are usually small in crowded suspensions [22]. In addition, the fluid is assumed to flow parallel to the vessel wall, while the fields v and G only change along r : (3.3b) Therefore, the microinertia i is also a function of r . Because of the assumption of incompressibility, the conservation of mass (Eq. (2.1)) reduces to the incompressibility condition: v · ∇ = 0.
Accordingly, the micro-incompressibility is assumed: Furthermore, the flow is assumed to be stationary, such that all partial derivatives with respect to time vanish. As a consequence of the assumptions, the convective term also vanishes: Finally, external body forces or moments of any kind are not considered. The constitutive equations (2.2) can be substituted in the balance equations (2.1), which allows to simplify the BVP, as derived in the following.

Boundary conditions
The boundary condition for the velocity field corresponds to the typical homogeneous Dirichlet or no-slip boundary condition on Γ 1 : For the gyration tensor, on the other hand, the postulation of boundary conditions is not that simple: Multiple different boundary conditions, such as (in-)homogeneous Dirichlet boundary conditions, see [7], or homogeneous Neumann boundary conditions for the micro-rotation are possible, see also [1,22]. However, here a coupling boundary condition is applied. This is motivated by the fact that the velocity gradient is the only field having a direct effect on the suspended particles. The vorticity will cause the particles to rotate near the wall, while the shearing of the fluid will induce a shearing and stretching of the particles. Formally, this effect is expressed as [22] The coupling parameters λ 2 and λ 3 can be obtained from a comparison with the linear and angular velocity of a rigid, spherical particle in a dilute suspension, respectively: This particle is assumed to roll and slip along the wall, depending on the concentration of the suspension, i.e., the hematocrit. The linear and angular velocity were derived using perturbation methods in [22]. For a crowded suspension of deformable particles, these findings are applied by relating the microshearing and microstretching to the derived linear particle velocity and adjusting the derived particle rotation by a crowding factor (1 − Ht) to account for frictional effects [22]: The remaining coupling parameter λ 1 is irrelevant in this case due to the incompressibility condition, but it is linked to the migration velocity and assumed to vanish anyway [22]. Finally, all fields are required to be regular on Ω, which will result in mathematical boundary conditions on Γ 2 , as shown during the derivation of the solution.

Constitutive equations and viscosity coefficients
Under the assumptions stated in Sect. 3.1, the constitutive equations simplify to [22]: where the new viscosity coefficients μ, ζ 1 , ζ 2 consist of combinations of viscosity coefficients introduced in Sect. 2.2, which are given by [22] μ : Substituting Eqs. (3.7) in the Clausius-Duhem inequality (2.7) yields multiple restrictions on the viscosity coefficients, of which the most important ones are given by [22, p. 145]: Furthermore, in the context of blood, several conclusions regarding the magnitudes of the viscosity coefficients and trends with respect to the initial hematocrit were made in [22]. They originate from a comparison with the stress tensor for dilute suspensions [14]. Since blood is not considered a dilute but a dense suspension, they cannot be adopted one-to-one. The findings do not include how crowded particles affect each other, but the remaining findings are still assumed to hold: μ increases with increasing hematocrit.
-|ζ 1 | and |ζ 2 | also increase with increasing hematocrit, but at a higher rate than μ, because they govern the influence of deformability on the resulting stress. Hence, larger concentration yields stronger interactions between particles. This effect will, for a dense suspension, exceed the effect of particles influencing each other via the surrounding fluid. -κ, on the other hand, represents the stress due to particle rotation, and thus will also increase at a higher rate than μ. It is of the same order of magnitude as ζ 2 , but cannot have a larger order of magnitude than μ. -The rotational viscosity coefficients γ i correspond, if normalized by μ, to a characteristic length of the particles and are therefore chosen such that √ γ i /μ matches the size of particles, i.e., around e-6m to e-5m.

Solution to the micromorphic pipe flow
The steady flow of a micromorphic fluid as described above is summarized by the following BVP in tensor representation: (3.10) The final solutions to this BVP are stated here, while a comprehensive derivation is given in "Appendix A." The solutions are stated as functions of the normalized radial coordinate ρ:=r/R, while the solutions themselves are normalized by The expression for the velocity field reads (3.11) and the micromotion is given by where G (i j) and G [i j] denote the components of sym(G) and skw(G), respectively. The constants C 1 and C 2 within the solutions are: Here, α and β are the positive solutions of the following equation: The α 2 i and K i are composed of viscosity coefficients as follows: 2μ + κ)). In the following, the obtained solutions are presented and discussed for some special cases. Here, Eqs. (3.12) and (3.11) correspond to Eqs. (5.8) from [22, p. 149], although, the constants C 1 and C 2 deviate from these used in [22], which will also be investigated below.

Results and discussion
By means of the solution obtained above, the applicability of a micromorphic fluid as a model for human blood can be discussed. The analysis focuses on pipes with small radius, which represent simplified blood vessels. First, the relationship between experimental data and analytical solutions is discussed to facilitate a meaningful discussion of the results. Second, the results are compared to corresponding results from the literature (cf. [22]). Third, a new set of viscosity coefficients is derived from the Fåhraeus-Lindqvist effect. Then, the resulting fluid model is compared to experimental data and, finally, the influence of several viscosity coefficients on the micromorphic fluid behavior is investigated.
The viscosity coefficients used in this section generally depend on the initial RBC hematocrit. A motivation for restrictions regarding their magnitude is given in [22] and briefly repeated in Sect. 3.3. Furthermore, the values for the coupling boundary condition in this section are specified in Eq. (3.6). Thus, the micromotion of the fluid at the wall is considered to be driven by the surrounding fluid motion, but also slowed down by particle-wall interaction and crowding effects.

Normalization of experimental data
As explained in Sect. 1, the viscosity μ of blood is not constant with respect to the vessel radius R, but decreases with decreasing radius for R < 300 µm, which corresponds to the well-known Fåhraeus-Lindqvist effect (cf. [13]). In order to conduct a meaningful discussion, it is necessary to clarify how the fields are normalized. To achieve this, the apparent viscosity μ app is reformulated in terms of the micromorphic viscosity coefficient μ from Eqs. (3.8) [22]: where Q denotes the total volume flow rate, given by By using ρ = r/R and dr = R dρ, the volume flow rate can be written in terms of the dimensionless velocity from Eq. (3.11): The remaining integral is denoted as Q, which yields an adapted, dimensionless expression for the volume flow rate [22, p. 152]: For a Newtonian fluid, the apparent viscosity is constant, because the Newtonian volume flow rate Q is always proportional to P R 4 . Consequently, Q = 1/4 will always hold. Hence, the apparent viscosity can be interpreted as a measure for the change of the normalized volume flow rate for a varying radius. A variable apparent viscosity indicates a deviation from Newtonian behavior. To normalize the solution properly, the micromorphic viscosity coefficient μ can be set to a reference viscosity for all R. For instance, a suitable reference viscosity is the asymptotic value of approximately 3.77mPa s obtained from experiments with large radii by [20]. By choosing v 0 as a function of a constant μ, the normalized volume flow rate Q for a micromorphic fluid will change with the vessel radius, because the normalized solution v(ρ)/v 0 still depends on R. The choice of a reference velocity ofṽ 0 = −P R 2 /(4μ app ) with an experimentally obtained value of μ app corresponds to a normalization by the experimental volumetrically averaged velocity v m , see Eqs. (4.1, 4.4): If the fluid model shows the same shear-thinning behavior as blood, the normalized volume flow rate will stay constant with varying radius. Thus, if the subsequently computed volume flow rate of v(ρ)/ṽ 0 does change with R or dimensionless experimental data does not fit the model, it is an indication for a flawed model or wrongly chosen parameters.

Comparison to Kang and Eringen
The solution derived above is compared to the reference solution from Kang and Eringen [22]. In order to enable this comparison, all model parameters are matching those from the reference. The experimentally obtained apparent viscosity of blood specified in [19,21] for R = 100 µm was chosen for the micromorphic viscosity coefficient μ to compare velocity profiles at different R with each other and with experimental data. As explained above, this is a rather arbitrary and counter-intuitive choice for the viscosity, since it neither corresponds to the asymptotic nor the actual apparent viscosity at R ∈ {20 µm, 35 µm}. Due to the normalization with v 0 and thereby with μ, the choice of μ does not affect the velocity and micromotion profiles. Therefore, the values from [22] are adopted for now which includes the following viscosity coefficients:  There are two sign errors in the boundary values given in [22]. The first one was already outlined in "Appendix A" in Eq. (A.21): According to the coupling boundary from Eq. (3.5), the boundary condition should read This formulation was also used for the derivation of our solution. However, after substituting the coefficients λ 2 and λ 3 from [22, p. 151] into Eqs. (5.6) from [22], the boundary conditions read: Comparison of Eqs. (4.7) and (4.8) reveals a sign error within the boundary condition for the microrotation in [22]. This sign error would result in a negative relation between the microrotation and the vorticity at the wall. Hence, the micro-rotation at the wall would not only be smaller than the corresponding vorticity due to friction but also in the opposite direction. This behavior is not physically plausible and the formulation from [22] is therefore considered to be faulty. The solutions to the BVP in Eq. (3.10) are exemplarily compared for R = 20 µm. This comparison includes our analytical solution, the solutions published in [22], and a numerical solution. Before correcting the flawed boundary condition as described above, all solutions are compared for the physically implausible boundary condition stated by Kang and Eringen. It will be shown that our derived analytical solution does not coincide with the analytical solution from Kang and Eringen. Therefore, a numerical solution is used to validate our derived solution. This numerical solution was computed using the NDSolve-function of the Mathematica software, in which the three scalar PDEs along with the boundary conditions were specified. Next to the mathematical expression of the analytical solution, the velocity and the microrotation field were also plotted in by Kang and Eringen in Fig. 3, see [22, p. 151], which are also compared to the analytical solutions. The fact that these plotted solutions deviate from both analytical solutions indicates misprints or errors in [22]. Accordingly, the coupling parameters λ 2 and λ 3 will be corrected such that the physically meaningful boundary condition is prescribed, and the same comparison is performed again. After the adjustment, our analytical solution for the microrotation coincides with the plotted microrotation solution from Kang and Eringen, whereas their corresponding solution in mathematical form differs. For the velocity solution on the other hand, our validated analytical solution does neither match the graphical nor the analytical solution published in [22]. This leads to the conclusion that the analytical solutions from Kang and Eringen contain, first, misprints because their analytical solutions do not match their plotted solutions, and second, errors, since both their graphical and analytical solution for the velocity field differ from the correct solution to the BVP. The derivation of this conclusion is given in more detail below.
For the original boundary conditions from [22], the different velocity profiles are shown in Fig. 2, and the micromotion is displayed in Fig. 3. It is noticeable that the numerical results for all fields coincide with our derived solutions, whereas the analytical solution from Kang and Eringen deviates from them and-in the case of velocity and microrotation-from their plotted solution. For the microshearing, no graphical solution was provided by Kang and Eringen. Hence, our analytical solution solves the BVP, unlike the solutions from Kang and Eringen. The small difference between the analytical and graphical solution for the velocity field in [22] could originate from inaccuracies in plotting, and thus, no statement regarding the source of the difference to our solution can be made yet. In order to determine if a misprint in the boundary condition contributed to these differences, the flawed boundary condition is corrected and the solutions are compared again.
After correcting the boundary condition, differences between our analytical solutions and those from Kang and Eringen still exist, as depicted in Figs. 4 and 5. However, for the microrotation in Fig. 5b, the situation is different: The correction of the boundary conditions causes our solution and the plotted solution from Kang and Eringen to match, while their analytical solution still differs. Minor differences between our solution and Kang and Eringen's graphical solution originate from the extraction of the graphical solution. It should also be noted that the plotted solution from [22] was negative for all ρ, whereas it should be positive for all ρ, at least when considering the correct boundary condition. For Fig. 5b, this error was treated by mirroring the graphical solution from [22]. The good accordance of our analytical with the graphical solution from Kang and Eringen indicates that the corrected coupling boundary condition indeed corresponds to the boundary condition they intended to use and that their analytical solution for the microrotation simply contains typographical errors. Regarding the velocity field, neither the mathematical expression nor the plot from [22] agrees with our validated solution. Thus, there is another error in the velocity solution that cannot be explained by a typographical error, although it could be due to the typographical error in the microrotation solution. For the microshearing it is similar: Albeit small, there is a difference between our validated analytical solution and the one from [22]. Consequently, the latter also contains errors. To summarize, the solutions to the micromorphic pipe flow from Kang and Eringen contain several errors, as already stated. One error regarded the misprints in the solution for the microrotation, while the exact origin of the remaining errors could not be identified. Moreover, the analysis showed that the flawed boundary condition from [22] merely corresponds to a typographical error in the publication and that the physically meaningful boundary condition from Eqs. (4.7) was also considered by Kang and Eringen.

Derivation of model parameters
Due to the counter-intuitive choice of μ and the flawed equations in [22] which might have led to the viscosity coefficients from Eqs. (4.6), a new set of model parameters is derived based on experimental data. To acquire this new set of parameters, the apparent viscosity of a micromorphic fluid is fitted to the actual apparent viscosity of human blood with an initial hematocrit of 45%. The reference data were collected and aggregated in [27]. It comprises the results of 18 different studies of RBCs suspended in a Newtonian fluid with hematocrit values from 40% to 45%, where results with hematocrit of less than 45% were extrapolated. Since the suspending fluid differed across the considered studies, an average for the relative apparent viscosity μ rel , i.e., the apparent viscosity divided by the viscosity of the suspending fluid μ s , was determined. From hereinafter the micromorphic viscosity μ denotes the asymptotic apparent viscosity as a function of μ s under consideration of the results from [27]: The fitting was conducted with the help of the curve_fit function from the python package scipy. optimize. This corresponds to applying the non-linear least-squares analysis to find a suitable set of parameters. The analysis was conducted for 19 data points in the range of radius from approximately 0.01mm to 2mm. For smaller vessel radii, the vessel dimension approximates the dimension of single RBCs, and the continuum approach ceases to be valid. The following bounds were imposed on the parameters according to the heuristics in [22]:   [27] and [19] The fitting resulted in the following set of parameters: (4.10) It should be noted that the values obtained for the aggregated gyration viscosity coefficients K i as well as the coefficient ζ 2 approximately lie on their respective bounds. This indicates that the optimal values for K i and ζ 2 lie beyond the imposed bounds. Nonetheless, for the following analysis the values from Eqs. 4.10 are considered.

Shear-thinning behavior of a micromorphic fluid
In this section, micromorphic fluid behavior is compared to the Fåhraeus-Lindqvist effect. The asymptotic relative viscosity is computed for R ∈ [10 µm, 1000 µm] and plotted against the experimental results from Pries et al. [27] and Haynes [19]. The former was used to determine the material parameters, while the latter served as reference data in the article from Kang and Eringen [22]. Since this analysis considers the relative instead of the absolute apparent viscosity, the results from [19] are normalized by the viscosity of the suspending fluid, i.e., Acid-Citrate-Dextrose (ACD). The different experimental results are not comparable one-to-one though, because the experiments in [19] were conducted for Ht ∈ {0%, 20%, 30%, 40%, 50%, 60%, 70%, 80%}. Therefore, the results were interpolated to 45% quadratically. The resulting apparent viscosity is shown in Fig. 6. Additionally, the viscosity of two experiments conducted by Bugliarello and Sevilla [4] is included in the plot which are discussed in Sect. 4.5. At first glance, it stands out how well the micromorphic relative apparent viscosity matches the actual relative apparent viscosity of RBC suspensions, as compiled in [27]. Clearly, this is due to the experimental data being used for the derivation of the model parameters, but the ability to reproduce the Fåhraeus-Lindqvist effect at such accuracy is remarkable. It indicates that, for sufficiently high shear rates, a micromorphic fluid with the parameters from Eqs. (4.10), will exhibit nearly the same flow resistance as human blood in glass capillaries. However, for the micromorphic fluid to be a viable model for blood in microcirculation, it is not only required to follow the Fåhraeus-Lindqvist effect, but also to resemble the actual motion of human blood. Consequently, the velocity distribution of a micromorphic fluid is investigated in the following section. Another interesting observation concerns the experimental results from Haynes [19]: At Ht = 40%, they coincide well with the data from Pries et al. [27], whereas they exceed the averaged values if interpolated to Ht = 45%. This strong variation across experimental data from literature was also pointed out in [27].

Comparison to experimental velocity profiles
In this section, several experimentally obtained velocity profiles of RBC suspensions in narrow glass capillaries are compared to those of a micromorphic fluid. This concerns data from Bugliarello and Sevilla [4] and Goldsmith and Marlow [15].
For R ∈ {20 µm, 35 µm}, a comparison with experimental data from Bugliarello and Sevilla [4] is possible. The experiments were conducted on a pressure-driven flow of RBCs suspended in ACD at Ht = 40%. Along with the absolute velocity, pressure-flow-diagrams were published, where the volume flow rate can be computed from the accompanying mean velocity v m . This allows to determine v 0 which, in turn, allows for a normalization of given data. Since the published v m does not agree with the respective velocity measurements, it is adjusted for the computation of v 0 (see "Appendix B"). Furthermore, the relative apparent viscosity was computed for both experiments using the ACD viscosity measured in [19] (see Fig. 6). Apparently, the relative apparent viscosity from the measured blood flows exceed the asymptotic relative apparent viscosity of human blood. Apart from noise, a possible explanation for the deviation is that the pseudo-shear rate v from Eq. (1.3) is too small for the apparent viscosity to tend to its asymptotic limit. Especially, for R = 35 µm with v ≈ 17 s −1 , this is a plausible explanation. For R = 20 µm with v > 100 s −1 on the other hand, the deviation should be negligible. In both cases, the normalization term v 0 was additionally shifted by the ratio of the different apparent viscosities. These two adjustments combined basically correspond to a normalization with respect to the volume flow rate. The resulting velocity profiles are shown in Fig. 7a, 7b.
Although the parameter optimization was only executed with respect to the apparent viscosity, the micromorphic velocity fields closely resemble the velocity fields of blood. This is especially noticeable if compared to a Newtonian velocity profile with the same volume flow rate: The micromorphic fluid exhibits the relatively blunted shape toward the vessel center, comparable to crowded RBCs forming a core region with smaller velocity differences. Near the vessel wall on the other hand, the velocity increases rapidly, similar to a friction-poor plasma layer.
Similar results are obtained from a comparison with experimental data from Goldsmith and Marlow [15], where the velocity fields of RBC and ghost cell suspensions were measured at R ≈ 40 µm and varying hematocrit. Since in this case only a volume flow rate but no pressure gradient was given, the data are normalized to meet the volume flow rate of the micromorphic fluid. The resulting comparison is displayed in Fig. 8. Although the pseudo-shear rate v is below 50 s −1 for all considered data sets, the general, blunted shape of the micromorphic and RBC suspension velocity profile is quite similar. Especially for an increasing pseudo-shear rate, the experimental data tends toward the micromorphic solution. However, it should also be noted that the volume flow rate of such low-shear experiments will most likely not agree with the micromorphic volume flow rate.
In conclusion, not only does the micromorphic fluid with adequately chosen parameters meet the Fåhraeus-Lindqvist effect, but experimental velocity profiles can also be reproduced satisfactorily in a circular pipe flow. As to why micromorphic dynamics coincide with hemodynamics, and which parameters influence this behavior, will be investigated below.  Using the derived material parameters from Eqs. (4.10), the micromorphic fluid motion resembles the motion of blood very well. Now, in order to get a better understanding of the associated material parameters, this section covers the influence of the viscosity coefficients contained in the stress tensor t, namely κ and ζ 1 , on the (micro-)motion and on the accompanying shear-thinning behavior. Of course, the viscosity μ is also included in t, but a variation of μ will only affect the apparent viscosity, while it does not implicate a change within the normalized (micro-)motion.
For this study, the difference between vorticity and microrotation as well as the microshearing are of particular interest: Microrotations differing strongly from the surrounding fluid vorticity produce stress causing the fluid resistance to increase, while the microshearing and microstretching along the fluid shearing, i.e., the strain rate, causes stress reduction within the system. Therefore, the viscosity coefficients are varied around their optimized values, and the resulting change in relative apparent viscosity μ rel , velocity v, microshearing sym(G), strain-rate d, and relative microrotation Ω is depicted. Here, the relative microrotation is defined as Ω:=skw(G) − w. The parameter study regarding the field quantities is carried out at R = 35 µm. For a variation of ζ 1 , the micromotion results are presented in Fig. 9, where the microshearing and relative micromotion were normalized by the strain-rate and vorticity, respectively.
Varying ζ 1 has an insignificant effect on the normalized microshearing: Irrespective of ζ 1 , it nearly vanishes in the RBC-rich core region, while stronger microshearing is induced in the high-shear area near the vessel wall. Thus, the ratio of strain-rate to micro-shearing is mostly independent of ζ 1 . In contrast, the normalized relative microrotation is strongly influenced by a variation of ζ 1 in the core region: In relation to the fluid vorticity, decreasing ζ 1 also yields a decreasing relative microrotation. As for the microshearing, a variation Fig. 10 Influence of ζ 1 on the relative microshearing, relative microrotation, and velocity field of a micromorphic fluid in the microcirculation Fig. 11 Influence of κ on the microshearing and relative microrotation field of a micromorphic fluid in the microcirculation of ζ 1 has a diminishing effect in the high-shear area close to the vessel wall. The concomitant change on the macro-scale fields is displayed in Fig. 10.
The strain-rate is simply increased in magnitude near the vessel wall for decreasing ζ 1 , whereas the strainrate in the inner vessel is only influenced marginally. Mathematically speaking, this is a consequence of ζ 1 directly influencing the Neumann boundary condition imposed on the velocity field. A possible interpretation can be inferred from the stress tensor t: ζ 1 governs the stress due to shearing on the micro-scale. In fact, Kang and Eringen state that higher |ζ 1 |, |ζ 2 | correspond to more stress exerted on the system due to deformation, and thus, a "more flexible substructure," [22, p. 146]. However, since t is linear in ζ 1 , changing the sign of ζ 1 also reverses the effect. For ζ 1 > 0, the stress in the system due to microshearing along the strain rate is increased, while ζ 1 < 0 induces stress relief as a consequence of microshearing along the strain rate. Hence, decreasing ζ 1 leads to a more flexible substructure. As a consequence, decreasing ζ 1 also leads to an increased velocity along the entire vessel. Figure 11 shows how κ influences the micromotion: The variation of κ barely affects the normalized microshearing, while increasing κ results in a lower normalized relative microrotation, especially in the low-shear area near the vessel center. In the high shear flow near the vessel wall on the other hand, the microrotation-tovorticity-ratio is mostly indifferent to a variation of κ due to the imposed boundary condition. The corresponding changes on the macroscopic field quantities are shown in Fig. 12. As a consequence of an increase in κ, the strain-rate declines, and the velocity diminishes. As visible in the formulation of the stress tensor t and as already stated in [22], large κ induce a high flow resistance due to particles rotating at a lower angular velocity than the surrounding vorticity field.
These phenomena can also be observed for the relative apparent viscosity, as shown in Fig. 13. An increase of ζ 1 leads to a less pronounced shear-thinning and even shear-thickening for ζ 1 > 0. In a similar manner, higher values for κ damp the shear-thinning.
In conclusion, especially the negative material parameter ζ 1 governs the micromorphic shear-thinning behavior through stress reduction induced by particle deformation, whereas κ controls the stress due to differences between particle rotation and fluid vorticity. Hence, the micromorphic shear-thinning behavior is   [22]. Particle flexibility also corresponds to the actual cause for shear-thinning in blood [16]. For the given configuration, both stress reduction and stress build-up due to micromotion are mostly observed near the vessel wall. Moreover, the micromotion is only activated in areas with sufficiently high strain-rates, where the ratio of micromotion to strain-rate is always similar, irrespective of the chosen κ and ζ 1 . This displays a major difference between micromorphic fluid behavior and hemodynamics: In actual blood, there is no particle deformation close to the wall, because there are no RBCs located close to the wall. However, Kang and Eringen motivate this discrepancy by an inverse relation between Ω and the RBC-concentration: High magnitudes of Ω, i.e., a strong deviation of vorticity and microrotation, indicate that a particle has a greater sphere of influence, and thus, the local hematocrit must be lower.

Influence of varying vessel diameter
Regarding the motion on a macroscopic scale, the most distinguishing property of vertebrate blood compared to a Newtonian fluid is its shear-thinning behavior, i.e., the decrease of apparent viscosity μ app along increasing strain-rate. As described in Sect. 1.1, the apparent viscosity of blood can be defined as a function of R, which, in turn, is satisfactorily modeled by a micromorphic fluid with our derived viscosity coefficients. This section is dedicated to analyzing the corresponding motion of a micromorphic fluid for varying vessel radii R. The considered vessel radii range from 10 µm to 500 µm to cover the full scope of microcirculation. In the former sections, it has already been shown that the velocity profile takes a more blunted shape than a Newtonian fluid. This evolution along varying R can be observed in Fig. 14 for the micromorphic fluid, where it was again put in context with a Newtonian fluid of the same normalized volume flow rate.
For large vessel radii, the micromorphic fluid almost fits the Newtonian fluid as expected, since here the inertial effects dominate the viscous effects of the material. From there, the micromorphic fluid profiles become more blunted until R = 35 µm, from where the trend retrieves until the velocity profile is almost parabolic again at 10 µm. It is not clear if this corresponds to the actual behavior of blood, or if the continuum approach ceases to be valid. More experimental data are required to validate the micromorphic approach as a model for blood in small vessels. In addition to the shape change, the normalized velocity profile also becomes steeper if the vessel diameter is decreased. The normalized volume flow rate Q is thereby increased. As a result, the volume flow rate Q is no longer proportional to P R 4 , as it would be for a Newtonian fluid. This behavior corresponds to shear-thinning, or in the context of blood to the Fåhraeus-Lindqvist effect. Since the origin of the micromorphic shear-thinning can be attributed to the micromotion, the evolution of the microshearing and the relative microrotation along varying vessel diameter is investigated below, see Fig. 15.
For decreasing vessel diameter, both the microshearing G (r z) and relative microrotation Ω r z tend to a linear function. For large vessel radii beyond the scope of microcirculation on the other side, both G (r z) and Ω r z vanish in the core region, but are increased in magnitude near the vessel wall abruptly. This agrees with the fact that blood can be considered a Newtonian fluid in large arteries. It also corresponds to the mechanism of the RBC-free, plasmatic layer in blood flow, whose thickness does not scale with the vessel radius [19,31]. The relative plasmatic layer thickness decreases along increasing vessel diameter, as does the micromotion's area of influence in Fig. 15. This behavior is also reflected in the strain-rate, as displayed in Fig. 16. For all radii but R = 10 µm, the micromorphic fluid can be divided into a low-viscosity and high-viscosity area, indicated by the the steep strain-rate profile at the wall and the rather shallow profile in the core region, respectively. Just like the plasmatic layer, the relative thickness of the low-viscosity part shrinks with increasing vessel diameter.

Conclusion
While the calculations made in [22] partly contained misprints and the analysis was methodically flawed, the conclusions reached regarding the velocity profile and apparent viscosity are still valid: Experimental data is reproducible, shear-thinning occurs and even follows the Fåhraeus-Lindqvist effect, while the reason for that, namely particle deformability, is even in accordance with the corresponding physical interpretation. Therefore, the so far studied effects indicate a high viability of a micromorphic fluid as a model for blood. However, there are more effects observed in blood, such as the evolution of the plasmatic layer size or the Fåhraeus effect. They depend on the concentration, which is not a field quantity in micromorphic fluid theory. In order to study these effects anyway, Kang and Eringen derived a formulation for the RBC-concentration as a field quantity. This formulation is based on energy dissipation due to particle rotation and fluid shear in the presence of particles (cf. [18,35]). Since the particles were considered to be rigid and spherical, this analysis is not error-free applicable to a fluid model that considers a deformable microstructure. Therefore, it is not included here. But the approach could be used to compute the concentration of a microfluid with a microstructure of rigid, spherical particles, such as Eringen's micropolar fluid [9].
Despite the good results, the theory of micromorphic fluids is complex in that it introduces nine additional kinematic degrees of freedom and several additional viscosity coefficients compared to a Newtonian fluid. Hence, it is mathematically viable, and the linear stress tensor is a great advantage toward many popular models for blood depending on a strain-rate measureγ = √ 2d ·· d, see [30] for a comprehensive review. Such models including a nonlinear stress tensor easily result in numerical complications and restrictions. However, simulating a micromorphic fluid numerically can also result in a large model, requiring a lot of processing power and memory. In addition, the determination of more than 20 material parameters is a complicated task that requires a lot of reference data. Especially for blood, acquiring a large amount of data proves to be difficult, because the suspended particles are not translucent. Thus, a first step would be the model reduction to identify and discard redundant and unnecessary terms, in particular with regard to the first stress moments tensor, which alone contains 15 gyration viscosity coefficients. Unfortunately, this task again requires a lot of experimental data. Therefore, a less complex model, such as a modified micropolar fluid model, is desired to accurately compute blood flow in micro-circulation.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Funding
Here the thermodynamic pressure was replaced by the hydrostatic pressure because of − p1 = t • , where the superscript circle denotes the spherical part of t, i.e., the thermodynamic pressure. In addition to the ODE in Eqs. (A.1c), the second equation reveals information about the pressure: Hence, the ODE in Eq. (A.1c) is solved by For the sake of readability, the pressure gradient in z will hereinafter be referred to as P:=∂ p/∂z. Substituting the r z-component of the stress tensor t, as displayed in Eq. (3.7a), in the ODE, yields where the constant C 1 must vanish due to the regularity requirement at r = 0. The obtained equation corresponds to Eq. (5.3a) from [22, p. 148]. Rewriting the expression leads to a function describing the velocity gradient in terms of the micro-rotation, the micro-stretching and the pressure gradient: where the different K i consists of the gyration viscosity coefficients while the α 2 i are composed of the summarized viscosity coefficients: The expressions for α 2 i are adopted from [22], and Eqs. (A.8) correspond to Eqs. (5.3b)-(5.3c) in [22, p. 148]. In the following, these derived equations are used for finding a solution for the gyration tensor as well as for the velocity field. This section concerns the solution of the above-derived BVP. The derived and partially decoupled field equations are solved below, starting with the solution for G, which yields the velocity solution after a simple integration.

A.3 Solution for the micro-motion
In order to solve Eqs. (A.8), the homogeneous problem is considered first, and the inhomogeneous problem is solved afterwards to obtain a general solution. Finally, the boundary conditions are applied to find the specific solution.

A.3.1 Homogeneous solution
The different expressions in Eqs. (A.8) are similar to the modified Bessel differential equation which is solved by the modified Bessel functions of the first and second kind. Therefore, the following ansatz is made: Here, I n (x) and K n (x) denote the Bessel functions of the first and second kind, respectively. Since K n (αr ) is not defined for r = 0, and the solution requires to be defined for all r ∈ Ω,Ĉ =D = 0 can be concluded. In order to solve the homogeneous problem, the ansatz function is substituted in Eq. (A.11) using the linear coordinate transformation x(r ) = αr with and analogously Since Eq. (A.17c) is a homogeneous system of linear equations, the rows must depend on each other, which is fulfilled if the determinant of A vanishes. Otherwise, the solution to the linear equation system would be the trivial solution, i.e., C = D = 0. Hence, the solution to the homogeneous system of differential equations is found by adding the two linear independent solutions G r z (r ) = C 1 I 1 (αr ) + C 2 I 1 (βr ), where α and β can be obtained from solving Finally, after plugging Eqs. (A.18) into Eq. (A.8a), one obtains the following relations between C i and D i :

A.3.2 Inhomogeneous solution
Since the right-hand side of Eqs. (A.8) also contains an inhomogeneous part, the inhomogeneous problem must be solved as well. The solution of the inhomogeneous problem is called particular solution, in the following denoted by the subscript p, whereas the homogeneous solution is denoted by the subscript h. The final solution is obtained by adding the homogeneous and particular solutions: By inspecting this inhomogeneous part of Eqs. (A.8), it can be deduced that linear terms in G r z,p = Ar and G zr,p = Br solve the inhomogeneous PDEs, because the second derivatives vanish, while the first derivatives are eliminated. This leaves the following system of equations to be solved: As a consequence of the definitions of α 2 i (Eqs. (A.10)), the following equations are obtained to solve the inhomogeneous BVP:

A.3.3 Specific solution
To find the specific solution of the system of PDEs (A.8), one needs to substitute the boundary conditions into Eqs. (A.20) and solve them for C 1 , C 2 . Hence, the equality must hold. Note that there is a deviation toward [22], where a mix-up regarding both signs within the boundary conditions seems to have happened. Furthermore, in [22] the right-hand side is missing the factor 1/2, because it is contained in λ 2 , λ 3  Again, the equations can be transformed into a system of linear equations which is solved by inverting the system matrix: This inversion is only possible if f 1 = f 2 , thus, if K 1 = K 2 . The constants C 1 , C 2 can therefore be expressed as In order to obtain a dimensionless solution, the constants can also be rewritten in terms of the normalization velocity v 0 := − P R 2 /(4μ), as suggested in [22]: where the dimensionless wall velocity gradient can be expressed using Hence, the normalized microrotation and microshearing are expressed as where ρ = r/R is the normalized spatial coordinate, while G (i j) and G [i j] denote the i j-component of sym(G) and skw(G), respectively. These equations are now substituted in Eq. (A.4) to obtain the solution of the velocity field.

A.4 Solution for the velocity field
The derived solutions for the micro-motion G r z and G zr are substituted in Eq. (A.4). This yields the velocity gradient: v (r ) = 1 μ + κ/2 This expression can be integrated, which yields the general solution for the velocity field v(r ) = 1 μ + κ/2 where C 4 is a constant of integration. Using the homogeneous Dirichlet boundary condition on Γ 1 , the constant C 4 is found: 2α Finally, C 4 is substituted in Eq. (A.29), and the velocity field for a micromorphic pipe flow reads v(ρ) v 0 = 1 − ρ 2

B Investigation of experimental data from Bugliarello and Sevilla [4]
The experimental data from Bugliarello and Sevilla originates from a pressure-driven flow, and therefore, the characteristic diagram linking the pressure gradient to the volume flow rate were published in [4]. With the help of the volumetric velocity v m = Q/(π R 2 ), the volume flow rate can be computed, and the corresponding pressure gradient can be extracted from said plots. The knowledge of Q and P enables the computation of v 0 and with it a proper normalization of the velocity profiles. In order to verify that the experimental data matches the volumetric velocity, a micromorphic fluid was fitted to the data and v m was computed for different v 0 , see Figs. 17 and 18. For both R = 20 µm and R = 35 µm, the normalization term v 0 was chosen such that a) v m corresponds to the published value (dotted line), b) the error between velocity profile and data is minimal (solid line), and c) the velocity profile matches the data better and there exists a data point in the characteristic diagram (dashed line). In the end, the reference velocity v 0 was computed using the pressure gradient yielded by the dashed velocity profiles. This corresponds to the data points at (5.25 × 10 5 kgm −2 s −2 , 1.17 × 10 −2 mm 3 s −1 ) and (7 × 10 4 kgm −2 s −2 , 1.45 × 10 −2 mm 3 s −1 ) for R = 20 µm and R = 35 µm, respectively. The viscosity coefficient μ in v 0 corresponds to the asymptotic viscosity for large radius taken from [19], because the test conditions of the two studies were alike: Both [19] and [4] investigated the behavior of RBCs suspended in ACD at approximately 25 • C. It should be noted that due to the inconsistency within the experimental data, it cannot be assured that this normalization reflects the actual experimental results accurately. However, since the Fåhraeus-Lindqvist effect is satisfactorily reproduced by the micromorphic fluid model, it is also assumed that the micromorphic fluid exhibits a similar volume flow rate as an actual RBC suspension in pressure-driven pipe flows at large shear rates.