Abstract
For narrow tubes, red blood cells concentrate in the core region, leaving an annular zone called cellfree layer. This has an impact on both the tube hematocrit level (Fåhraeus effect) and the apparent blood viscosity (Fåhraeus–Lindqvist effect). Blood flow, mass transfer across the microvessel membrane, and diffusion in the tissue affect the solute concentration profiles. The Krogh tissue cylinder concept, limiting mass transfer to a cylinder around each microvessel, and the marginal zone concept (introduced by Haynes to analyze blood flow dynamics in narrow tubes) are both used to develop a model for solute transfer from blood in microvessels to the surrounding tissues, based on fundamentals. The analysis accounts for advection and diffusion in each zone of the microvessel, solute transport in the microvessel membranes, and diffusion and reaction in the tissues. The present investigation provides an analytical solution. The approach can be extended to treat other kinetic models, while accounting for Fåhraeus and Fåhraeus–Lindqvist effects in blood microvessels. The model is validated against published results for glucose transport from blood to tissue.
1 Introduction
Transport of a solute from blood to the surrounding tissues involves transport in the microvessel, through the microvessel wall, and transport in the tissue. Blood flow in the microvessels affects the concentration profiles. For small blood vessels below 500 μm diameter, diameter affects both the tube hematocrit (Fåhraeus effect) and the apparent blood viscosity (Fåhraeus–Lindqvist effect) with red blood cells (RBCs) concentrating in the core region [1,2,3,4,5,6,7,8]. Haynes marginal zone theory considers blood as consisting of two layers: a marginal zone (cellfree layer) surrounding a core region where RBCs concentrate [7]. Both zones are treated as continuum fluids. The theory was used to investigate Fåhraeus and Fåhraeus–Lindqvist effects [9,10,11]. Other blood flow models include shearinduced models [12,13,14,15] using the shearinduced particle migration concept [16,17] and elastic stressinduced models [18,19,20,21,22]. Non Newtonian and viscoelastic models for blood are reviewed in Arciero et al. [23] and in the references therein. Mesoscale models account for the particulate aspect of blood [23,24].
The Krogh cylinder model was introduced by Krogh to solve for transport of oxygen in tissue [25]. The assumptions used and extensions of the model are discussed in studies [23,26,27,28] and the references therein. In the case of oxygen transport, accounting for enhanced oxygen transport by myoglobin and reaction in RBCs is required [26,28]. Different geometries for RBCs were proposed [28]. An analytical solution, assuming cylindrical RBCs to occupy the lumen of the blood vessel shows the intracapillary high resistance to oxygen transport [29]. Groebe [30] accounted for axisymmetric diffusion of oxygen through three layers assuming cylindrical RBCs with plasma filling the gap between them. An analytical solution and a discussion of the simplifying assumptions are provided in ref. [30]. Blood is assumed to be formed by spherical RBCs, in line, equidistant, and surrounded by plasma in ref. [31]. A frame of reference moving with the RBCs is selected [31]. Other shapes for RBCs were selected and assumed to be surrounded by plasma: cylindrical (using a frame of reference in which RBCs are fixed) [32] and parachute form [33,34] (with a reference system in which RBCs are fixed in ref. [33], while considering RBCs motion in the frame of reference in the study [34]). Transport of other species like NO in relation with oxygen transport are discussed in refs. [23,27,28,30]. A model for biofluids flow in slowly and successively contracting and expanding rectangular geometry with two weakly permeable walls (allowing for flow through the walls) is solved using the homotopy method [35].
A model for solute transfer accounting for advection and diffusion using Krogh model for transport in tissue is presented in Fournier [3]. For intracapillary transport, blood velocity is assumed uniform inside the capillary and an overall transport coefficient from blood to the outer surface of the capillary wall (external side of the capillary membrane) is used in the final expressions for solute concentration in ref. [3].
The present model uses Krogh concept for solute transport in the tissue [25] (with mass transfer limited to a cylinder around each microvessel) along with Haynes marginal zone concept for blood flow [7] (considering the microvessel as formed of a cellfree layer and a core region where RBCs concentrate). The model for solute transport from blood to the surrounding tissue accounts for advection and diffusion in both the core region and the cellfree layer in the microvessel, solute transport in the microvessel membrane, and diffusion and reaction in the tissue (Figure 1). The approach is based on the fundamentals of fluid mechanics and mass transfer, where the velocity profiles in the two zones, core and cellfree layer, are considered to account for advection in the mass transfer equations. An analytical approach combining the different steps in solute transport is provided.
2 Governing equations
Haynes theory considers blood as consisting of two phases: an annular cellfree layer and a core layer containing RBCs. The velocity profile is discussed in Fournier [3], Chebbi [11], and the references therein. Only boundary conditions and related results of the present work are presented here. A schematic diagram of a microvessel and its Krogh cylinder is shown in Figure 1. The cylindrical coordinates r and z are used.
For fully developed flow, the velocity profile is independent of z. The following boundary conditions express continuity of the velocity and shear stress at the two zones interface r = r _{i}, in addition to finite shear stress at the axis of the microvessel and noslip boundary condition at the inner surface of the microvessel wall
Applying the momentum balance for the two zones yields
where v denotes velocity, μ represents viscosity, r is the radial coordinate, R is the microvessel radius, p _{gr} is the absolute value of the pressure gradient, and subscripts a, c, and i denote the annular region, the core region, and the twozone interfacial area, respectively.
The velocity profiles in the core region and the cellfree layer (annular region) are given by
where
The average velocity v ^{*} is given by
where Q is the blood volumetric flow rate.
The core hematocrit H _{C} is assumed uniform. The conservation of hematocrit yields the following relations [3,11]:
where H
_{T} is the tube hematocrit, H
_{D} is the discharge hematocrit, and
The apparent viscosity μ _{app} is given by the Hagen–Poiseuille equation and is related to the core zone viscosity by:
The core zone viscosity is assumed to be a function of H _{C} (up to 0.6) according to ref. [3]:
where
and T is blood temperature in Kelvin.
The mass transfer model accounts for advectiondiffusion in the microvessel side, mass transfer across the microvessel membrane, and diffusion and zeroorder reaction in the surrounding tissues demarcated by the Krogh cylinder. For small concentrations and nearly constant temperature, the diffusion coefficients are constant, and the use of conservation of mass of solute yields the following partial differential equations assuming steady state (zero accumulation term):
where z is the axial cylindrical coordinate, c denotes the solute concentration, D represents the diffusion coefficient of the solute, and subscripts c and a refer to the core and annular zones, respectively. The effective diffusion coefficient is used for the core region as it consists of plasma and RBCs. In the case of oxygen transport for instance, a reaction term involving the solute consumption needs to be accounted for in the solute mass balance. The present model does not include the case where reaction involving the solute occurs in the blood side.
The following boundary conditions express continuity of concentration and fluxes at the interface between the core and annular zones, in addition to finite concentration gradient at the centerline:
The effective diffusion coefficient D _{c} can be obtained using Maxwell equation [3,36]:
where D _{RBC} is the solute diffusion coefficient in RBCs. Mass transfer through the microvessel membrane is mainly due to diffusion, with negligible transport by filtration [3]. Mass transfer in the membrane is considered in the present model as purely diffusional, and the boundary condition at the microvessel wall can be written at the inner side of the membrane as [3]:
and as
at the outer side of the microvessel membrane, where subscript T refers to the tissue, and P _{M} and t _{m} denote the solute permeability and membrane thickness, respectively.
At the limit of the Krogh tissue cylinder, r = r _{T}, we have [3]:
which restricts mass transfer to the Krogh tissue cylinder surrounding the specific microvessel. This condition prevents mass transfer between the different Krogh cylinders.
In case the critical length L _{cr}, at which the solute concentration is zero, is larger than the microvessel length, then at the Krogh radius r _{T}, Eq. (22) is still valid, but the reaction term in the model needs to be restricted to a smaller part of the Krogh cylinder where the solute is non zero.
In the tissue, we assume zeroorder reaction occurring at a molar rate of R _{0} per unit volume. Assuming constant D _{T} for the same arguments listed before for D _{a} and D _{c}, and using the constant density approximation, transport of solute can be easily shown to be governed by (Appendix B.11 [36]):
Solving for c _{T} using boundary condition 22 gives [3]:
where
In order to determine the concentration profile, an expression for
3 Solution
Assuming sufficiently large Péclet number and large microvessel length, we neglect axial diffusion in the microvessel side and in its Krogh tissue cylinder. Excluding a small region near z = 0, we seek a solution of the form
where
The use of boundary conditions 16 and 20 yields
Substituting for velocities into Eqs. (14) and (15) while using Eqs. (6), (7), and (25) gives:
Boundary conditions (16)–(18) and (20) yield:
Integrating Eqs. (26) and (27) while using Eq. (29) gives:
where K is an integration constant.
In order to find b, a solute mass balance is made considering a control volume including the Krogh cylinder, the membrane, and blood limited to the interval z – z + Δz:
Considering the Krogh cylinder as our selected control volume gives:
Rearranging the above equation yields
Using Eq. (36) shows that the solute concentration profile, given by Eq. (24), is consistent with boundary condition (21).
Using the following definitions:
The dimensionless variables are defined as
The dimensionless profiles are obtained using Eqs. (6), (7), (31), and (32).
Converting Eqs. (8) and (9) to a dimensionless form gives:
Using the relationship between concentrations across the microvessel membrane, Eq. (30), yields in a dimensionless form:
Applying Eqs. (43) and (30) at the inner side of the microvessel wall, we get:
Differentiating Eq. (41) yields:
Applying Eq. (45) at
Assuming an average concentration c _{0} at z = 0 and applying a mass balance yields:
Using Eq. (34) and substituting for b in Eq. (47) gives:
Splitting the integral term in Eq. (48) into two terms to account for the existence of two zones yields in dimensionless form:
Substituting for the velocity and concentration profiles in Eqs. (39)–(41), we get:
Integrating gives:
where J _{c} and J _{a} are given by:
Applying Eq. (48) at the interface between the core and annular zones yields:
where Λ is given by:
Combining Eqs. (49), (52), (53), and (56) yields:
Using Eq. (44) yields the following expression for the centerline value
The mass transfer coefficient can be determined using:
The average solute concentration in blood c _{av} is given by:
in which v and c represent v _{a} and c _{a} in the annular zone, and v _{c} and c _{c} in the core zone, respectively. Then, using Eqs. (25) and (48) yields:
Substituting Eqs. (62) and (25) into Eq. (61) gives:
Then, combining Eqs. (63) and (36) provides:
4 Computational procedure and results
Solving the problem requires a number of steps. The dimensionless profiles
The dimensionless profiles for
and an intermediate case: D _{c}/D _{a} = 1. The following table provides an application of the abovementioned procedure as an example.
Table 1 shows the impact of diffusion in RBCs on the concentrations inside the microvessel, with obviously less deviations at high diffusion coefficient D
_{c} as anticipated, as higher diffusion coefficient yields less resistance to mass transfer. The concentrations decrease with increase in r at constant z. This allows solute transport from blood to tissue. Some negative values for
Constant  Limiting low D _{c} value  Intermediate case: D _{c} = D _{a}  Limiting high D _{c} value 

