Multidimensional modelling for the carotid artery blood flow. (English) Zbl 1178.76395

Summary: A multidimensional 3D-1D FEM model of the whole arterial tree is implemented. It comprises a 3D compliant model of the carotid bifurcation coupled with a 1D model for the remaining part of the arterial tree. With this approach, difficulties arising from the treatment of boundary conditions for the 3D model are naturally handled. The Navier-Stokes equations are used as the governing equations for the blood flow while an elastic compliant model is implemented for the arterial wall. Also, the A.L.E. formulation is considered within the 3D blood regions taking into account the domain deformations produced by the moving artery wall. This 3D model is complemented with a 1D model of the entire arterial tree, in order to appropriately set inflow and outflow boundary conditions for the former. The reduced 1D model solves the momentum and continuity equations in compliant tubes so as to reproduce the propagation of the pressure pulse in the arterial network. Also, a volumetric flow rate is imposed at the inlet to model the systolic work of the heart. The peripheral arteriole beds are simulated with the well-known lumped Windkessel model. A standard geometry of the carotid bifurcation is discretized with P1 bubble-P1 tetrahedral elements. The obtained results properly reproduce the general flow patterns reported in the literature. Very good agreement between the outcomes of a pure 1D model and those of the combined multidimensional model was obtained. It is worth noting that this kind of model may provide useful information to gain understanding in the genesis and development of arterial diseases.


76Z05 Physiological flows
76M10 Finite element methods applied to problems in fluid mechanics
92C10 Biomechanics
74F10 Fluid-solid interactions (including aero- and hydro-elasticity, porosity, etc.)
Full Text: DOI


