Next Article in Journal
Convergence of Special Sequences of Semi-Exponential Operators
Previous Article in Journal
Enhancing the Transferability of Adversarial Examples with Feature Transformation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Extended Thermodynamics Model for Blood Flow

Department of Mathematical, Computer, Physical and Earth Sciences, University of Messina (Italy), V.le F. Stagno D’Alcontres 31, I-98166 Messina, Italy
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(16), 2977; https://doi.org/10.3390/math10162977
Submission received: 20 July 2022 / Revised: 12 August 2022 / Accepted: 17 August 2022 / Published: 18 August 2022
(This article belongs to the Section Mathematical Physics)

Abstract

:
A model for blood flow is introduced in the context of the Rational Extended Thermodynamics (RET). The balance equations are applied to the two-hierarchy structure recently introduced by Ruggeri and Sugiyama. The constitutive relations are derived with universal physical principles and the remaining constitutive functions are evaluated by use of the kinetic theory. The model herein obtained is a hyperbolic generalization of a classical blood flow model. Our equations by construction have the same physical proprieties of the classical system; in addition, owing to its hyperbolic structure, our model avoids the unphysical feature of instantaneous diffusive effects which is typical of parabolic systems. Furthermore we expect that our model, as all RET systems, can describe the physical phenomena better than the classical ones when the fields change rapidly or one has steep gradients.

1. Introduction

Gidaspow in 1994 [1] and Gidaspow and Huang in 2009 [2] employed and studied a model for blood flow based on the kinetic theory which assumes blood as a mixture of plasma and red blood cells (RBCs). This model allows to explain the Fahraeus–Lindqvist effect, the decrease of RBCs viscosity when the diameter of blood vessels decreases. The cause of this effect is the random RBCs kinetic energy per unit mass, known as granular temperature. When this temperature rises near the wall of the vessel, RBCs oscillate and they migrate from the wall to the center of the vessel. This mechanism is called the shear induced diffusion and it is a driving force connected to a gradient of shear. It is observed with the microscope that the smaller the diameter of the blood vessel, the stronger the shear diffusion. This creates a layer of free cells near the wall and consequently a collapse of the viscosity in the narrow vessels. The role of kinetic theory is to provide a theoretical expression for the blood viscosity. The results of Gidaspow’s paper [2] are in agreement with the experimental data collected for red blood cell concentrations in narrow vessels. Inspired by this work [2], many researchers have developed models on blood flow to study physiological and patalogical cases. Huang, in her PhD dissertation, applies Gidaspow’s theory to the transport of LDL and HDL lipoproteins and platelets. In [3], the pulsing flow in a coronary artery is studied using a typical cardiac waveform. The viscosity values of the red blood cells are numerically determined to appreciate the effect and to identify sites where arteriosclerosis could develop. Fedosov et al. [4] deal with the changes in the shape, mechanics and physical properties of red blood cells, determined by infectious diseases, environmental and hereditary factors. They evaluate how the biorheology of blood flow varies in subjects suffering from malaria and sickle cell anemia and they compare with healthy cases. Felicetti et al. [5] numerically treat the detection and the treatment of tumors by recognizing tumor markers on the blood cell surfaces. The information is conveyed on one side by nanomachines that have similar size to white blood cells and they can flow in vessels. On the other side, smart probes are used to detect information extrabody about the localization of tumors. Wu et al. [6] study benchmark problems including a sudden expansion and various driven slots and crevices through the vessels. Their results turned out to be a predictive model for thrombosis in blood-wetted devices. Lopes et al. [7] treat the flow of blood in the carotid vessels and compare models in which the vessels are idealized with rigid walls and elastic walls. They focus on the properties of the arteries that become more rigid with the appearance of arteriosclerosis. Furthemore, they speculate that the carotid sinus is the area where arteriosclerosis is most likely to occur, given its low values of Time-Averaged Wall Shear Stress. A 3-phase model, consisting of blood, red blood cells and white blood cells, is formulated by Melka et al. [8]. In particular, they analyze the aorta and its thoracic branches of a child of 8 years, affected by a congenital heart defect, the coarctation of the aorta. Another 3-phase model, this time consisting of plasma, platelets and red blood cells, is explained by Gidaspow and Chandra [9]. It is observed numerically that platelets and red blood cells move in opposite directions. The high shear on the walls of the vessel causes an increase in granular pressure and temperature that makes the red blood cells oscillate and push them towards the center. On the contrary, the platelets push towards the walls. From these articles, it is clear how expendable the Gispaspow’s theory is in scientific research.
In this paper, we start from this classical model by Gidaspow and Huang [2] in order to derive a hyperbolic system in the context of Rational Extended Thermodynamics [10,11,12]. Following this non-equilibrium theory, we consider as field variables not only the classical ones in [2], which are mass densities, velocities and RBCs temperature, but also the stress tensors, the RBG dynamic pressure and the heat flux. In order to describe this extended set of field variables, we need to introduce a new system of field equations. Following the idea of RET, we introduce for these fields a new set of balance laws, and, in particular, here we are inspired by a two-hierarchy structure model recently introduced in [11,12] for polyatomic gases. The balance equations needs closure and, following again the main idea of RET, we close them assuming local and instantaneous constitutive relations and using well-known universal physical principles, such as the entropy principle (existence of an entropy inequality and concavity of the entropy density) and the principle of relativity. The exploitation of the principle of relativity helps us in the determination of the velocity-dependence of the constitutive relations, while the principle of entropy implies some useful mathematical relations. The mathematical exploitation of the entropy principle is made by use of the Lagrange multipliers [13] (main field). The constitutive functions are all determined except for the thermal and caloric equation of state and the production terms.
We identify the remaining terms using the expressions of the RBCs pressure, viscosity and heat conductivity determined by use of the kinetic theory in [1]. In this way, we assure that our equations coincide with the classical GH model when we are in presence of slowly changing fields and small gradients. In addition, our model can describe other processes where rapid time changes or when a strong deviation from equilibrium occur. In fact, it has been shown that the field equations of RET can describe a range various non-equilibrium phenomena such as light shattering, sound waves, heat waves, structure of shock waves [10,11,12]. RET has been applied to monoatomic gases [10] and mixtures [10,11,14,15] with many interesting results. Recently, RET has been generalized to dense and rarefied polyatomic gases both in the classical [11,12,16,17,18] and in the relativistic framework [12,19,20], for metal electrons [21,22], to quantum systems [23] and also for biological systems [24,25,26,27], providing relevant results and good agreement with experimental data.
The model herein introduced, as all ET systems [28,29], is symmetric hyperbolic. Hyperbolicity guarantees finite speeds of propagation, while symmetric hyperbolic systems imply well-posedness of Cauchy problems (i.e., existence, uniqueness and continuous dependence on the data). Hyperbolic models overcome the paradox of infinite propagation speed of disturbances, and they are better suited to describe transient regimes. Moreover, they offer a richer scenario of dynamics which can be useful to describe biological problems, see for example [24,25,26,27] and the references therein. We expect that the model herein introduced can also describe phenomena in blood flow in which a rapid change in the blood pressure occurs, such as the sudden occlusion of a vessel or rapid initial changes in the fields due to external conditions.
Furthermore, from the results obtained for the heat transfer problem in a monoatomic gas or gas-mixtures, see for example [14,15,30,31] and the references therein, we also expect that our model can show more interesting spatial effects also for time-independent phenomena. Indeed, it has been shown that the field equations of RET are able to put in evidence more spatial effects than classical Navier–Stokes–Fourier equations, when a radial symmetry is present or a particular geometry is taken into account. Therefore, we expect that our model can describe more accurately the blood flow in a vessel and in particular in the case of change in the radial symmetry, like stenosis, or other defects of vessels.
Finally, it was recently shown that the field equations of RET are to describe the thermal-diffusion effect, contrary to the classical equations. In consideration of this RET result, we expect that our model is able to show some interesting additional features to the description of the Fahraeus–Lindqvist effect, which is anyway already present in the Gidaspow classical model [2].
The remainder of this article is the following: In Section 2, we recall the Gidaspow classical model, while in Section 3, we introduce the balance equations for our model. Section 4, Section 5 and Section 6 are devoted to the applications of the physical principles for the determinations of the constitutive relations. In particular, the Galilean invariance principle helps us to determine the dependence of the constitutive relations on the velocity, while the entropy principle implies some algebraic relations for the constitutive relations. In order to avoid too long calculations, we limit the study near the thermodynamic equilibrium by expanding the constitutive relations in terms of the non-equilibrium fluxes. In Section 7, a summary of the consequences of the principle is made, and in Section 8, the comparison with the Gidaspow model is shown with the identification of the last unknown functions. So all coefficients and functions here introduced can be obtained in terms of those calculated by Gidaspow using the kinetic theory. In Section 9, some simple calculations are presented in order to show with a simple numerical example the applicability of our model. Finally, in Section 10, some conclusions and final remarks are discussed.

