Abstract

Biomathematics is an interdisciplinary subject consisting of mathematics and biology, which is widely applicable for the analysis of biological problems. In this paper, we provide a mathematical model of two-phase hepatic blood flow in a jaundice patient’s artery. The blood flow is thought to be a two-phased process. The clinical data of a jaundice patient (blood pressure and hemoglobin) is gathered. To begin, hemoglobin is transformed into hematocrit, and blood pressure is turned to a decline in blood pressure. For the examination of hepatic arteries in Newtonian and non-Newtonian movements, a mathematical model is constructed. The relationship between two-phase blood flow flux and blood pressure reduction in the hepatic artery is established. For various hematocrit levels, the blood pressure decrease is determined. The patient’s states are defined by the slope of the linear relationship between computed blood pressure decrease and hematocrit.

1. Introduction

Due to globalization and the multidisciplinary development of life science and mathematical science, mathematical modeling becomes a topic of interest for the scientific community. Mathematical modeling plays a significant role to understand the intricacies of infectious diseases [13]. It helps to study the mechanisms responsible for observed epidemiological patterns, assesses the efficiency of control policies, and predicts epidemiological inclinations. Nowadays, human beings suffer from many disorders due to global warming, pressure, and stress to lifestyle [4]. Some physiological elements that occur during exposure to weightlessness may include alteration in blood flow to the liver. So many disorders can be seen like diabetes, blood pressure, and other diseases estimation of hepatic blood flow.

Bio-mathematics is an emerging and dynamic field. The term bio-mathematics refers to the use of quantities and mathematical methods to solve biological problems [1, 2]. The bio-mathematics is an interdisciplinary subject to understand the mechanical properties of living tissues, anatomy and physiology in health and illness, and introductory biomechanics: from cell to organisms, blood flow, and microcirculation [57]. Bio-fluid mechanics is a well-established branch of bio-mathematics with the help of its normal functions, changed due to alternation via mathematical analysis [7].

It is well known that the heart circulates blood in the body using elastic tubes: the arteries, capillaries, and veins. A proper flow of blood is necessary for good health as it circulates oxygen and other essential nutrients to various parts of the body. Two main factors that affect the blood flow are known as the blood vessels and the properties of blood. Therefore, the modeling of blood flow is very complex as compared to other fluid flows. Blood is a non-Newtonian fluid; hence its viscosity is not constant, resulting in nonlinearity of stress and strain rate relation [8]. However, this fact can be approximated by the power law or Casson fluid model [9]. The walls of blood vessels are elastic and disobey Hooke’s law. The curvature of the wall can also affect the properties of blood flow. Due to the composite elastic materials of the wall, a nonlinear stress-strain relation exists. Blood flow is pulsatile resulting in the periodic nature of pressure gradient and the fluid velocity [10]. It is caused due to the beating of the heart. The radii of the tube are variable and can be further divided into multiple veins. Most of the flow is not fully developed due to the high Reynolds’s numbers and only inlet flow occurs. Blood is a nonhomogeneous, unsteady, and pulsatile flow that flows from the elastic tube which is branched repeatedly.

The first attempt to understand the mechanism of blood flow through various types of vessels was made by Leonardo da Vinci and Descartes [11]. Later, Aird and Silverman explained the nature of the circulation of blood in the cardiovascular system [11, 12]. In 1899, the first mathematical model of propagation of arterial pulse waves is given by Frank [13]. Sharan and Popel [14] explained the execution of the two-phase model for the blood flow in slim tubes and explained how the powerful extended viscosity near the wall affects blood flow [14]. The viscosity of human blood concerning its performance ability in a high static magnetic field is discussed by Haik et al. [15]. Formaggia et al. [16] have studied the one-dimensional nonlinear system, which defines the blood as compliant with arteries [16]. They used the Navier–Stoke′s equation, which shows the algebraic relation between intramural pressure and vessel. Jiang et al. [17] concentrated on the physical demonstration and propelled reproductions of fuel-fluid two-stage stream streams in atomization and splashes [17]. It was proposed that the large blood vessels behave like a Newtonian. The proposed model calculates the parameters of blood flow have no restrictions on the vessel size [18]. The precarious and incompressible stream of non-Newtonian liquid through composite stenosis was also considered [19]. The two phases Bingham model was proposed by Ahnert et al. [20, 21]. They used the differential equation for solutions and solve the mathematical problems. They also used the drift-flux model which defines the behavior of the model like the resistance of the flow and also use the different types of parameters. Kumawat et al. [22] mathematically analyze two-phase blood flow through a stenosed curved artery with hematocrit and temperature dependent viscosity [22]. A patient-specific artery geometry in the presence of stenosis (plaque) was considered. In 2016, Achaba et al. [23] studied the blood flow in arteries through a non-Newtonian viscosity model [23]. They define the two major difficulties. The first is a constitutive equation. No one model accepted the behavior of blood viscosity, and the second is the highly nonlinear equation for blood motion. They used a two-dimensional equation. The study reveals that the power law model is better for non-Newtonian blood flow. Mekheimer and their research group did extensive work on the nanoparticles drug delivery to blood hemodynamics in diseased organs [2428]. Sharma et al. [29] proposed entropy analysis of thermally radiating MHD slip flow of hybrid nanoparticles (Au-/Blood) through a tapered multistenosed artery [29]. Bhatti and Abdelsalam [30] studied the peristaltic propulsion of hybrid nanofluid flow with Tantalum and Gold nanoparticles under magnetic effects [30].

