Cardiovascular disease is a major threat to human health. The study on the pathogenesis and prevention of cardiovascular disease has received special attention. In this paper, we have contributed to the derivation of a mathematical model for the nonlinear waves in an artery. From the Navier–Stokes equations and continuity equation, the vorticity equation satisfied by the blood flow is established. And based on the multiscale analysis and perturbation method, a new model of the Boussinesq equation with viscous term is derived to describe the propagation of a viscous fluid through a thin tube. In order to be more consistent with the flow of the fluid, the time-fractional Boussinesq equation with viscous term is deduced by employing the semi-inverse method and the fractional variational principle. Moreover, the approximate analytical solution of the fractional equation is obtained, and the effect of viscosity on the amplitude and width of the wave is studied. Finally, the effects of the fractional order parameters and vessel radius on blood flow volume are discussed and analyzed.

1. Introduction

In recent years, in the field of biological rheology, the rheology which is related to the blood, blood vessels, and heart that constitute the human blood circulation has been developed rapidly. Blood circulation is a complex system that can cause great damage to the whole body when a certain organ is diseased. In order to better understand the physiological and pathological behavior of the human cardiovascular system, it is necessary to deeply study the dynamics of blood flow in the arteries and the mechanical factors of blood flow.

Soliton phenomena exist in many fields [13]. Organism is a completely nonlinear complex medium, and both the blood composition and the structure of blood vessels show obvious nonlinear characteristics. The nonlinearity of blood flow has long been discovered by Womersley [4, 5] and McDonald [6, 7], which also provides a new direction and way for people to understand the law of life movement. Since then, many scholars have begun to study the field and made great progress. Ravindran et al. [8] derived the nonlinear Schödinger equation (the terms of pseudodifferential operators) governing the modulation of periodic waves. The KdV equation describing blood flow is obtained by Sigeo [9], to explain the steepness of pressure waves during propagation. Hashizume [10] analyzed the propagation of pressure waves from a theoretical perspective. Liu [11] combined arterial flow as a balanced flow with periodic small pulsatile flow and analyzed the effect of vascular elasticity on blood flow. Demiray [12, 13] considered propagation of wave through a viscous incompressible fluid contained in a prestressed thin elastic tube. Choy [14] deduced the mathematical model of nonlinear wave modulation of artery with stenosis.

Fractional derivative theory and methods [1518] are widely used in the study of nonequilibrium systems of various intermediate processes and critical phenomena in physics and mechanics, especially in nonlinear science [1922]. Fractional differential equations are transformed in a standard differential equation by replacing the time derivative or the space derivative with the fractional derivative. Compared with integer-order differential equations, the most important advantage of fractional derivative equations is that it can better fit some natural physical processes and dynamic system processes. Moreover, the study of solving partial differential equations also has a new exploration, such as the extended tanh method [23], the exp-function method [24], the variational-iteration method [15, 25], the Hirota bilinear method [19], and the -expansion method [26].

The structure of the full article is as follows: in Section 2, the Boussinesq equation with viscous term is derived by the multiscale analysis and perturbation method and used for the first time to describe blood flow. Based on the new model, we obtain the time-fractional Boussinesq equation with viscous term in Section 3. In Section 4, the approximate analytical solution of the above fraction equation and the viscous effect is discussed. Finally, the effects of fractional order, vascular radius, and blood flow velocity on stroke volume are analyzed and studied.

2. Derivation of the Boussinesq Equation

As we all know, the cardiovascular system is a complete closed conduit system. Thus, the blood can be considered as an incompressible non-Newtonian fluid. When dynamic equilibrium of the blood is disturbed by a pressure pulse generated by the motion of the heart, a harmonic wave type of motion will be developed in the blood. Although previous researchers have done some research about this question, the viscosity of fluid is often ignored. We consider blood vessel as cylindrical shape and adopt cylindrical coordinate system to describe the motion of blood. In this paper, we set as the central axis of a blood vessel and axial coordinate and set as the radial coordinates. Therefore, based on Womersley theory [4], basic equations of pulsatile flow in arteries are described as follows:where and are the axial and radial velocities of blood, respectively; is the atmospheric pressure; is the blood density; and represents the viscosity coefficient of blood. From (3), we can define the flow function of blood which satisfies

Substituting (4) into (1) and (2), we can gain the vorticity equation of blood flow with viscosity term:where . Through analysis, we assume that the basic law of blood flow can be expressed by the Poiseuille flow. Therefore, the speed of blood can be expressed aswhere is the maximum velocity of central axis of the blood vessel. represents the length of the blood fluid, and represents viscosity coefficient of blood. is the pressure difference between two ends of blood, and is the radius of the blood vessel. Next, when considering the viscosity of the vascular wall, we can get the revised form of the Poiseuille flow:

We set basic flow function of blood as . Therefore, it meetswhere we define disturbed flow function which represents disturbance of blood flow. Therefore, the flow function can be expressed as

And then, from (8) and (4), we can obtain

Substituting (9) and (10) into (5) and considering that there are other effects of dissipation which can offset , we obtainwhere and . The above equation is a vorticity equation of disturbed flow function which contains information about the effects on the heart when physiological or psychological condition changes dramatically.

In the previous studies, researchers have obtained the KdV equation from the Navier–Stokes equation and continuity equation of blood flow. However, with the development of nonlinear theory increasing, more equations which have stronger nonlinear properties are needed to describe problems in reality. Therefore, the Boussinesq equation model will be deduced in this paper. In order to obtain this equation, we firstly define coordinate transformations with long wave approximation and adopt the space-time transformation as follows:and we setwhere is a small parameter. Substituting (12) into (11) yields

Then, we extend perturbed flow function about as follows:

Due to the fact that the velocity of wave propagation in blood is much faster than that of blood flow on the axis of the blood vessel, we can obtain . Thus, substituting (15) into (14), we can yield all levels of approximation equations about :where .

For , it is easy to see that has a solution, and the form is ; then, we have

For , by analysis, we assume , where satisfies the following equation:

Let us consider the equation of order . By substituting , , (17), and (18) into the equation and multiplying the both sides of the order equation by , we have

And then integrating it with respect to from 0 to , we can have the following equation:where

Equation (20) is a new Boussinesq equation with viscosity effect which can be used to describe the state of the disturbance flow of blood with viscous properties. Comparing with the previous KdV model, the above model has stronger nonlinear properties. Besides, it can be concluded that there is disturbance flow in the form of solitary waves in the arteries.

3. Derivation of the Time-Fractional Boussinesq Equation

In this section, we seek for the time-fractional Boussinesq equation by using the semi-inverse method and the fractional variational principle. Firstly, we use a potential function , where gives the potential equation of (20) in the following form:

The functional of the potential equation (22) can be represented bywhere , , , and are the Lagrangian multipliers which can be obtained later by taking the variation of (23). and are considered as the fixed functions.

Making use of integrating (23) by parts and assuming , we can gain

Applying the variation optimum condition, we can derive the Euler equation as follows:

Comparing (25) with (22), we can obtain the following Lagrangian multipliers:

Therefore, the Lagrangian form of (20) can be given by

At this time, the Lagrangian form of the time-fractional Boussinesq equation can be represented aswhere based on the left Riemann–Liouville compression derivative, the fractional derivative is defined as

Therefore, the time-fractional form of (20) can be expressed by

Thus, the variation of functional (30) can be obtained:

As we all know, the fractional integration by parts is obtained by the following rule:based on the right Riemann–Liouville fractional derivative, is obtained by

Using (31) and (32), we get

By using variation principle, when , we obtain the optimization of the variation of functional. And then, we yield the Euler–Lagrange equation of (20) as the following equation:

Substituting (28) into (35) and letting , (35) is turned into the following form:

Equation (36) is a new model, which is gained by the time-fractional method. Therefore, it is named as time-fractional Boussinesq equation.

4. Solution of the Time-Fractional Boussinesq Equation

In the previous section, a new time-fractional Boussinesq equation had been gained. It can be used to describe the flow characteristics of blood with viscous properties. Due to this reason, we decide to yield the solution of the time-fractional Boussinesq equation to further characterize the flow state of blood.

Let we define the fractional complex transform as follows:where is the unknown constant. And then, we can yield the fractional derivatives into classical derivatives with (37):

It is obvious that (36) can be deduced to

Next, we consider to gain the solution of (39). Due to the fact that describes viscous properties of blood, we set as a small parameter. Assuming that and considering to get the properties of (39), two time scales are defined as follows:and is expanded as

Thus, we obtain all levels of approximation equations about :

For , it is obvious that the solution of the equation can be expressed aswhere and represents the maximum amplitude at initial moment. Then, for , we use the equation to determine the form of . Assuming that

By substituting (44) into the order one of equation from (42), the following equation is obtained:where . In order to describe the structure of , we consider to choose proper function which should be orthogonal to as

Multiplying both sides of (45) by , integrating them by parts and using and , we can see that satisfies

By assuming , the solution of (47) can be obtained as

Substituting (48) into (46) and using (45), we obtain

Therefore, we obtain the approximate solution of (39) as follows:where . The above equation is the solution of (39) and can be used to describe the flow state of blood more specifically.

Finally, the solution of the time-fractional Boussinesq equation iswhere .