2. Gidaspow Model

Gidaspow’s model [2] describes blood as a mixture of a fluid and the RBCs and uses as field equations the two conservation laws of mass for the two constituents, the two balance laws of momentum and a balance law for the total energy of RBCs. These equations, typical of classical thermodynamics, can be recast in the following form:
ρ f t + ρ f v k f x k = 0 , ρ s t + ρ s v k s x k = 0 , ρ f v i f t + ρ i k f + ρ f v i f v k f x k = G i f + I i f , ρ s v i s t + ρ i k s + ρ s v i s v k s x k = G i s + I i s , t ρ s e s + 1 2 ρ s v l s v l s + x k q k s + ρ k l s v l s + ρ s e s v k s + 1 2 ρ s v l s v l s v k s = v l s G l s + I s + v l s I l s ,
where ρ f and ρ s are respectively the plasma and the RBCs densities, v k f and v k s their velocities, ρ i k f and ρ i k s the stress tensors, e s is the RBCs specific interval energy and q k s is the RBCs heat flux. The field Equation (1) are closed by the classical Navier–Stokes and Fourier equations, which assume the stress tensor and the heat flux proportional to the gradient of the velocity or the temperature:
ρ i j s = 2 μ s v ( i s x j ) + 2 3 μ s ξ s v k s x k δ i j ρ i j f = 2 μ f v ( i f x j ) 2 3 v k f x k δ i j , q i s = K s θ s x i .
In the previous equations and through the paper, the round brackets in the indexes denote the symmetric parts of a tensor. The viscosities μ s , ξ s and μ f and the heat conductivity K s are calculated by Gidaspow by use of the kinetic theory and written explicitly in [2]. The right-hand side terms in Equation (1) are the body forces
G i f = ρ f g i + ε s i p f , G i s = ρ s g i ε s i p f ,
due to external specific force g i , the actions that the fluid exerts on the RBCs ε s i p f and vice versa ( ε s is the volume fraction of RBCs). Furthermore, the production terms take into account the interaction between plasma and RBCs, that are
I i f = I i s = β v i s v i f , I s = γ .
Moreover, the coefficients β and γ are evaluated in [2].
Equations (1)–(3) form the classical system of field equations, which consists of nine equations in the nine fields ρ α , v i α and θ α .

3. Balance Equations

