Skip to main content
Article Open Access

Comparison of stented bifurcation and straight vessel 3D-simulation with a prior simulated velocity profile inlet

  • EMAIL logo , , , and
Published/Copyright: September 30, 2016

Abstract

Coronary diseases are the main reason for death in the western world. Bio-fluid mechanical correlations with arterial diseases are in the focus of our research. To treat occluded vessels, stents are implanted. Stent implantations can be associated with blood flow disruptions leading to restenosis or thrombosis formation. Numerical flow simulation is a promising tool to evaluate the hemodynamic performance of cardiovascular implants, but is resource-intensive in time and computational power. Therefore, a reduction in grid size would be beneficial due to economic exploitation of computational cost. The purpose of this numerical investigation is to substitute the computational domain of a distal stented bifurcation with a stented straight vessel by using the right inlet condition. The deviation of the results of the two different methods to simulate the blood flow situation in a bifurcation is marginal. This inlet can be used for standardised simulations of bifurcations were lesions commonly occur.

1 Introduction

Coronary diseases are the main reason for death in the western world. In Germany one seventh of these cases of death are resulting from a myocardial infarction, which is caused by an occlusion of a coronary artery [1]. The affected arteries are treated with a so-called stent, a scaffold that widens the vessel to ensure an unimpeded blood flow.

Since the introduction of bare metal stents, the development goes to decreasing stent struts, so the blood flow is decreasingly influenced due to the stent segments.

Bioresorbable stents made of materials like polymers, show, as a result of their inferior material properties, higher struts as bare metal stents to ensure strength properties. These larger cross profiles build an obstruction and change the blood flow situation in the vessels. It is common knowledge that disturbances in the blood flow are correlated with diseases like the peripheral artery disease. Stent implantations can be associated with those disruptions that lead to very high and very low shear rates, which can lead to restenosis or thrombosis formation [2], [3]. The regions of interest can be found especially in bifurcations [4].

2 Aim of the simulation

The aim of this simulation is to compare the blood flow situation in a bifurcation and a tube, where a velocity profile taken out of the simulation of a bifurcation is used as an inlet. Thus there would be an enormous saving in time for preparing bifurcation models, generating meshes and in simulation time, because of the lower number of cells.

For this numerical study there will be the following operating points:

  • Modelling of three bifurcations with different angles between the side branches, the generic stent and a tube

  • Definition of the boundary conditions and execution of the simulations

  • Taking the velocity profile off the bifurcation simulation and using it as an inlet for the stented tube

  • Evaluation and comparison of the blood flow fields.

3 Material and methods

3.1 Construction of stented bifurcations and mesh generation

For the construction of the 3D bifurcation model CAD software is used. The inlet length until the junction of the two side branches is 17 mm. The diameter of the main branch is conical and goes from 4.71 mm to 3.5 mm and runs into the symmetrical side branches with a diameter of 3.5 mm and a length of 17 mm each [5]. Three bifurcations are constructed whereby they only differ in the constructed angle between the two side branches. The angles 70°, 90° and 110° are designed, since the average angle between LAD and LCx is about 82° ± 17° [6].

The generic stent is, as well as the bifurcation model, constructed with CAD software. The structure of the stent is helical, comparable to a spring, has a length of 8 mm and is located 12 mm ahead of the outlet in one of the two side branches. The cross section of the stent struts protruding into the lumen of the bifurcation is a semicircle with a height of 150 μm. The stent is cut out of the bifurcation.

The meshes of the three stented bifurcation models are made with open source software OpenFOAM and the tool snappyHexMesh. Each mesh contains 11–12 million volume cells.

To compare the blood flow conditions in the bifurcation model and in a tube with a velocity field taken out of the bifurcation model, a stented tube is needed. The tube is constructed and meshed in OpenFOAM and the stent is cut out with snappyHexMesh. It has about 3 million volume cells.

For each volume cell, the Navier-Stokes-Equation has to be solved (see eq. 1):

