Home Medicine Transient Euler-Lagrange/DEM simulation of stent thrombosis
Article Open Access

Transient Euler-Lagrange/DEM simulation of stent thrombosis

  • Michael Stiehm EMAIL logo , Carolin Wüstenhagen , Niels Grabow and Klaus-Peter Schmitz
Published/Copyright: September 30, 2016

Abstract

Stent implantation is the treatment of choice for cardiovascular diseases. By introduction of biodegradable thick strut stents investigations of thrombosis formation is one focus of research. This study deals with a transient Euler-Lagrange/DEM approach to simulate the flow field, platelet movement and clotting. The recirculation zones prolong particle residence time. As a result, the vicinity of stent struts shown a particularly higher risk for stent thrombosis.

1 Introduction

The treatment of choice for cardiovascular diseases initiated by atherosclerosis is the implantation of coronary stents. These implants open the affected artery to obtain the necessary blood flow rate.

Unfortunately the stent struts induce unphysiological flow conditions like recirculation zones, stasis and vortices resulting in post-operative complications. The most fatal one is the stent thrombosis which leads to myocardial infarction with high mortality rates. Especially with the introduction of biodegradable stents consisting of thick struts made of polymers or magnesium alloys thrombus formation is one focus of research again [1].

By means of numerical fluid dynamic simulations the understanding of thrombosis formation based on flow disturbances increases. Simulations give detailed information about flow distribution and platelet transport.

Blood is a suspension of viscoelastic particles like red blood cells and platelets in a Newtonian fluid, the blood plasma. The influence of the particles on fluid mechanical behaviour is often approximated by a non-Newtonian approach like the Carreau-Yasuda model [2]. With these models the shear-thinning behaviour of the whole-blood can be taken into account. Especially in regions with a shear rate below 100 s−1 this effect plays a significant role [3]. Nevertheless, the assumption of an ideal mixed and homogenous fluid neglects the transport mechanisms of thrombogenious particles and the adhesion of blood components.

Therefore Feng et al. simulated the viscous behaviour of blood in a stenosed two dimensional vessel with an advanced discrete particle dynamic (DPD) approach [4]. This model consists of multiple discrete particles interacting with each other. The results have shown a very good agreement of the DPD model and conventional computational fluid dynamic (CFD) method. Freund presented a review of further numerical models which also simulate the deformation of the red blood cells [5].

But all these discrete particle models request huge computational costs and so they are yet limited to capillary flow or simple cases.

Another numerical model, the so called Euler-Lagrange/DEM approach (discrete element model), combines the DPD and CFD method [6]. Therein the suspension is separated in a fluid and a solid phase. The fluid phase is treated as a homogeneous fluid with Newtonian or non-Newtonian flow behaviour and is calculated in the Eulerian framework as usual in common CFD simulations.

The solid phase consists of discrete zero dimensional particles and the movement is calculated in the Lagrangian framework based on Newton’s law. The interaction between both phases is considered by the source term in the momentum equation and the particle-particle interaction is often solved by a spring damper collision model. This kind of numerical model is successfully used by Chesnutt and Hanto simulate the process of stent thrombosis formation under steady state flow conditions [7].

The transient Euler-Lagrange/DEM model presented in this paper is based on the Lagrangian particle solver, which is implemented in the OpenFOAM framework [8].

2 Methods

The vessel is approximated by a rigid straight 2D channel with a height of h = 1.35 mm and an overall length of 10 h. The stent is modelled by 4 square struts with an edge length of 150 μm. Detailed information of stent design parameters can be seen in Figure 1.

Figure 1: Design parameters of the stent.
Figure 1:

Design parameters of the stent.

The particles are seeded in the red dotted box (Figure 1) every Δt = 0.04 s. In summation around 20,000 particles were modelled in one pulse (70 bpm). The particle count is 300,000 mm−3 [7]. This results in an average particle distance of 15 μm. The particle diameter is d = 2.4 μm.

