Simulations of time harmonic blood flow in the Mesenteric artery: comparing finite element and lattice Boltzmann methods
© Axner et al; licensee BioMed Central Ltd. 2009
Received: 27 April 2009
Accepted: 2 October 2009
Published: 2 October 2009
Systolic blood flow has been simulated in the abdominal aorta and the superior mesenteric artery. The simulations were carried out using two different computational hemodynamic methods: the finite element method to solve the Navier Stokes equations and the lattice Boltzmann method.
We have validated the lattice Boltzmann method for systolic flows by comparing the velocity and pressure profiles of simulated blood flow between methods. We have also analyzed flow-specific characteristics such as the formation of a vortex at curvatures and traces of flow.
The lattice Boltzmann Method is as accurate as a Navier Stokes solver for computing complex blood flows. As such it is a good alternative for computational hemodynamics, certainly in situation where coupling to other models is required.
Atherosclerosis is the most common cardiovascular disease that affects the arteries [1–6]. It commonly first develops at the origins of tributaries, bifurcations and at curvatures of the arteries, in areas which are associated with low or oscillatory wall shear stress. The complexity of the human vascular system and the time-harmonic character of blood flow make analysis and prediction of its behavior very difficult. To obtain an information of blood flow on the locations of these regions in the human cardiovascular system, numerical studies are essential.
The lattice Boltzmann method (LBM) was proposed two decades ago, and is now recognized as a well-established method in computational fluid dynamics [7, 8]. It has been used to simulate fluid flow in a wide range of complex geometries such as porous media [9–11] or geometries from medical applications [12–19]. Several studies comparing LB and FE methods for transient flows have shown good agreement between, where LBM is less expensive in terms of memory consumption and uses comparable computational times [20, 21]. However, in all these studies, whether in 2D or 3D geometries, the fluid flow was time-independent. Comparison of LB and FE methods for time-harmonic flows is lacking.
We have compared LB and FE simulation results for a systolic flow in a realistic 3D geometry. As an application, we have used blood flow in the abdominal aorta(AA) and one of the major abdominal branches, the superior mesenteric artery (SMA). The SMA takes origin from the anterior surface of the AA, distal to the root of the celiac trunk. Diffuse atherosclerotic disease is rare in the SMA its occlusion leads to the death of patients.
For quantitative comparison between the two methods we used the velocity and pressure profiles of the simulated flows. We also studied the details of flow characteristics such as vortex formation and traces of flow.
Lattice Boltzmann method
b is the number of directions.
Finite element method
As described in Jeays et al. , a transient CFD model was constructed by morphing a parametric mesh constructed from simple geometric primitives. The advantage of this process is that it is easy to control the element size distribution mapped onto the original geometry. It is robust in operation, and is ideally suited to the generation of dynamic CFD meshes of arterial systems that are free from major pathology. An unstructured triangular mesh was generated for the CFD analysis using the ANSYS pre-processor. Flow boundary conditions were determined based on phase contrast MRI velocity measurements. FLOTRAN (ANSYS Inc.) was used as a FEM solver for simulating the 3D Navier-Stokes equations as (see Jeays et al. ).
Results and Discussion
Experiments and validation
The geometry was obtained from an Magnetic Resonance Imaging scan obtained from a healthy volunteer with full ethical approval. A triangular mesh was generated for use in the FEM solver . The 3D mesh for FEM comprised 100465 nodes, with 650 nodes on the AA inlet, 295 on the AA outlet and 207 on the SMA outlet. From this mesh a voxel mesh was generated for LBM using a special 3D editing and mesh generation tool HemoSolve . In this tool the parts of interest of geometry can be selected and can be enhanced with inlet and outlet layers on its end-points. Hemosolve has sufficient functionality to mimic a real surgical procedure. The generated voxel mesh from FEM triangular mesh was comprised of a total of 74468 fluid nodes with 694 nodes on the AA inlet, 316 on the AA outlet and 224 on the SMA outlet. For these simulations the walls of geometries are considered rigid as the effect of the elasticity is proven to be less than 3% (see ref. ).
The following flow parameters were applied in FEM: the density of the blood ρ was 1000 kgm-3, Newtonian viscosity μ was 4 * 10-3 Pa.s, the cardiac cycle duration T was 0.86 s with 5 ms time-steps δt and the blood was considered incompressible. The maximum velocity u of 0.8 m/s at a peak systole in the aorta, combined with the diameter D of 0.0165 m, gave a maximum Reynolds number (Re) of 3300 . In order to apply the same flow condition in the LBM solver we converted all parameters into dimensionless numbers. From the FEM parameters we computed a Womersley number of α = 11.14. α is the ratio between the time-harmonic flow frequency and viscous effects and is defined as .
Using these two constraints (Re = 3300 and α = 11.14) and aiming at no more than 10% simulation error (see ) we chose the maximum velocity u max = 0.08. Next, we minimized the execution time by applying the constraint optimization scheme, using the same procedure as described in  where we assume that the parametrization used for pipe flow was also applicable to the current geometry. This gave T = 17200, e.g. in physical units a δt = 5 * 10-5 time step, while D = 30 lattice points. Finally from these parameters we derived ν = 0.00068.
We applied the same inlet/outlet boundary conditions as described in Jeays et al. . Velocities (created using Womersley's solution ) were specified at the proximal AA opening, the pressure waveform (created using Westerhof's model ) was specified at the distal AA opening and a free-flow boundary condition was applied at the outflow of the SMA. We first transformed these given velocities into dimensionless values and then applied the velocity and pressure boundary condition as described in Ref.  at the proximal AA and distal AA openings respectively.
In this study, for LBM a bounce back on links (BBL) boundary condition is applied on the walls. We did not use more accurate boundary condition such as Bouzidi boundary conditions (BBC) . These demand more elaborate computations especially for irregular geometries connected with the alpha parameter. From previous studies the order of both methods in terms of grid spacing are known [23, 25]. Moreover, many formulations of solid boundary conditions are known and their influence on the accuracy on the flow fields is well understood [12, 25]. We expect however that using BBC for LBM will only improve the results presented here, in the sense that either the accuracy will increase, or the choice of optimal simulation parameters changes drastically . That is however not the topic of this paper.
Comparison between two methods
After running simulations using both methods we compared the resulting velocity profiles at three different transverse planes (A, B, and C) as shown in Fig. 1.
There is very good agreement between velocity profiles obtained by FEM and LBM. Moreover, the position of peak velocity in Fig. 2 and Fig. 3, e.g. before the bifurcation, is skewed to the left, the posterior aspect of the vessel, while in Fig. 4, e.g. after bifurcation, it is towards the right, the anterior aspect of the vessel. This indicates a spiraling of the flow in the AA.
where a is either diameter D or time t. The greatest average difference is observed close to the vessel walls (Fig. 6) and is approximately 10% while the maximum difference is not more than 13%. In Fig. 7 we see that the greatest average difference as a function of time is observed after peak systole and is about 6% while the maximum is about 13%.
In Fig. 10 we also show the velocity profiles at peak systole in the last cardiac cycle executed, when beat to beat convergence had been achieved, along the cutting plane shown in Fig. 1 (right). Here the expected flow behavior and the formation of a vortex at the top of SMA just below the outer wall can be clearly seen. For FEM the vortex first appeared at 0.26 seconds and disappeared at 0.56 seconds while, for LBM, it appeared at time-step 5200 and disappears at time-step 11200. Samples were taken every 100 time-steps and inspected visually. These results are completely identical if the δt = 5 * 10-5 time-step is considered for the LBM simulation.
We have validated the LBM for time-harmonic flows by comparing the simulation results with those obtained from numerically solving the Navier Stokes equations with FEM. As an experimental geometry we used the geometry of the human vascular system, the AA together with one of its major abdominal branches, the SMA. We compared velocity and pressure profiles of simulated time-harmonic blood flow for both methods and demonstrated a very good agreement. The maximum differences for velocity profiles were greatest next to the vessel walls. These were less than 10%. The maximum difference, as a function of time, was 6%. Moreover, the spiraling of the flow profiles in the AA and the time of vortex formation in SMA coincided with both methods.
The main purpose of this paper is to demonstrate that the LBM is a good alternative for computational hemodynamics. It is self evident that validation of a this new method against existing established methodologies (that is, using Navier-Stokes solvers) is an important step. We have drawn the Navier Stokes solution data for comparison from previous work by Jeays et al. . We have chosen to focus our own efforts on demonstrating that the LB solver can reproduce the published Navier Stokes analysis by Jeays et al.  of a complex transient flow pattern in an arterial bifurcation, with the flow domain constructed from a medical image and the boundary conditions based on in vivo measurements. Some papers (e.g Ku et al , Ford et al.  have focused on the validation of Navier-Stokes solvers in the context of arterial flows, both by direct measurement in vivo and by measurement in anatomically-realistic experimental phantoms. In contrast, the emphasis in this manuscript is whether the numerical simulations produce an accurate solution of the problem posed. It is very well established that the most important determinant of any arterial flow, particularly in terms of the important local features such as flow separations, is the geometry. The whole point of the current work is to demonstrate that these local details are comparable between the LB solver and an established Navier Stokes solver. This paper confirms that the LBM represents a viable alternative to the Navier Stokes solvers for complex transient flow, in this case in an arterial bifurcation. It was quite remarkable that LBM stayed stable for such large Re = 3300 numbers. As for the computational expenses, LBM is less memory consuming and reported (see ref. ) to have computational times comparable to the FEM. With the modified simulator (see ref. ) we expect LBM to perform even faster.
With the increasing interest in using LBM for computational hemodynamics [15, 17, 18] this study shows that LBM can be considered to be an alternative as a solver for computational hemodynamics, producing results of equal quality to Navier Stokes solvers.
Another main difference between both methods lies in the mesh generation, which for LBM is obviously much easier. Moreover, suggested by Bernsdorf in his thesis , the LBM seems to be very well suited for the use in multiphysics models, for instance in blood clotting studies .
The work was funded by the Dutch National Science Foundation, NWO, Token 2000 Distributed Interactive Medical Exploratory for 3D Medical Images (DIME) project (634.000.024 ). Special thanks to Michael Scarpa from the University of Amsterdam for help with visualization of the geometries.
- Dampney RRA: Functional Organization of Central Pathways Regulating the Cardiovascular System. Physiological Reviews 1994,74(2):323–364.Google Scholar
- Wald NJ, Law MR: A Strategy to Reduce Cardiovascular Disease by more than 80%. BMJ Publishing Group 2003, 326: 1419–1423.Google Scholar
- Thubrikar MJ: Vascular Mechanics and Pathology. Spinger Science+Business Media; 2007.View ArticleGoogle Scholar
- Taylor CA, Hughes TJ, Zarins CK: Finite Element Modeling of Blood Flow in Arteries. Computer Methods in Applies Mechanics and Engineering 1998, 158: 155–196. 10.1016/S0045-7825(98)80008-XMathSciNetView ArticleGoogle Scholar
- Taylor CA, Draney MT: Experimental and Computational Methods in Cardiovascular Fluid Mechanics. Annual Review Fluid Mechanics 2004, 36: 197–231. 10.1146/annurev.fluid.36.050802.121944MathSciNetView ArticleGoogle Scholar
- Vorp DA: Biomechanics of abdominal aortic aneurysm. J of Biom 2007,44(9):1887–1902. 10.1016/j.jbiomech.2006.09.003View ArticleGoogle Scholar
- Succi S: The Lattice Boltzmann Equation for Fluid Dynamics and Beyond. New York:Oxford; 2001.Google Scholar
- Chen S, Doolen GD: Lattice Boltzmann Method for Fluid Flows. Annual Review Fluid Mechanics 1998, 30: 329–364. 10.1146/annurev.fluid.30.1.329MathSciNetView ArticleGoogle Scholar
- Koponen A, Kandhai B, Héllen E, Alava M, Hoekstra A, Kataja M, Niskanen K, Sloot P, Timonen J: Permeability of three-dimensional random fibre webs. Physical Review Letters 1998, 80: 716–719. 10.1103/PhysRevLett.80.716View ArticleGoogle Scholar
- Clague D, Kandhai D, Zhang R, Sloot PM: On the Hydraulic Permeability of (Un)Bounded Fibrous Media Using the Lattice-Boltzmann Method. Physical Review E 2000, 61: 616–625. 10.1103/PhysRevE.61.616View ArticleGoogle Scholar
- Kandhai D, Hlushkou D, Hoekstra AG, Sloot PM, van As H, Tallarek U: Influence of Stagnant Zones on Transient and Asymptotic Dispersion in Macroscopically Homogeneous Porous Media. Physical Review letters 2002,88(23):234501. 1–4. 10.1103/PhysRevLett.88.234501View ArticleGoogle Scholar
- Artoli AM, Hoekstra AG, Sloot PM: Mesoscopic Simulations of Systolic Flow in the Human Abdominal Aorta. Journal of Biomechanics 2006,39(5):873–884. 10.1016/j.jbiomech.2005.01.033View ArticleGoogle Scholar
- Artoli AM, Hoekstra AG, Sloot PM: Simulation of a Systolic Cycle in a Realistic Artery with the Lattice Boltzmann BGK Method. International Journal of Modern Physics B 2003,17(1–2):95–98. 10.1142/S0217979203017138View ArticleGoogle Scholar
- Fang H, Wamg Z, Lin Z, Liu M: Lattice Boltzmann Method for Simulating the Viscous Flow in Large Destensible Blood Vessels. Physical Review E 2002, 65: 051925–11. 10.1103/PhysRevE.65.051925View ArticleGoogle Scholar
- Harrison SE, Bernsdorf J, Hose DR, Lawford PV: Development of a Lattice Boltzmann Framework for Numerical Simulation of Thrombosis. Int J of Mod Phys C 2007,18(4):483–491. 10.1142/S0129183107010711View ArticleGoogle Scholar
- Hirabayashi M, Ohta M, Rüfenacht D, Chopard B: Numerical analysis of the flow pattern effect on the flow reduction performance in the cerebral aneurysm by stent implantation. Mathematics and Computers in Simulation 2006, 72: 128. 10.1016/j.matcom.2006.05.037MathSciNetView ArticleGoogle Scholar
- Chopard B, Ouared R, Ruefenacht D, Yilmaz H: Lattice Boltzmann Modeling of Thrombosis in Giant Aneurysms. Int J Mod Phys C 2007,18(4):712–721. 10.1142/S0129183107010978View ArticleGoogle Scholar
- Bernsdorf J, Harrison S, Hose R, Lawford P, Smith S: Numerical Simulation of Clotting Processes: A Lattice Boltzmann Application in Medical Physics. Mathematics and Computers in Simulation 2006,72(2–6):89–92. 10.1016/j.matcom.2006.05.008MathSciNetView ArticleGoogle Scholar
- Krafczyk M, Schulz M, Rank E, Cerrolaza M: Simulation and visualization of transient 3d-flow in arteries with an articial heart valve using lattice-boltzmann method. Journ of Biomech 1998,31(453):462.Google Scholar
- Kandhai B, Vidal DE, Hoekstra A, Hoefsloot H, Iedema P, Sloot P: Lattice-Boltzmann and Finite Element Simulations of Fluid Flow in a SMRX Mixer. International Journal of Numerical Methethods in Fluids 1999, 31: 1019–1033. 10.1002/(SICI)1097-0363(19991130)31:6<1019::AID-FLD915>3.0.CO;2-IView ArticleGoogle Scholar
- Geller S, Krafczyk M, Toelke J, Turek S, Hron J: Benchmark Computations Based on Lattice-Boltzmann, Finite Element and Finite Volume Methods for Laminar Flows. Computers and Fluids 2006, 35: 888–897. 10.1016/j.compfluid.2005.08.009View ArticleGoogle Scholar
- Qian YH, d'Humieres D, Lallemand P: Lattice BGK Models for Navier-Stokes Equation. Europhys Letters 1992, 17: 479–484. 10.1209/0295-5075/17/6/001View ArticleGoogle Scholar
- Jeays A, Lawford P, Gillott R, Spencer P, Barber D, Bardhan K, Hose D: Characterisation of the Haemodynamics of the Superior Mesenteric Artery. Journal of Biomechanics 2006,40(9):1916–1926. 10.1016/j.jbiomech.2006.09.009View ArticleGoogle Scholar
- Axner L, HA Marquering AH, Geest R, Reiber J, Sloot P: Towards Decision Support in Vascular Surgery through Computational Hemodynamics. Computers in Biology and Medicine 2008, in press.Google Scholar
- Axner L, Hoekstra A, Sloot P: Simulating Time Harmonic Flows with the Lattice Boltzmann Method. Physical Review E 2007,75(3):036709+7. 10.1103/PhysRevE.75.036709View ArticleGoogle Scholar
- Womersley J: Method for the Calculation of Velocity, Rate of Flow and Viscous Drag in Arteries when the Pressure Gradient is Known. Journal of Physiology 1955,127(3):553–563.View ArticleGoogle Scholar
- Westerhof N, Bosman F, de Vries C, Noordergraaf A: Analog Studies of the Human Systemic Arterial Tree. Journal of Biomechanics 1969, 2: 121–143. 10.1016/0021-9290(69)90024-4View ArticleGoogle Scholar
- Zou Q, He X: On pressure and velocity boundary conditions for the lattice Boltzmann BGK model. Phys of Fluids 1997,9(6):1591–1598. 10.1063/1.869307MathSciNetView ArticleGoogle Scholar
- Bouzidi M, Firdaouss M, Lallemand P: Momentum transfer of a lattice-Boltzmann fluid with boundaries. Phys Fluids 2001, 13: 3452. 10.1063/1.1399290View ArticleGoogle Scholar
- Ku J, Draney M, Arko F, Lee W, Chan F, Pelc N, Zarins C, Taylor C: In Vivo Validation of Numerical Prediction of Blood Flow in Arterial Grafts. Annals of Biomedical Engineering 2002,30(6):743–752. 10.1114/1.1496086View ArticleGoogle Scholar
- Ford M, Nikolov H, Lownie S, DeMont E, Kalata W, Loth F, Holdsworth D, Steinman D: PIV-Measured Versus CFD-Predicted Flow Dynamics in Anatomically-Realistic Cerebral Aneurysm Models. Journal of Biomechanical Engineering 2008,130(2):021015. 10.1115/1.2900724View ArticleGoogle Scholar
- Axner L, Bernsdorf J, Zeiser T, Lammers P, Linxweiler J, Hoekstra A: Performance Evaluation of a Parallel Sparse Lattice Boltzmann Solver. Journal of Computational Physics 2007,227(10):4895. 10.1016/j.jcp.2008.01.013MathSciNetView ArticleGoogle Scholar
- Bernsdorf J: Simulation of Complex Flows and Multi-Physics with the Lattice-Boltzmann Method. PhD thesis. Universiteit of Amsterdam; 2006.Google Scholar
- Bernsdorf J, Harrison S, Hose R, Lawford P, Smith S: Concurrent Numerical Simulation of Flow and Blood Clotting Using the Lattice Boltzmann Technique. IEEE Computer Society Press, Volume ICPADS 2005 2005, 336.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.