On One Nonlinear Mathematical Model of Blood Circulation with the Vessel Walls Reaction within the Hereditary Theory Petro Pukacha,c, Myroslava Vovka, Yurii Mylyanb, Halyna Bilushchaka and Pavlo Pukacha a Lviv Polytechnic National University, 12 Bandera str., Lviv, 79013, Ukraine b Danylo Halytsky Lviv National Medical University, 69 Pekarska str., Lviv, 79010, Ukraine c Hetman Petro Sahaidachnyi National Army Academy, 32 Heroes of Maidan str., Lviv, 79026, Ukraine Abstract The research demonstrates sufficient conditions of the existence and uniqueness for the solution in the oscillation mathematical model of the blood flow under nonlinear dissipative forces action within the theory of hereditary tube with biofactor. These results facilitated the application of different (explicit and implicit) numerical methods in further studies of the dynamical characteristics of solutions in the considered oscillation mathematical models. Numerical integration of the movement equations by Runge-Kutta 4th order method and Geer 2nd order method in a model case within this research enabled the estimation of the influence of different physical and mechanical factors on the amplitude and frequency of the oscillation process. The use of hybrid methods for the oscillation modeling in the nonlinear isotropic elastic environment on the example of a vessel enabled the formulation of the equation of an object’s mechanical state based on energy approaches and the theory of mechanical fields in the continuous environments. Keywords 1 Mathematical Model, Nonlinear Vibrations, Biofactor, Blood circulation, Vessel. 1. Introduction Modern social development trends such as the problems of human survival and healthy lifestyle preservation are closely interconnected with the general human problems. Under such circumstances, top priority tasks are to improve the quality of human life, to devise a formula for active longevity, and to raise the individual living standards. Physical and spiritual human self-awareness should be also considered. That is why the problem of an adequate mathematical modeling of the processes in living organisms is the problem of current interest for the modern healthcare and science in general. The spectrum of the considered problems is so wide that the whole investigation review is almost impossible. Brief but sufficiently capacious review of the mathematical models for many medical- biological processes, based on the well-known models and methods of the mechanics of continuous media is presented in [1]. More deepened analysis of medical-biological and mathematical aspects is proposed, particularly, in [2-10]. Numerical modeling of the biomechanical processes in the medical practice is realized using the models of mechanics of continuous media and numerical methods of solving the corresponding partial differential equations systems. Such modeling takes into account the development and realization of the numerical methods, adapted to the specified concrete tasks, development of the numerical method algorithm and its program package, visualization of the obtained results. To study some medical processes there is necessary to solve numerically the differential equations systems [7]. Biological and medical problems involving to numerical solutions of the partial differential equations are described The 3rd International Conference on Informatics & Data-Driven Medicine (IDDM 2020), November 19–21, 2020, Växjö, Sweden EMAIL: ppukach@gmail.com(Petro Pukach); mira.i.kopych@gmail.com (Myroslava Vovk); myp.ct2019@gmail.com (Yurii Mylyan); halyna.i.bilushchak@lpnu.ua (Halyna Bilushchak); pavlopukach@mail.com (Pavlo Pukach) ORCID: 0000-0002-0359-5025 (A. 1); 0000-0002-7818-7755 (A. 2); 0000-0002-5518-0956 (A. 3); 0000-0002-1226-8050 (A. 4); 0000- 0002-0488-6828 (A. 5) ©️ 2020 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR Workshop Proceedings (CEUR-WS.org) in the papers [6-10]. Rheological relationships for the biological continuous media are developed in [11]. The mechanical model of the heart was considered in [12, 13]. Description of the simplest mathematical models of the circulatory system and heart one can find in [14-19]. Circulatory system consisting from the large and lesser circulation, possesses very serious and different functions, that is why their modeling in the normal and pathological conditions, is very important task in medicine. For today the most adequacy to the real physical circulatory systems there are dynamical models of the pulsating flows of the incompressible fluids in the elastic tubes system. 2. Hemodynamics equation in the linear mathematical models of blood circulation Hemodynamics equation differs in every case [20]. 1. The linear elastic tube. Oscillations equation of the blood flow is the next: , where is blood consumption, is area of the tube cross-section, is the tube radius, is blood averaged flow rate, is blood density, is coefficient of kinematic viscosity. 2. The linear hereditary tube. Oscillations equation is: , where , is stress, is deformation in the corresponding is layer of a multilayer cylindrical tube (vessel); is some hereditary function. 3. The linear hereditary biofactor tube. Oscillations equation is , where ,  is time delay of the reaction, , , , is radial displacement of the wall for the multilayer package as a whole. Solutions of these equations can be found as finite sum of the main oscillation and the higher harmonics, using the harmonic analysis methods. 3. Problem statement. Mathematical models of blood circulation within the nonlinear theory Study of the wave propagation process in the deformed tubes with the liquid leaking through the tube is widespread applicated [21, 22]. These problems are actual, in particular, in case of blood circulation modeling in alive organisms. Problems of blood flows and oscillations propagation in large blood vessels are very important to understand the functioning, regulation and control the cardiovascular system [23]. As follows, diagnostics, surgery and prosthetics are bound up the hemodynamics [20, 24]. In the mathematical modeling of blood flow there is considered pulsed systaltic blood flow in the multilayer elastic or viscous elastic tube with the variable cross-section. More complicated mathematical models of blood circulation in the tubes that possess reaction on the external action (biofactor). This type models describe blood circulation in arteries and veins. To study complicated impulses, specific to the circulatory system, it is often necessary to consider the nonlinear models instead the linear ones. It is impossible to find the exact solution in this case. Background of the nonlinear mathematical models study is numerical (computer) modeling. Application of this approach also can’t be universal study method due to the other problems, for example, procedure convergence, numerical method’s stability, accuracy of computation. That’s why it is reasonable to develop hybrid methods to study the nonlinear oscillations mathematical models combining both the qualitative and numerical approaches [25, 26]. It is realized thorough qualitative description of the solution’s characteristics respectively to the problem. Usually, it is based on Galerkin method or on its different modifications. After that the numerical methods are applicated to find the approximate solutions. Herewith choice of the numerical method is of no principle from the theoretical aspect, and can be determined by effectiveness of the numerical realization only. Using all mentioned above, let’s consider the oscillations equation of blood flow in the nonlinear isotropic environment. It would be studied the problem of the blood flow oscillations equation within the mathematical model of the linear hereditary tube with biofactors and nonlinear external dissipative forces in the form , (1) where , is the linear blood density in the finite tube with fixed length , is uniformly distributed by tube length force, causing the initial oscillations, , are the coefficients of the external and internal dissipations of the environment respectively, . The nonlinear oscillations of the tube with the constant radius of the cross-section R in the case of the initial amplitude replacement and the initial zero rate would be considered. The case of rigidly fixed by length tube also would be considered. It is reasonable to study the equation (1) in the rectangle mixed problem with the initial conditions , (2) and the boundary conditions . (3) For the case the equation (1) is considered in the previous chapter within the linear elastic tube with the reaction (biofactor) model. In this paper the linear model is involved by the nonlinear factors. 4. Model case of the elastic tube with the reaction and computer modeling results For the mathematical modeling of the free nonlinear small transversal oscillations there is used long fixed on the endpoints tube under the force action on the unit of length . Free oscillations of the blood flow are described by the problem (1), (2), (3). Under these conditions the space discretization of the equation (1) is realized. Let be the quantity of the discretization components, and be difference interval over the coordinate . Solving of the problem (1), (2), (3) results in the numerical integration on the time interval of some difference equations system under some initial conditions. Computer modeling of the transient processes was realized on the model example of the analysis of the small transversal oscillations of the elastic thin tube m with the finite length. Tube length is 1 m, blood density , wall thickness is m. Tube is under the initial pertubation of the force applied to its center in the dilation direction (that means perpendicular to the line length). System parameters are the next: , x  0,1 m,   1 . Integration of the equations in the mechanical mode is realized via the explicit Runge-Kutta fourth order method and implicit second order Geer method. Numerical results almost coinside. Integration step of the explicit Runge-Kutta method is 1105 s, implicit Geer method is 1104 s. The nonlinear algebraic equations system on every step by the variable is solved via the simple iteration method. Three modes of the object were studied. The first mode presents the oscillations of the sufficiently small flows in the tube (the linear internal dissipation of the mechanical energy is present only, ). The second mode presents half-filled tube (the linear internal and linear external dissipations of the mechanical energy are present, , ). The third mode presents the oscillations of the filled tube (the linear internal and the nonlinear external dissipations of the mechanical energy are present, , ). Obviously, that these assumptions are adapted, but even in this case with the sufficient adequacy extent there are described the real physical processes in the object. Thus, there realized three experiments taking into consideration mentioned modes. To confirm the validity of the system model also were carried out two additional experiments consisting in study of the transitional processes with different tube thickness. The first experiment examined m (Fig. 1, Fig. 2), the second experiment examined m (Fig. 3). Figure 1: Transient motions of the tube central component ,( m): 1 is the first experiment, 2 is the second experiment, 3 is the third experiment Figure 2: Transient motions of the tube central component ( m): 1 is the second experiment, 2 is the third experiment at time span Figure 3: Transient motions of the tube central component ( m): 1 is the first experiment, 2 is the second experiment, 3 is the third experiment Analyzing the family of curves, one can see the essential influence of the internal dissipative processes on the tube oscillations. Partly filled tube can be treated as an object with the linear external dissipation in contradistinction to the fully filled tube. Turbulent processes of the blood circulation cause the nonlinear influence on the vessel walls. This influence depends on the vessel wall thickness. This fact is clearly fixed on the Fig. 3. Reducing of the tube thickness causes the increasing of the eigenoscillation frequency and contrariwise. This fact absolutely corresponds to the classical elasticity theory. The oscillations damping in the tube with the less thickness are more intensive being dependent on the internal processes in the vessel body. The nonlinear characteristics of the blood flow are manifested stronger in the vessels with the thin walls, that is well-understood from the physical point of view. 5. Conclusions The elaboration of the mathematical models of the physiological processes in the able-bodied organism, and also medical problems that follow in the sick mode of the patient, can be considered as mathematical modeling domain that is intensively developed. Hereby the qualitative study of the mathematical model and the next numerical modeling often are effective and available instrument of the biological and medical problems investigation. To confirm the problem correctness in the nonlinear mathematical model of the blood circulation in the paper are used the fundamental methods of the nonlinear boundary problems general theory. Basing on the results of the numerical modeling there is proved the sufficient adequacy of the obtained model to the real prototype. It is shown that the nonlinear medium promotes to the quicker oscillations damping and causes the inharmonic processes in the system. The well-known fact, that increasing of the vessel walls thickness causes the reducing eigenoscillations frequency of the system and contrariwise also is reaffirmed. 6. References 1. G.I. Marchuk, Mathematical models in immunology, Nauka, Moscow, 1985 [in Russian]. 2. E. Haier, G. Wiener, Solution of ordinary differential equations. Rigid and differential- hyperbolic problems, Mir, Moscow, 1999 [in Russian]. 3. I.M. Makarov, Informatics and medicine, Nauka, Moscow, 1997 [in Russian]. 4. O.M. Belotserkovsky, A.S. Kholodov, Computer models and medical progress, Nauka, Moscow, 2001 [in Russian]. 5. O.M. Belotserkovsky, Computer and brain. New technologies, Nauka, Moscow, 2005 [in Russian]. 6. A.A. Samarsky, Theory of difference schemes, Nauka, Moscow, 1977 [in Russian]. 7. O.M. Belotserkovsky, Numerical modeling in continuum mechanics, Fizmatlit, Moscow, 1994 [in Russian]. 8. K.M. Magomedov, A.S. Kholodov, Grid-characteristic methods, Nauka, Moscow, 1988 [in Russian]. 9. S.A. Regirer, Lectures on biological mechanics, Nauka, Moscow, 1980 [in Russian]. 10. V.I. Kondaurov, A.V. Nikitin, Final deformations of viscoelastic muscle tissues, Applied Mathematics and Mechanics 51(3) (1987) 443–452 [in Russian]. 11. M.B. Akhundov, Forced vibrations of elastic and viscoelastic rods in contact with a medium with the property of self-regulation, Herald of Baku State University 2 (2011) 63–72 [in Russian]. 12. K.D. Gosta, P.J. Hunter, J.M. Pogers, G.M. Gussione, L.K. Waldman, A.D. McCulloch, A three dimensional limite elements method for large elastic deformations of ventricular myocardium. Part I, J. Biomech. Eng. 118(4) (1996) 452–463. doi.org/10.1115/1.2796032 13. S.C. Panda, R. Natarajon, Finite-element method of stress analysis in the human left ventricular layered wull structure, Med. Biol. Eng. Comp. 15 (1977) 67–71. doi.org/10.1007/BF02441577 14. P.I. Begun, P.N. Afonin, Modeling in biomechanics, Higher School, Moscow (2004) [in Russian]. 15. M.V. Abakumov, I.V. Ashmetov, N.B. Eshkova, V.B. Koshelev, S.I. Mukhin, N.V. Sosnin, V.F. Tishkin, A.P. Favorsky, A.B. Khrumenko, Methods of mathematical modeling of the cardiovascular system, Mathematical modeling 12 (2) (2000) 106–117 [in Russian]. 16. I.V. Ashmetov, A.Ya. Bunicheva, S.I. Mukhin, T.V. Sokolova, N.V. Sosnin, A.P. Favorsky, Mathematical modeling of hemodynamics in the brain and in the great circle of blood circulation, in: Computer and brain. New technologies, Nauka, Moscow, 2005 [in Russian]. 17. A.S. Kholodov, Some dynamic models of external respiration and blood circulation taking into account their connectivity and transport, in: Computer models and the progress of medicine, Nauka, Moscow, 2001 [in Russian]. 18. A.V. Evdokimov, A.S. Kholodov, Quasi-stationary spatially distributed model of closed blood circulation of the human body, in: Computer models and the progress of medicine, Nauka, Moscow 2001 [in Russian]. 19. T.J. Peadley, R.C. Schroter, M.F. Sudllow, Energy loses and pressure drop in models of human airways, Respir. Physiol. 9 (1970) 371–386. doi.org/10.1016/0034-5687(70)90093-9 20. V.G. Nasibov, Waves in a multilayer tube of variable cross section with flowing viscous liquid, Proceedings of the IMM of the Azerbaijan Academy of Sciences VI (XIV) (1997) 245-258 [in Russian]. 21. P.Ya. Pukach, I.V. Kuzio, Z.M. Nytrebych, V.S. Ilkiv, Analytical methods for determining the effect of the dynamic process on the nonlinear flexural vibrations and the strength of compressed shaft, Naukovyi Visnyk Natsіonalnoho Hіrnychoho Unіversytetu 5 (2017) 69-76. URL: http://nbuv.gov.ua/UJRN/Nvngu_2017_5_13 22. P.Ya. Pukach, I.V. Kuzio, Z.M. Nytrebych, V.S. Ilkiv, Asymptotic method for investigating resonant regimes of nonlinear bending vibrations of elastic shaft, Naukovyi Visnyk Natsіonalnoho Hіrnychoho Unіversytetu 1 (2018) 68-73. doi.org/10.29202/nvngu/2018-1/9 23. Z. Nytrebych, V. Ilkiv, P. Pukach, O. Malanchuk, I. Kohut, A. Senyk, Analytical method to study a mathematical model of wave processes under twopoint time conditions, Eastern-European Journal of Enterprise Technologies 1(7(97)) (2019) 74-83. doi.org/10.15587/1729- 4061.2019.155148 24. V.G. Nasibov, Pulsating flow of a viscous liquid in a multilayer elastic tube of variable cross section with reaction, in: Collection of scientific works on mechanics, Az ISU, Baku, 1994, pp. 134-139 [in Russian]. 25. S.P. Lavrenyuk, P.Ya. Pukach, Mixed problem for a nonlinear hyperbolic equation in a domain unbounded with respect to space variables, Ukrainian Mathematical Journal 59(11) (2007) 1708- 1718. doi.org/10.1007/s11253-008-0020-0 26. P.Y. Pukach, Qualitative Methods for the Investigation of a Mathematical Model of Nonlinear Vibrations of a Conveyer Belt, Journal of Mathematical Sciences 198(1) (2014) 31-38. doi.org/10.1007/s10958-014-1770-x