The procedure of the Euler-Lagrange/DEM approach can be divided into two main steps.

  1. The Navier-Stokes equation (NSE) is solved for the whole domain in the Eulerian frame:

    U=0ρfUt+ρf(UU)=p+μf2USP

    where U is the velocity, ρ the density, p the pressure and μ the dynamic viscosity of the homogenous phase. The source term SP represents the forces on the solid particles and connects the Eulerian with the Lagrangian phase. Subscripts refer variables either to fluid phase or particles.

    For the inlet velocity condition a transient block profile is used. The mean Reynolds number is about Re = 250. The waveform was adopted from [9], see Figure 2.

  2. The equation of motion is solved for the inserted particles. The particle force is a result of the difference of the particle momentum in every control volume. The drag force FD influences the particle motion. By assuming a spherical particle, the drag force is calculated by:

    FD=12ρfCDUrel2AP

    where Urel = |U−Up|is the relative velocity and CD is the drag coefficient:

    CD=24Re

    Considering Strokes regime, which means the particle Reynolds number is below 0.1. The particle Reynolds number can be calculated from:

    ReP=ρfDPUrelμf

    These two steps for flow field and particle simulation can be done separately for a steady state flow or in one time step to regard a transient flow field as we used here.

Figure 2: Waveform used as inlet velocity condition based on [9].
Figure 2:

Waveform used as inlet velocity condition based on [9].

In this study only one way coupling is considered, so the source term in the momentum equation for the fluid phase is neglected. The particles are affected by the flow field via drag force but the particles do not influence the flow field.

The collision model implemented in the OpenFOAM solver is a non-linear viscoelastic model proposed by [10]. It’s based on the Hertz theory which means that particles are considered as elastic bodies, see Figure 3. Relevant simulation parameters can be found in Table 1.

Figure 3: Spring damper model used for particle collision.
Figure 3:

Spring damper model used for particle collision.

Table 1:

Simulation parameters and boundary conditions.

PropertyValue or condition
Artery radius or channel height1.35 mm
Strut height150 μm
Strut width150 μm
Mean fluid velocity0.31 m/s
Fluid density1060 kg/m3
Dynamic viscosity0.00355 kg/ms
Mean Reynolds number250
Inlet boundary conditionTransient velocity block profile
Outlet boundary conditionZero-pressure
Solid wall boundary conditionNo-slip
Centerline boundary conditionSymmetry
Number of particles1000 particles every Δt = 0.04 s
Particle diameter2.4 μm
Particle density1060 kg/m3

3 Results and discussion

3.1 Flow field

Numerous papers deal with the flow field in stented vessels. Therefore, only a short description of the velocity field is given here. As expected, recirculation zones can be identified up and down stream of the stent struts. Between both struts there exists one continuous recirculation zone, see Figure 4.

Figure 4 Velocity distribution around two stent struts during systole (t = 0.4 s).
Figure 4

Velocity distribution around two stent struts during systole (t = 0.4 s).

3.2 Formation of stent thrombosis

The thrombosis formation is initiated by the contact of separate particles due to spatial flow disturbances triggered by stent struts. At that moment two or more particles collide and merge the developed clot starts to tumble. Also the clot inertia causes a difference in velocity between separate particles and the clot, whereby more particles get in contact with the clot and fuse, see Figure 5.

Figure 5 Timeline of thrombosis formation upstream of a stent strut.
Figure 5

Timeline of thrombosis formation upstream of a stent strut.

The process of thrombosis formation in the flow field reinforces itself after ignition by the stent struts.

Additionally, the particle residence time increases enormous in the recirculation zones and so the risk of thrombosis increases. Figure 6 depicts particle positions and particle ages in the vicinity of stent struts. It can be seen, that a clot consists of particles with different ages. So, the particle laden flow isn’t layered anymore but mixed by the stent. Several computational and in vivo studies have shown that the vicinity of stent struts is particularly at risk [11].