Starting from Gidaspow’s model, we introduce a system of hyperbolic balance equations. First of all, we consider as field variables not only the densities, the velocities and the temperature but also the stress tensors and the heat flux. The fields equations for all these variables are the two conservation laws of mass for the two constituents, the two balance laws of momentum and the balance law for the total energy of RBCs, as in (1). Furthermore, instead of assuming the classical Navier–Stokes and Fourier expressions (2) for the stress tensors and for the heat flux, we introduce for these new variables supplementary balance equations. In particular, the system of balance equations that we assume has the following form:
ρ α t + ρ α v k α x k = 0 , ρ α v i α t + F i k α x k = G i α + I i α = H i α , M l l α t + M k l l α x k = 2 v l α H l α + 2 I α , F i j α t + F i j k α x k = 2 v ( i α H j ) α + S i j α , M i l l s t + M i k l l s x k = 3 v ( i s v l s H l ) s + 2 ρ i l s ρ s H l s + 2 e s H i s + 2 v i s I s + S i l l s ,
with α = f or s.
Equation (5) 1 , 2 and (5) 3 with α = s coincide exactly with the classical Gidaspow Equation (1) indicating with 2 M l l α total energy and 2 M i l l α the total energy flux. Equation (5) 3 with α = f is a supplementary equation here introduced for the total energy for the fluid. The other equations of (5) are introduced assuming that the fluxes in the first equations are the densities in the new equations. This structure was obtained by integration from the Boltzmann equation for a monatomic gasses [10] and recently this structure was modified for polyatomic gasses [11,12,17], where the two hierarchies structure (for the F’s and the M’s) is introduced by considerations from the kinetic theory.
System (5) represents a set of 24 equations for the 24 field variables that are the two densities ρ α , the two velocities v i α , the two temperatures θ α , the traceless (here and in the following, square brackets in the indexes denote the traceless part of a symmetric tensor) part stress tensors ρ < i j > α , the RBCs dynamic pressure Δ s and heat flux q i s . For sake of simplicity, we assume that q i f = 0 . The aim of this paper is to close system (5), that is to express all quantities in (5) in terms of these 24 fields.
First of all, in the spirit of Rational Extended Thermodynamics, we assume that the unknown quantities in (5) are local and instantaneous. This means that is at one time and one point they depend on the field variables at the same time and at the same point. So, indicating a generic unknown quantity with φ , it must hold
φ = φ ρ α , v i α , θ α , ρ < i j > α , Δ s , q i s .
Then, if the functions φ are known, we can insert them into the balance Equation (5) in order to obtain a closed system of 24 field equations in the 24 fields. The obtained system, for the assumption in (6), will be a set of quasi-linear first order differential equations, whose solution is called thermodynamic process.

4. Decomposition in Convective and Non Convective Terms

We start now to decompose the densities, fluxes and productions in (5) into convective and non-convective parts, according to the single component case. Indeed, assuming that the balance equations for each constituent of the mixture must be invariant under a Galilean transformation, it is possible to determine the velocity dependence of the quantities in (5) on the two velocities, so we have:
F i k α = ρ i k α + ρ α v i α v k α , F i j k α = ρ i j k α + 3 ρ ( i j α v k ) α + ρ α v i α v j α v k α , M l l α = 2 ρ α e α + ρ α v l α v l α , M i l l s = 2 q i s + 2 ρ i l s v l s + 2 ρ s e s v i s + ρ s v l s v l s v i s , M i k l l s = m i k l l s + 2 ρ i k l s v l s + 4 q ( i s v k ) s + 4 ρ l ( i s v k ) s v l s + ρ i k s v l s v l s + 2 ρ s e s v i s v k s + ρ s v i s v k s v l s v l s , S i j α = s i j α , S i l l s = s i l l s + 2 v l s s i j s .
The explicit calculations for expressions (7) are shown in [10] pp. 53, 54.
It is easy to recognize the decomposition (7) 1 of the moment fluxes in the stress tensors and velocity dependent parts, the decomposition (7) 3 of the total energies in internal and kinetic energies and the decomposition (7) 4 for the total energy flux M i l l s in the sum of the heat flux and the velocity dependent parts. These three decompositions are already present in (1). Then, we may say that ρ i k α in (7) are the stress tensors which, according to Gidaspow’s model, are expressed by
ρ i k s = ρ < i k > s + p s + Δ s δ i k , ρ i k f = ρ < i k > f + p f δ i k ,
while q i s is the RBCs heat flux.
Then, by insertion of the decomposition (7) into the balance Equation (5), one obtains
d α ρ α d t + ρ α v k α x k = 0 , ρ α d α v i α d t + ρ i k α x k = H i α , d α d t ρ α e α + ρ α e α v k α x k + ρ k l α v l α x k + q k α x k = I α , d α ρ i j α d t + ρ i j α v k α x k + 2 ρ k ( i α v j ) α x k + ρ k i j α x k = s i j α , 2 d s q i s d t + m i k l l s x k + 2 ρ i k l s v l s x k + 4 q ( i s v k ) s x k 2 ρ i l s ρ s ρ l k s x k 2 e s ρ i k s x k = s i l l s ,
where the two material derivatives d α / d t = / t + v k α / x k have been introduced. It can be clearly seen that these equations do not depend on the velocities but only on their derivatives and this is in agreement with the Galileian invariance applied to each constituent of the mixture.
With these decompositions, the unknown functions become ρ i j k α , m i k l l s , s i j α and s i l l s . These quantities must be determined in terms of the set ρ β , θ β , ρ < i j > β , Δ s , q i s and, in virtue of the principle of material frame indifference applied to the whole mixture ([10] p. 86, [11] p. 313, [12,32] pp. 6, 177, [33]), they cannot depend on the velocity of the two components v i β but only on their relative velocity V i s = v i s v i f . Equivalently, they can be considered as functions of the diffusion velocity
u i s = v i s v i where v i = ρ s v i s + ρ f v i f ρ s + ρ f
is the mixture velocity. Summarizing, the constitutive functions to be determined are
ρ i j k α = ρ i j k α ρ β , θ β , u i s , ρ < i j > β , Δ s , q i s , m i k l l s = m i k l l s ρ β , θ β , u i s , ρ < i j > β , Δ s , q i s , s i j α = s i j α ρ β , θ β , u i s , ρ < i j > β , Δ s , q i s , s i l l s = s i l l s ρ β , θ β , u i s , ρ < i j > β , Δ s , q i s .