In this article, we consider the hepatic blood flow in arterioles with respect to the nature of the hepatic circulatory system in humans. The Herschel–Bulkley non-Newtonian model in Bio-fluid physiological is investigated.

2. Important Formulations

In this section, we discuss the meaning of Newtonian and non-Newtonian flow. The mathematical formulation can be found in our previous articles [3135].

2.1. Covariant Vectors

A set of n function of n coordinates in a coordinate system are said to form the components of covariant vectors if they transform to another coordinate system according to the following rule [36]:

2.2. Contravariant Vectors

A set of n function of n coordinates in a coordinate system are said to form the components of covariant vectors if they transform to another coordinate system according to the following rule [36]:

2.3. Tensors of Second Order
2.3.1. Covariant Tensors of Order Two

A set of functions of coordinates in a coordinates system are said to form the components of a covariant tensor of order two if they transform to another coordinate system as follows [37]:

2.3.2. Contravariant Tensors of Order Two

A set of n2 functions of n coordinates in a coordinates system are said to form the components of a contravariant tensor of order two, if they transform to another coordinate system as follows [37]:

2.3.3. Mixed Tensors of Order Two

A set of n2 functions of n coordinates in a coordinates system are said to form the components of a mixed tensor of order two, if they transform to another coordinate system as follows [37]:

2.4. Christoffel’s Symbols
2.4.1. Christoffel Symbol of First Kind

The Christoffel 3 index symbol of the first kind are denoted by and defined by the following equation [37]:

2.4.2. Christoffel Symbol of Second Kind

The Christoffel 3 index symbol of the second kind is denoted by and defined by the following equation [37]:

3. Modeling of Blood Flow through Vessels

3.1. Numerical Analysis for Hepatic Arteries

The power law equation of continuity is written as follows [38]:

The motion equation is expressed as follows [39]:where is power law constitutive equation. The blood density equation iswhere is the volume ratio of blood cells and is hematocrit. The viscosity of a mixture of blood is expressed as follows:

In cylindrical form,

Tonsorial form in cylindrical coordinates,

The conjugate metric tensor is as follows:

The nonvanishing Christoffel’s symbols of 1st kind are as follows:

And, all other Christoffel’s symbols of 1st kind are zero. The nonvanishing Christoffel’s symbols of 2nd kind are as follows:

And, all other Christoffel’s symbols of 2nd kind are zero. Relation between contravariant and component of velocity of-blood flow is as follows:

The component of: . The components of shearing stress tensor,

The covariant derivative of is

3.2. Solution for Newtonian

The blood flow in the artery is symmetric with respect to the axis. Hence, , and do not depend upon , also . Since only one component of velocity, which is along the axis is effective. Now . Flow is steady, then obtain the following equation [39]:

Equation of continuity reduces to

Equation of motion reduces to

Z-component,

Let us assume that the pressure gradient of blood flow in arteries is .

Integration gives us,where are boundary conditions and is constant.

Integration of equation (26) is as follows:

Again, using the second boundary condition , then from equation (27), we get the constant of integration value B is

Putting the above-given value of B from equation (27), we obtain the following equation:where R is the radius of artery. If is the flux through the artery tube, then

Therefore,

Now,

Integration of both sides with limit initial to final,

3.3. Solution for Non-Newtonian

Now

The blood flow,

Equation of continuity reduces to,

Equation of motion reduces to,

Z-component,

Let us assume that the pressure gradient of blood flow in arteries is . Then,

We arrive at,

Let us consider boundary conditions; , and . We arrive at

After integration,

Under no slip boundary condition (, ), equation (41) can be written as follows:

By inserting the value of B in (41),

Equation (43) describes the velocity of blood flow in arteries. Total flux flow of blood () through a tube of arteries is defined as follows:

Thus,

Equation of motion,

Integration of equation (46) gives us the pressure drop of blood.

4. Analysis for Hepatic Arteries

4.1. For Newtonian Motion

If Q is the flux through the tube, then

Using (21), equation (48) can be written as follows:

Now,

Integration of both sides with limit initial to final,

Table 1 demonstrates the measured hemoglobin and standard BP of a jaundice patient (name: A, age: 43, and sex: F). Firstly, the hemoglobin is converted to hematocrit using the relation mentioned in Table 2. After that, the BP drop is assessed from standard BP using the following relation [31]:where is systolic and is diastolic. Later, the BP drop is converted to BP (in ). Let us consider that and blood pressure drop is 4332.9. Equation (11) gives us,

The relation between flow flux and blood pressure drop of two-phase blood flow in the hepatic artery is expressed as follows:

Using (10) & (63), the solution of equation (54) is obtained as follows:

Different values of gives us the blood pressure drop (Table 1). Figure 1 demonstrates the BP drop as a function of hematocrit. It is observed that the BP drop is increased as hematocrit enhances, which means that the hemoglobin of the patient is normal.

4.2. For Non-Newtonian Motion

The total flow-flux of blood through a tube of the arteries is Q defined by the following equation:

Both sides take power n,

We know that the for non-Newtonian motion. Hence,

Integration of both sides with limit initial to final,

Pressure drop of blood,

Let us consider that H = 0.02858 and BP drop = 4332.90 and

Using equation (11),

Again, by using equation (11), we arrive at,

Rearrangement of (59) provides the blood flow flux ,

Solution of equation (64) gives the value of ,

Inserting the values of , , , and in (60), we arrive at,

By inserting the value of in equation (66), we get a BP drop (Table 3).

Figure 2 demonstrates the relation between the calculated BP drop and hematocrit for non-Newtonian motion. It is observed that the BP drop increases as the hematocrit increases.

Figures 3(a) and 3(b). Mathematical data of BP drop (calculated using (66)) and BP (in ) as a function of Hematocrit for non-Newtonian motion. Solid lines represent the linear fitting.

Rheological properties of the Reiner–Rivlin fluid model for blood flow through a tapered artery with stenosis have been studied by Akbar et al. [45]. Elogail and Mekheimer [46] implemented a numerical study that simulates blood flow through a microvessel involving oxytactic microorganisms and nanoparticles [46]. The oxytactic microorganisms exhibit negative chemotaxis to gradients of oxygen (oxygen repellents). These microorganisms are to batter infected hypoxic tumor cells as drug carriers [46]. Awad et al. [47] studied the flow of a non-Newtonian fluid with nonzero yield stress [47]. Navier stokes equation is used to simulate this subject mathematically. The elasticity on the stenosis arterial walls is simulated by Rubinow and Keller model and the Mazumdar model [47]. Kumawat et al. [22] mathematically analyze two-phase blood flow through a stenosed curved artery with hematocrit and temperature dependent viscosity [22].

Last but not least, we developed mathematical for blood flow in arteries during jaundice. We collected the blood pressure data of jaundice patients from the hospital and later applied a mathematical model to analyze the data.

5. Conclusion

We developed a deterministic model of the link between hematocrit and blood pressure fluctuations in jaundice patients in this study. The mathematical analysis is validated on the clinical data and calculated BP values of the jaundice patient (name: A, age: 43, and sex: F). The role of blood flow and hemoglobin malfunction is investigated in this research. It is feasible to propose a patient for better therapy based on the trend line of the association between hematocrit and BP decline. With the use of model Newtonian motion and non-Newtonian motion, if the trend line exhibits positive slope, then the patient’s hemoglobin is normal. If the trend line indicates a negative slope, however, the management of jaundice patients should be modified. This study is useful to predict the hemoglobin status of jaundice patients on the basis of blood pressure measurement.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Conceptualization was done by A. S. and R. A. K.; data curation was done by A. S., R. A. K., and T. A.; formal analysis was done by A. S. and R. A. K.; investigation was done by A. S., S. K., and R. A. K.; methodology was done by A. S. and R. A. K.; project administration and validation was done by A. S. and R. A. K; visualization was done by A. S., T. A., and R. A. K; writing-original draft was done by A. S., S. K., T. A., and R. A. K.; writing-review and editing was done by A. S., S. K., and R. A. K. All authors read and approved the final manuscript.