(1)ρDuDt=grad p+div (2η(γ˙)D)

where ρ is density, p is pressure, u is velocity vector, t is time, η is absolute viscosity, γ˙ is shear rate and D the tensor of deformation.

3.2 Boundary conditions

Steady-state blood flow is simulated in the bifurcation as well as in the tube model. Furthermore stiff vessel walls and no slip conditions are assumed. At the inlet of the bifurcation model there is a mean velocity of v = 0.2 m/s in a block type profile. This matches a Reynolds number of 250, which is to be found in coronary vessels. The shear thinning behaviour of blood is approximated by the non-Newtonian Carreau model:

μf=μ+(μ0μ)(1+(λγ˙)2)(n1)/2

The used parameters are adopted from [7]. At the outlet a pressure of p = 0 is defined.

3.3 Velocity profile

For the inlet of the tube a velocity profile is taken out of the completely converged simulation of a bifurcation model which is not stented. The velocity profile is taken 2 mm ahead of the virtual beginning of the stent, as well as the inlet of the tube model is 2 mm ahead of the stent (see Figure 1).

Figure 1 (A) Stented bifurcation model and (B) bifurcation model and stented tube model; blue arrows show the position of the obtained velocity profile which is used as an inlet in further simulations.
Figure 1

(A) Stented bifurcation model and (B) bifurcation model and stented tube model; blue arrows show the position of the obtained velocity profile which is used as an inlet in further simulations.

To get the velocity profile all velocity vectors and their positions that lie on an imagined plain normal to the stented side branch are written out. A three dimensional illustration of the mapped velocity profile can be seen in Figure 2.

Figure 2 Mapped velocity profile of the bifurcation with 110 degree angle between the side branches.
Figure 2

Mapped velocity profile of the bifurcation with 110 degree angle between the side branches.

4 Results and discussion

4.1 Results and comparison of the blood flow situations

For comparison of the blood flow situation in the stented bifurcations with various angles and the stented tube models with a prior simulated inlet, there are severe criteria needed. Since restenosis is likely to occur at locations with a wall shear stress below 0.4 Pa, the area of the vessel walls with such shear stress is compared. Further on the maximal velocity and the maximal shear stress that are found in the stented vessel segments are compared.

The difference between the bifurcation and the tube having an area with a wall shear stress below 0.4 Pa is minor. The variation in percentage on all criteria is an average of 6.6%. A comparison of the wall shear stress distribution of the bifurcation with a 110 degree angle is shown in Figure 3.

Figure 3 Wall shear stress distribution: red area has a wall shear stress <0.4 Pa. (A) bifurcation simulation, (B) tube simulation.
Figure 3

Wall shear stress distribution: red area has a wall shear stress <0.4 Pa. (A) bifurcation simulation, (B) tube simulation.

Table 1

Percentage deviation of the stented tube simulation compared to the stented bifurcation simulation.

Area of wall shear stress < 0.4 Pa mm2Max. velocity (m/s)Max. shear stress (s−1)
Bifurcation 70 degree angle1.67 × 10−50.4817370
Tube 70 degree angle1.55 × 10−50.4866980
Bifurcation 90 degree angle1.70 × 10−50.4846460
Tube 90 degree angle1.47 × 10−50.4806470
Bifurcation 110 degree angle1.09 × 10−50.4846480
Tube 110 degree angle1.14 × 10−50.4927970

4.2 Discussion

The analysis showed that the maximum deviation can be found at the comparison between the simulations of the bifurcation of 110 degree angle. This deviation occurs since the 110 degree bifurcation obtains the most complex flux with its vortex.

The deviation of the results of the two different methods to simulate the blood flow situation in a bifurcation is marginal. These slight differences can be the reason of having not completely similar meshes in the stented region due to the automatic meshing tool snappyHexMesh.