5. Entropy Principle

The entropy principle asserts that it exists a quantity h, called entropy, which satisfies a balance equation of the form
h t + h k x k = Φ 0
for all thermodynamic processes, that is for all solutions of the field equations. In this paper, in line with the methods of Rational Extended Thermodynamics [10,11,12], we use this physical principle in order to determine the constitutive Equation (11).
First of all, as described before for the fields, the concave entropy density h, the entropy flux, h k = h v k + ϕ k , and the entropy production Φ must depend on the field variables through constitutive functions of the form
h = h ρ β , θ s , u i s , ρ < i j > β , Δ s , q i s , ϕ k = ϕ k ρ β , θ β , u i β , ρ < i j > β , Δ s , q i s , Φ = Φ ρ β , θ β , u i β , ρ < i j > β , Δ s , q i s .
Then, since the entropy principle must hold for all thermodynamic processes that is for all solutions of the field Equation (9), we use the field Equation (9) as constraints for the validity of the entropy principle, This can be mathematically done, by introduction [13] of the so-called Liu–Lagrange multipliers λ ’s, so the entropy inequality (12) becomes
h t + h v k + ϕ k x k α = f , s λ α d α ρ α d t + ρ α v k α x k + 1 ρ α ρ α λ i α 2 ρ i k s λ k n n s 3 ρ s θ s λ i n n s ρ α d α v i α d t + ρ i k α x k H i α + λ ^ l l α d α d t ρ α e α + ρ α e α v k α x k + ρ k l α v l α x k + q k α x k I α + λ i j α d α ρ i j α d t + ρ i j α v k α x k + 2 ρ k ( i α v j ) α x k + ρ k i j α x k s i j α λ i n n s 2 d s q i s d t + m i k l l s x k + 2 ρ i k l s v l s x k + 4 q ( i s v k ) s x k 2 ρ i l s ρ s ρ l k s x k 2 e s ρ i k s x k s i l l s   =   Ψ   0 ,
which must be valid for all values of ρ β , θ s , u i s , ρ < i j > β , Δ s and q i s .
We start now the mathematical evaluations of (14). We have to say that the introduction of the Lagrange multipliers adds new unknown quantities. In agreements with the other quantities, we assume that also the λ ’s must depend on the fields through constitutive relations of the form (11) and (13).
It is easy to see that in (14) there is only one term which depends on the external force H i α so, since the constitutive relations must be valid independently of the forces applied to the system, we can conclude that
ρ α λ i α 2 ρ i k s λ k n n s 3 ρ s θ s λ i n n s = 0
must hold.
Then, in equality (14) there are only few terms that depend only on the fields and not on their derivatives. These terms are the production terms in the balance equations and the entropy production, so we recover the so-called residual inequality, that yields the entropy production in terms of the other quantities
Φ = λ ^ l l s I s + α λ i j α s i j α + λ i n n s s i l l s 0
and whose sign must be verified when the other quantities are determined.
Another relation, that can be easily determined form (14), is the following:
α λ ^ v v s ρ i j s + 2 λ i v α ρ j v α + 2 λ v n n s ρ i j v s + q j s δ i v = h α ρ α λ α + 3 2 ρ s θ s λ l l s + λ k v α ρ k v α + 2 q v s λ v l l s δ i j .
This relation comes out eliminating all coefficients of the velocity v k from equality (14).
We can now insert the first relations (15)–(17) into the balance law (14) in order to obtain a simpler equation. A very accurate evaluation of the complete entropy principle requires a lot of complex calculations that we prefer not to present here. In the next section, we start by limiting our analysis to processes near the thermodynamic equilibrium. A more general case will be presented in [34].

6. Evaluation of the Quantities Near Equilibrium

As explained before, from now on we limit our attention to processes close to thermodynamic equilibrium, characterized by vanishing non-equilibrium fluxes ρ < i j > α , Δ s , u i α and q i α . This can be done by mathematically expanding the constitutive relations around the thermodynamic equilibrium and considering the terms which are at least quadratic in the non-equilibrium fluxes. In this way the entropy quantities become
h = α h 0 α + h 1 α ρ < i j > α ρ < i j > α + h 2 α q i α q i α + h 3 α Δ α Δ α + h 4 α u i α u i α + h 5 α q i α u i α + O 3 , ϕ i = α ϕ 1 α q i α + ϕ 2 α u i α + ϕ 3 α ρ < i j > α q j α + ϕ 4 α ρ < i j > α u j α + ϕ 5 α q i α Δ α + ϕ 6 α u i α Δ α + O 3 ,
the constitutive relations are given by
ρ i j k α = 3 β 1 α q ( i α δ j k ) + β 2 α u ( i α δ j k ) + O 2 , m i j l l s = γ 1 s ρ < i j > s + 1 3 γ 0 s + γ 2 s Δ s δ i j + O 2
and the Lagrange multipliers take the form
λ α = λ 0 α + λ 2 α Δ α + O 2 , λ ^ l l α = ν 0 α + ν 2 α Δ α + O 2 , λ i j α = r 1 α ρ < i j > α + 1 3 r 0 α + r 2 α Δ α δ i j + O 2 , λ i l l s = s 1 q i s + s 2 u i s + O 2 .
So the determination of the constitutive functions reduces to the calculations of all coefficients in (18)–(20), which depends only the equilibrium variables ρ α , θ α .
We then proceed by insertion of relations (15)–(20) into the entropy inequality (14) which becomes a very long but linear equation in the derivatives of the field variables ρ β , θ s , u i s , ρ < i j > β , Δ s and q i s . As already said, Equation (14) must be valid for all field variables; therefore, it must be valid also for all gradients of these variables. Then, setting the coefficients of the derivatives of ρ β , θ s , u i s , ρ < i j > β , Δ s and q i s equal to zero, we obtain a big set of algebraic equations for the determination of the functions in (18)–(20). We prefer not to present all these equations here, but an interested reader can find similar calculations in [10,11,12,34]. We present here only some equilibrium relations and the results.