μ _{c}/μ _{a}  1.947  1.947  1.947 
D _{c}/D _{a}  0.4771  1.000  3.192 
H _{C}  42.2%  42.2%  42.2% 
μ _{app}/μ _{a}  1.469  1.469  1.469 
H _{D}  40.0%  40.0%  40.0% 

−2.937  −2.937  −2.937 

−1.509  −1.509  −1.509 

0.5581  0.5581  0.5581 

0.2343  0.2343  0.2343 

0.2993  0.1433  0.04554 

−0.3266  −0.1553  −0.04800 

−0.3783  −0.2070  −0.09970 
I _{c}  9.132 × 10^{−3}  4.591 × 10^{−3}  1.747 × 10^{−3} 
I _{a}  −9.132 × 10^{−3}  −4.591 × 10^{−3}  −1.747 × 10^{−3} 
As an example, the model is applied to glucose transport from blood flowing in a capillary to the surrounding Krogh tissue cylinder. For the purpose of comparison, the data and values from ref. [3] are used and summarized in Table 2.
Property  Symbol  Value  Ref. [3] 

D  Diffusion coefficient in water  0.91 × 10^{−5} cm^{2}/s  p. 221 
D _{a}  Diffusion coefficient in plasma  0.63 D  Eq. (5.60) 
2R  Capillary diameter  10 μm  Table 5.1 
v*  Average blood velocity  0.05 cm/s  Table 5.1 
t _{m}  Wall thickness  0.5 μm  Table 5.1 
L  Length  1,000 μm  Table 5.1 
R _{0}  Rate of glucose consumption  0.01 μmol/cm^{3} s  p. 221 
P _{M}  Membrane permeability  5.76 × 10^{−5} cm/s  p. 221 
c _{0}  Average concentration at z = 0  5 μmol/cm^{3}  p. 221 
a  Glucose molecular radius  0.36 nm  p. 221 
D _{T}/D  Ratio of D _{T} to D  ≈0.9  Figure 5.12 
r _{T}  Krogh cylinder radius (assumed in [3])  40 μm  p. 229 
The dimensional concentration profiles are obtained from Eq. (25) as:
At the critical value z _{cr} given by:
As noticed from Figure 2, all glucose concentrations drop as z increases due to solute transport from blood to tissue, followed by diffusion and chemical reaction in the Krogh tissue surrounding the capillary. The results for the average glucose concentration in blood, c _{av} (z), and the glucose concentration at the outer surface of the capillary wall, c _{T}(R + t _{m}, z), are both very consistent with Fournier’s results. The results for the concentrations at the Krogh tissue radius, c _{T}(r _{T}, z), are very close to those at the outer surface of the capillary wall as seen from Figure 2, and from Figure 5.18 in ref. [3]. The profiles for the deviations of concentration in blood from the average concentration, c(r,z)−c _{av}(z), show a decrease in concentration as r gets closer to the capillary inner wall as glucose is transported from blood to tissue, with less deviations from the average concentration of glucose, at higher D _{c}/D _{a} (Figure 3). The drop in glucose concentration in the Krogh cylinder from its value at the outer surface of the wall, c _{T}(r, z) − c _{T} (R + t _{m}, z), is plotted in Figure 4. Considering the previous results in Figures 2–4 and Table 3, we conclude that the significant change in the rdirection occurs in the membrane, as noted by Fournier [3].
Constant  Limiting low D _{c} value  Intermediate case: D _{c} = D _{a}  Limiting high D _{c} value 

