Next Article in Journal
Recent Updates on the Use of Agro-Food Waste for Biogas Production
Previous Article in Journal
Probability Analysis of Hypertension-Related Symptoms Based on XGBoost and Clustering Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Study of the Blood Flow in a Deformable Human Aorta

1
Department of Radiological Sciences and Medical Imaging, College of Applied Medical Sciences, Majmaah University, Majmaah 11952, Saudi Arabia
2
Laboratory of Electronics and Microelectronics, Faculty of Science of Monastir, University of Monastir, Environment Boulevard, Monastir 5019, Tunisia
3
Department of Physics, College of Science, Majmaah University, Zulfi 11932, Saudi Arabia
4
Department of Computer Science and Information, College of Science, Majmaah University, Zulfi 11932, Saudi Arabia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(6), 1216; https://doi.org/10.3390/app9061216
Submission received: 17 February 2019 / Revised: 13 March 2019 / Accepted: 19 March 2019 / Published: 22 March 2019

Abstract

:

Featured Application

This article developed a computer simulation based on the finite element method to investigate the interaction between blood flow and arterial wall deformation. The numerical results show that the blood flow has an unsteady and complex behavior. The velocity field, pressure, and displacement distributions were investigated during a cardiac cycle. The method described in this study will be useful for investigating cardiovascular diseases and performing hemodynamics analyses.

Abstract

In this work, we present a numerical investigation of blood flow in a portion of the human vascular system. More precisely, the present work analyzed the blood flow in the upper portion of the aorta. The aorta and its ramified blood vessels are surrounded by the cardiac muscle. The blood flow generates pressure on the internal surfaces of the artery and its ramifications, thereby causing deformation of the cardiac muscle. The numerical analysis used the Navier–Stokes equations as the governing equations of blood flow for the calculation of the velocity field and pressure distribution in the blood. The neo-Hookean hyperelastic model was used for the description of the behavior of the vessel walls. The velocity and pressure distributions were analyzed. The deformation of the vessel was also investigated. The numerical results could be used to better understand and predict the factors that trigger cardiovascular diseases and distortions of the aorta and as a diagnostic tool in clinical applications.

1. Introduction

The cardiovascular system consists of the heart and vessels (arteries and veins). It has the function to distribute to the organs oxygen and nutrients essential for their life through the blood. The heart and its various organs can be affected by different pathologies, called cardiovascular diseases, namely, stenosis, aneurysm, occlusion, atherosclerosis, etc. According to recent statistics, cardiovascular diseases are the third leading cause of death in the western world [1]. The understanding of the vascular function requires knowledge of the behavior and the material properties of blood vessels [2,3]. In addition, the deep characterization of blood flow is crucial for the understanding of the behavior of the whole human system and for the improvement of human health [4,5]. Even though a few experimental works are available on the histological structure of arteries, they are still limited by the difficulty of making rapid analyses on the heart [6,7].
Numerical simulation is a widely used tool for understanding many phenomena and especially, it was confirmed to be useful in studying the blood flow in the cardiovascular system, to diagnose cardiovascular diseases [8,9,10,11,12], evaluate the hemodynamics of healthy and diseased blood vessels [13,14,15] and design and optimize vascular medical devices [16,17]. Over the last decades, various computational techniques have been used extensively by a large number of researchers seeking to simulate and analyze blood flow and to understand the relationship between vascular diseases and hemodynamics [18,19]. In addition, numerical simulation and computer models provide a fast, simple, safe way to study phenomena involved in vascular diseases.
Recently, there have been numerous studies of blood flow related to the aorta that use computer simulation to better investigate blood flow variables. Panta et al. [20] analyzed, through their numerical simulation, blood velocity as a function of the shape of human arteries, namely straight, bent, T-shaped, and the main arterial network models. Wood et al. [21] combined magnetic resonance imaging (MRI) and computational fluid dynamic (CFD) simulation of flow to analyze the human descending aorta. Their objective was to investigate the potential of combining CFD simulations with in vivo MRI scans. Tokuda et al. [22] simulated blood flow numerically using the finite element method to cognize the mechanism of stroke during cardio-pulmonary bypass (CPB). The location of turbulence, wall pressure and flow distribution within the aortic arch were investigated. They succeeded in simulating and visualizing the blood flow during CPB using CFD techniques. Torii et al. [23] developed a computer modeling technique based on the deforming-spatial-domain/stabilized space–time method to study fluid–structure interactions between flowing blood and arterial walls. They demonstrated that the wall shear stress is a very important factor in cardiovascular diseases, and their results proved that numerical modeling is an operative tool in hemodynamic studies of cardiovascular diseases.
Bazilevs et al. [24] presented a computational framework for the simulation of patient-specific cerebral aneurysm configurations and also described a new approach for the computation of blood vessel tissue prestress. They studied some variables, namely, hemodynamics, wall shear stress, and wall tension, in order to check the relevance of fluid–structure interactions when adapting the rigidity of the arterial wall and the relevance of a flexible wall. Crosetto et al. [25] performed a numerical algorithm based on finite elements to evaluate the mean velocity and temporal evolution of pressure and the distribution of wall shear stress on diverse sections of the aorta. More specifically, they reported the results of 3D simulations using a physiological geometry over several heartbeats. They changed the inlet condition of flux by pressure which caused a backward diastolic flow through the section and thereafter led to the closure of the aortic valve and the circulation in the coronary arteries. However, these observations were only predicted numerically. Garje et al. [26] studied numerically the changes in blood flow vessels and the growth of a blockage, which is an occlusion of the coronary arteries typical of coronary artery disease (CAD), due to an accumulation of cholesterol. Their work was focused on the size and position of the blockage. Priyadharshini and Ponalagusamy [27] studied numerically the effects of magneto-hydrodynamics on blood flow parameters in a stenotic artery under an applied external magnetic field and body acceleration.
In this work, we intended to simulate the interaction of blood flow with the artery wall. We simulated the pressure along the surface of the artery, the displacements and stresses of the walls of part of a child’s aorta and some of its ramifications embedded in the heart muscle when subjected to the pressure caused by heart beating.
In Section 2, the geometry of the aorta and its ramifications are presented. Then, the Navier–Stokes equations and the equation governing the displacement of the blood vessels are introduced. The initial and boundary conditions, as well as the numerical method, are also presented. The results are discussed in Section 3.