6.1. Gibbs Equation

We evaluate the coefficients in expansions (18)–(20), starting with the equilibrium terms. Indeed, from the derivatives of ρ α and θ α in (18)–(20) we obtain
h 0 ρ α = h 0 α ρ α = λ 0 α + ν 0 α ρ α ρ α e α , h 0 θ α = h 0 α θ α = ν 0 α θ α ρ α e α .
which can be recasted in the differential form
d h 0 α ρ α = ν 0 α e α θ α d θ α + λ 0 α ρ α + ν 0 α e α ρ α + ν 0 α e α ρ α h α ρ α 2 d ρ α .
This last equation can be compared with the well-known Gibbs relation, that expresses the variation of the equilibrium entropy density in terms of the thermodynamic variables. The expression of the Gibbs relation appropriate to each component of a mixture can be found for example in [10] p. 90 and it assumes the form
d h 0 α ρ α = 1 θ α d e α p α ρ α 2 d ρ α .
Then, by comparison of (22) and (23), taking also into account that e α = e α ρ α , θ α , it is possible to determine the equilibrium part of the Lagrange multipliers, i.e.,
ν 0 α = 1 θ α , λ 0 α = h α ρ α p α ρ α θ α e α θ α = μ α θ α , r 0 α = 0 ,
where μ α represents the chemical potential.

6.2. Evaluation of the Functions Outside Equilibrium

By use of the Lagrange multipliers at equilibrium (24) it is possible to evaluate from the entropy inequality all coefficients in the constitutive relations (14) and (18)–(20) as explained before. We have the entropy quantities
h = α = f , s ρ α μ α θ α + 1 θ α e α + p α ρ α ρ < i j > α ρ < i j > α 4 p α θ α + q i s q i s Ψ s + Δ α Δ α 2 θ α Γ α , ϕ i = q i s θ s + α = f , s ρ α μ α θ α + 1 θ α e α + p α ρ α u i α γ 1 s 2 e s + p s ρ s Ψ s ρ < i j > s q j s + 1 θ s Γ s p θ s ρ s c V s + 5 3 p α θ s γ 1 s 2 e s + p s ρ s Ψ s q i s Δ s ,
the first order terms in the constitutive relations
ρ i j k α = 3 β 1 α q ( i α δ j k ) , m i j l l s = γ 1 s ρ < i j > s + 1 3 γ 0 s + γ 2 s Δ s δ i j
and in the Lagrange multipliers
λ α = μ α θ α + 1 θ α Γ α p ρ α + p θ α ρ α c V α e α + ρ α e ρ α Δ α , λ i α = 2 Ψ α p α ρ α + e α q i α , λ ^ l l α = 1 θ α 1 θ α Γ α p θ α ρ α c V α Δ α , λ i j α = 1 2 θ α p α ρ < i j > α + 1 3 θ α Γ α Δ α δ i j , λ i l l s = 1 Ψ s q i s ,
with
γ 0 α ρ α = 6 e α + p α ρ α p ρ α , ρ α p α γ 1 α = 2 2 p α ρ α + e α p ρ α , γ 2 α = 3 Γ α 5 3 p α γ 1 α 2 e α + p α ρ α Γ α + 5 3 p α + Ψ α p θ α ρ α c V α θ α , Ψ α = T α 2 3 γ 0 α θ α + 2 θ α 2 e α + p α ρ α p θ α , Γ α = 5 3 p α + ρ α p ρ α + θ α p θ α 2 ρ α c V α , β 1 α = p α θ α Ψ α γ 1 α 2 p α ρ α + e α .
The symbol c V α represents the specific heat at constant volume, defined as c V α = e α / θ α . We underline that these constitutive relations are in agreement with relations (5.58) p. 123 of Ruggeri and Sugiyama [11] since the structure for RBCs equations is equivalent to their two hierarchies system. We conclude that all functions introduced in the present analysis are determined except for the pressure p α , the internal energy e α and the production terms.

6.3. Production Terms

In this subsection, we will limit the generality of the production terms indeed we assume also for the production terms linear constitutive functions and in particular we apply the BGK assumption [35]. Therefore, the production term in each equation is proportional to the density of the same equation minus its equilibrium value. Explicitly we have:
s < i j > α = ρ < i j > α τ σ α , s l l s = 3 Δ s τ Δ s s i l l s = 2 q i s τ q s ,
where τ σ α , τ Δ s and τ q s are the well-known relaxation times that are functions of ρ α and θ α . In the next section we identify these relaxation times by comparison with the Gidaspow’s model [2].

7. Summary of the Consequences of the Physical Principles

We summarize here all consequences of the Galilein and entropy principles for the determination of the field equations. Our field equations are given by system (9) with the stress tensors in (8). The constitutive functions, which we calculated here, explicitly (limiting our analysis to linear terms) are written in (26) with the coefficients (28). The production terms are instead summarized in (29). The terms due to the external forces H i α coincide with those introduced by Gidaspow (3); The interaction terms in H i α and I α are given by (4) in terms of the coefficients β and γ .
Therefore, all is explicit in terms of the fields, when the functions p α and e α , the relaxation terms τ σ α , τ Δ s and τ q s and the coefficients β and γ in the interaction terms are known in terms of the fields. In the next section, we evaluate these last coefficients by comparison of our model with that of Gidaspow, where he obtained these function with the help of the kinetic theory [1].

8. Comparison with Gidaspow’s Model