b, μmol/cm^{3} μm  1.256 × 10^{−3}  1.256 × 10^{−3}  1.256 × 10^{−3} 

2.738 × 10^{−2}  2.738 × 10^{−2}  2.738 × 10^{−2} 
c _{T}(R + t _{m},z)−c _{av}(z), μmol/cm^{3}  −2.7356  −2.7309  −2.7280 
k _{m}, cm/s  1.516 × 10^{−2}  2.769 × 10^{−2}  5.750 × 10^{−2} 
z _{cr}, μm  1,791  1,795  1,797 
The mass transfer coefficient results are reported in Table 3, showing higher values at higher D _{c}/D _{a}. Nevertheless, the significantly lower permeability P _{M} plays the predominant role in determining the overall mass transfer coefficient as also mentioned by Fournier [3]. The critical values z _{cr} are close in the three D _{c}/D _{a} cases (about 1.8 mm), and larger than the capillary length in Table 2.
5 Conclusion
The concentration profiles are determined based on solute conservation of mass, mass transport equations, and relevant boundary conditions. The steady state model accounts for convectiondiffusion in the microvessel, transport through the microvessel membrane, and diffusion and zeroorder reaction in the Krogh tissue cylinder. Velocity gradients in the microvessel are accounted for in mass transfer calculations using the marginal zone theory results in the cellfree layer and the core region where RBCs concentrate. The solution for mass transfer requires a detailed procedure starting with blood flow calculations, and ending with the determination of the concentrations at the core zonecellfree layer interface, and both sides of the microvessel membrane. The use of the concentration values at the boundary points yields different profiles in the microvessel and in the Krogh cylinder using the relevant equations discussed above. The solution excludes a small region near z equal to zero. Large Péclet number and microvessel length are required for neglecting the axial diffusional terms in the microvessel and tissue sides. The estimation of diffusion coefficients and permeability is not discussed here and can be found in Fournier [3]. The approach used here can be extended to treat different kinetic models, while accounting for Fåhraeus and Fåhraeus–Lindqvist effects in blood microvessels. The use of Maxwell equation to determine the effective diffusivity of the solute in the core region of the microvessel is not restricted to low volume fractions of the dispersed phase (RBCs in our case) as noted in ref. [36]. The model is validated against the results of Fournier [3] for the average glucose concentration in blood, and the glucose concentrations at the outer surface of the capillary wall and at the Krogh cylinder radius.

