Abstract
This article discusses a semi-analytical method for solving modelling equations for blood flow through compliant vessels. He's method, also called the variational iteration method (VIM), is used to find an approach to a solution in an explicit form. The chosen problem has been handled by an asymptotic reduction of the associated incompressible Navier–Stokes equations, to obtain a one-dimensional nonlocal blood flow model with source terms. Numerical examples were chosen based on experimentally valued parameters and were used to highlight and illustrate the benefits of VIM. The results obtained from the blood flow model were compared to a numerical reference method.
Keywords:
1. Introduction
Differential equations play a crucial role in modelling natural phenomena in reality. Scientists deal with many properties and actions (collisions, chemical bonding, solitary wave, viscosity, temperature, pressure, etc.) making the equations that govern these problems increasingly difficult to solve it analytically [Citation1–6]. For this purpose, researchers have resorted to finding solutions using many methods to approach it, as the direct algebraic methods [Citation7–9], rational expansion method [Citation10], reductive perturbation technique [Citation11], extended mapping [Citation12], Adomian decomposition method, homotopy perturbation method, differential transform method, and VIM [Citation13–18]. VIM, for example, is a semi-analytical technique among these mentioned methods. It was invented by He at the end of the last century [Citation19,Citation20]. Many authors have developed VIM in the last decades, as Wazwaz [Citation21,Citation22], which has proposed to solve various types of nonlinear problems. The VIM solution is treated by an iterative process leading to an accurate solution. The advantage of this method is that it can provide an approximate solution in an explicit form, unlike classical numerical methods which can only provide it by interpolation. In this study, we treat the problem of blood flow behaviour through compliant vessels using the VIM. It is therefore necessary to reduce the number of variables involved in this problem. This approach is based on an asymptotic reduction process applied to the associated incompressible Navier–Stokes equations. It should be noted that recently this subject has been treated using effective numerical methods by certain researchers, one can quote for example [Citation23,Citation24]. However, in this work, we present a way to use the advantages of VIM to find a solution approach in an explicit form without using a linearization or discretization procedure. We suggest applying this method to study the blood flow in the artery because this study allows us to know the blood flow by measuring the blood pressure and determining its flow in the arteries. Blood flow is assumed to be Newtonian and is governed by the equation of continuity and the equilibrium equation of momentum (which are known as the Navier–Stokes equations). We apply VIM to the equations governing blood flow in arteries in the presence of a source term. We proposed a comparison between the VIM solution with a numerical method based on the MATLAB PDEPE function to evaluate this proposed method.
2. Global one-dimensional blood flow modelling
In recent years, a lot of works have proposed mathematical and numerical modelling of the human cardiovascular system. These models are based on Navier–Stokes equations [Citation26,Citation27]. The 3D model describes the flow motion in compliant vessels and its interaction with wall displacement [Citation28,Citation29]. To simplify these equations, many authors have studied a monodirectional flow system [Citation25,Citation29,Citation31], where the 1D flow equations are generated from a single vessel system. This model was used to study systematic circulation since their activity is based on reasonable calculation costs. This model describes an incompressible viscous fluid flow in an elastic tube taking account the blood vessels and their interaction with wall displacement. Several approaches can be used to derive the equations by using the asymptotic analysis of incompressible Navier–Stokes equations in long and narrow channels [Citation29,Citation30]. Two approaches are distinguished. The first integrates the Navier–Stokes equations in a generic section. The second is based on adimensional underlying equations [Citation30,Citation31]. This approach seems more realist and gives global results.
2.1. Asymptotic reduction
In the asymptotic reduction, the equations will be treated with nondimensional variables, for more details see [Citation30]. Basing on that the artery or aorta is an almost cylindrical vessel and the blood flows in the axial direction. Then, the system of equations used must be in cylindrical coordinates , when the x-axis is aligned with the symmetry axis of the vessel. We denote by the velocity. We will first assume that the angular velocity is zero to obtain the following equations of motion: (1) (1) (2) (2) Translating the reduced nondimensional equations and using the incompressibility condition, we obtain (3) (3) We consider the characteristic quantities and that are the radial and axial characteristic velocity, L is the characteristic length, is the ship's inner radius, and the following appropriate dimensionless variables: By a simple calculation [Citation30], the first reduced motion equation will be as follows: (4) (4) The second reduced motion equation will be as follows: (5) (5) This equation implies that the pressure is constant across the cross-section of the vessel. Finally, the reduced incompressibility condition is as follows: (6) (6) In this step, we will express the equations in terms of the average quantity in the transverse zone. Let be the inner ship radius.
Therefore, (7) (7) (8) (8) where is the axial mean velocity, and α is the correction term.
We integrate the equations on to which expresses the terms of the average quantity by specifying the boundary condition with : For the incompressibility condition, we integrate Equations (Equation6(6) (6) ) and using the definition of , then For the first equation of motion, we integrate the first reduced motion equation, we have (9) (9) The dimensionless reduced and average equations are as follows: (10) (10) (11) (11) Finally, we define the average cross-sectional velocity U and the coefficient α, then we find their relation with dimensionless quantities (12) (12) where R the internal ship radius in dimensional terms.
The dimensional axial velocity is (13) (13) with even for (14) (14) Using these definitions, the reduced equations are transformed into their dimensional form to obtain (15) (15) (16) (16) To obtain the equations written according to the average magnitude, one can take a typical approximation of the speed profile [Citation30] (17) (17) Equations (Equation14(14) (14) ) and (Equation17(17) (17) ) involve the following relationship between the speed profile form determined by γ and the correction coefficient α: By posing and denote the volumetric flow, we obtain (18) (18) where is the mass source per unit length, and represents the acceleration of the flow due to several forces which are applied on the fluid (gravity, friction, etc.) [Citation29], and is the viscous flow resistance per unit length of the vessel.
3. VIM solution of blood flow problem
The variational iteration method [Citation21,Citation22], which is a modified and generalized method of the Lagrange multipliers [Citation19], demonstrates an effective, easy, and accurate resolution of a large class of nonlinear problems with approximations that converge quickly to a precise solution. To achieve our goal in the management of systems of nonlinear partial differential equations, we use VIM, and to illustrate this method we will consider the following system which writes in the form of an operator: (19) (19) where and are linear and nonlinear operators, respectively, h and g design the source terms. We can build the following functional correction: (20) (20) (21) (21) The Lagrange multipliers can be recognized using the variational theory, are appropriate initial functions, the index n indicates the nth-iteration , and are considered as restricted variations, i.e. [Citation14,Citation22]. The successive approximations solutions can be easily obtained.
For the system Equation (Equation18(18) (18) ), we have . Then, the VIM method to solve the blood flow problem can be written as follows: (22) (22)
4. Numerical results
We will present in this section some examples to show the effectiveness of the described method.
4.1. Example 1
The first example is chosen with linear source terms to illustrate the VIM for a 1D-blood model. To verify the accuracy of the proposed solution, we relate the numerical results to a numerical solution using MATLAB's PDEPE function. Let the following parameters: We consider the initial condition (23) (23) Let , and , then the problem of the blood flow model (Equation18(18) (18) ) is as follows: (24) (24) Using Equations (Equation22(22) (22) ), we get the following recursive process: (25) (25) The first approach of VIM is written by The second The approximation of averaged over the cross-section linear velocity is given by . It is a two-step VIM solution for the blood flow model.
Figure shows the VIM solution (solid line) with the numerical solution (solid line with circle markers). Both solutions are almost identical. We note that the maximum relative error is on the order of (Table ).
4.2. Example 2
In this example, we use a model that takes into account the elastic properties of the arteries described in [Citation33]. This model is characterized by the following pressure expression: (26) (26) where , and is the velocity of small disturbances propagation along the vessel.
Let the initial condition (27) (27) Using Equation (Equation27(27) (27) ), we get .
The first approximation is given The second approximation With
Figure illustrates the variation of the axial velocity for the time and the position of the cross-section , for a the blood flow model taking into account the elastic properties of the arteries. Figure shows the behaviour of the axial velocity as a function of time. Knowing that we took different fixed values of x for each solution to see the effect of the position of the artery cross-section on the axial velocity.
4.3. Example 3
We consider the mechanical model treated in [Citation29,Citation32,Citation34] and describes the elastic properties of both arteries and veins given by the following pressure: (28) (28) where depends of Poisson's ratio, Young's modulus, and the wall thickness to radius ratio.
Considering the following source terms: For the numerical study, we put and n=0.
Let the following initial condition: We get the first step of VIM The second step is given by With
Figure illustrates the change in axial velocity as a function of time t which varies in the interval and the position of the cross-section with .
5. Discussion
In Example 4.1, the performance of the VIM solution is compared to a numerical method using the MATLAB PDEPE function. Note that this example is entirely illustrative and based on simple expressions of pressure and source terms. The outcomes show that after only two iterations of VIM, we obtain a reasonable result. Examples 4.2 and 4.3 were chosen on the basis of parameters evaluated experimentally. We can distinguish the flexibility and ability of VIM to solve the blood flow problem through compliant vessels, in the presence of the complexity of the added terms. Despite these performances, we distinguish some limitations of VIM. In these examples, going from one iteration to another can lead to more complicated integrals in the computation of the solution.
6. Conclusion
This article examined a model based on the Navier–Stokes system. This system describes the behaviour of blood flow in the arteries. Thanks to source terms and pressure expressions that take into account the properties of the elastic walls of the arteries, the model can cover certain interactions of the blood with the surfaces of the arteries. To simplify the study of the problem, a process of asymptotic reduction of the incompressible Navier–Stokes equations was applied while keeping the problem its nonlocal character. Numerically, we have proposed an approximate solution to the blood flow problem in an explicit form, using the variation iteration method. The simplicity and usability of VIM is demonstrated despite the complexity and diversity of this topic. We also used a numerical method based on Matlab's PDEPE function to discuss and evaluate the VIM. The main advantage of the VIM is that the results obtained are reasonable in two iteration steps. The VIM is perfectly capable of providing a solution despite the complexity of the pressure and/or the expression of the source terms.
Disclosure statement
No potential conflict of interest was reported by the authors.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Additional information
Funding
References
- Seadawy A, Iqbal M, Dianchen L. Applications of propagation of long-wave with dissipation and dispersion in nonlinear media via solitary wave solutions of generalized Kadomtsev–Petviashvili modified equal width dynamical equation. Computers Math Appl. 2019;78:3620–3632.
- Seadawy AR, Iqbal M, Lu D. The nonlinear diffusion reaction dynamical system with quadratic and cubic nonlinearities with analytical investigations. Int J Modern Phys B. 2020;34(9):2050085.
- Seadawy AR, Iqbal M, Lu D. Nonlinear wave solutions of the Kudryashov–Sinelshchikov dynamical equation in mixtures liquid–gas bubbles under the consideration of heat transfer and viscosity. J Taibah University Sci. 2019;13(1):1060–1072.
- Nasreen N, Seadawy AR, Lu D, et al. Construction of modulation instability analysis and optical soliton solutions of pertubed nonlinear Schrödinger dynamical equation with power law nonlinearity in non-kerr medium. Res Phys. 2019;13:102263.
- Arshad M, Seadawy AR, Lu D, et al. Dispersive solitary wave solutions of strain wave dynamical model and its stability. Commun Theor Phys. 2019;71:1155–1162.
- Chamekh M, Brik N, Abid M. Mathematical and numerical study of the concentration effect of red cells in blood. J King Saud Univ Sci. 2019;31(4):1463–1470.
- Iqbal M, Seadawy AR, Lu D. Construction of solitary wave solutions to the nonlinear modified Kortewege–de Vries dynamical equation in unmagnetized plasma via mathematical methods. Modern Phys Lett. 2018;33(32):1850183.
- Iqbal M, Seadawy AR, Lu D. Dispersive solitary wave solutions of nonlinear further modified Korteweg–de vries dynamical equation in an unmagnetized dusty plasma. Modern Phys Lett. 2018;33(37):1850217.
- Iqbal M, Seadawy AR, Lu D. Applications of nonlinear longitudinal wave equation in a magneto-electro-elastic circular rod and new solitary wave solutions. Modern Phys Lett B. 2019;33(18):1950210.
- Iqbal M, Seadawy AR, Khalil OH, et al. Propagation of long internal waves in density stratified ocean for the (2 + 1)-dimensional nonlinear Nizhnik–Novikov–Vesselov dynamical equation. Res Phys. 2020;16:102838.
- Seadawy AR, Iqbal M, Lu D. Propagation of kink and anti-kink wave solitons for the nonlinear damped modified Korteweg–de Vries equation arising in ion-acoustic wave in an unmagnetized collisional dusty plasma. Physica A. 2020;544:123560.
- Arshad M, Seadawy AR, Lu D, et al. Elliptic function solutions, modulation instability and optical solitons analysis of the paraxial wave dynamical model with kerr media. Opt Quantum Electron. 2021;53:7.
- Arshad M, Lu D, Wang J. (N+1)-dimensional fractional reduced differential transform method for fractional order partial differential equations. Commun Nonlinear Sci Numerical Simul. 2017;48:509–519.
- Barati S, Ivaz K. Variational iteration method for solving systems of linear delay differential equations. Int J Comput Math Sci. 2012;6:132–135.
- Alderremy AA, Elzaki TM, Chamekh M. New transform iterative method for solving some Klein–Gordon equations. Res Phys. 2018;10:655–659.
- Elzaki MT, Chamekh M. Solving nonlinear fractional differential equations using a new decomposition method. Universal J Appl Math Comput. 2018;6:27–35.
- Chamekh M, Elzaki TM. Explicit solution for some generalized fluids in laminar flow with slip boundary conditions. J Math Computer Sci. 2018;18:72–281.
- Chamekh M, Latrach MA, Jday F. Multi-step semi-analytical solutions for a CHIKungunya virus system. J Umm Al-Qura Univ Appl Sci. 2023.
- He JH. Variational iteration method – a kind of non-linear analytical technique: some examples. The variational iteration method: a reliable analytic tool for solving linear and nonlinear wave equations. Int J Non-Linear Mech. 1999;34:699–708.
- He JH. Some asymptotic methods for strongly nonlinear equations. Int J Modern Phys B. 2006;20(10):1141–1199.
- Wazwaz AM. The variational iteration method for solving linear and nonlinear systems of PDEs. Computers Math Appl. 2007;54:895–902.
- Wazwaz AM. The variational iteration method: a reliable analytic tool for solving linear and nonlinear wave equations. Computers Math Appl. 2007;54:926–932.
- Bertaglia G, Caleffi V, Valiani A. Modeling blood flow in viscoelastic vessels: the 1D augmented fluid-structure interaction system. Comput Methods Appl Mech Eng. 2020;360:112772.
- Krivovichev GV. Computational analysis of one-dimensional models for simulation of blood flow in vascular networks. J Comput Sci. 2022;62:101705.
- Quarteroni A. Cardiovascular mathematics. In: Proceedings of the International Congress of Mathematicians, Madrid, Spain; 2006.
- Quarteroni A, Valli A. Domain decomposition methods for partial differential equations. Oxford: International Oxford University Press; 1999.
- Quarteroni A, Tuveri M, Veneziani A. Computational vascular fluid dynamics: problems, models and methods. Comput Vis Sci. 2000;2(4):163–197.
- Quarteroni A, Veneziani A. Analysis of a geometrical multiscale model based on the coupling of ODE and PDE for blood flow simulations. Multiscale Modeling Simul. 2003;1:173–195.
- Bessonov N, Sequeira A, Simakov S, et al. Methods of blood flow modelling. Math Model Nat Phenom. 2016;11:1–25.
- Čanić S, Kim EH. Mathematical analysis of the quasilinear effects in a hyperbolic model blood flow through compliant axi-symmetric vessels. Math Methods Appl Sci. 2003;26:1161–1186.
- Formaggia L, Veneziani A. Reduced and multiscale models for the human cardiovascular system. Lecture notes VKI lecture series, researchgate.net 2003.
- Brook BS, Falle SAEG, Pedley TJ. Numerical solutions for unsteady gravity-driven flows in collapsible tubes: evolution and roll-wave instability of a steady state. J Fluid Mech. 1999;396:223–256.
- Kholodov AS. Some dynamical models of external breathing and haemodynamics accounting for their coupling and substance transport. Computer models and medicine progress. Moscow: Nauka; 2001. p. 127–163 (Russian).
- Müller LO, Parés C, Toro E. Well-balanced high-order numerical schemes for one-dimensional blood flow in vessels with varying mechanical properties. J Comput Phys. 2013;242:53–85.