In this section, we finally determine the remaining functions by comparison with the Gidaspow–Huang model [2]. First of all, we identify the RBCs pressure, that by calculations [1] in kinetic theory is given by
p s = 1 + 2 χ ρ s θ s with χ = 1 + e ϵ s g 0 ,
where ϵ s is the RBCs volume fraction, g 0 their radial distribution function, that can be chosen according to the Bagnold’s equations, and e is the restitution coefficient, a measure of the elasticity of the particle collisions (see [1,2,36] for more details). We inserted in Table 1 some values for these quantities taken in [2]. Furthermore, we assume that the RBCs energy density is given by e s = 3 θ s , so that c V s = 3 .
By use of the expression for the RBCs pressure (30), it is possible to evaluate all functions in (28) in terms of the RBCs density and temperature, which are
γ 0 s = 12 2 + χ 1 + 2 χ ρ s θ s 2 , γ 1 s = 2 5 + 4 χ θ s , Ψ s = 4 2 + χ 1 + 2 χ ρ s θ s 3 , Γ α = 1 3 1 4 χ 2 ρ s θ s , γ 2 s = 6 5 + 4 χ θ s , β 1 s = 1 2 1 + 2 χ 2 + χ .
In order to evaluate the relaxation times of the model, we substitute the constitutive relations (26) and the production terms (29) into the balance Equation (9) 4 , 5 and we proceed with the first Maxwellian iterations [10]. This method is really simple to be used, but it is very useful: In practice, we neglect from the left-hand side of Equation (9) 4 , 5 the non-equilibrium terms, that is, we consider the case in which no-rapid changes in the fields occur. The obtained equations become
2 p α v < i α x j > = ρ < i j > α τ σ α , 5 3 p s v k s x k = Δ s τ Δ α , 1 3 γ 0 s θ s 2 p s ρ s + e s p s θ s θ s x k = 2 q k s τ q α ,
which have exactly the form of the stress tensor, the dynamic pressure and the heat flux predicted by the Navier–Stokes and the Fourier laws (2). This method is particularly useful; indeed, we can firstly evaluate, by comparison of (2) with (32), our relaxation times in terms of the Gidaspow viscosities and heat conductivity. Furthermore, we make sure that our equations predict the same results of the classical model when no-rapid changes in the fields occur. From the comparison and taking into account the expression of the viscosities and heat conductivity obtained in [2], we easily obtain
τ σ s = 1 1 + 2 χ 4 5 + 10 96 π ϵ s χ 2 1 + 4 5 χ 2 χ d p π θ s , τ σ f = ε f μ f p f , τ Δ s = 4 5 1 1 + 2 χ χ d p π θ s , τ q s = 1 1 + 2 χ 2 + χ 1 + 75 384 π ϵ s χ 2 1 + 6 5 χ 2 χ d p π θ s ,
where d p represents the RBCs diameter, while ϵ f represents the fluid volume fraction with ϵ s + ϵ f = 1 .
Finally, the coefficients β and γ in the expression of the interaction terms I i f and I s in (9), are taken from the analogous terms in [2] so we set
β = 150 ε s 2 μ f 1 ε s d p 2 i f ε f 0.8 , 0 i f ε f > 0.8 , and γ = 12 1 e 2 g 0 d p π ρ s ε s θ s 3 2 .
The first relation takes into account only the linear term of [2] that is more appropriate for the linear constitutive equations presented in this paper.
All coefficients of our model are evaluated in terms of the field variables; furthermore, with the obtained functions, h is a convex function and the residual inequality (16) is verified.
Summarizing, system (9) with the relations (8), (26), (30) and (31), the productions terms (29) and (33), the terms due to the external forces in (3) and the interaction terms (4) and (34) is our system of field equations. It is a closed system of 24 field equations for the 24 field variables ρ α , v i α , θ α , ρ < i j > α , Δ s and q i s .

9. A Simple Numerical Example

In order to show clearly the applicability of our model, we present here a simple numerical example. We stress that these results are only a first simple application of our equations and they are not exhaustive of all the applicability of our model. They are also part of another paper [37], so an interested reader can found there more details.
We use the field equations herein obtained in order to describe the stationary blood flow in a vessel, which is idealized as the gap between two infinite parallel plates. We can assume that the flow is stationary and the fields depend only on the x-coordinate orthogonal to the two walls, so the complete set of field equations becomes a more simple set of ODE. We also assume that the RBCs and fluid velocities have the direction of the flow and they are driven by an external pressure, considered as external force.
It order to recover easily the solution, we linearize the ODE around an equilibrium state characterized by constant ρ 0 α , θ 0 α and vanishing fluxes. We illustrate in Figure 1 and Figure 2 the solution obtained for ϵ s = 0.24 , with the parameters in Table 1, and in terms of the dimentionless variables
x ^ = x L , θ ^ α = θ s θ 0 α , ρ ^ < i j > α = ρ < i j > α ρ 0 s θ 0 s , v ^ 3 α = v 3 α θ 0 s .
The solutions show that the RBCs temperature decreases in the center of the vessel, where instead the RBCs density increases. This result is in agreement with the Fahraeus–Linqvist effect, as already shown in [2]. The two velocities are illustrated in Figure 2, showing that the fluid anticipates RBCs, as they are carried from the fluid self.
The solutions in Figure 1 and Figure 2 are somehow less general than those presented in [2], since here only the 1D plane case is considered, no-radial or 2D geometry, nor slip at the boundary as in [2] and furthermore the equations are linearized around the constant state. Anyway, integration of this very simple and linearized system yields (see for details [37]) already the presence of some non-vanishing stress components. These non-vanishing components do not affect visibly the classical fields in Figure 1 and Figure 2, due to the strong assumptions in this chapter, but they will play surely an important role in all fields when all simplifying assumptions will be removed. This will be shown clearly in further research.

10. Conclusions and Final Remarks

