Document Type: Research Paper


1 Department of Chemical Engineering, Jundi-Shapur University of Technology, Dezful, Iran

2 Department of Chemical Engineering, Amirkabir University of Technology, Tehran, Iran



One of the critical limitations of studies on cardiovascular blood flow simulation is to determine outlet boundary conditions accurately. In the present study, for the first time, pore network model is proposed as a useful technique to take into account interaction between blood flow and other body organs. Thus body organs are simulated by pore network model. Thanks to the method, pressure distribution among the porous medium of organ is determined and consequently the required boundary conditions are obtained for the simulation of arterial blood flow. The comparison between permeability resulted from developed model and experimental results shows that the difference is about 3% for the assumption of non-Newtonian blood flow through organ. This indicates the pore network model can accurately simulate velocity and pressure in the organs. Afterwards, a 3D patient-specific abdominal aorta was simulated under the proposed outlet boundary condition. The maximum deviation of predicted pressure from physiological data is 11.14% near the systole instant. Generally, the predicted pressure and velocity profiles are evident that the model can adequately simulate the blood flow through the arteries which feed main organs.

Graphical Abstract


Main Subjects

[1] A.D. Caballero and S. Lain, “A Review on Computational Fluid Dynamics Modelling in Human Thoracic Aorta”, Cardiovascular Engineering and Technology, Vol 4, No. 2, pp. 103-130, (2013).

[2] J.Y. Park, C.Y. Park, C.M. Hwang, K. Sun and B.G. Min, “Pseudo-organ boundary conditions applied to a computational fluid dynamics model of the human aorta”, Computers in Biology and Medicine, Vol 37, No. 8, pp. 1063-1072, (2007).

[3] A.C. Benim, A. Nahavandi, et al., “Simulation of blood flow in human aorta with emphasis on outlet boundary conditions”, Appl. Math. Modell., Vol 35, No. 7, pp. 31753188, (2011).

[4] E.O. Kung, A.S. Les, et al., “In vitro validation of finite element model of AAA hemodynamics incorporating realistic outlet boundary conditions”, J. Biomech. Eng., Vol. 133, No. 4, pp. 041003, (2011).

[5] I.E. Vignon-Clementel, C.A. Figueroa, K.E. Jansen and C. Taylor, “Outflow boundary conditions for three dimensional finite element modelling of blood flow and pressure in arterries”, Comput. Methods Appl. Mech. Eng., Vol. 195, No. 29-32, pp. 37763796, (2006).

[6] A.B. Deoghare, “Modelling of Human Abdominal Artery for Blood flow Analysis”, Materials Today: Proceedings, Vol 5, No. 5, pp. 12877-12885, (2018).

[7] M. Oshima, R. Torii, et al., “Finite element simulation of blood flow in the cerebral artery”, Comput. Methods Appl. Mech. Eng., Vol. 191, No. 6-7, pp. 661-671, (2001).

[8] C.A. Taylor, T.J.R. Hughes and C.K. Zarins, “Finite element modelling of blood flow in arteries”, Comput. Methods Appl. Mech. Eng., Vol. 158, No. 1-2, pp. 155-196, (1998).

[9] M. Borse, S. Bhushan, D. Walters, G. Burgreen, “Numerical simulations of flow pattern and particle trajectories in feline aorta for hypertrophic cardiomyopathy heart conditions”, Engineering Applications of Computational Fluid Mechanics, Vol.12, No. 1, pp. 57-73, (2018).

[10] J. Renner, R. Gardhagen, et al., “A method for subject specific estimation of aortic wall shear stress”, WSEAS Trans. Biol. Biomed., Vol. 6,  No. 3, pp. 49-57, (2009).

 [11] F.N.V.d. Vosse, J.d. Hart, et al., “Finite element based computational methods for cardiovascular fluid structure interaction”, J. Eng. Math., Vol. 47, No. 3-4, pp. 335-368, (2003).

[12] D. Gallo, G.D. Santis, et al. “On the use of in vivo measured flow rates as boundary conditions for image-based hemodynamic models of the human aorta: implications for indicators of abnormal flow”, Ann. Biomed. Eng., Vol. 40, No. 3, pp.  729-741, (2012). 

[13] N. Shahcheraghi, H.A. Dwyer, A.Y. Cheer, A.I. Barakat and T. Rutaganira, “Unsteady and three-dimensional simulation of blood flow in the human aortic arch”, J. Biomech. Eng., Vol. 124, No. 4, pp. 378-387, (2002). 

[14] E. Shalman, M. Rosenfeld, E. Dgany and S. Einav, “Numerical modeling of the flow in stenosed coronary artery. The relationship between main hemodynamic parameters”, Comput. Biol. Med., Vol. 32, No. 5, pp.  329-344, (2002).

[15] H.J. Kim, I.E. Vignon-Clementel, et al., “On coupling a lumped parameter heart model and a three-dimensional finite element aorta model”, Ann. Biomed. Eng., Vol. 37, No. 11, pp.  2153-2169, (2009). 