Figure 6 Particle position in the vicinity of stent struts, particles are coloured by their age.
Figure 6

Particle position in the vicinity of stent struts, particles are coloured by their age.

4 Conclusion

In this paper a transient Euler-Lagrange/DEM approach for simulation of thrombosis formation is introduced. This approach allows simulations of transient flow field and discrete particle interaction in one time step. Particle collision is taken into account by a non-linear viscoelastic collision model. Biochemical mechanisms like platelet activation are not considered.

Stent struts induce flow disturbances like recirculation zones promoting platelet collision. As a result thrombosis formation can be initiated. Furthermore, particle residence time is prolonged in recirculation zones up and down stream of stent struts leading to higher risk of thrombosis formation.

This approach goes beyond common CFD simulation and is beneficial for thrombosis forecast in stented coronary vessels. Transient simulations of three dimensional cases are planned for further research.

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] Byrne RA, Kastrati A. Bioresorbable drug-eluting stents - an immature technology in need of mature application. JACC: Cardiovasc Interv. 2015;8:198–200.10.1016/j.jcin.2014.11.010Search in Google Scholar

[2] Leuprecht A, Perktold K. Computer simulation of non-newtonian effects on blood flow in large arteries. Comput Meth. 2001;4:149–63.10.1080/10255840008908002Search in Google Scholar

[3] Fung YC. Biomechanics: Circulation. 2nd ed, New York: Springer, 1997. p. 571.10.1007/978-1-4757-2696-1Search in Google Scholar

[4] Feng R, Xenos M, Girdhar G, Kang W, Davenport JW, Deng Y, et al. Viscous flow simulation in a stenosis model using discrete particle dynamics: a comparison between DPD and CFD. Biomech Model Mechanobiol. 2012;11:119–29.10.1007/s10237-011-0297-zSearch in Google Scholar

[5] Freund JB. Numerical Simulation of Flowing Blood Cells. Annu Rev Fluid Mech. 2014;46:67–95.10.1146/annurev-fluid-010313-141349Search in Google Scholar

[6] Oevermann M, Gerber S, Behrendt F. Euler–Lagrange/DEM simulation of wood gasification in a bubbling fluidized bed reactor. Particuology. 2009;7:307–16.10.1016/j.partic.2009.04.004Search in Google Scholar

[7] Chesnutt JKW, Han HC. Simulation of the microscopic process during initiation of stent thrombosis. Comput Biol Med. 2015;56:182–91.10.1016/j.compbiomed.2014.11.006Search in Google Scholar

[8] Vallier A. Simulations of cavitation – from the large vapour structures to the small bubble dynamics. Dissertation. 2013Search in Google Scholar

[9] HuoY, Choy JS, Svendsen M, Sinha AK, Kassab GS. Effects of vessel compliance on flow pattern in porcine epicardial right coronary arterial tree. J Biomech. 2009;42:594–602.10.1016/j.jbiomech.2008.12.011Search in Google Scholar

[10] Tsuji Y, Tanaka T, Ishida T. Lagrangian numerical simulation of plug flow of cohesionless particles in a horizontal pipe. Powder Technol. 1992;71:239–50.10.1016/0032-5910(92)88030-LSearch in Google Scholar

[11] Kolandaivelu K, Swaminathan R, Gibson WJ, Kolachalama VB, Nguyen-Ehrenreich KL, Giddings VL et al. Stent thrombogenicity early in high-risk interventional settings is driven by stent design and deployment and protected by polymer-drug coatings. Circulation. 2011;123:1400–9.10.1161/CIRCULATIONAHA.110.003210Search in Google Scholar PubMed PubMed Central

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

©2016 Michael Stiehm 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 5.3.2026 from https://www.degruyterbrill.com/document/doi/10.1515/cdbme-2016-0066/html
Scroll to top button