In this paper, starting from a model for blood flow obtained by Gidaspow and Huang [2] in Classical Thermodynamics, we derived an hyperbolic model based on the methods proper of Rational Extended Thermodynamics [10,11,12]. The stress tensors and the heat flux are considered as new field variables and a new set of field equations is assumed. The introduction of new equations involves new quantities that are here determined as constitutive functions in terms of the fields. This was done using the Galileian and entropy principles. We determine the remaining functions with Gidaspow’s coefficients obtained in [1]. All parameters herein introduced are determined explicitly in terms of the fields and a simple numerical solution is shown in order to illustrate the applicability of our model.
The model herein obtained is hyperbolic, so it overcomes the paradox of infinite propagation speed of disturbances typical of parabolic ones. Our model can describe transient regimes better than the classical one and it offers a richer scenario of solutions as many other hyperbolic systems. Furthermore, by construction our model converges to the classical one [2] when the field gradients are negligible. On the contrary, we expect a better description of the phenomena in blood flow when strong changes in the blood occur. As seen for gases and mixture flow, we also expect that our Extended Thermodynamics model can describe more properly situations in with the domain presents more curvature of changes in the radial symmetry.
It must be said that, unfortunately, the integration of our field equations is more complicate than the classical ones, since our equations contain in practice the differential Equation (9) 4 , 5 instead of the classical Navier–Stokes and Fourier laws (2). This will imply a more complicate (normally numeric) integration of the whole set of equations and the introduction of initial values for stress tensors and heat flux.
Finally, our model has the same range of applicability of the classical model [2], so it can describe a great range of blood phenomena and cases. Some applications of our model are already under investigation and they will be part of further works.

Author Contributions

These authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Istituto Nazionale di Alta Matematica Francesco Severi: GNFM and Prin Progect: No. 2017YBKNCE “Multiscale phenomena in Continuum Mechanics: singular limits, off-equilibrium and transitions.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are grateful to the anonymous referees for their helpful remarks and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gidaspow, D. Multiphase Flow and Fluidization Continuum and Kinetic Theory Descriptions; Academic Press: Cambridge, MA, USA, 1994; ISBN 0-12-282470-9. [Google Scholar]
  2. Gidaspow, D.; Huang, J. Kinetic theory based model for blood flow and its viscosity. Ann. Biomed. Eng. 2009, 37, 1534–1545. [Google Scholar] [CrossRef] [PubMed]
  3. Huang, J.; Gidaspow, R.L.D. Pulsatile flow in a coronary artery using multiphase kinetic theory. J. Biomech. 2009, 42, 743–754. [Google Scholar] [CrossRef] [PubMed]
  4. Fedosov, D.; Dao, M.; Karniadakis, G.; Suresh, S. Computational Biorheology of Human Blood Flow in Health and Disease. Ann. Biomed. Eng. 2014, 42, 368–387. [Google Scholar] [CrossRef] [PubMed]
  5. Felicetti, L.; Femminella, M.; Reali, G.; Liò, P. A Molecular Communication System in Blood Vessels for Tumor Detection. In Proceedings of the NANOCOM’ 14: ACM The First Annual International Conference on Nanoscale Computing and Communication, Atlanta, GA, USA, 6–9 May 2014. [Google Scholar]
  6. Wu, W.; Yang, F.; Antaki, J.; Aubry, N.; Massoudi, M. Study of blood flow in several benchmark micro-channels using a two-fluid approach. Int. J. Eng. Sci. 2015, 95, 49–59. [Google Scholar] [CrossRef]
  7. Lopes, D.; Puga, H.; Teixeira, J.C.; Teixeira, S.F. Influence of arterial mechanical properties on carotid blood flow: Comparison of CFD and FSI, studies. Int. J. Mech. Sci. 2019, 160, 209–218. [Google Scholar] [CrossRef]
  8. Melka, B.; Gracka, M.; Adamczyk, W.; Rojczyk, M.; Golda, A.; Nowak, A.; Białecki, R.; Ostrowski, Z. Multiphase simulation of blood flow within main thoracic arteries of 8-year-old child with coarctation of the aorta. Heat Mass Transfer 2018, 54, 2405–2413. [Google Scholar] [CrossRef]
  9. Gidaspow, D.; Chandra, V. Unequal granular temperature model for motion of platelets to the wall and red blood cells to the center. Chem. Eng. Sci. 2014, 117, 107–113. [Google Scholar] [CrossRef]
  10. Müller, I.; Ruggeri, T. Rational Extended Thermodynamics; Springer: New York, NY, USA, 1998. [Google Scholar]
  11. Ruggeri, T.; Sugiyama, M. Rational Extended Thermodynamics Beyond the Monatomic Gas; Springer: New York, NY, USA, 2015. [Google Scholar]
  12. Ruggeri, T.; Sugiyama, M. Classical and Relativistic Rational Extended Thermodynamics of Gases; Springer: Cham, Switzerland, 2021. [Google Scholar]
  13. Liu, I.S. Method of Lagrange multipliers for exploitation of the entropy principle. Arch. Ration. Mech. Anal. 1972, 46, 131. [Google Scholar] [CrossRef]
  14. Barbera, E.; Brini, F. Heat transfer in gas mixtures: Advantages of an extended thermodynamics approach. Phys Lett. A 2011, 375, 827–831. [Google Scholar] [CrossRef]
  15. Barbera, E.; Brini, F. Frame dependence of stationary heat transfer in an inert mixture of ideal gases. Acta Mech. 2014, 225, 3285–3307. [Google Scholar] [CrossRef]
  16. Carrisi, M.C.; Montisci, S.; Pennisi, S. Entropy Principle and Galilean Relativity for Dense Gases, the General Solution without Approximations. Entropy 2013, 15, 1035–1056. [Google Scholar] [CrossRef]
  17. Arima, T.; Taniguchi, S.; Ruggeri, T.; Sugiyama, M. Extended thermodynamics of dense gases. Contin. Mech. Thermodyn. 2011, 24, 271. [Google Scholar] [CrossRef]
  18. Arima, T.; Carrisi, M.C.; Pennisi, S.; Ruggeri, T. Which moments are appropriate to describe gases with internal structure in Rational Extended Thermodynamics? Int. J.-Non-Linear Mech. 2021, 137, 103820. [Google Scholar] [CrossRef]
  19. Pennisi, S.; Ruggeri, T. Relativistic extended thermodynamics of rarefied polyatomic gases. Ann. Phys. 2017, 377, 414–445. [Google Scholar] [CrossRef]
  20. Arima, T.; Carrisi, M.C.; Pennisi, S.; Ruggeri, T. Relativistic Kinetic Theory of Polyatomic Gases: Classical Limit of a New Hierarchy of Moments and Qualitative Analysis. Partial. Differ. Equ. Appl. 2022, 3, 1–21. [Google Scholar] [CrossRef]
  21. Barbera, E.; Brini, F. A rational extended thermodynamics model for metal electrons in bounded domains. Ric. Mat. 2019, 68, 37–56. [Google Scholar] [CrossRef]
  22. Barbera, E.; Brini, F. New extended thermodynamics balance equations for an electron gas confined in a metallic body. Ric. Mat. 2021, 70, 181–194. [Google Scholar] [CrossRef]
  23. Trovato, M. Quantum maximum entropy principle and quantum statistics in extended thermodynamics. Acta Appl. Math. 2014, 132, 605–619. [Google Scholar] [CrossRef]
  24. Barbera, E.; Curró, C.; Valenti, G. On discontinuous travelling wave solutions for a class of hyperbolic reaction-diffusion models. Physical D 2015, 308, 116. [Google Scholar] [CrossRef]
  25. Consolo, G.; Curró, C.; Grifó, G.; Valenti, G. Oscillatory periodic pattern dynamics in hyperbolic reaction-advection-diffusion models. Phys. Rev. E 2022, 105, 034206. [Google Scholar] [CrossRef]
  26. Consolo, G.; Curró, C.; Valenti, G. Turing vegetation patterns in a generalized hyperbolic Klausmeier model. Math. Methods Appl. Sci. 2020, 43, 10474–10489. [Google Scholar] [CrossRef]
  27. Consolo, G.; Grifó, G. Eckhaus instability of stationary patterns in hyperbolic reaction–diffusion models on large finite domains. Partial. Differ. Equations Appl. 2022, 57. [Google Scholar] [CrossRef]
  28. Boillat, G. Sur l’existence et la recherche d’équations de conservation supplémentaires pour les systèmes hyperboliques. Compt. Rend. Hebd. Seances Acad. Sci. A 1974, 278, 909. [Google Scholar]
  29. Ruggeri, T.; Strumia, A. Main field and convex covariant density for quasi-linear hyperbolic systems. Relativistic fluid dynamics. Ann. Inst. Henri Poincaré 1981, 34, 65. [Google Scholar]
  30. Barbera, E.; Brini, F. Non-isothermal axial flow of a rarefied gas between two coaxial cylinders. AAPP 2017, 95, A1. [Google Scholar]
  31. Barbera, E.; Brini, F. Stationary heat transfer in helicoidal flows of a rarefied gas. EPL 2017, 120, 34001. [Google Scholar] [CrossRef]
  32. Müller, I. Thermodynamics; Pitman: London, UK, 1985; ISBN 0-273-08577-8. [Google Scholar]
  33. Kremer, G.M. Extended thermodynamics of mixtures of ideal gases. Int. J. Eng. Sci. 1987, 25, 95–115. [Google Scholar] [CrossRef]
  34. Pollino, A.; (University of Messina, Messina, Italy). Dissertation, 2024; in preparation.
  35. Bhatnagar, P.L.; Gross, E.P.; Krook, M. A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component systems. Phys. Rev. 1954, 94, 511–525. [Google Scholar] [CrossRef]
  36. Savage, S.B. Granular flow at high shear rates. In Theory of Dispersed Multiphase Flow; Meyer, R.E., Ed.; Academic Press: New York, NY, USA, 1983; pp. 339–358. [Google Scholar]
  37. Barbera, E.; Pollino, A.; (University of Messina, Messina, Italy). A Linear Simulation of Blood Flow in a Vessel, 2022. in preparation.
