BLOOD FLOW SIMULATION USING SPH METHOD IN LS-DYNA, ANALYSIS OF ADVANTAGES AND DISADVANTAGES

1st International Conference on Chemo and BioInformatics, ICCBIKG  2021, (255-258)

AUTHOR(S) / AUTOR(I): Marko D. Topalovic, Aleksandar V. Nikolic, Miroslav M. Zivkovic

E-ADRESS / E-ADRESA: topalovic@kg.ac.rs , dziga@kg.ac.rs, zile@kg.ac.rs

Download Full Pdf   

DOI: 10.46793/ICCBI21.255T

ABSTRACT / SAŽETAK:

The purpose of this research was to investigate the possibility of blood flow modelling in LS-DYNA using its SPH solver and SPH-FEM coupling. SPH and FEM methods are both based on the continuum mechanics, and SPH uses Lagrangian material framework, while FEM can use both Lagrangian for solid, and Eulerian formulation for fluid analysis. SPH implementation is mesh-free giving it the capability to model very large deformations without mesh distortions. However, this comes at a high computational price, so the number of SPH particles needs to be significantly lower in comparison to the number of FEM elements in the Eulerian analysis of the same fluid domain. In the case of combined SPH-FEM analysis, the blood vessel wall is modelled with FEM shell elements, while the blood inside is modelled with SPH particles. The contact between the two is done using nodes to surface algorithm, while if we use the SPH only, there is no need for the specific contact definition. The Lagrangian framework of the SPH method means that we need to generate particles at one end, and to destroy them on the other, in order to generate a continuous fluid flow. To do this we used activation and deactivation planes, which is a solution implemented in the commercial LS-Dyna SPH solver. In the results section of the paper, the velocity field of blood obtained by implementation of described modelling methodology is shown. SPH-FEM coupling offers greater possibilities to study the effects of wall deformations, tracking of movement of solid particle inclusion, or mixing two different fluids, but it requires elaborate contact definition, and prolonged analysis time in comparison to the FEM CFD analysis.

KEY WORDS / KLJUČNE REČI:

SPH, FEM, Lagrangian formulation, Eulerian formulation, Fluid structure interaction

REFERENCES / LITERATURA:

  • Liu, M. Liu., Smoothed Particle Hydrodynamics, World Scientific Publishing, (2003).
  • Vignjević, J. Campbell., Brief review of development of the smooth particle hydrodynamics (SPH) method, IConSSM 2011, Vlasina lake, Serbia, 5-8 July (2011) 24-43.
  • A. Gingold, J. J. Monaghan., Smoothed particle hydrodynamics: theory and application to non- spherical stars, Monthly Notices Royal Astronomical Society, 181 (1997) 375-389.
  • B. Lucy., A numerical approach to the testing of fusion process, Astronomical Journal, 88 (1997) 1013-1024.
  • J. Monaghan, H, Pongracic., Artificial viscosity for particle methods, Applied Numerical Mathematics, 1 (1985) 187-194.
  • D. Libersky, A. G. Petschek, T. C. Carney, J. R. Hipp, F. A. Allahdadi., High Strain Lagrangian Hydrodynamics: A Three-Dimensional SPH Code for Dynamic Material Response, Journal of Computational Physics, 109 (1993) 67-75.
  • Kojic, N. Filipovic, B. Stojanovic, N. Kojic., Computer Modeling in Bioengineering: Theoretical Background, Examples and Software, (2009) John Wiley & Sons.
  • O. Hallquist., LS-DYNA Theory Manual. (2006) Livermore Software Technology Corporation.