In order to study the influence of blood viscosity and fractional order on the evolution of blood flow, we change the value of , , , and to get related results. From Figure 1, we can see that when viscosity parameter , the amplitude of solitary waves is largest; at the same time, the width of solitary waves is the smallest. With the decrease of the amplitude, the width increases gradually. Figure 2 depicts the evolution of solitary waves in blood under different viscosity coefficient . It is evident that amplitude of the solitary waves becomes smaller with time from Figure 2(a), while the width of the solitary waves becomes larger with time . And we can know that the viscosity effect leads to the amplitude of the solitary waves to decrease and the width of the solitary waves to increase. In addition, in contrast with Figure 2(b), we find with the increasing viscosity coefficient , the rate of amplitude reduction has been accelerated and the direction of the wave changes from the positive direction of to the negative direction of it. Figure 3 shows the solitary waves under different fractional order . We find the peak value of amplitude is not linear with , and there is a minimum of the peak with the corresponding critical value of . According to the critical value, the peak value of amplitude increases firstly and then decreases.

5. Analysis and Discussion

In this section, we consider to calculate the value which can represent the state of blood by using the results from the previous sections. First, we search for the solution of (17) and (18). Taking and substituting them into (17) yield

And the solution of (52) is expressed aswhere are the constants and and are expressed as

Thus, the solution of (17) can be described as

Similarly, by substituting (55) into (18), we obtainwhere is a constant. We take the physiological parameters of the human body as follows:

And set , , and . Therefore, by using (21), we obtain the coefficients of (36) as follows:

Due to that the basic flow in a human artery can be assumed to be the Poiseuille flow, the stream function can be expressed as

So the stream function can be obtained as follows:where and following (12).

Because in the cylindrical coordinate system, the difference of the flow function values of any two points in the plane passing through the symmetry axis multiplied by is equal to the value, which is the flow rate of the rotating surface formed by the rotation of any curve passing through these two points around the symmetric axis, we decide to calculate blood flow by using the values of the outlet and entrance of blood vessels. We take the flow function of the central position of the vessel at the entrance of the aorta as , where and , and take the flow function at the junction of heart and aorta as , where and . Thus, we obtain

Thus, the blood flow at the outlet of the heart can be expressed as

In a cardiac cycle, the time-dependent image of the blood flow volume of the cardiac outlet is shown in Figure 4.

From Figure 4, we can see that, when is smaller, this model can be used to describe people with less blood flow or faster blood flow conduction. And when is large, this model can be used to describe people with slow blood flow in the body, or blood output from the heart is not easy to circulate throughout the body. These people may have lesions in some parts of the blood vessels, causing obstruction to blood flow. Or they are in good physical condition, but they are at high risk of cardiovascular and cerebrovascular diseases. Furthermore, we can find the appropriate value of to describe different populations by combining the actual monitoring data.

As we all know, stroke volume of healthy people is basically stable. When stroke volume is not stable for a long time, we can infer that there may be something wrong with that person’s cardiovascular system. If the amount of stroke volume is smaller, it will lead to insufficient oxygen supply to other tissues and organs. On the contrary, when stroke volume is too large, it will increase the load of human blood vessels, thus causing some cardiovascular and cerebrovascular diseases. According to the above introduction, we can find that stroke volume in each cycle of the heart is an important index which can reflect the health status of human cardiovascular system. Thus, we give the description of the stroke volume as follows:where represents the pulsatile cycle of the heart. By calculating, with and , we obtain the stroke volume in each cycle of the heart as which is close to the theoretical calculation value . Therefore, this model can well reflect the actual blood flow of the human body.

Considering that there is always a certain degree of difference in genetic characteristics of different people, it is easy to misjudge the disease condition if we use a unified health standard to measure the population with different characteristics. Thus, we consider to yield several different standards which can be used to measure the health of different groups of people by adjusting the value of fractional order . Figure 5 shows that the stroke volume increases with the increase of . has a larger growth rate in the range of value from 0.4 to 0.8, while it grows slowly in other value ranges. Thus, we can consider selecting the appropriate value of from 0.4 to 0.8 with the actual test data.

Because the radius of the blood vessel can influence stroke volume, we consider changing the value of to get several results as shown in Figure 6.

In Figure 6, the values of from bottom to top are 0.2, 0.4, 0.6, and 0.8. We can find that under the condition of the same vessel radius, stroke volume increases with the increase of the maximum velocity of blood vessel center. In addition, under the same , the value of is also related to the size of the vascular radius. As the vessel radius increases, stroke volume increases as well.

6. Conclusion

In this paper, we derive time-fractional Boussinesq equation with viscous term to describe blood flow and discuss the effects of fractional order and viscosity on blood flow volume. This indicates that local hemodynamic factors play an important role in the occurrence and development of cardiovascular disease. Therefore, it is of great significance to study the changes of blood flow in blood vessels to reveal the pathogenesis of major vascular diseases.

Data Availability

No data were used to support the findings of the study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.


This work was jointly sponsored by the National Key Research and Development Program of China (2017YFC1502102) and National Natural Science Foundation of China (41675102).