Figure 1. Dimensionless RBCs temperature (a) and dimensionless RBCs density (b) obtained in terms of x ^ with ε s = 0.24 .
Figure 1. Dimensionless RBCs temperature (a) and dimensionless RBCs density (b) obtained in terms of x ^ with ε s = 0.24 .
Mathematics 10 02977 g001
Figure 2. Dimensionless RBCs and fluid velocities in function of x ^ obtained for ε s = 0.24 .
Figure 2. Dimensionless RBCs and fluid velocities in function of x ^ obtained for ε s = 0.24 .
Mathematics 10 02977 g002
Table 1. Values of the parameters of the problems [2].
Table 1. Values of the parameters of the problems [2].
Tube diameter0.19 mm
Plasma density1020 kg m−3
Plasma viscosity0.0012 kg m−1 s−1
RBC size8 μm
Plasma density1092 kg m−3
Restitution coefficient0.95
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Barbera, E.; Pollino, A. An Extended Thermodynamics Model for Blood Flow. Mathematics 2022, 10, 2977. https://doi.org/10.3390/math10162977

AMA Style

Barbera E, Pollino A. An Extended Thermodynamics Model for Blood Flow. Mathematics. 2022; 10(16):2977. https://doi.org/10.3390/math10162977

Chicago/Turabian Style

Barbera, Elvira, and Annamaria Pollino. 2022. "An Extended Thermodynamics Model for Blood Flow" Mathematics 10, no. 16: 2977. https://doi.org/10.3390/math10162977

APA Style

Barbera, E., & Pollino, A. (2022). An Extended Thermodynamics Model for Blood Flow. Mathematics, 10(16), 2977. https://doi.org/10.3390/math10162977

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