In addition, only the mapped velocity profile is used as inlet condition for the simulations of the straight vessels. It could be beneficial to transfer the pressure field or the pressure gradient field as well. Here the described method has its limits. Numerical simulations admit only the velocity or the pressure field as an inlet.

5 Conclusion

This method is successful due to its possibility of saving time and computational resources. Using this inlet in standardised simulations can be a benefit of it. A lot of simulations are made with a block type profile, a Hagen-Poiseuille-profile or symmetric profiles in general. This asymmetric inlet simulates a bifurcation right before the stented vessel segment, were lesions commonly occur and can expose possible stent design optimizations. Further simulations could include the mapping of the gradient pressure field.

Acknowledgement

Financial support by the Federal Ministry of Education and Research (BMBF) within RESPONSE “Partnership for Innovation in Implant Technology” is gratefully acknowledged.

Author’s Statement

Research funding: The author state no funding involved. Conflict of interest: Authors state no conflict of interest. Material and Methods: Informed consent: Informed consent is not applicable. Ethical approval: The conducted research is not related to either human or animal use.

References

[1] Destatis Statistisches Bundesamt, Wiesbaden. Die 10 häufigsten Todesfälle durch Herz-Kreislauf-Erkrankungen. https://www.destatis.de/DE/ZahlenFakten/GesellschaftStaat/Gesundheit/Todesursachen/Tabellen/Sterbefaelle_HerzKreislaufErkrankungenInsgesamt.html. Accessed on April 18, 2016.Search in Google Scholar

[2] Holme PA, Ørvim U, Hamers M, Solum NO, Brosstad FR, Barstad RM, et al. Shear-induced platelet activation and platelet microparticle formation at blood flow conditions as in arteries with a severe stenosis. Arterioscler Thromb Vasc Biol. 1997;17:646–53.10.1161/01.ATV.17.4.646Search in Google Scholar

[3] Malek AM, Alper SL, Izumo S. Hemodynamic shear stress and its role in atherosclerosis. J Am Med Assoc. 1999;21:2035–42.10.1001/jama.282.21.2035Search in Google Scholar PubMed

[4] Asakura T, Karino T. Flow patterns and spatial distribution of atherosclerotic lesions in human coronary arteries. Circ Res. 1990;66:1045–66.10.1161/01.RES.66.4.1045Search in Google Scholar PubMed

[5] Lassen JF, Holm NR, Stankovic G, Lefèvre T, Chieffo A, Hildick-Smith D, et al. Percutaneous coronary intervention for coronary bifurcation disease: consensus from the first 10 years of the European Bifurcation Club meetings. EuroIntervention. 2014;10:545–60.10.4244/EIJV10I5A97Search in Google Scholar PubMed

[6] Malvè M, Gharib AM, Yazdani SK, Finet G, Martínez MA, Pettigrew R, et al. Tortuosity of Coronary bifurcation as a potential local risk factor for atherosclerosis: CFD steady state study based on in vivo dynamic CT measurements. Ann Biomed Eng. 2015;43:82–93.10.1007/s10439-014-1056-ySearch in Google Scholar PubMed PubMed Central

[7] Zhao W, Wang X, Jiang Y, Zhang J. Non-newtonian behavior of blood and arterial curvature influence variations of wall shear stress in stented artery. V European Conference on Computational Fluid Dynamics. Lisbon, Portugal, 14–17 June 2010.Search in Google Scholar

Published Online: 2016-9-30
Published in Print: 2016-9-1

©2016 Carolin Wüstenhagen et al., licensee De Gruyter.

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.