2. Methods

2.1. Simulated Geometry

Figure 1 shows the simulated geometry. The domain was divided into three parts: blood, aorta and its ramifications (branching vessels), and cardiac muscle. The cardiac muscle is very important because it protects the artery from the deformation due to the pressure. The inlet and outlet orifices are also shown in this figure.

2.2. Governing Equations

In this study, we aimed to compute the flow of blood in the aorta and to investigate the deformation of the cardiac muscle due to blood pressure and the behavior of the vessel walls. The numerical study consisted of two parts.
Part 1:
Fluid dynamics analysis which required to solve the Navier–Stokes equations. It included the calculation of the velocity field and the distribution of blood pressure (variable over time and in space).
Part 2:
Mechanical analysis of the deformation of the cardiac tissue and artery.
Navier–Stokes equations: The blood flow fluid dynamics was described using the time-dependent Navier–Stokes equations and the equations for the conservation of mass. The fluid is assumed to be incompressible and homogeneous, with a Newtonian behavior. The blood was assumed to be Newtonian in this study, as its viscosity is almost constant in arteries with relatively large diameters (5 mm); therefore, non-Newtonian effects were neglected [28].
The Navier–Stokes equations may be written as follows:
. u = 0 ,
ρ ( u t + u . u ) = σ  
where ρ is the blood density, u is the blood velocity, and σ is the stress tensor [29]. The relationship describing the stress tensor and the strain rate tensor is defined as:
σ = p I + μ ( u + ( u ) T )
where p is the blood pressure, μ is the dynamic viscosity of the blood, I is the identity tensor, and T denotes the transpose tensor.
The fluid properties such as the density and the dynamic viscosity were assumed to be constant. Their values were, respectively ρ b l o o d = 1060 kg/m3 and μ b l o o d = 0.005 Pa.s.
Structural mechanics: The motion of the vessel wall can be mathematically described by the following equation [30]:
ρ w 2 u s o l i d t 2 σ s = ρ w f s
where ρ w is the wall density, f s represents the external body forces acting on the solid, and σ s is the Cauchy stress tensor. The wall can be the artery or the cardiac muscle. The value of the wall density is ρ a r t e r y = 960 kg/m3 for the artery and ρ m u s c l e = 1200 kg/m3 for the cardiac muscle. It should be noted that u is a velocity (unit m/s); however, u s o l i d is a displacement (unit meter).
The artery is assumed to be subjected to large strains and undergo large deformations. Its behavior is non-linear. More precisely, it behaves as a hyperelastic neo-Hookean material. The constitutive law of this material involves the material shear modulus μ s and the bulk modulus κ s . The values of these parameters are provided in Table 1. In the present study, the material is assumed to be isothermal, and no effect of the temperature is considered.
For a three-dimensional, compressible, neo-Hookean material such as the artery, the strain energy density is given by [24]:
W = 1 2 μ s ( 1 J 2 / 3 t r ( C ) 3 ) + 1 2 κ s [ 1 2 ( J 2 1 ) ln ( J ) ]
where J = det F , F is the deformation gradient tensor, and C = F T F is the Cauchy–Green deformation gradient tensor. From the above constitutive law, the expressions of various tensors can be deduced by means of more or less arduous calculations.