[1] Bruinsma, P.; Arts, T.; Dankelman, J.; Spaan, J.A.E., Model of the coronary circulation based on pressure dependence of coronary resistance and compliance, Basic res. cardiol., 83, 510-524, (1988)
[2] Caro, C.G.; Fitz-Gerald, J.M.; Schroter, R.C., Atheroma and arterial wall shear dependent mass transfer mechanism for atherogenesis, Proc. roy. soc. lond. biol. B, 177, 109-159, (1971)
[3] Rindt, C.C.M.; Steenhoven, A.A.; Van Janssen, J.D.; Reneman, R.S.; Segal, A., A numerical analysis of steady flow in a three-dimensional model of the carotid artery bifurcation, J. biomech., 23, 461-473, (1990)
[4] Ku, D.N.; Giddens, D.P., Hemodynamics of the normal carotid bifurcation, Ultrasound med. biol., 11, 13-26, (1985)
[5] Ku, D.N.; Giddens, D.P., Laser Doppler anemometer measurements of pulsatile flow in a model carotid bifurcation, J. biomech., 20, 407-421, (1987)
[6] Ku, D.N.; Giddens, D.P.; Zarins, C.K.; Glagov, S., Pulsatile flow and atherosclerosis in the human carotid bifurcation, Arteriosclerosis, 5, 293-302, (1985)
[7] Botnar, R.H.; Rappitsch, G.; Scheidegger, M.B.; Liepsch, D.; Perktold, K.; Boesiger, P., Hemodynamics in the carotid artery bifurcation: a comparison between numerical simulations and in vitro MRI measurements, J. biomech., 33, 137-144, (2000)
[8] Sonneveld, P., Multigrid and conjugate gradient methods as convergence acceleration techniques, Multigrid methods in integral differential equations, (1985), Clarendon Press Oxford · Zbl 0577.65086
[9] Y. Saad, SPARSEKIT: a basic tool kit for sparse matrix computation version 2. Available from: <http//:www-users.cs.umn.edu/ saad/software/SPARSEKIT/sparsekit.html>, University of Illinois, 1985.
[10] Urquiza, S.A.; Vénere, M.J., An application framework architecture for FEM and other related solvers, (), 3099-3109
[11] S.A. Urquiza, A.L. Reutemann, M.J. Vénere, R.A. Feijóo, Acoplamiento de modelos unidimensionales y multidimensionales para la resolución de problemas hemodinámicos, in: Proc. ENIEF 2001, XII Congress on Numerical Methods and Their Applications, AMCA, Córdoba, 2002.
[12] S.A. Urquiza, M.J. Vénere, F.M. Clara, R.A. Feijóo, Finite element (one-dimensional) hemodynamic model of the human arterial system, in: E. Onate, G. Bugeda, B. Suarez, (Eds.), Proc. ECCOMAS 2000, European Community on Computational Methods in Applied Sciences and Engineering, Artes Gráficas Torres S.A., Barcelona, 2000.
[13] Giddens, D.P.; Zarins, C.K.; Glagov, S., The role of fluid mechanics in the localization and detection of atherosclerosis, J. biomech. engrg., 115, 588-594, (1993)
[14] Anayiotos, A.S.; Jones, S.A.; Giddens, D.P.; Glagov, S.; Zarins, C.K., Shear stress at a compliant model of the human carotid bifurcation, J. biomech. engng., 117, 98-106, (1994)
[15] Pironneau, O., Finite element methods for fluids, (1989), Wiley New York · Zbl 0665.73059
[16] Quarteroni, A., Modelling the cardiovascular system: a mathematical challenge, (), 961-972 · Zbl 1036.92013
[17] Formaggia, L.; Gerbeau, J.F.; Nobile, F.; Quarteroni, A., On the coupling of 3D and 1D navier – stokes equations for flow problems in compliant vessels, Comput. methods appl. mech. engrg., 191, 561-582, (2001) · Zbl 1007.74035
[18] Formaggia, L.; Gerbeau, J.F.; Nobile, F.; Quarteroni, A., Numerical treatment of defective boundary conditions for the navier – stokes equations, SIAM J. numer. anal., 40, 376-401, (2002) · Zbl 1020.35070
[19] R. Pietrabissa, A. Quarteroni, G. Dubini, A. Veneziani, F. Migliavacca, S. Ragni, From the global cardiovascular hemodynamics down to the local blood motion: preliminary applications of a multiscale approach, in: E. Oate, G. Bugeda, B. Suarez, (Eds.), Proc. ECCOMAS 2000, European Community on Computational Methods in Applied Sciences and Engineering, Artes Gráficas Torres S.A., Barcelona, 2000.
[20] Quarteroni, A.; Tuveri, M.; Veneziani, A., Computational vascular fluid dynamics: problems, models and methods, Comput. visual. sci., 2, 163-197, (2000) · Zbl 1096.76042
[21] Formaggia, L.; Nobile, F.; Quarteroni, A.; Veneziani, A., Multiscale modelling of the circulatory system: a preliminary analysis, Comput. visual. sci., 2, 75-83, (1999) · Zbl 1067.76624
[22] F. Nobile, Numerical approximation of fluid – structure interaction problems with application to haemodynamics, Ph.D. Thesis, Report Number 2458, EPFL, Lausanne, 2001.
[23] Guillard, H.; Farhat, C., On the significance of the geometric conservation law for computations on moving meshes, Comput. methods appl. mech. engrg., 190, 1467-1482, (2000) · Zbl 0993.76049
[24] Koobus, B.; Farhat, C., Second-order time-accurate and geometrically conservative implicit schemes for computations on unstructured dynamic meshes, Comput. methods appl. mech. engrg., 170, 103-129, (1999) · Zbl 0943.76055
[25] M. Lesoinne, C. Farhat, Geometric conservation laws for aeroelastic computations using unstructured dynamics meshes, in: 12th AIAA Computational Fluid Dynamics Conference (AIAA-95-1709, San Diego, June 1995). · Zbl 0896.76044
[26] Stettler, J.C.; Niederer, P.; Anliker, M., Theoretical analysis of arterial hemodynamics including the influence of bifurcations part I, Ann. biomed. engrg., 9, 145-164, (1981)
[27] Perktold, K.; Rappitsch, G., Computer simulation of local blood flow and vessel mechanics in a compliant carotid artery bifurcation model, J. biomech., 28, 845-856, (1995)
[28] Snyder, M.F.; Rideout, V.C.; Hillestad, R.J., Computer modelling of the human systemic arterial tree, J. biomech. engrg., 1, 341-353, (1968)
[29] Avolio, A.P., Multi-branched model of the human arterial system, Med. biol. engrg. comput., 18, 709-718, (1980)
[30] Stergiopulos, N.; Young, D.F.; Rogge, T.R., Computer simulation of arterial flow with applications to arterial and aortic stenoses, J. biomech., 25, 1477-1488, (1992)
[31] Schaaf, B.W.; Abrecht, P.H., Digital computer simulation of human systemic arterial pulse wave transmission: a nonlinear model, J. biomech. engrg., 15, 349-362, (1982)
[32] Hughes, T.J.R.; Liu, W.K.; Zimmermann, T.K., Lagrangian – eulerian finite element formulation for incompressible viscous flows, Comput. methods appl. mech. engrg., 29, 329-349, (1981) · Zbl 0482.76039
[33] T.J.R. Hughes, A study of the one-dimensional theory of arterial pulse propagation, Ph.D. Thesis, Report Number 74-13, Structures and Materials Research, Department of Civil Engineering, Division of Structural Engineering and Structural Mechanics, University of California, Berkeley, 1974.
[34] Hughes, T.J.R.; Brooks, A., A theoretical framework for petrov – galerkin methods with discontinuous functions: application to the stream-upwind procedure, (), 46-65
[35] Taylor, C.A.; Hughes, T.J.R.; Zarins, C.K., Finite element blood flow modeling in arteries, Comput. methods appl. mech. engrg., 158, 155-196, (1998) · Zbl 0953.76058
[36] Hughes, T.J.R.; Franca, L.P.; Hulbert, G.M., A new finite element formulation for computational fluid dynamics: VIII. the Galerkin/least squares method for advective – diffusive equations, Comput. methods appl. mech. engrg., 73, 173-189, (1989) · Zbl 0697.76100
[37] Bharadvaj, B.K.; Mabon, R.F.; Giddens, D.P., Steady flow in a model of the human carotid bifurcation—I. flow visualisation, J. biomech., 15, 349-362, (1982)
[38] Bharadvaj, B.K.; Mabon, R.F.; Giddens, D.P., Steady flow in a model of the human carotid bifurcation—II. laser-Doppler anemometer measurements, J. biomech., 15, 363-378, (1982)
[39] Richtmyer, R.D.; Morton, K.W., Difference methods for initial-value problems, (1967), Interscience Publishers New York · Zbl 0155.47502
[40] Courant, R.; Hilbert, D., Methods of mathematical physics, vol. II, (1962), Interscience Publisher New York · Zbl 0729.00007
[41] Carey, G.F.; Jiang, B.N., Least squares finite elements for first-order hyperbolic systems, Int. J. numer. methods fluids, 26, 81-93, (1988) · Zbl 0641.65080
[42] Young, B.N.; Carey, G.F., A stable least-squares FEM for non-linear hyperbolic problems, Int. J. numer. methods fluids, 8, 933-942, (1988) · Zbl 0666.76087
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.