[16] M. Dadvar, M. Sohrabi and M. Sahimi, “Pore network model of deactivation of immobilized glucose isomerase in packed-bed reactors I: Two-dimensional simulations at the particle level”, Chemical Engineering Science, Vol. 56, No. 8, pp. 2803-2819, (2001).

[17] M. Dadvar and M. Sahimi, “The effective diffusivities in porous media with and without nonlinear reactions”, Chemical Engineering Science, Vol. 62, No. 5, pp. 1466-1476, (2007). 

[18] M. Dadvar and M. Sahimi, “Pore network model of deactivation of immobilized glucose isomerase in packed-bed reactors. Part III: Multiscale modelling”, Chemical Engineering Science, Vol. 58, No. 22, pp. 4935-4951, (2003). 

[19] M. Dadvar and M. Sahimi, “Pore network model of deactivation of immobilized glucose isomerase in packed-bed reactors II: three-dimensional simulation at the particle level”, Chemical Engineering Science, Vol. 57, No. 6, pp. 939-952, (2002). 

[20] E.N. Marieb and K. Hoehn, Anatomy and physiology, 3rd ed., San Francisco, (2008). 

[21] C. Debbaut, J. Vierendeels, et al., “Perfusion characteristics of the human hepatic microcirculation based on 3D reconstructions and computational fluid dynamic analysis”, J Biomech Eng, Vol. 134, No. 1, pp. 011003-011010, (2012). 

[22] S. Karimi, M. Dabagh, P. Vasava, M. Dadvar, B. Dabir and P. Jalali, “Effect of rheological models on the hemodynamics within human aorta”, Journal of Non-Newtonian Fluid Mechanics, Vol. 207,  pp. 42-54, (2014).

[23] J. Boyd, J.M. Buick and S. Green, “Analysis of the Casson and Carreau-Yasuda non-Newtonian blood models in steady and oscillatory flows using the lattice Boltzmann method”, Phys. Fluids, Vol. 19, No. 9, pp. 093103, (2007).

[24] G.P. Biro, “Comparison of acute cardiovascular effects and oxygen-supply following haemodilution with dextran, stroma-free haemoglobin solution and fluorocarbon suspension”, Cardiovasc. Res., Vol. 16, No. 4, pp. 192-204, (1982). 

[25] E.W. Merrill, “Rheology of blood”, Physiol. Rev., Vol. 49, No. 4, pp. 683-888, (1969). 

[26] R. Skalak, S.R. Keller and T.W. Secomb, “ASME centennial historical perspective paper: mechanics of blood flow”, J. Biomech. Eng., Vol. 103, No. 2, pp. 102-115, (1981).

[27] R.B. Bird, W.E. Stewart, E.N. Lightfoot, Transport phenomena, John Wiley & Sons, (2007).

[28] J.E. Moore and D.N. Ku, “Pulsatile velocity measurements in the model of human abdominal aortic under simulated exercise and postprandial conditions”, Journal of Biomechanical Engineering, Vol. 116, No. 1, pp. 107-111, (1994).

[29] T. Shipkowitz, V.G. Rodgers, J. Lee and K.B. Chandran, “Numerical study on the effects of secondary flow in the human aorta on the local shear stresses in the abdominal aortic branches”, Journal of biomechanics, Vol. 33, No. 6, pp. 717-728, (2000).

[30] S.M.A. Khader, et al., “Haemodynamics Behaviour in Normal and Stenosed Renal Artery using Computational Fluid Dynamics”, Journal of Advanced Research in Fluid Mechanics and Thermal Sciences, Vol. 51, No. 1, pp. 80-90, (2018).

[31] C.C. Botar, Á.Á. Tóth, et al., “Dynamic simulation and Doppler Ultrasonography validation of blood flow behavior in Abdominal Aortic Aneurysm”, Physica Medica, Vol. 37, pp. 1-8, (2017).

 [32] D. Elad, S. Einav, Physical and flow properties of blood, McGraw-Hill Companies, (2009).

[33] C.J. Mills, I.T. Gabe, et al., “Pressure-flow relationships and vascular impedance in man”, Cardiovasc. Res., Vol. 4, No. 4, pp. 405-417, (1970).

[34] A.A. Soares, S. Gonzaga, C. Oliveira, A. Simões, A.I. Rouboa, “Computational fluid dynamics in abdominal aorta bifurcation: non-Newtonian versus Newtonian blood flow in a real case study”, Computer methods in biomechanics and biomedical engineering, Vol. 20, No. 8, pp. 822-831, (2017).

[35] F. Abdulrahim, “A Three-Dimensional Transient Numerical Study and Analysis of Blood Flow Dynamics in Renal Artery Aneurysms”, University of Central Oklahoma, (2017).

[36] D. Lee and J.Y. Chen, “Numerical simulation of steady flow fields in a model of abdomial aorta with its peripheral branches”, journal of biomechanics, Vol. 35, No. 8, pp. 1115-1122, (2002).

[37] T. Mabotuwana, L.K. Cheng and A.J. Pullan, “A model of blood flow in the mesentric arterial system”, Biomedical Engineering OnLine, Vol. 6, No. 1, pp. 17-29, (2007).