2.3. Initial and Boundary Conditions

As mentioned above, the model uses three domains, namely, the blood, the artery and the cardiac muscle. The numerical simulation considers the solution of the 3D Navier–Stokes equations and the fluid–structure interactions. The initial and boundary conditions required to solve the equations of the model are given hereafter.
The human heart beat is cyclically time-varying. Thus, the velocity at the inlet of the aorta is not set to be constant. It has a periodic profile as a function of time. During a cycle, the velocity varies between a minimal and a maximal value. The duration of each period is 0.8 s. We introduced the following function [31]:
f ( t ) = { 0.8   sin ( 2.857 π t )                                                                 0 t 0.35 0.6   sin ( 2.22 π ( t 0.35 ) )                             0.35 t 0.8  
The inlet velocity for the artery in this simulation was fixed to 0.3 [m/s]*f(t). The temporal evolution of the inlet velocity is represented in Figure 2. The outlet condition pressure was considered to be 125 [mmHg]*f(t). These conditions are similar to those adopted by Choudhari and Panse [31]. This function has two peaks which have different amplitudes and different periods.
On the aortic walls, a non-slip boundary condition is applied, and the wall is assumed to be impermeable since the fluid is supposed not to penetrate through the walls as it flows along the vessel.

2.4. Numerical Method

In this study, the finite element method was used to solve the set equations describing the problem [32]. The simulated geometry was meshed in tetrahedral discretization as shown in Figure 3a. Figure 3b depicts the radial discretization of the blood aorta. The central part stands for the blood. The artery and its thickness appear clearly in this figure. However, the muscle was considered as having no thickness, and only a boundary condition was applied to represent this surface. Indeed, the muscle was considered a free surface, with no external force applied on it.
The model used to solve the governing equations is two-way coupled. Actually, any deformation of blood vessel walls has an effect on the fluid flow and vice versa [33,34].
Mesh sizes are summarized in Table 2. The total number of elements were 195,050, and the total number of nodes was 33,506.

3. Results

The results concern mainly the temporal evolution of the velocity field, the pressure, and the wall displacement.
Figure 4 shows the blood speed distribution inside the vessel at different times. To discuss this figure, one should keep in mind the profile of the inlet velocity depicted in Figure 2. At t = 0.1 s, the fluid was mainly directed to the first inlet. The section between the first and second outlets started to move. However, the remaining part of the aorta was almost at rest. At t = 0.2 s, the fluid also reached the second and the third exits. Progressively, the blood flowed into the vessel accelerating or decelerating. The un-steady flow remained laminar.
From the results of Figure 4, one can also compute the flow rate in each output as well as the stroke volume.
Figure 5 exhibits the temporal evolution of the maximal velocity. The maximum value reached varied from 0.1 to 0.32 m/s depending on time. The relative curve is coarsely similar to that of the input velocity. It is constituted by two peaks that have different amplitudes and different frequencies. A least-square fit allowed to obtain the amplitude of these peaks. We obtained 0.35 m/s for the first peak and 0.26 m/s for the second one. The discrepancy between the computation results and the theoretical results obtained by the least-square fit was in a tolerable limit. The mean uncertainty was about 10%.
It was also interesting to analyze the pressure distribution versus time. Figure 6 exhibits the pressure distribution at t = 0.1 s and t = 0.7 s. It is clear that the pressure has a complex distribution. Furthermore, this distribution depends strongly on time, since the input conditions are unsteady.
The values of the wall shear stress applied to the blood vessel walls during a cardiac cycle vary strongly in time and position. This has an important effect on cardiovascular diseases.
Figure 7 shows the temporal evolution of the maximal pressure. This figure confirms the large variation in pressure during a cardiac cycle. Evidently, this has an important effect on the wall shear stress and consequently on cardiovascular diseases. Like the maximal velocity, the temporal evolution of the maximal pressure consisted of two pulses. The amplitudes of the pulses were 12.65 kPa and 9.62 kPa. The mean uncertainty between the computational results and the data obtained from the least-square fit was lesser than 8%.
Figure 8 depicts the displacement of the blood vessel wall at t = 0.1 s and t = 0.7 s. It is clear that a large deformation of the vessel occurred mainly in the vicinity of the embranchments. However, all the vessels may be subjected to large deformation.
The temporal evolution of maximal displacement during a cardiac cycle was investigated. Figure 9 presents the related results. Similarly to the maximal velocity and the maximal pressure, the temporal evolution of the maximal pressure consisted of two pulses. Using the least-square fit, we found that the amplitudes were 8.33 μm and 6.26 μm. The mean uncertainty due to the fit was about 3%.
It appears clearly that the maximal displacement depended on time and could reach 8 μm. These values are in accordance with the important shear stress on the vessel. This has an impact on cardiovascular diseases such as stenosis.
Our results have the same order of magnitude as previous works [31,35,36,37]. For instance, Choudhari and Panse obtained a radial deformation ranging from 1.5 to 2.2 µm and a maximal velocity about 0.6 m/s.
In this study, the pressure, the velocity, and the deformation were predicted with respect to a given blood flow velocity. It could be interesting, in future work, to analyze the inverse problem. A numerical method such as the iterative gradient search method is expected to be useful to perform this analysis [38,39].

4. Conclusions

In this work, we developed a computer simulation based on the finite element method to investigate the interaction between blood flow and arterial wall deformation. This simulation tool was applied to the upper part of a child’s aorta. The numerical results showed that the blood flow had an unsteady and complex behavior. The velocity field, pressure, and displacement distributions were investigated during a cardiac cycle. Further, the present study reveals that modeling the blood flow is helpful to the investigation of cardiovascular diseases and it can be an effective tool in hemodynamics studies of cardiovascular diseases.
Experimental validation on real human aorta and a silicon rubber model would be very useful. We intend to perform this validation in a future work.

Author Contributions

M.S. proposed the idea, contributed to data acquisition, performed the simulation. H.B. contributed to data analysis and result discussion. A.B. performed the algorithm construction.

Acknowledgments

The authors thank the Deanship of Scientific Research and the Deanship of Community Service at Majmaah University, Saudi Arabia, for supporting this work. The authors thank Shatha Mohammed Al-Mutairy, Danah Bander Al Mutairi, Reef Saleh Al-Aybani, Sarah Ibrahim Altwaijry for participating in this work.

Conflicts of Interest

Competing Interests: The authors declare no competing financial interests.

References

  1. Kroon, M.; Holzapfel, G.A. Modeling of saccular aneurysm growth in a human middle cerebral artery. J. Biomech. Eng. 2008, 130, 051012. [Google Scholar] [CrossRef] [PubMed]
  2. Kamalanand, K.; Srinivasan, S. Modelling and analysis of normal and atherosclerotic blood vessel mechanics using 3d finite element models. Intact J. Soft Comput. 2011, 2, 261–264. [Google Scholar] [CrossRef]
  3. Vito, R.P.; Dixon, S.A. Blood Vessel Constitutive models-1995–2002. Annu. Rev. Biomed. Eng. 2003, 5, 413–439. [Google Scholar] [CrossRef]
  4. Taylor, C.A.; Figueroa, C.A. Patient-specific model of cardiovascular mechanics. Annu. Rev. Biomed. Eng. 2009, 11, 109–134. [Google Scholar] [CrossRef] [PubMed]
  5. Gerbeau, J.-F.; Vidrascu, M. A quasi-newton algorithm based on a reduced model for fluid-structure interaction problems in blood flows. Math. Model. Numer. Anal. 2003, 37, 631–647. [Google Scholar] [CrossRef]
  6. Zulliger, M.A.; Fridez, P.; Hayashi, K.; Stergiopulos, N. A strain energy function for arteries accounting for wall composition and structure. J. Biomech. 2004, 37, 989–1000. [Google Scholar] [CrossRef] [PubMed]
  7. Spencer, A.J.M. Continuum Mechanics; Longman Scientific & Technical: Essex, UK, 1980. [Google Scholar]
  8. Ras, R.T.; Streppel, M.T.; Draijera, R.; Zock, P. Flow-mediated dilation and cardiovascular risk prediction: A systematic review with meta-analysis. Int. J. Cardiol. 2013, 20, 168, 344–351. [Google Scholar] [CrossRef] [PubMed]
  9. Bastien, M.; Poirier, P.; Lemieux, I.; Després, J.-P. Overview of epidemiology and contribution of obesity to cardiovascular disease. Prog. Cardiovasc. Dis. 2014, 56, 369–381. [Google Scholar] [CrossRef]
  10. Regitz-Zagrosek, V.; Oertelt-Prigione, S.; Prescott, E.; Franconi, F.; Gerdts, E.; Foryst-Ludwig, A.; Maas, A.H.; Kautzky-Willer, A.; Knappe-Wegner, D.; Kintscher, U.; et al. Gender in cardiovascular diseases: Impact on clinical manifestations, management, and outcomes. Eur. Heart J. 2016, 37, 24–34. [Google Scholar] [CrossRef]
  11. Babu, G.C.; Shantharajah, S.P. Optimal body mass index cutoff point for cardiovascular disease and high blood pressure. Neural Comput. Appl. 2018. [Google Scholar] [CrossRef]
  12. Narayan, S.; Sathiyamoorthy, E. A novel recommender system based on FFT with machine learning for predicting and identifying heart diseases. Neural Comput. Appl. 2018. [Google Scholar] [CrossRef]
  13. Leo, M.C.; Jayashree, S.; Darren, R.G.; Marsel, M.M.; Todd, B.P. Hemodynamic response changes in cerebrovascular disease: Implications for functional MR imaging. AJNR Am. J. Neuroradiol. 2002, 23, 1222–1228. [Google Scholar]
  14. Pang, C.L.; Alcock, R.; Pilkington, N.; Reis, T.; Roobottom, C. Determining the haemodynamic significance of arterial stenosis: The relationship between CT angiography, computational fluid dynamics, and non-invasive fractional flow reserve. Clin. Radiol. 2016, 71, 750–757. [Google Scholar] [CrossRef] [PubMed]
  15. Dirk, J.D.; Akos, K.; Daphne, M.; John, M.C.J. Regulation of coronary blood flow in health and ischemic heart disease. Prog. Cardiovasc. Dis. 2015, 57, 409–422. [Google Scholar]
  16. Mark, S.S.; Francis, D.P.; Joseph, G.R.; Leslie, W.M.; Benjamin, S.; Stuart, D.R.; Randall, C.S.; Leway, C.; Andrew, J.B.; Suzanne, C.; et al. Clinical management of continuous-flow left ventricular assist devices in advanced heart failure. J. Heart Lung Transplant. 2010, 29, S1–S39. [Google Scholar]
  17. Xenosa, M.; Girdhara, G.; Alemua, Y.; Jesty, J.; Slepian, M.; Einava, S.; Bluestein, D. Device Thrombogenicity Emulator (DTE)-Design optimization methodology for cardiovascular devices: A Study in two bileaflet MHV designs. J. Biomech. 2010, 43, 2400–2409. [Google Scholar] [CrossRef] [PubMed]
  18. Purkayastha, S.; Fadar, O.; Mehregan, A.; Salat, D.H.; Moscufo, N.; Meier, D.S.; Guttmann, C.R.G.; Fisher, N.D.L.; Lipsitz, L.A.; Sorond, F.A. Impaired cerebrovascular hemodynamics are associated withcerebral white matter damage. J. Cereb. Blood Flow Metab. 2014, 34, 228–234. [Google Scholar] [CrossRef] [PubMed]
  19. Lewis, G.D.; Bossone, E.; Naeije, R.; Grünig, E.; Saggar, R.; Lancellotti, P.; Ghio, S.; Varga, J.; Rajagopalan, S.; Oudiz, R.; et al. Pulmonary vascular hemodynamic response to exercise in cardiopulmonary diseases. Circulation 2013, 128, 1470–1479. [Google Scholar] [CrossRef] [PubMed]
  20. Panta, Y.M.; Harvey, M. Technical brief: Computational fluid dynamic (CFD) analysis of blood flow through human arteries. J. Comput. Simul. Model. 2012, 2, 27–29. [Google Scholar]
  21. Wood, N.B.; Weston, S.J.; Kilner, P.J.; Gosman, A.D.; Firmin, D.N. Combined MR imaging and CFD simulation of flow in the human descending aorta. J. Magn. Reson. Imaging 2001, 13, 699–713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Tokuda, Y.; Song, M.H.; Ueda, Y.; Usui, A.; Akita, T.; Yoneyama, S.; Maruyama, S. Three-dimensional numerical simulation of blood flow in the aortic arch during cardiopulmonary bypass. Eur. J. Cardiothorac. Surg. 2008, 33, 164–167. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Torii, R.; Oshima, M.; Kobayashi, T.; Takagi, K.; Tezduyard, T.E. Computer modeling of cardiovascular fluid–structure interactions with the deforming-spatial-domain/stabilized space–time formulation. Comput. Methods Appl. Mech. Eng. 2006, 195, 1885–1895. [Google Scholar] [CrossRef]
  24. Bazilevs, Y.; Hsu, M.C.; Zhang, Y.; Wang, W.; Kvamsdal, T.; Hentschel, S.; Isaksen, J.G. Computational vascular fluid–Structure interaction: Methodology and application to cerebral aneurysms. Biomech. Model. Mechanobiol. 2010, 9, 481–498. [Google Scholar] [CrossRef] [PubMed]
  25. Crosettoa, P.; Reymond, P.; Deparis, S.; Kontaxakis, D.; Stergiopulos, N.; Quarteroni, A. Fluid–Structure interaction simulation of aortic blood flow. Comput. Fluids 2011, 43, 46–57. [Google Scholar] [CrossRef]
  26. Garje, A.; Adhav, Y.G.; Bodas, D. Design and simulation of blocked blood vessel for early detection of heart diseases. In Proceedings of the 2015 2nd International Symposium on Physics and Technology of Sensors Pune, India, 7–10 March 2015; pp. 204–208. [Google Scholar]
  27. Priyadharshini, S.; Ponalagusamy, R. Mathematical modelling for pulsatile flow of Casson fluid along with magnetic nanoparticles in a stenosed artery under external magnetic field and body acceleration. Neural Comput. Appl. 2017. [Google Scholar] [CrossRef]
  28. McDonald, D.A. Blood Flow in Arteries, 2nd ed.; Edward Arnold: London, UK, 1974. [Google Scholar]
  29. Chuchard, P.; Puapansawat, T.; Siriapisith, T.; Wu, Y.H.; Wiwatanapataphee, B. Numerical simulation of blood flow through the system of coronary arteries with diseased left anterior descending. Int. J. Math. Comput. Simul. 2011, 5, 334–341. [Google Scholar]
  30. Torii, R.; Oshima, M.; Kobayashi, T.; Takagi, K.; Tezduyard, T.E. Influence of wall elasticity in patient-specific hemodynamic simulations. Comput. Fluids 2007, 36, 160–168. [Google Scholar] [CrossRef]
  31. Choudhari, P.; Panse, M.S. Finite Element Modeling and simulation of arteries in the human arm to study the aortic pulse wave propagation. Procedia Comput. Sci. 2016, 93, 721–727. [Google Scholar] [CrossRef]
  32. Rao, S.S.; Chen, Z. The Finite Element Methods in Engineering; Elsevier Science & Technology Books: Miami, FL, USA, 2004. [Google Scholar]
  33. Urquiza, S.A.; Blanco, P.J.; Venere, M.J.; Feijoo, R.A. Multidimensional modelling for the carotid artery blood flow. Comput. Methods Appl. Mech. Eng. 2006, 195, 4002–4017. [Google Scholar] [CrossRef]
  34. Wang, Y.; Quaini, A.; Čanić, S. A higher-order discontinuous Galerkin/arbitrary Lagrangian Eulerian partitioned approach to solving fluid–Structure interaction problems with incompressible, viscous fluids and elastic structures. J. Sci. Comput. 2018, 76, 481–520. [Google Scholar] [CrossRef]
  35. Kim, K.; Kim, J.U.; Beak, H.M.; Kim, S.K. Modeling of the artery tree in the human upper extremity and numerical simulation of blood flow in the artery tree. Trans. Korean Soc. Mech. Eng. B. 2016, 40, 221–226. [Google Scholar] [CrossRef]
  36. Wu, Q.; Liu, X.; Pan, L.; Tao, W.; Qian, C. Simulation analysis of blood flow in arteries of the human arm. Biomed. Eng. 2017, 29, 1750031. [Google Scholar] [CrossRef]
  37. Marie, W.; Valerie, L.; Emilie, M. Inlet boundary conditions for blood flow simulations in truncated arterial networks. J. Biomech. 2011, 44, 897–903. [Google Scholar]
  38. Jinqiang, N.; Steven, Y.L. Inverse identification of Johnson-Cook material constants based on modified chip formation model and iterative gradient search using temperature and force measurements. Int. J. Adv. Manuf. Technol. 2019. [Google Scholar] [CrossRef]
  39. Jinqiang, N.; Vinh, N.; Yong, H.; Karl, T.H.; Steven, Y.L. Inverse determination of Johnson–Cook model constants of ultra-fine-grained titanium based on chip formation model and iterative gradient search. Int. J. Adv. Manuf. Technol. 2018, 99, 1131–1140. [Google Scholar]
Figure 1. Schematic showing a view of the aorta, its ramifications (branching vessels), and the cardiac muscle.
Figure 1. Schematic showing a view of the aorta, its ramifications (branching vessels), and the cardiac muscle.
Applsci 09 01216 g001
Figure 2. Flow velocity waveform in the inlet used in the simulation for the flow analysis.
Figure 2. Flow velocity waveform in the inlet used in the simulation for the flow analysis.
Applsci 09 01216 g002
Figure 3. (a) Geometry mesh, (b) zoom of the mesh of the inlet/outlet surface.
Figure 3. (a) Geometry mesh, (b) zoom of the mesh of the inlet/outlet surface.
Applsci 09 01216 g003
Figure 4. Blood velocity field distribution inside the vessel at different times; (a) t = 0.1 s (b) t = 0.2 s (c) t = 0.4 s (d) t = 0.6 s (e) t = 0.7 s and (f) t = 0.8 s.
Figure 4. Blood velocity field distribution inside the vessel at different times; (a) t = 0.1 s (b) t = 0.2 s (c) t = 0.4 s (d) t = 0.6 s (e) t = 0.7 s and (f) t = 0.8 s.
Applsci 09 01216 g004
Figure 5. Temporal evolution of the maximal velocity.
Figure 5. Temporal evolution of the maximal velocity.
Applsci 09 01216 g005
Figure 6. Pressure distribution inside the vessel at different times; (a) t = 0.1 s and (b) t = 0.7 s.
Figure 6. Pressure distribution inside the vessel at different times; (a) t = 0.1 s and (b) t = 0.7 s.
Applsci 09 01216 g006
Figure 7. Temporal evolution of the maximal pressure.
Figure 7. Temporal evolution of the maximal pressure.
Applsci 09 01216 g007
Figure 8. Displacements in the blood vessel using a hyperelastic model; (a) t = 0.1 s and (b) t = 0.7 s.
Figure 8. Displacements in the blood vessel using a hyperelastic model; (a) t = 0.1 s and (b) t = 0.7 s.
Applsci 09 01216 g008
Figure 9. Temporal evolution of the maximal displacement.
Figure 9. Temporal evolution of the maximal displacement.
Applsci 09 01216 g009
Table 1. Mechanical parameters of the artery.
Table 1. Mechanical parameters of the artery.
ParameterValue
Density ρ a r t e r y = 960 kg/m3
Shear modulus μ s = 62.04 × 10 5 N/m2
Bulk modulus κ s = 124 × 10 5 N/m2
Table 2. Finite-element mesh sizes used in this simulation.
Table 2. Finite-element mesh sizes used in this simulation.
Fluid Elements Solid ElementsTotal ElementsTotal Nodes
85,083109,967195,05033,506

Share and Cite

MDPI and ACS Style

Selmi, M.; Belmabrouk, H.; Bajahzar, A. Numerical Study of the Blood Flow in a Deformable Human Aorta. Appl. Sci. 2019, 9, 1216. https://doi.org/10.3390/app9061216

AMA Style

Selmi M, Belmabrouk H, Bajahzar A. Numerical Study of the Blood Flow in a Deformable Human Aorta. Applied Sciences. 2019; 9(6):1216. https://doi.org/10.3390/app9061216

Chicago/Turabian Style

Selmi, Marwa, Hafedh Belmabrouk, and Abdullah Bajahzar. 2019. "Numerical Study of the Blood Flow in a Deformable Human Aorta" Applied Sciences 9, no. 6: 1216. https://doi.org/10.3390/app9061216

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