Articles in the same Issue

  1. Synthesis and characterization of PIL/pNIPAAm hybrid hydrogels
  2. Novel blood protein based scaffolds for cardiovascular tissue engineering
  3. Cell adhesion and viability of human endothelial cells on electrospun polymer scaffolds
  4. Effects of heat treatment and welding process on superelastic behaviour and microstructure of micro electron beam welded NiTi
  5. Long-term stable modifications of silicone elastomer for improved hemocompatibility
  6. The effect of thermal treatment on the mechanical properties of PLLA tubular specimens
  7. Biocompatible wear-resistant thick ceramic coating
  8. Protection of active implant electronics with organosilicon open air plasma coating for plastic overmolding
  9. Examination of dielectric strength of thin Parylene C films under various conditions
  10. Open air plasma deposited antimicrobial SiOx/TiOx composite films for biomedical applications
  11. Systemic analysis about residual chloroform in PLLA films
  12. A macrophage model of osseointegration
  13. Towards in silico prognosis using big data
  14. Technical concept and evaluation of a novel shoulder simulator with adaptive muscle force generation and free motion
  15. Usability evaluation of a locomotor therapy device considering different strategies
  16. Hypoxia-on-a-chip
  17. Integration of a semi-automatic in-vitro RFA procedure into an experimental setup
  18. Fabrication of MEMS-based 3D-μECoG-MEAs
  19. High speed digital interfacing for a neural data acquisition system
  20. Bionic forceps for the handling of sensitive tissue
  21. Experimental studies on 3D printing of barium titanate ceramics for medical applications
  22. Patient specific root-analogue dental implants – additive manufacturing and finite element analysis
  23. 3D printing – a key technology for tailored biomedical cell culture lab ware
  24. 3D printing of hydrogels in a temperature controlled environment with high spatial resolution
  25. Biocompatibility of photopolymers for additive manufacturing
  26. Biochemical piezoresistive sensors based on pH- and glucose-sensitive hydrogels for medical applications
  27. Novel wireless measurement system of pressure dedicated to in vivo studies
  28. Portable auricular device for real-time swallow and chew detection
  29. Detection of miRNA using a surface plasmon resonance biosensor and antibody amplification
  30. Simulation and evaluation of stimulation scenarios for targeted vestibular nerve excitation
  31. Deep brain stimulation: increasing efficiency by alternative waveforms
  32. Prediction of immediately occurring microsleep events from brain electric signals
  33. Determining cardiac vagal threshold from short term heart rate complexity
  34. Classification of cardiac excitation patterns during atrial fibrillation
  35. An algorithm to automatically determine the cycle length coverage to identify rotational activity during atrial fibrillation – a simulation study
  36. Deriving respiration from high resolution 12-channel-ECG during cycling exercise
  37. Reducing of gradient induced artifacts on the ECG signal during MRI examinations using Wilcoxon filter
  38. Automatic detection and mapping of double potentials in intracardiac electrograms
  39. Modeling the pelvic region for non-invasive pelvic intraoperative neuromonitoring
  40. Postprocessing algorithm for automated analysis of pelvic intraoperative neuromonitoring signals
  41. Best practice: surgeon driven application in pelvic operations
  42. Vasomotor assessment by camera-based photoplethysmography
  43. Classification of morphologic changes in photoplethysmographic waveforms
  44. Novel computation of pulse transit time from multi-channel PPG signals by wavelet transform
  45. Efficient design of FIR filter based low-pass differentiators for biomedical signal processing
  46. Nonlinear causal influences assessed by mutual compression entropy
  47. Comparative study of methods for solving the correspondence problem in EMD applications
  48. fNIRS for future use in auditory diagnostics
  49. Semi-automated detection of fractional shortening in zebrafish embryo heart videos
  50. Blood pressure measurement on the cheek
  51. Derivation of the respiratory rate from directly and indirectly measured respiratory signals using autocorrelation
  52. Left cardiac atrioventricular delay and inter-ventricular delay in cardiac resynchronization therapy responder and non-responder
  53. An automatic systolic peak detector of blood pressure waveforms using 4th order cumulants
  54. Real-time QRS detection using integrated variance for ECG gated cardiac MRI
  55. Preprocessing of unipolar signals acquired by a novel intracardiac mapping system
  56. In-vitro experiments to characterize ventricular electromechanics
  57. Continuous non-invasive monitoring of blood pressure in the operating room: a cuffless optical technology at the fingertip
  58. Application of microwave sensor technology in cardiovascular disease for plaque detection
  59. Artificial blood circulatory and special Ultrasound Doppler probes for detecting and sizing gaseous embolism
  60. Detection of microsleep events in a car driving simulation study using electrocardiographic features
  61. A method to determine the kink resistance of stents and stent delivery systems according to international standards
  62. Comparison of stented bifurcation and straight vessel 3D-simulation with a prior simulated velocity profile inlet
  63. Transient Euler-Lagrange/DEM simulation of stent thrombosis
  64. Automated control of the laser welding process of heart valve scaffolds
  65. Automation of a test bench for accessing the bendability of electrospun vascular grafts
  66. Influence of storage conditions on the release of growth factors in platelet-rich blood derivatives
  67. Cryopreservation of cells using defined serum-free cryoprotective agents
  68. New bioreactor vessel for tissue engineering of human nasal septal chondrocytes
  69. Determination of the membrane hydraulic permeability of MSCs
  70. Climate retainment in carbon dioxide incubators
  71. Multiple factors influencing OR ventilation system effectiveness
  72. Evaluation of an app-based stress protocol
  73. Medication process in Styrian hospitals
  74. Control tower to surgical theater
  75. Development of a skull phantom for the assessment of implant X-ray visibility
  76. Surgical navigation with QR codes
  77. Investigation of the pressure gradient of embolic protection devices
  78. Computer assistance in femoral derotation osteotomy: a bottom-up approach
  79. Automatic depth scanning system for 3D infrared thermography
  80. A service for monitoring the quality of intraoperative cone beam CT images
  81. Resectoscope with an easy to use twist mechanism for improved handling
  82. In vitro simulation of distribution processes following intramuscular injection
  83. Adjusting inkjet printhead parameters to deposit drugs into micro-sized reservoirs
  84. A flexible standalone system with integrated sensor feedback for multi-pad electrode FES of the hand
  85. Smart control for functional electrical stimulation with optimal pulse intensity
  86. Tactile display on the remaining hand for unilateral hand amputees
  87. Effects of sustained electrical stimulation on spasticity assessed by the pendulum test
  88. An improved tracking framework for ultrasound probe localization in image-guided radiosurgery
  89. Improvement of a subviral particle tracker by the use of a LAP-Kalman-algorithm
  90. Learning discriminative classification models for grading anal intraepithelial neoplasia
  91. Regularization of EIT reconstruction based on multi-scales wavelet transforms
  92. Assessing MRI susceptibility artefact through an indicator of image distortion
  93. EyeGuidance – a computer controlled system to guide eye movements
  94. A framework for feedback-based segmentation of 3D image stacks
  95. Doppler optical coherence tomography as a promising tool for detecting fluid in the human middle ear
  96. 3D Local in vivo Environment (LivE) imaging for single cell protein analysis of bone tissue
  97. Inside-Out access strategy using new trans-vascular catheter approach
  98. US/MRI fusion with new optical tracking and marker approach for interventional procedures inside the MRI suite
  99. Impact of different registration methods in MEG source analysis
  100. 3D segmentation of thyroid ultrasound images using active contours
  101. Designing a compact MRI motion phantom
  102. Cerebral cortex classification by conditional random fields applied to intraoperative thermal imaging
  103. Classification of indirect immunofluorescence images using thresholded local binary count features
  104. Analysis of muscle fatigue conditions using time-frequency images and GLCM features
  105. Numerical evaluation of image parameters of ETR-1
  106. Fabrication of a compliant phantom of the human aortic arch for use in Particle Image Velocimetry (PIV) experimentation
  107. Effect of the number of electrodes on the reconstructed lung shape in electrical impedance tomography
  108. Hardware dependencies of GPU-accelerated beamformer performances for microwave breast cancer detection
  109. Computer assisted assessment of progressing osteoradionecrosis of the jaw for clinical diagnosis and treatment
  110. Evaluation of reconstruction parameters of electrical impedance tomography on aorta detection during saline bolus injection
  111. Evaluation of open-source software for the lung segmentation
  112. Automatic determination of lung features of CF patients in CT scans
  113. Image analysis of self-organized multicellular patterns
  114. Effect of key parameters on synthesis of superparamagnetic nanoparticles (SPIONs)
  115. Radiopacity assessment of neurovascular implants
  116. Development of a desiccant based dielectric for monitoring humidity conditions in miniaturized hermetic implantable packages
  117. Development of an artifact-free aneurysm clip
  118. Enhancing the regeneration of bone defects by alkalizing the peri-implant zone – an in vitro approach
  119. Rapid prototyping of replica knee implants for in vitro testing
  120. Protecting ultra- and hyperhydrophilic implant surfaces in dry state from loss of wettability
  121. Advanced wettability analysis of implant surfaces
  122. Patient-specific hip prostheses designed by surgeons
  123. Plasma treatment on novel carbon fiber reinforced PEEK cages to enhance bioactivity
  124. Wear of a total intervertebral disc prosthesis
  125. Digital health and digital biomarkers – enabling value chains on health data
  126. Usability in the lifecycle of medical software development
  127. Influence of different test gases in a non-destructive 100% quality control system for medical devices
  128. Device development guided by user satisfaction survey on auricular vagus nerve stimulation
  129. Empirical assessment of the time course of innovation in biomedical engineering: first results of a comparative approach
  130. Effect of left atrial hypertrophy on P-wave morphology in a computational model
  131. Simulation of intracardiac electrograms around acute ablation lesions
  132. Parametrization of activation based cardiac electrophysiology models using bidomain model simulations
  133. Assessment of nasal resistance using computational fluid dynamics
  134. Resistance in a non-linear autoregressive model of pulmonary mechanics
  135. Inspiratory and expiratory elastance in a non-linear autoregressive model of pulmonary mechanics
  136. Determination of regional lung function in cystic fibrosis using electrical impedance tomography
  137. Development of parietal bone surrogates for parietal graft lift training
  138. Numerical simulation of mechanically stimulated bone remodelling
  139. Conversion of engineering stresses to Cauchy stresses in tensile and compression tests of thermoplastic polymers
  140. Numerical examinations of simplified spondylodesis models concerning energy absorption in magnetic resonance imaging
  141. Principle study on the signal connection at transabdominal fetal pulse oximetry
  142. Influence of Siluron® insertion on model drug distribution in the simulated vitreous body
  143. Evaluating different approaches to identify a three parameter gas exchange model
  144. Effects of fibrosis on the extracellular potential based on 3D reconstructions from histological sections of heart tissue
  145. From imaging to hemodynamics – how reconstruction kernels influence the blood flow predictions in intracranial aneurysms
  146. Flow optimised design of a novel point-of-care diagnostic device for the detection of disease specific biomarkers
  147. Improved FPGA controlled artificial vascular system for plethysmographic measurements
  148. Minimally spaced electrode positions for multi-functional chest sensors: ECG and respiratory signal estimation
  149. Automated detection of alveolar arches for nasoalveolar molding in cleft lip and palate treatment
  150. Control scheme selection in human-machine- interfaces by analysis of activity signals
  151. Event-based sampling for reducing communication load in realtime human motion analysis by wireless inertial sensor networks
  152. Automatic pairing of inertial sensors to lower limb segments – a plug-and-play approach
  153. Contactless respiratory monitoring system for magnetic resonance imaging applications using a laser range sensor
  154. Interactive monitoring system for visual respiratory biofeedback
  155. Development of a low-cost senor based aid for visually impaired people
  156. Patient assistive system for the shoulder joint
  157. A passive beating heart setup for interventional cardiology training
Downloaded on 21.4.2026 from https://www.degruyterbrill.com/document/doi/10.1515/cdbme-2016-0065/html
Scroll to top button