Funding information: The authors state no funding involved.

Author contributions: All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

Conflict of interest: The authors state no conflict of interest.
References
[1] Fåhraeus R. The suspension stability of blood. Physiol Rev 1929;9:241–74.10.1152/physrev.1929.9.2.241Search in Google Scholar
[2] Fåhraeus R, Lindqvist T. The viscosity of the blood in narrow capillary tubes. Am J Physiol. 1931;96:562–8.10.1152/ajplegacy.1931.96.3.562Search in Google Scholar
[3] Fournier RL. Basic transport phenomena in biomedical engineering. Boca Raton: CRC Press; 2012.10.1201/b14885Search in Google Scholar
[4] Goldsmith HL, Cokelet GR, Gaehtgens P. Robin Fåhræus: evolution of his concepts in cardiovascular physiology. Am J Physiol Heart Circ Physiol. 1989;257:H1005–15.10.1152/ajpheart.1989.257.3.H1005Search in Google Scholar PubMed
[5] Toksvang LN, Berg RMG. Using a classic paper by Robin Fåhraeus and Torsten Lindqvist to teach basic hemorheology. Adv Physiol Educ. 2013;37:129–33.10.1152/advan.00009.2013Search in Google Scholar PubMed
[6] Secomb TW, Pries AR. Blood viscosity in microvessels: experiment and theory. Comptes Rendus Phys. 2013;14:470–8.10.1016/j.crhy.2013.04.002Search in Google Scholar PubMed PubMed Central
[7] Haynes RF. Physical basis of the dependence of blood viscosity on tube radius. Am J Physiol. 1960;198:1193–200.10.1152/ajplegacy.1960.198.6.1193Search in Google Scholar PubMed
[8] Pries AR, Neuhaus D, Gaehtgens P. Blood viscosity in tube flow: dependence on diameter and hematocrit. Am J Physiol Heart Circ Phys. 1992;263:H1770–8.10.1152/ajpheart.1992.263.6.H1770Search in Google Scholar PubMed
[9] Sharan M, Popel AS. A twophase model for flow of blood in narrow tubes with increased effective viscosity near the wall. Biorheology. 2001;38:415–28.Search in Google Scholar
[10] Sriram K, Intaglietta M, Tartakovsky DM. NonNewtonian flow of blood in arterioles: consequences for wall shear stress measurements. Microcirculation. 2014;21(7):628–39.10.1111/micc.12141Search in Google Scholar PubMed PubMed Central
[11] Chebbi R. Dynamics of blood flow: modeling of the Fåhræus–Lindqvist effect. J Biol Phys. 2015;41(3):313–26.10.1007/s1086701593761Search in Google Scholar PubMed PubMed Central
[12] Weert KV. Numerical and experimental analysis of shearinduced migration in suspension flow. (A thesis for the degree of Master). Eindhoven, Netherlands: Eindhoven University; 2005.Search in Google Scholar
[13] Mansour MH, Bressloff NW, Shearman CP. Red blood cell migration in microvessels. Biorheology. 2010;47:73–93.10.3233/BIR20100560Search in Google Scholar PubMed
[14] Chebbi R., Dynamics of blood flow: modeling of Fåhraeus and FåhraeusLindqvist effects using a shearinduced red blood cell migration model. J Biol Phys. 2018;44:591–603.10.1007/s1086701895085Search in Google Scholar PubMed PubMed Central
[15] Chebbi R. A twozone shearinduced red blood cell migration model for blood flow in microvessels. Front Phys. 2019;7:206.10.3389/fphy.2019.00206Search in Google Scholar
[16] Leighton DT, Acrivos A. The shearinduced migration of particles in concentrated suspension. J Fluid Mech. 1987;181:415–39.10.1017/S0022112087002155Search in Google Scholar
[17] Phillips RJ, Armstrong RC, Brown RA. A constitutive equation for concentrated suspensions that accounts for shearinduced particle migration. Phys Fluids. 1992;4:30–40.10.1063/1.858498Search in Google Scholar
[18] MoyersGonzalez M, Owens RG, Fang J. A nonhomogeneous constitutive model for human blood. Part 1. Model derivation and steady flow. J Fluid Mech. 2008;617:327–54.10.1017/S002211200800428XSearch in Google Scholar
[19] MoyersGonzalez MA, Owens RG. Mathematical modelling of the celldepleted peripheral layer in the steady flow of blood in a tube. Biorheology. 2010;47(1):39–71.10.3233/BIR20100558Search in Google Scholar PubMed
[20] Dimakopoulos Y, Kelesidis G, Tsouka S, Georgiou GC, Tsamopoulos J. Hemodynamics in stenotic vessels of small diameter under steady state conditions: Effect of viscoelasticity and migration of red blood cells. Biorheology. 2015;52(3):183–210.10.3233/BIR14033Search in Google Scholar PubMed
[21] Mavrantzas VG, Beris AN. Modelling the rheology and the flowinduced concentration changes in polymer solutions. Phys Rev Lett. 1992;69:273–6. Errata. 1993;70,2659.10.1103/PhysRevLett.69.273Search in Google Scholar PubMed
[22] Tsouka S, Dimakopoulos Y, Mavrantzas V, Tsamopoulos J. Stressgradient induced migration of polymers in corrugated channels. J Rheol. 2014;58(4):911–47.10.1122/1.4880245Search in Google Scholar
[23] Arciero JC, Causin P, Malgaroli F. Mathematical methods for modeling the microcirculation. AIMS Biophys. 2017;4:362–99.10.3934/biophy.2017.3.362Search in Google Scholar
[24] Bessonov N, Sequeira A, Simakov S, Vassilevski Yu, Volpert V. Methods of blood flow modelling. Math Model Nat Phenom. 2016;11:1–25.10.1051/mmnp/201611101Search in Google Scholar
[25] Krogh A. The number and distribution of capillaries in muscles with calculations of the oxygen pressure head necessary for supplying the tissue. J Physiol. 1919;52:409–15.10.1113/jphysiol.1919.sp001839Search in Google Scholar
[26] Popel AS. Theory of oxygen transport to tissue. Crit Rev Biomed Eng. 1989;17:257–321.Search in Google Scholar
[27] Goldman D. Theoretical models of microvascular oxygen transport to tissue. Microcirculation 2008;15:795–811.10.1080/10739680801938289Search in Google Scholar
[28] Truskey GA, Yuan F, Katz DF. Transport phenomena in biological systems. New Jersey: Pearson; 2010.Search in Google Scholar
[29] Hellums JD. The resistance to oxygen transport in the capillaries relative to that in the surrounding tissue. Microvasc Res. 1977;13:131–6.10.1016/00262862(77)901224Search in Google Scholar
[30] Groebe K. An easytouse model for O2 supply to red muscle. Validity of assumptions, sensitivity to errors in data. Biophys J. 1995;68:1246–69.10.1016/S00063495(95)803000Search in Google Scholar
[31] Federspiel WJ, Popel AS. A theoretical analysis of the effect of the particulate nature of blood on oxygen release in capillaries. Microvasc Res. 1986;32:164–89.10.1016/00262862(86)90052XSearch in Google Scholar
[32] Eggleton CD, Vadapalli A, Roy TK, Popel AS. Calculations of intracapillary oxygen tension distributions in muscle. Math Biosci. 2000;167:123–43.10.1016/S00255564(00)000389Search in Google Scholar
[33] Vadapalli A, Goldman D, Popel AS. Calculations of oxygen transport by red blood cells and hemoglobin solutions in capillaries. Artif Cell Blood Substit Immobil Biotechnol. 2002;30:157–88.10.1081/BIO120004338Search in Google Scholar
[34] Lucker A, Weber B, Jenny PA. Dynamic model of oxygen transport from capillaries to tissue with moving red blood cells. Am J Physiol Heart Circ Physiol. 2015;308:H206–H216.10.1152/ajpheart.00447.2014Search in Google Scholar PubMed
[35] Dinarvand S, Rashidi MM, Doosthoseini A. Analytical approximate solutions for twodimensional viscous flow through expanding or contracting gaps with permeable walls. Cent Eur J Phys. 2009;7(4):791–9.10.2478/s115340090024xSearch in Google Scholar
[36] Bird RB, Stewart WE, Lightfoot EN. Transport phenomena. New York: Wiley; 2007.Search in Google Scholar
© 2022 Rachid Chebbi, published by De Gruyter
This work is licensed under the Creative Commons Attribution 4.0 International License.